Background: The interrelationships between air pollution, lung function and the incidence of childhood asthma have yet to be established. A study was undertaken to determine whether lung function is associated with new onset asthma and whether this relationship varies by exposure to ambient air pollutants.
Methods: A cohort of children aged 9–10 years without asthma or wheeze at study entry were identified from the Children’s Health Study and followed for 8 years. The participants resided in 12 communities with a wide range of ambient air pollutants that were measured continuously. Spirometric testing was performed and a medical diagnosis of asthma was ascertained annually. Proportional hazard regression models were fitted to investigate the relationship between lung function at study entry and the subsequent development of asthma and to determine whether air pollutants modify these associations.
Results: The level of airway flow was associated with new onset asthma. Over the 10th–90th percentile range of forced expiratory flow over the mid-range of expiration (FEF25–75, 57.1%), the hazard ratio (HR) of new onset asthma was 0.50 (95% CI 0.35 to 0.71). This protective effect of better lung function was reduced in children exposed to higher levels of particulate matter with an aerodynamic diameter <2.5 μm (PM2.5). Over the 10th–90th percentile range of FEF25–75, the HR of new onset asthma was 0.34 (95% CI 0.21 to 0.56) in communities with low PM2.5 (<13.7 μg/m3) and 0.76 (95% CI 0.45 to 1.26) in communities with high PM2.5 (⩾13.7 μg/m3). A similar pattern was observed for forced expiratory volume in 1 s. Little variation in HR was observed for ozone.
Conclusion: Exposure to high levels of PM2.5 attenuates the protective effect of better lung function against new onset asthma.
Statistics from Altmetric.com
Asthma is a major cause of childhood morbidity and has high social and economic costs.1 Current research suggests a complex aetiological pathway for asthma,2 and there is emerging evidence that air pollution is one environmental factor that may be involved in the pathogenesis of the disease.3 Both indoor pollution4 and traffic exposure5 have been associated with an increased risk of asthma in children. In the Children’s Health Study (CHS),6 7 a longitudinal cohort study of air pollution, genetics and respiratory health, we have reported that the genetic risk factors for asthma depend upon exposure to ambient air pollution.8 For example, we found that the protective effect of a polymorphism in the promoter region of tumour necrosis factor α was attenuated by ambient ozone levels. This suggests that air pollution may overwhelm the protective effects of certain alleles and increase the risk of asthma in a subset of children.
Based on these observations, we suggest that air pollution modifies the risk of asthma associated with other constitutive factors such as lung function which may reflect an integrative index of respiratory health. The evidence supports a protective role for high lung function in asthma development. Flow rates during the first year of life have been associated with the subsequent onset of asthma and persistent wheeze in some9 10 but not all studies.11 However, if different environmental exposures such as air pollution attenuate the protective effect of better lung function in some studies, these differences might be expected to occur based on pathophysiology. One class of pollutants of interest is ambient particulate matter which has been associated with small airway remodelling (a hallmark of asthma),12 low airway flow rates13 14 and impaired lung function growth.15 These studies suggest that there may be joint effects of air pollution and lung function on asthma.
To assess the hypothesis that higher lung function is associated with reduced risk for childhood asthma, but that ambient air pollution attenuates this effect, we examined air pollution, lung function and health data from the CHS. Among children without asthma at study entry, we examined the relationship between spirometric measurements of lung volumes and airway flows and the subsequent risk of new onset asthma in communities with different levels of air pollution.
Study design and cohort
The design, methods and characteristics of the CHS have been described previously.6 7 Briefly, school children were enrolled into the study in 1993 from 12 communities in southern California (selected on the basis of different ambient pollution levels). At study entry parents or guardians of the children completed a questionnaire regarding children’s medical history and sociodemographic factors (for details see Methods in online supplement available at http://thorax.bmj.com/supplemental). This study was approved by the Institutional Review Board of University of Southern California.
To provide a targeted age range (9–10 years of age at baseline) in which to assess the effect of lung function on the risk of asthma from childhood through adolescence, the analysis was restricted to fourth grade cohorts who (1) were free of any wheezing or physician-diagnosed asthma at study entry, (2) could perform adequate lung function manoeuvres at study entry and (3) had at least one follow-up assessment (n = 2057).
New onset asthma
Children who were disease-free at baseline and reported physician-diagnosed asthma at annual follow-up testing were classified as having new onset asthma. As the follow-up occurred annually, the date of onset was assigned to the midpoint of the interval between the interview date when asthma diagnosis was first reported and the previous interview date. During the interview the children also provided information regarding recent use of any inhaled asthma medication.
Air pollution data
Ambient levels of ozone (O3), nitrogen dioxide (NO2), particulate matter with an aerodynamic diameter <10 µm (PM10) and <2.5 µm (PM2.5), acid vapour and elemental and organic carbon in each of the 12 communities were measured at air monitoring sites from 1994 onwards. Long-term mean pollutant levels (from 1994 to the end of 2003) were calculated for use in the statistical analysis. Owing to a high correlation between PM2.5, NO2, acid vapour, PM10 and elemental and organic carbon (but not ozone), these can be considered as a correlated non-ozone “package” of pollutants with a similar pattern relative to each other across the 12 communities. Communities defined as “high” or “low” based on any of the non-ozone “package” of pollutants were the same for all of the correlated pollutants (for details see Methods section in online supplement available at http://thorax.bmj.com/supplemental).
Lung function measurement
Data on children’s lung function were collected by trained field technicians who visited the schools annually. Maximal effort spirometric tests, standing height and weight were measured. Details regarding the lung function testing protocol have been published previously.6 These activities covered the period 1993–2004. Three measures of lung function were analysed: forced vital capacity (FVC), forced expiratory volume in 1 s (FEV1) and forced expiratory flow over the mid-range of expiration (FEF25–75). Sex-specific percentage predicted lung function values were calculated using linear regression (for detail see Methods section in online supplement available at http://thorax.bmj.com/supplemental).
The percentage predicted lung function values were scaled to the 10th–90th percentile range of the corresponding lung function, and also categorised into three groups based on cut-off points that provided an adequate distribution for analyses. The categories were <90%, 90–100% and >100% of predicted for FVC and FEV1 and <100%, 100–120% and >120% for FEF25–75.
Sociodemographic and medical history information
Personal information such as ethnicity, birth weight, premature birth, maternal smoking during pregnancy and allergy histories was collected at study entry. A family history of asthma was defined as asthma in any of the biological parents. Body mass index (BMI) was categorised into age- and sex-specific percentiles based on the Centers for Disease Control (CDC) BMI growth charts using 1-month age intervals.16 Participants with BMI at or above the 85th percentile were classified as overweight. These personal characteristics and household and indoor exposures (pets, pests, humidifier use and household smoking) were considered as potential effect modifiers as well as confounders in the analysis.
We fitted Cox proportional hazards models with sex- and age-specific (age defined as integer age at study entry) baseline hazards to investigate the association between new onset asthma and lung function at study entry. Initially, we fitted lung function as both categorical and continuous terms. In the absence of any non-linear association between lung function and new onset asthma (tested by adding quadratic terms for the lung functions in the continuous models as well as by fitting models with categorical lung function terms) based on likelihood ratio tests, we report results using lung function as a continuous term. The hazard ratio (HR) can be interpreted as the change in risk of new onset asthma as the lung function increases over the 10th–90th percentile range of the corresponding lung function. We also fitted proportional hazard models treating annual lung function as 1- or 2-year lagged time-dependent covariates. All models were adjusted for community and race/ethnicity. Additional covariates were considered for inclusion in the model based on whether their inclusion changed the lung function effect estimate by more than 10%. Heterogeneity of associations among subgroups was assessed by comparing appropriate models with and without interaction terms.
To assess the effect of ambient air pollution on the relationship between lung function and new onset asthma, we estimated the heterogeneity of association using community levels of air pollutants measured at one monitor in each community. To address this issue we fitted hierarchical two-stage models to these time-dependent data (for details see Methods section in online supplement available at http://thorax.bmj.com/supplemental).17
All analyses were conducted using SAS software (SAS Institute, Cary, NC, USA) Version 9.1. All hypothesis testing was conducted assuming a 0.05 significance level and a two-sided alternative hypothesis.
Most of the 2057 disease-free children included in this analysis were 10 years of age or less at study entry and non-Hispanic white (table 1). Most children had health insurance and parents with at least a high school education. The median and mean years of follow-up were 8 and 6.5 years, respectively, for the 1993 cohort and 8 and 6.2 years, respectively, for the 1996 cohort (8 years maximum possible for both cohorts). The completeness of the follow-up data was calculated as the ratio of observed and maximum possible years of follow-up. Overall, the children were followed for 79% of the possible time of observation over the 8-year period of the study. The completeness of follow-up did not vary substantially across any of the subgroups (table 1). Based on telephone interviews with the families of subjects with asthma or wheeze at study entry who left the schools, the loss to follow-up was mainly due to employment-related moves of families out of the school catchment area (data not shown). Thus, the loss to follow-up was random and not related to any exposure or disease status.
There was no difference between the children included in the study and those who were excluded owing to insufficient lung function assessments (n = 521) with respect to age, sex, in utero and environmental exposure to smoking, parental education, parental asthma status and health insurance. However, those included in the study were slightly more likely to be non-Hispanic white, to have dogs at home and a higher parental income than those with insufficient lung function assessments (see table E3 in online supplement available at http://thorax.bmj.com/supplemental).
Lung function and asthma
There were 212 new cases of asthma resulting in an incidence rate of 16.1/1000 person-years. The incidence of newly diagnosed asthma was inversely associated with measures of airflow at study entry (see table E4 in online supplement available at http://thorax.bmj.com/supplemental). For example, the incidence rate of newly diagnosed asthma increased from 9.5/1000 person-years for children with percentage predicted FEF25–75 values ⩾120% to 20.4/1000 person-years for children with FEF25–75 values ⩽100%. A similar inverse relationship was observed for FEV1 but was less clear for FVC.
Over the 10th–90th percentile range for FEF25–75 (57.1) the hazard ratio (HR) of new onset asthma was 0.50 (95% CI 0.35 to 0.71), table 2. A similar significant inverse association was observed for FEV1. The results from lagged models were very similar to the baseline analysis, except for FVC where a statistically significant inverse relationship was also observed between FVC and new onset asthma. These inverse associations were not confounded or modified by birth weight, premature birth, parental history of asthma, history of allergy, second hand smoke, current and in utero maternal smoking, personal smoking, physical activity, health insurance, overweight or by parental income or education.
Interrelationship between ambient air pollution, lung function and new onset asthma
The incidence rates across lung function categories differed markedly between “high” (13.7–29.5 μg/m3) and “low” (5.7–8.5 μg/m3) PM2.5 communities (table 3). The incidence rate of asthma for FEF25–75 ⩾120% in the “high” PM2.5 communities was 15.9/1000 person-years compared with 6.4/1000 person-years in “low” PM2.5 communities. However, there was little difference in the incidence rate of asthma for FEF25–75 80–100% in the “high” and “low” PM2.5 communities. In addition, loss of protection by high lung function against new onset asthma in the “high” PM2.5 communities was observed for all of the lung function measures. Over the 10th–90th percentile range of FEV25–75, for example, the HR of new onset asthma was 0.34 (95% CI 0.21 to 0.56) in the “low” PM2.5 communities, whereas this protective effect was reduced to 0.76 (95% CI 0.45 to 1.26) in the “high” PM2.5 communities (table 4). There were no substantial differences in the effect of lung function between “high” and “low” ozone communities.
Using a hierarchical model, we evaluated the effect of individual air pollutants (NO2, PM10, PM2.5, acid vapour, ozone and elemental and organic carbon) on the association between lung function and asthma (fig 1). The loss of the protective effect from better lung function can be appreciated from the graphs. The modifying effect of PM2.5, PM10 and organic carbon was statistically significant (p⩽0.05) and that of NO2, elemental carbon and acid vapour was marginally significant (p⩽0.08). Of all the pollutants, PM2.5 appeared to have the strongest modifying effect on the association between lung function with asthma as it had the highest R2 value (0.42). For each 57.1% (10th–90th range) change in FEF25–75, the risk of asthma decreased by 0.35-fold in a community with ambient PM2.5 levels of 10 μg/m3 whereas this protective effect was attenuated to 0.90-fold in a community with ambient PM2.5 levels of 25 μg/m3 (fig 1).
To assess the effects of potential misclassification of new onset asthma we limited the definition of incident asthma cases to those who also reported recent use of inhalers. The observed associations were similar to the primary analysis (see tables E5 and E6, model 1, in online supplement available at http://thorax.bmj.com/supplemental). Restricting the analysis to children aged 10 years or less at study entry (see tables E5 and E6, model 2, in online supplement available at http://thorax.bmj.com/supplemental) or children without extreme values of lung function (5th–95th percentile) (see tables E5 and E6, model 3, in online supplement available at http://thorax.bmj.com/supplemental) also did not alter our initial findings. The association also remained unchanged after adjusting for ambient PM2.5 (see table E5, model 4, in online supplement available at http://thorax.bmj.com/supplemental) and risk factors of low lung function after birth, ie, preterm birth, birth weight, maternal smoking during pregnancy and parental history of asthma (tables E5 and E6, model 4, in online supplement available at http://thorax.bmj.com/supplemental).
The joint effect of air pollution and lung function in the development of asthma during adolescence has not been characterised previously. The results of this study show that better airflow, characterised by higher FEF25–75 and FEV1 during childhood, was associated with a decreased risk of new onset asthma during adolescence. However, exposure to high levels of PM2.5 (or PM10, NO2, acid vapour, elemental carbon or organic carbon) attenuated this protective association of lung function on the occurrence of asthma.
Although the protective effect of better lung function on the development of asthma has been observed in other studies, the attenuation of this effect with increasing levels of particulate and related non-ozone pollutants is a new observation. The explanation of this pattern of effects is unclear. However, chronic inflammation and resulting airway remodelling is a central pathophysiological characteristic of asthma18 and, in one necroscopic study, chronic exposure to high levels of particulate air pollution was associated with small airway remodelling characterised by increase in fibrous tissue and smooth muscle in respiratory bronchioles.12 The authors concluded that those changes could result in chronic airflow obstruction. Furthermore, exposure to a high level of particulate matter and ozone was associated with radiological evidence of bronchiolar disease and mild bronchial wall thickening on CT scans in children 5–13 years of age.19 In our study it is possible that chronic inflammation in the distal airways induced by air pollution led to remodelling of the airways that modified the protection conferred by better lung function (or associated risk factors) on the subsequent development of asthma.
We have interpreted our results as a protective effect of better lung function that was attenuated by the effects of air pollution. Our data do not support an alternative hypothesis that particulate air pollution increased rates of incident asthma among children with poor lung function at study entry because rates among those with poor lung function were similar in both low and high pollution communities (table 3). This interpretation is also consistent with previous analyses in this cohort in which we found that the reduced risk of asthma associated with the TNFα-308 GG genotype (a promoter variant of tumour necrosis factor) was attenuated by ambient ozone levels.8 It is therefore likely that evolutionary selection has resulted in lung structure and other characteristics that promote better respiratory health. Better lung function may be a marker for lower susceptibility to airway pathophysiology.
Higher flow rates are well correlated with a lower prevalence of airway hyper-reactivity,20–22 a phenotype characteristic of asthma which has common genetic determinants with the disease.23 Furthermore, it has been observed that the age-related decline in the prevalence of airway hyper-reactivity is larger among those with higher FEV1%.24
Our results are consistent with the limited number of previous cohort studies that have found that higher flow rates in school children were associated with a lower risk of asthma and wheezing in adults.25 26 Because our cohort was not recruited at birth, we could not evaluate the effect of lung function very early in life on the risk of asthma. Studies that measured lung function in the first month10 or first year11 of life showed that lung function tracked throughout childhood and adolescence. However, the association between early life lung function and subsequent development of wheeze and asthma was not consistent in these studies. In a US study, lung function during the first year of life was associated with subsequent late onset or persistent wheeze but not with transient wheeze.27 An Australian study reported that lung function in the first month of life was associated with persistent wheeze but not with transient or late onset wheeze.10 It is possible that different constitutive factors such as genetic predisposition28 and environmental exposures such as maternal smoking or other air pollutants28 or sensitivity to different allergens influenced the pattern of wheeze associated with early deficits of lung function.10 29 A recent study of a Norwegian birth cohort found that better lung function within days of birth was associated with a lower prevalence of asthma (16.2%) by age 10 compared with those with worse lung function (24.3%) after birth. However, the authors concluded that lung function ascertained within a few days of birth was not a good predictor of asthma at age 10 as the positive predictive value for asthma ranged from 24.3 to 31.3.9
The strength of our study was the 8-year follow-up of children in the fourth grade and the collection of their yearly lung function, asthma diagnosis and air pollution exposure in a consistent manner. The observed incidence rate of physician-diagnosed asthma in the present study (16.1 cases/1000 person-years) was similar to the increasing occurrence of asthma in recent decades in children.30–33 The retention rate in this cohort is similar to other longitudinal studies involving children with lung function assessments,27 34 with an observed follow-up of 79% of the total possible person-years. Furthermore, the 521 children excluded from the study because of insufficient follow-up data did not differ from those included in this analysis with respect to covariates pertinent to the incidence of asthma (see table E3 in online supplement available at http://thorax.bmj.com/supplemental).
As this is not a birth cohort study, the temporal relationship between lung function and asthma diagnosis remains a potential concern. Because low lung function tracks over time and is associated with different patterns of wheeze,10 27 34 it is possible that the cases of new onset asthma in our study were undiagnosed cases of asthma and had low lung function at study entry. This seems unlikely because children with any history of wheeze or asthma as well as children without a history of wheezing at study entry were excluded (∼40% of the cohort). Furthermore, the flow rates at study entry in children who developed new onset asthma were significantly higher than those of children with early transient or early persistent wheeze (see Lung function by wheeze status in online supplement available at http://thorax.bmj.com/supplemental) who could have low lung function10 11 27 and might be misclassified as “non-asthmatic” at study entry by forgetful parents. The observed associations also remained essentially unchanged even after removing those with extreme values (see tables E5 and E6, model 3, in online supplement available at http://thorax.bmj.com/supplemental). Furthermore, if the observed association between flow rates and new onset asthma were mediated through early low lung function, then adjustment for risk factors for low lung function after birth and ambient PM2.5 should have attenuated the effects of lung flow on the risk of asthma. However, we did not observe any such effect (see tables E5 and E6, model 4, in online supplement available at http://thorax.bmj.com/supplemental). Tracking of lung function, undiagnosed pre-existing asthma and low lung function at study entry due to ambient PM2.5 level are therefore unlikely to explain our results.
Another potential limitation of our study could be the accuracy of new onset asthma as it was based on personal interview of the children; however, a recent study noted that children as young as 7 years can provide information regarding their asthma with an acceptable level of validity and reliability.35 In our analysis we addressed the possibility that some children might have pre-existing undiagnosed asthma by excluding incident asthma cases diagnosed within the first 2 years of follow-up in the lagged time-dependent analysis. The similarity in the associations between the time-dependent lagged and baseline models makes it unlikely that the observed associations were only due to misclassification of asthma status at cohort entry or during follow-up.
Another potential concern may be that the results of this study are due to chance as we have used three different lung function measurements and tested interaction with ozone and the six correlated non-ozone pollutants. However, we do not think this is the case as the choice of the lung function measurements and the pollutants are based on a priori hypothesis15 and all the results show a consistent pattern. Furthermore, the lung function parameters and the non-ozone group of pollutants are correlated so a Bonferroni type of adjustment is not appropriate in this setting.
In conclusion, our data suggest that better lung function during childhood provides protection against the development of asthma during adolescence. However, this protective effect of better lung function was attenuated by exposure to higher levels of ambient particulate matter (PM2.5 and PM10), NO2, acid vapour and elemental and organic carbon. Further research is needed to clarify the mechanism by which air pollution modifies the protective effect of better lung function on the development of asthma.
Further details are given in the online supplement available at http://thorax.bmj.com/supplemental.
We acknowledge all the participants and their parents of the CHS Study, as well as all the teachers and staff of the participating schools. Without their active participation we would not be able to address these intriguing and important aspects of respiratory health in children and adolescents. We are also grateful to the members of our experienced field team whose conscientious and courteous work ethic enabled us to gather this immense data efficiently. We also thank our colleagues who are involved in the maintenance and quality control of the database.
This work was supported by the Southern California Environmental Health Sciences Center (grant number 5P30ES007048) funded by the National Institute of Environmental Health Sciences; the Children’s Environmental Health Center (grant numbers 5P01ES009581, R826708-01 and RD831861-01) funded by the National Institute of Environmental Health Sciences and the Environmental Protection Agency; the National Institute of Environmental Health Sciences (grant number 5P01ES011627); the National Heart, Lung and Blood Institute (grant numbers 5R01HL61768 and 1R01HL76647); the California Air Resources Board (contract number 94-331); and the Hastings Foundation. The study sponsors had no role in the study design; the collection, analysis and interpretation of data; the writing of the report; or in the decision to submit the paper for publication.
Competing interests: None.
- forced expiratory flow over the mid-range of expiration
- forced expiratory volume in 1 s
- forced vital capacity
- nitrogen dioxide
- PM2.5, particulate matter with an aerodynamic diameter <10 µm and <2.5 µm
If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.