Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Effect of Temperature on Cystic Fibrosis Lung Disease and Infections: A Replicated Cohort Study

  • Joseph M. Collaco ,

    mcollac1@jhmi.edu

    Affiliation Eudowood Division of Pediatric Respiratory Sciences, Johns Hopkins University School of Medicine, Baltimore, Maryland, United States of America

  • John McGready,

    Affiliation Department of Biostatistics, Johns Hopkins Bloomberg School of Public Health, Baltimore, Maryland, United States of America

  • Deanna M. Green,

    Affiliation Eudowood Division of Pediatric Respiratory Sciences, Johns Hopkins University School of Medicine, Baltimore, Maryland, United States of America

  • Kathleen M. Naughton,

    Affiliation McKusick-Nathans Institute of Genetic Medicine, Johns Hopkins University School of Medicine, Baltimore, Maryland, United States of America

  • Christopher P. Watson,

    Affiliation McKusick-Nathans Institute of Genetic Medicine, Johns Hopkins University School of Medicine, Baltimore, Maryland, United States of America

  • Timothy Shields,

    Affiliation Department of Molecular Microbiology and Immunology, Johns Hopkins Bloomberg School of Public Health, Baltimore, Maryland, United States of America

  • Scott C. Bell,

    Affiliation The Prince Charles Hospital and University of Queensland, Brisbane, Queensland, Australia

  • Claire E. Wainwright,

    Affiliation Royal Children's Hospital and University of Queensland, Brisbane, Queensland, Australia

  • for the ACFBAL Study Group, Garry R. Cutting3

    Affiliation McKusick-Nathans Institute of Genetic Medicine, Johns Hopkins University School of Medicine, Baltimore, Maryland, United States of America

Abstract

Background

Progressive lung disease accounts for the majority of morbidity and mortality observed in cystic fibrosis (CF). Beyond secondhand smoke exposure and socio-economic status, the effect of specific environmental factors on CF lung function is largely unknown.

Methods

Multivariate regression was used to assess correlation between specific environmental factors, the presence of pulmonary pathogens, and variation in lung function using subjects enrolled in the U.S. CF Twin and Sibling Study (CFTSS: n = 1378). Significant associations were tested for replication in the U.S. CF Foundation Patient Registry (CFF: n = 16439), the Australian CF Data Registry (ACFDR: n = 1801), and prospectively ascertained subjects from Australia/New Zealand (ACFBAL: n = 167).

Results

In CFTSS subjects, the presence of Pseudomonas aeruginosa (OR = 1.06 per °F; p<0.001) was associated with warmer annual ambient temperatures. This finding was independently replicated in the CFF (1.02; p<0.001), ACFDR (1.05; p = 0.002), and ACFBAL (1.09; p = 0.003) subjects. Warmer temperatures (−0.34 points per °F; p = 0.005) and public insurance (−6.43 points; p<0.001) were associated with lower lung function in the CFTSS subjects. These findings were replicated in the CFF subjects (temperature: −0.31; p<0.001; insurance: −9.11; p<0.001) and similar in the ACFDR subjects (temperature: −0.23; p = 0.057). The association between temperature and lung function was minimally influenced by P. aeruginosa. Similarly, the association between temperature and P. aeruginosa was largely independent of lung function.

Conclusions

Ambient temperature is associated with prevalence of P. aeruginosa and lung function in four independent samples of CF patients from two continents.

Introduction

Cystic fibrosis (CF) is an autosomal recessive disorder caused by mutations in the CFTR gene. Progressive obstructive lung disease and recurrent respiratory infections account for the majority of morbidity and mortality associated with CF. Significant variation in CF lung disease exists, even among individuals with identical mutations [1][3], and approximately half of this variation is attributable to environmental and/or stochastic (random) factors [4], [5]. To improve outcomes for this life-limiting disorder, it is important to identify specific environmental factors that ameliorate or exacerbate disease. Prior CF studies have demonstrated associations between lung function and local environmental factors, including secondhand smoke exposure [6][10], air pollution [11], household income [12][14], maternal education [13], and insurance status [13], [15][17]. However, with the exception of air pollution [11], there has been limited work examining environmental factors operating on a geographic scale, such as climate.

