Postinspiratory and preBötzinger complexes contribute to respiratory-sympathetic coupling in mice before and after chronic intermittent hypoxia

The sympathetic nervous system modulates arterial blood pressure. Individuals with obstructive sleep apnea (OSA) experience numerous nightly hypoxic episodes and exhibit elevated sympathetic activity to the cardiovascular system leading to hypertension. This suggests that OSA disrupts normal respiratory-sympathetic coupling. This study investigates the role of the postinspiratory complex (PiCo) and preBötzinger complex (preBötC) in respiratory-sympathetic coupling under control conditions and following exposure to chronic intermittent hypoxia (CIH) for 21 days (5% O₂–80 bouts/day). The surface of the ventral brainstem was exposed in urethane (1.5 g/kg) anesthetized, spontaneously breathing adult mice. Cholinergic (ChAT), glutamatergic (Vglut2), and neurons that co-express ChAT and Vglut2 at PiCo, as well as Dbx1 and Vglut2 neurons at preBötC, were optogenetically stimulated while recording activity from the diaphragm (DIA), vagus nerve (cVN), and cervical sympathetic nerve (cSN). Following CIH exposure, baseline cSN activity increased, breathing frequency increased, and expiratory time decreased. In control mice, stimulating PiCo specific cholinergic-glutamatergic neurons caused a sympathetic burst during all phases of the respiratory cycle, whereas optogenetic activation of cholinergic-glutamatergic PiCo neurons in CIH mice increased sympathetic activity only during postinspiration and late expiration. Stimulation of glutamatergic PiCo neurons increased cSN activity during the postinspiratory phase in control and CIH mice. Optogenetic stimulation of ChAT containing neurons in the PiCo area did not affect sympathetic activity under control or CIH conditions. Stimulating Dbx1 or Vglut2 neurons in preBötC evoked an inspiration and a concomitant cSN burst under control and CIH conditions. Taken together, these results suggest that PiCo and preBötC contribute to respiratory-sympathetic coupling, which is altered by CIH, and may contribute to the hypertension observed in patients with OSA.

1 Introduction

