The daily dynamics of cystic fibrosis airway microbiota during clinical stability and at exacerbation
© Carmody et al.; licensee BioMed Central. 2015
Received: 21 October 2014
Accepted: 24 February 2015
Published: 1 April 2015
Recent work indicates that the airways of persons with cystic fibrosis (CF) typically harbor complex bacterial communities. However, the day-to-day stability of these communities is unknown. Further, airway community dynamics during the days corresponding to the onset of symptoms of respiratory exacerbation have not been studied.
Using 16S rRNA amplicon sequencing of 95 daily sputum specimens collected from four adults with CF, we observed varying degrees of day-to-day stability in airway bacterial community structures during periods of clinical stability. Differences were observed between study subjects with respect to the degree of community changes at the onset of exacerbation. Decreases in the relative abundance of dominant taxa were observed in three subjects at exacerbation. We observed no relationship between total bacterial load and clinical status and detected no viruses by multiplex PCR.
CF airway microbial communities are relatively stable during periods of clinical stability. Changes in microbial community structure are associated with some, but not all, pulmonary exacerbations, supporting previous observations suggesting that distinct types of exacerbations occur in CF. Decreased abundance of species that are dominant at baseline suggests a role for less abundant taxa in some exacerbations. Daily sampling revealed patterns of change in microbial community structures that may prove useful in the prediction and management of CF pulmonary exacerbations.
KeywordsCystic fibrosis Respiratory exacerbation Lung microbiome Airway microbiome
Chronic pulmonary infection and inflammation culminating in respiratory failure is the primary cause of death in people with cystic fibrosis (CF). Studies employing culture-independent analyses have shown that CF airways typically harbor complex microbial communities comprised of numerous bacterial species, particularly during the early and intermediate stages of lung disease [1-9]. While these studies have made important contributions toward understanding the dynamics of these communities, they have been limited by cross-sectional design and/or their reliance on analysis of a very small number (often only one or a few) of respiratory samples taken during periods of interest. Samples representing the periodic exacerbations of respiratory symptoms that characterize CF most often have been obtained after the initiation of antibiotics for exacerbation treatment. A more complete understanding of the dynamics of airway microbiota, in particular with respect to changes occurring at the time of pulmonary exacerbation, would benefit from an analysis of serial samples obtained in the days preceding and coincident with exacerbation onset, but prior to the initiation of antibiotic therapy for exacerbation. A meaningful assessment of airway microbial community changes during this crucial period can only be made, however, in the context of an understanding of the day-to-day variability of airway communities during periods of clinical stability. Rigorous analyses of the daily dynamics of the airway microbiota around these critical periods are lacking.
We obtained daily sputum samples from persons with CF in order to characterize changes in the bacterial and viral airway community during periods of clinical stability and during the onset of exacerbation of clinical symptoms, but prior to the initiation of antibiotics for exacerbation. We quantified community stability to assess the day-to-day variation in community structures and examine the degree of inter-individual variation in microbiome stability. We also describe the changes in community structure at the onset of symptoms of exacerbation in these four individuals.
Demographic and sample collection information for four study subjects
% FEV 1 a
Most abundant OTU at baseline
Number of samples
Day-to-day community variation when clinically stable
We next asked whether the degree of day-to-day change in baseline bacterial community structure varied between subjects and/or within subjects over time. The mean BC dissimilarity among day-to-day pairwise comparisons of all baseline samples varied between the four subjects, with subject A having the highest mean (0.255) and subject D having the lowest (0.144) (Figure 2). The mean BC dissimilarity during the second week of baseline samples from subject A was greater than that observed during the prior week (0.388 vs. 0.057, respectively). Further, a very strong negative correlation was observed between day-to-day BC dissimilarity and days to start of symptoms for subject A, indicating that the community became more variable as this subject approached exacerbation (P = −0.842). In contrast, strong correlations between day-to-day BC dissimilarity and days to the start of symptoms were not observed for the other three subjects (P < 0.6).
Homogeneity of molecular variance (HOMOVA), a nonparametric analog to Bartlett's test for homogeneity of variance, was used to determine if subjects differed significantly from each other with respect to overall variability in community structure at baseline. Significantly greater variation in community structure was observed in subject A, compared to either subject C or subject D (P = 0.002 and 0.003, respectively), but not in comparison to subject B (P = 0.452). In contrast, no statistically significant differences were observed in the degree of variation in community structure at baseline between subjects B, C, and D.
Community movement during exacerbation
Study subjects recorded symptoms on a daily questionnaire based on the Respiratory Symptoms Domain of the Cystic Fibrosis Questionnaire-Revised [10,11], and a symptom score was calculated by summing the total number of positive symptoms each day. We observed no strong correlations between day-to-day BC dissimilarity and either symptom score or change in symptom score from the previous day (P < 0.6). When grouped into baseline and exacerbation categories, the average day-to-day BC dissimilarity was higher for exacerbation samples compared to baseline samples for subject B (0.361 and 0.213, respectively), and lower for exacerbation samples compared to baseline samples for subject D (0.047 and 0.144, respectively). (It should be noted, however, that only three exacerbation samples were available for subject D prior to the initiation of antibiotics). No large differences in average day-to-day BC dissimilarity were observed when comparing day-to-day BC dissimilarity of exacerbation samples to baseline samples for subject A (0.246 and 0.255, respectively) and subject C (0.177 and 0.196, respectively).
After exploring global community changes, we looked for patterns of change in specific genera at exacerbation onset. A strong negative correlation was observed between collection day and the relative abundance of Staphylococcus for subject A (P = −0.659), while a moderate correlation was observed between collection day and the relative abundance of Streptococcus for subject C (P = −0.426). The average relative abundance of these operational taxonomic units (OTUs) was lower in exacerbation samples compared to baseline samples in both subjects (A 0.607 vs. 0.788, respectively; C 0.290 vs. 0.342, respectively). The relative abundance of the dominant OTU for subject B (Burkholderia) also decreased at exacerbation compared to baseline (average relative abundance 0.525 vs. 0.387, respectively), although we observed only a weak negative correlation between collection day and Burkholderia relative abundance (P = −0.253). During exacerbation, the community structure of subject B's samples was quite variable; large day-to-day shifts in relative abundances of the top OTUs were observed (Figure 3, Additional file 1: Table S1). No large changes in relative abundance of subject D's dominant OTU were observed at exacerbation, nor were strong correlations observed between absolute abundance of the dominant OTUs (calculated as the product of total bacterial load and the relative abundance of the OTU) and collection day for any subject.
Total bacterial load was not markedly different between baseline and exacerbation samples for any of the four subjects, nor were strong correlations observed between bacterial load and symptom score (P < 0.6). While total bacterial load for subjects B, C, and D fluctuated within a 2-log range over the entire period of observation, the total load for subject A was more variable, fluctuating over a greater than 3-log range (Additional file 3: Figure S2). None of the 31 samples assessed for the presence of respiratory viruses tested positive by the FilmArray Respiratory Panel.
Defining airway community stability
Bacterial community stability was quantified across the entire period of observation by defining, for each subject, a zone for changes in community structure, within which variations in community composition could be explained by stochastic events (‘stable’) as described in Methods (Figure 4). When communities move outside this zone, we interpret this to be due to a perturbation that changes the community structure (‘unstable’). Communities in subject A moved outside the stability zone 7 days prior to symptom onset, before becoming relatively stable again after 4 days of symptoms. Communities in subject B also moved outside the stability zone, but unlike subject A, this occurred after the onset of symptoms and continued until at least the beginning of antibiotic treatment. Communities in subject C remained within the stability zone for the entire period of observation, indicating relatively stable communities before and during exacerbation onset, while subject D moved outside the stability zone briefly in the week prior to the onset of symptoms.
Repeated exacerbations of pulmonary symptoms in persons with CF are associated with a progressive decline in lung function. Although clinical management based on antibiotic treatment and intensified airway clearance is usually at least partially effective [10-12], patients fail to return to 100% of baseline lung function in greater than 40% of exacerbations treated with intravenous antibiotics . Unfortunately, relatively little is known about the pathophysiologic basis of exacerbations. The identification of microbial signals that correlate with exacerbation would aid efforts to prevent or mitigate the effects of these events and would inform studies investigating the mechanisms by which these factors may drive exacerbation. A prerequisite to understanding how CF airway communities change at exacerbation, however, is a better appreciation of the day-to-day variability in these communities during periods of clinical stability. Unfortunately, no previous studies have investigated CF airway microbial community variability with this level of granularity. To this end, we analyzed the daily dynamics of the airway bacterial communities in persons with CF during periods of clinical stability and at the onset of pulmonary exacerbations, prior to the initiation of antibiotic therapy for exacerbation.
In the absence of clinical symptoms, bacterial airway communities remained relatively stable. Subject A had overall greater variability at baseline than subjects C or D, due to destabilization of subject A's community during the week immediately preceding the onset of exacerbation symptoms. The airway community in subject D showed more subtle instability during this time. Because the baseline samples from all four subjects were obtained 1 to 3 weeks immediately preceding the onset of symptoms, it is possible that community variability was already influenced by impending exacerbation even in the absence of clinical symptoms. Analyses of daily samples obtained during periods of clinical stability further removed from exacerbation are warranted to better ascertain the degree of day-to-day variability in bacterial airway communities in CF in general. For the purposes of our study, however, we believed it is important to assess changes in community structure during the period of clinical stability immediately preceding the onset of exacerbation symptoms.
The consistency in these subject's airway communities over the course of several days has bearing on questions pertaining to bacterial biogeography within diseased lungs and on how well sputum samples represent airway communities in microbiome studies in CF. Recent studies have demonstrated that microbial communities vary over the spatial landscape of the diseased lung [14-17]. The remarkable level of community stability observed in our study subjects, particularly subject D over the first 10 days of observation (Figures 2 and 4), would not be expected if serial sputum specimens were sampling distinct regions of diseased lungs from 1 day to the next. This consistency also suggests that measures of lung microbial community structure derived from expectorated sputum are not substantially impacted by ‘contamination’ from oropharyngeal microbiota. In a recent comparison of multiple body sites, the temporal variability in oral cavity microbiomes was found to be significantly greater than those in vaginal, gut, skin, and anterior nares sites . If sputum samples were significantly contaminated with oral microbiota during expectoration, we would expect to see variability reflecting that observed in serial oral cavity samples. Finer scale sampling of oral sites, analyzed alongside paired sputum samples from the same subject, would further illuminate this issue.
With the onset of symptoms of exacerbation, we observed distinct patterns of bacterial community change among the four subjects. The community in subject A showed considerable movement several days before the onset of symptoms, then shifted to an alternate stable state different in structure from the previous week. This change can be visualized by NMDS (Figure 5), wherein subject A's baseline samples occupy a different space in the ordination plot than the exacerbation samples, the result, in part, of a decrease in the relative abundance of Staphylococcus and increases in the relative abundances of Pseudomonas and Prevotella (Figure 3, Table S1).
The airway bacterial community in subject B destabilized coincident with the onset of exacerbation symptoms (Figure 4). However, in contrast to subject A, the community did not shift to a new structure (that is, the centroid of subject B's exacerbation samples in the NMDS plot remained near that of the baseline samples, Figure 5). However, approximately 1 week after the onset of symptoms, the community shifted intermittently from one dominated by Burkholderia to one dominated by Pseudomonas.
The communities in subjects C and D were relatively stable through the onset of symptoms (Figure 4). This result was unexpected for subject C, who had the greatest baseline community diversity among the subjects. Previous work suggested that a high degree of baseline community diversity predisposes to greater change at exacerbation . Subject D was prescribed antibiotic therapy very soon (3 days) after the onset of exacerbation symptoms. The three daily samples immediately prior to antibiotic administration showed exceptionally stable communities, with a maximum BC dissimilarity between communities of only 0.032, which is within the range of the replicate controls. Whether more movement of this subject's airway community would have been observed in a larger set of samples before the initiation of treatment is unknown.
Despite the inter-subject differences in community dynamics with exacerbation onset in these four subjects, common features are also apparent. In three of four subjects, the dominant taxon (Staphylococcus, Burkholderia, and Streptococcus in subjects A, B, and C, respectively) decreased in relative abundance around the time of exacerbation. Although this change did not reflect a decrease in the absolute density of these species, due to day-to-day variability in total bacterial density, this finding is consistent with previous work that demonstrated decreases in the relative abundance of the dominant taxa at the time of exacerbation in another set of subjects . The pathophysiologic implications of this intriguing finding require further elucidation, but suggest a role for less abundant taxa in some pulmonary exacerbations.
In all four subjects, despite subtle day-to-day variability in total bacterial load, we observed no overall changes in total bacterial density between baseline and exacerbation samples. This finding supports studies that similarly failed to show significant changes in overall bacterial density with exacerbation [3,7,9]. These findings suggest that shifts in the relative abundance of bacterial community members, rather than changes in total bacterial density, are more likely associated with alterations in clinical state.
Although respiratory viral infection is associated with a significant minority of pulmonary exacerbations in CF [19-21], we detected no respiratory viruses in samples obtained around the time of exacerbation onset in any of the four subjects. This analysis was limited to the 17 respiratory viruses included in the commercial viral detection system we employed. Nevertheless, this panel would be expected to detect the great majority of viruses responsible for viral respiratory infection in CF based on recent epidemiologic surveys [19,20].
Our study design allowed us to assess the correlation between subject-reported symptoms and airway community structure. In this regard, we observed no strong correlations between daily symptom scores and day-to-day community variability. However, it is possible that a larger scale study and/or a more refined symptom scoring system might detect such relationships, if they indeed exist. Weighting of symptoms, rather than merely summing the number of symptoms on each day, could reveal correlations with airway community movement that were not detected in this study. Similarly, consideration of other patient-specific features (for example, disease stage) and/or analysis of the temporal relationship between specific symptoms and microbiomic measures may prove worthwhile. Our study was limited to four subjects at a single center. Although samples from three of four subjects contained Pseudomonas, none were dominated by this taxon as is often seen in CF samples. The applicability of these findings to the larger CF population is thus unclear. Finally, we restricted our analyses to predominantly descriptive reports of trends within and between subjects, given the small sample size and correlated nature of the data. Our results indicate the potential for significance given increased power from an expanded pool of subjects.
Rapid advances in technical capability for microbiome analyses will provide opportunities to study larger sample sets to more completely characterize the spectrum of microbiome changes associated with pulmonary exacerbation. Transcriptomic and/or metabolomic analyses would complement these efforts in enhancing our understanding of how changes in the activity of airway microbiota relate to changes in the host's clinical state.
Daily airway bacterial community structures in CF are relatively stable during periods of clinical stability. The onset of symptoms of exacerbation may be heralded by marked shifts in these communities, even in the absence of viral infection or antibiotic therapy for the treatment of exacerbation. Whether such shifts only occur in a minority of exacerbations, as observed here, is unclear. Monitoring of airway microbial community structures therefore has the potential to identify signatures that could be useful in predicting CF respiratory exacerbations, which, in turn, may enable more prompt and/or appropriate therapy of exacerbation. A better understanding of the changes in airway communities around the time of exacerbations will also inform studies to elucidate the pathophysiologic mechanisms of exacerbation in CF.
Sputum specimens and metadata
With approval of the Institutional Review Board of the University of Michigan Medical School, sputum specimens were collected daily by four persons with CF and stored at 4°C for up to 27 days before being sent on ice to the investigator's laboratory for processing and storage at −80°C. Subjects were part of a larger cohort of CF patients participating in long-term daily collection of sputum for microbiome analysis. The first four exacerbations captured from this cohort were chosen for analysis in this study. Study subjects recorded symptoms and antibiotic use on a daily questionnaire based on the Respiratory Symptoms Domain of the Cystic Fibrosis Questionnaire-Revised [22,23]. Symptoms recorded included increased cough, change in sputum, shortness of breath, fever, decreased energy, decreased appetite, and missed work/activities. Evidence of hemoptysis in the specimen was recorded by the laboratory. A symptom score was calculated by summing the total number of positive symptoms, plus hemoptysis, each day. Baseline samples were defined as those obtained during a period of consecutive symptom-free days that culminated in the onset of symptoms of exacerbation. ‘Exacerbation’ samples were defined as those obtained from the start of consistent symptoms to the day prior to the start of antibiotic treatment for exacerbation . Six consecutive daily sputum samples obtained from subject A served as technical control samples. DNA was prepared from each of these six samples and each DNA samples was sequenced in duplicate on the same sequencing run. These six pairs of replicate control samples served to assess the degree of background variability between and within DNA sequencing runs (Figure 1).
DNA was prepared from frozen sputum as previously described . Briefly, samples were treated with Sputolysin (EMD Chemicals, Gibbstown, NJ, USA) and subjected to mechanical disruption by bead beating before DNA extraction by the MagNA Pure nucleic acid purification platform (Roche Diagnostics Corp., Indianapolis, IN, USA). Total bacterial load was measured by quantitative PCR using the universal primer/probe set of Nadkarni et al.  targeting the bacterial 16S rRNA gene as previously described .
Samples were selected for viral testing whenever a subject's clinical state changed (that is, from no symptoms to having symptoms or vice versa) or at least weekly when clinical state was unchanged for more than 7 days. The FilmArray Respiratory Panel (BioFire Diagnostics, Salt Lake City, UT, USA) was used to test sputum samples for the presence of 17 respiratory viruses [26,27]. Briefly, 1 ml of sputum was liquefied with SnotBuster SL Solution (Copan Diagnostics Inc., Murrieta, GA, USA) by mixing for 30 s at 2,000 rpm. After incubating at room temperature for 15 min, 300 μl of the solution was transferred to sample buffer and inoculated into the FilmArray pouch, which was loaded onto the FilmArray instrument according to the manufacturer's instructions.
DNA sequencing and community ecological metrics
Bar-coded pyrosequencing of a portion of the 16S rRNA V3-V5 hypervariable region was performed by the Human Genome Sequencing Center at Baylor College of Medicine using protocols developed for the Human Microbiome Project (http://www.hmpdacc.org/resources/tools_protocols.php) as previously described , with the following modifications. The software package mothur v1.29  was used to process sequences, which were assigned to OTUs using an average neighbor algorithm with a 0.03 dissimilarity cutoff as previously described . The total number of reads for each community was rarefied to 1,807, the smallest number of reads obtained in the sample set, to control for differences in sequencing depth before alpha (nonparametric Shannon index ) and community dissimilarity measures were calculated. BC dissimilarity was used to measure the difference in community structures between paired samples . Dominant OTU was defined as the most abundant OTU in the majority of baseline samples.
Measure of community stability
Community structure stability was quantified by defining, for each subject, a ‘stability zone’ within which variations in community composition can be explained by stochastic events (stable). Movement of the community outside of this zone is likely due to a perturbation that overrides these stochastic events (unstable) . For each time point, the similarity in community composition was calculated between that time point and the next time point as [1 − Bray-Curtis dissimilarity]. This similarity score, used as a proxy for stability, ranges from 0 to 1. A value of 0 indicates no similarity between communities (implying maximum instability), while a value of 1 indicates no differences between communities (implying perfect stability). We calculated the expected variability in this score due to stochastic processes by randomly permuting counts across the OTU-by-time point matrix for each patient 1,000 times. This generated 1,000 random matrices with neither sample nor temporal selection, while controlling for subject specificity, sampling depth, and global taxon abundance. For each permutation, the stability score [1 − Bray-Curtis dissimilarity] was calculated at all time points. Finally, 95% confidence intervals based on these 1,000 simulated values at each time point was used to create the stability zone.
Spearman's rank correlation coefficient P was calculated to study the associations between community measures, patient metadata, and time. A correlation coefficient P ≥ 0.8 was considered a very strong correlation, P = 0.60 to 0.79 a strong correlation, P = 0.4 to 0.59 a moderate correlation, P = 0.20 to 0.39 a weak correlation, and P ≤ 0.19 a negligible correlation. HOMOVA [32,33] analysis was used to test whether groups of baseline samples from each subject had equal variances in community structure. BC-based NMDS was used to visualize changes in community composition over time.
Availability of supporting data
Sequences have been deposited at NCBI SRA (SRA: SRP051730; BioProject: PRNJA271691).
nonmetric multidimensional scaling
operational taxonomic unit
homogeneity of molecular variance
The authors thank the dedication of the study subjects whose generous participation made this study possible. This work was supported by the U.S. National Institutes of Health NHLBI grant 1RC1HL100809-01 and CTSA grant UL1RR024986. JJL also was supported by the Cystic Fibrosis Foundation. AV was supported by a T32 grant from National Institutes of Health NHLBI (HL007749). Additional support was provided by the Charles Woodson Pediatric Research Fund and the Nesbitt Program for Cystic Fibrosis Research.
- Rogers GB, Carroll MP, Serisier DJ, Hockey PM, Jones G, Bruce KD. Characterization of bacterial community diversity in cystic fibrosis lung infections by use of 16 s ribosomal DNA terminal restriction fragment length polymorphism profiling. J Clin Microbiol. 2004;42:5176–83.View ArticlePubMed CentralPubMedGoogle Scholar
- Sibley CD, Parkins MD, Rabin HR, Duan K, Norgaard JC, Surette MG. A polymicrobial perspective of pulmonary infections exposes an enigmatic pathogen in cystic fibrosis patients. Proc Natl Acad Sci U S A. 2008;105:15070–5.View ArticlePubMed CentralPubMedGoogle Scholar
- Stressmann FA, Rogers GB, Marsh P, Lilley AK, Daniels TW, Carroll MP, et al. Does bacterial density in cystic fibrosis sputum increase prior to pulmonary exacerbation? J Cyst Fibros. 2011;10:357–65.View ArticlePubMedGoogle Scholar
- Tunney MM, Klem ER, Fodor AA, Gilpin DF, Moriarty TF, McGrath SJ, et al. Use of culture and molecular analysis to determine the effect of antibiotic treatment on microbial community diversity and abundance during exacerbation in patients with cystic fibrosis. Thorax. 2011;66:579–84.View ArticlePubMedGoogle Scholar
- Fodor AA, Klem ER, Gilpin DF, Elborn JS, Boucher RC, Tunney MM, et al. The adult cystic fibrosis airway microbiota is stable over time and infection type, and highly resilient to antibiotic treatment of exacerbations. PLoS One. 2012;7:e45001.View ArticlePubMed CentralPubMedGoogle Scholar
- Stressmann FA, Rogers GB, van der Gast CJ, Marsh P, Vermeer LS, Carroll MP, et al. Long-term cultivation-independent microbial diversity analysis demonstrates that bacterial communities infecting the adult cystic fibrosis lung show stability and resilience. Thorax. 2012;67:867–73.View ArticlePubMedGoogle Scholar
- Carmody LA, Zhao J, Schloss PD, Petrosino JF, Murray S, Young VB, et al. Changes in cystic fibrosis airway microbiota at pulmonary exacerbation. Ann Am Thorac Soc. 2013;10:179–87.View ArticlePubMed CentralPubMedGoogle Scholar
- Daniels TW, Rogers GB, Stressmann FA, van der Gast CJ, Bruce KD, Jones GR, et al. Impact of antibiotic treatment for pulmonary exacerbations on bacterial diversity in cystic fibrosis. J Cyst Fibros. 2013;12(1):22–8.View ArticlePubMedGoogle Scholar
- Price KE, Hampton TH, Gifford AH, Dolben EL, Hogan DA, Morrison HG, et al. Unique microbial communities persist in individual cystic fibrosis patients throughout a clinical exacerbation. Microbiome. 2013;1:27.View ArticlePubMed CentralPubMedGoogle Scholar
- Smith AL, Fiel SB, Mayer-Hamblett N, Ramsey B, Burns JL. Susceptibility testing of Pseudomonas aeruginosa isolates and clinical response to parenteral antibiotic administration: lack of association in cystic fibrosis. Chest. 2003;123:1495–502.View ArticlePubMedGoogle Scholar
- Cunningham S, McColm JR, Mallinson A, Boyd I, Marshall TG. Duration of effect of intravenous antibiotics on spirometry and sputum cytokines in children with cystic fibrosis. Pediatr Pulmonol. 2003;36:43–8.View ArticlePubMedGoogle Scholar
- Konstan MW, Morgan WJ, Butler SM, Pasta DJ, Craib ML, Silva SJ, et al. Risk factors for rate of decline in forced expiratory volume in one second in children and adolescents with cystic fibrosis. J Pediatr. 2007;151:134–9.View ArticlePubMedGoogle Scholar
- Sanders DB, Hoffman LR, Emerson J, Gibson RL, Rosenfeld M, Redding GJ, et al. Return of FEV1 after pulmonary exacerbation in children with cystic fibrosis. Pediatr Pulmonol. 2010;45:127–34.View ArticlePubMedGoogle Scholar
- Erb-Downward JR, Thompson DL, Han MK, Freeman CM, McCloskey L, Schmidt LA, et al. Analysis of the lung microbiome in the "healthy" smoker and in COPD. PLoS One. 2011;6:e16384.View ArticlePubMed CentralPubMedGoogle Scholar
- Willner D, Haynes MR, Furlan M, Schmieder R, Lim YW, Rainey PB, et al. Spatial distribution of microbial communities in the cystic fibrosis lung. ISME J. 2012;6:471–4.View ArticlePubMed CentralPubMedGoogle Scholar
- Willner D, Haynes MR, Furlan M, Hanson N, Kirby B, Lim YW, et al. Case studies of the spatial heterogeneity of DNA viruses in the cystic fibrosis lung. Am J Respir Cell Mol Biol. 2012;46:127–31.View ArticlePubMed CentralPubMedGoogle Scholar
- Brown PS, Pope CE, Marsh RL, Qin X, McNamara S, Gibson R, Burns JL, Deutsch G, Hoffman LR. Directly sampling the lung of a young child with cystic fibrosis reveals diverse microbiota. Ann Am Thorac Soc. 2014 Jul 29. [Epub ahead of print].
- Ding T, Schloss PD. Dynamics and associations of microbial community types across the human body. Nature. 2014;509:357–60.View ArticlePubMedGoogle Scholar
- Hoek RA, Paats MS, Pas SD, Bakker M, Hoogsteden HC, Boucher CA, et al. Incidence of viral respiratory pathogens causing exacerbations in adult cystic fibrosis patients. Scand J Infect Dis. 2013;45:65–9.View ArticlePubMedGoogle Scholar
- Flight WG, Bright-Thomas RJ, Tilston P, Mutton KJ, Guiver M, Morris J, et al. Incidence and clinical impact of respiratory viruses in adults with cystic fibrosis. Thorax. 2013;0:1–7.Google Scholar
- Etherington C, Naseer R, Conway SP, Whitaker P, Denton M, Peckham DG. The role of respiratory viruses in adult patients with cystic fibrosis receiving intravenous antibiotics for a pulmonary exacerbation. J Cyst Fibros. 2014;13:49–55.View ArticlePubMedGoogle Scholar
- Modi AC, Quittner AL. Validation of a disease-specific measure of health-related quality of life for children with cystic fibrosis. J Pediatr Psychol. 2003;28:535–46.View ArticlePubMedGoogle Scholar
- Goss CH, Quittner AL. Patient-reported outcomes in cystic fibrosis. Proc Am Thorac Soc. 2007;4:378–86.View ArticlePubMed CentralPubMedGoogle Scholar
- Zhao J, Schloss PD, Kalikin LM, Carmody LA, Foster BK, Petrosino JF, et al. Decade-long bacterial community dynamics in cystic fibrosis airways. Proc Natl Acad Sci U S A. 2012;109:5809–14.View ArticlePubMed CentralPubMedGoogle Scholar
- Nadkarni MA, Martin FE, Jacques NA, Hunter N. Determination of bacterial load by real-time pcr using a broad-range (universal) probe and primers set. Microbiology. 2002;148:257–66.PubMedGoogle Scholar
- Branche AR, Walsh EE, Formica MA, Falsey AR. Detection of respiratory viruses in sputum from adults by use of automated multiplex PCR. J Clin Microbiol. 2014;52:3590–6.View ArticlePubMedGoogle Scholar
- Ruggiero P, McMillen T, Tang YW, Babady NE. Evaluation of the BioFire FilmArray respiratory panel and the GenMark eSensor respiratory viral panel on lower respiratory tract specimens. J Clin Microbiol. 2014;52:288–90.View ArticlePubMed CentralPubMedGoogle Scholar
- Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75:7537–41.View ArticlePubMed CentralPubMedGoogle Scholar
- Chao A, Shen T-J. Nonparametric estimation of Shannon's index of diversity when there are unseen species in sample. Environ Ecol Stat. 2003;10:429–43.View ArticleGoogle Scholar
- Bray JR, Curtis JT. An ordination of the upland forest communities of southern Wisconsin. Ecol Monogr. 1957;27:326–49.View ArticleGoogle Scholar
- Sul WJ, Oliver TA, Ducklow HW, Amaral-Zettler LA, Sogin ML. Marine bacteria exhibit a bipolar distribution. Proc Natl Acad Sci U S A. 2013;110:2342–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Stewart CN, Excoffier L. Assessing population genetic structure and variability with RAPD data: application to Vaccinium macrocarpon (American Cranberry). J Evol Biol. 1996;9:153–71.View ArticleGoogle Scholar
- Schloss PD. Evaluating different approaches that test whether microbial communities have the same structure. ISME J. 2008;2:265–75.View ArticlePubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.