Using data from the U.S. CF Twin and Sibling Study, we sought to determine (i) whether selected geographic factors (factors mapped using geospatial techniques) were associated with the presence of significant respiratory pathogens in CF, (ii) whether geographic factors were associated with variation in lung function, and (iii) the relative impact of geographic factors compared to other environmental factors (demographic or household) in multivariate analyses. We attempted to replicate key findings in subjects from two national CF patient registries (United States and Australia). The role of infection was also examined in a prospective sample of CF patients from Australia and New Zealand.

Methods

Ethics Statement

Written informed consent was obtained from CFTSS subjects. This study, including the CFF, ACFDR, and ACFBAL data downloads, were specifically approved by the Johns Hopkins University Institutional Review Board (NA_00035659, NA_00019677). Participants: U.S. CF Twin and Sibling Study (CFTSS): Subjects were recruited from CF centers based on having a twin or sibling also with CF (n = 1658 in 817 families) [18]. Data were collected between 10/27/00 and 9/25/09 with data supplementation from the U.S. CF Foundation Patient Registry through 12/31/08. U.S. CF Foundation Patient Registry (CFF): Anonymized data from the calendar year 2007 was provided (n = 24,799). Australian Cystic Fibrosis Data Registry (ACFDR): Anonymized data from 1998–2008 was provided (n = 3789). Australian Cystic Fibrosis Bronchoalveolar Lavage Trial (ACFBAL): Anonymized data was provided for 168 infants diagnosed with CF in Australia and New Zealand between 6/10/99 and 1/18/05 who were prospectively followed for microbial acquisition until 5 years of age. Subjects were excluded if lung function (CFTSS, CFF, ACFDR), respiratory culture, or residential postal/zip code data were not available, or if actively smoking (CFTSS, CFF)(Figures S1 and S2). CFTSS subjects enrolled in the CFF Patient Registry (n = 1435) were excluded from the CFF sample.

Household Variables

Secondhand smoke exposure was defined as any reported home exposure (Table 1). Maternal education was characterized as 1 = less than a high school education; 2 = high school graduate; 3 = some college education; 4 = college graduate. Household income was estimated using zip code and 2000 U.S. Census data. Insurance status was determined from the most recent year of data. Household density was the number of persons residing within the household. Household data were not available for Australian subjects.

Geographic Variables

U.S. subjects were mapped to the center of their most recent known residential zip code (finest resolution available), to derive measures for humidity, temperature, air pollution (Fine particulate matter: PM2.5), elevation, and distance to the CF care center of record using ArcGIS 9.3 (ESRI; Redlands, CA) and ERDAS Imagine 9.3 (ERDAS; Atlanta, GA). Source data were selected on the basis the time period available that matched the clinical data and included 1961–1990 mean monthly and annual relative humidity and temperatures (The Climate Source; Corvallis, OR; 2 km resolution), 2006 PM2.5 measures (U.S. Environmental Protection Agency), elevations (U.S. Geological Survey), and population densities (2000 U.S. Census). PM2.5 measures were derived for only subjects living within 30 miles of a pollution monitor based on previous work by Goss et al. [11]. Australian subjects were mapped to the center of their residential postal code to derive measures for mean annual temperature using System for Automated Geoscientific Analyses (SAGA User Group Association; Hamburg, Germany), GDAL libraries and Python bindings (Open Source Geospatial Foundation; Vancouver, BC), and Shapely v1.2 (Python Software Foundation; Wolfeboro Falls, NH). Source data included 1961–1990 and 1971–2000 mean annual temperatures (Australian Bureau of Meteorology and the [New Zealand] National Climate Database, respectively).

Outcome Variables

