{'fraction_non_alphanumeric': 0.04199440074656712, 'fraction_numerical': 0.02737412789405857, 'mean_word_length': 4.768777236605999, 'pattern_counts': {'":': 0, '<': 6, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 7, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Excitatory neurons in paraventricular hypothalamus contributed to the mechanism underlying acupuncture regulating the swallowing function

Paraventricular hypothalamus (PVH) is demonstrated to regulate stress, feeding behaviors, and other related homeostatic processes. However, no direct evidence has been investigated for the role of PVH in swallowing function. Acupuncture therapy at Lianquan (CV23) acupoint has been reported to improve the swallowing function in clinical trials, but its underlying mechanism still needs to be uncovered. Thus, we aimed to explore whether PVH involved the acupuncture mediated regulating swallowing function. Chemogenetics, electromyography (EMG) recording, and immunofluorescence staining methods were combined to demonstrate that neurons in PVH could be activated by electroacupuncture (EA) stimulation at CV23, and this neuronal cluster was represented as excitatory neurons. Furthermore, we mapped both the inputs and outputs of PVH neurons using viral tracing. The neurons in PVH projected with the brain regions, including parabrachial nucleus (PBN) and the solitary tract nucleus (NTS), which both participated in the swallowing process. The EA function regulating the swallowing was attenuated after inhibiting the neurons in PVH in the post stroke dysphagia. In conclusion, this study suggested that EA at CV23 could regulate swallowing function involving the excitatory neurons in PVH.

www.nature.com/scientificreports/ PVH Paraventricular hypothalamus S1 The primary sensory cortex ST36 Zusanli VDB Nucleus of the vertical limb of the diagonal band VLM Ventrolateral medulla Swallowing is a critical process for allowing food and fluid to be ingested safely and efficiently, and thereby maintaining physiological and biochemical function. Dysphagia, which commonly occurs after stroke, is associated with the risk of death, occurrence of pneumonia, poor quality of life, and longer hospital stay 1 . Electroacupuncture (EA) therapy is an important and traditional intervention in dysphagia in China [2][3][4] . The treatment of dysphagia after stroke with acupuncture at Lianquan (CV23) has been documented in ancient Chinese medicine books, and its efficacy has been practiced in both pre-clinical experiments and the clinical trials 3,5,6 . However, the mechanism underlying EA at CV23 acupoint on swallowing function was still elusive. The process of food intake depends on counts of chews and swallows 7 . It is known that solid food is transported to the pharynx actively in parallel to it being crushed by chewing and mixed with saliva in the oral cavity, which has to be mixed for a food bolus to reach the swallowing threshold 8 . Feeding behavior is regulated by neural circuits containing the hypothalamus and hindbrain in the mammals 9 . Paraventricular hypothalamus (PVH), which is located in the ventral diencephalon adjacent to the third ventricle, is essential to food intake 10 . Varela et al. suggested that the central lateral parabrachial nucleus (cLPBN) is highly innervated by fibers from the PVH, and the cLPBN-PVH circuit is necessary to induce feeding behaviors 11 . Recent studies indicated that arcuate nucleus (ARC)-PVH and ventral part of the lateral septal nucleus (LSv)-PVH circuits are also important in regulating the food intake behavior 12,13 . Roman et al. reported that the stimulation of the nucleus of solitary tract nucleus (NTS)-PVH circuit signals induces appetitive behavior 14 . Thus, we hypothesized that neurons in PVH might also regulate the swallowing function.
Our previous studies have shown that EA treatment is effective for swallowing disorders. It suggested that EA at CV23 could increase the motor conduction velocity in hypoglossal nerve, enhance the electromyography (EMG) of mylohyoid muscle, and promote release of the substance P in the local acupoint 15 . Besides that, EA at CV23 could activate the swallowing related neurons in the primary motor cortex (M1), parabrachial nucleus (PBN), NTS, nucleus ambiguous (NA) and ventrolateral medulla (VLM) 2,15 , and increase the expression of brain-derived neurotrophic factor 5 . However, no evidence showed the role of PVH in the acupuncture improving the swallowing function.
In the present study, we suggested that the excitatory neurons in PVH could be selectively activated by EA at CV23. The function of PVH is crucial for swallowing process and EA mediated efficacy in both physiological and pathological conditions. This study provided the fundamental and important insight of the EA at CV23 to treat swallowing related disorders for further clinical application.

Results
Effects of electroacupuncture at CV23 on c-Fos expression in PVH. c-Fos staining was widely used to label the neurons newly activated 16 . Electroacupuncture can induce c-Fos expression in brain regions (Fig. 1A). The results showed that c-Fos positive neurons in M1, the primary sensory cortex (S1), PBN and NTS increased by the EA at CV23 (Fig. 1B), which was consistent with previous studies suggesting these brain regions related to the swallowing process [17][18][19] . Besides that, we found that EA at CV23 induced the c-Fos expression in PVH (Fig. 1C). Furthermore, low c-Fos expression were observed in caudate putamen (striatum) (CPu), nucleus of the vertical limb of the diagonal band (VDB), medial lemniscus (ml) and pontine reticular nucleus, ventral part (Pnv), which has been reported to be few related with the swallowing function ( Fig. S1A) 20,21 . The density of c-Fos neurons in the PVH in the CV23 group were significantly increased, comparing to that in the Zusanli (ST36) group (Control: 6580 ± 387.8, EA at CV23 group: 15,441 ± 1730, cells/mm 3 , p < 0.001). There was no statistical difference between the control group and the ST36 group (Fig. 1D). These results suggested that PVH specifically responds to stimulation at CV23 acupoint. Compared with the responses to ST36 acupoint linking to the anti-inflammatory, pains and gastrointestinal disturbs 22,23 , PVH was more sensitive to CV23 acupoint stimulation, which was thought to regulate the swallowing function (Fig. 1D). The anatomical studies have demonstrated the presence of glutamatergic neurons within the PVH, as well as occasional GABAergic neurons 24 . Double staining showed that c-Fos neurons in PVH activated by EA was co-labelled with NeuN + , a marker of neurons, while with few GABAergic neurons in the GAD67-GFP mice (Fig. 1E). Thus, our results suggested that the excitatory neurons in the PVH could respond to EA stimulation at CV23.

Inputs and outputs of PVH neurons.
To further explore the other brain regions connecting to excitatory neurons in PVH, we mapped the monosynaptic outputs of PVH neurons using mCherry with CaMKIIα promoter ( Fig. 2A,B). The mCherry + fibers were observed in the NTS and PBN (Fig. 2C), which suggested neurons in NTS and PBN could receive input from the excitatory neurons in PVH. Notably, NTS contains the generator neurons involved in the triggering, shaping, and timing of the sequential or rhythmic swallowing patterns 25 , while PBN participates in mediating the dive reflex, laryngeal adductor control, swallowing function and upper airway tone, which acts as a mediator that relate primarily to laryngeal closure, upper airway tone and swallowing 26 .
We next aimed to explore how the neuronal population in PVH receives inputs from other brain regions. To solve this problem, we injected retrograde tracer CTB-488 into PVH to screen the monosynaptic projection inputs ( Fig. 2A,D www.nature.com/scientificreports/ transition area (Apir), anteroventral thalamic nucleus (AV), central gray of the pons (CGPn), PBN and NTS. The activation in contralateral cingulate cortex was associated with a better motor recovery 27 . LSv is associated with food intake 13 . NTS and PBN is involved in the swallowing related process 17,26 . All together, these results suggested the structural basis for PVH to play a role in function associated with the swallowing process. www.nature.com/scientificreports/ www.nature.com/scientificreports/ ( Fig. 3A-C). To further confirm the efficacy, the sham acupoint adjacent to CV23 region was set (Fig. S2A), and the results showed EA at sham acupoint didn't change the density of c-Fos neurons in the stroke condition (Fig. S2B&C). The c-Fos expression results indicated that neurons in the PVH could be activated by stroke induction, and further activated by EA stimulation at CV23.

Inhibition of PVH neurons prevented the EA mediated efficacy for improving the swallowing function.
To further verify the role of excitatory neuronal cluster in the PVH in swallowing, we expressed an inhibitory (hM4D) designer receptor exclusively activated by designer drugs (DREADD) in excitatory neurons in the PVH was expressed, and this can effectively inhibit the activity of excitatory neurons (Fig. 4A,B). Pharyngeal swallows were associated with transient increases in pharyngeal pressure, contraction of the mylohyoid muscle as recorded by EMG. Treatment of mice with clozapine N-oxide (CNO) significantly attenuated swallowing response (Fig. 4C). As shown in Fig. 4D, focal cortical ischemia significantly attenuated area under the curve (AUC) of EMG responses (Control: 0.1151 ± 0.01043, PSD: 0.06066 ± 0.004932, p < 0.05). These results suggested that focal cortical ischemia induced dysphagia. AUC of EMG responses was also significantly elevated by EA at CV23 (PSD: 0.06066 ± 0.004932, PSD + EA: 0.09912 ± 0.01325, p < 0.05). However, AUC of EMG responses wasn't changed by EA at sham acupoint (Fig. S3A&B). The result indicated that electroacupuncture at CV23 not sham acupoint improved swallowing function. The model groups with chemoinhibition the excitatory neurons in PVH injected with CNO (i.p.) significantly decreased the swallowing function compared to the model groups (PSD: 0.06066 ± 0.004932, PSD + HM4DI + CNO: 0.01521 ± 0.001004, p < 0.001) (Fig. 4E). However, the model groups with chemoinhibition the excitatory neurons in PVH injected with saline (i.p.) didn't affect swallowing function (PSD + HM4DI + CNO: 0.01521 ± 0.001004, PSD + HM4DI + Saline: 0.04962 ± 0.01118, p < 0.01). While the inhibition of PVH didn't elevate swallowing by EA at CV23 (PSD + EA:0.09912 ± 0.01325, PSD + EA + HM4DI + CNO: 0.03645 ± 0.007895, p < 0.01) (Fig. 4F). Overall, these results demonstrated the excitatory neurons in PVH involved the process of the swallowing both in the physiological and pathological condition, and also mediated the efficacy of EA at CV23. www.nature.com/scientificreports/

Discussion
The present study aimed to detect the effect of EA at CV23 on swallowing and investigated the role of excitatory neurons in PVH during this process. Our results showed that the EA at CV23 could activate excitatory neurons in PVH. The neurons in PVH both projected to and received input from many brain regions associated with food intake, including LSv, PBN and NTS [11][12][13][14] , which provides a structural basis for PVH to play a role in swallowing function. Functionally, swallowing related response was attenuated by chemoinhibition of PVH in the PSD model. EA at CV23 was used to treat for kinds of diseases, including acute dental pain 29 , essential tremor 30 , cranial nerve palsy 31 , stroke 32 and swallowing disorders 2, 6, 15, 33 . Acupuncture at CV23 has been evidenced to activate specific brain regions, including NA 15, 33 , VLM 2, 15 , NTS 15 , striatum 34 , central amygdala (CeA) 35   www.nature.com/scientificreports/ brainstem 41 . The swallowing CPG includes two main groups of neurons located within the medulla oblongata: a dorsal swallowing group located within NTS and a ventral swallowing group located in the VLM adjacent to NA 42 . The coordination of swallowing reflex with rhythmic jaw movements could be regulated by the CeA 43 . The striatum has a prominent role in selecting which motor program should be called into action 44 . These brain regions play different roles in the swallowing process. Our results showed that the neurons in PVH could be activated by EA stimulation at CV23, and this neuronal cluster was demonstrated as the excitatory neurons. The PVH was suggested to be composed of three main types of neurons, including magnocellular, parvocellular, and long-projecting neurons, which plays imperative roles in the regulation of energy balance and various endocrinological activities 45 . The arginine vasopressin positive neurons in the PVH participate in the regulation of feeding behaviors 46 . However, the excitatory neuron in this study belongs to which neuronal type above-mentioned is elusive. In the PVH, melanocortin-4 receptor is expressed on most glutamatergic neurons that project to PBN and dorsal motor nucleus of the vagus nerve (DMV) in the brainstem 47,48 . Further studies would explore the specific neuronal type and clarify the role of neuronal type involved in swallowing.
The monosynaptic input and output neurons were screened throughout the whole brain. The results demonstrated that cingulate cortex, AcbSh and LSv showed strong projection to PVH (Fig. 2E). Neuronal variability in the posterior cingulate cortex (PCC) was significantly reduced after stroke 49 . The dominance of activation in contralateral cingulate cortex was associated with a better motor recovery in the early stage after stroke 27 . AcbSh sends substantial projections to at least two brain regions known to play a role in the control of feeding: the lateral hypothalamus and medial ventral pallidum 50 . LSv plays a critical role in emotionality, social behavior and feeding processes, through neural connections with the hippocampus and hypothalamus 51 . Glucagon-like peptide 1 receptors (GLP-1R) are expressed in the LSv, and endogenous GLP-1 affects feeding and motivation for food 52 . The anterior cingulate cortex (ACC) mediates the food foraging-related behaviors 53 , which is a key region underlying neural processing of social decision-making, specifically tending to compete for foraging high predictive reward food 54 . While PVH neurons both project to the PBN, they synoptically engage distinct efferent nodes, the pre-locus coeruleus and cLPBN, sufficient and necessary to control food intake 55 . The hindbrain, including NTS, PBN, and DMV, is an essential region which could integrate the feeding information from the PVH, other brain regions, and peripheral tissue 56 . Swallowing is an indispensable part during food intake. We hypothesized that these brain regions are involved in swallowing, and all these provide the basis for PVH to play a role in swallowing function.
This study provided the compelling evidence on how the excitatory neurons in PVH involved the process of EA stimulation at CV23 for regulating swallowing function. On the one hand, this result emphasized the role of excitatory neurons in PVH during the swallowing process, on the other hand, this represented new insight for understanding the mechanism of acupuncture treatment at CV23, which could guide the further clinical application of acupuncture therapy.

Materials and methods
Animals. C57BL/6 J mice were purchased from the Animal Laboratory Center of Guangzhou University of Chinese Medicine (license No. SCXK (Yue) 2018-0034). GAD67-GFP knock-in mice were from Professor Yongjun Chen 57 . In total, 130 mice have been used in this study as follows. Firstly, 20 mice were used in the preliminary experiments for exploring PSD molding parameters and virus injection sites. Secondly, 27 mice were used in the immunofluorescence, dividing into control group (N = 6), EA at CV23 group (N = 7), EA at ST36 group (N = 3), PSD group (N = 5) and PSD + EA group (N = 6). Thirdly, 8 mice were used in whole-brain mapping. Fourth, 45 mice were used in electromyography recording, divided into control group (N = 3), PSD group (N =9 , PSD + EA group (N = 9), PSD + HM4DI + Saline group (N = 7), PSD + HM4DI + CNO group (N = 8), and PSD + HM4DI + CNO + EA group (N = 9). We used 100 mice, presented results were obtained from 60 mice. 20 mice were excluded due to the death of mice after modeling and remaining 20 mice were used for preliminary experiments. At last, 30 mice were used to supplement the experiment for sham EA experiments, control group (N = 8), PSD group (N = 11) and PSD + Sham EA group (N = 11). Due to death of the 7 mice after modeling, presented results were obtained from 23 mice. The mice were housed in cages (5 mice/cage) with an ambient temperature of 25 ± 2 °C, a 12-h light-dark cycle, and ad libitum access to water and food. Animal care and experimental manipulation were conducted in accordance with the National Institutes of Health Guide for the Care and Use of Laboratory Animals and approved by the Committee for Care and Use of Research Animals of Guangzhou University of Chinese Medicine (approval No. 20170303), March 3, 2017. The study was carried out in compliance with the ARRIVE guidelines. PSD modeling. Mice were anesthetized with intraperitoneally (i.p.) injected Tribromoethanol (1.25%, 125 mg/kg Sigma-Aldrich, Saint Louis, MO, USA). Ten minutes prior to illumination, Rose Bengal solution (1.5%, Sigma-Aldrich, Saint Louis, MO, USA, 100 mg/kg, i.p.) was administered. Mice were fixed on stereotaxic apparatus (RWD Biotechnology, Shenzhen, China). The skin was cut, following the skull was exposed, and a fiber optic cable delivering laser light was placed onto the skull using a stereotaxic frame. The laser light was centered 0.13 mm anterior to the bregma and 1.1 mm lateral from the midline, targeting the right motor cortex (AP: -0.13, ML: 1.1). A laser (wave length: 530 nm; power: 6.5 mw) was used to irradiate an area of ~ 2mm 2 . After 7 min of irradiation, the scalp was sutured, and mice were placed back in the home cage to recover from anesthesia. Swallowing-related muscle functions were evaluated in awake mice. EA treatment. EA was applied at CV23, ST36 and sham acupoint (located at the right side of cervical region) for 15 min, respectively, and the control group did not take any treatment. Anesthesia was induced www.nature.com/scientificreports/ using 4% isoflurane, and maintained with 2% isoflurane using a mask. After routinely sterilized neck skin, we inserted an acupuncture needle to the upper margin of the midline of the mandible, where CV23 acupoint located. Another needle was inserted 2 mm adjacent to CV23. The needling depth of CV23 is 5 mm. In mouse, ST36, located at 3 mm below the fibular head, in the posterolateral part of the knew. Perfusion and tissue sectioning. c-Fos immunoreactivity was investigated in GAD67-GFP and C57 mice following EA treatment. One hour after the treatment, mice were anesthetized with injected (i.p.) Tribromoethanol (125 mg/kg) and perfused transcardially with 0.9% cold saline (Macklin; China) followed by 4% paraformaldehyde (PFA) (mubiotech, China) in PBS (Thermo Fisher Scientific, America). The brains were extracted, postfixed overnight in 4% PFA at 4 °C and cryopreserved in 15% and 30% sucrose in PBS, and then embedded in Tissue-Tek OCT (Sakura, Japan) compound. Each brain (40 μm) was cut in the transverse plane using a freezing microtome (Thermo, Micro International GmbH, Germany). Confocal microscopy. Confocal fluorescence images (Nikon, Japan) were acquired on a Nikon scanning laser microscope using × 20 and × 40 air objectives. Image analysis was performed using either ImageJ (1.52a, National Institutes of Health, America) or NIS-Viewer v.4.5 (Nikon, Japan).

Quantification and analysis.
Expression of c-Fos was manually quantified at various brain regions by an observer blind to the mouse's experimental condition, using a brain atlas for guidance. Area was measured using ImageJ software (1.52a, National Institutes of Health, America). Next, the PVH (median eminence not included) was contoured and the area was measured on all images based on the labelling for c-Fos. The reference volume was determined as the sum of the area multiplied by the distance between sampled sections (40 μm). The c-Fos positive cell density was calculated for each slice by dividing the total sampled cell numbers by the total volume of the region. A schematic diagram of the selected brain regions is shown in Figs. 1, 3, S2.
Electromyography recording in the Mylohyoid Muscle In Vivo. Anesthesia was induced with 4% isoflurane (Sigma, America). Mice were fixed supine on stereotaxic apparatus. A tube (0.85*0.42 mm PE tube) was adjusted for placement under the tongue. A recording electrode was inserted into the mylohyoid muscle. The reference electrode was inserted into the masseter muscle. The EMG activity of the mylohyoid muscle was evoked and acquired with the Spike2 software (CED, Cambridge, UK) when mice were given water (2 μl/s, 10 s) by a micro-injection pump (HARVARD, USA) after waking up from anesthesia. Data was digitized at 5 kHz with a Power 1401 digitizer (CED, UK) and band-pass filtered at 0.1-1 kHz with a 1902 differential AC amplifier (CED, UK). Evaluated parameters of EMG response is the area under the curve (AUC, represented by mV* msec) of integrated 58 . The area under the curve was calculated in a 10 s window from the onset of water delivery in EMG analysis.
Stereotaxic surgery. Mice were anesthetized with injected (i.p.) Tribromoethanol (125 mg/kg) and fixed on stereotaxic apparatus. After incision, the skull was exposed and holes were drilled over the targeted areas at the following coordinates from bregma. Virus or CTB was injected at a rate of 30 nl/min. After 10 min, the needle was retracted and the incision was closed with a sterile suture. For chemogenomic experiments, rAAV-CaMKIIα-hM4D(Gi)-EGFP (300 nl) was bilateral injected in PVH for inhibition. rAAV-CaMKII-mCherry (50 nl) was unilateral injected in PVH for tracing. Mice were given about 3 weeks to allow expression before next manipulation. CTB-488 (50 nl) was unilateral injected in PVH for afferent trace. CTB was injected 5 days prior to sacrifice for retrograde transport. The viruses or CTB-488 were injected according to the following coordinates: PVH (AP: -0.82; ML: ± 0.2; DV: -4.9). Coordinate mapping and injection was performed using stereotaxic apparatus (RWD Biotechnology, Shenzhen, China).

CNO-induced inhibition.
DREADD-related experiments were performed after waiting 3 weeks for viral expression. Mice were briefly anesthetized with isoflurane then fixed supine on stereotaxic apparatus. We