PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Pediatr Pulmonol. Author manuscript; available in PMC 2012 April 22.
Published in final edited form as:
PMCID: PMC3332105
NIHMSID: NIHMS360607

High-Flow Nasal Cannula: Impact on Oxygenation and Ventilation in an Acute Lung Model

Meg Frizzola, DO,*,1 Thomas L. Miller, PhD,2,3 M.E. Rodriguez, MD,1 Yan Zhu, MD,1 Jorge Rojas, MD,4 Anne Hesek,1 A. Stump,1 Thomas H. Shaffer, MSE, PhD,1,2,5 and Kevin Dysart, MD1,2

SUMMARY

Introduction

High-flow nasal cannula therapy (HFNC) has been shown to be more effective than continuous positive airway pressure (CPAP) in reducing intubations and ventilator days. HFNC likely provides mechanisms to support respiratory efficiency beyond application of distending pressure. We reason that HFNC washout of nasopharyngeal dead space impacts CO2 removal along with oxygenation. The aim of this study was to demonstrate the flow dependence of CO2 reduction and improved oxygenation during HFNC and the dependence on leak around the nasal prongs.

Materials and Methods

Neonatal piglets (n=13; 2-6kg) were injured with IV oleic acid and supported with HFNC at 2 through 8 L/minute. High and low leak around the nasal prongs was accomplished by using single and double prong cannulae, respectively. Measurement of hemodynamic, respiratory and blood gas parameters were made at each setting following 10 minutes for physiologic equilibration. Tracheal pressures were recorded by transmural catheters.

Results

With HFNC, CO2 trended downward in a flow dependent manner independent of leak. Oxygenation and tracheal pressures increased in a flow dependent manner with the greatest effect during double prong. At 8L/minute, tracheal pressures did not exceed 6±1 cmH2O.

Conclusions

HFNC improves gas exchange in a flow dependent manner; double prong had greater impact on O2; single prong had greater impact on CO2 elimination.

Keywords: high-flow nasal cannula, animal model, gas exchange, dead space, lung mechanics

INTRODUCTION

Continuous positive airway pressure (CPAP) is often a preferred intervention over mechanical ventilation with neonatal respiratory distress syndrome in that oxygenation is improved and lung volume is recruited while supporting spontaneous ventilatory efforts.1,2 Some essential clinical criteria to remain on CPAP are effective spontaneous respiratory effort and CO2 elimination. Hypercapnia, or apnea that may be secondary to hypercapnia, are cited as some of the more common reasons for progressing to more invasive forms of ventilatory support. Therefore, it is thought that if CO2 retention during CPAP can be reduced or eliminated, many infants can be spared potential lung injury and subsequent chronic lung disease induced by mechanical ventilation.

High-flow nasal cannula (HFNC) (>2 liters per minute, (lpm)) has been proposed by some physicians as an alternative to CPAP in maintaining infants without the need for intubation and mechanical ventilation.3 The use of HFNC in place of conventional methods has been controversial and as yet not well documented in the literature as being more effective than CPAP. Given the paucity of translational research investigating the mechanism of action of HFNC, we sought to investigate this mechanism further. We hypothesized that the mechanism of gas exchange with HFNC is associated with the washout of nasopharyngeal anatomical dead space in addition to the potential impact of positive airway pressure generation.

In a recent study, we used tracheal gas insufflation (TGI) with CPAP in lung injured piglets using a specially designed endotracheal tube to allow for insufflation of fresh gas flow (No. 6501.30; Vygon, Ecouen, France).4 The results showed asignificant reduction in CO2 retention associated with prosthetic dead space washout in spontaneously breathing piglets with oleic acid induced lung injury. Such a therapy, in non-intubated patients, may increase the number of infants that can be sustained without intubation or mechanical ventilation and these infants may be weaned sooner.

Given our TGI results, we reasoned that the beneficial clinical responses to HFNC may be related to the high-flow rates effectively flushing the nasopharyngeal cavity of expiratory gas. Thus, we proposed that the HFNC may provide a practical means to wash out nasopharyngeal dead space during spontaneous breathing without the need for intubation. The present study was designed to demonstrate this washout effect and the flow dependence of CO2 elimination and oxygenation during HFNC. In addition, the study was conducted to demonstrate the independent effects of airway pressure under different leak conditions. We also hypothesized that HFNC therapy would affect CO2 retention in a flow dependent manner without markedly increasing intra-tracheal pressure and improve respiratory parameters compared to a control group of animals treated with mask CPAP (minimal leak).