Raw FEV1 (liters) measurements were converted into CF-specific percentiles [19], excluding measurements obtained before 6 years of age and after lung transplantation. Lung function was defined as the best percentile in the most recent year of data. The higher lung function measured in the CFTSS subjects vs. other subjects likely reflects the younger age of CFTSS subjects (Table 1; Table S1). The lower lung function observed in the ACFDR subjects may represent the limitation of applying of a U.S.-based phenotype to an Australian population. Using all available respiratory culture data, subjects were considered to be positive for an organism (Pseudomonas aeruginosa, mucoid P. aeruginosa, or Burkholderia cepacia complex) if they had any cultures positive for that organism. The higher prevalence of infectious organisms in the CFTSS and ACFDR subjects vs. the CFF subjects reflects the multi-year culture data compared to the single year data, respectively (Table 1). Age of acquisition for an organism was defined as the date of the first positive culture following at least one prior negative culture (CFTSS, ACFBAL) [20]. Means of ascertainment likely account for age of acquisition differences between the retrospective CFTSS and prospective ACFBAL samples.

Data Analysis

Regressions clustered by family (CFTSS) or by CF care center (CFF), ANOVA, chi square, and student's t-tests were performed using Stata 10 (StataCorp LP; College Station, TX). Center-level data were not available for Australian subjects. CFTSS subjects served as the primary population being the best characterized of the study populations. The CFTSS preliminary multivariate (logistic for presence of pathogens; linear for lung function) regression models included all predictor variables significant in univariate modeling. The final CFTSS models were generated by dropping non-significant (p≥0.05) predictors from the preliminary models in a stepwise manner. The final CFTSS models were used in replication regressions for the other samples. For U.S. and Australian subjects, temperature quartiles were derived from all CFTSS subjects (n = 1557) and ACFDR subjects (n = 3635) with temperature data, respectively. Kaplan-Meier plots were constructed for the age of acquisition for P. aeruginosa by temperature quartiles for the CFTSS and ACFBAL subjects. The effect of mediation by infection on lung function was estimated by dividing the difference of regression coefficients of temperature adjusted and unadjusted for the mediator by coefficient of temperature unadjusted for the mediator [21], [22]. The mediation effect by lung function on infection was estimated similarly using the logarithms of the temperature odds ratios.

Results

Mean Annual Temperature is Associated with Prevalence of and Age at P. Aeruginosa Infection

Given the importance of bacterial infection in CF lung function decline [23], [24], we tested whether selected environmental factors were associated with specific organisms using multivariate logistic regression. Among CFTSS subjects, ambient temperature (Table 2: OR = 1.06 per °F; p<0.001) was associated with the prevalence of P. aeruginosa. CFTR genotype (based on the number of the most common mutation (F508del) present, OR = 2.05 per F508del mutation; p<0.001), age at the last respiratory culture (OR = 1.19 per year; p<0.001), and age at diagnosis (OR = 0.84 per year; p<0.001) were also associated with P. aeruginosa prevalence, as previously reported (Complete regression results in Table S2) [3], [20], [25]. In contrast to previous studies [26], [27], none of the household factors tested were associated with the presence of P. aeruginosa. The observed association was not accounted for by increased culturing frequency in warmer regions, as a higher frequency of cultures was associated with colder temperatures (Regression coefficient p<0.001; n = 1370). Also, the association was not accounted for by a temperature-humidity interaction (Interaction term p = 0.45).

thumbnail
Table 2. Annual Ambient Temperature is a Predictor of P. aeruginosa Infection.1

https://doi.org/10.1371/journal.pone.0027784.t002

The association of warmer temperatures and an increased P. aeruginosa prevalence was also seen in the CFF subjects (OR = 1.02; p<0.001). As both the CFTSS and CFF subjects are U.S.-based, this association could be subject to regional biases in socio-economic status, culture, or clinical care patterns. Australian subjects are ideal for replication of geographic modifiers of health outcomes in CF as, like the U.S., the continent has a European immigration-derived CF population and a wide range of climatic conditions, including temperature. The association of warmer temperatures and the prevalence of P. aeruginosa was replicated in both the ACFDR (OR: 1.05; p = 0.002) and ACFBAL subjects (OR: 1.09; p = 0.003). Only older age was associated with the presence of mucoid P. aeruginosa (Table S3: CFTSS, CFF) or B. cepacia complex (Table S4: CFTSS, CFF, ACFDR), which corresponds to age-specific rates of prevalence reported in CF [24].

