Lack of evidence for ectopic sprouting of genetically labeled Aβ touch afferents in inflammatory and neuropathic trigeminal pain
© Zhang et al.; licensee BioMed Central. 2015
Received: 27 February 2015
Accepted: 30 March 2015
Published: 10 April 2015
Mechanical and in particular tactile allodynia is a hallmark of chronic pain in which innocuous touch becomes painful. Previous cholera toxin B (CTB)-based neural tracing experiments and electrophysiology studies had suggested that aberrant axon sprouting from touch sensory afferents into pain-processing laminae after injury is a possible anatomical substrate underlying mechanical allodynia. This hypothesis was later challenged by experiments using intra-axonal labeling of A-fiber neurons, as well as single-neuron labeling of electrophysiologically identified sensory neurons. However, no studies have used genetically labeled neurons to examine this issue, and most studies were performed on spinal but not trigeminal sensory neurons which are the relevant neurons for orofacial pain, where allodynia oftentimes plays a dominant clinical role.
We recently discovered that parvalbumin::Cre (Pv::Cre) labels two types of Aβ touch neurons in trigeminal ganglion. Using a Pv::CreER driver and a Cre-dependent reporter mouse, we specifically labeled these Aβ trigeminal touch afferents by timed taxomifen injection prior to inflammation or infraorbital nerve injury (ION transection). We then examined the peripheral and central projections of labeled axons into the brainstem caudalis nucleus after injuries vs controls. We found no evidence for ectopic sprouting of Pv::CreER labeled trigeminal Aβ axons into the superficial trigeminal noci-receptive laminae. Furthermore, there was also no evidence for peripheral sprouting.
CreER-based labeling prior to injury precluded the issue of phenotypic changes of neurons after injury. Our results suggest that touch allodynia in chronic orofacial pain is unlikely caused by ectopic sprouting of Aβ trigeminal afferents.
Mechanical allodynia, in which light touch is perceived as painful , is a hallmark of many chronic pain conditions. Numerous previous studies reported changes that were interpreted as ectopic collateral sprouting of Aβ touch afferents into the “pain-processing” superficial laminae I-II after nerve injury [2-6]. The abnormal sprouting was observed initially using bulk transport of cholera toxin B (CTB), and later using single-fiber labeling and electrophysiology [2,5,6]. The use of CTB as a specific marker taken up selectively by large-diameter Aβ neurons in the periphery was later found to be problematic [7-10]. Furthermore, intra-axonal recording followed by intra-axonal injection of the neural tracer into the identified Aβ fibers revealed a complete lack of evidence for sprouting of putative Aβ touch afferents (low-threshold mechanoreceptors, or LTMRs) into the superficial laminae . Another evidence argue against Aβ fibers sprouting came from the finding that the plexus of VGLUT1-immunoreactive boutons in laminae IIi-III of the lumbar dorsal horn (most of which originate from Aβ LTMR afferents) did not extend more dorsally after sciatic nerve transection . More recently, Koerber and colleagues used an ex vivo prep to record from the soma of previously axotomized peripheral neurons, then labeled the recorded neurons with intrasomal staining. This approach also demonstrated that individually labeled Aβ LTMR afferents do not sprout into “pain” laminae after injury, but instead, they discovered that myelinated nociceptive neurons (high threshold mechanoreceptors HTMRs) form flame-shaped terminal arbors that project throughout lamina I-V. They concluded that these myelinated HTMRs could account for the morphological findings previously assigned to sprouting LTMRs .
These previous studies were performed on dorsal root ganglion (DRG) sensory neurons. We want to examine trigeminal ganglion (TG) sensory neurons that innervate the head, face and its adjacent sentient structures such as teeth, sinuses, dura mater, cornea, temporo-mandibular joint (TMJ). TG sensory neurons are most relevant for orofacial pain and headaches, forms of pain that share a clinical hallmark of significant mechanical allodynia, e.g. in trigeminal neuralgia, atypical face pain, dental pulpitis, keratitis, TMJ disorder, postherpetic trigeminal nerve pain and migraine [14-21]. We therefore felt compelled to probe the neuro-anatomical basis of mechanical allodynia in the TG system. In animal models, TG nerve injury shares similarities with its DRG counterpart, but there are also significant differences . Previously, trigeminal nerve transection studies had reported changes in central axon arborizations after transection, but these transection studies were performed in fetal or neonatal mice [23,24]. Similar experiments in adult animals had not been performed. Furthermore, we wanted to use a genetic method to label Aβ TG touch neurons of defined morphological types before injury, in order to determine for these precisely defined Aβ touch afferent whether they sprout ectopically into noci-receptive laminae in the brainstem after injury.
Painful orofacial cues are detected by nociceptive TG neurons and relayed to the brainstem trigeminal spinal nucleus, the spinal caudalis (SpC). Similarly to the spinal cord dorsal horn, C-fiber peptidergic and non-peptidergic TG nociceptive neurons project to the superficial layers of SpC including lamina I and lamina II , whereas TG LTMR neurons project to deeper laminae in SpC (III –V). Importantly, our recent experiments revealed that in the mouse TG system, parvalbumin::Cre (Pv::Cre) selectively labeled two types of Aβ touch (LTMR) neurons: primarily slowly adapting Merkel-ending neurons and a small number of rapidly adapting longitudinal lanceolate-ending neurons . In the current study we also use Pv::Cre and Pv::CreER (in which tamoxifen-inducible CreER is knocked into the Pv locus)  mouse lines and subject them to two orofacial pain models in order to probe peripheral and central axonal projections of TG Aβ afferents in pathological conditions.
Results and discussion
Using neonatal tamoxifen injection in Pv-CreER; Ai14 mice to specifically label individual TG Aβ touch afferents
No evidence of ectopic central sprouting of TG Aβ touch afferents in a CFA-induced trigeminal inflammatory pain model
We first carried out experiments using an inflammatory pain model. Numerous previous studies showed that injection of complete Freund’s adjuvant (CFA) into the whisker pad (which are innervated by the infraorbital branch of the maxillary nerve) resulted in prolonged mechanical allodynia that lasted 10 days or longer [29-35]. We repeated the CFA inflammation experiments using the “Pv:CreER; Ai14”- mice after neonatal tamoxifen injection at P5. After CFA injection into the whisker pad, we used von Frey hairs to examine mechanical sensitivity at different time points. In keeping with previous experiences, we observed that 3 ~ 4 days following initial inflammatory injury, mechanical allodynia was most pronounced, and it lasted more than 10 days (Figure 1A). Thus, we chose to examine the axon projections at day 3 and day 10 post CFA injection, to determine whether there would be any transient ectopic sprouting (at day 3) or sprouting at late stage (at day 10) (Figure 1B and F, schematic of experimental procedure). CFA was injected into the left whisker pad of adult “Pv-CreER; Ai14”-mice (that had tamoxifen injection at P5). The whisker pad on the right side was injected with IFA and served as internal control. We used isolectin B4 (IB4) to mark the inner lamina of the “pain-processing” layers in SpC. IB4 is a marker for non-peptidergic nociceptive fibers that project to lamina II of SpC, same as for spinal sensory neurons . Figure 1C-E shows representative images of SpC sections 3 days after CFA injection. Figure 1C is an overview of the SpC. The left side of the SpC is CFA-injected and the right side is control (Figure 1C). High magnification images from the regions receiving inputs from the whisker pad (boxed area) are shown in Figure 1D-E. Other than the main incoming axons from labeled TG en route to deep layers, we did not observe any tomato-labeled Aβ axon collaterals or boutons sprouting into either the IB4 positive lamina II or the more superficial lamina I, where peptidergic nociceptive fibers terminate. This was observed in both the CFA-treated and the control side (Figure 1C-E, N = 5 mice). In samples collected 10 days after CFA injection (Figure 1F), we also failed to observe any ectopic sprouting of tomato-labeled axons into the noci-receptive layers of the SpC (Figure 1G-I, N = 5). Therefore, CFA-induced inflammatory pain did not cause abnormal central sprouting of Pv-expressing LTMR TG afferent axons into the superficial laminae.
No evidence of ectopic central sprouting of TG Aβ touch afferents in a trigeminal neuropathic pain model
P5 “Pv::CreER; Ai14”-mice were injected with tamoxifen to specifically label individual TG Aβ afferents before injury. Subsequently, ION transection was performed on the left side, and sham operation was carried out on the right side in 4-week adult mice. 3, 14, 22 days after injury, mice were perfused and the SpC region of the samples were serial sectioned and stained with IB4-488. At all time points examined, we did not observe any labeled axon terminals sprouting into the IB4 positive or more superficial layers of the SpC (Figures 2 and 3, N = 5 mice for each group). High magnification images of the whisker representing regions in SpC are shown, and there were no differences observed between the control side and the injured side (Figures 2 and 3).
Lack of peripheral sprouting of TG Aβ touch afferents in both CFA-induced inflammatory and neuropathic pain models
Previous studies using intra-axonal labeling of DRG neurons, combined with electrophysiology did not support the hypothesis that Aβ afferents ectopically sprout into the superficial noci-receptive laminae of the spinal cord [11,13]. However, whether there was transient Aβ sprouting remains unknown because the analyses were performed at one time point in the later stages after nerve injury. There is also a possibility that some LTMR Aβ neurons undergo a complete phenoptypic change into HTMRs that project throughout lamina I-V after injury. These two important issues can be resolved by genetically pre-labeling of Aβ afferents and examining the central projections of labeled neurons after injury.
In this study, we selectively labeled a subset of Aβ low-threshold trigeminal touch neurons prior to inflammatory or neuropathic injuries, using a genetic labeling strategy relying on the Pv:CreER mouse line. Our genetic method allowed clear visualization of central axonal arborizations from the Aβ LTMR touch neurons. Subsequently, the labeled animals were examined at several time points either after inflammatory or nerve ligation-mediated injury. We found no evidence of ectopic peripheral or central axonal sprouting, either transiently or in more chronic conditions in all animals examined. Thus, the long lasting mechanical allodynia in our orofacial pain models is unlikely caused by ectopic Aβ axonal sprouting into “pain-processing” layers of the trigeminal spinal nucleus. Since we only labeled Pv-expressing Aβ TG neurons, we cannot rule out the possibility that other types of Aβ LTMR touch neurons that were not labeled in our study sprout axons after injury. However, our results together with previous studies strongly argue against the ectopic sprouting hypothesis. Instead, we favor the idea that gene expression changes in primary sensory neurons and/or sensitization or disinhibition of central processing circuits contribute more significantly to the tactile hypersensitivity [40,41].
Materials and methods
Pv::CreER; Ai14 (Cre-dependent TdTomato reporter) and Pv-Cre; RΦAP adult male mice of at least 4-week old were used for injury models and subsequent axon projection studies. Both alleles are heterozygous. Pv::Cre mice (stock number 008069) and Ai14 (stock number 007914) were purchased from the Jackson Laboratories. RΦAP and Pv::CreER mice have been described previously [26,42]. All animal procedures were approved by The Duke University Institutional Animal Care and Use Committee.
The tamoxifen solution dissolved in corn oil (0.05 mg/g; body weight) was injected subcutaneously in Pv-CreER; Ai14 pups at P5.
Complete Freund’s Adjuvant (CFA) injection
All mice were briefly anesthetized with 2% isoflurane and injections were performed using a 30-G needle on a 10 μl Hamilton syringe. 10 μl of complete Freund adjuvant (CFA, 5 mg/ml; Chondrex, Redmond, WA) were injected into the left (ipsilateral) whisker pad. Control sides were injected with the same volume of incomplete Freund’s adjuvant (IFA, Chondrex). Following CFA treatment, the mice were subjected for behavior tests, and were sacrificed either at day 3 and 10 respectively for histology analysis.
Infraorbital Nerve (ION) injury model
Mice were anesthetized with Ketamine/Xylazine (Sigma, 40 mg/3 mg/kg, i.p.), and all surgeries were performed in sterile conditions under a surgical microscope [4,9]. The infraorbital muscle was gently dissected from the bone until the orbit could be gently retracted. The infraorbital nerve on the left was dissected free from the bone at its most rostral extent in the orbital cavity. The nerve was then ligated with chromic gut suture (6–0, Angiotech) and transected at the point just distal to the ligature. For the sham operation on the right side, only skin incision and muscle dissection were performed, and the nerve was not touched and no chromic gut suture was inserted. All skin incisions were sutured with 5–0 nylon non-absorbable monofilament (Angiotech), and mice were subjected for behavior tests after surgery, and were sacrificed on day 3, 14 and 22 respectively for histology analysis.
Mice after CFA injection or ION transection injury were subject to von Frey hair testing. Briefly, two different von Frey filaments (0.4 and 2.0 g) were used with the smallest hair used first. The response to each of the filaments is scored as follows: score 0, no response; score 1, detection - the mice detected and explored the von Frey hair; score 2, head slowly withdrawal reaction; score 3, escape/attack; score 4, asymmetric face grooming – at least three face wash strokes directed toward the stimulated facial area . Three applications were performed for each filament and the mean of the measurements is used. All data were expressed as mean ± SEM. Differences between groups were evaluated using 2-tailed Student’s t test (experimental against sham control), or in the case of multiple groups, one-way ANOVA followed by post-hoc Bonferroni test. The criterion for statistical significance was p < 0.05.
IB4-staing and placental alkaline phosphatase staining (PLAP staining) were performed using standard procedures. For PLAP staining, the tissue sections were incubated in PBS at 65°C for 3 hours, and developed with 1:200 NBT/BCIP (Roche) in staining solution which contains 0.1 M Tris–HCl, pH 9.5, 0.1 M NaCl and 5 mM MgCl2. GS-IB4 Alexa 488-conjugated (Invitrogen) was used at 1:1000.
This study was supported by NIH grants DE19440 (FW), DE19440S1 (FW and WL), DE018549 (WL), F33DE024668 (YC), support from the Duke Institute for Brain Science (DIBS) to WL and FW, a Leon Goldberg Fellowship (Nicholas School for the Environment, Duke University) to YC. We thank Bao-Xia Han for technical assistance.
- Campbell JN, Raja SN, Meyer RA, Mackinnon SE. Myelinated afferents signal the hyperalgesia associated with nerve injury. Pain. 1988;32(1):89–94.View ArticlePubMedGoogle Scholar
- Kohama I, Ishikawa K, Kocsis JD. Synaptic reorganization in the substantia gelatinosa after peripheral nerve neuroma formation: aberrant innervation of lamina II neurons by Abeta afferents. J Neurosci. 2000;20(4):1538–49.PubMed CentralPubMedGoogle Scholar
- Shortland P, Woolf CJ. Chronic peripheral nerve section results in a rearrangement of the central axonal arborizations of axotomized A beta primary afferent neurons in the rat spinal cord. J Comp Neurol. 1993;330(1):65–82.View ArticlePubMedGoogle Scholar
- Woolf CJ, Shortland P, Coggeshall RE. Peripheral nerve injury triggers central sprouting of myelinated afferents. Nature. 1992;355(6355):75–8.View ArticlePubMedGoogle Scholar
- Koerber HR, Mirnics K, Brown PB, Mendell LM. Central sprouting and functional plasticity of regenerated primary afferents. J Neurosci. 1994;14(6):3655–71.PubMedGoogle Scholar
- Koerber HR, Mirnics K, Kavookjian AM, Light AR. Ultrastructural analysis of ectopic synaptic boutons arising from peripherally regenerated primary afferent fibers. J Neurophysiol. 1999;81(4):1636–44.PubMedGoogle Scholar
- Shehab SAS, Spike RC, Todd AJ. Evidence against cholera toxin B subunit as a reliable tracer for sprouting of primary afferents following peripheral nerve injury. Brain Res. 2003;964(2):218–27.View ArticlePubMedGoogle Scholar
- Tong YG, Wang HF, Ju G, Grant G, Hokfelt T, Zhang X. Increased uptake and transport of cholera toxin B-subunit in dorsal root ganglion neurons after peripheral axotomy: possible implications for sensory sprouting. J Comp Neurol. 1999;404(2):143–58.View ArticlePubMedGoogle Scholar
- Bao L, Wang HF, Cai HJ, Tong YG, Jin SX, Lu YJ, et al. Peripheral axotomy induces only very limited sprouting of coarse myelinated afferents into inner lamina II of rat spinal cord. Eur J Neurosci. 2002;16(2):175–85.View ArticlePubMedGoogle Scholar
- Santha P, Jancso G. Transganglionic transport of choleragenoid by capsaicin-sensitive C-fibre afferents to the substantia gelatinosa of the spinal dorsal horn after peripheral nerve section. Neuroscience. 2003;116(3):621–7.View ArticlePubMedGoogle Scholar
- Hughes DI, Scott DT, Todd AJ, Riddell JS. Lack of evidence for sprouting of A beta afferents into the superficial laminas of the spinal cord dorsal horn after nerve section. J Neurosci. 2003;23(29):9491–9.PubMedGoogle Scholar
- Hughes DI, Polgar E, Shehab SA, Todd AJ. Peripheral axotomy induces depletion of the vesicular glutamate transporter VGLUT1 in central terminals of myelinated afferent fibres in the rat spinal cord. Brain Res. 2004;1017(1–2):69–76.View ArticlePubMedGoogle Scholar
- Woodbury CJ, Kullmann FA, McIlwrath SL, Koerber HR. Identity of myelinated cutaneous sensory neurons projecting to nocireceptive laminae following nerve injury in adult mice. J Comp Neurol. 2008;508(3):500–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Benemei S, De Cesaris F, Fusi C, Rossi E, Lupi C, Geppetti P. TRPA1 and other TRP channels in migraine. J Headache Pain. 2013;14:71.View ArticlePubMed CentralPubMedGoogle Scholar
- Oaklander AL. Mechanisms of pain and itch caused by herpes zoster (shingles). J Pain. 2008;9(1 Suppl 1):S10–8.View ArticlePubMedGoogle Scholar
- Noseda R, Burstein R. Migraine pathophysiology: anatomy of the trigeminovascular pathway and associated neurological symptoms, cortical spreading depression, sensitization, and modulation of pain. Pain. 2013;154 Suppl 1:S44–53.View ArticlePubMedGoogle Scholar
- Ashkenazi A. Allodynia in cluster headache. Curr Pain Headache Rep. 2010;14(2):140–4.View ArticlePubMedGoogle Scholar
- Schurks M, Diener HC. Migraine, allodynia, and implications for treatment. Eur J Neurol. 2008;15(12):1279–85.View ArticlePubMedGoogle Scholar
- Oshinsky ML. Insights from experimental studies into allodynia and its treatment. Curr Pain Headache Rep. 2006;10(3):225–30.View ArticlePubMedGoogle Scholar
- Moskowitz MA. Neurogenic inflammation in the pathophysiology and treatment of migraine. Neurology. 1993;43(6 Suppl 3):S16–20.PubMedGoogle Scholar
- Rosenthal P, Borsook D. The corneal pain system. Part I: the missing piece of the dry eye puzzle. Ocul Surf. 2012;10(1):2–14.View ArticlePubMedGoogle Scholar
- Fried K, Bongenhielm U, Boissonade FM, Robinson PP. Nerve injury-induced pain in the trigeminal system. Neuroscientist. 2001;7(2):155–65.View ArticlePubMedGoogle Scholar
- Rhoades RW, Chiaia NL, Macdonald GJ, Jacquin MF. Effect of fetal infraorbital nerve transection upon trigeminal primary afferent-projections in the Rat. J Comp Neurol. 1989;287(1):82–97.View ArticlePubMedGoogle Scholar
- Waite PME, dePermentier PJ. Effect of neonatal capsaicin and infraorbital nerve section on whisker-related patterns in the rat trigeminal nucleus. J Comp Neurol. 1997;385(4):599–615.View ArticlePubMedGoogle Scholar
- Bae YC, Oh JM, Hwang SJ, Shigenaga Y, Valtschanoff JG. Expression of vanilloid receptor TRPV1 in the rat trigeminal sensory nuclei. J Comp Neurol. 2004;478(1):62–71.View ArticlePubMedGoogle Scholar
- Sakurai K, Akiyama M, Cai B, Scott A, Han BX, Takatoh J, et al. The organization of submodality-specific touch afferent inputs in the vibrissa column. Cell Rep. 2013;5(1):87–98.View ArticlePubMedGoogle Scholar
- Madisen L, Zwingman TA, Sunkin SM, Oh SW, Zariwala HA, Gu H, et al. A robust and high-throughput Cre reporting and characterization system for the whole mouse brain. Nat Neurosci. 2010;13(1):133–40.View ArticlePubMed CentralPubMedGoogle Scholar
- Bennett-Clarke CA, Chiaia NL, Jacquin MF, Rhoades RW. Parvalbumin and calbindin immunocytochemistry reveal functionally distinct cell groups and vibrissa-related patterns in the trigeminal brainstem complex of the adult rat. J Comp Neurol. 1992;320(3):323–38.View ArticlePubMedGoogle Scholar
- Krzyzanowska A, Pittolo S, Cabrerizo M, Sanchez-Lopez J, Krishnasamy S, Venero C, et al. Assessing nociceptive sensitivity in mouse models of inflammatory and neuropathic trigeminal pain. J Neurosci Methods. 2011;201(1):46–54.View ArticlePubMedGoogle Scholar
- Benoliel R, Wilensky A, Tal M, Eliav E. Application of a pro-inflammatory agent to the orbital portion of the rat infraorbital nerve induces changes indicative of ongoing trigeminal pain. Pain. 2002;99(3):567–78.View ArticlePubMedGoogle Scholar
- Chen Y, Williams SH, McNulty AL, Hong JH, Lee SH, Rothfusz NE, et al. Temporomandibular joint pain: a critical role for Trpv4 in the trigeminal ganglion. Pain. 2013;154(8):1295–304.View ArticlePubMed CentralPubMedGoogle Scholar
- Morgan JR, Gebhart GF. Characterization of a model of chronic orofacial hyperalgesia in the rat: contribution of NA(V) 1.8. J Pain. 2008;9(6):522–31.View ArticlePubMedGoogle Scholar
- Okumura M, Iwata K, Yasuda K, Inoue K, Shinoda M, Honda K, et al. Alternation of gene expression in trigeminal ganglion neurons following complete Freund’s adjuvant or capsaicin injection into the Rat face. J Mol Neurosci. 2010;42(2):200–9.View ArticlePubMedGoogle Scholar
- Garrett FG, Hawkins JL, Overmyer AE, Hayden JB, Durham PL. Validation of a novel Rat-holding device for studying heat- and mechanical-evoked trigeminal nocifensive behavioral responses. J Orofac Pain. 2012;26(4):337–44.PubMedGoogle Scholar
- Simonic-Kocijan S, Zhao XH, Liu W, Wu YW, Uhac I, Wang KW. TRPV1 channel-mediated bilateral allodynia induced by unilateral masseter muscle inflammation in rats. Mol Pain. 2013;9. doi:10.1186/1744-8069-9-68.Google Scholar
- Wang H, Rivero-Melian C, Robertson B, Grant G. Transganglionic transport and binding of the isolectin B4 from Griffonia simplicifolia I in rat primary sensory neurons. Neuroscience. 1994;62(2):539–51.View ArticlePubMedGoogle Scholar
- Kernisant M, Gear RW, Jasmin L, Vit JP, Ohara PT. Chronic constriction injury of the infraorbital nerve in the rat using modified syringe needle. J Neurosci Methods. 2008;172(1):43–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Alvarez P, Dieb W, Hafidi A, Voisin DL, Dallel R. Insular cortex representation of dynamic mechanical allodynia in trigeminal neuropathic rats. Neurobiol Dis. 2009;33(1):89–95.View ArticlePubMedGoogle Scholar
- Lim EJ, Jeon HJ, Yang GY, Lee MK, Ju JS, Han SR, et al. Intracisternal administration of mitogen-activated protein kinase inhibitors reduced mechanical allodynia following chronic constriction injury of infraorbital nerve in rats. Prog Neuro-Psychoph. 2007;31(6):1322–9.View ArticleGoogle Scholar
- Woolf CJ. Central sensitization: implications for the diagnosis and treatment of pain. Pain. 2011;152(3 Suppl):S2–15.View ArticlePubMed CentralPubMedGoogle Scholar
- Kuner R. Central mechanisms of pathological pain. Nat Med. 2010;16(11):1258–66.View ArticlePubMedGoogle Scholar
- Que J, Wilm B, Hasegawa H, Wang F, Bader D, Hogan BL. Mesothelium contributes to vascular smooth muscle and mesenchyme during lung development. Proc Natl Acad Sci U S A. 2008;105(43):16626–30.View ArticlePubMed CentralPubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.