MATERIALS AND METHODS

Animal Preparation

Thirteen spontaneously breathing newborn piglets (15-20 days old, weight 2-6 kg) were instrumented, placed on mask CPAP of 5 cmH2O and injured by intravenous administration of oleic acid. All thirteen piglets were treated with CPAP (minimal leak), single prong (SP) HFNC with a high degree of leak around the nasal prongs (HI LEAK), and double prong (DP) HFNC with a low degree of leak around the nasal prongs (LOW LEAK). The three treatment conditions were administered in a randomized, repeated measure, crossover design described below. Following experimentation, animals were sacrificed with pentobarbital (50 mg/kg) and saturated potassium chloride (2 mEq). At baseline, injury and at each treatment condition (following 10 minutes for physiologic equilibration) measurements were made of arterial blood chemistry, pulmonary function and intra-tracheal pressure. All procedures for the animal protocol in this study were approved by the Life Science Center Animal Care and Use Committee at Nemours Children’s Clinic - Wilmington of The Nemours Foundation.

Instrumentation and Injury Protocol

The piglets were anesthetized with an anesthetic solution (ketamine: 23 mg/kg; azepromazine: 0.58 mg/kg; and xylazine: 0.8 mg/kg [KAX]) given as two 1 ml/kg intramuscular injections separated by 10 minutes. The skin and soft tissues were locally infiltrated with 0.5% lidocaine HCl (4 mg/kg) and instrumented with 5.0 F catheters placed in the jugular vein and carotid artery. Additionally, a 5.0 F fluid filled catheter was placed directly into the mid trachea approximately 2 cm below the cricoid ring at a perpendicular orientation for measurements of dynamic tracheal pressure. The tracheal catheter insertion site was sealed with surgical glue to eliminate leak around the catheter. Subsequent anesthesia was maintained with an intravenous infusion of KAX at 0.4 ml/kg/hr, along with diazepam every two hours and additional KAX every one half an hour. Maintenance fluid was provided by a continuous infusion of 5% dextrose solution at a rate of 6 ml/kg/hr. Arterial blood pressure was monitored by attaching the arterial catheter to a standard pressure transducer via bedside patient monitor (Model M1175A, Hewlett Packard). ECG electrodes were placed for monitoring. Throughout the protocol the animal’s rectal temperature was monitored and maintained at 37-38°C on a radiant warmer bed (Resuscitaire®; Hill-Rom Air-Shields, Hatboro, PA).

Following instrumentation, baseline measurements were made and then lung injury was initiated. Oleic acid (0.08 ml/kg; Sigma) was emulsified in 0.5 ml/kg blood and infused into the central venous line in four equal doses with ample time to allow stabilization between doses. The injury model has been previously described4 and the oleic acid dose was refined in a pilot study to achieve a 50% reduction in PaO2 and respiratory compliance while preserving spontaneous respiration.

Treatment Protocol

Following the injury procedure and a 30 minute stabilization period, all thirteen animals were exposed to HI LEAK, LOW LEAK and CPAP (minimal leak) in randomly assigned order. For each treatment scenario, flow rate (HFNC) or pressure (CPAP) was increased incrementally by 2 (lpm or cmH2O) from least to greatest magnitude and then incrementally back to the least; the two sets of data for each increment were then averaged to filter out the residual effects of the previous level of treatment. To be consistent with current neonatal clinical practice, HFNC flow rates ranged from 2 to 8 lpm and CPAP pressures used ranged from 2 to 6 cmH2O flow was held constant at 8 l/min. Inspired oxygen fraction was always held at 1.0 for consistency and so that PaO2 reflects the a/A ratio.

HFNC therapy was administered using a Vapotherm 2000i (Vapotherm, Inc, Stevensville, MD) with a pediatric nasal cannula. Given the small diameter of the piglet nares relative to the human, the pediatric cannula fully occluded both nares of the piglet, thus representing the LOW LEAK condition. To achieve the HIGH LEAK condition where only half the area of the nares was obstructed (as per manufacturer’s recommendation), one of the prongs of a pediatric cannula was removed and the hole sealed, such that one nare of the piglet was occluded with the cannula prong and the other nare was open to exhaust.