To examine whether the relationship between temperature and P. aeruginosa was non-linear (as regression analysis assumes linearity), we examined the prevalence of P. aeruginosa by temperature quartiles (Table S5). In all four samples of CF subjects, the prevalence of P. aeruginosa was higher in warmer temperature quartiles (p values<0.001–0.005) paralleling the results of the regression analysis. In addition, age at acquisition of P. aeruginosa in the retrospective CFTSS and prospective ACFBAL samples were 15 and 9 months earlier, respectively, in the warmest temperature quartile compared to the coldest (CFTSS p = 0.04; Figure S3: log rank test p = 0.006; ACFBAL p = 0.04; log rank test p = 0.001). There were no replicated associations between temperature quartile and prevalence or age of acquisition of mucoid P. aeruginosa or B. cepacia complex.

Mean Annual Temperature is Associated with Lung Function

We examined the role of environmental factors on lung function among CFTSS subjects and found that warmer temperatures (−0.34 per °F; p = 0.005) were associated with lower lung function (Table 3; complete regression results in Table S1). Although the lung function measure used does account for age [19], older age was associated with lower lung function (−0.40 per year; p<0.001), likely due to cohort or survival effect. Public insurance (vs. private insurance) was also associated with lower lung function (−6.43; p<0.001), consistent with prior studies [13], [15], [16].

thumbnail
Table 3. Annual Ambient Temperature and Insurance Status are Predictors of Lung Function (CF-specific FEV1).1

https://doi.org/10.1371/journal.pone.0027784.t003

Temperature also was associated with lung function in the CFF subjects (−0.31; p<0.001), with a similar magnitude to the CFTSS (−0.34). This corresponds to ∼3 percentile point drop in FEV1 for each 10°F increase in mean annual temperature. In addition, warmer temperatures tended to be associated with lower lung function in the ACFDR subjects (−0.23; p = 0.057). Mean annual temperatures in Australia are higher than in the U.S., which may influence the magnitude of the coefficient and its significance. The lower lung function observed in the ACFDR subjects could be a function of overall warmer temperatures in Australia. To assess whether this association was independent of CFTR genotype, we examined white patients in the largest sample of CF subjects (CFF) who had identical CFTR genotypes (F508del homozygotes) and found that warmer temperatures remained associated with lower lung function with a similar co-efficient (−0.35; n = 6367; p<0.001) to that from the entire CFF sample (−0.31).

Seasonal Difference in Temperature is also associated with Lung Function

To examine the possibility of seasonal effects, lung function data in both January and July was examined in both U.S. samples of subjects (Table 4). In the larger CFF sample, lung function was higher in January than in July (n = 2145; p<0.001) and there was a trend towards higher lung function in January in all four quartiles. Thus, there is a suggestion that individual patients have higher lung function in January than in July, which may reflect temperature differences between the two seasons. Of note, lung function was higher in colder temperature quartiles in both CFTSS and CFF subjects, regardless of whether sampled in January or July, suggesting the predominance of average annual temperatures over seasonal fluctuations.

thumbnail
Table 4. Comparison of Lung Function (CF-specific FEV1) by Temperature Quartile and Seasonal Extremes.

https://doi.org/10.1371/journal.pone.0027784.t004

Temperature may act on Lung Function and P. aeruginosa through Independent Mechanisms

Having observed temperature to be associated with lung function and P. aeruginosa, we sought to determine whether temperature acts through the same or different mechanisms on these outcomes. One pathway might be that temperature affects infection rates, which in turn alters lung function. To test this, we included the presence of P. aeruginosa as a predictor of lung function in a mediation analysis, and found that the association between temperature and lung function decreased (Table S6: CFTSS: −0.34 to −0.29 per °F; CFF: −0.31 to −0.29; ACFDR: −0.23 to −0.18). Using these regression co-efficients, infection with P. aeruginosa is estimated to account for 15% of the association between temperature and lung function among CFTSS subjects (e.g., (−0.34 – −0.29)/ −0.34), 6% among CFF subjects and 22% among ACFDR subjects. The alternate pathway is that temperature affects lung function, which in turn alters the acquisition of infection. For this pathway, lung function is estimated to account for 7% of the association between temperature and P. aeruginosa among CFTSS subjects (e.g., (ln(1.057456)-ln(1.053488))/ln(1.057456)), 20% among CFF subjects, and 6% among ACFDR subjects (Table S7). Thus, there is limited overlap between the effects of temperature on the two clinical outcomes, suggesting that mechanisms through which temperature may act on lung function vs. on P. aeruginosa infection are largely different.

