PreBötzinger complex neurons drive respiratory modulation of blood pressure and heart rate
Abstract
Heart rate and blood pressure oscillate in phase with respiratory activity. A component of these oscillations is generated centrally, with respiratory neurons entraining the activity of pre-sympathetic and parasympathetic cardiovascular neurons. Using a combination of optogenetic inhibition and excitation in vivo and in situ in rats, as well as neuronal tracing, we demonstrate that preBötzinger Complex (preBötC) neurons, which form the kernel for inspiratory rhythm generation, directly modulate cardiovascular activity. Specifically, inhibitory preBötC neurons modulate cardiac parasympathetic neuron activity whilst excitatory preBötC neurons modulate sympathetic vasomotor neuron activity, generating heart rate and blood pressure oscillations in phase with respiration. Our data reveal yet more functions entrained to the activity of the preBötC, with a role in generating cardiorespiratory oscillations. The findings have implications for cardiovascular pathologies, such as hypertension and heart failure, where respiratory entrainment of heart rate is diminished and respiratory entrainment of blood pressure exaggerated.
Introduction
The respiratory and cardiovascular systems act in synergy to maintain blood gas homeostasis in animals. Beyond their respective separate roles, these systems are also functionally coupled, with oscillations in blood pressure (BP), called Traube-Hering waves, and in heart rate (HR), called respiratory sinus arrhythmia (RSA), that are in phase with the respiratory cycle. This respiratory entrainment of cardiovascular activity is a highly conserved physiological property present in vertebrates, including mammals, fish, amphibians and reptiles (Taylor et al., 2010; Taylor et al., 1999), and was recently found in the early air-breathing primitive lungfish (Monteiro et al., 2018). When within physiological range, Traube-Hering waves are important contributors to vascular tone changes (Briant et al., 2015), RSA reduces cardiac energetic cost (Ben-Tal et al., 2012), and both are believed to optimize blood tissue perfusion and gas exchange. However, an alteration in the amplitude of these oscillations, typically with enlargement of Traube-Hering waves and concurrent abolishment of RSA, occurs in common cardiovascular diseases (e.g. hypertension and heart failure) and contributes to their pathogeneses (Machado et al., 2017; Menuet et al., 2017; Palatini and Julius, 2009; Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology, 1996). The pathological relevance of these oscillations extends to numerous neurological diseases, not superficially linked to cardiovascular function, including depression (Brush et al., 2019) and amyotrophic lateral sclerosis (Merico and Cavinato, 2011).
Traube-Hering waves and RSA are generated by multiple mechanisms, including changes in thoracic pressure and viscerosensory afferent feedback, but their origin has a central component, as they remain in heart- and lung-denervated humans, and in paralyzed and deafferented animal preparations (Dick et al., 2014; Farmer et al., 2016; Menuet et al., 2017; Taha et al., 1995). Traube-Hering waves are generated by bursts of vasomotor sympathetic nerve activity that occur in phase with the respiratory cycle (RespSNA), spanning most of inspiration and the very beginning of post-inspiration (the first phase of expiration). As a result of the slow vascular contraction dynamics, Traube-Hering waves peak during post-inspiration. RespSNA is the result of respiratory modulation of the activity of pre-sympathetic neurons in the rostral ventrolateral medulla oblongata, and in particular of the excitatory inspiratory modulation of adrenergic C1 neurons (Menuet et al., 2017; Moraes et al., 2013). RSA is generated by a cardio-inhibitory vagal parasympathetic tone that decreases during the second part of expiration and into inspiration, resulting in HR oscillations that peak at the end of inspiration. Cardiac vagal preganglionic neurons of the nucleus ambiguus are weakly inhibited during the second part of expiration, and strongly inhibited during inspiration, by both GABAergic and glycinergic inputs (Gilbey et al., 1984; Neff et al., 2003). However, the specific cellular components of the central respiratory pattern generating circuit that are responsible for modulation of both vasomotor pre-sympathetic and cardiac vagal preganglionic neurons are not fully elucidated.
Respiratory and cardiovascular neurons are present in adjacent and slightly intermingled cell columns in the ventral medulla oblongata. The respiratory neurons are classified by the relationship between their activity pattern and phase of the respiratory cycle – for example inspiratory neurons fire action potentials during inspiration. Core to the respiratory column is a small, yet heterogeneous, group of neurons in the preBötzinger Complex (preBötC), which consists of excitatory (glutamatergic) inspiratory neurons, and inhibitory (GABAergic and glycinergic) neurons that can be inspiratory or expiratory (Baertsch et al., 2018; Koizumi et al., 2013; Morgado-Valle et al., 2010; Sherman et al., 2015). The preBötC is critical for respiratory rhythm generation, and also acts as a master oscillator entraining the activity of other oscillators, like the sniffing and whisking oscillators, in a unidirectional fashion (Deschênes et al., 2016; Moore et al., 2013). Our recent work on the pre-synaptic connectome of C1 neurons shows that preBötC neurons have the anatomical connections to provide the inspiratory input to these pre-sympathetic neurons (Dempsey et al., 2017; Menuet et al., 2017). It is also proposed that the strong inspiratory inhibition of cardiac vagal preganglionic neurons comes from the preBötC (Frank et al., 2009).
Using optogenetics, which provides the spatial resolution to restrict opsin expression to a small group of neurons such as those in the preBötC, and temporal control to excite or inhibit neurons during specific phases of the respiratory cycle, we show that the preBötC is a cardiorespiratory hub that drives RespSNA, Traube-Hering waves and RSA.
Results
Injections of lentivirus expressing the light-activated anion channel GtACR2 are targeted to the preBötC
To induce rapid, reversible and powerful inhibition of neural activity in the preBötC, we injected a lentivirus inducing expression of GtACR2 (Govorunova et al., 2015) bilaterally. Injections were directed by electrophysiological recordings to identify inspiratory neurons in the rostral-most part of the ventral inspiratory cell column in the medulla oblongata. Post hoc histological analysis revealed the GtACR2 expression was targeted to the gap in the rostrocaudal column of parvalbumin-expressing neurons (Figure 1A; Alheid et al., 2002), dorsal to tyrosine hydroxylase (TH)-expressing A1-C1 neurons, ventral to the caudal tip of the compact part of the nucleus ambiguus, and at the level containing the highest density of neurokinin type 1 receptor (NK1R)-positive neurons (Figure 1B). Together these markers characterize the preBötC area. GtACR2 was expressed in both excitatory (glutamatergic) and inhibitory (GABAergic and glycinergic) neurons in the preBötC (Figure 1C–E). A heatmap representation of GtACR2 expression in all the 12 animals used for the in situ Working Heart-Brainstem Preparation (WHBP) experiments shows that intersectional GtACR2 expression across all animals is restricted to the preBötC area (Figure 1F).
The preBötC is necessary for inspiratory rhythm generation
The preBötC is mostly studied for its critical role in inspiratory rhythm generation. Prolonged GtACR2-mediated preBötC photoinhibition, either in vivo in anesthetized rats (n = 5) or in the in situ WHBP (n = 12), produced immediate cessation of inspiratory activity and a long-lasting apnea (Figure 2A–C; Figure 3A–C). Following the apnea, whilst preBötC photoinhibition remained activated, inspiratory activity resumed, but with an altered pattern of intercostal muscle (intC) electromyograph (EMG) activity in vivo (Figure 2A–B) or phrenic nerve activity in situ (PNA; Figure 3A–C and H; Figure 3—figure supplement 1A–B). In the WHBP, the longest apnea was produced by 10 and 20 Hz photoinhibition (Figure 3A–C). When PNA resumed, 5–10 Hz photoinhibition induced shorter and smaller PNA bursts, whilst PNA burst characteristics were not altered by 20–50 Hz photoinhibition.
A single 10 ms light pulse at any time during inspiration was sufficient to immediately stop the PNA burst and trigger post-inspiration (Figure 7A–C). The duration of the following expiratory period was correlated to the duration and amplitude of the stopped inspiration, mostly by adjusting the post-inspiratory phase duration. Photoinhibition of the PreBötC only during the expiratory period prolonged the expiratory time for the duration of the photoinhibition, and PNA activity resumed following a short (~1 s) post-photoinhibition delay (Figure 7E).
In whole-cell recordings of GtACR2-expressing neurons (n = 7) in brainstem slices of the dorsal medulla, we observed strong, and sustained, hyperpolarisation with GtACR2 photoactivation (Figure 2F), induced by large inhibitory conductances (Figure 2E). This hyperpolarisation continued for the entire period of light delivery, and in spontaneously active neurons, complete blockade of action potential firing was obtained with 10 ms duration light pulses at 20 Hz. On the basis of these recordings, it is not clear why inspiratory activity resumed during prolonged photoinhibition.
Prolonged preBötC photoinhibition induces decreased vasomotor sympathetic activity and increased cardiac parasympathetic activity
Due to the profound apnea, which would alter blood gas composition and sensory afferent input, detailed analysis of the cardiorespiratory response to preBötC photoinhibition was not feasible in vivo and hence detailed analyses were performed in the WHBP, an oxygenated, perfused preparation devoid of the depressive effects of anaesthesia and where central respiratory perturbations do not influence peripheral blood gases or pH (Menuet et al., 2017; Paton, 1996). During prolonged preBötC photoinhibition, thoracic sympathetic nerve activity (tSNA) showed initial inhibition during apnea, and then a return toward control levels as central respiratory activity resumed (Figure 3A and F; Figure 3—figure supplement 1E and G; Figure 3—figure supplement 2A). This was reflected in an initial decrease in perfusion pressure (PP), followed by an acute increase (Figure 3A and G; Figure 3—figure supplement 1F–G) that were maximal with 20 Hz photoinhibition. We also observed a strong increase in vagal nerve activity (VNA) and decrease in HR (bradycardia) (Figure 2A–B and D; Figure 3A–B and D–E; Figure 3—figure supplement 1C–D). Unlike the sympathetic vasomotor effect, the bradycardia lasted the entire photoinhibition period and was not correlated to the apnea duration (Figure 3—figure supplement 1G). The bradycardia was due to increased cardiac parasympathetic drive (measured as non-resp VNA, Figure 3A and D, Figure 3—figure supplement 2B) as it was abolished by bath application of the muscarinic antagonist, atropine (Figure 4). The biphasic tSNA and PP responses to preBötC photoinhibition were maintained during atropine (Figure 4C).
The preBötC is sympatho-excitatory and parasympatho-inhibitory
To avoid secondary circuit effects arising from sustained photoinhibition, we used single light-pulse preBötC photoinhibition to test direct functional connectivity with sympathetic and parasympathetic activities. Laser-triggered averaging of tSNA and VNA following low-frequency (1 Hz) preBötC photoinhibition confirmed opposite sympathetic and parasympathetic effects (Figure 5; Figure 5—figure supplement 1). PreBötC photoinhibition induced a small sympatho-excitation that peaked ~65 ms following the light pulse, and mainly a strong sympatho-inhibition that was maximal at ~135 ms. PreBötC photoinhibition induced dramatic parasympathetic excitation, which could be seen on VNA recordings in response to individual light pulses (Figure 5A–B). This parasympatho-excitation was very fast, peaking at ~7 ms (Figure 5C).
We conclude that neurons within the preBötC exert direct excitatory effects on vasomotor sympathetic pathways and inhibitory effects on cardiac parasympathic pathways. We next investigated the role played by PreBötC neurons in the respiratory modulation of tSNA (RespSNA), PP (Traube-Hering waves) and HR (RSA).
PreBötC activity contributes to the generation of RespSNA and Traube-Hering waves
When respiration resumed during prolonged preBötC photoinhibition, the return of RespSNA and Traube-Hering waves mirrored the return of PNA inspiratory activity (Figure 3A–C and Figure 6A). Like PNA, both RespSNA and Traube-Hering wave amplitudes were strongly decreased during 10 Hz preBötC photoinhibition, but unaffected during 20 Hz photoinhibition (Figure 6A–B; Figure 6—figure supplement 1A and 6S1D), and overall individual changes in RespSNA area under the curve (AUC) and Traube-Hering wave amplitude correlated closely with changes in PNA amplitude (Figure 6C). Decreased RespSNA duration during 10 Hz photoinhibition was specifically due to a shorter inspiratory component (Figure 6D; Figure 6—figure supplement 1A and C). In contrast, mean tonic tSNA was not altered during preBötC photoinhibition (Figure 6B; Figure 6—figure supplement 1B).
Single light pulses delivered during early inspiration arrested PNA bursts (AUC reduced by ~75%) and almost abolished the associated RespSNA and Traube-Hering wave (Figure 7A–D; Figure 7—figure supplement 1A). PreBötC photoinhibition during early post-inspiration had no impact on PNA or RespSNA and induced a small increase in Traube-Hering wave amplitude. Tetanic preBötC photoinhibition during expiration induced exaggerated RespSNA and Traube-Hering wave amplitude, which is likely due to the increased pulse pressure associated with the bradycardia (Figure 7E).
Together, these results show that the preBötC exerts a direct phasic inspiratory drive to vasomotor sympathetic activity, with a major involvement in the generation of RespSNA and Traube-Hering waves.
PreBötC activity is strongly involved in the generation of RSA
RSA was almost abolished during both 10 Hz and 20 Hz preBötC photoinhibition (Figure 6A–B; Figure 6—figure supplement 1D), and changes in RSA were independent from changes in PNA amplitude (Figure 6C). The loss of baseline RSA evoked by systemic atropine demonstrates that RSA is generated via parasympathetic activity (Figure 4A–B). These results indicate that the preBötC is a source of ongoing respiratory phase-locked inhibitory input to cardiac parasympathetic preganglionic neurons.
Baseline RSA is characterized by peaks in HR at the end of inspiration, decreases during the first part of post-inspiration, and increases continuously during the second part of post-inspiration, late-expiration (E2) and inspiration (Figure 7—figure supplement 1B). Single light pulse preBötC photoinhibition during inspiration did not affect the amplitude of the RSA during the respiratory cycle in which stimulation occurred (Figure 7A–D; Figure 7—figure supplement 1), but rather decreased RSA amplitude in the following respiratory cycle due to a smaller post-photoinhibition HR decrease that didn’t reach baseline level before HR increased again. Photoinhibition of preBötC during early post-inspiration induced a larger decrease in HR which increased the amplitude of the following cycle of RSA. Prolonged photoinhibition during expiration strongly increased VNA, decreased HR and prevented the tachycardic component of RSA until photoinhibition was stopped, whereupon HR increased abruptly prior to the next inspiration (Figure 7E).
Together, these results show that the preBötC exerts a complex cardiac parasympathetic influence that spans both inspiratory and expiratory periods, with a major involvement in the generation of RSA.
PreBötC photoexcitation decreases the respiratory command amplitude, resulting in decreased RespSNA, Traube-Hering waves and RSA
The same injection protocol described for the GtACR2 expression was used to produce restricted expression of channelrhodopsin-2 (ChR2) in preBötC bilaterally (Figure 9A). Strong ChR2-positive terminals were found in the facial motor nucleus (Figure 9B), which is characteristic of neurons in the preBötC region, and not those in the inspiratory rostral ventral respiratory group (Deschênes et al., 2016). ChR2-mediated preBötC photoexcitation (n = 8) increased respiratory frequency without altering inspiratory duration, and decreased PNA inspiratory and VNA post-inspiratory amplitudes (Figure 8A, C and D), therefore reducing each respiratory cycle intensity. We focused on 20 Hz photoexcitation, as this produced maximal effects across all animals.
PreBötC photoexcitation did not alter overall mean tSNA amplitude (which includes both RespSNA and tonic tSNA) (Figure 8F), but disrupted respiratory entrainment of SNA leading to decreased RespSNA amplitude and duration, with shorter inspiratory and post-inspiratory components (Figure 8I–J). As a consequence both Traube-Hering wave amplitude (Figure 8H) and mean PP (Figure 8G) were decreased. Changes in RespSNA magnitude and Traube-Hering wave amplitude induced by preBötC photoexcitation both correlate with changes in PNA inspiratory amplitude (Figure 8K). On the other hand, preBötC photoexcitation increased tonic (non-respiratory) tSNA (Figure 8I), and laser-triggered averaging of 1 Hz preBötC photoexcitation revealed a principal sympatho-excitatory effect (Figure 8L–M).
Photoexcitation of preBötC did not alter mean HR (Figure 8E) but decreased RSA amplitude (Figure 8H), due to loss of the post-inspiratory/expiratory bradycardic component. Changes in RSA induced by preBötC photoexcitation were independent of changes in PNA inspiratory amplitude (Figure 8K). Laser-triggered averaging of VNA during continuous low-frequency (1 Hz) preBötC photoexcitation revealed a small short latency parasympatho-inhibition.
PreBötC neurons project directly to putative pre-sympathetic vasomotor and parasympathetic neurons
To map axonal projections from preBötC neurons we performed small (5 nl), unilateral virus injections under electrophysiological guidance to induce restricted expression of tdTomato in neurons in the preBötC core (n = 2). Successful targeting of respiratory preBötC neurons was indicated by substantial numbers of transduced fibres crossing the midline and terminal fields within the contralateral preBötC (Figure 9C), defining features of functionally identified preBötC neurons (Bouvier et al., 2010; Koizumi et al., 2013). TdTomato-positive preBötC terminals were found in close apposition to both RVLM catecholaminergic C1 neurons, which include pre-sympathetic vasomotor neurons (Figure 9D), and nucleus ambiguus neurons, which include cardiac parasympathetic preganglionic neurons (Figure 9E).
Discussion
Respiratory modulation of sympathetic and parasympathetic activity, leading to Traube-Hering waves and RSA, is a well characterized physiological feature but with little mechanistic understanding. Here, we demonstrate that this respiratory-cardiovascular entrainment is coded at the core of the brainstem respiratory network. In addition to its role in inspiratory rhythm generation, the preBötC provides inspiratory-locked excitatory drive to sympathetic activity and subsequent generation of Traube-Hering waves, and inhibitory drive to parasympathetic outputs, to generate RSA. Photoinhibition of preBötC exerts profound effects on HR and BP in both anesthetized in vivo and in situ rat models, highlighting the importance of this ongoing activity under basal conditions. These data reinforce other recent observations (Deschênes et al., 2016; Moore et al., 2013) supporting the view that oscillatory activity arising from the preBötC provides widespread influences extending well past circuits related to breathing.
Targeting the preBötC
The preBötC was defined as a group of excitatory inspiratory neurons with intrinsic rhythmogenic properties, located in the ventral respiratory column (Peña et al., 2004; Smith et al., 1991). It remains classified as a region with relatively ambiguous anatomic boundaries that contains a heterogeneous cell population composed of rhythmogenic, mostly inspiratory, excitatory and inhibitory neurons (Baertsch et al., 2019; Baertsch et al., 2018; Koizumi et al., 2013; Kuwana et al., 2006; Morgado-Valle et al., 2010; Sherman et al., 2015). No single and exclusive marker of the preBötC has been found so far. We studied neurons in the preBötC using the pan-neuronal CBA promoter that targets neurons irrespective of phenotype. Our anatomical validation of correct transduction (Figures 1 and 9) is supported functionally, as photoinhibition caused immediate cessation of inspiratory activity and a prolonged apnea both in vivo and in situ (Figures 2 and 3), whereas photoexcitation increased respiratory frequency while decreasing inspiratory amplitude (Figure 8). These effects on breathing are consistent with previous optogenetic perturbations of the preBötC (Alsahafi et al., 2015; Cui et al., 2016; Koizumi et al., 2016; Vann et al., 2018). Adjacent to the preBötC are caudal ventrolateral medulla (CVLM) GABAergic neurons, which form part of the sympathetic vasomotor circuit and inhibit rostral ventrolateral medulla pre-sympathetic neurons (Guyenet, 2006), as well as some cardiac parasympathetic preganglionic neurons of the nucleus ambiguus (Dergacheva et al., 2010; Gourine et al., 2016). Some of these neurons could have been transduced with our viral approach; however, our physiological observations do not support this being a substantial influence, as GtACR2-mediated photoinhibition caused inverse effects to those expected if these neurons were transduced (which would be sympathoexcitation/increase in blood pressure and parasympatho-inhibition/tachycardia) (Figures 3 and 5).
When inspiration resumed during photoinhibition, we observed a crenelated activity pattern in inspiratory intC EMG and PNA bursts (Figures 2 and 3). This is likely to result from photoinhibition of some rostral ventral respiratory group phrenic premotor neurons, which intermingle with the caudal pole of the preBötC. Importantly, these neurons are not thought to project to respiratory/cardiovascular neurons other than phrenic motoneurons and lateral reticular nucleus (Wu et al., 2017), so they are unlikely to mediate the cardiac parasympathetic and vasomotor sympathetic effects found in this study.
Caveats related to optically-activated chloride channels
To interrogate interactions between preBötC respiratory circuits and autonomic outputs we used the high temporal resolution of the natural light-activated anion channel GtACR2, which induces large chloride photocurrents with faster kinetics than other inhibitory opsins (Govorunova et al., 2015; Figure 2). Interpretation of our results requires careful consideration, as the chloride electrochemical gradient is heterogeneous among subcellular compartments, and axonal activation of chloride channels can lead to terminal neurotransmitter release (Messier et al., 2018). This phenomenon is unlikely to explain the current results, which are consistent with neural inhibition and were opposed by optogenetic excitation, highlighted by the strong, perfectly inverse sympathetic (tSNA) and parasympathetic (VNA) effects observed with single light pulse ChR2-mediated photoexcitation (Figures 5 and 8). We conclude that in this study and with this cell group, GtACR2 reliably inhibited neuronal activity.
The preBötC provides powerful, ongoing inhibition of cardiac parasympathetic activity
Cardiac parasympathetic preganglionic neurons are located in the nucleus ambiguus and the dorsal motor nucleus of the vagus nerve (Gourine et al., 2016). Among these, only nucleus ambiguus neurons show respiratory-related phasic patterns of activity, with an inhibition of their activity during inspiration, an excitation during post-inspiration and a weak inhibition during late expiration (E2) (Gilbey et al., 1984; Neff et al., 2003). Respiratory inputs into nucleus ambiguus cardiac preganglionic neurons could be widely distributed, although the strongest evidence indicates they come from the pontine Kölliker-Fuse nucleus for the post-inspiratory excitation (Farmer et al., 2016), and from GABAergic and glycinergic neurons that could be located in the preBötC for the inspiratory inhibition (Frank et al., 2009; Neff et al., 2003).
Single light pulse preBötC photoinhibition induced intense vagal activation (Figure 7), whilst prolonged photoinhibition, both in vivo and in situ, caused profound bradycardia, with up to a 70% decrease in HR (Figures 2 and 3), which was blocked by of the parasympatholytic agent atropine (Figure 4). The simplest pathway that could explain this result is a direct inhibitory connection between preBötC and cardiac vagal motoneurons, a scheme supported by our tracing studies (Figure 9) and the short latencies of vagal responses to optical excitation and inhibition (Figures 5 and 8). The preBötC contains both GABAergic and glycinergic neurons (Koizumi et al., 2013; Kuwana et al., 2006; Morgado-Valle et al., 2010) that display four activity patterns: tonic, inspiratory, pre-inspiratory and expiratory (Baertsch et al., 2018; Morgado-Valle et al., 2010; Sherman et al., 2015). It is likely that most of these functional subgroups are involved in cardiac parasympatho-modulation.
Bradycardia occurs independent of when in the respiratory cycle photoinhibition commences, showing that tonic neurons must be involved. Most interestingly, RSA, which is mediated by a progressive decrease of vagal parasympathetic drive during expiration that peaks during inspiration (Bouairi et al., 2004; Dergacheva et al., 2010; Gilbey et al., 1984; Neff et al., 2003), was nearly abolished by preBötC photoinhibition, even after the re-commencement of PNA during prolonged photoinhibition. RSA was also abolished with prolonged preBötC photoexcitation, while there was a tendency for an increase in HR (Figure 8). This shows that prolonged photostimuli of the preBötC, that are likely to suppress phasic activities in targeted neurons and render them either inactive (photoinhibition) or tonically active (photoexcitation), are able to suppress RSA. This could be explained by suppression of the phasic activity of inhibitory preBötC neurons, composed of a small (~10%) proportion of inhibitory expiratory neurons (Baertsch et al., 2018; Carroll et al., 2013), and a larger proportion of inhibitory pre-inspiratory and inspiratory neurons. As a result, phasic late expiratory and inspiratory inhibition of cardiac vagal parasympathetic activity would be suppressed, therefore suppressing RSA. The preBötC is therefore critically involved in tonic regulation of HR as well as the generation of RSA, through a combination of phasic inhibitory neurons.
The decrease in HR just after inspiration, during the post-inspiratory phase, likely involves an excitatory drive to cardiac parasympathetic preganglionic neurons external to the preBötC, most likely coming from the pontine Kölliker-Fuse nucleus (Farmer et al., 2016; Gilbey et al., 1984). Still, preBötC neurons exert an indirect influence on this, as a single photoinhibitory light pulse during inspiration modulated both the inspiratory and post-inspiratory phases, as well as amplitude of the following cycle of RSA, by reducing the post-inspiratory decrease in HR (Figure 7).
The preBötC provides excitatory drive to vasomotor sympathetic activity
Respiratory modulation of neurons involved in generating vasomotor sympathetic activity is found from spinal sympathetic preganglionic neurons of the intermediolateral cell column (Zhou and Gilbey, 1992) to medullary bulbospinal pre-sympathetic neurons of the RVLM (McAllen, 1987; Moraes et al., 2013) and caudal raphe (Gilbey et al., 1995), as well as in baroactivated CVLM neurons (Mandel and Schreihofer, 2006). Neurons in each of these regions display three main types of respiratory modulation, inspiratory activated, inspiratory inhibited, and post-inspiratory activated, which could be acquired mono-synaptically or poly-synaptically from respiratory neurons.
Prolonged preBötC photoinhibition reduced SNA and PP, but only for the apneic period (Figure 3). Single-pulse preBötC photoinhibition revealed a fast inhibition of SNA, peaking ~140 ms after the light pulse (Figure 5). This effect on SNA is delayed compared to the latencies of responses recorded in VNA, and cannot be attributed to the relatively slow conduction velocity of sympathetic pathways, as SNA responses to photostimulation of sympathetic premotor neurons are also lower than those seen here (Abbott et al., 2009; Menuet et al., 2014). These observations are consistent with a polysynaptic pathway prior to the recorded SNA, which could potentially include recruitment of neurons in the intermediate reticular nucleus (Toor et al., 2019). Also, inhibitory, GABAergic CVLM neurons are respiratory modulated (Mandel and Schreihofer, 2006) and could mediate inspiratory modulation of RVLM pre-sympathetic neurons arising from inhibitory or excitatory inspiratory preBötC neurons. Microinjection of the GABAA receptor agonist muscimol into the CVLM/preBötC abolishes RespSNA and respiratory modulation of RVLM neurons, but it is not clear whether this is a direct effect to CVLM inhibition, or due to the associated total loss of respiratory drive and consequent network effects (Koshiya and Guyenet, 1996). Microinjections of GABAA receptor antagonists in the RVLM enhance RespSNA (Guyenet et al., 1990; Menuet et al., 2017), which, together with our single-pulse preBötC photoinhibition experiments, favor the hypothesis that any respiratory phasic input from the CVLM is likely to inhibit presympathetic RVLM neurons.
Direct recordings show that inspiratory modulation of RVLM C1 neurons is mediated by excitatory post-synaptic potentials (Moraes et al., 2013), and we found evidence that neurons in the preBötC, albeit not definitively demonstrated to be inspiratory-modulated neurons, send direct projections onto C1 neurons (Figure 9). This confirms our previous data showing retrograde trans-synaptic labeling of preBötC NK1R neurons from pre-sympathetic C1 neurons (Dempsey et al., 2017; Menuet et al., 2017). Upon resumption of PNA, but during continued photoinhibition, the amplitude of RespSNA and Traube-Hering waves were decreased, strongly correlating with PNA amplitude (Figure 6). Single-pulse preBötC photoinhibition during inspiration decreased the concomitant RespSNA burst and subsequent Traube-Hering wave (Figure 7). During preBötC photoexcitation, inspiratory amplitude was decreased and this also correlated with decreased RespSNA and Traube-Hering waves amplitudes (Figure 8). We therefore conclude the preBötC provides ongoing excitatory drive to vasomotor sympathetic activity, most likely via pre-sympathetic C1 neurons whose activity correlates with inspiration (Menuet et al., 2017; Moraes et al., 2013), and drives RespSNA and Traube-Hering waves (Marina et al., 2011; Menuet et al., 2017; Moraes et al., 2017).
Single-pulse preBötC photoinhibition had no effect on RespSNA when applied during post-inspiration (Figure 7). Part of the post-inspiratory component of RespSNA originates from post-inspiratory neurons in the intermediate reticular nucleus (Toor et al., 2019), and similarly to RSA, the dorsolateral pons could also be involved (Baekey et al., 2008; Dick et al., 2009). Indeed, ponto-medullary transections strongly decreased RespSNA and Traube-Hering waves. Yet, this also caused a switch from eupnea to an apneustic breathing pattern, with prolonged and square shaped phrenic nerve activity. This highlights the difficulty to interrogate functional circuitry entrained by respiratory neurons, as alteration of the activity of the respiratory network will create rebound network effects and alterations of breathing, making it difficult to disentangle the direct versus indirect mechanisms inlvolving the neurons studied. Ideally, future experiments will enable modulation of individual synapses between respiratory neurons and neurons regulating sympathetic (or parasympathetic) activity.
Different preBötC respiratory neurons affect SNA compared to cardiac vagal activity
The current view is that excitatory preBötC neurons generate the inspiratory rhythm (Bouvier et al., 2010; Vann et al., 2018; Wang et al., 2014), while inhibitory preBötC neurons shape its pattern and frequency (Baertsch et al., 2018; Cregg et al., 2017; Sherman et al., 2015). Here, prolonged preBötC photoinhibition induced apnea but not for the entire photoinhibition period. Previous investigators have reported prolonged apnea upon preBötC inhibition induced by pharmacogenetic (Tan et al., 2008) or (brief) optogenetic neural inhibition (Vann et al., 2018). Like the inspiratory activity, tSNA also recovered respiratory modulation when inspiration resumed. By stark contrast, vagal activation and bradycardia remained affected for the entire period of photoinhibition. This observation provides strong support for affecting two different cell populations, an excitatory one involved in inspiratory rhythmogenesis and SNA modulation, and an inhibitory one involved in inspiratory patterning and cardiac vagal modulation. The most likely explanation for the return of respiratory activity during continued photoinhibition is that rhythmogenic, excitatory inspiratory neurons of the preBötC escaped GtACR2-mediated photoinhibition. The mechanism responsible for the rhythmogenic activity of preBötC is still debated, but involves excitatory pre-inspiratory and inspiratory neurons with cellular and/or micro-network electrophysiological properties particularly prone to burst generation, which could confer particular propensity to escape prolonged inhibition. If this is the mechanism, it appears to be a property of the excitatory neurons.
Respiratory cardiovascular interactions in the clinical setting
Assessment of the short-term variability of HR, and to a lesser extent BP, is common in clinical evaluation of cardiovascular health. The main spectral density found in both HR and BP variability is in phase with the respiratory frequency, in the high-frequency band, and is actually a quantification of RSA and Traube-Hering waves, respectively. These measurements provide an invaluable means for non-invasive assessment of parasympathetic and sympathetic drives to the cardiovascular system and have attracted general community interested due to the connection between HR variability and ‘healthiness’ scores. However, views on the analysis and interpretation of these data vary substantially (Hayano and Yuda, 2019). At the core of uncertainty regarding the value of HR and BP variability lies a lack of understanding of the mechanisms driving RSA and Traube-Hering waves.
Regardless of how it is measured, substantial evidence links altered respiratory modulation of autonomic activity with cardiovascular diseases. Exaggerated inspiratory modulation of sympathetic activity leads to increased Traube-Hering waves and drives the development of hypertension (Menuet et al., 2017). This study identifies a key source of this excitatory inspiratory modulation, and while modifying preBötC activity is not conceivable given its critical role, targeted therapeutics that only impact the preBötC to C1 neuron connection could be developed to prevent hypertension development. A common feature of many cardiovascular diseases, and prognostic indicator of poor survival, is reduced RSA, often in concert with increased Traube-Hering waves (Palatini and Julius, 2009; Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology, 1996). Our study also identifies the preBötC as the likely source of this dual alteration of cardiac parasympathetic and vasomotor sympathetic respiratory drives.
Conclusion
Whilst initially characterized as a key oscillator for driving inspiratory activity, the preBötC has already been shown to have more widespread functions, acting as a master oscillator driving other activities, such as sniffing and whisking (Deschênes et al., 2016; Moore et al., 2013). Redundancy is usually a characteristic feature of neuronal networks, as it provides robustness to perturbations (Li et al., 2016). It is therefore surprising to find that the preBötC is so critical for respiratory, cardiovascular and other oscillations. Our study provides evidence for another layer of influence for the preBötC as a hub oscillator that synchronizes cardio-respiratory function. This coordination is a fundamental physiological phenomenon observed across many different phyla that is likely to be of key importance for optimal perfusion, and function, of tissues.
Materials and methods
Animal experiments
Request a detailed protocolExperiments were conducted in accordance with the Australian National Health and Medical Research Council ‘Code of Practice for the Care and Use of Animals for Scientific Purposes’ and were approved by the University of Melbourne Animal Research Ethics and Biosafety Committees (ethics ID #1413273, #1614009, #1814599 and Florey 16-040). All experiments were performed on male Sprague Dawley (SD) rats. Animals were housed with a 12-h light-dark cycle (06H30 to 18H30), at a constant temperature (22 ± 1°C) with ad libitum access to standard rat chow and water.
Viruses
To perform optogenetic inhibition or excitation of preBötC neurons, we used viral-mediated transgenesis to drive the expression of two opsins, the guillardia theta anion channelrhodopsin 2 (GtACR2) and the humanized ChR2 H134R, fused to the enhanced yellow fluorescent protein (YFP) and tdTomato reporter genes, respectively. To trace preBötC terminals, a virus driving only tdTomato expression was used. We used the pan-neuronal promoter chicken ß-actin (CBA), either alone, or as part of the larger synthetic promoter CAG. Indeed, CAG induces the same expression profile as CBA, since CBA is the promoter component (A) of CAG (C being the cytomegalovirus early enhancer element, G being the splice acceptor of the rabbit ß-globulin gene). The viral vectors used were: Lv-CBA-GtACR2-YFP-WPRE (titre 6.95x109 IU/ml), AAV-CBA-GtACR2-YFP-WPRE (titre 9.65x1011 VP/ml), AAV-CAG-hChR2-H134R-tdTomato (titre 2.23x1010 VP/ml), AAV-CBA-tdTomato (6.32x1012 VP/ml). The GtACR2 constructs were cloned in-house from the pFUGW-hGtACR2-EYFP plasmid, which was a gift from John Spudich (Addgene plasmid # 67877; http://n2t.net/addgene:67877; RRID:Addgene_67877) (Govorunova et al., 2015), by replacing the UbiquitinC promoter with CBA. The ChR2-expressing virus was made from the pAAV-CAG-hChR2-H134R-tdTomato plasmid, which was a gift from Karel Svoboda (Addgene plasmid #28017; http://n2t.net/addgene:28017; RRID:Addgene_28017) (Mao et al., 2011). The tdTomato-expressing virus was made from pAM.DCA.spe.tdTomato plasmid.
Brainstem virus microinjections
Request a detailed protocolThirty to sixty minutes prior to surgery, 21-day-old (P21) rats were injected with a non-steroidal, anti-inflammatory drug (meloxicam, 1mg/kg, s.c., Metacam, Boehringer Ingelheim, Sydney, Australia) and then lightly anesthetized by inhalation of isoflurane in an induction box (Rhodia Australia Pty. Ltd., Notting Hill, Australia) prior to intramuscular injection of a mixture of ketamine (60mg/kg, i.m.; Lyppard, Dingley, Australia) and medetomidine (250 μg/kg, i.m.; Pfizer Animal Health, West Ryde, Australia). Surgery was initiated once a deep surgical level of anesthesia was obtained, as evidenced by loss of the pedal withdrawal and corneal reflexes. Body temperature was maintained at 37.5°C with a heat pad (TC-1000 Temperature controller, CWE Inc, USA) that was covered with an autoclaved absorbent pad. Rats were then placed in a stereotaxic frame with the nose ventro-flexed (40°). A midline incision was made over the occipital bone, portions of the bone overlying the cerebellum were removed with a dental drill, and the atlanto-occipital membrane was opened to expose the dorsal brainstem. Glass micropipettes (~30 µm tip diameter) were back-filled with injectate and connected via a silver wire to an electrophysiology amplifier (NeuroLog system, Digitimer Ltd., Welwyn Garden City, UK) to enable recording of neuronal activity. The pipette (20° angled in the rostro-caudal axis, tip pointing forward) was descended into the brainstem 1.5 mm lateral to the midline and extracellular recordings of multiunit activity used to functionally map the preBötC, defined as the most rostral point at which vigorous inspiratory-locked activity was isolated (typically 0.2 ± 0.3 mm rostral, 1.4 ± 0.3 mm ventral to the calamus scriptorius). Bilateral injections were performed for GtACR2 and ChR2 expressions, unilateral injection for tdTomato expression. For in vitro slice electrophysiology, GtACR2 expression was induced in the solitary tract nucleus (NTS). NTS injection coordinates, relative to calamus scriptorius for antero-posterior and medio-lateral coordinates, with all injections performed 0.4mm ventral to the dorsal surface of the brainstem at each injection coordinate: 0.4mm and 0.1mm caudal (midline); 0.1mm rostral, 0.25mm lateral (both sides); 0.4mm rostral, 0.4mm lateral (both sides). A picospritzer (World Precision Instruments, Sarasota, USA) was used to microinject viral vectors driving GtACR2-YFP (50 nl per injection site over 5 min), ChR2-tdTomato (100 nl per injection site over 1 min) or tdTomato alone (5 nl per injection site over 1 min). Pipettes were left in place for 5 min after injection, the pipette withdrawn, and wounds closed with sterile sutures, and anesthesia reversed with atipamazole (1 mg/kg, i.m., Antisedan, Pfizer Animal Health, West Ryde, Australia). Physiological experiments (working heart and brainstem preparation (WHBP), n = 12 GtACR2 and n = 8 ChR2; in vivo electrophysiology, n = 5 Lv-GtACR2 and n=3 AAV-GtACR2; and in vitro whole cell recordings, n = 7) were performed 10 days after viral injections, while animals used for anterograde tracing (n = 2) were sacrificed 20 days after vector injection.
In vivo electrophysiology
Request a detailed protocolLv-GtACR2 injected rats (n = 5), Figure 2A-D: 30 to 60 min prior to surgery, rats were injected with a non-steroidal, anti-inflammatory drug (meloxicam, 1mg/kg, s.c., Metacam, Boehringer Ingelheim, Sydney, Australia) and then lightly anesthetized by inhalation of isoflurane in an induction box (Rhodia Australia Pty. Ltd., Notting Hill, Australia) prior to intramuscular injection of a mixture of ketamine (60mg/kg, i.m.; Lyppard, Dingley, Australia) and medetomidine (250 μg/kg, i.m.; Pfizer Animal Health, West Ryde, Australia). Surgery was initiated once a deep surgical level of anesthesia was obtained, as evidenced by loss of the pedal withdrawal and corneal reflexes. Body temperature was maintained at 37.5°C with a heat pad (TC-1000 Temperature controller, CWE Inc, USA) that was covered with an autoclaved absorbent pad. Rats were placed in a stereotaxic frame with the nose slightly dorso-flexed (mouth bar set at 0). A small incision was made over the dorso-lateral thorax and bipolar stainless steel hook electrodes were placed in consecutive intercostal muscles with a ground electrode positioned under the skin. Recordings of inspiratory intercostal activity were amplified, filtered (50–1500 Hz band-pass, NeuroLog system, Digitimer Ltd., Welwyn Garden City, UK), sampled (5 kHz) and digitized (Power 1401, Cambridge Electrical Design, Cambridge, UK), and rectified and integrated using Spike2 software (Cambridge Electrical Design, Cambridge, UK). Heart rate (HR) was derived from the same biopotential using a window discriminator to trigger from the R-wave of the electrocardiogram.
AAV-GtACR2 injected rats (n = 3), Figure 2—figure supplement 1: 30 to 60 min prior to surgery, rats were injected with a non-steroidal, anti-inflammatory drug (meloxicam, 1mg/kg, s.c., Metacam, Boehringer Ingelheim, Sydney, Australia) and then lightly anesthetized by inhalation of isoflurane in an induction box (Rhodia Australia Pty. Ltd., Notting Hill, Australia). The rats were transferred to a surgical table and the nose placed in a nose cone delivering 3-3.5% isoflurane in 100% oxygen at a flow rate of 0.8 L/min. Surgery was initiated once a deep surgical level of anesthesia was obtained, as evidenced by loss of the pedal withdrawal and corneal reflexes. Body temperature was maintained at 37.5°C with a heat pad (TC-1000 Temperature controller, CWE Inc, USA) that was covered with an autoclaved absorbent pad. The femoral artery and femoral vein were cannulated for measurement of arterial pressure and drug administration, respectively (PE10 tubing (ID 0.28 x OD 0.61 mm) connected to PE50 (0.7 x 1.45 mm)). For diaphragm electromyography (DiaEMG) recordings, a lateral abdominal incision was made and two nylon-insulated stainless-steel wire electrodes ending with suture pads were placed in the costal diaphragm. Isoflurane was slowly replaced by urethane (1.2 mg/kg i.v., SigmaAldrich, St Louis, USA), following which rats were tracheotomized, and oxygen (100%) directed over the tracheotomy cannula. Rats were placed in a stereotaxic frame with the nose slightly dorso-flexed (mouth bar set at 0). The arterial catheter was connected to a pressure transducer (P23 Db, Statham Gould, USA), and recordings of DiaEMG were amplified and filtered (500-1000 Hz band-pass, NeuroLog System, Digitimer Ltd., UK). Both signals were sampled (5 kHz) and digitized (Power 1401, Cambridge Electrical Design, Cambridge, UK), and monitored using Spike2 software (Cambridge Electrical Design, Cambridge, UK).
Then, for both groups: A midline incision was made over the occipital bone, portions of the bone overlying the cerebellum were removed with a dental drill, and bilateral optical fibers (200 µm diameter; Doric Lenses, Quebec, Canada) connected to a 473nm DPSS Laser (Shanghai Laser & Optics Century Co., Shanghai, China) were descended vertically. Optical fibers were positioned 1.5 mm lateral (lambda level), and descended slowly while photostimulating at 20 Hz (10 ms pulses) until maximal respiratory and bradycardic effects were observed. The experimental protocol was then started after a 15 min period of acclimatization.
In vitro slice electrophysiology
Request a detailed protocolRats were deeply anesthetized with 5% isoflurane and the medulla removed, blocked and rapidly cooled in artificial cerebrospinal fluid (aCSF, 2 °C). The medulla was sectioned into 250 µm horizontal slices containing the transduced neurons (Leica VT1200s). Slices were continuously perfused in aCSF solution, containing (in mM): NaCl, 125; KCl, 3; KH2PO4, 1.2; MgSO4, 1.2; NaHCO3, 25; dextrose, 10; CaCl2, 2; (300 mOsmol), bubbled with 95% O2, 5% CO2, at 32°C. Recording pipettes (3.5 – 4.5 MΩ) contained a low Cl- internal solution (10 mM Cl-) containing (in mM): NaCl, 6; NaOH, 4; KOH, potassium gluconate, 130; EGTA, 11; CaCl2, 1; HEPES, 10; MgCl2, 1; 0.1% biocytin (pH 7.3, 290 mOsmol). As a consequence ECl- was -68.8 mV. Pipettes were visually guided to NTS transduced YFP-neurons using a fixed stage scope (Zeiss Examiner) and camera (Rolera EM-C2, Q-Imaging). Whole-cell recordings were made in either voltage clamp (holding voltage -60 mV or -40 mV) or current clamp mode (MultiClamp 700B and pClamp 10.3, Molecular Devices). Signals were sampled at 20 kHz and filtered at 10 kHz. Liquid junction potentials were not corrected (-6.2 mV at 32°C). Pulses of light (465 nm, 10 mW at fibre tip) from a light emitting diode (LED, Plexon) were delivered via fibre optic affixed to a micro-manipulator and positioned to illuminate the recorded neurons.
Working Heart-Brainstem preparation
Request a detailed protocolRats were anesthetized deeply with isoflurane until loss of the pedal withdrawal reflex, bisected below the diaphragm, exsanguinated, cooled in Ringer solution on ice (composition in mM: 125 NaCl, 24 NaHCO3, 5 KCl, 2.5 CaCl2, 1.25 MgSO4, 1.25 KH2PO4 and 10 dextrose, pH 7.3 after saturation with carbogen gas (chemicals were purchased from Sigma-Aldrich, Australia)) and decerebrated precollicularly. The lungs were removed and the descending aorta was isolated and cleaned. Retrograde perfusion of the thorax and head was achieved via a double-lumen catheter (ø 1.25 mm, DLR-4, Braintree Scientific, Braintree, USA) inserted into the descending aorta. The perfusate was Ringer solution containing Ficoll (1.25%), vasopressin acetate (0.5 nM, V9879, Sigma-Aldrich, Australia) and vecuronium bromide (2–4 μg.ml−1, Organon Teknika, Cambridge, UK) warmed to 31°C and gassed with carbogen (closed loop reperfusion circuit, 25–30 ml/min). The second lumen of the cannula was connected to a transducer to monitor perfusion pressure (PP) in the aorta. The head of the preparation was fixed with ear and mouth bars with the nose ventro-flexed (40°). Simultaneous recordings of phrenic nerve activity (PNA), cervical vagus nerve (left side) activity (VNA) and thoracic sympathetic chain activity (tSNA, between T8-10) were obtained using glass suction electrodes. Some WHBPs for both GtACR2 and ChR2 experiments were performed without VNA recordings, to preserve a fully functional vagal cardiac parasympathetic modulation by keeping both vagus nerves intact. No differences in HR and respiratory sinus arrhythmia were found between WHBPs with either unilateral or bilateral intact vagus nerves, therefore experiments were pooled for analysis. Neurograms were amplified, filtered (50–1500 Hz band-pass, NeuroLog system, Digitimer Ltd., Welwyn Garden City, UK), rectified and integrated using Spike2 (Cambridge Electrical Design, Cambridge, UK). HR was derived by using a window discriminator to trigger from the R-wave of the electrocardiogram recorded simultaneously through the phrenic nerve electrode. Bilateral optical fibers (200 µm diameter; Doric Lenses, Quebec, Canada) connected to a 473 nm DPSS Laser (Shanghai Laser and Optics Century Co., Shanghai, China) were descended, 20° angled in the rostro-caudal axis, tip pointing forward, to a position 400 µm dorsal to the injection coordinates. Pilot mapping experiments indicated that light delivery to these coordinates was associated with the largest physiological effects.
Optogenetic stimulation
Request a detailed protocolA light meter (PM20A fiber power meter, Thorlabs, Newton, NJ, USA) was used to calibrate light output from each optical fiber to 10 mW as this was found to induce maximal effects in our experimental set-up without any adverse effect on the tissue (Menuet et al., 2014). Photostimulation consisted of: i) frequency-response trials to compare physiological effects induced by GtACR2 photoactivation in the WHBP (5, 10, 20 and 50 Hz, 10 ms pulses, with respectively 300, 500, 1000 and 1500 pulses at each frequency) and at the single-cell level in vitro (1, 10, 20 and 50 Hz; 1, 10 and 100 ms pulses, 30s each). The results of these trials indicated that maximal physiological effects were evoked by 10 ms light pulses delivered at 20 Hz; in subsequent experiments in vivo or in the WHBP we used tetanic (20 Hz, 300-600 pulses), intermittent (1 Hz, 300 pulses), or single triggered 10ms pulses delivered at specific phases of the respiratory cycle. All photostimulation protocols were repeated multiple times, in a randomized manner, in each animal (technical replication). We found strong intra-preperation reproducibility for all protocols.
Histology
Request a detailed protocolAt the end of WHBP experiments, brainstems were removed, fixed by immersion in 4% formaldehyde in 0.1M sodium phosphate buffer (FA), for 24 hr at 4°C and cryoprotected in 20% sucrose. At the end of the in vivo and tracing experiments, animals were perfused transcardially with 200 ml phosphate-buffered saline followed by 400 ml of 4% FA. Brainstems were removed, post-fixed for 12 hr in 4% FA at 4°C and cryoprotected in 20% sucrose. Brainstems were mostly cut coronally (40 µm sections), except for parvalbumin labeling where they were cut sagitally (50 µm sections), using a cryostat. Fluorescence immunohistochemistry was performed as previously described (Chen et al., 2010; Menuet et al., 2014; Sevigny et al., 2012). Primary antibodies used were rabbit anti-TH (1:5000, Merck-Millipore Bioscience Research Reagents, Bayswater, VIC, Australia, AB152), chicken anti-GFP (1:5000, Abcam, Melbourne, VIC, Australia, AB13970), mouse anti-parvalbumin (1:10,000, Merck-Millipore Bioscience Research Reagents, Bayswater, VIC, Australia, MAB1572), Rabbit anti-NK1R (1:5000, Merck-Sigma-Aldrich, Castle Hill, NSW, Australia S8305), Rabbit anti-DsRed (1:5000, Takara Bio, Mountain View, CA, USA, 632496) and Goat anti-ChAT (1:1000; Merck-Millipore Bioscience Research Reagents, Bayswater, VIC, Australia, AB144P). The secondary antibodies used were Cy3-conjugated donkey anti-rabbit (1:500, Jackson ImmunoResearch Laboratories, West Grove, USA), Cy3-conjugated donkey anti-mouse (1:500, Jackson ImmunoResearch Laboratories, West Grove, USA), AlexaFluor-488 donkey anti-rabbit (1:500, Jackson ImmunoResearch Laboratories, West Grove, USA), AlexaFluor-488 donkey anti-mouse (1:500, Jackson ImmunoResearch Laboratories, West Grove, USA) and AlexaFluor-488 donkey anti-chicken (1:500, Jackson ImmunoResearch Laboratories, West Grove, USA).
In situ hybridization cRNA probes for glycine transporter 2 (GlyT2) and glutamic acid decarboxylase 67 (GAD67) were synthesized as previously described, permitting identification of glycinergic and GABAergic neurons respectively by in situ hybridization (Bowman et al., 2013; Le et al., 2016). A digoxigenin-labeled cRNA probe for vesicular glutamate transporter 2 (VGluT2), a marker of glutamate synthesis, was synthesized as follows: cDNA template (886 bp, Genebank reference sequence NM_053427) was amplified using the PCR primers (lower case) with T7 and Sp6 RNA polymerase promotors (upper case) attached F:GGATCCATTTAGGTGACACTATAGAAGtcaatgaaatccaacgtcca; R:GAATTCTAATACGACTCACTATAGGGAGAcaagagcacaggacaccaaa.
Following purification of the DNA template with gel extraction, antisense riboprobes were in vitro transcribed (Epicentre Technologies, Madison, WI) incorporating digoxigenin-11-UTP (Roche Applied Science) and validated as previously described (Kumar et al., 2012).
Sections containing GtACR2-YFP-labeled neurons were incubated in pre-hybridization buffer (50% formamide, 100 µg/ml heparin, 5 x SSC, pH 7.0, 1 x Denhardt’s solution, 250 µg/ml herring sperm DNA, 100 µg/ml yeast tRNA, 5% dextran sulphate, 0.1% Tween-20, reagents obtained from Sigma, Australia unless otherwise indicated) at 37°C (30 min), then 58°C (1 hr) before hybridization with gentle agitation with GlyT2, GAD67 or VGluT2 riboprobe (1000 ng/ml) at 58°C (12–18 hr). Sections were washed in 2 x SSC buffer with 0.1% Tween-20, followed by 0.2 x SSC buffer with 0.1% Tween-20, and then maleic acid buffer (0.1 M maleic acid, 0.15 M NaCl, 0.1% Tween-20). The tissue was then blocked in maleic acid buffer containing 2% Boehringer blocking reagent (Roche Applied Science, Manheim, Germany) and 10% normal horse serum.
Rabbit anti-GFP (1:500, Life Technologies, Scoresby, Australia, A-6455) was added to the blocking buffer and incubated at 4°C (24 hr), then at room temperature (4 hr). Sections were washed in TPBS (Tris–HCl 10 mM, sodium phosphate buffer 10 mM, 0.9% NaCl, pH 7.4, 3 × 30 min) and incubated overnight with AlexaFluor-488 donkey anti-rabbit IgG secondary antibody (1:500, Life Technologies, Scoresby, Australia, A21206) with 2% normal horse serum. DIG-labeled neurons were revealed by incubation in NTMT (0.1 M NaCl, 0.1 M Tris-HCl pH 9.5, 0.1 M MgCl2, 0.1% Tween-20, 2 mM tetramisole HCl) containing nitro blue tetrazolium (Roche Applied Science) and 5-bromo-4-chloro-3-indolyl phosphate salts (Roche Applied Science). The reaction was stopped by washing (0.1M Tris, 1 mM EDTA, pH 8.5, 3 × 15 min) when DIG-labelling was intense with minimal background staining. No labeling was seen when the sense probe was substituted for the anti-sense probe.
Mapping of reporter expression
Request a detailed protocolA complete map of the distribution of GtACR2-YFP expression from the 12 WHBP experiments was generated from 1:4 series of sections co-labeled with ChAT, TH and NK1R to assist with alignment of the sections into the Paxinos-Watson atlas (Paxinos and Watson, 2004). Heatmaps of pooled GtACR2-YFP expression, imaged using Zeiss Axio Imager.M2 microscope with an Apocromat 20x/0.8 air objective (Carl-Zeiss, NSW, Australia) and captured using Axiocam 506 mono camera and ZEN 2.6 (blue edition) imaging software, were made using Affinity Designer software (Serif Ltd., UK). Sections from four representative animals used in the WHBP were processed to examine colocalization of preBötC GtACR2-YFP with GAD67, VGluT2 and GlyT2 mRNA. Sections were imaged under epifluoresence and brightfield (Zeiss Z1, 10x objective, 0.3 NA). The expression profile of ChR2-tdTomato and GtACR2-YFP used for in vivo experiments were similar to that of GtACR2-YFP from WHBP experiments, so they were mapped qualitatively rather than quantitatively. The brainstems from two animals from the GtACR2-YFP in vivo experiments were cut sagitally for better visualization of rostro-caudal GtACR2-YFP expression.
For the tracing study, tdTomato expression at the injection site was examined under Zeiss Axio Imager.M2 microscope with an Apocromat 20x/0.8 air objective (Carl-Zeiss, NSW, Australia) and captured using Axiocam 506 mono camera and ZEN 2.6 (blue edition) imaging software. TdTomato terminals were imaged with a Zeiss LSM880, Axio Imager 2 confocal microscope, with a Plan-Apochromat 63x/1.4 oil DIC M27 objective, in a Z-stack array, acquiring up to 8 consecutive overlapping slices. All channels were acquired at a consistent pinhole setting of 1.0 μm. Each optical section was assessed for close apposition between axons labelled by tdTomato and ChAT-immunoreactive nucleus ambiguus vagal preganglionic neurons or TH-immunoreactive C1 pre-sympathetic neurons. Maximum intensity projections of 4–6 image stacks were produced using the Image J software (NIH, USA) to better represent the distribution of preBötC terminal labeling.
Data analysis and statistics
Request a detailed protocolNerve signals were rectified and integrated with a 50 ms time constant. For low-frequency photostimulation (1 Hz) laser triggered averaging of tSNA and VNA, the signals were rectified and integrated with a 5 ms time constant. For triggered photostimulations, PNA signals were rectified and integrated with a 0.5 ms DC remove time constant and a 100 ms smoothing time constant, to eliminate the ECG background. Phrenic-triggered (end of inspiratory burst) averaging of tSNA, systolic PP and HR were used to quantify respiratory modulation of tSNA (RespSNA), Traube-Hering waves (mmHg) and respiratory sinus arrhythmia (RSA, bpm), respectively. Tonic levels of tSNA were measured using the phrenic-triggered averaging of tSNA, as the average of tSNA during the expiratory period, therefore outside of the burst of respiratory modulation of SNA. Phrenic-triggered (end of inspiratory burst) averaging of VNA was used to quantify the peak post-inspiratory and non-respiratory components of VNA. Details regarding the period analysed are shown on Figure 3—figure supplement 2. For manual photostimulations, to avoid averaging distortion, phrenic-triggered averagings were performed on both consecutive respiratory-related activities to be analysed, and then traces were overlapped and merged at their mean intersection for further analysis. Statistical analysis was done on raw values, and then mean deltas in percentages were quantified (photostimulation vs pre-photostimulation period). Statistics shown on figures presenting only mean delta percentages (for presentation clarity, e.g. Figure 3) were done on the raw values as shown on the associated supplementary figures (e.g. Figure 3—figure supplement 1 and associated data source). Group data are presented as mean ± SEM. Normal distribution of data was verified with the Shapiro-Wilk test, and equal variance of data was tested (SigmaPlot v12, Systat Software Inc, Erkrath, Germany). Data were analyzed using one-way repeated measures or two-way repeated-measures ANOVA, followed by post hoc Holm-Sidak multiple-comparison tests, or paired student’s t tests, for data with normal distribution and equal variance. Data with failed normal distribution or failed equal variance were analyzed using Friedman repeated measures analysis of variance on ranks followed with pairwise multiple comparison Tukey test, or Wilcoxon signed rank test. Linear regression and correlation analysis was used to analyze parameters relationships (Prism v2.0, GraphPad, La Jolla, USA). When missing values were present, for instance because one photostimulation frequency could not be performed or some of the associated data were not appropriate for analysis (noise etc.), it was dealt with accordingly depending on the test used, mostly by excluding the sample from the whole analysis (cf. source data and statistics). All source data and detailed statistics are presented as supplementary information along with each figure. Differences were considered significant at p<0.05.
Data availability
All data generated or analysed during this study are included in the manuscript and supporting files. Source data and statistics files have been provided for Figures 2, 3, 4, 5, 6, 7 and 8.
References
-
Parvalbumin in respiratory neurons of the ventrolateral medulla of the adult ratJournal of Neurocytology 31:693–717.https://doi.org/10.1023/a:1025799830302
-
Optogenetic excitation of preBötzinger complex neurons potently drives inspiratory activity in vivoThe Journal of Physiology 593:3673–3692.https://doi.org/10.1113/JP270471
-
Respiratory sinus arrhythmia in freely moving and anesthetized ratsJournal of Applied Physiology 97:1431–1436.https://doi.org/10.1152/japplphysiol.00277.2004
-
Hindbrain interneurons and axon guidance signaling critical for breathingNature Neuroscience 13:1066–1074.https://doi.org/10.1038/nn.2622
-
Brain sources of inhibitory input to the rat rostral ventrolateral medullaJournal of Comparative Neurology 521:213–232.https://doi.org/10.1002/cne.23175
-
Respiratory modulated sympathetic activity: a putative mechanism for developing vascular resistance?The Journal of Physiology 593:5341–5360.https://doi.org/10.1113/JP271253
-
Patterns of inspiratory phase-dependent activity in the in vitro respiratory networkJournal of Neurophysiology 109:285–295.https://doi.org/10.1152/jn.00619.2012
-
Respiratory modulation of premotor cardiac vagal neurons in the brainstemRespiratory Physiology & Neurobiology 174:102–110.https://doi.org/10.1016/j.resp.2010.05.005
-
Cardio-respiratory coupling depends on the ponsRespiratory Physiology & Neurobiology 168:76–85.https://doi.org/10.1016/j.resp.2009.07.009
-
Cardiorespiratory coupling: common rhythms in cardiac, sympathetic, and respiratory activitiesProgress in Brain Research 209:191–205.https://doi.org/10.1016/B978-0-444-63274-6.00010-2
-
Brainstem sources of cardiac vagal tone and respiratory sinus arrhythmiaThe Journal of Physiology 594:7249–7265.https://doi.org/10.1113/JP273164
-
Respiratory-related discharge patterns of caudal raphe neurones projecting to the upper thoracic spinal cord in the ratJournal of the Autonomic Nervous System 50:263–273.https://doi.org/10.1016/0165-1838(94)00097-4
-
Cardiac vagal preganglionic neurones: an updateAutonomic Neuroscience 199:24–28.https://doi.org/10.1016/j.autneu.2016.06.003
-
Rostral ventrolateral medulla and sympathorespiratory integration in ratsAmerican Journal of Physiology-Regulatory, Integrative and Comparative Physiology 259:R1063–R1074.https://doi.org/10.1152/ajpregu.1990.259.5.R1063
-
The sympathetic control of blood pressureNature Reviews Neuroscience 7:335–346.https://doi.org/10.1038/nrn1902
-
Pitfalls of assessment of autonomic function by heart rate variabilityJournal of Physiological Anthropology 38:3.https://doi.org/10.1186/s40101-019-0193-2
-
Tonic sympathetic chemoreflex after blockade of respiratory rhythmogenesis in the ratThe Journal of Physiology 491 Pt 3:859–869.https://doi.org/10.1113/jphysiol.1996.sp021263
-
BookCombined In Situ Hybridization and Immunohistochemistry in Rat Brain Tissue Using Digoxigenin-Labeled RiboprobesIn: Badoer E, editors. Visualization Techniques: From Immunohistochemistry to Magnetic Resonance Imaging, Neuromethods. Totowa, NJ: Humana Press. pp. 31–52.
-
Electrophysiological and morphological characteristics of GABAergic respiratory neurons in the mouse pre-Bötzinger complexThe European Journal of Neuroscience 23:667–674.https://doi.org/10.1111/j.1460-9568.2006.04591.x
-
Somatostatin 2a receptors are not expressed on functionally identified respiratory neurons in the ventral respiratory column of the ratJournal of Comparative Neurology 524:1384–1398.https://doi.org/10.1002/cne.23912
-
Neurogenic hypertension and the secrets of respirationAmerican Journal of Physiology-Regulatory, Integrative and Comparative Physiology 312:R864–R872.https://doi.org/10.1152/ajpregu.00505.2016
-
Central respiratory modulation of barosensitive neurones in rat caudal ventrolateral medullaThe Journal of Physiology 572:881–896.https://doi.org/10.1113/jphysiol.2005.103622
-
Central respiratory modulation of subretrofacial bulbospinal neurones in the catThe Journal of Physiology 388:533–545.https://doi.org/10.1113/jphysiol.1987.sp016630
-
Catecholaminergic C3 neurons are sympathoexcitatory and involved in glucose homeostasisJournal of Neuroscience 34:15110–15122.https://doi.org/10.1523/JNEUROSCI.3179-14.2014
-
Autonomic dysfunction in the early stage of ALS with bulbar involvementAmyotrophic Lateral Sclerosis 12:363–367.https://doi.org/10.3109/17482968.2011.584628
-
Glycinergic pacemaker neurons in preBötzinger complex of neonatal mouseJournal of Neuroscience 30:3634–3639.https://doi.org/10.1523/JNEUROSCI.3040-09.2010
-
The role of cardiac autonomic function in hypertension and cardiovascular diseaseCurrent Hypertension Reports 11:199–205.https://doi.org/10.1007/s11906-009-0035-4
-
A working heart-brainstem preparation of the mouseJournal of Neuroscience Methods 65:63–68.https://doi.org/10.1016/0165-0270(95)00147-6
-
BookThe Rat Brain in Stereotaxic Coordinates - the New Coronal Set - (5th Edition)Elsevier. ISBN 978-0-08-047412-0.
-
Efferent projections of C3 adrenergic neurons in the rat central nervous systemThe Journal of Comparative Neurology 520:2352–2368.https://doi.org/10.1002/cne.23041
-
Respiratory sinus arrhythmia in humans: an obligatory role for vagal feedback from the lungsJournal of Applied Physiology 78:638–645.https://doi.org/10.1152/jappl.1995.78.2.638
-
Autonomic control of cardiorespiratory interactions in fish, amphibians and reptilesBrazilian Journal of Medical and Biological Research 43:600–610.https://doi.org/10.1590/S0100-879X2010007500044
-
A V0 core neuronal circuit for inspirationNature Communications 8:544.https://doi.org/10.1038/s41467-017-00589-2
Article and author information
Author details
Funding
National Health and Medical Research Council (App #1120477)
- Clément Menuet
- Simon McMullan
- Andrew M Allen
National Health and Medical Research Council (App #1156727)
- Clément Menuet
- Simon McMullan
- Andrew M Allen
Australian Research Council (DP120100920)
- Simon McMullan
- Andrew M Allen
Australian Research Council (DP170104582)
- Andrew M Allen
University of Melbourne (- McKenzie Research Fellowship)
- Clément Menuet
Fondation pour la Recherche Médicale (Fellowship ARF20160936221)
- Clément Menuet
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
The authors acknowledge the facilities and technical assistance of the Biological Optical Microscopy Platform (University of Melbourne) for the confocal microscopy images.
Ethics
Animal experimentation: Experiments were conducted in accordance with the Australian National Health and Medical Research Council 'Code of Practice for the Care and Use of Animals for Scientific Purposes' and were approved by the University of Melbourne Animal Research Ethics and Biosafety Committees (ethics ID #1413273, #1614009, #1814599 and Florey 16-040).
Copyright
© 2020, Menuet et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 5,981
- views
-
- 618
- downloads
-
- 61
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Neuroscience
Granule cells of the cerebellum make up to 175,000 excitatory synapses on a single Purkinje cell, encoding the wide variety of information from the mossy fibre inputs into the cerebellar cortex. The granule cell axon is made of an ascending portion and a long parallel fibre extending at right angles, an architecture suggesting that synapses formed by the two segments of the axon could encode different information. There are controversial indications that ascending axon (AA) and parallel fibre (PF) synapse properties and modalities of plasticity are different. We tested the hypothesis that AA and PF synapses encode different information, and that the association of these distinct inputs to Purkinje cells might be relevant to the circuit and trigger plasticity, similar to the coincident activation of PF and climbing fibre inputs. Here, by recording synaptic currents in Purkinje cells from either proximal or distal granule cells (mostly AA and PF synapses, respectively), we describe a new form of associative plasticity between these two distinct granule cell inputs. We show for the first time that synchronous AA and PF repetitive train stimulation, with inhibition intact, triggers long-term potentiation (LTP) at AA synapses specifically. Furthermore, the timing of the presentation of the two inputs controls the outcome of plasticity and induction requires NMDAR and mGluR1 activation. The long length of the PFs allows us to preferentially activate the two inputs independently, and despite a lack of morphological reconstruction of the connections, these observations reinforce the suggestion that AA and PF synapses have different coding capabilities and plasticity that is associative, enabling effective association of information transmitted via granule cells.
-
- Neuroscience
Sour taste, which is elicited by low pH, may serve to help animals distinguish appetitive from potentially harmful food sources. In all species studied to date, the attractiveness of oral acids is contingent on concentration. Many carboxylic acids are attractive at ecologically relevant concentrations but become aversive beyond some maximal concentration. Recent work found that Drosophila ionotropic receptors IR25a and IR76b expressed by sweet-responsive gustatory receptor neurons (GRNs) in the labellum, a peripheral gustatory organ, mediate appetitive feeding behaviors toward dilute carboxylic acids. Here, we disclose the existence of pharyngeal sensors in Drosophila melanogaster that detect ingested carboxylic acids and are also involved in the appetitive responses to carboxylic acids. These pharyngeal sensors rely on IR51b, IR94a, and IR94h, together with IR25a and IR76b, to drive responses to carboxylic acids. We then demonstrate that optogenetic activation of either Ir94a+ or Ir94h+ GRNs promotes an appetitive feeding response, confirming their contributions to appetitive feeding behavior. Our discovery of internal pharyngeal sour taste receptors opens up new avenues for investigating the internal sensation of tastants in insects.