Obstructive sleep apnea (OSA) represents a significant public health burden and is increasing in prevalence (Franklin and Lindberg, 2015). Cardiopulmonary sequelae of OSA include arterial hypertension, increased ventilatory response to hypercapnia, endothelial dysfunction, atherosclerosis, and disordered swallow-breathing coordination (Fletcher, 2001; Prabhakar et al., 2005; Schulz et al., 2005; Nanduri et al., 2012; Dumitrascu et al., 2013; Huff et al., 2022). Chronic intermittent hypoxia (CIH) is used as an experimental paradigm to model the pathophysiological changes associated with OSA. In rats, CIH increases sympathetic activity and disrupts respiratory-sympathetic coupling resulting in hypertension (Zoccal et al., 2008; Simms et al., 2009; Zoccal et al., 2009; Molkov et al., 2011). Furthermore, when exposure to CIH occurs postnatally, elevated arterial pressure levels and baroreflex dysfunctions can persist into adulthood (Soukhova-O'Hare et al., 2006; Karlen-Amarante et al., 2023). The cardiopulmonary effects of CIH have been attributed to the overexcitation of rostral ventrolateral medulla (RVLM) presympathetic neurons (Karlen-Amarante et al., 2023) via excitatory and inhibitory inputs from respiratory rhythmogenic networks (Lipski et al., 1996; Molkov et al., 2011; Moraes et al., 2012, 2013). However, the exact neuronal mechanisms underlying respiratory modulation of sympathetic activity remain to be characterized.

To investigate the role of central respiratory networks in sympathetic control and how their function is altered by CIH in mice, we focused on two respiratory networks located within the medulla: the postinspiratory complex (PiCo) and the pre-Bötzinger complex (preBötC). The preBötC generates inspiratory activity (Anderson and Ramirez, 2017), and contains excitatory, glutamatergic inspiratory neurons (vesicular glutamate transporter 2; Vglut2) that are primarily derived from progenitor cells characterized by the transcription factor Dbx1 (Bouvier et al., 2010; Gray et al., 2010; Cui et al., 2016). It is well-known that preBötC also regulates cardiorespiratory oscillations in sympathetic activity by providing excitatory input to presympathetic catecholaminergic neurons in the RVLM (Moraes et al., 2013; Menuet et al., 2017, 2020). In contrast, the extent to which PiCo contributes to respiratory-sympathetic coupling in mice has not been investigated. PiCo contains neurons that co-express markers for cholinergic (ChAT) and glutamatergic (Vglut2) neurotransmitters (Anderson et al., 2016) and is important for the control of postinspiratory activity (Anderson et al., 2016) and swallow-breathing coordination (Huff et al., 2023). Chemical inhibition of PiCo has been shown to reduce sympathetic output in rats (Toor et al., 2019). Thus, the present study tested the hypothesis that PiCo and preBötC modulate respiratory-sympathetic coupling in mice, which is altered by CIH.

2 Methods

2.1 Ethical approval and transgenic mice

All experiments and procedures were approved by the Seattle Children’s Research Institute (SCRI) Institutional Animal Care and Use Committee (IACUC) and conducted in accordance with the National Institutes of Health and ARRIVE guidelines (Percie du Sert et al., 2020). Adult (P42 to P152), female (n = 32) and male (n = 34) mice used in this study were bred at SCRI and housed in a temperature (22 ± 1°C) and humidity (25–75%) controlled room with a 12:12 h light–dark cycle. Mice were provided with a standard diet and water ad libitum.

Although many studies of the cardiopulmonary effects of CIH have utilized rats, the present study leverages the genetic homology mice have with humans (O'Brien et al., 1999) to advance understanding of the pathogenesis of OSA (Nanduri et al., 2012; Barnes et al., 2022). All the following Cre mice were crossed with homozygous mice containing a floxed STOP channelrhodopsin-2 (ChR2) fused to a EYFP (Ai32) reporter sequence from Jackson Laboratories (stock number 024109). Vglut2-IRES-cre (Vglut2), Vgat-IRES-cre (Vgat), and ChAT-IRES-cre (ChAT) homozygous breeder lines were obtained from Jackson Laboratories (stock numbers 028863; 016962, and 031661 respectively). Heterozygous Dbx1creERT2 mice were donated by Dr. Christopher Del Negro (College of William and Mary, VA) and a homozygous breeder line was generated at SCRI. Homozygous Dbx1creERT2 dams were plug checked and injected intraperitoneally (IP) at E10.5 with tamoxifen (24 mg/kg) to target preBötC Dbx1 neurons (Baertsch et al., 2019; Huff et al., 2022). ChAT-IRES-cre(Δneo), technically known as 129S-Chat^tm1(cre)Lowl/MwarJ strain, and Vglut2-IRES2-FlpO-D, technically known as 129S-Slc17a6^{tm1.1(flpo)Hze}/J were obtained from Jackson Laboratories (stock numbers 031661 and 030212, respectively). To generate double-transgenic mice, the ChATcre and Vglut2FlpO strains were interbred to generate compound homozygotes, referred to as ChATcre:Vglut2FlpO, in which neurons that have a developmental history of expressing both cholinergic and glutamatergic neurotransmitters express both Cre and Flp recombinases. Experimental groups are referred to as Vglut2cre:Ai32, Vgatcre:Ai32, ChATcre:Ai32, Dbx1cre:Ai32, and ChATcre:Vglut2FlpO:ChR2. All mice used in this study were on the C57BL/6 J (B6) background or back crossed with B6 mice (Vgatcre:Ai32 and Dbx1cre:Ai32).

2.2 Brainstem injection of adeno-associated virus (AAV)

All surgical procedures were performed under aseptic conditions. For targeted adeno-associated virus (AAV) injection, the mice (P42–131) were anesthetized with isoflurane (2%). Adequate anesthesia was confirmed by the absence of hind-paw withdrawal reflex. The hair over the skull was removed and the skin was disinfected. The mice were then placed prone on a stereotaxic apparatus (bite bar set at −3.5 mm; David Kopf Instruments, Tujunga, CA, USA). A bilateral 0.5 mm diameter craniotomy was drilled into the occipital bone caudal to lambda. Viral solutions were loaded into a 1.2 mm internal diameter glass pipette with a 20 μm tip (external diameter). To target the PiCo region with ChR2, the pipette was inserted into the brainstem at the following coordinates: 1.6 mm caudal to lambda, 1.1 mm lateral to the midline, and 4.8 mm below the dorsal surface of the cerebellum. Bilateral 150 nL injections of pAAV-hSyn Con/Fon hChR2(H134R)-EYFP adenovirus vector (cat# 55645-AAV8; AddGene, USA; abbreviated as AAV8-ConFon-ChR2-EYFP) were made at a rate of 50 nL/min using a glass micropipette and an automatic nanoliter injector (Nanoject II, Drummond Scientific Co., Broomall, PA). This AAV is a Cre-on/Flp-on ChR2-EYFP under the synapsin promoter and encodes the photoactivatable cation channel channelrhodopsin-2 (ChR2, H134R) fused to EYFP. The vector was diluted to a final titer of 1 × 10¹³ viral particles/ml with sterile phosphate-buffered saline. Afterward, the skin was sutured, and the animals were treated with ketoprofen (5 mg/kg/day for 2 days). No abnormal behaviors were observed during the recovery period. Mice were given three days to recover from intracranial viral injections before being randomly assigned to CIH or normoxia protocols.

2.3 Chronic intermittent hypoxia (CIH) exposure

Offspring from optogenetic lines were randomly assigned to CIH or normoxia (control) protocols. The CIH protocol used in the present study has been previously described (Huff et al., 2022). Briefly, mouse cages were kept inside a custom-built dual oxygen/nitrogen cycling chamber (Oxycycler, Huff Technologies Inc., Morganfield, KY). One chamber was used for the control group and another for the CIH group. Mice in the CIH chamber were exposed to intermittent episodes of hypoxia consisting of 1 min of nitrogen (N2) injection to reduce the percentage of inspired oxygen (O₂) inside the chamber from 21% to 4.5–5%. Then, compressed room air was injected into the chamber to return the percentage of O₂ to 21% for 5 min before starting a new cycle of hypoxia. The entire 6 min CIH protocol was repeated 80 times per day over 8 h during the mouse’s subjective night. For the remaining 16 h, CIH mice were kept under normoxic conditions (21% O₂). The injection of N2 and compressed air into the chambers was automated by custom computer software (Huff Technologies Inc.) and monitored by O₂ sensors inside the chambers. The gas injection occurred at the top of the chamber to avoid potential stress on the mice. Directly below the CIH chamber, control mice were kept in an identical chamber under normoxic conditions 24 h a day for 21 days. All mice had ad libitum access to food and water during both protocols.

2.4 In-vivo preparation

After CIH or normoxia exposure, mice from both experimental groups were anesthetized with urethane (1.5 g/kg, IP) and placed supine on a custom surgical table with a heating pad to keep mouse body temperature constant at 37°C. Adequate depth of anesthesia was confirmed every 15 min via heart and breathing rates and absence of hind-paw withdrawal reflex. The trachea was exposed through a midline incision and cannulated caudal to the larynx with a curved (180 degree) tracheal tube (24 G). A line delivering 100% O₂ was attached to the cannulated trachea to provide O₂ throughout the experiment. All the recorded nerves were isolated unilaterally, cut distally, and recorded via a unipolar suction electrode connected to a fire-polished, pulled glass pipette filled with aCSF. The corresponding contralateral nerve was preserved and remained intact. The cervical sympathetic nerve (cSN) was isolated and cut caudal to the superior cervical sympathetic ganglion. The cervical vagus (cVN) nerve was isolated, and the trachea and esophagus were cut and detached at the rostral end. The occipital bone was removed, and the ventral medullary surface was bathed with warm (~36°C) artificial cerebrospinal fluid (aCSF; in mM: 118 NaCl, 3 KCl, 25 NaHCO₃, 1 NaH₂PO₄, 1 MgCl₂, 1.5 CaCl₂, 30 D-glucose) equilibrated with carbogen (95% O₂, 5% CO₂) by a peristaltic pump (Dynamax RP-1, Rainin Instrument Co, Emeryville CA, USA). A bipolar electromyogram (EMG) electrode (A-M systems, Carlsborg, USA) was placed in the costal diaphragm muscle (DIA) to record the inspiratory activity and monitor anesthetic depth. The cSN and cVN nerve efferences were recorded using glass suction unipolar electrodes secured in a 3D micromanipulator (Narishige, Tokyo, Japan). All EMG and electroneurogram (ENG) signals were acquired with a A/D converter, (CED micro 1401, Cambridge Electronic Design, Cambridge, UK) amplified (AC Amplifier, model 1700, A-M System, Carlsborg, USA), and band-pass filtered (100–3 kHz). Using Spike 2 software (5 KHz, CED, Cambridge, UK), data was further processed using a band pass filter (200–700 Hz, 40 Hz transition gap), rectified, and smoothed (20 ms).

2.5 Optogenetic stimulation and data analysis

Optical fibers with a diameter of 200 μm were connected to a blue (470 nm) high-powered LED (ThorLabs Inc., Newton, NJ, USA) and positioned bilaterally on the ventral surface above preBötC or PiCo by utilizing the vertebral and basilar arteries as landmarks as described in previous studies (Anderson et al., 2016; Baertsch et al., 2018; Huff et al., 2022, 2023). Power was set at 0.75 mW/mm². To determine the probability of optogenetically evoking diaphragm bursts and sympathetic activity responses, 200 ms light pulses were triggered every 10 s for 4 min (total 25 stims) to stimulate ChAT, Vglut2, and neurons that co-express ChAT and Vglut2 in PiCo or Dbx1 and Vglut2 neurons in preBötC. This resulted in the stimulations occurring randomly across the respiratory cycle. To avoid artificially inflating degrees of freedom, stimulations were averaged within each respiratory phase within each mouse. In Vgat:Ai32 mice, continuous 10 s light pulses were triggered every 60 s to stimulate GABAergic neurons in preBötC.

The frequency of DIA discharge, or respiratory rate (RR), was determined by measuring the time interval over 10 consecutive bursts (expressed as bursts per minute). We quantified the DIA burst duration (time of inspiration, s) and interval between inspiratory discharges (time of expiration, s). Heart rate (HR) was derived from the DIA EMG. Spike 2 was used to detect and analyze the number of spikes with a half-width of 0.5 ms present in the cSN raw signal. In addition, the distribution of sympathetic discharge throughout the respiratory cycle was calculated to quantify the phase-specificity of respiratory-sympathetic coupling. The inspiratory phase was determined by the duration of the DIA burst. Expiration was divided into two phases based on the vagus activity: postinspiration and late expiration.

The cSN activity was recorded and analyzed in absolute units (volts, V). Noise levels were quantified at the end of each experiment and subtracted from the cSN activity. The integrated sympathetic activity signal was used to calculate the amplitude of the responses to the optogenetic stimulation. The amplitude of the cSN recordings were normalized to allow for comparisons between mice. During each mouse’s respective baseline, the peak cSN amplitudes during 10 respiratory cycles were standardized to the highest peak of normal physiological cSN discharge during the respiratory cycles and expressed as a percentage of the maximum cSN peak (% MAX). Then, optogenetically evoked responses in cSN activity were expressed as % MAX of baseline cSN maximum amplitude during inspiration. We have previously shown that stimulation of PiCo neurons has a high probability of evoking a swallow in a respiratory-phase dependent manner (Huff et al., 2023). To evaluate the PiCo contribution to respiratory-sympathetic coupling, only PiCo stimulations that did not trigger a swallow were included in this study.

2.6 Histology

At the end of the experiments, all animals were deeply anesthetized with 5% isoflurane in 100% oxygen and perfused through the ascending aorta with 20 mL of phosphate buffered saline (PB; pH 7.4) followed by 4% phosphate-buffered (0.1 M; pH 7.4; 20 mL) paraformaldehyde (Electron Microscopy Sciences, Fort Washington, PA). The brains were removed and stored in the perfusion fixative for 4 h at 4°C followed by 20% sucrose for 8 h. Coronal brain sections (25 μm) were serially cut using a cryostat and stored in cryoprotectant solution at −20°C (20% glycerol plus 30% ethylene glycol in 50 mL phosphate buffer, pH 7.4) prior to histological processing. All histochemical procedures were done using free-floating sections.

ChAT was detected using a polyclonal goat anti-ChAT antibody (AB144P; Millipore; 1:100) and EYFP was detected using a polyclonal mouse anti-GFP (06–896, Millipore; 1:1000) diluted in PB containing 2% normal donkey serum (017-000-121, Jackson Immuno Research Laboratories) and 0.3% Triton X-100, incubated for 24 h. Sections were subsequently rinsed in PB and incubated for 2 h in an Alexa 488 donkey anti-mouse antibody (715-545-150; 1:250; Jackson Immuno Research Laboratories) and Alexa 594 donkey anti-goat antibody (705-585-003; 1:250; Jackson Immuno Research Laboratories). For all secondary antibodies, control experiments confirmed that no labeling was observed when primary antibodies were omitted. The sections were mounted on slides in rostrocaudal sequential order, dried, and covered with fluoromount (00-4,958-02; Thermo Fisher). Coverslips were affixed with nail polish.

2.7 Cell counting, imaging, and data analysis

A VS120-S6-W Virtual Slide Scanner (Olympus; Tokyo, Japan) was used to scan all the sections. Images were taken with a color camera (Nikon DS-Fi3; Tokyo, Japan). To reduce bias, the photomicrography and counting were performed by a blinded researcher. Image J (version 1.41; National Institutes of Health, Bethesda, MD) was used for cell counting and Canvas software (ACD Systems, Victoria, Canada, v. 9.0) was used for line drawings. A one-in-two series of 25-μm brain sections was used per mouse, such that each section analyzed was 50 μm apart. The area analyzed was delimited (mean of 5,423 μm²) based on previous data (Anderson et al., 2016). Sections were examined to confirm the number of transfected cells. The number of labeled cells were bilaterally counted and reported as median ± interquartile range (IQR). Section alignment was relative to a reference section, as previously described (Anderson et al., 2016) and based on the Paxinos and Franklin atlas (Paxinos and Franklin, 2012).

2.8 Statistical analysis

Statistical analyses were performed using GraphPad Prism 8 (GraphPad, San Diego, CA). All data passed the Shapiro–Wilk test for normality. Two-way analysis of variance (ANOVA) followed by Bonferroni post hoc tests were used to test for phase-relationship responses. Unpaired, two-tailed t-tests were used to compare the baseline activity between control and CIH groups. The comparisons between baseline and optogenetic stimulation within the same group were made via repeated measures, two-way ANOVA or a mixed effects analysis (missing data). Differences were considered significant at p < 0.05. Investigators were not blinded during the analyses. Sample sizes were chosen based on previous studies (Anderson et al., 2016; Huff et al., 2022, 2023). Raw data summaries are presented in Tables 1–5. ANOVA tables are presented in Supplementary Table 1.

Table 1

Table 1. Phase-dependent sympathetic parameters during baseline conditions and during optogenetic stimulation of cholinergic-glutamatergic PiCo neurons in control and CIH mice.

Table 2

Table 2. Phase-dependent sympathetic parameters during baseline conditions and during optogenetic stimulation of glutamatergic PiCo neurons in control and CIH mice.

Table 3

Table 3. Phase-dependent sympathetic parameters during baseline conditions and during optogenetic stimulation of cholinergic PiCo neurons in control and CIH mice.

Table 4

Table 4. Sympathetic parameters during baseline conditions and during optogenetic stimulation of Dbx1 preBötC neurons in control and CIH mice.

Table 5

Table 5. Sympathetic parameters during baseline conditions and optogenetic stimulation of glutamatergic preBötC neurons in control and CIH mice.

3 Results

3.1 Stimulating cholinergic-glutamatergic neurons at PiCo increases sympathetic activity in a respiratory phase-dependent manner

Previous studies have established that PiCo contains neurons that co-express the genetic markers for ChAT and Vglut2 (Anderson et al., 2016; Huff et al., 2022, 2023). To evaluate the role of PiCo specific neurons in respiratory-sympathetic coupling, neurons co-expressing ChAT and Vglut2 were targeted for optogenetic stimulation. Figure 1 shows neurons from transfected cells located just dorsal to the nucleus ambiguus near −6.84 mm posterior to Bregma, ~1.1 mm from the midline, and ~ 0.7 mm above the marginal layer according to the Paxinos and Franklin mouse atlas (2012). Neurons were targeted using an intersectional genetic strategy via bilateral injections of the vector AAV8-ConFon-ChR2-EYFP into PiCo region of ChATcre:Vglut2FlpO mice (Figure 1A₁) left, (Figure 1A₂) right. These cells were all ChAT-immunoreactive (ChAT-ir, red) and EYFP positive (green), indicating successful expression of channelrhodopsin. Nucleus ambiguus cholinergic neurons were not transfected (red, Figures 1A_1,2). The median number of transfected neurons in both hemispheres was not significantly different between the control and CIH groups (Figure 1A₃). Repeated measures two-way ANOVA revealed that in-vivo optogenetic stimulation (blue vertical bar) of ChAT-Vglut2 PiCo neurons produced significant (F(1, 4) = 9.981; p = 0.0342) increases in cSN activity (green trace vs. gray baseline trace) in a respiratory phase-dependent manner (F(2, 8) = 34.06; p < 0.0001). Bonferroni post hoc multiple comparisons test confirmed that optogenetic stimulation evoked a significant increase in sympathetic activity during inspiration (p < 0.0001; Figure 1B₁), postinspiration (p < 0.0001; Figure 1B₂), and late expiration (p < 0.0001; Figure 1B₃). Heart rate was unaffected (Table 1).

Figure 1

Figure 1. Stimulating cholinergic-glutamatergic neurons at PiCo increases sympathetic activity in-vivo. Representative transverse hemisections in a ChATcre:Vglut2FlpO:ChR2 mouse at −6.8 mm posterior to bregma illustrates selective bilateral transfection of cholinergic-glutamatergic neurons of PiCo left (A₁) and right (A₂) via the pAAv-hSyn Con/Fon hChR2 (H134R)-EYFP vector. Violin plots denote median (solid line) and interquartile range (IQR; dashed lines) and individual mouse data is illustrated via dots. Data distributions are denoted by the edges of the violins. Unpaired, two-tailed t-test indicated number of transfected cells was not different between control (CON) and chronic intermittent hypoxia (CIH) groups. Representative traces of raw and integrated (∫) diaphragm (DIA), cervical sympathetic nerve (cSN), and vagus nerve (cVN) are shown for a control ChATcre:Vglut2FlpO:ChR2 mouse during inspiration (1), postinspiration (2), and late expiration (3). The vertical shaded blue area in the trace shows the 200 ms LED activation, and the green trace illustrates the response to the optogenetic stimulation relative to the gray baseline trace within the same representative mouse. Repeated measures, two-way ANOVA with a Bonferroni post hoc test revealed optogenetic stimulation of cholinergic-glutamatergic neurons at PiCo caused significant increases in sympathetic nerve activity during inspiration (B₁), postinspiration (B₂), and late expiration phases (B₃). *p < 0.05; n = 5 (2 male, 3 female).

3.2 Optogenetic stimulation of glutamatergic, but not cholinergic, neurons at PiCo causes phase-dependent increases in sympathetic activity

To compare the contribution of glutamatergic versus cholinergic neurons in the PiCo area, we optogenetically stimulated Vglut2 or ChAT containing neurons in control mice at random times in the respiratory cycle using previously determined landmarks (Anderson et al., 2016). Figure 2 shows representative data during the in-vivo LED stimulations (blue vertical bar) with the colored trace compared to the baseline trace (gray). Repeated measures, two-way ANOVA showed sympathetic responses to optogenetic stimulation (F(1, 5) = 11.31; p = 0.02) of glutamatergic neurons at PiCo were respiratory phase-dependent (F(2, 10) = 4.382; p = 0.043). Bonferroni post hoc test confirmed that stimulating Vglut2cre:Ai32 neurons at PiCo in control mice evoked significant increases in sympathetic activity during postinspiration (p = 0.0035; Figure 2A₂), but not during late expiration (Figure 2A₃) or inspiration (Figure 2A₁). This resulted in significantly decreased discharge duration (Table 2) but did not affect heart rate. By contrast, stimulation of cholinergic neurons did not cause a significant increase in sympathetic activity during any phase of the respiratory cycle (Figure 2B₁–B₃; Table 3).

Figure 2

Figure 2. Stimulation of glutamatergic, but not cholinergic, neurons at PiCo causes phase-dependent increases in sympathetic activity. Each panel contains representative traces of raw and integrated (∫) diaphragm (DIA), cervical sympathetic nerve (cSN), and vagus nerve (cVN) during baseline (gray) and optogenetic stimulation (colored) from a control Vglut2cre:Ai32 mouse (A), orange; n = 6 (2 male, 4 female) or a control ChATcre:Ai32 mouse (B), blue; n = 6 (4 male, 2 female) during each respiratory phase. The shaded light blue area in the trace represents the 200 ms LED activation at PiCo during inspiration (1), postinspiration (2), and late expiration (3). Violin plots illustrate median (solid line) and interquartile range (IQR; dashed lines) while each dot represents data from one mouse. The edges correspond to the data distributions. Repeated measures, two-way ANOVA with a Bonferroni post hoc test revealed optogenetic stimulation of glutamatergic neurons (orange trace) at PiCo evoked significant increases in sympathetic activity when the stimulation occurred during postinspiration (A₂), but not inspiration (A₁) or late expiration (A₃). Optogenetic stimulation of cholinergic neurons (blue trace) at PiCo did not evoke significant increases in sympathetic activity (B₁–B₃). *p < 0.05.

3.3 CIH increases respiratory rate and sympathetic discharge in mice

This study utilized CIH to model the cardiopulmonary effects of obstructive sleep apnea in mice. Figure 3 shows representative traces and summary data from mice that underwent control or CIH conditions. Figures 3A_1,2 have the respiratory cycle delineated into inspiration (I, orange), postinspiration (Post-I, gray), and late expiration (E2, green) to illustrate how the data were analyzed. The arrows in the CIH representative trace are pointing to the increased sympathetic nerve activity. Both control and CIH mice presented eupneic-like respiratory-sympathetic patterns as shown in representative traces in Figures 3A_1,2. Mice that underwent the CIH protocol exhibited higher respiratory rate (RR; p = 0.025; Figure 3B) and a reduced expiratory duration (p = 0.0130; Figure 3C). Furthermore, Figure 3D (p = 0.004) and 3E (p = 0.001) illustrates that cSN activity was significantly increased across the respiratory cycle following CIH. Compared to control mice, CIH mice exhibited significantly elevated sympathetic activity during inspiration (p = 0.0037; Figure 3F), postinspiration (p = 0.0019; Figure 3G), and late expiration (p = 0.0020; Figure 3H). This was confirmed by the increased postinspiration to inspiration cSN spike ratio (F(1, 64) = 83.09; p < 0.05) via repeated measures two-way ANOVA and a Bonferroni post hoc test (p = 0.039; Figure 3I).

Figure 3

Figure 3. CIH increases respiratory rate and sympathetic discharge in mice. Representative traces of raw and integrated (∫) diaphragm (DIA), cervical sympathetic nerve (cSN), and vagus nerve (cVN) from a urethane anesthetized mouse after control (CON, n = 32, A₁) or chronic intermittent hypoxia (CIH, n = 34, A₂) protocols. The phases of the respiratory cycle were determined via the DIA and cVN activity. The orange shaded area in the recording represents inspiration (I). The gray shaded area represents postinspiration (PI). The green shaded area represents the late expiratory phase of the respiratory cycle (E2). The spike distribution is plotted in the histogram and the arrows are pointing to the increased sympathetic nerve activity after CIH. Violin plots show median (solid line) and interquartile range (IQR; dashed lines) while the edges follow the distribution of the data collected from mice that received normoxia or CIH. Unpaired, two-tailed t-tests revealed CIH significantly increased respiratory rate (B) and decreased expiratory duration (C), while increasing overall cSN voltage (D) and firing rate (E). CIH significantly increased number of cSN spikes during inspiration (F), postinspiration (G), and late expiratory phase (H). Repeated measures, two-way ANOVA revealed the ratio of cSN spikes during postinspiration and inspiration was significantly different while the ratio of spikes during postinspiration and late expiratory phase was not (I). *p < 0.05; **p < 0.01 different from the control group. CON (n = 32, 16 male, 16 female); CIH (n = 34, 18 male, 16 female).

3.4 Optogenetic stimulation of glutamatergic neurons at PiCo causes increases in sympathetic activity following CIH exposure

Optogenetic stimulation of excitatory neurons at PiCo was repeated in mice that were exposed to CIH to investigate how CIH affects the role of PiCo in respiratory-sympathetic coupling. Figure 4 illustrates the diaphragm and sympathetic responses to LED activation of excitatory neurons at PiCo area in CIH mice. Following CIH exposure, mixed effects analysis and Bonferroni post hoc tests revealed optogenetic stimulation of neurons co-expressing ChAT and Vglut2 increased (F(1, 4) = 14.41; p = 0.019) the peak cSN amplitude when the stimulus occurred during postinspiration (p = 0.021; 4A₂) and late expiration (p = 0.006; 4A₃), but not inspiration (p = 0.318; 4A₁). The area under the curve and peak amplitudes of the cSN were also significantly increased in response to the LED activation during late expiration (Table 1). Heart rate remained unchanged (Table 1). Repeated measures, two-way ANOVA showed optogenetic stimulation of glutamatergic PiCo neurons evoked significant (F(1, 5) = 7.911; p = 0.037) increases in cSN activity only during postinspiration (p = 0.019; Figure 4B₂) without affecting heart rate (Table 2). Similar to the control mice (Figure 3B), stimulation of cholinergic PiCo neurons in CIH mice did not significantly increase cSN activity during any phase of the respiratory cycle (Figure 4C₁–C₃; Table 3). The maximum amplitudes of cSN responses to LED stimulation of cholinergic, glutamatergic, and cholinergic-glutamatergic neurons were not significantly different between the control and CIH groups.

Figure 4

Figure 4. Stimulation of glutamatergic neurons at PiCo causes increases in sympathetic activity after CIH exposure. For each panel, traces of raw and integrated (∫) diaphragm (DIA), cervical sympathetic nerve (cSN), and vagus nerve (cVN) are shown for one ChATcre:Vglut2FlpO:ChR2 (A); n = 5 (2 male, 3 female), Vglut2cre:Ai32 (B); n = 6 (1 male, 5 female), or ChATcre:Ai32 (C); n = 5 (5 male) mouse that underwent the CIH protocol. The 200 ms LED triggered during inspiration (1), postinspiration (2), and late expiration (3) is illustrated by the vertical shaded blue area, with the colored traces representing the response to the optogenetic stimulation relative to the gray baseline trace. Group data are displayed via violin plots with median (solid line), interquartile range (IQR; dashed lines), and data distributions (violin edges). Each dot represents data from one mouse. Mixed effects analysis and Bonferroni post hoc test showed optogenetic stimulation of neurons that co-express ChAT and Vglut2 at PiCo (red trace) evoked significant increases in sympathetic activity when the stimulation occurred during postinspiration (A₂) and late expiration (A₃), but not when it occurred during inspiration (A₁). Optogenetic stimulation of glutamatergic neurons of Vglut2cre:Ai32 mice at PiCo (blue trace) evoked significant increases in sympathetic activity only during postinspiration (B₂) as shown via repeated measures, two-way ANOVA and Bonferroni post hoc test. No significant differences in sympathetic activity were observed following optogenetic stimulations in the ChATcre:Ai32 mice (pink trace) exposed to CIH (C₁–C₃). *p < 0.05.

3.5 Optogenetic stimulation of Dbx1 neurons at preBötC evokes inspiration and sympathetic nerve activity

Since optogenetic stimulation of excitatory preBötC neurons always evoked an inspiration, it was impossible to stimulate randomly across the respiratory cycle as was done during the PiCo stimulations. Representative traces show that LED activation (blue vertical bar) of Dbx1cre:Ai32 neurons at preBötC increased sympathetic activity in both control mice (Figure 5A) and mice exposed to CIH (Figure 5B). Optogenetic stimulations of 200 ms increased the area under the curve of the integrated cSN activity (p = 0.0088; Figure 5A₁) in control mice. It also increased the maximum amplitude (p = 0.0004; Figure 5A₂) and absolute peak amplitude (p = 0.0336; Table 4) without changing cSN discharge duration (Table 4). In mice exposed to CIH, stimulation by LED of Dbx1 neurons significantly increased sympathetic activity (p = 0.0044; Figure 5B₁), maximum amplitude (p = 0.0031; Figure 5B₂), and absolute peak amplitude (p = 0.0009; Table 4). No change in cSN discharge duration was observed (Table 4). The overall median cSN responses to optogenetic stimulation were not different between the control and CIH groups (Figure 5C), however, the maximum peak amplitude of cSN following Dbx1 neuron activation at preBötC was increased in mice that were exposed to CIH (p = 0.0153; Figure 5D). We observed similar baseline levels in heart rate between control and CIH groups, and optogenetic stimulation did not affect heart rate (Table 4).

Figure 5

Figure 5. Optogenetic stimulation of Dbx1 neurons at preBötC evokes inspiration and sympathetic nerve activity. Representative data traces of raw and integrated (∫) diaphragm (DIA), cervical sympathetic (cSN), and vagus (cVN) nerve recordings from in-vivo preparations of a Dbx1 mouse maintained under normoxia CON, (A); n = 7 (6 male, 1 female) and exposed to chronic intermittent hypoxia CIH, (B); n = 6 (1 male, 5 female). The shaded blue area represents the 200 ms optogenetic stimulation. Violin plots illustrate median (solid line) and interquartile range (IQR; dashed lines) while each dot represents data from one mouse. Distribution of the data is shown by the plot edges. Optogenetic stimulation evoked significant increases in sympathetic activity for control (A₁,A₂) and CIH mice (B₁,B₂) as confirmed by paired, two-tailed t-tests. The change in sympathetic activity caused by optogenetic stimulation was not significantly different between control and CIH mice (C), however, the percent max amplitude was higher in the CIH group compared to controls (D). *p < 0.05; **p < 0.01; ***p < 0.0001 different from the control group.

3.6 Optogenetic stimulation of glutamatergic neurons at preBötC evokes inspiration and sympathetic nerve activity

Figure 6 displays responses to optogenetic stimulation (vertical blue bar) of Vglut2cre:Ai32 neurons at preBötC in control mice (Figure 6A) and CIH mice (Figure 6B). In control mice, LED stimulation significantly increased the cSN area under the curve (p = 0.0359; Figure 6A₁), maximum amplitude (p = 0.0032; Figure 6A₂), and absolute peak amplitude (p = 0.0190; Table 2) without altering cSN discharge duration (Table 5). After CIH exposure, stimulation of glutamatergic neurons also significantly increased sympathetic activity (p = 0.0064; Figure 6B₁), cSN maximum amplitude (p = 0.0005; Figure 6B₂) and absolute peak amplitude (p = 0.0399; Table 5). Duration of cSN discharge and heart rate were both unaffected (Table 5). Similar to Dbx1 stimulation, the maximal peak amplitude responses to Vglut2 neuronal stimulation were significantly higher in CIH mice compared to control mice (p = 0.0357; Figure 6D), although the magnitude of the responses were not (Figure 6C).

Figure 6

Figure 6. Optogenetic activation of glutamatergic neurons at preBötC evokes inspiration and sympathetic nerve activity. Representative raw and integrated (∫) diaphragm (DIA), cervical sympathetic (cSN), and vagus (cVN) nerve traces recorded from a control CON, (A); n = 8 (5 male, 3 female) and chronic intermittent hypoxia (CIH) Vglut2cre:Ai32 mouse (B); n = 7 (3 male, 4 female). Optogenetic stimulations of 200 ms are shown via the blue shaded region. Violin plots indicate median (solid line) and interquartile range (IQR; dashed lines) and individual mouse data are indicated by the dots. Distribution of the data is shown by the edges of the plots. Paired, two-tailed t-tests show that optogenetic stimulation caused significant increases in sympathetic activity for control (A₁,A₂) and CIH mice (B₁,B₂). The optogenetic evoked changes in sympathetic activity were not significantly different between control and CIH mice (C), however, the percent max amplitude of the responses was higher in the CIH group compared to the control group (D). *p < 0.05; **p < 0.01; ***p < 0.001 different from the control group.

3.7 Optogenetic stimulation of inhibitory neurons at preBötC decreases respiratory rate but does not alter sympathetic activity

To test for preBötC inhibitory control over sympathetic output, this study optogenetically stimulated Vgatcre:Ai32 neurons for 10 s at preBötC area. Figure S1 shows representative data from control (Supplementary Figure S1A) and CIH (Supplementary Figure S1C) mice. LED activation of inhibitory neurons decreased respiratory rate without altering the sympathetic discharge in both control (Supplementary Figure S1B) and CIH groups (Supplementary Figure S1D).

4 Discussion

Exposure to CIH altered the respiratory modulation of sympathetic discharge in urethane anesthetized, freely breathing mice. Mice exposed to CIH exhibited abnormal breathing patterns, disrupted respiratory-sympathetic coupling, and significantly elevated sympathetic activity. Here we provide new evidence that PiCo contributes to respiratory-sympathetic coupling in control mice and mice exposed to CIH. Optogenetic activation of glutamatergic-cholinergic neurons in PiCo region evoked increases in sympathetic activity throughout the respiratory cycle. For CIH mice, sympathetic activity was increased in response to LED stimulation of glutamatergic-cholinergic neurons during postinspiration and late expiration phases. Thus, optogenetic activation of PiCo specific neurons increased sympathetic activity in both control and CIH mice. To differentiate the effects of stimulating glutamatergic versus cholinergic neurons within PiCo, we found that stimulating glutamatergic neurons in PiCo evoked increases in sympathetic activity only during the postinspiratory phase in control mice. Following CIH, optogenetic stimulation of glutamatergic PiCo neurons also only evoked significant increases in sympathetic activity during the postinspiratory phase. Interestingly, stimulating only cholinergic PiCo neurons did not evoke an elevated sympathetic response compared to baseline in control or CIH mice. Stimulating Dbx1 or glutamatergic neurons in preBötC of control and CIH mice evoked inspiratory activity concomitant with increased sympathetic activity. Overall, sympathetic activity was differentially increased by phasic inputs from optogenetically stimulated excitatory respiratory neurons in the ventrolateral medulla. These findings, along with limitations of this study, will be discussed.

4.1 Chronic intermittent hypoxia alters respiratory-sympathetic coupling

CIH models the frequent hypoxic episodes that occur in humans with obstructive sleep apnea (Barnes et al., 2022) and other disorders of dysautonomia (Ramirez et al., 2023). In the present study, mice that underwent the CIH protocol exhibited higher baseline respiratory rate with reduced expiratory duration, as well as increased cervical sympathetic nerve activity (Figure 1). Furthermore, this sympathetic overactivity occurred throughout the respiratory cycle, indicating a decoupling of normal respiratory-sympathetic modulation (Figure 1). These results are consistent with the CIH-induced cardiopulmonary changes observed in rats (Zoccal et al., 2008; Molkov et al., 2011; Moraes et al., 2012; Karlen-Amarante et al., 2023) and increases in blood pressure seen in mice (Campen et al., 2005). These findings support the interpretation that CIH secondary to obstructive sleep apnea is a risk factor for developing high blood pressure due to increased sympathetic drive.

The mechanisms underlying the disruption of normal respiratory-sympathetic coupling following CIH remain unclear. Early life exposure to CIH elicits higher levels of blood pressure and a greater carotid body sensory response during subsequent exposure to hypoxia (Nanduri et al., 2012). This has been attributed to dysregulation of hypoxia-inducible factor (HIF) 1α and 2α which are oxygen sensitive transcriptional factors that regulate the homeostatic responses to hypoxia (Prabhakar and Semenza, 2012). Consistent with the results of the present study, presympathetic neurons in the RVLM of rats exposed to CIH show increased resting activity and firing variability (Karlen-Amarante et al., 2023). This increased firing activity of RVLM presympathetic neurons is not due to changes in the intrinsic electrophysiological properties induced by CIH (Almado et al., 2014). Rather, it is due to a plethora of factors that collectively alter presympathetic neuron activity and disrupt respiratory-sympathetic coupling. These factors include overexpression of HIF-1α (Karlen-Amarante et al., 2023), sensitization of central (Molkov et al., 2011) and peripheral chemoreceptors (Prabhakar et al., 2005, 2023), and dysregulation of astrocyte homeostatic control of blood pressure (Marina et al., 2020). However, more studies are needed to elucidate the complex mechanisms underlying sympathetic overactivation caused by CIH (Guyenet et al., 2020; Barnes et al., 2022).

4.2 Optogenetic activation of excitatory neurons in respiratory networks in-vivo increases sympathetic activity

In addition to sympathetic overactivity, this study also demonstrates that phasic optogenetic stimulation of excitatory neurons in PiCo and preBötC increased sympathetic activity. Mice that received the CIH protocol showed larger inspiratory related peak sympathetic bursts during optogenetic stimulation of Dbx1 (Figure 5D) and Vglut2 (Figure 6D) neurons at preBötC. This suggests that CIH caused an increase in excitatory tone to the RVLM presympathetic neurons from preBötC which altered respiratory-sympathetic modulation. The increased sympathetic activity following stimulation of excitatory neurons in preBötC complements existing evidence that neurons in the preBötC send direct projections to presympathetic C1 neurons (Moraes et al., 2013; Dempsey et al., 2017; Menuet et al., 2017, 2020) where they provide excitatory and inhibitory drive to vasomotor sympathetic activity. This tonic drive likely occurs via presympathetic C1 neurons whose activity positively correlates with inspiration (Marina et al., 2011; Moraes et al., 2013; Menuet et al., 2017). In rats, prolonged optogenetic stimulation of inhibitory preBötC neurons decreases sympathetic activity (Menuet et al., 2020) via suppression of the phasic activity of inhibitory expiratory preBötC neurons (Carroll et al., 2013; Baertsch et al., 2018), indicating the preBötC is critically involved in tonic regulation of the presympathetic neuronal discharge through a combination of excitatory and inhibitory neuronal mechanisms. However, in the present study, stimulation of GABAergic neurons in preBötC of mice did not yield any significant decreases in sympathetic activity compared to baseline (Supplementary Figure S1). This contrasting result may be attributed to differences between species, stimulation durations, or inhibitory neuronal genetic markers targeted.

In addition to preBötC, PiCo may also excite presympathetic neurons via its projections to the RVLM (Dempsey et al., 2017). Stimulating glutamatergic or cholinergic-glutamatergic neurons in the PiCo area had a more pronounced effect on sympathetic activity than stimulating excitatory neurons in preBötC. For example, excitatory preBötC stimulations caused a spike in sympathetic activity and a concurrent inspiration (Figures 5, 6), whereas only stimulating glutamatergic or more specifically, the cholinergic-glutamatergic PiCo neurons caused a disassociation of normal respiratory-sympathetic coupling by evoking a spike in sympathetic activity independent of inspiration (Figures 2–4). Moreover, optogenetic stimulation of cholinergic-glutamatergic PiCo neurons significantly increased sympathetic activity at all phases of the respiratory cycle and not only during the postinspiratory phase. However, since PiCo is typically only active in phase with postinspiration (Anderson et al., 2016), this region likely contributes to the postinspiratory component of respiratory-sympathetic coupling under control conditions. We recently demonstrated via respiratory phase reset analysis that activation of glutamatergic-cholinergic PiCo neurons during inspiration or at the beginning of postinspiration is highly likely to reset the respiratory rhythm and trigger a swallow, suggesting that PiCo plays a role in swallow-breathing coordination (See Huff et al., 2023; Figure 3A). In the present study, similar swallow responses and corresponding cSN activity were observed during optogenetic stimulation of PiCo, however, these responses were not included in the analysis, allowing the characterization of PiCo contribution to cardiorespiratory coupling separately from swallow behavior.

Following CIH, previous studies reported distinct firing patterns of RVLM neuronal activity during different phases of the respiratory cycle (Moraes et al., 2013; Zoccal et al., 2019) and overall increases in RVLM activity throughout the respiratory cycle (Zoccal et al., 2008; Moraes et al., 2013; Zoccal et al., 2019), leading to high blood pressure (Zoccal et al., 2007a,b; Souza et al., 2016; Bazilio et al., 2019). Since optogenetic stimulation of excitatory neurons in both preBötC and PiCo increased sympathetic activity, the activity from these two rhythmogenic networks may differentially modulate sympathetic activity during inspiration and postinspiration, respectively. PiCo is not always active throughout the respiratory cycle (Anderson et al., 2016), and postinspiratory motor activity is often missing in phrenic nerve recording of rodents (Dutschmann and Paton, 2002), which may explain why sympathetic activity is dominated by diaphragm burst activity under baseline conditions. The present data show that PiCo and preBötC both evoke sympathetic burst activity, but further studies are necessary to unravel their underlying connectivity to the RVLM sympathetic neurons. RVLM sympathetic premotor neurons include both C1 and non-C1 neurons, which differ in their functional and neurochemical phenotypes and could play differing roles in the generation of sympathetic nerve activity (Schreihofer and Guyenet, 1997; Stornetta et al., 2001).

4.3 Limitations and conclusions

There are several limitations to the present study that merit consideration. First, there are differences in healthy sympathetic patterns between rodents and humans. Specifically, sympathetic activity dominates during expiration instead of inspiration in humans (Zoccal et al., 2008, 2009; Karlen-Amarante et al., 2023). It is conceivable that, in humans, inspiratory driven sympathetic activity is largely inhibited by vagal afferents, as suggested by experimentally induced lung inflations which inhibit pre-and postganglionic sympathetic nerve activity in anesthetized or decerebrate, vagally intact laboratory animals (Seals et al., 1990). During quiet resting breathing, muscle sympathetic nerve activity (MSNA) is highest at end-expiration and lowest at end-inspiration (Hagbarth and Vallbo, 1968; Katayama et al., 2021). However, the respiratory modulation of MSNA is independent of changes in blood pressure and baroreceptor activity (Seals et al., 1990; Macefield and Wallin, 1995).

Second, measures of respiratory and sympathetic activity contained a high amount of variability which may have increased the likelihood of committing a type-1 error. Due to the nature of the methods by which the data must be collected, a high amount of variability is expected, as mentioned in the “Optogenetic stimulation and data analysis” section. However, every effort was made to quantify differences in dependent measures without unnecessary data transformations. No outliers were removed, and any missing data was missing at random. Despite this variability, significant effects due to optogenetic stimulations or CIH were observed, underscoring the robust nature of the phenomena.

Third, due to technical limitations associated with the exposed ventral brainstem in-vivo preparation, blood pressure was not included as a dependent measure in this study. Despite the absence of blood pressure data, we did not see any changes in the HR associated with the optogenetic stimulations. This is likely because the latency of the cSN response to the optogenetic stimulations was too short to be triggered by periphery circuits, suggesting a central origin of the sympathetic response. Previous studies have correlated sympathetic overactivity with the incidence of hypertension following CIH exposure in rodents (Greenberg et al., 1985; Fletcher, 2001; Zoccal et al., 2007b; Prabhakar and Kumar, 2010; Moraes et al., 2012) and humans (Somers et al., 1995), and these findings contribute to the interpretation of the present data. Most studies of sympathetic activity are performed using the thoracic, renal, or splanchnic sympathetic branches which all present a phasic activity modulation that increases during the end of inspiration and decreases throughout the respiratory cycle. Although the cSN activity in the present study has a monophasic discharge, this difference is likely due to the differential channeling of A- and C-fiber inputs to distinct populations of vasomotor neurons (McMullan et al., 2008).

Forth, it is a well-known limitation of the optogenetic technique that collateral activation of neurons outside of the target area may occur due to light scatter through tissue. It is possible other nearby glutamatergic neurons are stimulated when targeting the PiCo or preBötC areas, since glutamatergic neurons are distributed throughout the entire VRC (Bush and Ramirez, 2024). It is also an issue when stimulating only cholinergic PiCo neurons. Adjacent cholinergic neurons located in the Nucleus ambiguus will also be stimulated during stimulation of the cholinergic-glutamatergic interneurons in the PiCo region. Precautions were taken to circumvent this possibility such as titrating the LED power and using established, previously published landmarks for the position of the LED optrodes (Anderson et al., 2016; Huff et al., 2023). Furthermore, stimulations of cholinergic-glutamatergic neurons were performed to specifically limit stimulation only to PiCo neurons, as only PiCo neurons have been shown to express both ChAT and Vglut2 excitatory markers in that brainstem region (Anderson et al., 2016; Huff et al., 2023). PiCo is a relatively new region of interest, and future studies are necessary to fully elucidate the exact function, nomenclature, anatomical connections, and contribution of these neurons to respiratory rhythm, swallowing, and laryngeal behaviors as well other autonomic functions.

This study demonstrates that PiCo and preBötC contribute to respiratory-sympathetic coupling and that exposure to CIH alters this coupling. CIH exposure disrupted breathing, increased baseline sympathetic tone, and altered respiratory phase-dependent modulation of sympathetic activity. Optogenetically stimulating preBötC and PiCo glutamatergic neurons and PiCo specific cholinergic-glutamatergic neurons significantly increased sympathetic activity. Taken together, the present study suggests that the heightened sympathetic activity to the cardiovascular system in mice exposed to CIH may be driven by increased excitatory tone from the respiratory networks. Specifically, the current results support the interpretation that CIH causes a cascade of events which begin with an overactive carotid body (Nanduri et al., 2012) and lead to RLVM hyperactivity that is exacerbated by an increased excitatory drive from PiCo and preBötC. However, it is unlikely that these respiratory complexes are the only contributing factors to the RVLM overactivity, given that CIH is known to cause a variety of other changes in the respiratory network (Zoccal et al., 2008; Molkov et al., 2011; Moraes et al., 2012, 2013; Bazilio et al., 2019; Zoccal et al., 2019; Bittencourt-Silva et al., 2020; Karlen-Amarante et al., 2023). These findings are consistent with the pathogenesis of obstructive sleep apnea related hypertension, which is a risk factor for developing life-threatening cardiovascular diseases.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The animal study was approved by Seattle Childrens Research Institute Institutional Animal Care and Use Committee. The study was conducted in accordance with the local legislation and institutional requirements.

Author contributions

MK-A: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Software, Validation, Writing – original draft, Writing – review & editing. ZG: Conceptualization, Data curation, Formal analysis, Methodology, Software, Validation, Visualization, Writing – original draft, Writing – review & editing. AH: Conceptualization, Data curation, Investigation, Methodology, Software, Writing – original draft, Writing – review & editing. LO: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Software, Writing – original draft, Writing – review & editing. J-MR: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research was funded by NIH grants P01 HL144454 and Project 2 (awarded to J-MR) R01 HL144801 (awarded to J-MR), R01 HL 151389 (awarded to J-MR), F32 HL160102-01 (awarded to AH).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2024.1386737/full#supplementary-material

References

Almado, C. E., Leao, R. M., and Machado, B. H. (2014). Intrinsic properties of rostral ventrolateral medulla presympathetic and bulbospinal respiratory neurons of juvenile rats are not affected by chronic intermittent hypoxia. Exp. Physiol. 99, 937–950. doi: 10.1113/expphysiol.2013.077800

PubMed Abstract | Crossref Full Text | Google Scholar

Anderson, T. M., Garcia, A. J. 3rd, Baertsch, N. A., Pollak, J., Bloom, J. C., Wei, A. D., et al. (2016). A novel excitatory network for the control of breathing. Nature 536, 76–80. doi: 10.1038/nature18944

PubMed Abstract | Crossref Full Text | Google Scholar

Anderson, T. M., and Ramirez, J. M. (2017). Respiratory rhythm generation: triple oscillator hypothesis. F1000Res 6:139. doi: 10.12688/f1000research.10193.1

PubMed Abstract | Crossref Full Text | Google Scholar

Baertsch, N. A., Baertsch, H. C., and Ramirez, J. M. (2018). The interdependence of excitation and inhibition for the control of dynamic breathing rhythms. Nat. Commun. 9:843. doi: 10.1038/s41467-018-03223-x

PubMed Abstract | Crossref Full Text | Google Scholar

Baertsch, N. A., Severs, L. J., Anderson, T. M., and Ramirez, J. M. (2019). A spatially dynamic network underlies the generation of inspiratory behaviors. Proc. Natl. Acad. Sci. USA 116, 7493–7502. doi: 10.1073/pnas.1900523116

PubMed Abstract | Crossref Full Text | Google Scholar

Barnes, L. A., Mesarwi, O. A., and Sanchez-Azofra, A. (2022). The cardiovascular and metabolic effects of chronic hypoxia in animal models: a Mini-review. Front. Physiol. 13:873522. doi: 10.3389/fphys.2022.873522

PubMed Abstract | Crossref Full Text | Google Scholar

Bazilio, D. S., Bonagamba, L. G. H., Moraes, D. J. A., and Machado, B. H. (2019). Cardiovascular and respiratory profiles during the sleep-wake cycle of rats previously submitted to chronic intermittent hypoxia. Exp. Physiol. 104, 1408–1419. doi: 10.1113/EP087784

PubMed Abstract | Crossref Full Text | Google Scholar

Bittencourt-Silva, P. G., Menezes, M. F., Mendonca-Junior, B. A., Karlen-Amarante, M., and Zoccal, D. B. (2020). Postnatal intermittent hypoxia enhances phrenic and reduces vagal upper airway motor activities in rats by epigenetic mechanisms. Exp. Physiol. 105, 148–159. doi: 10.1113/EP087928

PubMed Abstract | Crossref Full Text | Google Scholar

Bouvier, J., Thoby-Brisson, M., Renier, N., Dubreuil, V., Ericson, J., Champagnat, J., et al. (2010). Hindbrain interneurons and axon guidance signaling critical for breathing. Nat. Neurosci. 13, 1066–1074. doi: 10.1038/nn.2622

PubMed Abstract | Crossref Full Text | Google Scholar

Bush, N. E., and Ramirez, J. M. (2024). Latent neural population dynamics underlying breathing, opioid-induced respiratory depression and gasping. Nat. Neurosci. 27, 259–271. doi: 10.1038/s41593-023-01520-3

PubMed Abstract | Crossref Full Text | Google Scholar

Campen, M. J., Shimoda, L. A., and O’Donnell, C. P. (2005). Acute and chronic cardiovascular effects of intermittent hypoxia in C57BL/6J mice. J. Appl. Physiol. 99, 2028–2035. doi: 10.1152/japplphysiol.00411.2005

Crossref Full Text | Google Scholar

Carroll, M. S., Viemari, J. C., and Ramirez, J. M. (2013). Patterns of inspiratory phase-dependent activity in the in vitro respiratory network. J Neurophysiol 109, 285–295.

Google Scholar

Cui, Y., Kam, K., Sherman, D., Janczewski, W. A., Zheng, Y., and Feldman, J. L. (2016). Defining preBotzinger complex rhythm- and pattern-generating neural microcircuits in vivo. Neuron 91, 602–614. doi: 10.1016/j.neuron.2016.07.003

PubMed Abstract | Crossref Full Text | Google Scholar

Dempsey, B., Le, S., Turner, A., Bokiniec, P., Ramadas, R., Bjaalie, J. G., et al. (2017). Mapping and analysis of the connectome of sympathetic premotor neurons in the rostral ventrolateral medulla of the rat using a volumetric brain atlas. Front Neural Circuits 11:9. doi: 10.3389/fncir.2017.00009

Crossref Full Text | Google Scholar

Dumitrascu, R., Heitmann, J., Seeger, W., Weissmann, N., and Schulz, R. (2013). Obstructive sleep apnea, oxidative stress and cardiovascular disease: lessons from animal studies. Oxidative Med. Cell. Longev. 2013:234631, 1–7. doi: 10.1155/2013/234631

Crossref Full Text | Google Scholar

Dutschmann, M., and Paton, J. F. (2002). Glycinergic inhibition is essential for co-ordinating cranial and spinal respiratory motor outputs in the neonatal rat. J. Physiol. 543, 643–653. doi: 10.1113/jphysiol.2001.013466

PubMed Abstract | Crossref Full Text | Google Scholar

Fletcher, E. C. (2001). Invited review: physiological consequences of intermittent hypoxia: systemic blood pressure. J. Appl. Physiol. 90, 1600–1605. doi: 10.1152/jappl.2001.90.4.1600

Crossref Full Text | Google Scholar

Franklin, K. A., and Lindberg, E. (2015). Obstructive sleep apnea is a common disorder in the population—a review on the epidemiology of sleep apnea. J. Thorac. Dis. 7:1311. doi: 10.3978/j.issn.2072-1439.2015.06.11

Crossref Full Text | Google Scholar

Gray, P. A., Hayes, J. A., Ling, G. Y., Llona, I., Tupal, S., Picardo, M. C., et al. (2010). Developmental origin of preBotzinger complex respiratory neurons. J. Neurosci. 30, 14883–14895. doi: 10.1523/JNEUROSCI.4031-10.2010

PubMed Abstract | Crossref Full Text | Google Scholar

Greenberg, H. E., Sica, A., Batson, D., and Scharf, S. M. (1985). Chronic intermittent hypoxia increases sympathetic responsiveness to hypoxia and hypercapnia. J. Appl. Physiol. 86, 298–305. doi: 10.1152/jappl.1999.86.1.298

Crossref Full Text | Google Scholar

Guyenet, P. G., Stornetta, R. L., Souza, G., Abbott, S. B. G., and Brooks, V. L. (2020). Neuronal networks in hypertension: recent advances. Hypertension 76, 300–311. doi: 10.1161/HYPERTENSIONAHA.120.14521

PubMed Abstract | Crossref Full Text | Google Scholar

Hagbarth, K. E., and Vallbo, A. B. (1968). Pulse and respiratory grouping of sympathetic impulses in human muscle-nerves. Acta Physiol. Scand. 74, 96–108. doi: 10.1111/j.1365-201X.1968.tb10904.x

PubMed Abstract | Crossref Full Text | Google Scholar

Huff, A., Karlen-Amarante, M., Oliveira, L. M., and Ramirez, J. M. (2023). Role of the postinspiratory complex in regulating swallow-breathing coordination and other laryngeal behaviors. eLife 12:e86103. doi: 10.7554/eLife.86103

PubMed Abstract | Crossref Full Text | Google Scholar

Huff, A., Karlen-Amarante, M., Pitts, T., and Ramirez, J. M. (2022). Optogenetic stimulation of pre-Botzinger complex reveals novel circuit interactions in swallowing-breathing coordination. Proc. Natl. Acad. Sci. USA 119:e2121095119. doi: 10.1073/pnas.2121095119

Crossref Full Text | Google Scholar

Karlen-Amarante, M., Ip, I., Katayama, P. L., Colombari, E., Bittencourt-Silva, P. G., Menezes, M. F., et al. (2023). Sympathetic dysregulation induced by postnatal intermittent hypoxia. Sleep 46:zsad055. doi: 10.1093/sleep/zsad055

PubMed Abstract | Crossref Full Text | Google Scholar

Katayama, K., Dominelli, P. B., Foster, G. E., Kipp, S., Leahy, M. G., Ishida, K., et al. (2021). Respiratory modulation of sympathetic vasomotor outflow during graded leg cycling. J. Appl. Physiol. 131, 858–867. doi: 10.1152/japplphysiol.00118.2021

Crossref Full Text | Google Scholar

Lipski, J., Kanjhan, R., Kruszewska, B., and Rong, W. (1996). Properties of presympathetic neurones in the rostral ventrolateral medulla in the rat: an intracellular study in vivo. J. Physiol. 490, 729–744. doi: 10.1113/jphysiol.1996.sp021181

PubMed Abstract | Crossref Full Text | Google Scholar

Macefield, V. G., and Wallin, B. G. (1995). Modulation of muscle sympathetic activity during spontaneous and artificial ventilation and apnoea in humans. J. Auton. Nerv. Syst. 53, 137–147. doi: 10.1016/0165-1838(94)00173-H

Crossref Full Text | Google Scholar

Marina, N., Abdala, A. P., Korsak, A., Simms, A. E., Allen, A. M., Paton, J. F., et al. (2011). Control of sympathetic vasomotor tone by catecholaminergic C1 neurones of the rostral ventrolateral medulla oblongata. Cardiovasc. Res. 91, 703–710. doi: 10.1093/cvr/cvr128

Crossref Full Text | Google Scholar

Marina, N., Christie, I. N., Korsak, A., Doronin, M., Brazhe, A., Hosford, P. S., et al. (2020). Astrocytes monitor cerebral perfusion and control systemic circulation to maintain brain blood flow. Nat. Commun. 11:131. doi: 10.1038/s41467-019-13956-y

PubMed Abstract | Crossref Full Text | Google Scholar

McMullan, S., Pathmanandavel, K., Pilowsky, P. M., and Goodchild, A. K. (2008). Somatic nerve stimulation evokes qualitatively different somatosympathetic responses in the cervical and splanchnic sympathetic nerves in the rat. Brain Res. 1217, 139–147. doi: 10.1016/j.brainres.2008.04.034

PubMed Abstract | Crossref Full Text | Google Scholar

Menuet, C., Connelly, A. A., Bassi, J. K., Melo, M. R., Le, S., Kamar, J., et al. (2020). PreBötzinger complex neurons drive respiratory modulation of blood pressure and heart rate. eLife 9:e57288. doi: 10.7554/eLife.57288

PubMed Abstract | Crossref Full Text | Google Scholar

Menuet, C., Le, S., Dempsey, B., Connelly, A. A., Kamar, J. L., Jancovski, N., et al. (2017). Excessive respiratory modulation of blood pressure triggers hypertension. Cell Metab. 25, 739–748. doi: 10.1016/j.cmet.2017.01.019

PubMed Abstract | Crossref Full Text | Google Scholar

Molkov, Y. I., Zoccal, D. B., Moraes, D. J., Paton, J. F., Machado, B. H., and Rybak, I. A. (2011). Intermittent hypoxia-induced sensitization of central chemoreceptors contributes to sympathetic nerve activity during late expiration in rats. J. Neurophysiol. 105, 3080–3091. doi: 10.1152/jn.00070.2011

PubMed Abstract | Crossref Full Text | Google Scholar

Moraes, D. J., da Silva, M. P., Bonagamba, L. G., Mecawi, A. S., Zoccal, D. B., Antunes-Rodrigues, J., et al. (2013). Electrophysiological properties of rostral ventrolateral medulla presympathetic neurons modulated by the respiratory network in rats. J. Neurosci. 33, 19223–19237. doi: 10.1523/JNEUROSCI.3041-13.2013

Crossref Full Text | Google Scholar

Moraes, D. J., Zoccal, D. B., and Machado, B. H. (2012). Medullary respiratory network drives sympathetic overactivity and hypertension in rats submitted to chronic intermittent hypoxia. Hypertension 60, 1374–1380. doi: 10.1161/HYPERTENSIONAHA.111.189332

PubMed Abstract | Crossref Full Text | Google Scholar

Nanduri, J., Makarenko, V., Reddy, V. D., Yuan, G., Pawar, A., Wang, N., et al. (2012). Epigenetic regulation of hypoxic sensing disrupts cardiorespiratory homeostasis. Proc. Natl. Acad. Sci. USA 109, 2515–2520. doi: 10.1073/pnas.1120600109

Crossref Full Text | Google Scholar

O'Brien, S. J., Menotti-Raymond, M., Murphy, W. J., Nash, W. G., Wienberg, J., Stanyon, R., et al. (1999). The promise of comparative genomics in mammals. Science 286, 458–481. doi: 10.1126/science.286.5439.458

Crossref Full Text | Google Scholar

Paxinos, G., and Franklin, F. (2012). The mouse brain in estereotaxic coordinates. Elsevier, vol. 1. 1st Edn. Sydney: Academic Press. Elsevier.

Google Scholar

Percie du Sert, N., Hurst, V., Ahluwalia, A., Alam, S., Avey, M. T., Baker, M., et al. (2020). The ARRIVE guidelines 2.0: updated guidelines for reporting animal research. J. Cereb. Blood Flow Metab. 40, 1769–1777. doi: 10.1177/0271678X20943823

PubMed Abstract | Crossref Full Text | Google Scholar

Prabhakar, N. R., and Kumar, G. K. (2010). Mechanisms of sympathetic activation and blood pressure elevation by intermittent hypoxia. Respir. Physiol. Neurobiol. 174, 156–161. doi: 10.1016/j.resp.2010.08.021

PubMed Abstract | Crossref Full Text | Google Scholar

Prabhakar, N. R., Peng, Y. J., Jacono, F. J., Kumar, G. K., and Dick, T. E. (2005). Cardiovascular alterations by chronic intermittent hypoxia: importance of carotid body chemoreflexes. Clin. Exp. Pharmacol. Physiol. 32, 447–449. doi: 10.1111/j.1440-1681.2005.04209.x

PubMed Abstract | Crossref Full Text | Google Scholar

Prabhakar, N. R., Peng, Y. J., and Nanduri, J. (2023). Carotid body hypersensitivity in intermittent hypoxia and obtructive sleep apnoea. J. Physiol. 601, 5481–5494. doi: 10.1113/JP284111

Crossref Full Text | Google Scholar

Prabhakar, N. R., and Semenza, G. L. (2012). Adaptive and maladaptive cardiorespiratory responses to continuous and intermittent hypoxia mediated by hypoxia-inducible factors 1 and 2. Physiol. Rev. 92, 967–1003. doi: 10.1152/physrev.00030.2011

PubMed Abstract | Crossref Full Text | Google Scholar

Ramirez, J. M., Carroll, M. S., Burgraff, N., Rand, C. M., and Weese-Mayer, D. E. (2023). A narrative review of the mechanisms and consequences of intermittent hypoxia and the role of advanced analytic techniques in pediatric autonomic disorders. Clin. Auton. Res. 33, 287–300. doi: 10.1007/s10286-023-00958-6

Crossref Full Text | Google Scholar

Schreihofer, A. M., and Guyenet, P. G. (1997). Identification of C1 presympathetic neurons in rat rostral ventrolateral medulla by juxtacellular labeling in vivo. J. Comp. Neurol. 387, 524–536. doi: 10.1002/(SICI)1096-9861(19971103)387:4<524::AID-CNE4>3.0.CO;2-4

Crossref Full Text | Google Scholar

Schulz, R., Seeger, W., Fegbeutel, C., Husken, H., Bodeker, R. H., Tillmanns, H., et al. (2005). Changes in extracranial arteries in obstructive sleep apnoea. Eur. Respir. J. 25, 69–74. doi: 10.1183/09031936.04.00030004

PubMed Abstract | Crossref Full Text | Google Scholar

Seals, D. R., Suwarno, N. O., and Dempsey, J. A. (1990). Influence of lung volume on sympathetic nerve discharge in normal humans. Circ. Res. 67, 130–141. doi: 10.1161/01.RES.67.1.130

PubMed Abstract | Crossref Full Text | Google Scholar

Simms, A. E., Paton, J. F., Pickering, A. E., and Allen, A. M. (2009). Amplified respiratory-sympathetic coupling in the spontaneously hypertensive rat: does it contribute to hypertension? J. Physiol. 587, 597–610. doi: 10.1113/jphysiol.2008.165902

PubMed Abstract | Crossref Full Text | Google Scholar

Somers, V. K., Dyken, M. E., Clary, M. P., and Abboud, F. M. (1995). Sympathetic neural mechanisms in obstructive sleep apnea. J. Clin. Invest. 96, 1897–1904. doi: 10.1172/JCI118235

PubMed Abstract | Crossref Full Text | Google Scholar

Soukhova-O'Hare, G. K., Cheng, Z. J., Roberts, A. M., and Gozal, D. (2006). Postnatal intermittent hypoxia alters baroreflex function in adult rats. Am. J. Physiol. Heart Circ. Physiol. 290, H1157–H1164. doi: 10.1152/ajpheart.00767.2005

Crossref Full Text | Google Scholar

Souza, G. M., Bonagamba, L. G., Amorim, M. R., Moraes, D. J., and Machado, B. H. (2016). Inspiratory modulation of sympathetic activity is increased in female rats exposed to chronic intermittent hypoxia. Exp. Physiol. 101, 1345–1358. doi: 10.1113/EP085850

PubMed Abstract | Crossref Full Text | Google Scholar

Stornetta, R. L., Schreihofer, A. M., Pelaez, N. M., Sevigny, C. P., and Guyenet, P. G. (2001). Preproenkephalin mRNA is expressed by C1 and non-C1 barosensitive bulbospinal neurons in the rostral ventrolateral medulla of the rat. J. Comp. Neurol. 435, 111–126. doi: 10.1002/cne.1196

PubMed Abstract | Crossref Full Text | Google Scholar

Toor, R. U. A. S., Sun, Q.-J., Kumar, N. N., Le, S., Hildreth, C. M., Phillips, J. K., et al. (2019). Neurons in the intermediate reticular nucleus coordinate postinspiratory activity, swallowing, and respiratory-sympathetic coupling in the rat. J. Neurosci. 39, 9757–9766. doi: 10.1523/JNEUROSCI.0502-19.2019

PubMed Abstract | Crossref Full Text | Google Scholar

Zoccal, D. B., Bonagamba, L. G., Antunes-Rodrigues, J., and Machado, B. H. (2007a). Plasma corticosterone levels is elevated in rats submitted to chronic intermittent hypoxia. Auton. Neurosci. 134, 115–117. doi: 10.1016/j.autneu.2007.01.004

PubMed Abstract | Crossref Full Text | Google Scholar

Zoccal, D. B., Bonagamba, L. G., Oliveira, F. R., Antunes-Rodrigues, J., and Machado, B. H. (2007b). Increased sympathetic activity in rats submitted to chronic intermittent hypoxia. Exp. Physiol. 92, 79–85. doi: 10.1113/expphysiol.2006.035501

PubMed Abstract | Crossref Full Text | Google Scholar

Zoccal, D. B., Colombari, D. S. A., Colombari, E., Flor, K. C., da Silva, M. P., Costa-Silva, J. H., et al. (2019). Centrally acting adrenomedullin in the long-term potentiation of sympathetic vasoconstrictor activity induced by intermittent hypoxia in rats. Exp. Physiol. 104, 1371–1383. doi: 10.1113/EP087613

Crossref Full Text | Google Scholar

Zoccal, D. B., Paton, J. F., and Machado, B. H. (2009). Do changes in the coupling between respiratory and sympathetic activities contribute to neurogenic hypertension? Clin. Exp. Pharmacol. Physiol. 36, 1188–1196. doi: 10.1111/j.1440-1681.2009.05202.x

PubMed Abstract | Crossref Full Text | Google Scholar

Zoccal, D. B., Simms, A. E., Bonagamba, L. G., Braga, V. A., Pickering, A. E., Paton, J. F., et al. (2008). Increased sympathetic outflow in juvenile rats submitted to chronic intermittent hypoxia correlates with enhanced expiratory activity. J. Physiol. 586, 3253–3265. doi: 10.1113/jphysiol.2008.154187

PubMed Abstract | Crossref Full Text | Google Scholar