Discussion

Ambient temperature has been demonstrated to affect the prevalence of both infectious and non-infectious diseases [28]. Our analysis demonstrates that CF patients living in areas with warmer annual temperatures have a higher prevalence of and an earlier age of acquisition of P. aeruginosa. Given a constant temperature (37°C) within the airways, the association of P. aeruginosa and temperature is likely mediated outside the host. Our findings may be due to increased prevalence of P. aeruginosa in the environment secondary to more favorable conditions for the organism in warmer climates. Regional P. aeruginosa biodiversity may also play a role as genotypes of soil isolates of P. aeruginosa differ by geographic distance [29]. Genetic differences that alter the adhesion or biofilm properties of P. aeruginosa may interact with local environmental features to affect the likelihood of acquisition by a CF patient [30][32]. Indeed, higher ambient temperatures (30°C vs. 15°C) have been shown to alter the capacity of P. aeruginosa to adhere to surfaces [33]. It is possible that regional differences in neonatal screening, and P. aeruginosa detection and prevention may lead to differences in the longitudinal development of infection [27], [34].

Regarding lung function, we estimate that patients residing in the warmest regions in the U.S. would have CF-specific lung function 10 percentile points lower than if they had resided in the coldest regions of the U.S. where mean temperatures are ∼30 degrees (°F) lower. Translated to more clinically familiar measures (NHANES FEV1 percentages) [35], a hypothetical 18 year old white male with CF (Height: 175cm) with an FEV1 of 73.5% percent living in a cold climate would be expected to have an FEV1 of 66.1% had he resided in a 30 degree (°F) warmer climate. Furthermore, lung function is likely associated with season, as patients tend to have higher lung function in the colder month of January compared to July. Although temperature is associated with lung function, based on the univariate regression r-value (r = 0.0069), temperature accounts for only 0.7% of the variation seen in CF lung function. Also, the mechanisms through which temperature acts on lung function are unclear, but may include infectious agents, aeroallergens, air pollution, and socio-economic status.

Our study supports the concept that infectious agents may mediate the association between temperature and lung function, at least through P. aeruginosa, which accounts for 6–22% of the association between temperature and lung function. Further investigation into viruses and other CF pathogens may reveal other mediators. Temperature is a predictor of pollen loads as well as mold production [36], although differences in lung function have not been seen with allergy sensitive in CF (n = 55) [37]. Higher temperatures may exacerbate the effects of air pollution [38], [39], and air pollution has already been demonstrated to be associated with lower CF lung function [11]; thus, temperature and air pollution may interact to worsen CF lung function. Geographic variation in temperature does not appear to serve as a proxy for socio-economic status in our study as the association between lung function and temperature remained robust after adjusting for significant socio-economic factors (Table S1).

Study limitations include that subjects may reside in multiple locations over a lifetime. Lung function from the most recent year of data and each subject's last known postal or zip code were used to minimize this uncertainty. Furthermore, for the U.S. subjects, 79.5% (CFTSS) and 77.6% (CFF) live within the state in which they were born with another 8% (CFTSS) and 7.3% (CFF) living within an adjacent state. A limitation of cross-sectional data is the uncertainty in measuring infection status, and hence accurate ascertainment of infection prevalence. To address this issue, longitudinal respiratory culture data were obtained in the CFTSS and ACFDR samples; it should also be noted that this longitudinal assessment may lead to temporal effects as cultures and lung function data may have been obtained for different subjects as much as 10 years apart. Another limitation is our largely retrospective culture ascertainment, thus we sought out a prospectively ascertained sample (ACFBAL). There are other confounding factors that we were unable to assess, such as co-morbidities, family support, physical activity, etc. Finally, ACFDR lung function results are subject to the caveat of applying a U.S-based FEV1 phenotype to an Australian population.