CPAP was administered by a VIP Bird infant ventilator (Bird Products Corp., Palm Springs, CA) set to the CPAP mode. The ventilator circuit was attached to the piglet using a semi-cone shaped face mask lined with silicone that conformed to the piglet’s head shape to minimize the addition of prosthetic dead space. Furthermore, the base of the face mask (around the piglets head) allowed for minimal gas leak to allow washout of any remaining prosthetic dead space from the mask. Under all conditions, the ventilator was able to generate the prescribed levels of CPAP as indicated by ventilator circuit and pig tracheal pressure measures.

Pulmonary Function Assessment

Arterial blood gases and chemistry were measured by a standard patient blood analyzer (Stat Profile®, Nova Biomedical, Waltham, MA). Blood gas parameters measured included pH, PaCO2, PaO2, hemoglobin, hematocrit, HCO3, and base excess values. The blood gas measurements were obtained 10 minutes after each change in respiratory support to ensure adequate physiologic equilibration given that the respiratory system compensates for changes in acid-base status within this time period.5 Respiratory parameters were assessed by inductive plethysmography (SomnoStar PT, SensorMedics, Yorba Linda, CA), including relative tidal volume, respiratory rate and thoracoabdominal synchrony. In all animals, relative tidal volumes, respiratory rates and minute ventilations were determined prior to lung injury during manual bagging (via mask CPAP) with low (2-4 ml/kg) and high (6-8 ml/kg) tidal volume strategies in order to establish relative calibration standards for inductive plethysmography.

Data Analysis

Data was analyzed using linear regression or analysis of variance (ANOVA) to assess differences associated with type and/or magnitude of therapy. Post-hoc analyses were done, where appropriate, using Bonferroni comparisons. Significance was accepted with p values < 0.05.

RESULTS

Thirteen piglets were utilized for the study and demonstrated a stable injury; age 13 ± 8 days and weight 6.0 ± 0.2 kg. The injury was associated with a substantial increase in PaCO2 (37.3 mmHg ± 6.6 vs 49.1 mmHg ± 6.1; p < 0.01) despite the transition from room air to FiO2 of 1.0 and an increase in A-a gradient (23.3 +/− 12.1 vs 254.3 +/− 88.9; p < 0.01). These findings are demonstrated in Figure 1.

Figure 1
Baseline and Injury

Tracheal Pressure

As shown in Figure 2, it was possible to record direct temporal changes in tracheal pressure at each specific therapeutic intervention. A data acquisition recorder (DASH 8Xe) was used for the tracheal measurements. As such, these traces provided evaluation of pressure response to increments in CPAP or flow during HFNC perturbations. Based on summarized data shown in Figure 3A and 3B, there was a direct linear relationship between each prescribed CPAP setting and tracheal pressure (slope=0.6;r2=0.99;p=0.02). Similarly, a direct linear relationship existed between flow rate and tracheal pressure under the LOW LEAK condition (slope=0.21;r2=0.82; p =0.04). Under the HIGH LEAK condition, the linear relationship between flow and resultant tracheal pressure showed a significant rise in tracheal pressure between 2 and 8 lpm of flow (r2=0.81;p= 0.04). Overall, mean tracheal pressures were higher during the low leak condition compared to high leak flow conditions.

Figure 2
Real-time Tracheal Pressures
Figure 3
Tracheal Pressures

Pulmonary and Physiologic Parameters

