Much remains unknown about the consequences of very low birth weight (VLBW) and bronchopulmonary dysplasia (BPD) on adult lungs. We hypothesized that VLBW adults would have impaired lung function compared with controls, and those with a history of BPD would have worse lung function than those without.
At age 26 to 30 years, 226 VLBW survivors of the New Zealand VLBW cohort and 100 term controls born in 1986 underwent lung function tests including spirometry, plethysmographic lung volumes, diffusing capacity of the lung for carbon monoxide, and single-breath nitrogen washout (SBN2).
An obstructive spirometry pattern was identified in 35% VLBW subjects versus 14% controls, with the majority showing mild obstruction. Compared with controls, VLBW survivors demonstrated significantly lower forced expiratory volume in 1 second (FEV1), FEV1/forced vital capacity (FVC) ratio (FEV1/FVC), forced expiratory flow at 25% to 75% of FVC and higher residual volume (RV), RV/total lung capacity (TLC) ratio (RV/TLC), decreased diffusing capacity of the lung for carbon monoxide, and increased phase III slope for SBN2. The differences persisted after adjustment for sex and smoking status. Within the VLBW group, subjects with BPD showed significant reduction in FEV1, FEV1/FVC, and forced expiratory flow at 25% to 75% of FVC, and increase in RV, RV/TLC, and phase III slope for SBN2, versus subjects without. The differences remained after adjustment for confounders.
Adult VLBW survivors showed a higher incidence of airflow obstruction, gas trapping, reduced gas exchange, and increased ventilatory inhomogeneity versus controls. The findings suggest pulmonary effects due to VLBW persist into adulthood, and BPD is a further insult on small airway function.
Ongoing airflow obstruction presents in survivors born low birth weight or very preterm in childhood through early adulthood. Studies assessing other aspects of lung function reported conflicting results and largely focused on subjects at young ages.
The findings suggest airflow obstruction, gas exchange inefficiency, and increased ventilatory inhomogeneity due to very low birth weight persist into adulthood, and bronchopulmonary dysplasia is a further insult on small airway function.
Survival rates of very low birth weight (VLBW) (birth weight <1500 g) or very preterm (VP) (gestation <32 weeks) infants have increased considerably from the late 1970s onward in highly developed countries. As these individuals reach adulthood, there is increasing interest in their health and welfare.1,2 In addition to respiratory care being a major focus in the neonatal period, respiratory problems are among the most common causes of morbidity in early and middle childhood among survivors,3 and lung function has been 1 of the most frequently studied health outcomes in adolescence.4 However, there are few population-based studies in which authors have investigated lung function in VLBW and/or VP survivors compared with term-born controls beyond late adolescence.4–7 Most studies have focused on airflow obstruction with few assessing other aspects of lung function.4,8,9 In addition, there are no New Zealand data in this crucial area.
Our aim with the current study was to compare the results of a range of lung function tests in the prospectively enrolled New Zealand 1986 VLBW national cohort with those in term-born controls.10 We hypothesized that VLBW adults would have impaired lung function compared with controls, and that within the VLBW group, those with a diagnosis of bronchopulmonary dysplasia (BPD) would have worse lung function than those without BPD.
Methods
The New Zealand 1986 VLBW National Cohort
The New Zealand 1986 VLBW cohort consists of survivors of the 413 VLBW infants who were live born and admitted to a neonatal unit in New Zealand that year, all of whom were included in a prospective audit of retinopathy of prematurity.10 Of these infants, 338 (82%) survived to be discharged from the hospital. Between March 2013 and October 2016, 250 cohort members (77% of the 323 known survivors), now aged 26 to 30 years, participated in a follow-up study, and 229 VLBW adults, together with 100 term controls, underwent assessment in Christchurch over 2 days, covering physical and mental health, physiologic assessments, and cognitive, neuropsychological, and social functioning (Supplemental Fig 2). Available perinatal data for the VLBW cohort included receipt of antenatal corticosteroids (ANS), days of assisted ventilation via an endotracheal tube, and a diagnosis of BPD (defined as receiving supplementary oxygen at 36 weeks’ postmenstrual age). Exogenous surfactant was not available in 1986. Participants in the control group, first recruited at 22 to 23 years, were born at term in New Zealand in 1986 and not admitted to a neonatal unit.10 Ethics approval was granted by the Upper South B Regional Ethics Committee, New Zealand. Written informed consent was given by all participants.
Lung Function Tests
The assessments included anthropometric measurements (height was measured by using a Harpenden stadiometer; weight was measured by using electronic scales) and a comprehensive questionnaire of health and welfare, which included smoking history and other items that were based on the European Community Respiratory Health Survey II.11 Lung function tests were undertaken at the Respiratory Physiology Laboratory, Christchurch Hospital, which is accredited by the Thoracic Society of Australia and New Zealand. All tests were performed on 1 day in the following order: spirometry, plethysmographic lung volumes, single-breath diffusing capacity of the lung for carbon monoxide (DLCO) (all with Autobox V62J; CareFusion, Yorba Linda, CA), and single-breath nitrogen washout (SBN2) (SensorMedics Vmax Spectra exercise test station; CareFusion). The tests were performed with the subject sitting and with a nose clip in place; respiratory physiologists were blinded as to the participants’ group. Testing protocols adhered to the American Thoracic Society and European Respiratory Society pulmonary function standards12–15 and interpretation.16 DLCO gas mixture contained 0.3% carbon monoxide and 0.3% methane with balance air. DLCO values were corrected for venous blood hemoglobin. Only technically acceptable and repeatable data were included in data analyses. Values for spirometry, lung volumes, and DLCO were presented as z scores and percentages of the predicted (% pred) on the basis of the reference values from the Global Lung Function Initiative (GLI) (2012),17 Quanjer et al (1993)18 and GLI (2017),19 respectively. The SBN2 results were reported as the mean of the measured phase III slopes.
Statistical Analysis
Between-group comparisons were conducted by using the independent samples t test for differences in means and the χ2 test for comparison of percentages. Mean between-group differences in lung function were adjusted for covariate factors by using multiple linear regression methods. With a nominal P value of .05, the study had 80% power to detect a mean difference in lung function between VLBW participants and controls in excess of 0.34 SD. For subgroup comparisons (eg, by BPD status within the VLBW cohort), the corresponding effect size estimates ranged between 0.40 and 0.45 SD. A Bonferroni correction was applied to account for multiple significance testing across the lung function measures, taking into account the average correlation between outcomes. The Bonferroni adjusted P value was P < .009. Statistical analyses were conducted by using SAS 9.4 (SAS Institute, Inc, Cary, NC), and graphing of outcomes were conducted by using GraphPad Prism version 7.
Results
Characteristics of Participants
Of the 323 known VLBW survivors, there were few differences between those assessed and not assessed in demographic and perinatal variables, except those assessed had a lower mean birth weight but were less likely to have moderate or severe disability at 7 to 8 years (Supplemental Table 6). A total of 228 VLBW survivors and 100 controls underwent lung function testing with valid data obtained from 226 VLBW survivors and all controls. The demographic and anthropomorphic data for VLBW adults and controls are shown in Table 1. Compared with controls, both men and women in the VLBW group were shorter (P < .001), but only women had lower weight (P < .001) and BMI (P = .022). Participants in the VLBW group were more likely to have ever smoked (P = .034). Among the VLBW group, 128 (57%) had received ANS, 72 (32%) were small for gestational age (SGA) (birth weight <10% percentile), and 46 (20%) were diagnosed with BPD. Self-reported respiratory symptoms are shown in Table 2, with the only statistically significant difference between the groups being a higher incidence of any wheeze in the past 12 months in the VLBW cohort.
Demographic and Perinatal Characteristics of VLBW Adults and Controls
Measure . | VLBW (n = 226) . | Control (n = 100) . | P . |
---|---|---|---|
Sex, n (%) | |||
Male | 99 (44) | 37 (37) | — |
Female | 127 (56) | 63 (63) | .25 |
Age at assessment, y | 28.5 (1.1) | 28.3 (0.9) | .13 |
Height, cm | 167.2 (8.9) | 171.2 (9.1) | <.001 |
Male | 174.5 (5.9) | 179.7 (7.0) | <.001 |
Female | 161.3 (6.1) | 166.2 (6.0) | <.001 |
Wt, kg | 75.4 (19.4) | 82.5 (18.4) | .002 |
Male | 81.2 (18.1) | 85.4 (15.8) | .21 |
Female | 70.8 (19.3) | 80.7 (19.6) | .001 |
BMI | 26.6 (6.5) | 28.2 (6.4) | .049 |
Male | 26.6 (5.5) | 26.4 (4.7) | .87 |
Female | 26.7 (7.2) | 29.3 (7.0) | .022 |
Ethnicity, n (%) | |||
Māori or Pacific Islander | 72 (32) | 24 (24) | — |
White or other | 154 (68) | 76 (76) | .15 |
Smoking status, n (%) | |||
Never smoked | 123 (54) | 67 (67) | — |
Ever smoked | 103 (46) | 33 (33) | .034 |
Current smoker | 71 (31) | 21 (21) | .054 |
Birth wt, g | 1132 (238) | 3372 (565) | <.001 |
Gestational age, wk | 29.2 (2.5) | — | — |
BPD, n (%) | 46 (20) | — | — |
Ventilation, n (%) | 147 (65) | — | — |
Ventilation, d | 16.6 (16.4) | — | — |
ANS, n (%) | 128 (57) | — | — |
SGA, n (% | 72 (32) | — | — |
RDS, n (%) | 124 (55) | — | — |
Maternal smoking in pregnancy, n (%) | 91 (41) | — | — |
Measure . | VLBW (n = 226) . | Control (n = 100) . | P . |
---|---|---|---|
Sex, n (%) | |||
Male | 99 (44) | 37 (37) | — |
Female | 127 (56) | 63 (63) | .25 |
Age at assessment, y | 28.5 (1.1) | 28.3 (0.9) | .13 |
Height, cm | 167.2 (8.9) | 171.2 (9.1) | <.001 |
Male | 174.5 (5.9) | 179.7 (7.0) | <.001 |
Female | 161.3 (6.1) | 166.2 (6.0) | <.001 |
Wt, kg | 75.4 (19.4) | 82.5 (18.4) | .002 |
Male | 81.2 (18.1) | 85.4 (15.8) | .21 |
Female | 70.8 (19.3) | 80.7 (19.6) | .001 |
BMI | 26.6 (6.5) | 28.2 (6.4) | .049 |
Male | 26.6 (5.5) | 26.4 (4.7) | .87 |
Female | 26.7 (7.2) | 29.3 (7.0) | .022 |
Ethnicity, n (%) | |||
Māori or Pacific Islander | 72 (32) | 24 (24) | — |
White or other | 154 (68) | 76 (76) | .15 |
Smoking status, n (%) | |||
Never smoked | 123 (54) | 67 (67) | — |
Ever smoked | 103 (46) | 33 (33) | .034 |
Current smoker | 71 (31) | 21 (21) | .054 |
Birth wt, g | 1132 (238) | 3372 (565) | <.001 |
Gestational age, wk | 29.2 (2.5) | — | — |
BPD, n (%) | 46 (20) | — | — |
Ventilation, n (%) | 147 (65) | — | — |
Ventilation, d | 16.6 (16.4) | — | — |
ANS, n (%) | 128 (57) | — | — |
SGA, n (% | 72 (32) | — | — |
RDS, n (%) | 124 (55) | — | — |
Maternal smoking in pregnancy, n (%) | 91 (41) | — | — |
Data are presented as mean (SD) unless otherwise indicated as n (%). RDS, respiratory distress syndrome; —, not applicable.
Self-reported Respiratory Problems in VLBW Adults and Controls
Measure . | VLBW, % (n = 224) . | Controls, % (n = 100) . | P . |
---|---|---|---|
Any wheeze past 12 mo | 34.8 | 23.0 | .034 |
Asthma ever | 43.3 | 33.0 | .08 |
Asthma past 12 mo | 7.6 | 7.0 | .85 |
Current medication for asthma | 10.7 | 13.0 | .55 |
Inhaled medicines to help breathing past 12 mo | 16.5 | 18.0 | .74 |
Troubled by shortness of breath when hurrying on level ground or walking up a slight hill | 16.4 | 9.2 | .09 |
Measure . | VLBW, % (n = 224) . | Controls, % (n = 100) . | P . |
---|---|---|---|
Any wheeze past 12 mo | 34.8 | 23.0 | .034 |
Asthma ever | 43.3 | 33.0 | .08 |
Asthma past 12 mo | 7.6 | 7.0 | .85 |
Current medication for asthma | 10.7 | 13.0 | .55 |
Inhaled medicines to help breathing past 12 mo | 16.5 | 18.0 | .74 |
Troubled by shortness of breath when hurrying on level ground or walking up a slight hill | 16.4 | 9.2 | .09 |
Lung Function Between the VLBW and Control Groups
Compared with controls, VLBW adults showed significantly lower z scores for forced expiratory volume in 1 second (FEV1), FEV1/forced vital capacity (FVC) ratio (FEV1/FVC), and forced expiratory flow at 25% to 75% of forced vital capacity (FEF25–75), whereas the difference in FVC z scores was not statistically significant (Table 3, Fig 1). The spirometric pattern is shown in Table 4; 63% of VLBW participants had normal spirometry versus 82% of controls (P = .0004). The predominant abnormality in VLBW participants was obstruction (35% vs 14% of controls), which was mild in 66 (30%) and more severe in 13 (5.5%) compared with 2% of control participants.
Pulmonary Function Comparisons Between VLBW Adults and Controls
Measure . | VLBW . | Control . | Mean Difference (95% CI)a . | P . | Adjusted Mean Difference (95% CI)b . | P . |
---|---|---|---|---|---|---|
Spirometry | ||||||
n | 224 | 100 | — | — | — | — |
FEV1 z score | −0.67 (1.20) | −0.13 (1.17) | −0.54 (−0.83 to −0.26) | .0002c | −0.56 (−0.84 to −0.27) | .0002c |
FVC z score | 0.18 (0.99) | 0.31 (1.03) | −0.13 (−0.36 to 0.11) | .30 | −0.15 (−0.39 to 0.09) | .23 |
FEV1/FVC z score | −1.21 (1.13) | 0.65 (1.01) | −0.56 (−0.82 to −0.30) | <.0001c | −0.54 (−0.80 to −0.27) | <.0001c |
FEF25–75 z score | −1.29 (1.28) | −0.52 (1.18) | −0.77 (−1.07 to −0.47) | <.0001c | −0.76 (−1.06 to −0.47) | <.0001c |
Lung volumes | ||||||
n | 226 | 100 | — | — | — | — |
TLC z score | 0.01 (1.01) | −0.01 (1.00) | 0.03 (−0.21 to 0.26) | .82 | 0.03 (−0.37 to 0.14) | .83 |
FRC z score | −0.56 (1.10) | −0.73 (1.16) | 0.18 (−0.09 to 0.44) | .19 | 0.13 (−0.13 to 0.39) | .32 |
RV z score | −0.80 (0.98) | −1.22 (0.90) | 0.42 (0.20 to 0.65) | .0003c | 0.42 (0.19 to 0.64) | .0004c |
RV/TLC z score | −1.02 (0.98) | −1.49 (0.83) | 0.46 (0.23 to 0.68) | <.0001c | 0.45 (0.22 to 0.67) | .0001c |
FRC/TLC z score | −0.67 (1.22) | −1.00 (1.22) | 0.33 (0.04 to 0.61) | .026 | 0.33 (0.04 to 0.62) | .026 |
DLCO | ||||||
n | 224 | 97 | — | — | — | — |
DLCO z score | −0.71 (1.01) | −0.11 (0.79) | −0.61 (−0.83 to −0.38) | <.0001c | −0.53 (−0.75 to −0.31) | <.0001c |
VA z score | −0.08 (1.01) | −0.14 (0.96) | 0.06 (−0.18 to 0.30) | .61 | 0.03 (−0.21 to 0.26) | .82 |
KCO z score | −0.61 (1.02) | −0.03 (0.98) | −0.58 (−0.82 to −0.34) | <.0001c | −0.49 (−0.72 to −0.26) | <.0001c |
SBN2 | ||||||
n | 201 | 95 | — | — | — | — |
Phase III slope (% N2/L) | 1.13 (0.62) | 0.93 (0.41) | 0.20 (0.07 to 0.34) | .004c | 0.21 (0.07 to 0.34) | .002c |
Measure . | VLBW . | Control . | Mean Difference (95% CI)a . | P . | Adjusted Mean Difference (95% CI)b . | P . |
---|---|---|---|---|---|---|
Spirometry | ||||||
n | 224 | 100 | — | — | — | — |
FEV1 z score | −0.67 (1.20) | −0.13 (1.17) | −0.54 (−0.83 to −0.26) | .0002c | −0.56 (−0.84 to −0.27) | .0002c |
FVC z score | 0.18 (0.99) | 0.31 (1.03) | −0.13 (−0.36 to 0.11) | .30 | −0.15 (−0.39 to 0.09) | .23 |
FEV1/FVC z score | −1.21 (1.13) | 0.65 (1.01) | −0.56 (−0.82 to −0.30) | <.0001c | −0.54 (−0.80 to −0.27) | <.0001c |
FEF25–75 z score | −1.29 (1.28) | −0.52 (1.18) | −0.77 (−1.07 to −0.47) | <.0001c | −0.76 (−1.06 to −0.47) | <.0001c |
Lung volumes | ||||||
n | 226 | 100 | — | — | — | — |
TLC z score | 0.01 (1.01) | −0.01 (1.00) | 0.03 (−0.21 to 0.26) | .82 | 0.03 (−0.37 to 0.14) | .83 |
FRC z score | −0.56 (1.10) | −0.73 (1.16) | 0.18 (−0.09 to 0.44) | .19 | 0.13 (−0.13 to 0.39) | .32 |
RV z score | −0.80 (0.98) | −1.22 (0.90) | 0.42 (0.20 to 0.65) | .0003c | 0.42 (0.19 to 0.64) | .0004c |
RV/TLC z score | −1.02 (0.98) | −1.49 (0.83) | 0.46 (0.23 to 0.68) | <.0001c | 0.45 (0.22 to 0.67) | .0001c |
FRC/TLC z score | −0.67 (1.22) | −1.00 (1.22) | 0.33 (0.04 to 0.61) | .026 | 0.33 (0.04 to 0.62) | .026 |
DLCO | ||||||
n | 224 | 97 | — | — | — | — |
DLCO z score | −0.71 (1.01) | −0.11 (0.79) | −0.61 (−0.83 to −0.38) | <.0001c | −0.53 (−0.75 to −0.31) | <.0001c |
VA z score | −0.08 (1.01) | −0.14 (0.96) | 0.06 (−0.18 to 0.30) | .61 | 0.03 (−0.21 to 0.26) | .82 |
KCO z score | −0.61 (1.02) | −0.03 (0.98) | −0.58 (−0.82 to −0.34) | <.0001c | −0.49 (−0.72 to −0.26) | <.0001c |
SBN2 | ||||||
n | 201 | 95 | — | — | — | — |
Phase III slope (% N2/L) | 1.13 (0.62) | 0.93 (0.41) | 0.20 (0.07 to 0.34) | .004c | 0.21 (0.07 to 0.34) | .002c |
Data are presented as mean (SD) and mean difference (95% confidence interval). Severity of obstruction was determined on the basis of FEV1 % predicted value. —, not applicable.
Mean VLBW – controls difference.
Mean difference adjusted for sex and ever smoker.
Mean difference is statistically significant after correction for multiple comparisons (Bonferroni corrected P = .009).
Pulmonary function test results in VLBW adults and controls. The solid horizontal bars indicate the mean values. For derivation of z scores, see the main text.
Pulmonary function test results in VLBW adults and controls. The solid horizontal bars indicate the mean values. For derivation of z scores, see the main text.
Spirometric Pattern in VLBW Adults and Controls
Spirometric pattern . | VLBW, n (%) . | Controls, n (%) . |
---|---|---|
Normal spirometry | 140 (63) | 82 (82) |
Restrictive pattern | 5 (2) | 4 (4) |
Obstructive pattern | 79 (35) | 14 (14) |
Mild obstruction | 66 (30) | 12 (12) |
Moderate obstruction | 7 (3) | 1 (1) |
Moderately severe obstruction | 5 (2) | — |
Severe obstruction | 1 (0.5) | 1 (1) |
Spirometric pattern . | VLBW, n (%) . | Controls, n (%) . |
---|---|---|
Normal spirometry | 140 (63) | 82 (82) |
Restrictive pattern | 5 (2) | 4 (4) |
Obstructive pattern | 79 (35) | 14 (14) |
Mild obstruction | 66 (30) | 12 (12) |
Moderate obstruction | 7 (3) | 1 (1) |
Moderately severe obstruction | 5 (2) | — |
Severe obstruction | 1 (0.5) | 1 (1) |
Obstructive pattern: FEV1/FVC below the lower limit of normal (eg, below the fifth percentile of the predicted value; z score < −1.64); severity of obstruction was determined on the basis of the FEV1 % predicted value; mild: FEV1 ≥70%; moderate: 60% to 69%; moderately severe: 50% to 59%; severe: 35% to 49%; very severe: <35%; restrictive pattern: normal FEV1/FVC with FVC below the lower limit of normal.16 —, not applicable.
For lung volumes, there were no significant differences in the z scores for total lung capacity (TLC) or functional residual capacity (FRC) between the groups. However, higher z scores in residual volume (RV) and RV/TLC ratio were found in the VLBW group compared with controls (Table 3, Fig 1).
The VLBW group had significantly lower z scores for DLCO and transfer coefficient for carbon monoxide (KCO) than controls. The z scores for alveolar volume (VA) were not significantly different between the groups. The phase III slope in SBN2 in the VLBW group was significantly higher than controls. All of these lung function differences persisted when adjusted for sex and smoking status (Table 3, Fig 1). Between-group differences in lung function revealed the same statistical significance when variables were analyzed as percent predicted values (Supplemental Table 7).
Lung Function in Those With and Without BPD Within the VLBW Group
Within the VLBW group, z scores for FEV1, FEV1/FVC, and FEF25–75 were all significantly lower in subjects with BPD versus subjects without BPD (Table 5). There was no between-group difference in z score for TLC, FRC, or FRC/TLC after correction for multiple comparisons, but RV and RV/TLC were higher in subjects with BPD than in subjects without BPD. DLCO was lower in subjects with BPD but not statistically significant on Bonferroni correction; VA and KCO were not significantly different. The phase III slope was significantly higher in subjects with BPD than in subjects without BPD. These effects also persisted when adjusted for sex, smoking status, birth weight (<1000 g vs ≥1000 g), gestation (<28 weeks vs ≥28 weeks), SGA status, and maternal smoking during pregnancy. The patterns remained the same when data were analyzed as percent predicted values (Supplemental Table 8). Duration of ventilation was significantly correlated with BPD (r = 0.51). Comparing ventilated and nonventilated subjects revealed broadly similar results to the BPD and non-BPD comparisons (Supplemental Table 9).
Pulmonary Function Comparisons for VLBW Adults With and Without BPD
Measure . | VLBW BPD . | VLBW Non-BPD . | Mean Difference (95% CI)a . | P . | Adjusted Mean Difference (95% CI)b . | P . |
---|---|---|---|---|---|---|
Spirometry | ||||||
n | 46 | 178 | — | — | — | — |
FEV1 z score | −1.34 (1.41) | −0.50 (1.08) | −0.85 (−1.23 to −0.47) | <.0001c | −0.75 (−1.15 to −0.34) | .0003c |
FVC z score | −0.14 (1.10) | 0.26 (0.95) | −0.40 (−0.72 to −0.08) | .014 | −0.33 (−0.68 to 0.01) | .06 |
FEV1/FVC z score | −1.67 (1.28) | −1.09 (1.06) | −0.58 (−0.94 to −0.22) | .0016c | −0.55 (−0.93 to −0.17) | .005c |
FEF25–75 z score | −1.90 (1.45) | −1.13 (1.19) | −0.77 (−1.17 to −0.36) | .0002c | −0.69 (−1.12 to −0.26) | .002c |
Lung volumes | ||||||
n | 46 | 180 | — | — | — | — |
TLC z score | −0.01 (1.01) | 0.02 (1.01) | −0.03 (−0.35 to 0.30) | .88 | 0.04 (−0.31 to 0.38) | .83 |
FRC z score | −0.33 (1.18) | −0.62 (1.08) | 0.29 (−0.07 to 0.65) | .11 | 0.43 (0.05 to 0.79) | .027 |
RV z score | −0.33 (1.23) | −0.92 (0.87) | 0.59 (0.29 to 0.90) | .0002c | 0.58 (0.25 to 0.91) | .0006c |
RV/TLC z score | −0.49 (1.23) | −1.17 (0.86) | 0.68 (0.37 to 0.99) | <.0001c | 0.63 (0.30 to 0.96) | .0002c |
FRC/TLC z score | −0.33 (1.40) | −0.76 (1.16) | 0.43 (0.03 to 0.82) | .034 | 0.55 (0.13 to 0.96) | .01 |
DLCO | ||||||
n | 46 | 178 | — | — | — | — |
DLCO z score | −0.98 (1.06) | −0.64 (0.99) | −0.33 (−0.66 to −0.01) | .045 | −0.26 (−0.59 to 0.07) | .13 |
VA z score | −0.28 (1.01) | −0.02 (1.01) | −0.26 (−0.59 to 0.07) | .12 | 0.23 (−0.58 to 0.11) | .18 |
KCO z score | −0.66 (1.06) | −0.60 (1.01) | −0.07 (−0.40 to 0.27) | .69 | −0.02 (−0.35 to 0.32) | .93 |
SBN2 | ||||||
n | 37 | 164 | — | — | — | — |
Phase III slope (% N2/L) | 1.44 (0.81) | 1.06 (0.54) | 0.38 (0.16 to 0.59) | .0006c | 0.39 (0.17 to 0.60) | .0005c |
Measure . | VLBW BPD . | VLBW Non-BPD . | Mean Difference (95% CI)a . | P . | Adjusted Mean Difference (95% CI)b . | P . |
---|---|---|---|---|---|---|
Spirometry | ||||||
n | 46 | 178 | — | — | — | — |
FEV1 z score | −1.34 (1.41) | −0.50 (1.08) | −0.85 (−1.23 to −0.47) | <.0001c | −0.75 (−1.15 to −0.34) | .0003c |
FVC z score | −0.14 (1.10) | 0.26 (0.95) | −0.40 (−0.72 to −0.08) | .014 | −0.33 (−0.68 to 0.01) | .06 |
FEV1/FVC z score | −1.67 (1.28) | −1.09 (1.06) | −0.58 (−0.94 to −0.22) | .0016c | −0.55 (−0.93 to −0.17) | .005c |
FEF25–75 z score | −1.90 (1.45) | −1.13 (1.19) | −0.77 (−1.17 to −0.36) | .0002c | −0.69 (−1.12 to −0.26) | .002c |
Lung volumes | ||||||
n | 46 | 180 | — | — | — | — |
TLC z score | −0.01 (1.01) | 0.02 (1.01) | −0.03 (−0.35 to 0.30) | .88 | 0.04 (−0.31 to 0.38) | .83 |
FRC z score | −0.33 (1.18) | −0.62 (1.08) | 0.29 (−0.07 to 0.65) | .11 | 0.43 (0.05 to 0.79) | .027 |
RV z score | −0.33 (1.23) | −0.92 (0.87) | 0.59 (0.29 to 0.90) | .0002c | 0.58 (0.25 to 0.91) | .0006c |
RV/TLC z score | −0.49 (1.23) | −1.17 (0.86) | 0.68 (0.37 to 0.99) | <.0001c | 0.63 (0.30 to 0.96) | .0002c |
FRC/TLC z score | −0.33 (1.40) | −0.76 (1.16) | 0.43 (0.03 to 0.82) | .034 | 0.55 (0.13 to 0.96) | .01 |
DLCO | ||||||
n | 46 | 178 | — | — | — | — |
DLCO z score | −0.98 (1.06) | −0.64 (0.99) | −0.33 (−0.66 to −0.01) | .045 | −0.26 (−0.59 to 0.07) | .13 |
VA z score | −0.28 (1.01) | −0.02 (1.01) | −0.26 (−0.59 to 0.07) | .12 | 0.23 (−0.58 to 0.11) | .18 |
KCO z score | −0.66 (1.06) | −0.60 (1.01) | −0.07 (−0.40 to 0.27) | .69 | −0.02 (−0.35 to 0.32) | .93 |
SBN2 | ||||||
n | 37 | 164 | — | — | — | — |
Phase III slope (% N2/L) | 1.44 (0.81) | 1.06 (0.54) | 0.38 (0.16 to 0.59) | .0006c | 0.39 (0.17 to 0.60) | .0005c |
Data are presented as mean (SD) and mean difference (95% confidence interval). —, not applicable.
Mean VLBW BPD – non-BPD difference.
Mean difference adjusted for sex, smoking, birth weight (<1000, ≥1000 g), gestation (<28, ≥28 wk), SGA, and maternal smoking in pregnancy.
Mean difference is statistically significant after correction for multiple comparisons (Bonferroni corrected P = .009).
Pulmonary function comparisons between SGA and non-SGA VLBW adults revealed nominally significant differences of lower FEV1/FVC and FEF25–75 z scores and higher phase III slope in the SGA group, but none remained significant on Bonferroni correction (Supplemental Table 10). Restricting the analysis to those <30 weeks’ gestation to reduce any sampling bias did not alter these results.
Differences in lung function between extremely low birth weight (ELBW) (birth weight <1000 g) subjects and other VLBW participants were small and did not reach statistical significance when adjusted for multiple comparisons (Supplemental Table 11).
There were no differences in spirometry and lung volumes between extremely preterm (EP) (gestation <28 weeks) and other VLBW participants (Supplemental Table 11). DLCO and KCO were significantly lower in EP survivors, but there was no difference in VA or phase III slope.
There was no statistically significant difference in any lung function parameter between those VLBW survivors who received ANSs and those who did not (Supplemental Table 9).
Discussion
In this comprehensive assessment of static lung function in a national cohort of VLBW survivors at age of 26 to 30 years, we found that compared with term-born controls, VLBW adults showed a higher incidence of airflow obstruction, gas trapping, reduced gas exchange, and higher ventilatory inhomogeneity. The majority of VLBW survivors did however have normal spirometry. Within the VLBW group, participants with a history of BPD had further reduced airflow, increased gas trapping, and ventilatory inhomogeneity compared with those without BPD. These findings extend and enhance our understanding of lung function in adult VLBW survivors by incorporating lung volumes and diffusing capacity, which are rarely reported, and measurement of ventilatory inhomogeneity by SBN2, which is novel.
Numerous published spirometry data,4 including a recent individual patient data meta-analysis,7 reveal ongoing expiratory airflow obstruction in survivors born VLBW or VP in childhood through early adulthood. Airflow obstruction in VP children is associated with abnormal lung structure changes, including subpleural opacities, bronchial wall thickening, and hypoattenuated lung areas, as detected by computed tomography scan.20 Our results are broadly consistent with the literature, with significantly lower z scores for parameters reflecting airflow (FEV1, FEV1/FVC, and FEF25–75) in adult VLBW survivors; and these values were further reduced in those with a history of BPD. Similar results in those with SGA compared with non-SGA participants did not remain significant after Bonferroni correction. The finding of reduced airflow with stable FVC indicates reduced airway caliber,20 and reduced FEF25–75 suggests that obstruction presents in distal small airways. In our VLBW cohort, 63% had normal spirometry; airflow obstruction was generally mild with only a small proportion having more severe airflow obstruction. We found no association between either male sex or ANS and spirometry results, in contrast to a recently published meta-analysis of late adolescence and adult cohort studies.7 Male sex is associated with increased mortality in VLBW and VP infants21 and with increased BPD in survivors,22 whereas ANS increase survival but have little impact on rates of BPD.23 The effects on survival may well mask any longer-term consequences.
Studies assessing lung volumes in former VLBW and/or VP subjects in childhood through late adolescence have provided conflicting results.20,24–27 Reports in adulthood, however, are rare. Our data revealed increased gas trapping, demonstrated by significantly elevated RV and RV/TLC ratio in VLBW adult survivors compared with controls, and in those with BPD compared with those without. These findings are physiologically consistent with our spirometry data, which reveal more airflow obstruction in these participants than in their counterparts, and with published lung volume data in VP children.27,28
Although there are several studies reporting on gas exchange (DLCO) in former VLBW and/or VP infants in childhood, with conflicting results, there have been few studies in adults.29 Narang et al30 reported preterm subjects (<38 weeks) had significantly reduced DLCO compared with term controls at 21 years of age. Landry et al31 reported preterm subjects (<37 weeks) with BPD had significantly reduced DLCO than other preterm subjects or term controls at 22 years of age. Our data revealed reduced DLCO and KCO in VLBW adults compared with controls, although significance was no longer found in those with BPD compared with no BPD after Bonferroni correction. These findings support the hypothesis that damaged alveolar compartments and pulmonary vasculature resulting in reduced gas exchange, as previously reported in VLBW children, persist into adulthood, with more evidence required to assess the importance of BPD.
Using the lung clearance index, a few researchers have investigated ventilatory inhomogeneity in preterm infants and children with inconsistent results.29 To our knowledge, this study is the first used to report SBN2 in VLBW adults. We have demonstrated evidence of increased ventilatory inhomogeneity, as determined by significantly elevated phase III slope, in VLBW survivors and in those with BPD.
Taken together with previous studies, there are several implications from our results. By the end of the third decade of life, these young adults will have already passed their peak lung function,4,8,17 and anything that may hasten respiratory function decline should be avoided. Chief among those is smoking, and it is disappointing that 31% of the VLBW cohort (and 21% of controls) were current smokers. Neonatologists and pediatricians engaged in follow-up clinics should emphasize the importance of VLBW and/or VP adults not smoking. Annual influenza immunization should be recommended8 and not just for individuals with a history of BPD. The arguments that everybody should be free of exposure to secondhand smoke and air pollution are particularly pertinent for VP adults.8,9 There might also be a case for routine lung function testing of VLBW and/or VP survivors in middle childhood with targeted messages for those identified as already having abnormal lung function. With increasing evidence of long-term health consequences after VP birth, it is important that clinicians in all settings caring for adolescents and young adults ask about gestational age8,32 ; yet a recent survey suggested that respiratory specialists did so infrequently.33 Taking an appropriate history should alert practitioners, and VP adults with respiratory symptoms should be referred for detailed lung function tests.
Although neonatal care has continued to advance with improved survival rates for VP infants compared to 1986, BPD is still a common outcome, albeit with perhaps a slightly different pathologic basis and milder manifestation in the neonatal period34 and beyond.35 Lung function tests in midchildhood in survivors with this “new BPD” revealed increased airway obstruction compared with EP survivors without BPD and with term controls, similar to follow-up studies from the presurfactant era.36 Although the longer-term outcome of these individuals is unclear, it seems probable that the results from our cohort will have relevance for contemporary neonatal care survivors.37
There is a need for ongoing monitoring of lung function and functional status as this and other similar cohorts age to assess the pattern of physiologic and functional decline. Of interest will be whether lung function decline is uniform or whether 2 distinct groups emerge. If the rate of decline in lung function is accelerated in the group identified as having abnormality in childhood, this allows targeted interventions in this group and reassurance to others. Currently, only a minority of VLBW adults in this cohort have any (generally mild) respiratory symptoms. But the literature suggests that in addition to obstructive airways disease, which is not always reversible, structural lung changes of emphysema may lead to future symptoms.8 Our finding of increased ventilatory inhomogeneity is consistent with this scenario, and there may be more information available from cardiopulmonary exercise testing. Additionally, the lung changes associated with VP birth coupled with the effects on cardiac and pulmonary vascular development could increase the risks of pulmonary hypertension in adult survivors.9,38,39
Strengths of our study include the fact that this is a national population-based cohort prospectively enrolled at birth in a single year and assessed in 1 center with a high participation rate, and is thus able to give a more complete picture of the respiratory phenotype of VP adults than hospital-based studies or those focused only on graduates with BPD. We performed a range of lung function tests using the latest GLI reference equations. Although born before the availability of surfactant, more than half of the cohort had been exposed to ANS, and both assisted ventilation and parenteral nutrition were routinely available. Weaknesses include the fact that we did not assess bronchial responsiveness to direct stimulation (eg, methacholine) or the response to bronchodilators. High-resolution computed tomography scans would be desirable to further delineate structural lung abnormalities.8 The cohort was recruited on the basis of birth weight, which was usual in the 1980s, but the inclusion of more mature infants of 32 weeks’ gestation and greater, many of whom were SGA, could have affected the overall results. Although we assessed the cohort at 7 to 8 years of age, when questionnaires included evidence of wheeze and a diagnosis of asthma (52% VLBW versus 27% contemporary controls),40 we did not conduct lung function testing at that time and therefore cannot report longitudinal data for these outcomes. Lastly, although this is a national cohort and the largest data set reported so far,7 our numbers are still relatively small.
Conclusions
The current study reveals a higher incidence of airflow obstruction, gas trapping, reduced gas exchange, and increased ventilatory inhomogeneity in adult VLBW survivors versus age-matched controls; subjects with BPD showed further worsened airflow obstruction compared with those without BPD. The findings suggest the pulmonary effects due to VLBW, although mild, persist into adulthood, and BPD is a further insult on small airway function. Our data highlight the importance of avoiding factors that may accelerate lung function decline, particularly smoking, in VLBW survivors. Routine lung function testing and ongoing monitoring are imperative in VLBW and/or VP survivors, with targeted intervention for those with abnormal lung function.
Acknowledgments
We are grateful to all participants in this study for their enthusiastic support of the assessment days. In addition, we acknowledge the support of scientific and medical staff in the Respiratory Physiology Laboratory at Christchurch Hospital. We thank the study project manager, Julia Martin, for all her hard work and organizational skills.
This article is dedicated to the memories of Dr Maureen Swanney and Josh Stanton, senior members of the research team and dearly loved respiratory physiologist colleagues, both of whom died during the course of the study.
†The late Dr Swanney and Mr Stanton were key members of the research team and contributed to writing abstract reports on partial data from the study but died before full preparation of this article.
Mr Yang and Ms Kingsford contributed to the study design, collected the data, conducted data analysis, drafted the initial manuscript, and reviewed and revised the manuscript; Prof Horwood contributed to the study design, conducted data analysis, and critically reviewed and revised the manuscript; Dr Epton contributed to the study design and data analysis and critically reviewed and revised the manuscript; Dr Swanney and Mr Stanton contributed to the study design, collected the data, conducted initial data analysis, and drafted abstract reports on partial data: both were deceased before the present manuscript was completed; and Prof Darlow conceived and designed the overall study, obtained funding, supervised data collection, and critically reviewed and revised the manuscript.
This trial has been registered with the Australian New Zealand Clinical Trials Registry (https://www.anzctr.org.au/) (identifier ACTRN12612000995875).
FUNDING: Funded by a project grant from the New Zealand Health Research Council (12-129) with additional funding from Cure Kids, a project grant from the Child Health Research Foundation (CHRF 5040), and the Canterbury District Health Board laboratory operational budget.
- ANS
antenatal corticosteroids
- BPD
bronchopulmonary dysplasia
- DLCO
single-breath diffusing capacity of the lung for carbon monoxide
- ELBW
extremely low birth weight
- EP
extremely preterm
- FEF25–75
forced expiratory flow at 25% to 75% of forced vital capacity
- FEV1
forced expiratory volume in 1 second
- FRC
functional residual capacity
- FVC
forced vital capacity
- GLI
Global Lung Function Initiative
- KCO
transfer coefficient for carbon monoxide
- RV
residual volume
- SBN2
single-breath nitrogen washout
- SGA
small for gestational age
- TLC
total lung capacity
- VA
alveolar volume
- VLBW
very low birth weight
- VP
very preterm
- % pred
percentage of the predicted
References
Competing Interests
POTENTIAL CONFLICT OF INTEREST: The authors have indicated they have no potential conflicts of interest to disclose.
FINANCIAL DISCLOSURE: The authors have indicated they have no financial relationships relevant to this article to disclose.
Comments