Our findings suggest that accounting for temperature should be considered in the design of both epidemiological studies of infection and/or lung function and clinical trials that encompass broad geographic areas. As the effects of temperature upon lung function and P. aeruginosa are not within patients' and clinical providers' control, geography may need to be considered when comparing the performance of CF Care Centers.

Supporting Information

Figure S1.

Derivation of Cystic Fibrosis Twin-Sibling Study (CFTSS) sample and Cystic Fibrosis Foundation Patient Registry (CFF) sample outlining exclusions.

https://doi.org/10.1371/journal.pone.0027784.s001

(DOC)

Figure S2.

Derivation of Australian Cystic Fibrosis Data Registry (ACFDR) sample and Australian Cystic Fibrosis BAL Study (ACFBAL) sample outlining exclusions.

https://doi.org/10.1371/journal.pone.0027784.s002

(DOC)

Figure S3.

Surivial analysis for first positive respiratory culture for Pseudomonas aeruginosa in the CFTSS and ACFBAL samples by temperature quartile. In both study samples, the warmer temperatures were associated with earlier acquisition of P. aeruginosa.

https://doi.org/10.1371/journal.pone.0027784.s003

(DOC)

Table S1.

Complete Regression Analyses for Predictors of Lung Function (CF-specific FEV1).

https://doi.org/10.1371/journal.pone.0027784.s004

(DOC)

Table S2.

Complete Logistic Regression Analyses for Predictors of P. aeruginosa Infection.

https://doi.org/10.1371/journal.pone.0027784.s005

(DOC)

Table S3.

Complete Logistic Regression Analyses for Predictors of Mucoid P. aeruginosa Infection.

https://doi.org/10.1371/journal.pone.0027784.s006

(DOC)

Table S4.

Complete Logistic Regression Analyses for Predictors of B. cepacia Complex Infection.

https://doi.org/10.1371/journal.pone.0027784.s007

(DOC)

Table S5.

Study Outcomes by Temperature Quartile.

https://doi.org/10.1371/journal.pone.0027784.s008

(DOC)

Table S6.

Regression Analyses for Lung Function: Assessing the Mediation Effect of P. aeruginosa.

https://doi.org/10.1371/journal.pone.0027784.s009

(DOC)

Table S7.

Regression Analyses for P. aeruginosa : Assessing the Mediation Effect of Lung Function.

https://doi.org/10.1371/journal.pone.0027784.s010

(DOC)

Acknowledgments

The authors thank the Cystic Fibrosis Foundation Patient Registry, especially Ase Sewall and Monica Brooks; the Australian Cystic Fibrosis Data Registry, especially Geoff Sims, Terry Stewart, and its advisory committee; the Australian Cystic Fibrosis Bronchoalveolar Lavage Study, especially Joyce Cheney; the support of the Office of Medical Research, Queensland Health; Scott Blackman and Beryl Rosenstein for helpful discussions, and most importantly the patients with cystic fibrosis and their families, research coordinators, nurses, and physicians who are participating in the U.S. Cystic Fibrosis Twin and Sibling Study. Joseph Collaco and Garry Cutting assume responsibility for the contents of this manuscript.

Author Contributions

Conceived and designed the experiments: JC JM GC. Performed the experiments: JC DG KN C.Watson TS. Analyzed the data: JC GC. Contributed reagents/materials/analysis tools: C.Watson TS SB C.Wainwright GC. Wrote the paper: JC GC. Reviewed and edited the manuscript: JC JM DG KN C.Watson TS SB C.Wainwright GC.