No differences were seen in respiratory rate and relative tidal volume across therapies or orders of magnitude (i.e., incremental pressures (CPAP) or flow rates (HFNC). The relative minute ventilation values determined through integration of rates and relative tidal volumes (inductive plethysmography) also showed no differences as demonstrated in Table 1.

Table 1
Minute Ventilation Values (AU) as a Function Flow

Physiologic parameters were measured for all piglets under every condition. Measurements included rectal temperature, systolic blood pressure, diastolic blood pressure, mean arterial blood pressure, heart rate, respiratory rate and SaO2 values. The means and standard deviations for each parameter under all conditions are depicted in Table 2. There was no significance found for any of these values.

Table 2
Physiologic Parameters as a function of CPAP, SP, and DP across experimental conditions.

Ventilation Indices

Figure 4 shows the effects of incremental increases in CPAP and HFNC therapies on PaCO2. With initiation of CPAP (2 cmH2O), arterial carbon dioxide tension (PaCO2) decreased from the basal injury level (p<0.05); however, no further decease in PaCO2 was seen with the incremental increases in CPAP pressure. With HFNC, in the LOW LEAK condition a significant relationship was found between flow rate and partial pressure of carbon dioxide (slope=-0.73; r2 = 0.92;p<0.01) where PaCO2 decreased with incremental increases in flow rate in a sigmoidal fashion. In the HIGH LEAK condition, PaCO2 was restored to baseline levels immediately (2 lpm of flow) and therefore did not decrease further with incremental increases in flow rate.

Figure 4
Partial Pressure Carbon Dioxide

Oxygenation

Oxygenation versus CPAP/flow rate data is shown in Figure 5. With CPAP administration, arterial oxygen tension (PaO2) rose to over 400 mmHg but did not change with incremental increases in CPAP pressure. However, PaO2 showed a flow dependent increase for both LOW LEAK and HIGH LEAK conditions. In the LOW LEAK condition, PaO2 responded in a significant linear relationship with throughout the range of flow rates 2-8 lpm (slope=30.6; r2 = 0.85;p = 0.03), and there was a 400 mmHg plateau from 8-10 lpm. Under the HIGH LEAK condition, PaO2 rapidly increased from 2 to 6 lpm before reaching a 400 mmHg plateau above 6 lpm.

Figure 5
Partial Pressure Oxygen

DISCUSSION

This study demonstrated two important findings. First, with HFNC tracheal pressures were comparable to CPAP pressures at the same flow range. Second, washout of nasopharyngeal dead space is associated with improved gas exchange using HFNC and may allow patients to breathe more comfortably as represented clinically.6 As our results demonstrate, tracheal pressures in the neonatal piglets were no greater than with traditional CPAP ventilation and pressures. This is consistent with HFNC clinical studies that report lower esophageal pressures compared to CPAP at 6 cmH2O.7-11

With respect to gas exchange, the impact of increasing flow on ventilation and oxygenation occurred independent of tracheal pressure generation alone. With increasing CPAP pressure in this injured lung model, neither PaCO2 nor PaO2 demonstrated progressive change. However, with HFNC under both leak conditions PaCO2 and PaO2 improved in a somewhat flow dependent manner reflected by saturation curves, (PaCO2 decreased with increasing flow until saturation and PaO2 increased with increasing flow until saturation). These saturation relationships are consistent with nasopharyngeal dead space washout (like) as demonstrated in the literature from tracheal gas insufflation.

Tracheal gas insufflation is a means of flushing prosthetic dead space using a catheter inserted into an endotracheal tube, or a specially designed endotracheal tube for gas insufflation.6,12-15 With this technique, insufflation of fresh gas is fed through the catheter at low flows (typically 0.5-1.0 lpm). This method supplies fresh gas during inspiration, and also flushes the prosthetic tube dead space between the catheter tip and the oral/nasal openings during expiration and the end expiratory pause. By flushing prosthetic dead space during the expiratory phase of the respiratory cycle, the subsequent inspiratory breath contains more fresh gas and less residual end-expiratory gas. Therefore, alveolar gas fractions are moved toward fresh gas values; i.e., less carbon dioxide and more oxygen, thereby improving respiratory efficiency.16-18 However, the limitation of tracheal gas insufflation is that it requires intubation.4

With gas insufflation, at any given liter flow the volume of dead space to be flushed is determined by the time allotted by end-expiration and the expiratory pause.6,12-15 Therefore, greater gas flows will flush more dead space until a flow rate is reached that can flush all of the available anatomical/prosthetic dead space in the allotted time; beyond this liter flow, increased flow rates accomplish no further effect, generating a saturation cure for alveolar and thus arterial gas composition.

In the current study, tracheal pressures with HFNC rose in a linear fashion and did not exceed 6 cm H2O up to flows of 8 lpm, while both ventilation and oxygenation responses reached saturation near 6-8 lpm flow. Therefore, we determined that nasopharyngeal dead space washout was a more predominant factor in determining gas exchange response to incremental flow as compared to incremental pressure increases. As compared to CPAP and LOW LEAK, the partial pressure of carbon dioxide was lower for flow rates (< 6 lpm) in the high leak condition. We attribute this phenomenon to a better washout of the nasopharyngeal cavity with HIGH LEAK, making a high leak situation more desirable since effective gas exchange can be achieved at a lower tracheal pressure.

Although several recent articles have reported estimates of tracheal pressure, positive end distending pressure and continuous positive airway pressure generated during the use of high frequency nasal cannula therapy, these studies are confounded by several measurement limitations. The results to date have been collected by means of nasal pharyngeal pressure monitoring, oral cavity pressure monitoring or esophageal pressure monitoring. Although these methods are all minimally invasive; they provide only an estimation of actual airway pressure. The major complication with these forms of measurement are the following: inaccurate pressure assessment due to catheter placement or catheter design, catheter orientation with respect to the gas flow, frequency response of the measurement system, as well as variable gas leak from the oral cavity. In our study, we used precise, accurate, and responsive dynamic tracheal pressure monitoring/recording via a fluid filled catheter surgically placed mid trachea with a perpendicular orientation (with respect to gas flow). Using this technique, as described, we did not find significantly elevated tracheal pressures over the flow ranges tested. This is the first time this method has been used to directly measure tracheal pressures during high flow nasal cannula therapy.3, 7, 19-21

Taking into account the possibility of other physiologic adjustments in minute ventilation impacting alveolar and therefore blood gas compositions, we measured relative minute ventilation using inductive plethysmography. No changes in this relative minute ventilation were noted, indicating that changes in blood gas values were directly associated with the change in therapy paradigm. In addition, since therapeutic interventions were randomized during the study, there were no alterations in lung injury or alveolar dead space, which could have accounted for changes in gas exchange. Taken together, these data support the conclusion that HFNC may augment carbon dioxide elimination by means of nasopharyngeal dead space washout.

As in all animal model studies, a limitation to this study is the translation of data to clinical application. The anatomy of a piglet, especially the nasopharyngeal anatomy is very different from a human. In this regard, we believe that the cavernous human nasopharynx, as opposed to the elongated pig nasopharynx would flush more easily without as much pressure development. Nonetheless, we contend that the anatomy is similar enough to study physiologic mechanisms, and the piglet as model for treatment of lung injury has been well established.4,22

The current data suggest that HFNC therapy may be used as a first line option over traditional CPAP support. This study shows that HFNC provides improved levels of arterial CO2, whereas CPAP does not, independent of changes in minute ventilation. Therefore, HFNC can be viewed as a truly non-invasive means of ventilatory support. In addition, by way of nasopharyngeal dead space elimination, HFNC can serve as a means of oxygenation support independent of, or in addition to, supplemental oxygen administration. Clinical studies are now warranted to refine flow needs in the human scenario to optimize these effects demonstrated in this preclinical model.

ACKNOWLEDGEMENTS

The authors gratefully acknowledge The Nemours Foundation, National Institutes of Health COBRE grant and Vapotherm, Inc. (Stevensville, MD) in part for financial support of the study and provision of equipment.

Footnotes

This work was supported in part by a NIH COBRE P20 RR 020173 grant, a Nemours Foundation Grant and an unrestricted research/educational grant from VapoTherm, Inc. The authors declare no bias on behalf of the authors or Vapotherm, Inc., who supported this study.

Data from this manuscript was presented at the 2008 Annual Meeting of the Pediatric Academic Societies.

REFERENCES

1. Gregory GA, Kitterman JA, Phibbs RH, Tooley WH, Hamilton WK. Treatment of the idiopathic respiratory-distress syndrome with continuous positive airway pressure. N Engl J Med. 1971;284:1333–1340. [PubMed]
2. Saunders RA, Milner AD, Hopkin IE. The effects of continuous positive airway pressure on lung mechanics and lung volumes in the neonate. Biol Neonate. 1976;29:178–186. [PubMed]
3. Shoemaker MT, Pierce MR, Yoder BA, DiGeronimo RJ. High flow nasal cannula versus nasal CPAP for neonatal respiratory disease: a retrospective study. J Perinatol. 2007;27:85–91. [PubMed]
4. Miller TL, Blackson TJ, Shaffer TH, Touch SM. Tracheal gas insufflation-augmented continuous positive airway pressure in a spontaneously breathing model of neonatal respiratory distress. Pediatr Pulmonol. 2004;38:386–395. [PubMed]
5. Hall J. Regulation of Acid Base Balance. In: Guyton AC, Hall JE, editors. Guyton and Hall Textbook of Medical Physiology. 11th edition Elsevier Saunders; Philadelphia: 2006. pp. 383–401.
6. Danan C, Dassieu G, Janaud JC, Brochard L. Efficacy of dead-space washout in mechanically ventilated premature newborns. Am J Respir Crit Care Med. 1996;153:1571–1576. [PubMed]
7. Kubicka ZJ, Limauro J, Darnall RA. Heated, humidified high flow nasal cannula therapy: yet another way to deliver continuous positive airway pressure? Pediatrics. 2008;121:82–88. [PubMed]
8. Spence KL, Murphy D, Kilian C, McGonigle R, Kilani RA. High flow nasal cannula as a device to provide continuous positive airway pressure in infants. J Perinatol. 2007;27:772–775. [PubMed]
9. Wilkinson DJ, Andersen CC, Smith K, Holberton J. Pharyngeal pressure with high-flow nasal cannulae in premature infants. J Perinatol. 2008;28:42–47. [PubMed]
10. Saslow JG, Aghai ZH, Nakhla TA, Hart JJ, Lawrysh R, Stahl GE, Pyon KH. Work of breathing using high-flow nasal cannula in preterm infants. J Perinatol. 2006;26:476–480. [PubMed]
11. Kahn DJ, Courtney SE, Steele AM, Habib RH. Unpredictability of delivered bubble nasal continuous positive airway pressure: role of bias flow magnitude and nares-prong air leaks. Pediatric Res. 2007;62:343–347. [PubMed]
12. Dassieu G, Brochard L, Agudze E, Patkaï J, Janaud JC, Danan C. Continuous tracheal gas insufflation enables a volume reduction strategy in hyaline membrane disease: technical aspects and clinical results. Intensive Care Med. 1998;24:1076–1082. [PubMed]
13. Claure N, D’Ugard C, Bancalari E. Elimination of ventilator dead space during synchronized ventilation in premature infants. J Pediatr. 2003;143:315–320. [PubMed]
14. Burke WC, Nahum A, Ravenscraft SA, Nakos G, Adams AB, Marcy TW, Marini JJ. Modes of tracheal gas insufflation. Comparison of continuous and phase-specific gas injection in normal dogs. Am Rev Respir Dis. 1993;148:562–568. [PubMed]
15. Bernath MA, Henning R. Tracheal gas insufflation reduces requirements for mechanical ventilation in a rabbit model of respiratory distress syndrome. Anaesth Intensive Care. 1997;25:15–22. [PubMed]
16. Nahum A. Animal and lung model studies of tracheal gas insufflation. Respir Care. 2001;46:149–157. [PubMed]
17. Dassieu G, Brochard L, Benani M, Avenel S, Danan C. Continuous tracheal gas insufflation in preterm infants with hyaline membrane disease. A prospective randomized trial. Am J Respir Crit Care Med. 2000;162:826–831. [PubMed]
18. Nakos G, Lachana A, Prekates A, Pneumatikos J, Guillaume M, Pappas K, Tsagaris H. Respiratory effects of tracheal gas insufflation in spontaneously breathing COPD patients. Intensive Care Med. 1995;21:904–912. [PubMed]
19. Parke R, McGuinness S, Eccleston M. Nasal high-flow therapy delivers low level positive airway pressure. Br J of Anaesth. 2009;103:886–890. [PMC free article] [PubMed]
20. Lampland A, Plumm B, Meyers P, Worwa C, Mammel MC. Observational study of humidified high-flow nasal cannula compared with nasal continuous positive airway pressure. J Pediatr. 2009;154:177–182. [PubMed]
21. Locke RG, Wolfson MR, Shaffer TH, Rubenstein SD, Greenspan JS. Pediatrics. 1993;91:135–138. [PubMed]
22. Miller TL, Singhaus CJ, Sherman TI, Greenspan JS, Shaffer TH. Physiologic implications of helium as a carrier gas for inhaled nitric oxide in a neonatal model of Bethanecol-induced bronchoconstriction. Pediatr Crit Care Med. 2006;7:159–164. [PubMed]