Abstract
The nodose and jugular vagal ganglia supply sensory innervation to the airways and lungs. Jugular vagal airway sensory neurons wire into a brainstem circuit with ascending projections into the submedius thalamic nucleus (SubM) and ventrolateral orbital cortex (VLO), regions known to regulate the endogenous analgesia system. Here we investigate whether the SubM–VLO circuit exerts descending regulation over airway vagal reflexes in male and female rats using a range of neuroanatomical tracing, reflex physiology, and chemogenetic techniques. Anterograde and retrograde neuroanatomical tracing confirmed the connectivity of the SubM and VLO. Laryngeal stimulation in anesthetized rats reduced respiration, a reflex that was potently inhibited by activation of SubM. Conversely, inhibition of SubM potentiated laryngeal reflex responses, while prior lesions of VLO abolished the effects of SubM stimulation. In conscious rats, selective chemogenetic activation of SubM neurons specifically projecting to VLO significantly inhibited respiratory responses evoked by inhalation of the nociceptor stimulant capsaicin. Jugular vagal inputs to SubM via the medullary paratrigeminal nucleus were confirmed using anterograde transsynaptic conditional herpes viral tracing. Respiratory responses evoked by microinjections of capsaicin into the paratrigeminal nucleus were significantly attenuated by SubM stimulation, whereas those evoked via the nucleus of the solitary tract were unaltered. These data suggest that jugular vagal sensory pathways input to a nociceptive thalamocortical circuit capable of regulating jugular sensory processing in the medulla. This circuit organization suggests an intersection between vagal sensory pathways and the endogenous analgesia system, potentially important for understanding vagal sensory processing in health and mechanisms of hypersensitivity in disease.
SIGNIFICANCE STATEMENT Jugular vagal sensory pathways are increasingly recognized for their important role in defensive respiratory responses evoked from the airways. Jugular ganglia neurons wire into a central circuit that is notable for overlapping with somatosensory processing networks in the brain rather than the viscerosensory circuits in receipt of inputs from the nodose vagal ganglia. Here we demonstrate a novel and functionally relevant example of intersection between vagal and somatosensory processing in the brain. The findings of the study offer new insights into interactions between vagal and spinal sensory processing, including the medullary targets of the endogenous analgesia system, and offer new insights into the central processes involved in airway defense in health and disease.
- airway innervation
- jugular vagal ganglia
- respiratory
- submedius thalamic nucleus
- vagus nerve
- ventrolateral orbital cortex
Introduction
The mammalian airways and lungs are innervated by populations of sensory nerve fibers that are mostly vagal in origin. Two different vagal ganglia, named the nodose and jugular ganglia, supply these sensory fibers and when activated by airway stimuli they initiate a range of sensations and motor responses that regulate breathing, airway clearance, and airway patency. Importantly, nodose and jugular vagal ganglia neurons are embryologically distinct, with the former representing the quintessential vagal visceral sensors and the latter resembling more the somatosensory neurons of the dorsal root ganglia (Nassenstein et al., 2010; McGovern et al., 2015a; Kupari et al., 2019; Mazzone et al., 2020). Within the airway tree, a transition exists from a substantial jugular innervation present in the larynx and proximal airways to the predominantly nodose innervation present in the intrapulmonary airways and lung parenchyma. This transition of vagal sensory innervation is suggestive of a shift from somatosensory processes contributing to laryngeal control to more traditional viscerosensory regulation of the lower airways and lungs.
Responses evoked through vagal sensory innervation of the airways are regulated via reasonably well defined reflex circuits in the brainstem and less well defined networks in the higher brain. Importantly, we have previously shown that the brainstem pathways contributing to nodose and jugular vagal sensory-evoked responses are not the same. Nodose sensory neurons innervating the airways, like most other vagal visceral afferents, terminate centrally in the nucleus of the solitary tract (NTS), whereas neurons of the jugular vagal ganglia terminate within the trigeminal nuclei, notably in and around the paratrigeminal nucleus (Pa5) in the caudal medulla (Driessen et al., 2015; McGovern et al., 2015b; Bautista et al., 2019; Farrell et al., 2020; Kim et al., 2020). These termination sites, and the resultant networks in receipt of outputs from NTS and Pa5 neurons, provide further evidence of the dual viscerosensory and somatosensory controls at play in the airways and lungs.
The behavioral responses that accompany airway stimulation involve network activity in subcortical and cortical brain regions that indirectly receive airway vagal sensory inputs. Our human brain imaging studies demonstrate multiple networks for discriminative and affective processing of airway sensory inputs, but additionally have identified higher brain processes capable of modulating airway vagal sensory processing, including modulation through placebo conditioning and heterotopic painful stimuli (conditioned pain modulation; Mazzone et al., 2009; Farrell et al., 2012; Leech et al., 2012, 2013; Abubakar et al., 2020). Similar inhibitory controls have been described modulating noxious spinal somatosensory processing leading to the perception of pain. Indeed, a well described endogenous neural system, referred to as the “endogenous analgesia system” or the “descending antinociceptive/analgesic pathway,” can modulate sensory inputs to second-order neurons at the level of the spinal dorsal horn. The network organization of this system is relatively well described, and the higher-order controls include multiple subcortical and prefrontal pathways that ultimately regulate descending inputs to the spinal dorsal horn to tune spinal sensory processing (Basbaum and Fields, 1984; Dubner and Ren, 2013). One pathway regulating the descending analgesia system involves the activation of the submedius thalamic nucleus (SubM) and its connectivity to the ventrolateral orbital cortex (VLO; Coffield et al., 1992; Yoshida et al., 1992; Zhang et al., 1998, 1999; Tang et al., 2009; Kuramoto et al., 2017). Notably, we have previously shown that airway-specific vagal sensory inputs to the brainstem reach the SubM via connectivity specifically within in the medullary Pa5 (and not the NTS), suggestive that jugular vagal pathways may also intersect with the descending analgesia system (McGovern et al., 2015b). Given these combined observations, we set out in the present study to investigate whether the SubM is capable of orchestrating regulatory controls over laryngeal vagal sensory processing in the brainstem.
Materials and Methods
Animals
Experiments were approved by an accredited institutional Animal Ethics Committee and were conducted on adult Sprague Dawley rats of either sex (weight range, 300–500 g; N = 109). Animals were housed on a 12 h light/dark cycle and given ad libitum access to food and water. No statistical differences were observed between the measured responses of male and female rats in initial experiments, and, in turn, the data were pooled for all subsequent statistical comparisons.
Physiologic studies in anesthetized rats
Rats were anesthetized with urethane (1.5 g/kg, i.p.), the level of which was determined to be adequate by assessing the palpebral and limb withdrawal reflexes. Animals were placed supine on a thermostatically controlled heating pad set at 36.5°C, and body temperature was maintained via feedback from a rectal probe. A midline incision was made along the ventral surface of the neck to expose the larynx, trachea and underlying nerves and blood vessels. The left carotid artery was cannulated using polyethylene tubing (internal diameter, 0.5 mm; outer diameter, 0.9 mm) attached to a pressure transducer filled with heparinized saline (50 U/ml; Sigma-Aldrich) to measure arterial blood pressure (ABP). The distal extrathoracic trachea was cannulated and connected via a side port to a pressure transducer to measure the tracheal pressure and changes associated with spontaneous respiration. Output from pressure transducers was digitized and recorded (LabScribe, iWorx). A midline incision was made through the larynx, and a platinum bipolar stimulating electrode was placed on the mucosal surface (Fig. 1A). After a 20 min stabilization period, electrical stimulation (S48 square pulse stimulator, Grass Instruments) of the larynx was performed to activate the peripheral terminals of vagal afferents innervating this site (Driessen et al., 2015, 2020). A voltage response curve was initially conducted by delivering increasing voltages (0.5–5 V; 0.5 V increments at 2 min intervals) to the larynx, at a constant maximum stimulating frequency (32 Hz), pulse duration (1 ms), and train duration (10 s). In each preparation, the optimum voltage to elicit complete apneic responses, mimicking the laryngeal chemoreflex evoked by aspiration or chemical activation of vagal laryngeal afferents (Xia et al., 2008), was used to conduct a frequency response analysis by sequentially increasing the stimulation frequency (1–32 Hz). Across all experiments, the average voltage used was 1.85 ± 0.12 V (range, 1–3 V), which is expected to activate all subtypes of sensory terminals innervating the larynx (McAllen et al., 2018). The specificity of vagal sensory pathways mediating respiratory effects was confirmed after unilateral and bilateral vagotomy. Following this initial baseline electrical stimulation of the larynx, rats were placed onto a stereotaxic frame (model 940, David Kopf Instruments), the skull exposed, and burr holes drilled at the corresponding stereotaxic coordinates for bilateral microinjections into the SubM [anteroposterior (AP), −2.75 mm; mediolateral (ML), ±2.00 mm (12° angle); dorsoventral (DV), −7.10 mm, relative to bregma] of vehicle (0.9% saline; 250 nl), 5-HT (5 µg of 5-HT-HCl, 150 nl; Sigma-Aldrich) to stimulate analgesic promoting SubM neurons via 5-HT1A and 5-HT2 receptors, or the generic neuronal inhibitor muscimol (10 ng, 150 nl; Sigma-Aldrich) to inhibit SubM-mediated analgesic neurons (Xiao et al., 1999, 2005a,b; Jia et al., 2004). Microinjections of pharmacological agents were performed using glass-pulled micropipettes (tip diameter, ∼20 µm) connected to a microprocessor-controlled picopump (model PV820, World Precision Instruments). Solutions were combined with red fluorescent beads (Lumofluor) to allow easy visualization of injection sites. Injections outside of the SubM target site, into adjacent nuclei, served as control experiments to ensure that evoked responses were mediated by the SubM. In some experiments, an electrolytic lesion of the ventrolateral orbital cortex (AP, +3.50 mm; ML, ±2.50 mm; DV, −5.10 mm) was performed using an LMD Lesion Making Device (15 s, 1 mA, DC constant current; Ugo Basile), with sham lesions involving the insertion of the electrode but without delivering any current. After microinjections, rats were immediately returned to a supine position and the frequency response to electrical stimulation of the larynx reassessed in a paired experimental design. In all instances, the second series of laryngeal stimuli was completed within a period of 30 min from the first microinjection.
In some experiments, bilateral microinjections of the nociceptor stimulant capsaicin were performed into the medullary NTS and Pa5, instead of electrical stimulations of the larynx, to activate the central terminals of vagal afferent neurons projecting to these discrete sites (McGovern et al., 2012, 2015a,b; Driessen et al., 2015; Mazzone and Undem, 2016; Kim et al., 2020). In these experiments, rats were set up to record respiration and blood pressure parameters, following which they were placed onto a stereotaxic frame and either 5-HT or vehicle was bilaterally microinjected into the SubM (as above). To access the medulla, the head was angled at 45° and a midline incision was made through the skin of the animal, the posterior neck muscles, and the dura mater between the occipital bone and C1 vertebra. Using obex as a reference point (∼0.5 mm rostral to the calamus scriptorius), microinjections of capsaicin (0.5 μm, 150 nl; Sigma-Aldrich) were made into either the caudal NTS (0 mm level to obex, 0.3 mm lateral to the midline, and 0.5 mm below the medullary surface) or Pa5 (0.5 mm rostral to obex, 2–3 mm lateral to the midline, and 0.5 mm below the medullary surface). Respiration and blood pressure parameters were monitored for 30 min after capsaicin microinjection.
Neuronal tracer and viral vector injections
Rats were anesthetized with isoflurane (4% induction and 2% maintenance in oxygen) via a nose cone and placed onto a stereotaxic frame as described above. The adequacy of anesthesia and body temperature was monitored as described above. Microinjections were made using the same stereotaxic coordinates described above. SubM–VLO connectivity was assessed by retrograde tracing from the VLO using green fluorescent retrobeads (N = 4; 250 nl; Lumofluor) or anterograde tracing from the SubM using biotin dextran amine (BDA; 10% w/v, 150 nl; N = 4; Sigma-Aldrich). For conscious behavioral experiments, we opted to target specifically the VLO-projecting SubM neurons by bilaterally microinjecting adeno-associated virus (AAV) retrograde-pkg-Cre (250 nl; catalog #24593, Addgene) into the VLO and either AAV8-hSyn-DiO-hM3D(Gq)-mCherry (N = 10; 150 nl; catalog #44361, Addgene) or control vector AAV8-hSyn-DiO-mCherry (N = 9; 150 nl; catalog #50459, Addgene) bilaterally into the SubM. For anterograde transsynaptic neuronal tracing, 15 µl of the conditional HSV1-H129floxed (1 × 109 pfu/ml) was injected into the airway lumen of rats (N = 6) that had previously received bilateral microinjections of the AAV2 pkg-Cre (250 nl; 1 × 1011 viral genomes/ml) into the Pa5 (see above for coordinates; see Fig. 4E). Following the completion of microinjections, sterile bone wax was applied to cover the holes in the top of the skull, incisions were sutured, and rats received subcutaneous injections of meloxicam (1 mg/kg) and 1 ml of 0.9% sterile saline. Rats were allowed to recover for 7 d for conventional tracing studies, 5–6 d for transsynaptic viral tracing studies (McGovern et al., 2015a,b), and 28 d for conscious respiratory behavioral studies.
Conscious respiratory studies
Conscious unrestrained rats were exposed to an aerosolized respiratory irritant in a four-site Buxco whole-body plethysmography system. The experimental design consisted of a 29 d protocol whereby rats were microinjected with either the hM3D(Gq) or control vector (see above; day 1), were allowed to recover for 23 d, then acclimatized to the plethysmography chambers (1 h/d) without aerosolized challenge (days 23, 25, and 27). On challenge day (day 29), all rats were administered the hM3D(Gq) ligand clozapine N-oxide (CNO; 3 mg/kg, i.p.), following which they were placed in the chambers for 30 min before receiving an initial aerosolization of saline (vehicle) followed by increasing challenge doses of capsaicin (0.3, 3, 30, 100, and 300 μm). Challenge substances were nebulized (Buxco Control-10 Aerosol Delivery Unit 10 LPM, Aerogen nebulizer unit AG-AL1000 with filter cap of 3.1 µm particle size) as 1 ml solutions over 5 min, at an air flow velocity of 5 L/min, and at 60% nebulizer duty cycle. Challenge substances were nebulized for 5 min, followed by a 5 min postchallenge response period. A further 5 min was allowed between each challenge dose to allow respiratory parameters and behaviors to return to baseline. Rats were left in the chambers for 30 min after the completion of the capsaicin challenges to allow for adequate time for induction of c-Fos post CNO. Two additional control groups of rats, no surgery + CNO (N = 6) and no surgery and no CNO (N = 6), underwent the same respiratory challenge protocol.
Tissue harvest and histologic analyses
Following completion of conscious, physiological, or tracing studies, rats were overdosed with sodium pentobarbital (100 mg/kg, i.p.) and transcardially perfused with 200 ml of 0.1 m PBS, pH 7.4, and 200 ml of 4% paraformaldehyde (PFA) in 0.1 m PBS. Brains were dissected and postfixed overnight in 4% PFA before being cryoprotected in 25% sucrose at 4°C. Tissues were frozen in OCT (optimal cutting temperature) embedding compound, and 50 µm cryostat cut sections encompassing the areas of interest were collected serially into 0.1 m PBS. Sections from the physiological studies were mounted directly onto gelatin-coated slides and coverslipped to assess SubM injection sites and VLO lesions. All fluorescent immunohistochemistry was performed at room temperature. Brain sections were rinsed then blocked in 10% donkey serum in 0.1 m PBS for 1 h before incubation for 48 h in the primary antibody of interest: rabbit anti-calbindin D-28k (1:1000; catalog #CB-38a, Swant); mouse anti-NeuN (1:1000; catalog #MAB377, Merck Millipore); rabbit anti-c-Fos (1:2000; catalog #26209, Immunostar). All antibodies were diluted in antibody diluent (2% donkey serum and 0.3% Triton X-100 in 0.1 m PBS). After the required incubation, sections were washed three times with 0.1 m PBS for 20 min each then incubated for 1 h in appropriate fluorescently conjugated secondary antibodies (1:500; Thermo Fisher Scientific). BDA traced sections were incubated in streptavidin-Alexa Fluor 594 conjugate (1:500; catalog #S32356, Thermo Fisher Scientific) for 1 h. Immunostained sections were kept in order, mounted onto gelatin-coated slides, and coverslipped with an antifade mounting media (Fluoroshield, catalog #F6182, Sigma-Aldrich). All tissues were visualized under the appropriate filter cube using a Leica DM6B LED Microscope, and images were captured using a Leica DFC7000T camera. Representative photomicrographs were assembled in Adobe Photoshop CS6.
Experimental design and statistical analyses
All statistical comparisons were made in Prism software (version 7.03; GraphPad Software). For Student's t test and ANOVA all T and p, and F and p values for the interaction effects and individual treatment effects, respectively, are reported either in the text or in corresponding figure legends. Degrees of freedom for each F value are designated by subscripts as per convention. Data are reported as the mean ± SD, unless stated otherwise; and differences were considered significant at p < 0.05.
Anesthetized physiological respiratory studies
Respiratory rate, blood pressure, and heart rate were calculated from the raw chart recordings over a 10 s period, 2 min before any stimulation (baseline) and then again during each laryngeal stimulation (evoked response), and multiplied by 6 to determine the equivalent per minute values, and respiration was plotted as the mean ± SD for all groups. For capsaicin medullary microinjection experiments, respiratory rate, blood pressure, and heart rate were collected from the raw chart recordings at 60 s intervals over a period of 10 min after the first microinjection of capsaicin. Baseline measurements were determined as the 60 s period before capsaicin microinjection with data presented as the percentage change relative to baseline and plotted as the mean ± SD. Peak blood pressure was normalized to baseline levels and plotted as the mean ± SD for all groups. Normalization was performed to remove confounding factors such as interanimal variation in baseline parameters. A D'Agostino–Pearson normality test was used, then stimulus response curve data were compared using a two-way ANOVA for repeated-measures followed by Bonferroni's multiple-comparisons test, while a Student's two-tailed paired/unpaired t test was performed to compare the maximal reduction in respiratory rate (Emax) and the electrical frequency required to elicit 50% of the maximum response (EF50).
Conscious respiratory studies
Respiratory frequency (breaths/min), tidal volume (milliliters per kilogram), minute volume (milliliters per minute), inspiratory and expiratory time (seconds), peak inspiratory and expiratory flow (milliliters per minute), and enhanced pause (Penh, a measure of airflow obstruction) were automatically sampled from the calibrated Buxco flow trace, corrected for animal body weight and ambient chamber temperature, and calculated in real time by the plethysmography software (Buxco Finepointe Software, version 2.1.0.9). Real-time data collected every 2 s were used to calculate the average per 5 min for each cohort at baseline (5 min before vehicle nebulization) and during each nebulized challenge (5 min nebulization time). Data were normalized to vehicle nebulization values for each rat and then plotted as the mean ± SD for all groups. Normalization was performed to remove confounding factors such as interanimal variation in baseline parameters. Data were first tested for normality using a D'Agostino–Pearson test, and the response to each dose of capsaicin was compared between groups using a two-way ANOVA for repeated measures followed by Bonferroni's multiple-comparisons test.
Neuroanatomical data
The efficacy of AAV transduction of the hM3D(Gq) and control vector were first quantified after conscious physiological studies by manually counting the number of mCherry-positive neurons within the entire rostrocaudal extent of the SubM. Counts from both left and right sides were pooled across the two groups and reported as the mean ± SD at every 100 µm. Student's two-tailed unpaired t test was used to compare the hM3D(Gq) and control AAV groups. In addition, tile-scan images at three representative rostrocaudal levels (distance from bregma: −2.20, −2.80, and −3.30 mm) were acquired from each hM3D(Gq)-injected rat with every mCherry-positive neuron (using the region of interest tool on the Leica X acquisition software) transposed to traces of the SubM to demonstrate the spatial specificity of the injection sites (left and right sides were kept separate). To quantify neuronal activation in transduced SubM neurons following CNO administration, the number of neurons expressing mCherry and/or c-Fos were counted at every 200 µm in both hM3D(Gq)-injected and control AAV-injected rats. These counts were presented as a percentage of the total number of mCherry and/or c-Fos neurons in each animal.
A similar approach was used for the anterograde transsynaptic tracing experiments with all H129 tdT-positive neurons counted within the entire rostrocaudal extent of the SubM and reported as the mean ± SD at every 100 µm. To demonstrate spatial specificity, traces of individual H129 tdT-positive SubM neurons were obtained from the peak level of infection (distance from bregma, −2.70 mm).
Results
Pharmacological modulation of the SubM regulates laryngeal-evoked respiratory reflexes via connectivity to the VLO in anesthetized rats
At optimum stimulus intensities, increasing stimulus frequency applied to the laryngeal mucosa evoked a frequency-dependent reduction in respiratory rate as well as an increase in mean ABP in urethane-anesthetized rats (Fig. 1B,D,G), while no effect was seen on heart rate [462.6 ± 36.0 beats per minute (bpm) at baseline compared with 473.1 ± 34.0 bpm at the highest stimulus frequency tested]. It must be noted that baseline heart rate was elevated compared with nonanesthetized rats, perhaps reflective of elevated sympathetic tone commonly observed with urethane anesthesia, and this may have masked any observable change in heart rate during stimulation. Evoked responses during electrical stimulation were substantially reduced following unilateral vagotomy (Fig. 1C′) and were absent following bilateral vagotomy confirming the specificity of the stimulus for vagal afferent pathways (Fig. 1C′′). Bilateral microinjection of vehicle (0.9% saline) into the SubM (N = 9) did not alter the stimulation frequency-dependent fall in respiratory rate (EF50 and Emax: paired, t(8) = 1.03, p = 0.332, and t(8) = 0.595, p = 0.569, respectively; Fig. 1D), or blood pressure (EF50 and Emax: paired, t(9) = 0.4771, p = 0.6460, and t(9) = 0.0423, p = 0.9674, respectively; Fig. 1G). Bilateral microinjection of 5-HT into the SubM to excite neurons (N = 9) had no effect on basal respiratory rate and blood pressure, but significantly attenuated the respiratory response (EF50 and Emax: paired, t(8) = 10.00, p < 0.0001, and t(8) = 8.13, p < 0.0001, respectively; Fig. 1E) and modestly reduced the blood pressure response associated with laryngeal stimulation (EF50 and Emax: paired, t(8) = 0.7242, p = 0.4962, and t(8) = 2.329, p = 0.0587, respectively; Fig. 1H). Bilateral microinjection of muscimol to inhibit neurons into the SubM (N = 6) also had no effect on basal respiratory rate and blood pressure and did not modify blood pressure responses evoked by laryngeal stimulation (EF50 and Emax: paired, t(6) = 0.2354, p = 1.000, and t(6) = 0.5410, p = 0.6261, respectively; Fig. 1I). Strikingly, muscimol had the opposite effect of 5-HT, significantly augmenting the respiratory response evoked by laryngeal stimulation (EF50 and Emax: paired, t(5) = 3.02, p = 0.029, and t(5) = 2.71, p = 0.042, respectively; Fig. 1F). Post hoc histologic assessment of the injection sites showed that 9 of 14 rats had successfully received bilateral SubM microinjection “hits,” whereas injections in the remaining rats were restricted to unilateral hits. Nevertheless, this variance did not appear to have a significant impact on response magnitudes evoked by 5-HT or muscimol (data not shown). In addition, to determine that the physiological responses observed were SubM specific, six rats received microinjections of 5-HT into brain regions surrounding the SubM (Fig. 2). This resulted in no significant change to respiratory rate (Fig. 2A–F), blood pressure, and heart rate (data not shown) after 5-HT microinjection.
Pharmacological modulation of the submedius thalamic nucleus regulates laryngeal-evoked respiratory reflexes. A, Anesthetized respiratory physiology setup and example of a SubM microinjection site. B, Electrical stimulation of the larynx results in a frequency-dependent slowing of respiration accompanied by hypertension. C–C′′, This respiratory effect (C) is partially reduced by unilateral vagotomy (C′) and completely abolished following bilateral vagotomy (C′′). D–F, Stimulus frequency-dependent responses are assessed before (Before) bilateral microinjection of pharmacological agents into the SubM and repeated afterward (After). D, Stimulation-evoked changes in respiratory rate before and after microinjection of vehicle (N = 9) into the SubM. E, Stimulation-evoked changes in respiratory rate before and after microinjection of 5-HT (N = 9) into the SubM. F, Stimulation-evoked changes in respiratory rate before and after microinjection of muscimol (N = 6) into the SubM. G–I, Blood pressure stimulus–response curves before and after microinjection of vehicle (N = 9; G), microinjection of 5-HT (N = 9; H), and microinjection of muscimol (N = 6) into the SubM (I). *p < 0.05, ****p < 0.0001, two-tailed paired t test (F–I). Scale bar, A, 300 µm. mt, Mamillary tract. Gray box (B, C) denotes a 10 s period of maximum electrical stimulus intensity.
Schematic representations of microinjections of 5-HT into brain regions surrounding the SubM. Each coronal schematic represents an individual rat, with the SubM highlighted in yellow shading and the “missed” injection sites highlighted in gray shading. A–F, Microinjection sites encompassed rat 1, ventrolateral and paracentral thalamic nuclei (A); rat 2, central medial thalamic nucleus (B); rat 3, central medial, paracentral, and ventral reuniens thalamic nuclei (C); rat 4, reuniens and ventromedial thalamic nuclei (D); rat 5, dorsal hypothalamic area (E); and rat 6, rhomboid and ventromedial thalamic nuclei (F). No significant effect on respiration was observed for each brain region (see graph insert for each brain image). Each coronal slice represents bregma level from −2.70 to −2.80 mm.
Microinjection of the anterograde tracer BDA into the SubM location responsive to 5-HT and muscimol above resulted in dense SubM neuron terminal field labeling in the VLO, spanning from bregma +3.20 to +4.00 mm. Many of the labeled terminals were concentrated in cortical layer 4, and less so in layer 5 (Fig. 3A). Furthermore, retrograde tracing from this VLO site using green fluorescent retrobeads identified a large population of neurons in the SubM (Fig. 3B), confirming the neuroanatomical connection. In functional experiments (Fig. 3C), bilateral sham lesions (N = 3) or electrolytic lesions of the VLO (N = 5) followed by bilateral microinjection of vehicle into the SubM had no effect on basal respiratory and blood pressure parameters (Fig. 3D,F). However, VLO lesions significantly potentiated the Emax evoked by laryngeal stimulation, when compared with the sham lesion rats (EF50 and Emax: unpaired, t(6) = 0.562, p = 0.594, and t(6) = 4.875, p = 0.0028, respectively; Fig. 3D). Although the potentiating effect was not as substantial as that seen following SubM muscimol microinjection (Fig. 1F), the data are nonetheless consistent with the VLO lesion unmasking a tonic inhibitory influence of the SubM–VLO circuit on respiratory reflex responses evoked by laryngeal stimulation. No effect of VLO lesion on the laryngeal-evoked blood pressure response was observed in this experiment (EF50 and Emax: unpaired, t(8) = 0.2471, p = 0.8120, and t(8) = 1.132, p = 0.2947, respectively; Fig. 3F). Bilateral microinjection of 5-HT into the SubM of sham (N = 4) or VLO lesioned (N = 6) rats did not alter basal respiratory or blood pressure parameters but characteristically reduced the magnitude of laryngeal-evoked respiratory slowing in sham-lesioned rats (Fig. 3E,G). Strikingly, this effect of 5-HT was blocked entirely in rats receiving bilateral VLO lesions (EF50 and Emax: unpaired, t(8) = 9.061, p < 0.0001, and t(8) = 7.369, p < 0.0001, respectively; Fig. 3E). A small difference was also seen in the blood pressure response associated with laryngeal stimulation between 5-HT-treated sham-lesioned and VLO-lesioned rats, testing significant for the Emax response (EF50 and Emax: unpaired, t(8) = 0.2939, p = 0.7787, and t(8) = 2.651, p = 0.0380, respectively; Fig. 3G). This finding is entirely consistent with outputs from the SubM to the VLO playing a central role in SubM-mediated regulation of laryngeal respiratory reflexes.
The submedius thalamic nucleus regulates laryngeal-evoked respiratory reflexes in anesthetized rats via connectivity to the ventrolateral orbital cortex. A, B, The SubM has dense projections to the VLO, as revealed by anterograde (A) and retrograde (B) neuroanatomical tracing. C, Schematic summary of the physiological experiment design to test the effect of bilateral electrolytic lesioning of the VLO (insert photomicrograph demonstrates example lesion site) or sham lesions on SubM-mediated modulation of laryngeal reflexes. D, Laryngeal stimulation-evoked changes in respiratory rate for lesion VLO (N = 5) and sham VLO (N = 3) plus microinjection of vehicle into the SubM. E, Laryngeal stimulation-evoked changes in respiratory rate for lesion VLO (N = 6) and sham VLO (N = 4) plus microinjection of 5-HT into the SubM. F, Blood pressure stimulus–response curves for lesion VLO (N = 5) and sham VLO (N = 3) plus microinjection of vehicle into the SubM. G, Blood pressure response curves for lesion VLO (N = 6) and sham VLO (N = 4) plus microinjection of 5-HT into the SubM. *p < 0.05, **p < 0.01, ****p < 0.0001, two-tailed unpaired t test (D–G). AI, Agranular insula; CI, claustrum; mt, mamillary tract; rf, rhinal fissure; VM, ventromedial thalamus.
The SubM–VLO circuit regulates nocifensive respiratory reflexes in conscious rats
Functional targeting of the SubM–VLO circuit in conscious rats was achieved using DREADD (designer receptor exclusively activated by designer drugs) technology introduced using a combinatorial viral vector approach. In conscious, freely behaving rats undergoing whole-body plethysmography, no change in baseline respiratory parameters were observed 30 min following administration of CNO between any of the experimental groups (Table 1). In naive rats (i.e., not treated with AAV vectors or CNO; N = 6), increasing concentrations of aerosolized capsaicin produced a dose-dependent increase in respiratory rate and tidal volume (respiratory rate, repeated-measures two-way ANOVA: interaction between group and dose, F(4) = 0.21, p = 0.9338; dose effect, F(4) = 3.19, p = 0.0207; group effect, F(1) = 0.51, p = 0.4770; Fig. 4A; tidal volume, repeated-measures two-way ANOVA: interaction between group and dose, F(4) = 0.91, p = 0.4656; dose effect, F(4) = 4.83, p = 0.0023; group effect, F(1) = 2.48, p = 0.1216; Fig. 4D), the former largely driven by a progressive, albeit variable, reduction in inspiratory time (repeated-measures two-way ANOVA: interaction between group and dose, F(4) = 0.72, p = 0.5792; dose effect, F(4) = 3.57, p = 0.0123; group effect, F(1) = 16.60, p = 0.0002; Fig. 4B). No effect was observed on expiratory time (repeated-measures two-way ANOVA: interaction between group and dose, F(4) = 0.50, p = 0.7354; dose effect, F(4) = 0.50, p = 0.7351; group effect, F(1) = 1.52, p = 0.2231; Fig. 4C). Importantly, CNO did not significantly modify any respiratory parameter response evoked by any concentration of inhaled capsaicin in control experiments involving naive rats (Fig. 4, Table 2), arguing that CNO has no observable off-target biological activity in this experimental model. Capsaicin inhalation in CNO-treated control AAV rats (N = 9) also produced respiratory responses that were largely indistinguishable from those seen in naive control rats, characterized by a dose-dependent increase in respiratory rate, tidal volume, and a progressive shortening of inspiratory time (Fig. 5C–E, Table 2).
Baseline respiratory measurements after CNO administration in naive (N = 6), naive + 3 mg/kg CNO (N = 6), control AAV + 3 mg/kg CNO (N = 9), and hM3D(Gq) + 3 mg/kg CNO rats (N = 10) rats
Respiratory parameters measured during conscious capsaicin challenge in in naive (N = 6), naive + 3 mg/kg CNO (N = 6), control AAV + 3 mg/kg CNO (N = 9), and hM3D(Gq) + 3 mg/kg CNO rats (N = 10) rats
Respiratory parameters measured during conscious capsaicin challenge are unaltered by CNO in the absence of hM3D(Gq) DREADD expression. A–D, Respiratory rate (RR; A), inspiratory time (Ti; B), expiratory time (Te; C), and tidal volume (Tv; D) in naive (N = 6) or naive + 3 mg/kg CNO (N = 6) rats.
Activation of the SubM–VLO circuit in conscious unrestrained rats regulates nocifensive respiratory reflexes induced by inhalation of capsaicin. A, Microinjection of AAVrg pkg-Cre into the VLO to achieve Cre-dependent expression of either hM3D(Gq) mCherry or control mCherry in SubM–VLO projecting neurons only. B, Schematic outlining the conscious respiratory challenge protocol. C–F, Respiratory parameters measured during conscious capsaicin challenge in hM3D(Gq) + CNO (N = 10) or control AAV + CNO (N = 9) rats: respiratory rate (RR; C), inspiratory time (Ti; D), expiratory time (Te; E), Tidal volume (Tv; F). G, H, Representative images of c-Fos (green) expression in SubM neurons expressing either hM3D(Gq) mCherry (G) or control AAV mCherry (H) following administration of CNO 2 h before rats were killed. Arrowhead indicates colabeled neuron in box insert. I, Quantification of the number of SubM neurons expressing mCherry and/or c-Fos represented as a percentage of the total number of mCherry and/or c-Fos neurons counted [hM3D(Gq), 529 neurons; control AAV, 580 neurons]. *p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001, hM3D(Gq) versus control AAV by Bonferroni's multiple comparisons (C–F). Scale bars: G, H, 75 µm.
By contrast, CNO administration to rats expressing the hM3D(Gq) in VLO-projecting SubM neurons (N = 10) significantly blunted the dose-dependent capsaicin-evoked increases in respiratory rate (repeated-measures two-way ANOVA: interaction between group and dose, F(4) = 0.24, p = 0.917; dose effect, F(4) = 25.02, p < 0.0001; group effect, F(1) = 40.72, p < 0.0001; Fig. 5C) and the accompanying reduction in inspiratory time (repeated-measures two-way ANOVA: interaction between group and dose, F(4) = 3.58, p = 0.011; dose effect, F(4) = 24.08, p < 0.0001; group effect, F(1) = 19.20, p = 0.0005; Fig. 5D). Notably, these rats also displayed a significant increase in expiratory time (repeated-measures two-way ANOVA: interaction between group and dose, F(4) = 0.55, p = 0.696; dose effect, F(4) = 5.15, p = 0.0011; group effect, F(1) = 20.99, p = 0.0003; Fig. 5E), which resulted in overall respiratory rate falling significantly below baseline levels at capsaicin challenge doses of 0.3–30 μm (Fig. 5C). The consistent increase in tidal volume evoked by capsaicin was largely unaffected by CNO in hM3D(Gq) rats (repeated-measures two-way ANOVA: interaction between group and dose, F(4) = 2.30, p = 0.0674; dose effect, F(4) = 35.38, p < 0.0001; group effect, F(1) = 0.76, p = 0.396; Fig. 5F), suggesting that breathing rate and depth are independently regulated by capsaicin inhalation. Overall, the substantial increase in minute ventilation evoked by capsaicin in naive and control AAV-transfected animals was abrogated by the activation of hM3D(Gq) expressed by VLO-projecting SubM neurons and converted to a hypoventilatory response at lower capsaicin challenge concentrations (Table 2). As expected, CNO treatment in these experiments induced the expression of c-Fos in the majority of SubM–VLO-projecting neurons virally transfected with hM3D(Gq) with only a very small number of transfected neurons expressing c-Fos in the control AAV rats (Fig. 5G–I), confirming the specificity of CNO for hM3D(Gq)-transduced neurons.
The combination of AAVrg pkg-Cre injected into the VLO and AAV8 hSyn-DiO-hM3D(Gq)-mCherry (or control vector) injected into the SubM resulted in robust expression of vector reporter specifically in the SubM neurons with output connectivity to the VLO (Fig. 6A). Transduced neurons spanned a rostrocaudal distribution within the SubM, which ranged from −2.00 to −3.40 mm relative to bregma, with the largest density of neurons observed at approximately −2.70 to −2.90 mm. More often these neurons were located within the dorsal half of the nucleus and especially neighboring the dorsal boundary (Fig. 6B–D). No differences in the spread or positioning of transduced neurons were seen between the hM3D(Gq) and control AAV groups (Fig. 6C). We found evidence for transduced neurons within other brain regions located in central and medial thalamus; however, the number of transduced neurons in each region was significantly lower than that observed in the SubM of hM3D(Gq) rats (one-way ANOVA: F(8) = 15.87, p < 0.0001; Fig. 6A).
Neuroanatomical organization of neurons in the submedius thalamic nucleus involved in modulating vagal reflexes. A, Total number of transduced neurons in the brain of hM3D(Gq) mCherry rats (N = 10). Data are presented as the mean ± SD for each brain region. B, Schematic representation of the distribution of VLO-projecting hM3D(Gq)-transduced neurons within the SubM. Each coronal slice summarizes the location of individual hM3D(Gq) mCherry neurons (red dots, from all hM3D(Gq) rats, N = 10) at three different bregma levels spanning the entire rostrocaudal extent of the SubM. C, Quantification of the rostrocaudal spread (at 100 µm intervals) of neurons transduced by hM3D(Gq) and control AAV (N = 9 rats). Data are represented as the mean ± SD. D–E, Representative images of hM3D(Gq) mCherry-transduced neurons (D) with examples of high-power individual neurons (D′) and control AAV mCherry-transduced neurons (E) in the SubM. F, Schematic representation of the anterograde transsynaptic HSV1-H129floxed viral tracer (permanently switches to tdTomato (tdT) fluorescence in the presence of Cre-recombinase; McGovern et al., 2015a) injected into the airways of rats that had previously been microinjected with the AAV pkg Cre construct into the medullary Pa5. G, Rostrocaudal spread of SubM H129 tdT-infected neurons in receipt of Pa5 projecting airway neurons at 100 µm intervals (N = 6). Data are represented as the mean ± SD. Coronal slice summarizes the location of individual tdT neurons (red dots, from all rats) at the peak bregma level of H129 tdT infection. H, Representative image of H129 tdT-infected neurons in the SubM. ****p < 0.0001, SubM, determined by Bonferroni's multiple comparisons (A). VM, Ventromedial thalamic nucleus; CM, central medial thalamic nucleus; PC, paracentral thalamic nucleus; MD, mediodorsal thalamic nucleus; Rh, rhomboid thalamic nucleus; AM, anteromedial thalamic nucleus; IAM, interanteromedial thalamic nucleus; AMV, anteromedial ventral thalamic nucleus; mt, mamillary tract.
Specificity of the SubM–VLO circuit for vagal processing via the Pa5
Airway afferents that terminate in the medullary Pa5 are known to project to the SubM (McGovern et al., 2015a). Using the same anterograde transsynaptic tracing approach with our conditional HSV1 H129floxed virus, we determined the spatial distribution of H129 tdT-infected neurons within the SubM that are in receipt of Pa5-projecting airway afferents. Infected neurons spanned a rostrocaudal distribution within the SubM that ranged from −2.30 to −3.20 mm, with the largest density observed at approximately −2.70 mm (Fig. 6G). Interestingly, the majority of these neurons was located in the dorsal half of the SubM nucleus (Fig. 6G,H), consistent with what we observed with the transduction patterns of both the hM3D(Gq) and control AAV experiments. Activation of vagal afferent central terminals in both the NTS and Pa5 was achieved by microinjecting capsaicin into these two brain sites. Following SubM treatment with vehicle, microinjection of capsaicin (0.5 μm) into the medullary caudal NTS (N = 8; Fig. 7A) produced a considerable decrease in respiratory frequency (Fig. 7B,C) and a decrease in blood pressure (Fig. 7B,D). No change was observed in these parameters following activation of the SubM with 5-HT (repeated-measures two-way ANOVA: interaction between time and treatment, F(10) = 0.51, p = 0.8723; time effect, F(10) = 10.64, p < 0.0001; treatment effect, F(1) = 0.08, p = 0.7881; Fig. 7C,D). Microinjection of capsaicin into the Pa5 (N = 8; Fig. 7E) following SubM vehicle treatment decreased respiratory rate (Fig. 7F,G) and produced a substantial increase in blood pressure (Fig. 7F,H). Activation of the SubM with 5-HT was able to block the decrease in respiratory rate (repeated-measures two-way ANOVA: interaction between time and treatment, F(10) = 5.56, p < 0.0001; time effect, F(10) = 1.11, p = 0.3691; treatment effect, F(1) = 208.57, p < 0.0001; Fig. 7G) observed following capsaicin microinjection into the Pa5; however, no effect was seen in blood pressure (Fig. 7H).
Differential effects on stimulation of bulbar vagal integration nuclei mediated by the submedius thalamic nucleus. A, E, Schematic of experimental design showing bilateral microinjections of vehicle or 5-HT into the SubM, preceding microinjection of capsaicin into either the Pa5 (N = 8) or NTS (N = 8) to activate the central terminals of primary nociceptors. B–D, F–H, Representative chart traces show respiratory and blood pressure movements before and after bilateral microinjection of capsaicin (red line) into either the NTS (B) or Pa5 (F) following vehicle microinjection into the SubM. NTS capsaicin-evoked changes in respiratory frequency (C) and blood pressure (D) were unaffected, whereas Pa5 capsaicin-evoked changes in respiratory frequency (G), but not blood pressure (H), were significantly blunted by SubM 5-HT microinjection in anesthetized rats. Data are represented as the mean ± SD. ****p < 0.0001, SubM 5-HT versus SubM vehicle by Bonferroni's multiple-comparisons test (B).
Discussion
We provide compelling evidence that vagal airway nociceptive inputs to the brainstem are regulated by a central descending control system, capable of substantially blunting airway nociceptor reflexes. Essential to this system is circuitry between the SubM and VLO, confirmed in the present study using a combination of conventional and viral neuroanatomical tracing, and conventional and chemogenetic functional studies. Descending inhibition of nociceptive sensory processing orchestrated by the SubM–VLO circuit was demonstrated at the level of vagal sensory processing sites within the brainstem Pa5, where jugular afferents terminate, but not in the caudal NTS, where nodose afferents project. These observations resemble the thalamocortical network regulating the endogenous analgesia system, which exerts comparable regulation of spinal dorsal horn nociceptive processing (Tang et al., 2009). The potential intersection of the endogenous analgesia system with jugular vagal sensory processing in the brainstem offers new insight into the regulation of airway nociceptive reflexes and respiratory sensations.
A SubM–VLO circuit provides both tonic and recruitable top-down control over vagal somatosensory processing
The somatosensory innervation of the nose and pharynx (supplied by trigeminal and glossopharyngeal nerves) is partially preserved in the larynx and proximal airways via innervation from the jugular vagal ganglia, before transitioning to predominantly viscerosensory nodose innervation throughout the remainder of the airways and lungs (McGovern et al., 2015b; Mazzone and Undem, 2016). A somatosensory representation of the airways in the brain has important functional implications, including for the perception of airway irritation. For example, studies of noxious spinal processing have shown that although graded sensory stimuli generally result in graded action potential responses in dorsal horn neurons, the level of sensation experienced may not correlate with these measures as “top-down” neural pathways are capable of modulating nociceptive processing in the spinal dorsal horn (Millan, 2002; Schweinhardt and Bushnell, 2010; Ossipov et al., 2014). One prominent pathway involves descending projections from the prefrontal cortex (e.g., rostral agranular insular cortex and VLO) to the midbrain periaqueductal gray (PAG), a region central to the endogenous analgesia system. Indeed, stimulation of these prefrontal cortex regions in animals, or the dorsolateral prefrontal cortex in humans, induces analgesia (Hardy, 1975; Lee et al., 2015; Cheriyan and Sheets, 2018; Ong et al., 2019; Wu et al., 2020). Although a wide variety of inputs can regulate prefrontal cortical activity, inputs from the SubM have been shown to be a major contributor to prefrontal-mediated analgesia (Zhang et al., 1998; Tang et al., 2009). The present study suggests that this same network regulates aspects of vagal sensory processing in the brainstem.
We used reflex measures of breathing and autonomic tone as indicators of central vagal afferent circuit activity. As previously demonstrated (Driessen et al., 2015, 2020), electrical stimulation of the larynx evoked a vagally mediated cardiorespiratory response that mimicked the laryngeal chemoreflex. Our data confirm previous reports of dense neuroanatomical connectivity between the SubM and VLO (Coffield et al., 1992; Yoshida et al., 1992; Kuramoto et al., 2017). When activated by 5-HT in anesthetized rats, this SubM–VLO circuit produced potent inhibition of breathing (and to a lesser extent cardiovascular) reflexes evoked by directly stimulating the larynx. This functional effect is unlikely to be an artifact unmasked in anesthetized animals as targeted chemogenetic stimulation of the SubM–VLO circuit in conscious rats potently inhibited respiratory responses evoked by inhaling the nociceptive agent capsaicin.
Two observations from our physiological studies are of particular interest. First, no intervention performed in the SubM or VLO modified basal respiratory or cardiovascular measures, suggesting that SubM–VLO circuitry is devoid of regulatory inputs directly to cardiorespiratory premotor/motor neuron pools in the brainstem and spinal cord. The data may also argue that ongoing pulmonary stretch receptor feedback (nodose derived) to the brainstem is also not regulated by SubM–VLO activity, as alterations in this feedback would be expected to modify basal breathing (Nonomura et al., 2017). However, an action mediated at these sites cannot be categorially ruled out without detailed electrophysiological examination. Second, both SubM inhibition by muscimol and VLO lesion demonstrated a level of SubM–VLO circuit-dependent inhibitory tone over vagal sensory processing in the brainstem. This observation argues for the existence of central neuronal mechanisms capable of bidirectionally altering the gain of laryngeal-evoked respiratory reflexes.
Precisely how the SubM–VLO circuit induces the inhibition of vagal nociceptive respiratory responses is unclear. In spinal nociceptive processing, prefrontal neurons descend on PAG neurons, which in turn regulate spinally projecting brainstem neuronal pools in the rostral ventromedial medulla, locus coeruleus, and previously (Millan, 2002; Kim et al., 2018). A similar circuit (or components of it) may be involved in vagal sensory regulation in the brainstem. As we have previously reported, respiratory reflexes can be evoked by microinjections of capsaicin into either the NTS or Pa5, responses that are likely mediated in part by airway vagal afferents terminating in these regions (Mazzone and Canning, 2002; Kubin et al., 2006; Chang et al., 2015; Kim et al., 2020). Concomitant activation of the SubM significantly blunted capsaicin-evoked responses from the Pa5, without having an effect on those mediated by the NTS. This strongly indicates that SubM–VLO descending inhibitory influences are restricted to jugular somatosensory vagal afferent pathways. This observation is perhaps consistent with the disproportionate impact on breathing reflexes over the accompanying cardiovascular responses, of which the latter are dominated by the activation nodose and NTS circuits (Kalia and Richter, 1988; Bonham and McCrimmon, 1990; Kubin et al., 2006; Mazzone and Undem, 2016). Furthermore, in conscious rats, capsaicin inhalation evoked a paradoxical hypoventilation after activation of the SubM–VLO circuit, rather than simply preventing the increase in breathing observed in control animals, conceivably representing an underlying or alternate reflex response not sensitive to SubM–VLO activation.
Vagal inputs may regulate thalamocortical control of the descending analgesia system
The SubM contains glutamatergic neurons that project to the VLO and local interneurons responsive to noxious stimuli from a range of peripheral tissues (Zhang et al., 1998; Huo et al., 2009; Kuramoto et al., 2017). Previous studies have shown that it receives direct nociceptive inputs from laminar one neurons of the trigeminal nucleus and spinal dorsal horn (Craig and Burton, 1981; Ma et al., 1988; Dado and Giesler, 1990; Yoshida et al., 1991; Blomqvist et al., 1992). Using novel transsynaptic herpes viral tracing strategies, we confirm that airway vagal afferent pathways reach the SubM and VLO, terminating in regions that overlap extensively with other defined nociceptive inputs (McGovern et al., 2012, 2015a,b). The anatomic route of these vagal inputs is solely via the Pa5, indicative of only afferents of jugular vagal ganglia origin contributing to nociceptive regulation of the SubM–VLO circuit. Indeed, we previously reported that projections into the SubM and VLO are strikingly absent from nodose–NTS circuits (McGovern et al., 2015a). Visceral inputs via spinal routes to this thalamocortical circuit are well described. For example, both colorectal distension and noxious somatic stimulation have been shown to modulate the activity of the same neurons in the SubM and VLO (Kawakita et al., 1993, 1997; Fu et al., 2002; Barkova et al., 2005). Our data suggest that an alternative input from internal organs may arise from tissues supplied by the jugular vagal ganglia. The circumstances under which vagal input to the SubM–VLO circuit is recruited is not known. Conceivably, as the cervical vagus nerve and especially its proximal branches (e.g., the auricular nerve) contains axons arising from jugular ganglia neurons, our anatomic and functional data linking the jugular sensory pathway with the endogenous analgesia system offers an alternative circuitry via which vagus nerve stimulation protocols modulate pain and spinal nociceptive transmission (Randich and Gebhart, 1992; Ren et al., 1993; Multon and Schoenen, 2005; Straube et al., 2015; Mahadi et al., 2019).
Conclusions and implications for disorders of airway sensory processing
Our results further support the notion that airway tissues have a vagal somatosensory representation in the brain. Our extensive past and present animal and human data argue that this central somatosensory circuitry arises from jugular vagal afferents and involves relays via the Pa5 and ventrobasal thalamus (McGovern et al., 2015a,b; Farrell et al., 2020). Cough is a significant and difficult to treat clinical problem, often ascribed to disordered functioning of nociceptive vagal afferents of the nodose ganglia. However, mounting evidence is challenging this dogma. Animal and human studies support the notion that jugular, not nodose vagal nociceptive, pathways evoke coughing when stimulated (Driessen et al., 2016; Hewitt et al., 2016; Chou et al., 2018; Driessen et al., 2020). Intriguingly, both behavioral studies and clinical trials of new cough suppressants show profound inhibitory effects of placebo conditioning (up to 40% reduction) on cough and measures of airway irritant perception (Abdulqawi et al., 2015; Morice et al., 2019). We have shown in humans that placebo cough suppression is proportional to enhanced prefrontal cortex activity (Leech et al., 2013), in regions also involved in placebo analgesia and the recruitment of the endogenous analgesia system. When considered alongside the present data, these observations point toward a clinically relevant role of jugular vagal afferent pathways and descending control mechanisms in chronic cough. Given that altered central inhibition is common in patients with chronic pain, it will be intriguing to understand how this system functions in disorders impacting the airways. Chronic cough patients and long-term cigarette smokers both show signs of altered central inhibition, significant changes in PAG activity and cough control (Ando et al., 2016, 2019), suggesting that therapies aimed at reducing jugular vagal sensory processing or enhancing central inhibition may prove clinically useful for treating excessive coughing. Our findings may equally offer new insights into other pathophysiological conditions, such as migraine and nociceptive hypersensitivities in head, neck, and thoracic tissues known to receive jugular vagal afferent innervation.
Footnotes
A.E.M. is currently funded by National Health and Medical Research Council (NHMRC) of Australia Fellowship APP1121376. This research was supported by NHMRC Grant 1078943 to S.B.M.
The authors declare no competing financial interests.
- Correspondence should be addressed to Alice E. McGovern at alice.mcgovern{at}unimelb.edu.au.
