INTRODUCTION {#S5}
============

Bronchopulmonary Dysplasia (BPD) is the major pulmonary morbidity of prematurity, affecting up to 10,000 US infants annually ([@R1]). Increased airway reactivity and wheezing disorders are among the most common long-term pulmonary consequences of prematurity and BPD and worsen with increasing degree of prematurity ([@R2],[@R3]). Studies of former preterm infants demonstrate a combination of increased airway resistance, decreased lung compliance, and/or increased health care utilization (hospitalizations and asthma medications) from term corrected age through adulthood ([@R4]-[@R8]). Physicians in the Neonatal Intensive Care Unit (NICU) are challenged with maintaining adequate levels of O~2~ saturation while avoiding excessive periods of hypoxia and hyperoxia. Recurrent apnea and respiratory insufficiency associated with apnea of prematurity (AOP) lead to frequent episodes of O~2~ desaturation despite supplemental O~2~ therapy ([@R9]). Recent studies in preterm infants using pulse oximetry have well characterized the timing and incidence of these events in the first 2 months of postnatal life ([@R10]). The incidence of desaturations (\<80% saturation for 10-180 seconds) increases dramatically during the first 2 postnatal weeks followed by a steady decline thereafter with extremely preterm infants experiencing over 100 desaturation events per day ([@R11]). In the clinical setting, these events are often resolved by a transient increase in the fraction of inspired oxygen (FIO~2~), most commonly titrated by the bedside nurse. A consequence of increasing FIO~2~ to resolve the desaturation event, however, is an ensuing \"overshoot\" and a transient period of relative hyperoxia above baseline ([@R12]). Clinical studies demonstrate that the majority of desaturation events treated with O~2~ result in higher O~2~ saturations compared to an infant\'s target range ([@R13]). Furthermore, there is often a period of several minutes where the delivered FIO~2~ is greater than baseline after the event resolves ([@R13]). These episodes of reoxygenation following an O~2~ desaturation event may be expected to exacerbate oxidative stress, but their effects on long term airway reactivity and wheezing disorders in premature infants are unknown.

Murine lung development at birth is in the saccular stage of alveolarization, similar to the premature human lung at 24-36 weeks gestation, making it a suitable candidate to model lung disease secondary to prematurity ([@R14]). This model has been validated in numerous animal models of BPD which frequently utilize sustained and continuous hyperoxia exposure to replicate the BPD phenotype ([@R15]). Animal models of chronic intermittent hypoxia have similarly been shown to worsen allergen-induced pulmonary inflammation and have deleterious effects on respiratory function ([@R16]). Chronic intermittent hypoxia is also associated with increased levels of NOX2 and reactive oxygen species (ROS) production ([@R17]). Furthermore, re-oxygenation with 50% FIO~2~ after acute hypoxia exposure resulted in a transient surge in superoxide anion formation ([@R18]). This would suggest that hypoxia, followed by re-oxygenation with hyperoxia may be a potent stimulus of ROS formation. Despite these findings, the specific patterns that predispose to later lung airway injury have not been well characterized. In this study, we developed a novel neonatal mouse model of chronic intermittent hypoxia exposure with hyperoxia overshoots to determine if they have a deleterious effect on subsequent airway hyper-reactivity and respiratory mechanics. We hypothesized that neonatal chronic intermittent hypoxia with superimposed intermittent hyperoxia causes an increase in airway responsiveness to methacholine challenge and is associated with changes in pro-oxidant enzymes.

RESULTS {#S6}
=======

Body Weight {#S7}
-----------

Average body weight was not significantly different between groups at P21 (N~X~: 8.9 ± 1.9g; CIH~E~: 8.1 ± 1.5g; CIH~O/E~: 8.4 ± 1.4g; and CIH~O~: 8.4 ± 1.9g; p=0.25).

Airway Responsiveness to Methacholine {#S8}
-------------------------------------

Within each group there was a significant dose-dependent increase in Rrs ([Figure 1a](#F1){ref-type="fig"}) and decrease in Crs ([Figure 1b](#F1){ref-type="fig"}) when challenged with methacholine. CIH~O/E~ and CIH~E~ groups demonstrated increased respiratory system resistance at higher methacholine doses (50-200 mg/mL) compared to animals raised in room air ([Figure 1a](#F1){ref-type="fig"}; p\<0.05). Airway reactivity, defined as the difference between maximal contraction at the highest dose of methacholine and baseline, did not increase in the CIH~O/E~ or CIH~E~ groups compared to N~X~ (p=0.15 and p=0.53 respectively). There was a consistent decrease in specific respiratory system compliance at all methacholine doses among the CIH~O/E~ and CIH~E~ groups compared to normoxia ([Figure 1b](#F1){ref-type="fig"}; p\<0.05). There were no differences in resistance or compliance responses observed in the CIH~O~ group compared to room air controls.

Baseline Respiratory Mechanics {#S9}
------------------------------

Baseline respiratory mechanics at P21 are shown in [Figure 2](#F2){ref-type="fig"}. No differences in Rrs were observed between groups ([Figure 2a](#F2){ref-type="fig"}) at baseline. Baseline sCrs was decreased and specific elastance (sH) increased in the CIH~O/E~ and CIH~E~ groups compared to room air controls (both p\<0.05; [Figure 2b and 2c](#F2){ref-type="fig"}). Similarly, specific tissue damping (sG) was increased in the CIH~O/E~ and CIH~E~ groups compared to room air controls (p\<0.05, [Figure 2d](#F2){ref-type="fig"}). The CIH~O/E~ group demonstrated a further increase in sG compared to the CIH~E~ subjects (p\<0.05, [Figure 2d](#F2){ref-type="fig"}). No significant baseline differences were seen in the CIH~O~ group compared to room air controls.

Real Time PCR {#S10}
-------------

In an effort to describe changes in prooxidants and antioxidants gene expression, real time PCR of lung homogenates was performed at P21 for various superoxide dismutases (SOD) and NADPH oxidases (NOX). The antioxidants SOD1, SOD2, and SOD3 and their relative fold-change from normoxia are shown ([Figure 3a](#F3){ref-type="fig"}). SOD1 and SOD3 were decreased in CIH~E~ animals compared to normoxia but unchanged among other groups. SOD2 was unchanged in all groups. The prooxidants NOX1, NOX2, as well as the transcription factor and antioxidant-regulating protein nuclear factor (erythroid-derived-2)-like2 (NRF2) were also evaluated ([Figure 3b](#F3){ref-type="fig"}). NOX1 was decreased in CIH~O~ compared to all other groups (p\<0.05). NOX2 was increased in CIH~O/E~ compared to all other groups (p\<0.05). NRF2 was decreased in CIH~O~ compared to CIH~O/E~.

Immunohistochemistry {#S11}
--------------------

To assess long-term effects of intermittent hypoxia and hyperoxia on alveolarization, we quantified radial alveolar counts (RAC) at P21 ([Figure 4](#F4){ref-type="fig"}). There were no significant differences in RAC between groups. There were also no differences in the amount of airway α-smooth muscle actin between treatment groups. Representative images and staining for α-smooth muscle actin are demonstrated in [Figure 5](#F5){ref-type="fig"}.

DISCUSSION {#S12}
==========

Prematurity is an independent risk factor for wheezing disorders among extremely and moderately preterm infants, attributed to structural lung immaturity, pulmonary surfactant deficiency, mechanical ventilation, and supplemental oxygen therapy ([@R19]). Furthermore, emerging evidence implicates even late preterm birth as a risk factor for impaired respiratory function attributable to differences in lung and chest wall compliance, airway tethering, and susceptibility to postnatal infections such as RSV ([@R20]). Among infants born preterm, small airway obstruction is evident with increased odds of childhood asthma and decreased forced expiratory volume (FEV~1~) among infants with and without BPD ([@R21],[@R22]). We pursued this study to determine if a mouse model of FIO~2~ fluctuations, and delivery of specific patterns of oxygen exposure, contribute to pulmonary morbidity by altering respiratory mechanics, lung architecture, airway smooth muscle, and changes in prooxidant and antioxidant signaling.

In these experiments, CIH~O/E~ increased airway resistance at higher methacholine doses, decreased compliance, increased elastance, and increased tissue damping in pulmonary function testing. Specific tissue damping (sG) characterizes energy dissipation in the lung and markedly increases when there are regional differences in lung ventilation ([@R23]). sG increased in CIH~O/E~ compared to CIH~E~ suggesting that hypoxia immediately followed by hyperoxia may have some synergistic negative effect on regional respiratory mechanics. Chronic intermittent hypoxia alone had no effect on pulmonary mechanics, lung architecture, or airway smooth muscle compared to N~X~ controls. We did not, however, observe any corresponding changes in radial alveolar counts, airway smooth muscle hypertrophy, and mRNA expression for antioxidant and prooxidant markers except NOX2 in our CIH~O/E~ treated mice. Animals exposed to CIH~E~ had changes almost identical to CIH~O/E~ with the exception that increased NOX2 mRNA was unique to the latter group. Although these data demonstrate complex changes in lung mRNA for various pro- and anti-oxidant enzymes, the most significant observation was the increased Rrs to methacholine in the CIH~O/E~ group was associated with an increase in mRNA for the pro-oxidant NOX2. Taken together, these findings suggest that minute-to-minute swings in FIO~2~ above and below 21% are relevant stimuli that contribute to altered pulmonary mechanics.

The nature of intermittent hypoxia and hyperoxia events in preterm infants changes with advancing postnatal age in both timing and frequency ([@R11]). Altering these patterns places different stressors on the developing airway, potentially altering signaling pathways and disrupting lung development. Our study utilized a pattern of exposure with FIO~2~ cycling between 10% and 50% in 10-minute intervals for 24hrs/day, differentiating this exposure paradigm from other conventional exposures of sustained continuous hyperoxia that have shown airway hyperreactivity, airway smooth muscle hypertrophy, and decreased bronchiolar-alveolar attachments ([@R24],[@R25]). In addition, our model addresses the intermittent hypoxia/hyperoxia fluctuations that take place in neonatal intensive care units on the order of minutes, not hours and mimics the frequency of events (up to 100/day) documented in the NICU ([@R11]). We demonstrated deleterious effects on pulmonary function with an average FIO~2~, calculated as the \"area under the curve\", of 30-32% in the CIH~O/E~ group. These physiologic changes persist at P21, roughly correlating to the pulmonary stage in development for a toddler-aged child ([@R14]), a common point where respiratory morbidity is described from chronic lung disease due to prematurity ([@R26]). Chronic intermittent hypoxia in the absence of hyperoxia overshoot did not induce long term airway morbidity. Therefore, our characterization of at-risk patterns of intermittent hypoxia and hyperoxia may aid in guiding bedside management of labile preterm infants experiencing frequent apnea and desaturations to reduce hyperoxia overshoots and potentially their long-term risk of pulmonary morbidity.

Studies with different timing, frequency, and severity of intermittent hypoxia/hyperoxia compared to CIH~O/E~ have observed dysfunction in angiogenic and alveolar signaling pathways. One mouse model with single, 2-hour cycles of hypoxia (12%) with hyperoxia (50%) for 7 days demonstrated downregulation of HIF-1α and VEGF, implicating altered angiogenesis and thus impaired alveolar development ([@R27]). A different paradigm of clustered intermittent hypoxia superimposed on sustained hyperoxia showed alterations in VEGF levels normalized by treatment with a superoxide dismutase mimetic, further implicating altered angiogenesis as a factor ([@R28]). It is clear that the amount, severity, timing, and pattern of these swings in oxygen exposure can have wide-ranging effects on several signaling pathways, making it an intriguing area for further study.

The mechanism by which intermittent hypoxia with hyperoxia contributes to deleterious lung function is still not well established. In contrast to previously published studies of intermittent hypoxia/hyperoxia at one week ([@R27]), alveolarization was not affected in any study group following room air recovery at P21. We proceeded to investigate changes in signaling pathways that might be implicated in functional changes in the lung. RT-PCR for both prooxidant and antioxidant markers demonstrated increased NOX2 in only the CIH~O/E~ group. NOX2 was the first protein in the family of NADPH oxidase proteins to be discovered and is present in many different cell types with highest expression in phagocytes ([@R29]). It produces superoxide, peroxynitrites, and other reactive oxygen species that are involved in host defense ([@R29]). Counts for alveolar macrophages ([@R30]) were not different between the CIH~O/E~ group and other groups (data not shown), suggesting the increase in NOX2 signaling is independent of the presence of inflammatory cells in the lung parenchyma. Further investigation is needed to investigate what role NOX2 is playing in this CIH~O/E~ paradigm and how it interacts with other physiologic stressors.

We have evaluated changes in both airway smooth muscle and lung parenchyma as both may contribute to increased airway reactivity. The absence of airway smooth muscle hypertrophy is consistent with no change in baseline Rrs and decreased compliance might be associated with decreased alveolarization although we did not detect changes in alveolar counts. Therefore, we speculate that the pattern of intermittent exposures is having an effect on lung function that may not be anatomically detectable at this time point. These findings contrast with other models of sustained hyperoxia demonstrating increased airway smooth muscle hypertrophy with 40% sustained FIO~2~ ([@R25]) underscoring the importance of studying fluctuations in inspired FIO~2~ as they relate to pulmonary function and airway reactivity.

In conclusion, the CIH~O/E~ paradigm of intermittent hypoxia with hyperoxia overshoots for 7 days has long-term deleterious effects on lung function at P21 after 2-week recovery in room air but no detectable effect on alveolarization or smooth muscle hypertrophy. This study, like others underscores the negative effect specific patterns of oxygen exposure can have on lung development and lung function. Piecing together the cumulative dose of oxygen with the pattern in which it is delivered may provide valuable insight into how we respond to these events at the bedside. Increased oxidative stress may be implicated, however, the intermittent nature of these exposures may elicit a different signature of oxygen toxicity compared to models of sustained hyperoxia. These experiments are a first step in characterizing one pattern of oxygen exposure that may be implicated in long-term respiratory morbidity in preterm infants.

METHODS {#S13}
=======

Animals and Exposures {#S14}
---------------------

Animal protocols were approved by the Institutional Animal Care and Use Committee at Case Western Reserve University. Timed pregnant C57BL/6 female mice (Jackson Laboratories, Bar Harbor, ME) were purchased and housed in an AAALAC accredited animal facility on 12- hour light-dark cycles with ad libitum standard food and water. Newborn mouse pups were pooled within 24 hours of birth and randomly redistributed into treatment groups. Each litter and the nursing dam was placed in a modified exposure chamber designed to permit rapid exchange of chamber gas composition according to the desired exposure paradigm. Gas flow into the chamber was modified using individual mass flow controllers (SDPROC, AALBORG Instruments & Controls, Inc., Orangeburg, New York), programmed with a command module and corresponding software. Each flow control was independently controlled to allow airflow of gases (air, nitrogen or oxygen) into the chamber to achieve the desired exposure paradigm. An O~2~ sensor was placed into the exposure chamber to continuously record achieved FIO~2~. A schematic of the experimental setup is shown in [Figure 6](#F6){ref-type="fig"}.

To simulate a desaturation with hyperoxia overshoot, one group of mice were exposed to repeating cycles of hypoxia (10% O~2~) immediately followed by rebound hyperoxia (50% O~2~). After the chamber reached 50% O~2~, room air (21% oxygen) was purged into the chamber allowing gradual return of chamber conditions to normoxia within 5 minutes ([Figure 6](#F6){ref-type="fig"}). This exposure was repeated at 10 minute intervals, 24hrs/day for 7 days starting at P0. This treatment group that received intermittent hypoxia and hyperoxia overshoot was designated CIH~O/E~. Additional groups of mice were used to assess the separate effects of chronic intermittent hypoxia (CIH~O~) or hyperoxia (CIH~E~). All groups were compared to control mice raised in room air for the entire exposure period. After 7 days of respective exposures, animals were returned to room air until respiratory function testing or tissue harvesting at P21.

Respiratory Mechanics {#S15}
---------------------

Respiratory function testing was performed at P21 for all groups (N=7-11 per group). Each animal was placed under general anesthesia (intraperitoneal ketamine 200mg/kg and xylazine 20mg/kg), tracheostomized via a 19-gauge blunt tip cannula, and placed on a commercially available small mammal ventilator (flexiVent, SCIREQ, Montreal, Canada). Each subject was then paralyzed (pancuronium bromide 10mg/kg) and mechanically ventilated with a tidal volume of 10ml/kg, 150 breaths/min, PEEP of 3cm H~2~O, and FIO~2~ of 50%. Animals equilibrated for five minutes on the ventilator, then two recruitment breaths of sustained inspiration up to a pressure of 30 cm H~2~O for 3 seconds were administered followed by measurement of baseline respiratory mechanics. A dose-response curve was generated by administering aerosolized (Aeroneb, SCIREQ) methacholine (Sigma Aldrich, St. Louis, MO) for 10 seconds at increasing doses (saline control, 12.5, 25, 50, 100, and 200mg/ml). Respiratory System Resistance (R~rs~) and Specific Respiratory System Compliance (sC~rs~) were calculated by the ventilator software (Flexiware Version 7.2 Service Pack 2 Build 728, SCIREQ, Montreal, Canada) using a 1.2 second, 2.5 Hz single-frequency forced oscillation maneuver ([@R31]). Specific Tissue Damping (sG) and Specific Elastance (sH) were calculated using a low-frequency forced oscillation technique imparting a quasi-sinusoidal maneuver containing a wide range of frequencies. Input impedance was then analyzed through the previously described Constant Phase Model([@R23]). Specific values for respiratory function parameters were determined by dividing the result by the inspiratory capacity of each animal at baseline.

Real Time PCR {#S16}
-------------

RNA was extracted from snap-frozen lungs of 21 day old mice (N=7-16 per group) using TriZol reagent (Life Technologies, Carlsbad, CA) and quantified by Nanodrop spectroscopy (Thermo Scientific, Waltham, MA). cDNA was generated from one microgram of RNA by reverse transcription using qScript cDNA synthesis kit (Quanta Biosciences, Gaithersburg, MD). Real-time quantitative PCR was performed on a StepOne PCR system (Applied Biosystems, Foster City, CA) using TaqMan probes (Life Technologies) for SOD1 (Mm01344233_g1), SOD2 (Mm01313000_m1), SOD3 (Mm01213380_s1), NOX1 (Mm00549170_m1), NOX2 (Mm01287743_m1), NRF2 (Mm00801891_m1) compared to GAPDH (\#4351309) with PerfeCTa qPCR FastMix, UNG, ROX (Quanta Biosciences, Beverly, MA). Fold-changes are reported utilizing 2\^-ddCT method and StepOne software v2.3 (Applied Biosystems).

Histology and Morphometric Analysis {#S17}
-----------------------------------

Animals (N=7-12 per group) were sacrificed using anesthetic overdose of ketamine/xylazine and lungs inflated with 4% paraformaldehyde at 25cm H~2~O pressure on postnatal day 21. Paraffin-embedded sections (5μm) were prepared and stained with hematoxylin-eosin. Randomly chosen hematoxylin-eosin areas were photographed with 10X magnification from 7-10 animals per group. All of the morphometric and image analysis measurements were performed by a blinded investigator. The extent of alveolarization was determined by the radial alveolar count (RAC) method as previously described ([@R32]). Six terminal bronchioles per pup were analyzed for radial alveolar count, determined by counting the number of alveolar septa transected by a perpendicular line drawn from the terminal bronchiole to the nearest connective tissue septum.

Immunofluorescence and Analysis for Airway Smooth Muscle Actin {#S18}
--------------------------------------------------------------

The trachea was cannulated and the lungs (N=8-10 per group) were inflation fixed (25 cmH~2~O) for 10 min with 10% neutral-buffered formalin. The left lung was removed, prepared for immuno-staining, and post fixed for 2 days at 4°C. The fixed lung was embedded in paraffin and sectioned (5μm). Sections were de-paraffinized using xylene and then dehydrated in graded alcohol. Deparaffinized tissue sections were pretreated for antigen retrieval, and incubated with a monoclonal antibody against mouse anti α-smooth muscle actin (α-ASM; 1:400 dilution, Sigma-Aldrich, St. Louis, MO) overnight at 4°C. Immune complexes were captured with FITC-conjugated donkey anti-mouse (1:500, Alexa Fluor-488, Invitrogen, Carlsbad, CA). Appropriate negative controls were run by omitting the primary antibody to confirm nonspecific staining, which was not observed. Immunostained sections were cover slipped with Vectashield mounting medium (Vector Laboratories, Burlingame, CA) and visualized with a fluorescence microscope. Images of the immunostained sections were captured with a Rolera XR CCD camera (Q-Imaging, Surrey, Canada) mounted on a microscope (Leica Microsystems, Wetzlar, Germany). Five random images at 20X magnifications were captured from five to seven animals per group by using digital image analysis software with settings for color and size identification (*Image j*, NIH, Bathesda, MD). The area of airways and the green fluorescent areas α-ASM from five airways per animal were measured using the circular method ([@R33]). The amount of α-ASM area was expressed as α-ASM staining divided by airway basement membrane perimeter (*P*~bm~).

Statistical Analysis {#S19}
--------------------

Data are expressed as mean ± standard error of the mean with exception of animal weights (mean ± SD). Respiratory mechanics at baseline and for each dose of methacholine exposure were performed using two-way ANOVA with repeated measures. Analysis of mRNA expression was performed with a 1-way ANOVA on the ΔCT. Radial alveolar counts and alpha smooth muscle were performed with 1-way ANOVA. Statistical significance was defined as p \< 0.05.

**Statement of Financial Support:** This work was supported by the American Academy of Pediatrics Marshal Klaus Perinatal Research Award (Elk Grove Village, IL) and an endowment from the Rainbow Babies and Children's Foundation Fellowship Research Award Program (Cleveland, OH). Additional support provided by the National Institutes of Health T32-HD060537 training grant (Bethesda, MD). Dr. Dylag is the recipient of the National Institutes of Health Loan Repayment Program award (Bethesda, MD).

**Disclosure Statement:** The authors report no conflict of interest.

**Category of Study:** Translational

![Respiratory System Resistance (Rrs) and Specific Respiratory System Compliance (sCrs) and responses to methacholine in mechanically ventilated mice at P21. (a) Significant increases in Rrs were observed at higher methacholine doses of both groups who received intermittent hyperoxia (CIH~O/E~, dashed line closed triangles and CIH~E~, dashed line open circles) compared to room air controls (N~X~, solid line closed circles). (b) Significant decrease in specific respiratory system compliance in CIH~O/E~ and CIH~E~ at all doses compared to room air controls. No changes were observed between CIH~O~ (solid line, open triangles) and Nx. N=7-11 per group. p \< 0.05: \* CIH~O/E~ vs N~X~, †CIH~E~ vs N~X~.](nihms-829105-f0001){#F1}

![Baseline Respiratory Mechanics at P21. (a) Respiratory System Resistance (Rrs) (b) Specific Respiratory System Compliance (sCrs) (c) Specific Elastance (sH). (d) Specific Tissue Damping (sG). N=7-11 per group. \* p \< 0.05](nihms-829105-f0002){#F2}

![Real-Time PCR of Lung Homogenates. RNA was extracted from whole lung homogenates of pups at P21. Fold-change of mRNA expression are demonstrated for groups CIH~E~ (second from left, black), CIH~O/E~ (third from left, gray) and CIH~O~ (right, hatched) normalized to Nx (left, white) mice. (a) Antioxidant markers SOD1, SOD2, and SOD3 (b) Prooxidant makers of NOX1, NOX2, and NRF2. N=7-16 per group. \*p \< 0.05.](nihms-829105-f0003){#F3}

![Radial Alveolar Counts for all groups at P21. (a) Graph of counts for all four groups (b) Representative images of each group. No differences were seen between groups. Magnification bar 50μm. N=7-12 per group.](nihms-829105-f0004){#F4}

![Alpha-Actin in airway smooth muscle at P21. (a) Graph of smooth muscle alpha-actin expression in all four groups (b) Representative images from each group. No significant differences were observed between groups. Magnification bar 50μm. N=8-10 per group.](nihms-829105-f0005){#F5}

![(a) Schematic of Experimental Setup. Independent gas inputs of Nitrogen (N~2~), Oxygen (O~2~) and Room Air (RA) pass through computer-programmed flow controllers each connected to a command module and operated by a computer. Gas inputs feed into the animal cage attached to oximeter. Oxygen exposure paradigm for each experimental group with continuously monitored levels: (b) Chronic Intermittent Hypoxia/Hyperoxia (c) Chronic Intermittent Hyperoxia (d) Chronic Intermittent Hypoxia. Room air controls not illustrated. Two ten-minute repeating cycles illustrated. Cycles repeat 24hrs/day.](nihms-829105-f0006){#F6}