References

  1. 1. Kerem E, Corey M, Kerem B-S, Rommens J, Markiewicz D, et al. (1990) The relation between genotype and phenotype in cystic fibrosis--analysis of the most common mutation (deltaF508). N Engl J Med 323: 1517–1522.
  2. 2. The Cystic Fibrosis Genotype-Phenotype Consortium (1993) Correlation between genotype and phenotype in patients with cystic fibrosis. N Engl J Med 329: 1308–1313.
  3. 3. Koch C, Cuppens H, Rainisio M, Madessani U, Harms H, et al. (2001) European Epidemiologic Registry of Cystic Fibrosis (ERCF): comparison of major disease manifestations between patients with different classes of mutations. Pediatr Pulmonol 31: 1–12.
  4. 4. Collaco JM, Blackman SM, McGready J, Naughton KM, Cutting GR (2010) Quantification of the relative contribution of environmental and genetic factors to variation in cystic fibrosis lung function. J Pediatr 157: 802–807.
  5. 5. Stanke F, Becker T, Kumar V, Hedtfeld S, Becker C, et al. (2011) Genes that determine immunology and inflammation modify the basic defect of impaired ion conductance in cystic fibrosis epithelia. J Med Genet 48: 24–31.
  6. 6. Collaco JM, Vanscoy L, Bremer L, McDougal K, Blackman SM, et al. (2008) Interactions between secondhand smoke and genes that affect cystic fibrosis lung disease. JAMA 299: 417–424.
  7. 7. Beydon N, Amsallem F, Bellet M, Boule M, Chaussain M, et al. (2002) Pulmonary function tests in preschool children with cystic fibrosis. Am J Respir Crit Care Med 166: 1099–1104.
  8. 8. Campbell PW III, Parker RA, Roberts BT, Krishnamani MR, Phillips JA, III (1992) Association of poor clinical status and heavy exposure to tobacco smoke in patients with cystic fibrosis who are homozygous for the F508 deletion. J Pediatr 120: 261–264.
  9. 9. Rubin BK (1990) Exposure of children with cystic fibrosis to environmental tobacco smoke. N Engl J Med 323: 782–788.
  10. 10. Smyth A, O'Hea U, Williams G, Smyth R, Heaf D (1994) Passive smoking and impaired lung function in cystic fibrosis. Arch Dis Child 71: 353–354.
  11. 11. Goss CH, Newsom SA, Schildcrout JS, Sheppard L, Kaufman JD (2004) Effect of ambient air pollution on pulmonary exacerbations and lung function in cystic fibrosis. Am J Respir Crit Care Med 169: 816–821.
  12. 12. O'Connor GT, Quinton HB, Kneeland T, Kahn R, Lever T, et al. (2003) Median household income and mortality rate in cystic fibrosis. Pediatrics 111: e333–e339.
  13. 13. Schechter MS, McColley SA, Silva S, Haselkorn T, Konstan MW, et al. (2009) Association of socioeconomic status with the use of chronic therapies and healthcare utilization in children with cystic fibrosis. J Pediatr 155: 634–639.
  14. 14. Stephenson A, Hux J, Tullis E, Austin PC, Corey M, et al. (2011) Socioeconomic status and risk of hospitalization among individuals with cystic fibrosis in Ontario, Canada. Pediatr Pulmonol 46: 376–384.
  15. 15. Schechter MS, Shelton BJ, Margolis PA, Fitzsimmons SC (2001) The association of socioeconomic status with outcomes in cystic fibrosis patients in the United States. Am J Respir Crit Care Med 163: 1331–1337.
  16. 16. Schechter MS, Margolis PA (1998) Relationship between socioeconomic status and disease severity in cystic fibrosis. J Pediatr 132: 260–264.
  17. 17. Curtis JR, Burke W, Kassner AW, Aitken ML (1997) Absence of health insurance is associated with decreased life expectancy in patients with cystic fibrosis. Am J Respir Crit Care Med 155: 1921–1924.
  18. 18. Vanscoy LL, Blackman SM, Collaco JM, Bowers A, Lai T, et al. (2007) Heritability of lung disease severity in cystic fibrosis. Am J Respir Crit Care Med 175: 1036–1043.
  19. 19. Kulich M, Rosenfeld M, Campbell J, Kronmal R, Gibson RL, et al. (2005) Disease-specific reference equations for lung function in patients with cystic fibrosis. Am J Respir Crit Care Med 172: 885–891.
  20. 20. Green DM, McDougal KE, Blackman SM, Sosnay PR, Henderson LB, et al. (2010) Mutations that permit residual CFTR function delay acquisition of multiple respiratory pathogens in CF patients. Respir Res 11: 140.
  21. 21. Baron RM, Kenny DA (1986) The moderator-mediator variable distinction in social psychological research: conceptual, strategic, and statistical considerations. J Pers Soc Psychol 51: 1173–1182.
  22. 22. Kenny DA (2011) Mediation. 2011): http://davidakenny.net/cm/mediate.htm. (Last updated 08/23/2011; last accessed 09/03/.
  23. 23. Lipuma JJ (2005) Update on the Burkholderia cepacia complex. Curr Opin Pulm Med 11: 528–533.
  24. 24. Lipuma JJ (2010) The changing microbial epidemiology in cystic fibrosis. Clin Microbiol Rev 23: 299–323.
  25. 25. McKone EF, Emerson SS, Edwards KL, Aitken ML (2003) Effect of genotype on phenotype and mortality in cystic fibrosis: a retrospective cohort study. Lancet 361: 1671–1676.
  26. 26. Kosorok MR, Jalaluddin M, Farrell PM, Shen G, Colby CE, et al. (1998) Comprehensive analysis of risk factors for acquisition of Pseudomonas aeruginosa in young children with cystic fibrosis. Pediatr Pulmonol 26: 81–88.
  27. 27. Rosenfeld M, Emerson J, McNamara S, Joubran K, Retsch-Bogart G, et al. (2010) Baseline characteristics and factors associated with nutritional and pulmonary status at enrollment in the cystic fibrosis EPIC observational cohort. Pediatr Pulmonol 45: 934–944.
  28. 28. Patz JA, Campbell-Lendrum D, Holloway T, Foley JA (2005) Impact of regional climate change on human health. Nature 438: 310–317.
  29. 29. Cho JC, Tiedje JM (2000) Biogeography and degree of endemicity of fluorescent Pseudomonas strains in soil. Appl Environ Microbiol 66: 5448–5456.
  30. 30. Musken M, Di FS, Dotsch A, Fischer R, Haussler S (2010) Genetic determinants of Pseudomonas aeruginosa biofilm establishment. Microbiology 156: 431–441.
  31. 31. Conibear TC, Collins SL, Webb JS (2009) Role of mutation in Pseudomonas aeruginosa biofilm development. PLoS One 4: e6289.
  32. 32. Fonseca AP, Correia P, Sousa JC, Tenreiro R (2007) Association patterns of Pseudomonas aeruginosa clinical isolates as revealed by virulence traits, antibiotic resistance, serotype and genotype. FEMS Immunol Med Microbiol 51: 505–516.
  33. 33. Cappello S, Guglielmino SP (2006) Effects of growth temperature on polystyrene adhesion of Pseudomonas aeruginosa ATCC 27853. Brazilian Journal of Microbiology 37: 205–207.
  34. 34. Li Z, Kosorok MR, Farrell PM, Laxova A, West SE, et al. (2005) Longitudinal development of mucoid Pseudomonas aeruginosa infection and lung disease progression in children with cystic fibrosis. JAMA 293: 581–588.
  35. 35. Hankinson JL, Odencrantz JR, Fedan KB (1999) Spirometric reference values from a sample of the general U.S. population. Am J Respir Crit Care Med 159: 179–187.
  36. 36. U.S.Environmental Protection Agency (2008) A review of the impacts of climate variability and change on aeroallergens and their associated effects. US EPA/600/R-06/ 164F:
  37. 37. Hallstrand TS, Calenoff E, Becker JW, Henderson WR Jr, Aitken ML (2004) The role of allergy in manifestations of respiratory disease in adult cystic fibrosis. Ann Allergy Asthma Immunol 92: 228–233.
  38. 38. Roberts S (2004) Interactions between particulate air pollution and temperature in air pollution mortality time series studies. Environ Res 96: 328–337.
  39. 39. Bloomer BJ, Stehr JW, Piety CA, Salawitch RJ, Dickerson RR (2009) Observed relationships of ozone air pollution with temperature and emissions. Geophysical Research Letters 36: L09803: