- Jiangsu Province Key Laboratory of Anesthesiology, Xuzhou Medical University, Xuzhou, China
Objective: To investigate whether the CSF-contacting nucleus receives brainstem and spinal cord projections and to understand the functional significance of these connections.
Methods: The retrograde tracer cholera toxin B subunit (CB) was injected into the CSF-contacting nucleus in Sprague-Dawley rats according the previously reported stereotaxic coordinates. After 7–10 days, these rats were perfused and their brainstem and spinal cord were sliced (thickness, 40 μm) using a freezing microtome. All the sections were subjected to CB immunofluorescence staining. The distribution of CB-positive neuron in different brainstem and spinal cord areas was observed under fluorescence microscope.
Results: The retrograde labeled CB-positive neurons were found in the midbrain, pons, medulla oblongata, and spinal cord. Four functional areas including one hundred and twelve sub-regions have projections to the CSF-contacting nucleus. However, the density of CB-positive neuron distribution ranged from sparse to dense.
Conclusion: Based on the connectivity patterns of the CSF-contacting nucleus receives anatomical inputs from the brainstem and spinal cord, we preliminarily conclude and summarize that the CSF-contacting nucleus participates in pain, visceral activity, sleep and arousal, emotion, and drug addiction. The present study firstly illustrates the broad projections of the CSF-contacting nucleus from the brainstem and spinal cord, which implies the complicated functions of the nucleus especially for the unique roles of coordination in neural and body fluids regulation.
Introduction
The cerebrospinal fluid (CSF)-contacting nucleus is a unique nucleus in the brain. It is located within the ventral gray of the lower portion of the aqueduct (Aq) and upper portion of the fourth ventricle (4V) floor (Song et al., 2019). The outstanding feature of this nucleus is that the neural somata are located in the brain parenchyma with their axons stretching into the CSF (Song and Zhang, 2018; Song et al., 2019). The morphological connections of CSF-contacting nucleus with non-CSF-contacting neurons, glial cells, and blood vessels have been confirmed using electron microscopy (Zhang et al., 2003). The unique characteristic feature of the CSF-contacting nucleus is that it may be a key structure bridging the nerves and fluids (CSF and plasma). The study of this nucleus is about 30 years after we discovered and named this nucleus. The basic biological characteristics of the CSF-contacting nucleus such as specific labeling method (Lu et al., 2008), location and stereotaxic coordinates (Song et al., 2019), receptor, neurotransmitter, and ion channel distributions (Lu et al., 2011; Wang et al., 2014; Liu et al., 2017), and their relationships with morphine dependence and withdrawal, stress, sodium appetite, and pain (Lu et al., 2011; Wu et al., 2015; Xing et al., 2015; Zhou et al., 2017) have been revealed. However, the pathways and mechanisms associated with such biological activities have not been elucidated yet.
The CSF-contacting nucleus is also located in the brainstem and participates in pain and other behaviors (Lu et al., 2011; Xing et al., 2015; Liu et al., 2017). To demonstrate the functions and neural networks of the CSF-contacting nucleus in the central nervous system, it is import to identify the existence of connections from the brainstem and spinal cord to the CSF-contacting nucleus and to elucidate the connections from different functional regions. Therefore, we planned to inject retrograde tracer cholera toxin B subunit (CB) into the CSF-contacting nucleus. By using immunofluorescence technique, the projections from the brainstem and spinal cord to the CSF-contacting nucleus can be revealed. The possible functional significance can be speculated according to the projection relationships, which will deepen the understanding and lay a foundation for future research.
Materials and Methods
Experimental Animals
Specific pathogen-free (SPF) grade Sprague-Dawley rats weighing 250 ± 50 g were acquired from the Experimental Animal Center of Xuzhou Medical University. Rats successfully injected with the tracer into the CSF-contacting nucleus were used for observation and analysis (n = 6). All the animals were housed under a 12/12 light/dark cycle at the temperature of 22–25°C, with ad libitum access to food and water.
Tracer Administration
Rats were anesthetized with pentobarbital sodium (40 mg/kg, i.p.), and their heads were fixed on the stereotaxic instrument (Stoelting 51700, United States). Retrograde tracer CB solution (0.2 μl, 1%; Sigma, United States) was injected into the core of the CSF-contacting nucleus (Posterior to Bregma: 8.24 ± 0.18 mm, Lateral: 0.09 ± 0.01 mm, Depth: 6.45 ± 0.11 mm) (Song et al., 2019) by using Hamilton syringe with a 33G needle tip (Hamilton Company, Switzerland). The injection was applied under the microinfusion pump (KD Scientific, United States) for over 30 min periods. At the end of the injection, the microsyringe was left in place for 10–15 min before retraction.
Sampling and Sectioning
After 7–10 days of retrograde tracer CB injections, rats were perfused transcardially and sacrificed. The whole brain and spinal cord were isolated and sectioned (40 μm thickness) coronally using a cryostat (Leica CM1900, Germany). All the sections were placed in sequence and numbered. In this study, only the brainstem and spinal cord were analyzed.
Tracer Staining and CB-Positive Neuron Enumeration
After CB immunofluorescent staining (rabbit anti-CB primary antibody diluted in 1:600, Abcam; donkey anti-rabbit Alexa Fluor 488 secondary antibody diluted in 1:200, Life Technologies), the sections were sequentially mounted on slides, counterstained with DAPI, and coverslipped. The sections of the brainstem and spinal cord were visualized using fluorescence (Leica DM6, Germany) and confocal laser microscopes (Zeiss LSM 880, Germany) for different positions at the uniform standard. The cell density of CB-positive neurons (cell number/0.2 mm2 area) in each brain region was calculated using Image-Pro Plus 7.0 software. The density of CB-positive neurons was classified according to the densities: <5, 6–10, and >10, which correspond to sparse, moderate, and dense distributions, respectively.
Verification of the Retrograde Results by Anterograde Tracing
The recombinant adenoassociated virus vector (serotype 2/9) expressed enhanced GFP (EGFP) under control of a human synapsin I promoter was used for anterograde tracing to verify the retrograde tracing results. 0.2 μl rAAV2/9-hSyn-EGFP (1.2 × 1012 gene copies; Brainvta Company, Wuhan, China) was injected into the raphe magnus nucleus (RMg) (Posterior to Bregma: 10.92 ± 0.2 mm, Lateral: 0 mm, Depth: 10.5 ± 0.1 mm) (Paxinos and Watson, 2007). 1 μl (1 μg) Cholera toxin subunit B Alexa Fluor 555 conjugate (CB-555) (Brainvta Company, Wuhan, China) was microinjected into the lateral ventricle (Posterior to Bregma: 1.0 ± 0.1 mm, Lateral: 1.6 ± 0.1 mm, Depth: 4.0 ± 0.2 mm) (Paxinos and Watson, 2007) to label the CSF-contacting nucleus. The injection parameters were the same as above. Rats were perfused after 14–21 days of anterograde tracer rAAV2/9 injections. The axon terminals from the RMg to the CSF-contacting nucleus were captured under confocal laser microscopes (Zeiss LSM 880, Germany).
3D Reconstruction of the Brainstem and Spinal Cord Connections
The CB-positive neurons were aligned, segmented, and registered into rat common reference atlas (Paxinos and Watson, 2007). The 3D connections were reconstructed by using Imaris software version 8.4.1 (Bitplane, United States).
Statistics
SPSS 13.0 software was used for data analysis in the present study. Data are presented as mean ± SD.
Results
Injection of Retrograde Tracer Cholera Toxin B Subunit Into the CSF-Contacting Nucleus
Injections of the CB tracer produced dense immunofluorescence staining (green). The tracer was confined within the boundary of the CSF-contacting nucleus (Figures 1A,B), where the microsyringe needle tract can be seen to be located at the core of the CSF-contacting nucleus. Representative sections of the CSF-contacting nucleus are shown in Figures 1C,D.
Figure 1. Photographs showing injection of the tracer CB into the CSF-contacting nucleus. (A) The site of injection of tracer CB into the CSF-contacting nucleus. Green fluorescence labeling of CB can be seen in the wholly CSF-contacting nucleus. (B) Higher magnification of the boxed area in (A). Downwards arrow (↓) shows the passage of the injected needle. (C) The representative section of the CSF-contacting nucleus in the brain. (D) Higher magnification of the boxed area in C. Aq, aqueduct. Scale bar = 1 mm in (A), 100 μm in (B), 70 μm in (C), and 40 μm in (D).
Cellular Morphology of Brainstem and Spinal Cord Connections
After the retrograde tracer was injected into the CSF-contacting nucleus, it was transported in a retrograde manner along the axons; consequently, the neural somata projected from the brainstem and spinal cord can be detected.
The CB-positive neurons were round or fusiform in shape. Sizes of these neurons are different, and the processes are obvious (Figure 2A). In the reticular structure of the brainstem, the neurons are large and appear in polygonal shape; they have many processes, and the processes have many branches (Figures 2B,C). In the spinal cord, the neurons are fusiform, and the processes are obvious and longer (Figure 2D).
Figure 2. Cellular morphology of CB-positive neurons. Representative photographs from the midbrain (A), pons (B), medulla oblongata (C), and spinal cord (D). Scale bar = 40 μm.
Connection Sites of the Brainstem and Spinal Cord
The brainstem and spinal cord projections to the CSF-contacting nucleus can be identified by CB labeling.
In the midbrain, CB-positive neuron can be found in the precommissural nucleus (PrC), periaqueductal gray (PAG), nucleus of Darkschewitsch (DK), medial accessory oculomotor nucleus (MA3), pre-Edinger-Westphal nucleus (PrEW), Edinger-Westphal nucleus (EW), rostral linear nucleus of the raphe (RLi), interfascicular nucleus (IF), interstitial nucleus of Cajal (InC), red nucleus (R), ventral tegmental area (VTA), substantia nigra (SN), medial pretectal nucleus (MPT), posterior pretectal nucleus (PPT), olivary pretectal nucleus (OPT), nucleus of the optic tract (OT), posterior thalamic nuclear group (Po), interpeduncular nucleus (IP), supraoculomotor cap (Su3C), supraoculomotor periaqueductal gray (Su3), oculomotor nucleus (3N), trochlear nucleus (4N), mesencephalic reticular formation (mRt), superior colliculus (SC), inferior colliculus (IC), caudal linear nucleus of the raphe (CLi), subbrachial nucleus (SubB), and retrorubral field (RRF) (Figures 3, 4).
Figure 3. The distribution of neurons labeled by CB in the midbrain Part I (A–D). PrC, precommissural nucleus; VTA, ventral tegmental area; SN, substantia nigra; mRt, mesencephalic reticular formation; SubB, subbrachial nucleus; MA3, medial accessory oculomotor nucleus; PrEW, pre-Edinger-Westphal nucleus. Scale bar = 100 μm.
Figure 4. The distribution of neurons labeled by CB in the midbrain Part II (A–D). EW, Edinger-Westphal nucleus; Su3, supraoculomotor periaqueductal gray; Su3C, supraoculomotor cap; 3N, oculomotor nucleus; 4N, trochlear nucleus; IP, interpeduncular nucleus; PAG, periaqueductal gray. Scale bar = 100 μm.
Among them, the PrC, PAG, MA3, IP, Su3C, Su3, 3N, 4N, and SubB have dense projections; the PrEW, EW, IF, SN, MPT, PPT, OPT, OT, and mRt send moderate projections; and the DK, RLi, InC, R, VTA, Po, SC, IC, CLi, and RRF send sparse projections to the CSF-contacting nucleus (Figures 3, 4).
In the pons, CB-positive neuron can be found in the pontine nuclei (Pn), dorsal raphe nucleus (DR), median raphe nucleus (MnR), paramedian raphe nucleus (PMnR), isthmic reticular formation (isRt), pedunculopontine tegmental nucleus (PTg), subpeduncular tegmental nucleus (SPTg), microcellular tegmental nucleus (MiTg), precuneiform area (PrCnF), cuneiform nucleus (CnF), laterodorsal tegmental nucleus (LDTg), medial paralemniscial nucleus (MPL), nucleus of the lateral lemniscus (LL), A7 noradrenaline cells (A7), central gray (CG), pontine reticular nucleus oral part (PnO), pontine reticular nucleus caudal part (PnC), pontine reticular nucleus ventral part (PnV), Barrington’s nucleus (Bar), subcoeruleus nucleus (SubC), lateral parabrachial nucleus (LPB), medial parabrachial nucleus (MPB), Kolliker Fuse nucleus (KF), supratrigeminal nucleus (Su5), motor nucleus of trigeminal nerve (Mo5), medioventral periolivary nucleus (MVPO), lateroventral periolivary nucleus (LVPO), locus coeruleus (LC), superior vestibular nucleus (SuVe), medial vestibular nucleus parvicellular part (MVePC), medial vestibular nucleus magnocellular part (MVeMC), lateral vestibular nucleus (LVe), spinal vestibular nucleus (SpVe), nucleus of origin of efferents of the vestibular nerve (EVe), supragenual nucleus (SGe), ventral cochlear nucleus (VC), dorsal cochlear nucleus (DC), lateral superior olive (LSO), superior paraolivary nucleus (SPO), nucleus of the trapezoid body (Tz), abducens nucleus (6N), paraabducens nucleus (Pa6), raphe interpositus nucleus (RIP), raphe pallidus nucleus (RPa), dorsal paragigantocellular nucleus (DPGi), gigantocellular reticular nucleus (Gi), gigantocellular reticular nucleus alpha part (GiA), gigantocellular reticular nucleus ventral part (GiV), raphe magnus nucleus (RMg), facial nucleus (7N), A5 noradrenaline cells (A5), intermediate reticular nucleus alpha part (IRtA), parvicellular reticular nucleus alpha part (PCRtA), parapyramidal nucleus (PPy), lateral paragigantocellular nucleus alpha part (LPGiA), and lateral paragigantocellular nucleus external part (LPGiE) (Figure 5).
Figure 5. The distribution of neurons labeled by CB in the pons (A–F). DR, dorsal raphe nucleus; CnF, cuneiform nucleus; MnR, median raphe nucleus; PMnR, paramedian raphe nucleus; Pn, pontine nuclei; DTg, dorsal tegmental nucleus; LDTg, laterodorsal tegmental nucleus; Bar, Barrington’s nucleus; scp, superior cerebellar peduncle; LPB, lateral parabrachial nucleus; MPB, medial parabrachial nucleus; KF, Kolliker Fuse nucleus. Scale bar = 100 μm.
Among them, the Pn, DR, MnR, CnF, LDTg, CG, Bar, LPB, KF, LC, MVePC, EVe, SGe, DC, and 6N have strong projections; the PMnR, PTg, MPB, MVPO, SuVe, MVeMC, VC, TZ, RPa, RMg, LPGiA, and LPGiE have moderate projections; and the isRt, SPTg, MiTg, PrCnF, MPL, LL, A7, PnO, PnC, PnV, SubC, Su5, Mo5, LVPO, LVe, LSO, SPO, Pa6, RIP, DPGi, Gi, GiA, GiV, 7N, A5, IRtA, PCRtA, and PPy have sparse projections to the CSF-contacting nucleus (Figure 5).
In the medulla oblongata and spinal cord, CB-positive neurons can be found in the prepositus nucleus (Pr), solitary nucleus (Sol), raphe obscurus nucleus (Rob), epifascicular nucleus (EF), rostroventrolateral reticular nucleus (RVL), caudoventrolateral reticular nucleus (CVL), spinal trigeminal nucleus (Sp5), ambiguus nucleus (Amb), Botzinger complex (Bo), nucleus of Roller (Ro), external cuneate nucleus (ECu), cuneate nucleus (Cu), inferior olive (IO), dorsal motor nucleus of vagus (10N), hypoglossal nucleus (12N), rostral ventral respiratory group (RVRG), paratrigeminal nucleus (Pa5), lateral reticular nucleus (LRt), gracile nucleus (Gr), area postrema (AP), A1 noradrenaline cells (A1), A2 noradrenaline cells (A2), accessory nerve nucleus (11N), intermediate reticular nucleus (IRt), parvicellular reticular nucleus (PCRt), medullary reticular nucleus dorsal part (MdD), medullary reticular nucleus ventral part (MdV), and spinal cord (Figures 6–10).
Figure 6. The distribution of neurons labeled by CB in the medulla oblongata Part I (A–D). 6N, abducens nucleus; Pa6, paraabducens nucleus; DC, dorsal cochlear nucleus; Pr, prepositus nucleus; MVePC, medial vestibular nucleus parvicellular part; MVeMC, medial vestibular nucleus magnocellular part; Sol, solitary nucleus; Sp5, spinal trigeminal nucleus. Scale bar = 100 μm.
Figure 7. The distribution of neurons labeled by CB in the medulla oblongata Part II (A–D). Gi, gigantocellular reticular nucleus; GiA, gigantocellular reticular nucleus alpha part; RMg, raphe magnus nucleus; RPa, raphe pallidus nucleus; py, pyramidal tract; 10N, dorsal motor nucleus of vagus; 12N, hypoglossal nucleus; Sol, solitary nucleus. Scale bar = 100 μm.
Figure 8. The distribution of neurons labeled by CB in the medulla oblongata Part III (A–D). 10N, dorsal motor nucleus of vagus; 12N, hypoglossal nucleus; Sol, solitary nucleus; AP, area postrema. Scale bar = 100 μm.
Figure 9. The distribution of neurons labeled by CB in the medulla oblongata Part IV (A–E). Rob, raphe obscurus nucleus; IO, inferior olive; LRt, lateral reticular nucleus; Amb, ambiguus nucleus; CVL, caudoventrolateral reticular nucleus; 11N, accessory nerve nucleus. Scale bar = 100 μm.
Among them, the Pr, Sol, RVL, CVL, Amb, 10N, 12N, AP, and A1 have strong projections; the Rob, Sp5, Bo, IO, LRt, A2, 11N, and IRt have moderate projections; and EF, Ro, ECu, Cu, RVRG, Pa5, Gr, PCRt, MdD, MdV, and spinal cord have sparse projections to the CSF-contacting nucleus (Figures 6–10).
In summary, CB-positive neurons were distributed in 4 functional areas including 112 sub-regions in the brainstem and spinal cord. However, their distribution ranged from sparse to dense in each functional region. Other parts did not contain CB-positive neurons. In addition, the retrograde results were further verified by the anterograde method. After the anterograde tracer rAAV2/9 injection into the RMg, large number of axon terminals can be detected in the CSF-contacting nucleus (Figure 11).
Figure 11. The verification of the RMg projection to the CSF-contacting nucleus by anterograde tracing. (A) The injection site of the anterograde tracer rAAV2/9-hSyn-EGFP into the RMg. (B) The axon terminals of the RMg (green) make close appositions with the neurons in the CSF-contacting nucleus (red). RMg, raphe magnus nucleus; py, pyramidal tract; ml, medial lemniscus; Aq, aqueduct. Scale bar = 100 μm in A; Scale bar = 50 μm in B.
3D Reconstruction of the CB-Positive Neurons From the Brainstem and Spinal Cord
The distributions of CB-positive neurons throughout the brainstem and spinal cord were three dimensionally reconstructed. The densities of the connections were obvious in the 3D view. The red areas indicated dense connections (PrC, PAG, MA3, IP, Su3C, Su3, 3N, 4M, SubB, Pn, DR, MnR, CnF, LDTg, CG, Bar, LPB, KF, LC, MVePC, EVe, SGe, DC, 6N, Pr, Sol, RVL, CVL, Amb, 10N, 12N, AP, and A1); green areas indicated moderate connections (PrEW, EW, IF, SN, MPT, PPT, OPT, OT, mRt, PMnR, PTg, MPB, MVPO, SuVe, MVeMC, VC, Tz, RPa, RMg, LPGiA, LPGiE, ROb, Sp5, Bo, IO, LRt, A2, 11N, and IRt); and blue areas indicated sparse connections (DK, RLi, InC, R, VTA, Po, SC, IC, CLi, RRF, isRt, SPTg, MiTg, MPL, LL, A7, PnO, PnC, PnV, SubC, Su5, Mo5, LVPO, LVe, SpVe, LSO, SPO, Pa6, RIP, DPGi, Gi, GiA, GiV, 7N, A5, IRtA, PCRtA, PPy, EF, Ro, ECu, Cu, RVRG, Pa5, Gr, PCRt, MdD, MdV, and spinal cord) to the CSF-contacting nucleus (Figure 12).
Figure 12. 3D view of the patterns of the brainstem and spinal cord connections to the CSF-contacting nucleus. (A) Dorsal view, (B) posterior view, (C) lateral view, (D) ventral view. The red areas indicate strong connections; green areas indicate moderate connections; and blue areas indicate weak connections to the CSF-contacting nucleus.
The Amount of Projection From the Brainstem and Spinal Cord to the CSF-Contacting Nucleus
In the brainstem and spinal cord, the CB-positive neurons were found in 112 sub-regions. The amount of projection from these regions to the CSF-contacting nucleus is shown in Figure 13.
Figure 13. Whole brainstem and spinal cord statistics of CB-positive cell input to the CSF-contacting nucleus (mean ± SD, n = 6).
Discussion
The CSF-contacting nucleus is a unique nucleus in the brain. It has non-synaptic connections with blood vessels and CSF via the CSF-contacting neurons; this connection plays an important role in the regulation of body fluids. Additionally, it has synaptic connections with non-CSF-contacting neurons via the CSF-contacting neurons; this connection facilitates neuron-neuron crosstalk in the brain. The unique anatomical feature of the CSF-contacting nucleus is that it may be a key structure bridging the nervous and humoral regulating systems. It has been uncovered that this nucleus receives the projection from central nervous system, although the bidirectional synapses between CSF-contacting and non-CSF-contacting neurons were observed in the parenchyma by electron microscopy (Liang et al., 2007). Herein, we include a systematic discussion of the projections from different functional zones in the brainstem and spinal cord to the CSF-contacting nucleus. As the fibers of passage may have been taken up by the retrograde tracer CB, we also verify the results by the anterograde method.
Our results indicate that the CSF-contacting nucleus received extensive projections from 4 functional areas including 112 sub-regions of the brainstem and spinal cord (Figure 14). Among them, the CSF-contacting nucleus receives the connections from the brainstem cranial nerve nuclei, such as trigeminal nuclei and solitary nucleus. The CSF-contacting nucleus is also located in the brainstem, but the processes stretch into the CSF. The nucleus is more sensitive to the peripheral tracer CB-HRP or CB which implies that the CSF-contacting nucleus has the characteristics of the cranial nerve nuclei. According to the connection patterns of the CSF-contacting nucleus with the functional areas of the brainstem and spinal cord, the biological functions of the CSF-contacting nucleus can be predicted. Whilst the function of the CSF-contacting nucleus is not known, the following are the functional implications of these connections.
Figure 14. The schematic diagram of projections from functional areas in the brainstem and spinal cord to the CSF-contacting nucleus. Among them, midbrain contains 28 sub-regions, pons contains 56 sub-regions, and medulla oblongata and spinal cord contain 28 sub-regions.
Pain
The pain descending inhibitory system is located in the brainstem. This system can send fibers in a top-down manner to the dorsal spinal horn to inhibit pain transmission (Millan, 2002). The system comprises of the PAG, DR, rostral-ventromedial medulla (RVM), and LC (Reynolds, 1969; Francois et al., 2017; Torralba et al., 2018; Zhang et al., 2018). The system also sends dense connections to the CSF-contacting nucleus, which might be a novel pathway in pain modulation.
Apart from the descending inhibitory system, other regions (CnF, PTg, A5, A7, LRt, MdD, and spinal cord) (Zemlan and Behbehani, 1988; Dugast et al., 2003; Mravec et al., 2012; Huma et al., 2015; Francois et al., 2017; de Oliveira et al., 2018) can also regulate nociception and send projections to the CSF-contacting nucleus. For example, deep brain stimulation of the CnF neurons is known to produce antinociceptive effect (Zemlan and Behbehani, 1988). Lesion in the PTg impair the antinociceptive effects (de Oliveira et al., 2018).
Visceral Activity
In the brainstem, the Bar, CnF, Gi, LPB, MPB, PTg, 10N, A1, Amb, Bo, RVL, RVRG, and Sol are important for visceral activity modulation; they have extensive connections with the CSF-contacting nucleus. The Bar engages in modulating the bladder functions (Blanco et al., 2014). Manipulation of the neuronal functions in the CnF, Gi, PTg, A1, Amb, and RVL has great effects on cardiovascular responses (Chan and Chan, 1983; Elliott et al., 1985; Balan Junior et al., 2004; Topchiy et al., 2010; Shafei et al., 2012; Brailoiu et al., 2017). The Bo and RVRG are responsible for respiratory control. The Bo has been described to be involved in the compensatory responses to hypoxia (Alberto et al., 2014), while neurons of the RVRG can directly project to the phrenic motor neurons and discharge an inspiratory pattern (Hayakawa et al., 2004).
The LPB and MPB are located around the superior cerebellar peduncle (scp). They can receive various visceral inputs from the Sol and spinal cord and can relay the visceral information to other regions of the brain (Palmiter, 2018). The Sol is one of the most important sites involved in visceral modulation, and it is regarded as the entry point of afferent signals from various visceral mechanosensors, chemosensors, and C-fibers (Kotmanova et al., 2018; You et al., 2018). The 10N can integrate information and control the visceral activities (Browning et al., 2005; Cruz et al., 2007).
Sleep and Arousal
The CSF-contacting nucleus receives broad projections along the brainstem reticular formation—a region of the brain that plays a key role in regulating sleep and arousal (Moruzzi and Magoun, 1949). The PnO is a part of the ascending reticular activating system, and it regulates rapid eye movement sleep (REMS) and wakefulness (Nunez et al., 2010; Silvia et al., 2013; Vanini and Baghdoyan, 2013). The Gi has projections to the midbrain reticular formation and intralaminar thalamic nuclei, which help modulate arousal (Martin et al., 2011). During the REMS, activity of the LRt neurons changes significantly (Stettner et al., 2013).
Apart from the brainstem reticular formation, other regions of the brain that innervate the CSF-contacting nucleus, such as the DR and LC, can also modulate sleep and arousal. Pharmacological intervention the DR can significantly change REMS (Monti and Jantos, 2018). The LC is a key arousal node, and it is pivotal for modulating consciousness and effects of general anesthetics (Du et al., 2018).
Emotion
The CSF-contacting nucleus receives projections from the EW, IP, PAG, DR, LC, and MnR, which might be associated with emotion. The EW is a region implicated in stress responses and different aspects of stress-related behaviors (Xu et al., 2010, 2014). Lesion-based studies demonstrated that the IP is involved in regulating anxiety and fear (Mclaughlin et al., 2017). The PAG is regarded as an integrated site for coordinating rapid responses to aversive stimuli, such as those provoking panic, anxiety, or fear (Lowery-Gionta et al., 2018). The DR and MnR contain the major serotonergic (5-HT) populations (Pollak Dorocic et al., 2014). Both of the nuclei and 5-HT are involved in anxiety and depression, and they are crucial and central drug targets (Lanzenberger et al., 2012; An et al., 2013, 2016; Lopez Hill et al., 2013; Campos et al., 2018; de Souza et al., 2018; Donner et al., 2018). The LC is the major noradrenergic system with regard to stress and major depressive disorder (Chandley and Ordway, 2012; Fan et al., 2018).
Drug Addiction
The CSF-contacting nucleus receives input from the EW, IP, LDTg, and LPGi, which might be associated with drug addiction. After oral self-administration of ethanol, the EW shows significantly elevated c-Fos expression (Bachtell et al., 1999; Giardino et al., 2011). Pharmacological studies show that the IP mediates behavioral responses to nicotine reinforcement, aversion, and withdrawal (Ramiro et al., 2009; Antolin-Fontes et al., 2015; Morton et al., 2018). The LDTg is involved in drug reward processing and reinforcement in the development of addiction and abuse (Kaneda, 2018). The LPGi has been shown to play an important critical role in the development of opiate tolerance and dependence (Ahmadi-Soleimani et al., 2017).
In this study, we used tract-tracing method to reveal the CSF-contacting nucleus input patterns from the brainstem and spinal cord. The unique morphological feature of the CSF-contacting nucleus is that the somata are located in the brain parenchyma, which can receive inputs from the above-mentioned areas; the processes can form synaptic and non-synaptic connections with non-CSF-contacting neurons, CSF, or plasma. For one aspect, it might form the brainstem and spinal cord→CSF-contacting nucleus→non-CSF-contacting neurons circuits and participate in the regulation of life activities via neuron-neuron crosstalk. For another aspect, it might be involved in the brainstem and spinal cord→CSF-contacting nucleus→CSF/plasma circuit and modulate physiological functions via neuron-fluid interactions. According to the connection regularities of CSF-contacting nucleus from the brainstem and spinal cord, we preliminary conclude that the CSF-contacting nucleus involved in modulation pain, visceral activity, sleep and arousal, emotion, and drug addiction, etc. The current study provides morphological evidence for unveiling the significance of the CSF-contacting nucleus.
Data Availability Statement
The data that support the findings of this study are available from the corresponding author, upon reasonable request.
Ethics Statement
All animal experiments were approved by and performed in accordance with the guidelines of the Committee for Ethical Use of Laboratory Animals of Xuzhou Medical University.
Author Contributions
S-YS and L-CZ designed the study and prepared the manuscript. S-YS, YL, X-MZ, Y-HL, C-YB, C-JS, JH, and J-LC conducted the studies. All authors read and approved the manuscript.
Funding
This research is supported by the National Natural Science Foundation of China (Grant Nos. 81371243 and 81901131), the Natural Science Foundation of Jiangsu Province (Grant No. BK20190987) and Chinese Postdoctoral Science Foundation (Grant No. 2018M642328).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank all members of L-CZ group at Xuzhou Medical University.
References
Ahmadi-Soleimani, S. M., Azizi, H., Gompf, H. S., and Semnanian, S. (2017). Role of orexin type-1 receptors in paragiganto-coerulear modulation of opioid withdrawal and tolerance: a site specific focus. Neuropharmacology 126, 25–37. doi: 10.1016/j.neuropharm.2017.08.024
Alberto, F. D., Soriano, R. N., Sabino, J. P., and Branco, L. G. S. (2014). Endogenous hydrogen sulfide in the rostral ventrolateral medulla/Botzinger complex downregulates ventilatory responses to hypoxia. Respir. Physiol. Neurobiol. 200, 97–104. doi: 10.1016/j.resp.2014.06.007
An, Y., Inoue, T., Kitaichi, Y., Chen, C., Nakagawa, S., Wang, C., et al. (2016). Combined treatment with subchronic lithium and acute intracerebral mirtazapine microinjection into the median raphe nucleus exerted an anxiolytic-like effect synergistically. Eur. J. Pharmacol. 783, 112–116. doi: 10.1016/j.ejphar.2016.05.004
An, Y., Inoue, T., Kitaichi, Y., Izumi, T., Nakagawa, S., Song, N., et al. (2013). Anxiolytic-like effect of mirtazapine mediates its effect in the median raphe nucleus. Eur. J. Pharmacol. 720, 192–197. doi: 10.1016/j.ejphar.2013.09.078
Antolin-Fontes, B., Ables, J. L., Gorlich, A., and Ibanez-Tallon, I. (2015). The habenulo-interpeduncular pathway in nicotine aversion and withdrawal. Neuropharmacology 96, 213–222. doi: 10.1016/j.neuropharm.2014.11.019
Bachtell, R. K., Wang, Y. M., Freeman, P., Risinger, F. O., and Ryabinin, A. E. (1999). Alcohol drinking produces brain region-selective changes in expression of inducible transcription factors. Brain Res. 847, 157–165. doi: 10.1016/s0006-8993(99)02019-3
Balan Junior, A., Caous, C. A., Yu, Y. G., and Lindsey, C. J. (2004). Barosensitive neurons in the rat tractus solitarius and paratrigeminal nucleus: a new model for medullary, cardiovascular reflex regulation. Can. J. Physiol. Pharmacol. 82, 474–484. doi: 10.1139/y04-054
Blanco, L., Ros, C. M., Tarragon, E., Fernandez-Villalba, E., and Herrero, M. T. (2014). Functional role of Barrington’s nucleus in the micturition reflex: relevance in the surgical treatment of Parkinson’s disease. Neuroscience 266, 150–161. doi: 10.1016/j.neuroscience.2014.02.002
Brailoiu, E., McGuire, M., Shuler, S. A., Deliu, E., Barr, J. L., Abood, M. E., et al. (2017). Modulation of cardiac vagal tone by bradykinin acting on nucleus ambiguus. Neuroscience 365, 23–32. doi: 10.1016/j.neuroscience.2017.09.034
Browning, K. N., Coleman, F. H., and Travagli, R. A. (2005). Characterization of pancreas-projecting rat dorsal motor nucleus of vagus neurons. Am. J. Physiol. Gastrointestinal Liver Physiol. 288, G950–G955.
Campos, G. V., De, S. N., De, A. S., Lima, P. M., Abreu, A. R., Chianca, J. D., et al. (2018). Estrogen receptor β activation within dorsal raphe nucleus reverses anxiety-like behavior induced by food restriction in female rats. Behav. Brain Res. 357–358, 57–64. doi: 10.1016/j.bbr.2018.03.029
Chan, S. H., and Chan, J. Y. (1983). Correlated effects of clonidine on single-neuron activities in the gigantocellular reticular nucleus, arterial pressure and heart rate in the cat. Neurosci. Lett. 40, 139–143. doi: 10.1016/0304-3940(83)90292-6
Chandley, M. J., and Ordway, G. A. (2012). “Noradrenergic dysfunction in depression and suicide,” in The Neurobiological Basis of Suicide, ed. Y. Dwivedi (Boca Raton, FL: CRC Press).
Cruz, M. T., Murphy, E. C., Sahibzada, N., Verbalis, J. G., and Gillis, R. A. (2007). A reevaluation of the effects of stimulation of the dorsal motor nucleus of the vagus on gastric motility in the rat. Am. J. Physiol. Regul. Integr. Compar. Physiol. 292, R291–R307.
de Oliveira, R. C., de Oliveira, R., Falconi-Sobrinho, L. L., Biagioni, A. F., Almada, R. C., Dos Anjos-Garcia, T., et al. (2018). Neurotoxic lesions of the pedunculopontine tegmental nucleus impair the elaboration of postictal antinociception. Physiol. Behav. 194, 162–169. doi: 10.1016/j.physbeh.2018.05.011
de Souza, M. M., Silote, G. P., Herbst, L. S., Funck, V. R., Joca, S. R. L., and Beijamini, V. (2018). The antidepressant-like effect of galanin in the dorsal raphe nucleus of rats involves GAL2 receptors. Neurosci. Lett. 681, 26–30. doi: 10.1016/j.neulet.2018.05.029
Donner, N. C., Kubala, K. H., Hassell, J. E. Jr., Lieb, M. W., Nguyen, K. T., Heinze, J. D., et al. (2018). Two models of inescapable stress increase tph2 mRNA expression in the anxiety-related dorsomedial part of the dorsal raphe nucleus. Neurobiol. Stress 8, 68–81. doi: 10.1016/j.ynstr.2018.01.003
Du, W. J., Zhang, R. W., Li, J., Zhang, B. B., Peng, X. L., Cao, S., et al. (2018). The locus coeruleus modulates intravenous general anesthesia of zebrafish via a cooperative mechanism. Cell Rep. 24. 3146–3155.e3143.
Dugast, C., Almeida, A., and Lima, D. (2003). The medullary dorsal reticular nucleus enhances the responsiveness of spinal nociceptive neurons to peripheral stimulation in the rat. Eur. J. Neurosci. 18, 580–588. doi: 10.1046/j.1460-9568.2003.02782.x
Elliott, J. M., Stead, B. H., West, M. J., and Chalmers, J. (1985). Cardiovascular effects of intracisternal 6-hydroxydopamine and of subsequent lesions of the ventrolateral medulla coinciding with the Al group of noradrenaline cells in the rabbit. J. Auton. Nerv. Syst. 12, 117–130. doi: 10.1016/0165-1838(85)90055-4
Fan, Y., Chen, P., Raza, M. U., Szebeni, A., Szebeni, K., Ordway, G. A., et al. (2018). Altered expression of Phox2 transcription factors in the locus coeruleus in major depressive disorder mimicked by chronic stress and corticosterone treatment in vivo and in vitro. Neuroscience 393, 123–137. doi: 10.1016/j.neuroscience.2018.09.038
Francois, A., Low, S. A., Sypek, E. I., Christensen, A. J., Sotoudeh, C., Beier, K. T., et al. (2017). A brainstem-spinal cord inhibitory circuit for mechanical pain modulation by GABA and enkephalins. Neuron 93, 822–839.e826.
Giardino, W. J., Cocking, D. L., Kaur, S., Cunningham, C. L., and Ryabinin, A. E. (2011). Urocortin-1 within the centrally-projecting edinger-westphal nucleus is critical for ethanol preference. PLoS ONE 6:e26997. doi: 10.1371/journal.pone.0026997
Hayakawa, T., Takanaga, A., Tanaka, K., Maeda, S., and Seki, M. (2004). Ultrastructure of the rostral ventral respiratory group neurons in the ventrolateral medulla of the rat. Brain Res. 1027, 94–102. doi: 10.1016/j.brainres.2004.08.059
Huma, Z., Ireland, K., and Maxwell, D. J. (2015). The spino-bulbar-cerebellar pathway: activation of neurons projecting to the lateral reticular nucleus in the rat in response to noxious mechanical stimuli. Neurosci. Lett. 591, 197–201. doi: 10.1016/j.neulet.2015.02.047
Kaneda, K. (2018). Neuroplasticity in cholinergic neurons of the laterodorsal tegmental nucleus contributes to the development of cocaine addiction. Eur. J. Neurosci. 50, 2239–2246. doi: 10.1111/ejn.13962
Kotmanova, Z., Simera, M., Veternik, M., Martvon, L., Misek, J., Jakus, J., et al. (2018). GABA-ergic neurotransmission in the nucleus of the solitary tract modulates cough in the cat. Respir. Physiol. Neurobiol. 257, 100–106. doi: 10.1016/j.resp.2018.02.009
Lanzenberger, R., Kranz, G. S., Haeusler, D., Akimova, E., Savli, M., Hahn, A., et al. (2012). Prediction of SSRI treatment response in major depression based on serotonin transporter interplay between median raphe nucleus and projection areas. NeuroImage 63, 874–881. doi: 10.1016/j.neuroimage.2012.07.023
Liang, D., Zhang, L., Qin, C., and Zeng, Y. (2007). The ultrastructure of the distal CSF-contacting neurons in the dorsal raphe and the relationships with their surrounding tissues. Acta Anatom. Sin. 38, 410–413.
Liu, P. F., Fang, H. Z., Yang, Y., Zhang, Q. Q., Zhou, Q. Q., Chen, S. S., et al. (2017). Activation of P2X3 receptors in the cerebrospinal fluid-contacting nucleus neurons reduces formalin-induced pain behavior via PAG in a rat model. Neuroscience 358, 93–102. doi: 10.1016/j.neuroscience.2017.06.036
Lopez Hill, X., Pascovich, C., Urbanavicius, J., Torterolo, P., and Scorza, M. C. (2013). The median raphe nucleus participates in the depressive-like behavior induced by MCH: differences with the dorsal raphe nucleus. Peptides 50, 96–99. doi: 10.1016/j.peptides.2013.10.002
Lowery-Gionta, E. G., Diberto, J., Mazzone, C. M., and Kash, T. L. (2018). GABA neurons of the ventral periaqueductal gray area modulate behaviors associated with anxiety and conditioned fear. Brain Struct. Funct. 223, 3787–3799. doi: 10.1007/s00429-018-1724-z
Lu, X., Geng, X., Zhang, L. C., and Zeng, Y. (2008). The methodology for labeling the distal cerebrospinal fluid-contacting neurons in rats. J. Neurosci. Methods 168, 98–103. doi: 10.1016/j.jneumeth.2007.09.033
Lu, X. F., Li, Y. Y., Wang, C. G., Wei, J. Q., Ye, Y., Zhang, L. C., et al. (2011). Substance P in the cerebrospinal fluid-contacting nucleus contributes to morphine physical dependence in rats. Neurosci. Lett. 488, 188–192. doi: 10.1016/j.neulet.2010.11.026
Martin, E. M., Devidze, N., Shelley, D. N., Westberg, L., Fontaine, C., and Pfaff, D. W. (2011). Molecular and neuroanatomical characterization of single neurons in the mouse medullary gigantocellular reticular nucleus. J. Comp. Neurol. 519, 2574–2593. doi: 10.1002/cne.22639
Mclaughlin, I., Dani, J. A., and De, B. M. (2017). The medial habenula and interpeduncular nucleus circuitry is critical in addiction, anxiety, and mood regulation. J. Neurochem. 142(Suppl. 2), 130–143. doi: 10.1111/jnc.14008
Monti, J. M., and Jantos, H. (2018). The effects of local microinjection of selective dopamine D1 and D2 receptor agonists and antagonists into the dorsal raphe nucleus on sleep and wakefulness in the rat. Behav. Brain Res. 339, 11–18. doi: 10.1016/j.bbr.2017.11.006
Morton, G., Nasirova, N., Sparks, D. W., Brodsky, M., Sivakumaran, S., Lambe, E. K., et al. (2018). Chrna5-expressing neurons in the interpeduncular nucleus mediate aversion primed by prior stimulation or nicotine exposure. J. Neurosci. 38, 6900–6920. doi: 10.1523/jneurosci.0023-18.2018
Moruzzi, G., and Magoun, H. W. (1949). Brain stem reticular formation and activation of the EEG. Electroencephalogr. Clin. Neurophysiol. 1, 455–473. doi: 10.1016/0013-4694(49)90219-9
Mravec, B., Bodnar, I., Uhereczky, G., Kvetnansky, R., and Palkovits, M. (2012). Effect of lesions of A5 or A7 noradrenergic cell group or surgical transection of brainstem catecholamine pathways on plasma catecholamine levels in rats injected subcutaneously by formalin. Gen. Physiol. Biophys. 31, 247–254. doi: 10.4149/gpb_2012_029
Nunez, A., Garzon, M., and De-Andres, I. (2010). Relationship between the perifornical hypothalamic area and oral pontine reticular nucleus in the rat. Possible implication of the hypocretinergic projection in the control of rapid eye movement sleep. Eur. J. Neurosci. 24, 2834–2842. doi: 10.1111/j.1460-9568.2006.05159.x
Palmiter, R. D. (2018). The parabrachial nucleus: CGRP neurons function as a general alarm. Trends Neurosci. 41, 280–293. doi: 10.1016/j.tins.2018.03.007
Paxinos, G., and Watson, C. (2007). The Rat Brain in Stereotaxic Coordinates, 6th Edn. Amsterdam: Academic Press.
Pollak Dorocic, I., Furth, D., Xuan, Y., Johansson, Y., Pozzi, L., Silberberg, G., et al. (2014). A whole-brain atlas of inputs to serotonergic neurons of the dorsal and median raphe nuclei. Neuron 83, 663–678. doi: 10.1016/j.neuron.2014.07.002
Ramiro, S., Renea, S., Jim, B., and Mariella, D. B. (2009). Nicotinic receptors in the habenulo-interpeduncular system are necessary for nicotine withdrawal in mice. J. Neurosci. 29, 3014–3018. doi: 10.1523/jneurosci.4934-08.2009
Reynolds, D. V. (1969). Surgery in the rat during electrical analgesia induced by focal brain stimulation. Science 164, 444–445. doi: 10.1126/science.164.3878.444
Shafei, M. N., Nasimi, A., Alaei, H., Pourshanazari, A. A., and Hosseini, M. (2012). Role of cuneiform nucleus in regulation of sympathetic vasomotor tone in rats. Pathophysiology 19, 151–155. doi: 10.1016/j.pathophys.2011.11.001
Silvia, T., Rodrigo-Angulo, M. L., Angel, N. E., and Miguel, G. (2013). Synaptic interactions between perifornical lateral hypothalamic area, locus coeruleus nucleus and the oral pontine reticular nucleus are implicated in the stage succession during sleep-wakefulness cycle. Front. Neurosci. 7:216. doi: 10.3389/fnins.2013.00216
Song, S. Y., Li, Y. H., Bao, C. Y., Li, Y., Yin, P. C., Hong, J., et al. (2019). Stereotaxic coordinates and morphological characterization of a unique nucleus (CSF-contacting nucleus) in rat. Front. Neuroan. 13:47. doi: 10.3389/fnana.2019.00047
Song, S. Y., and Zhang, L. C. (2018). The establishment of a CSF-contacting nucleus “Knockout”, model animal. Front. Neuroan. 12:22. doi: 10.3389/fnana.2018.00022
Stettner, G. M., Lei, Y., Kate, B. H., and Leszek, K. (2013). Evidence that adrenergic ventrolateral medullary cells are activated whereas precerebellar lateral reticular nucleus neurons are suppressed during REM sleep. PLoS ONE 8:e62410. doi: 10.1371/journal.pone.0062410
Topchiy, I., Waxman, J., Radulovacki, M., and Carley, D. W. (2010). Functional topography of respiratory, cardiovascular and pontine-wave responses to glutamate microstimulation of the pedunculopontine tegmentum of the rat. Respir. Physiol. Neurobiol. 173, 64–70. doi: 10.1016/j.resp.2010.06.006
Torralba, M. L., Pilar-Cuellar, F., Bravo, L., Bruzos-Cidon, C., and Berrocoso, E. (2018). Opioid activity in the locus coeruleus is modulated by chronicneuropathic pain. Mol. Neurobiol. 56, 4135–4150. doi: 10.1007/s12035-018-1361-9
Vanini, G., and Baghdoyan, H. A. (2013). Extrasynaptic GABAA receptors in rat pontine reticular formation increase wakefulness. Sleep 36, 337–343. doi: 10.5665/sleep.2444
Wang, X. Y., Yan, W. W., Zhang, X. L., Liu, H., and Zhang, L. C. (2014). ASIC3 in the cerebrospinal fluid-contacting nucleus of brain parenchyma contributes to inflammatory pain in rats. Neurol. Res. 36, 270–275. doi: 10.1179/1743132813y.0000000297
Wu, Y. H., Song, S. Y., Liu, H., Xing, D., Wang, X., Fei, Y., et al. (2015). Role of adrenomedullin in the cerebrospinal fluid-contacting nucleus in the modulation of immobilization stress. Neuropeptides 51, 43–54. doi: 10.1016/j.npep.2015.03.007
Xing, D., Wu, Y., Li, G., Song, S., Liu, Y., Liu, H., et al. (2015). Role of cerebrospinal fluid-contacting nucleus in sodium sensing and sodium appetite. Physiol. Behav. 147, 291–299. doi: 10.1016/j.physbeh.2015.04.034
Xu, L., Bloem, B., Gaszner, B., Roubos, E. W., and Kozicz, T. (2010). Stress-related changes in the activity of cocaine- and amphetamine-regulated transcript and nesfatin neurons in the midbrain non-preganglionic Edinger-Westphal nucleus in the rat. Neuroscience 170, 478–488. doi: 10.1016/j.neuroscience.2010.07.001
Xu, L., Janssen, D., van der Knaap, N., Roubos, E. W., Leshan, R. L., and Myers, M. G. Jr., et al. (2014). Integration of stress and leptin signaling by CART producing neurons in the rodent midbrain centrally projecting Edinger-Westphal nucleus. Front. Neuroan. 8:8. doi: 10.3389/fnana.2014.00008
You, H., Hu, S., Ye, Q. P., Wang, F., Huang, J. P., Yang, Y. L., et al. (2018). Role of 5-HT1A in the nucleus of the solitary tract in the regulation of swallowing activities evoked by electroacupuncture in anesthetized rats. Neurosci. Lett. 687, 308–312. doi: 10.1016/j.neulet.2018.06.018
Zemlan, F. P., and Behbehani, M. M. (1988). Nucleus cuneiformis and pain modulation: anatomy and behavioral pharmacology. Brain Res. 453, 89–102. doi: 10.1016/0006-8993(88)90146-1
Zhang, H. A., Sang, N., Ge, X., Huang, Q., Li, X. L., and Sha, J. (2018). Nesfatin-1 in the dorsal raphe nucleus influences visceral sensitivity via 5-HT neurons in male maternally separated rats. Sci. Rep. 8:9334.
Zhang, L. C., Zeng, Y. M., Ting, J., Cao, J. P., and Wang, M. S. (2003). The distributions and signaling directions of the cerebrospinal fluid contacting neurons in the parenchyma of a rat brain. Brain Res. 989, 1–8. doi: 10.1016/s0006-8993(03)03123-8
Keywords: CSF-contacting nucleus, brainstem, spinal cord, projection, retrograde trace
Citation: Song S-Y, Li Y, Zhai X-M, Li Y-H, Bao C-Y, Shan C-J, Hong J, Cao J-L and Zhang L-C (2020) Connection Input Mapping and 3D Reconstruction of the Brainstem and Spinal Cord Projections to the CSF-Contacting Nucleus. Front. Neural Circuits 14:11. doi: 10.3389/fncir.2020.00011
Received: 18 December 2019; Accepted: 10 March 2020;
Published: 31 March 2020.
Edited by:
Manuel S. Malmierca, University of Salamanca, SpainReviewed by:
Daniel Llano, University of Illinois at Urbana-Champaign, United StatesNatasha N. Kumar, University of New South Wales, Australia
Copyright © 2020 Song, Li, Zhai, Li, Bao, Shan, Hong, Cao and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Li-Cai Zhang, licaizhang001@163.com