The Journal of Allergy and Clinical Immunology
Volume 127, Issue 2 , Pages 372-381.e3, February 2011

Airway microbiota and bronchial hyperresponsiveness in patients with suboptimally controlled asthma

  • Yvonne J. Huang, MD

      Affiliations

    • Division of Pulmonary and Critical Care Medicine, Department of Medicine, University of California, San Francisco, Calif
  • ,
  • Craig E. Nelson, PhD

      Affiliations

    • Marine Science Institute, University of California, Santa Barbara, Calif
  • ,
  • Eoin L. Brodie, PhD

      Affiliations

    • Earth Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, Calif
  • ,
  • Todd Z. DeSantis, MS

      Affiliations

    • Earth Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, Calif
  • ,
  • Marshall S. Baek, BS

      Affiliations

    • Department of Anesthesia and Perioperative Care, University of California, San Francisco, Calif
  • ,
  • Jane Liu, MS

      Affiliations

    • Division of Pulmonary and Critical Care Medicine, Department of Medicine, University of California, San Francisco, Calif
  • ,
  • Tanja Woyke, PhD

      Affiliations

    • US Department of Energy, Joint Genome Institute, Walnut Creek, Calif
  • ,
  • Martin Allgaier, PhD

      Affiliations

    • Department of Anesthesia and Perioperative Care, University of California, San Francisco, Calif
  • ,
  • Jim Bristow, MD

      Affiliations

    • US Department of Energy, Joint Genome Institute, Walnut Creek, Calif
  • ,
  • Jeanine P. Wiener-Kronish, MD

      Affiliations

    • Department of Anesthesia and Perioperative Care, University of California, San Francisco, Calif
  • ,
  • E. Rand Sutherland, MD, MPH

      Affiliations

    • Division of Pulmonary and Critical Care Medicine, Department of Medicine, National Jewish Health, Denver, Colo
  • ,
  • Tonya S. King, PhD

      Affiliations

    • Division of Biostatistics, Department of Public Health Sciences, Pennsylvania State University, Hershey, Pa
  • ,
  • Nikolina Icitovic, MAS

      Affiliations

    • Division of Biostatistics, Department of Public Health Sciences, Pennsylvania State University, Hershey, Pa
  • ,
  • Richard J. Martin, MD

      Affiliations

    • Division of Pulmonary and Critical Care Medicine, Department of Medicine, National Jewish Health, Denver, Colo
  • ,
  • William J. Calhoun, MD

      Affiliations

    • Division of Allergy, Pulmonary, Immunology, Critical Care and Sleep, Department of Internal Medicine, University of Texas Medical Branch at Galveston, Galveston, Tex
  • ,
  • Mario Castro, MD

      Affiliations

    • Division of Pulmonary and Critical Care Medicine, Department of Medicine, Washington University, St Louis, Mo
  • ,
  • Loren C. Denlinger, MD, PhD

      Affiliations

    • Division of Pulmonary and Critical Care Medicine, Department of Medicine, University of Wisconsin Schools of Medicine and Public Health, Madison, Wis
  • ,
  • Emily DiMango, MD

      Affiliations

    • Division of Pulmonary, Allergy, and Critical Care Medicine, Department of Medicine, Columbia University College of Physicians and Surgeons, New York, NY
  • ,
  • Monica Kraft, MD

      Affiliations

    • Division of Pulmonary, Allergy, and Critical Care Medicine, Department of Medicine, Duke University, Durham, NC
  • ,
  • Stephen P. Peters, MD, PhD

      Affiliations

    • Section on Pulmonary, Critical Care, Allergy and Immunological Diseases, Department of Internal Medicine, Wake Forest University Health Sciences, Winston-Salem, NC
  • ,
  • Stephen I. Wasserman, MD

      Affiliations

    • Allergy and Immunology Section, Department of Medicine, University of California San Diego, San Diego, Calif
  • ,
  • Michael E. Wechsler, MD

      Affiliations

    • Division of Pulmonary and Critical Care Medicine, Department of Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, Mass
  • ,
  • Homer A. Boushey, MD

      Affiliations

    • Division of Pulmonary and Critical Care Medicine, Department of Medicine, University of California, San Francisco, Calif
  • ,
  • Susan V. Lynch, PhD

      Affiliations

    • Division of Gastroenterology, Department of Medicine, University of California, San Francisco, Calif
    • Corresponding Author InformationReprint requests: Susan V. Lynch, PhD, Colitis and Crohn’s Disease Center, Division of Gastroenterology, Department of Medicine, Box 0538, University of California, San Francisco, CA 94143.
  • ,
  • National Heart, Lung, and Blood Institute’s Asthma Clinical Research Network

      Affiliations

    • Investigators of the Asthma Clinical Research Network are listed in Appendix E1 in this article’s Online Repository at www.jacionline.org.

Received 17 July 2010; received in revised form 8 October 2010; accepted 27 October 2010. published online 03 January 2011.

Article Outline

Background

Improvement in lung function after macrolide antibiotic therapy has been attributed to reduction in bronchial infection by specific bacteria. However, the airway might be populated by a more diverse microbiota, and clinical features of asthma might be associated with characteristics of the airway microbiota present.

Objective

We sought to determine whether relationships exist between the composition of the airway bacterial microbiota and clinical features of asthma using culture-independent tools capable of detecting the presence and relative abundance of most known bacteria.

Methods

In this pilot study bronchial epithelial brushings were collected from 65 adults with suboptimally controlled asthma participating in a multicenter study of the effects of clarithromycin on asthma control and 10 healthy control subjects. A combination of high-density 16S ribosomal RNA microarray and parallel clone library-sequencing analysis was used to profile the microbiota and examine relationships with clinical measurements.

Results

Compared with control subjects, 16S ribosomal RNA amplicon concentrations (a proxy for bacterial burden) and bacterial diversity were significantly higher among asthmatic patients. In multivariate analyses airway microbiota composition and diversity were significantly correlated with bronchial hyperresponsiveness. Specifically, the relative abundance of particular phylotypes, including members of the Comamonadaceae, Sphingomonadaceae, Oxalobacteraceae, and other bacterial families were highly correlated with the degree of bronchial hyperresponsiveness.

Conclusion: The composition of bronchial airway microbiota is associated with the degree of bronchial hyperresponsiveness among patients with suboptimally controlled asthma. These findings support the need for further functional studies to examine the potential contribution of members of the airway microbiota in asthma pathogenesis.

Key words: Microbiome, bacteria, asthma, 16S ribosomal RNA, PhyloChip

Abbreviations used: ICS, Inhaled corticosteroid, iNKT, Invariant natural killer T, MIA, Macrolides in Asthma, MRPP, Multiresponse Permutation Procedure and Indicator Species Analysis, NMDS, Nonmetric multidimensional scaling, Q-PCR, Quantitative PCR, 16S rRNA, 16S ribosomal RNA

 

Discuss this article on the JACI Journal Club blog: www.jaci-online.blogspot.com.

Interest in the potential role of bacterial colonization or infection of the bronchial mucosa in the pathogenesis of asthma has been raised by several recent reports: associations between bacterial colonization of neonatal airways and subsequent asthma development,1 evidence of bronchial infection by specific intracellular bacteria among adult asthmatic patients,2 and the efficacy of prolonged macrolide antibiotic treatment in subsets of asthmatic patients.3, 4, 5 Although the role of microbial infection in asthma pathogenesis remains unclear, studies to date have primarily applied culture-based or targeted molecular approaches to identify specific bacterial species of interest in the airways. However, recent developments in defining the human microbiome have demonstrated that the composition of bacterial communities colonizing mucosal surfaces, rather than simply the presence of individual species, can be important in defining states of health or disease.6, 7, 8, 9, 10, 11

Culture-independent microbiota profiling based on sequence polymorphisms in the 16S ribosomal RNA (16S rRNA) gene, which is present in all bacteria, has been widely applied in environmental ecologic studies12, 13, 14 and in studies of the human “superorganism.”9, 10, 11 16S rRNA–based phylogenetic analysis through sequencing or phylogenetic arrays permits detection of uncultured or fastidious bacteria, providing information on community composition without a priori knowledge of those present.15, 16, 17 One such platform, the 16S rRNA PhyloChip, has been applied in a number of environmental12, 14, 15, 18 and clinical6, 19, 20 studies. This high-density array contains approximately 500,000 probes that can differentiate approximately 8,500 bacterial taxa (defined as groups of organisms sharing ≥97% 16S rRNA sequence homology)15, 21 and demonstrates greater resolution of complex bacterial communities than traditional 16S rRNA clone library-sequencing approaches.6, 18, 21 It provides an ideal tool to examine relationships between microbiota composition, including the presence and relative abundance of members of the “rare biosphere,”22 and clinical characteristics of the disease.

Recognizing the frequency with which the tracheobronchial tree is exposed to the external environment and to secretions from the oropharynx or upper gastrointestinal tract, we hypothesized that complex bacterial communities might colonize asthmatic airway mucosa and exhibit relationships with clinical features of disease. Therefore in conjunction with a prospective study of the effects of extended clarithromycin therapy in adults with suboptimally controlled asthma,23 we conducted a pilot study of the airway microbiota in 65 adult asthmatic patients and 10 healthy control subjects using the PhyloChip and parallel clone library sequencing. Bronchial airway samples from asthmatic patients possessed greater bacterial burden and diversity than those from healthy subjects. Furthermore, the degree of bronchial hyperresponsiveness exhibited by subjects was related to community composition and the relative abundance of specific bacterial families comprising the airway microbiota. Portions of these results have previously been presented in abstract form.24, 25

Back to Article Outline

Methods 

Subjects 

Bronchial epithelial samples for microbial analysis were obtained from a subset of subjects enrolled in the Macrolides in Asthma (MIA) study23 (NCT00318708, clinicaltrials.gov) conducted by the National Heart, Lung, and Blood Institute–sponsored Asthma Clinical Research Network. Briefly, adults with clinically stable but suboptimally controlled asthma, defined as persistent symptoms on the Asthma Control Questionnaire26 after 4 weeks of standardized treatment with inhaled fluticasone, were studied. Relevant inclusion/exclusion criteria and study procedures are described in the Methods section in this article’s Online Repository at www.jacionline.org, including exclusion for symptoms of respiratory tract infection or asthma exacerbation within 6 weeks of the bronchoscopy visit. Bronchoscopy was performed before randomization to clarithromycin or placebo therapy. Three protected specimen bronchial brushings for microbiome analysis were obtained and sent to the University of California, San Francisco, where 10 healthy, nonsmoking, nonatopic adults without asthma were also enrolled. All protocols were approved by the institutional review board at each center.

Sample processing 

All samples first underwent screening by using PCR for the presence of the 16S rRNA gene with extracted total DNA (100 ng) and the universal 16S rRNA primers Bact-27F and Bact-1492R27 in a 40-cycle reaction. Samples were deemed positive or negative for bacteria based on the presence or absence of a visible 16S rRNA PCR product. Samples with any evidence of a 16S rRNA PCR product were subsequently subjected to 8 PCR reactions (per sample) by using the same primer set across a temperature gradient (48°C to 58°C) to maximize bacterial diversity captured. Amplicons from each sample were then pooled, purified, and gel quantified with E-gels (Invitrogen, Carlsbad, Calif) before hybridization of a standardized quantity to the array for each sample, as previously described.6, 15

Quantitative PCR analysis 

Quantitative PCR (Q-PCR) was performed for 21 samples to determine the total 16S rRNA copy number normalized to β-actin copy numbers. Further details are provided in the Methods section of this article’s Online Repository.

Data analysis 

Relationships between microbiota structure and study variables were analyzed by using multivariate statistics in combination with the Multi-response Permutation Procedure and Indicator Species Analysis (MRPP) and nonmetric multidimensional scaling (NMDS) ordination.28 Data were log transformed before analysis, where appropriate. Sorensen distance-based dissimilarity matrices were calculated and used for subsequent NMDS analysis to examine the relatedness of samples based on their microbiota composition. Three samples from asthmatic patients with mean relative Sorensen distances of greater than 2 deviations from the pooled mean were excluded for this analysis. Pearson correlations were performed to determine relationships between the relative abundance of all detected phylotypes and study variables with q value corrections for false discovery.29 Significance was conservatively based on a false discovery rate of less than 3% and a probability of no more than 2 taxa being falsely identified as significant for each test. Community diversity was estimated by using Shannon indices30 based on the number of taxa present (community “richness”) and their relative abundance (community “evenness”). Statistical analyses were conducted with PC-ORD version 5,31 JMP version 7 (SAS Institute, Inc, Chicago, Ill), and R software version 2.7.1 (http://www.r-project.org).

16S rRNA clone library sequencing 

Clone libraries for a subset of patients were constructed (by using the same 16S rRNA amplicon pool hybridized to the array) and sequenced at the Department of Energy’s Joint Genome Institute (299-374 clones per library). Quality-filtered sequences were aligned and chimera checked with NAST32 and Bellerophon,33 respectively, and then classified by means of G2-chip taxonomy (www.greengenes.lbl.gov).34

Phylogenetic tree construction 

Representative 16S rRNA sequences for taxa of interest were exported from the Greengenes database. A neighbor-joining tree with nearest-neighbor interchange was produced with FastTree.35 Trees were annotated by using the Interactive Tree of Life (http://www.itol.embl.de/).36

Back to Article Outline

Results 

Airway bacterial burden 

Bronchial brushings were obtained from 75 subjects: 65 asthmatic patients and 10 healthy control subjects. Subjects’ characteristics are summarized in Table I. All asthmatic patients fulfilled entry criteria for the parent MIA trial,23 and samples for this pilot study were processed as they were received with no additional selection criteria applied. On initial screening, 54 (83%) of 65 asthmatic patients and 8 (80%) of 10 control subjects exhibited a visible 16S rRNA PCR product. However, not all PCR-positive samples subsequently produced sufficient 16S rRNA amplicon for PhyloChip analysis. Although very small amplicon concentrations can be hybridized to the array, in our experience application of 250 ng provides good characterization of the respiratory microbiota.8, 19 A greater proportion of asthmatic patients (37/65 [56.9%]) than healthy subjects (3/10 [30%]) produced 250 ng of amplicon, although this difference was not significant (P = .17). To increase the number of samples for comparative analyses, we analyzed the community composition detected by using the PhyloChip in 6 paired samples with both 100 and 250 ng of amplicon. Finding no significant difference in the community composition detected in paired samples (PMRPP > 0.1), 7 additional subjects (5 asthmatic patients and 2 healthy subjects) were assayed by using PhyloChip with 100 ng of amplicon. Therefore array data from a total of 47 subjects (42 asthmatic patients and 5 healthy control subjects) were obtained for the subsequent analyses (see Fig E1 in this article’s Online Repository at www.jacionline.org).

Table I. Baseline characteristics of subjects
Asthmatic patients (n = 65)Healthy control subjects (n = 10)P value
Age (y)39.3 ± 11.235.5 ± 12.9NS
Sex (% female/male)52/4860/40NS
FEV1 (% predicted)74.9 ± 15.694.1 ± 12.2<.001
PC20 (mg · mL−1)1.3 (2.9)>16
Sputum eosinophils (%)2.0 ± 5.2NA
Sputum neutrophils (%)33.6 ± 23.7NA
Asthma Control Questionnaire score1.8 ± 0.7NA
Total IgE (IU · mL−1)128.4 (0.03)NA
Asthmatic group
Analyzed by PhyloChip (n = 42)Unable to analyze by
PhyloChip (n = 23)
P value
Age (y)40.4 ± 10.737.3 ± 12.1.3
Sex (% female/male)45/5552/48.4
FEV1 (% predicted)75.5 ± 13.677 ± 14.8.7
PC20 (mg · mL−1)1.2 (3.2)1.8 (2.2).2
Sputum eosinophils (%)2.2 ± 5.91.5 ± 3.3.5
Sputum neutrophils (%)30.0 ± 24.040.3 ± 22.3.1
Asthma Control Questionnaire score1.8 ± 0.81.8 ± 0.5.6
Positive allergen skin test result (reaction to ≥1 allergen)79% (33/42)74% (17/23)
Total IgE (IU · mL−1)114.6 (0.03)158.3 (0.03).4
History of last oral corticosteroid use (no. of subjects)
None since age 12 y155
Age 12 y to up to 2 years before study enrollment137.40
Within last 2 years before study enrollment1411
Exacerbation events during study (no. of subjects)
Within 1 wk of bronchoscopy52
>1 wk after bronchoscopy32.90
No exacerbations during study3419

Values are presented as means ± SDs, except where noted.

NA, Not applicable; NS, not significant.

Geometric mean (coefficient of variation) reported.

Samples with a negative 16S rRNA PCR screen result or less than 100 ng of total 16S rRNA PCR product for hybridization were not analyzed by using PhyloChip.

Fisher exact P value.

To confirm that 16S rRNA amplicon concentrations (in nanograms per microliter) reported by means of gel quantification accurately reflected bacterial burden, we performed total 16S rRNA Q-PCR and regression analysis on these independent measurements, which demonstrated strong concordance (r = 0.69, P < .001). These data also confirmed that samples with sufficient amplicon for array analysis had significantly greater 16S rRNA copy numbers compared with those unable to be analyzed by using PhyloChip (P = .002). Therefore gel-quantified 16S rRNA amplicon concentrations were used as a proxy for bacterial burden to examine potential relationships between airway bacterial burden and clinical variables. Of the samples analyzed by using PhyloChip, bacterial burden was significantly higher in asthmatic samples compared with that seen in control samples (P = .008, Fig 1). Bacterial burden appeared to decrease with increasing methacholine PC20 measurements, although this did not reach statistical significance (r = −0.21, P = .09), suggesting that community composition rather than actual bacterial burden might be more related to this clinical feature.

  • View full-size image.
  • Fig 1. 

    16S rRNA bacterial burden in bronchial brushings from asthmatic patients and healthy subjects. Of subjects analyzed by using PhyloChip (total n = 47), asthmatic patients exhibited greater bacterial burden (measured concentration of 16S rRNA PCR product) than healthy subjects (P = .008, Welch t test with log-transformed data). Median and interquartile ranges are noted.

Airway microbiota structure and clinical features of asthma 

Phylogenetic analysis identified 1941 bacterial taxa (detected in at least 1 subject) representing 161 bacterial families (see Table E1 in this article’s Online Repository at www.jacionline.org). As a form of validation, we compared the phylogenetic distribution of airway microbiota detected by using the array with reference trees for airway bacteria reported in the Human Microbiome Project (http://www.hmpdacc.org/reference_genomes.php; accessed May 13, 2010)37 and observed good agreement in the types of bacteria reported by both (data not shown). Community richness across samples from asthmatic patients was variable and ranged from 48 to 1240 taxa. Of the 5 healthy subjects, 3 exhibited low community richness (200-253 taxa), whereas 2 demonstrated greater richness (844 and 1121 taxa), suggesting that variability in airway microbiota composition was not exclusive to asthmatic patients.

We found no association between study center and variation in microbiota composition, indicating that sample collection by different operators in distinct geographic regions was not a major influence on community composition in this study. Multivariate analyses also showed no association between microbiota composition and spirometric measurements (FEV1 and forced vital capacity), sputum eosinophil or neutrophil percentages, history of exacerbations, or systemic corticosteroid treatment within the previous 2 years. However, microbiota structure differed significantly based on asthmatic or healthy group assignment (PMRPP = .028), and NMDS analysis demonstrated that variability in airway microbiota composition was correlated with both the degree of airway hyperresponsiveness (methacholine PC20 concentrations; r = 0.56, P < .001; Fig 2, A) and bacterial burden (r = −0.36, P = .015; Fig 2, B). This suggests that the types and relative abundance of bacteria present in the airways are associated with these 2 variables.

  • View full-size image.
  • Fig 2. 

    Relationships between variability in airway microbiota composition and community diversity with bronchial hyperresponsiveness. Each circle represents the microbiota present in a single sample (n = 44). A and B, NMDS ordination, based on Sorensen dissimilarity matrices, demonstrates that variability in community composition is most strongly related to bronchial hyperresponsiveness and bacterial burden. The NMDS axes have no inherent units and indicate spatial relationships of samples based on their phylogenetic (dis)similarity. Three samples from asthmatic patients were excluded due to extreme mean relative Sorensen distances from the pooled mean. C, Increased diversity (higher Shannon indices) is correlated with greater bronchial hyperresponsiveness (Pearson correlations with log-transformed data).

Airway bacterial diversity, community composition, and bronchial hyperresponsiveness 

Asthmatic patients possessed greater airway microbiota diversity than control subjects (P = .012). Methacholine PC20 concentrations were inversely correlated with diversity (r = −0.46, P < .003; Fig 2, C), suggesting a strong relationship between increased airway bacterial diversity and greater bronchial hyperresponsiveness.

Relationships between PC20 and relative abundance of all taxa detected were analyzed to determine which of the airway microbiota were most associated with bronchial hyperresponsiveness. After false discovery correction, approximately 100 taxa demonstrated a significant (P < .01, q ≤ 0.015) linear relationship between increasing relative abundance and greater bronchial hyperresponsiveness (Fig 3 and see Table E2 in this article’s Online Repository at www.jacionline.org). These taxa, primarily belonging to the Proteobacteria, represented 31 bacterial families, such as the Comamonadaceae, Sphingomonadaceae, Nitrosomonadaceae, Oxalobacteraceae, and Pseudomonadaceae, and included a number of potentially interesting organisms. Although samples were collected in a manner to minimize oral secretion contamination (eg, use of triple-lumen protected specimen brushes and sample collection early in the procedure preferentially from upper lobes), we recognized the potential for organisms associated with the oropharynx to reside in the lower airways. Thus we investigated whether these approximately 100 phylotypes have previously been identified in the oral cavity by interrogating both the extensive human oral microbiome sequence database38 and PubMed using species or genus names and the search terms “oral,” “oropharynx,” “mouth,” or “tongue.” The majority (87.5%) of these phylotypes were not found referenced to the oral cavity in either database.

  • View full-size image.
  • Fig 3. 

    Phylogenetic tree based on 16S rRNA representative gene sequences (GenBank accession numbers) of the approximately 100 bacterial taxa highly correlated with greater bronchial hyperresponsiveness (P < .01, q < 0.015). Colors represent different bacterial families. Asterisks denote taxa with member species previously associated with clinical disease or possessing notable functional features (see also Table E2 in this article’s Online Repository).

Posttreatment bronchoscopy was not performed in the parent study, and therefore examination of relationships between clarithromycin-induced changes in airway microbiota composition and treatment outcomes was not possible. Instead, we explored whether relationships existed between pretreatment airway microbiota and posttreatment clinical outcomes. Airway bacterial diversity was significantly greater in asthmatic subjects who demonstrated a significant reduction in bronchial reactivity in response to clarithromycin treatment (defined as at least a doubling in methacholine PC20 dose before and after treatment) than those with a smaller reduction (P = .03, Fig 4). This pattern was not evident in patients who received placebo intervention (P = .64).

  • View full-size image.
  • Fig 4. 

    Baseline airway microbiota diversity and change in bronchial hyperresponsiveness with clarithromycin treatment. Sixteen asthmatic patients with pretreatment samples analyzed by means of microarray received clarithromycin therapy. Subjects with significant improvement in bronchial hyperresponsiveness after clarithromycin (n = 9) possessed higher pretreatment airway bacterial diversity than those with no significant response (n = 3). Four subjects with incomplete data because of study discontinuation were not included.

16S rRNA clone library-sequencing validation 

The PhyloChip, like other array technologies, is potentially susceptible to cross-hybridization, despite features of the array design intended to minimize this. 16S rRNA clone libraries were constructed with the same amplicon pool hybridized to the array for a subset of asthmatic patients (n = 6) to validate array-based findings. An average of 333 clones per sample were analyzed, yielding approximately 2000 nonchimeric sequences representing 197 taxa, the majority of which were present in more than 1 sample. The majority of bacterial subfamilies identified by means of sequencing were also present in the corresponding array dataset, validating large numbers of array-based positive hits (Table II). Two subfamilies (the Fusobacteriaceae subfamily 3 and Prevotellaceae subfamily 1) were detected by means of sequencing only in 2 samples (1 subfamily per sample). However, taxa within these subfamilies were very close (positive-fraction value = 0.84 and 0.75, respectively) to the conservative arbitrary cutoff (positive-fraction value ≥ 0.9) used to define array-based presence or absence.15 Clone library analysis also indicated that 2 samples (AB01 and AB04) were each dominated by a specific bacterial family. The majority of sequences for sample AB01 were classified as Pasteurellaceae subfamily 1 (this sample also exhibited the lowest array-reported community richness [48 taxa]), whereas for subject AB04, a majority of the sequences belonged to the Cellulomonadaceae.

Table II. Bacterial families identified by using both 16S rRNA clone library sequencing and 16S rRNA PhyloChip
AB01 (374)AB02 (301)AB03 (299)
PhylumFamilySubfamilyPhylumFamilySubfamilyPhylumFamilySubfamily
FirmicutesStreptococcaceaesf-1ActinobacteriaCellulomonadaceaesf-1ActinobacteriaCellulomonadaceaesf-1
ProteobacteriaPasteurellaceaesf-1BacteroidetesFlavobacteriaceaesf-1BacteroidetesPrevotellaceaesf-1
BacteroidetesPrevotellaceaesf-1BacteroidetesPorphyromonadaceaesf-1FirmicutesEnterococcaceaesf-1
BacteroidetesPrevotellaceaesf-1FirmicutesPeptococcaceae/Acidaminococcaceaesf-11
FirmicutesLachnospiraceaesf-5FirmicutesPeptostreptococcaceaesf-5
FirmicutesMycoplasmataceaesf-1FirmicutesStaphylococcaceaesf-1
FirmicutesPeptococcaceae/Acidaminococcaceaesf-11FirmicutesStreptococcaceaesf-1
FirmicutesPeptostreptococcaceaesf-5FusobacteriaFusobacteriaceaesf-1
FirmicutesStaphylococcaceaesf-1FusobacteriaFusobacteriaceaesf-3
FirmicutesStreptococcaceaesf-1ProteobacteriaCampylobacteraceaesf-3
FusobacteriaFusobacteriaceaesf-1ProteobacteriaPasteurellaceaesf-1
FusobacteriaFusobacteriaceaesf-3TM7Unclassifiedsf-1
ProteobacteriaCampylobacteraceaesf-3
ProteobacteriaPasteurellaceaesf-1
SpirochaetesSpirochaetaceaesf-1
TM7Unclassifiedsf-1
AB04 (344)AB05 (332)AB32 (349)
PhylumFamilySubfamilyPhylumFamilySubfamilyPhylumFamilySubfamily
ActinobacteriaCellulomonadaceaesf-1ActinobacteriaMicrococcaceaesf-1ActinobacteriaCorynebacteriaceaesf-1
BacteroidetesPorphyromonadaceaesf-1BacteroidetesFlavobacteriaceaesf-1BacteroidetesPorphyromonadaceaesf-1
BacteroidetesPrevotellaceaesf-1BacteroidetesPorphyromonadaceaesf-1BacteroidetesPrevotellaceaesf-1
FirmicutesAerococcaceaesf-1BacteroidetesPrevotellaceaesf-1FirmicutesAerococcaceaesf-1
FirmicutesLachnospiraceaesf-5FirmicutesAerococcaceaesf-1FirmicutesEnterococcaceaesf-1
FirmicutesPeptococcaceae/Acidaminococcaceaesf-11FirmicutesEnterococcaceaesf-1FirmicutesLachnospiraceaesf-5
FirmicutesPeptostreptococcaceaesf-5FirmicutesErysipelotrichaceaesf-3FirmicutesMycoplasmataceaesf-1
FirmicutesStreptococcaceaesf-1FirmicutesLachnospiraceaesf-5FirmicutesPeptococcaceae/Acidaminococcaceaesf-11
FusobacteriaFusobacteriaceaesf-3FirmicutesPeptococcaceae/Acidaminococcaceaesf-11FirmicutesPeptostreptococcaceaesf-5
ProteobacteriaCampylobacteraceaesf-3FirmicutesPeptostreptococcaceaesf-5FirmicutesStaphylococcaceaesf-1
ProteobacteriaPasteurellaceaesf-1FirmicutesStaphylococcaceaesf-1FirmicutesStreptococcaceaesf-1
ProteobacteriaPseudomonadaceaesf-1FirmicutesStreptococcaceaesf-1FusobacteriaFusobacteriaceaesf-1
FusobacteriaFusobacteriaceaesf-1FusobacteriaFusobacteriaceaesf-3
ProteobacteriaCampylobacteraceaesf-3ProteobacteriaCampylobacteraceaesf-3
ProteobacteriaPasteurellaceaesf-1ProteobacteriaMoraxellaceaesf-3
ProteobacteriaPseudomonadaceaesf-1ProteobacteriaPasteurellaceaesf-1
SpirochaetesSpirochaetaceaesf-1ProteobacteriaPseudomonadaceaesf-1
TM7Unclassifiedsf-1TM7Unclassifiedsf-1

Phylotypes reported at the subfamily level classified by using Greengenes G2 chip taxonomy (number of clones sequenced per library indicated in parentheses). Subfamilies in boldface were undetected by means of array based on conservative criteria for determining taxon presence.

Back to Article Outline

Discussion 

To our knowledge, this is the first study to examine whether specific aspects of the airway microbiota are related to relevant clinical or physiologic features of asthma. Our findings suggest that bacterial diversity, variations in community composition, and the relative abundance of specific phylotypes are associated with the degree of bronchial hyperresponsiveness in asthmatic patients administered inhaled corticosteroids (ICSs). This culture-independent study also expands the repertoire of organisms of potential relevance to asthma pathogenesis within this cohort beyond those previously implicated.1, 2 These data, together with recent studies in other asthmatic patients,7 suggest that the microbiome of the airways, as with other discrete host niches, might be an important contributor to asthma pathophysiology and to the heterogeneity of disease.39, 40, 41

Overall, these findings highlight 2 key aspects in defining airway health. The existence of a bronchial microbiota suggests that the bronchial tree is unlikely to be completely sterile and that the composition of its microbiome might directly or indirectly modulate airway function. The role of airway microbiota in airway disease is underinvestigated compared with other human niches, such as the intestinal microbiota, which has a demonstrated role in chronic inflammatory gastrointestinal diseases.10, 42 That a bronchial epithelium–associated microbiota exists is not entirely unexpected, particularly in the setting of airway disease, ICS therapy, or both. Specific bacterial communities persist in the airways despite antimicrobial treatment in patients with chronic obstructive pulmonary disease and cystic fibrosis.8, 19 Recently, Hilty et al7 also described airway microbiota in patients with mild-to-moderate asthma (all receiving ICS therapy), patients with chronic obstructive pulmonary disease, and healthy subjects. Those with airway disease had an increased prevalence of pathogenic Proteobacteria compared with control subjects, although relationships to clinical measurements were not examined in their study.

Complex microbial communities are now recognized to reside in various host mucosal sites43 where perturbations of community structure are associated with disease.6, 8, 10, 42 Moreover, the airway epithelium is increasingly recognized as important in immunologic responses to environmental and microbial exposures.44 Thus the airway microbiome could potentially influence presentations of asthma. As such, our findings of specific community relationships with bronchial hyperresponsiveness in the context of a well-characterized asthmatic cohort receiving standardized baseline treatment implicate additional bacterial groups in asthma. These include a Nitrosomonas species possessing a functional nitric oxide reductase.45 Their relative abundance might be a microbial indicator of airway concentrations of nitric oxide that, although not assessed in this study, reflect airway inflammation46 and correlate with measures of airway hyperresponsiveness.47 Also notable is Oxalobacter formigenes (Oxalobacteraceae), an anaerobic bacterium known to be susceptible to macrolide antibiotics.48 Other bacteria of interest include the Comamonadaceae and Sphingomonadaceae. Members of the Comamonadaceae, previously identified in patients with cystic fibrosis,49 possess steroid-responsive degradation pathways,50, 51 and their presence could plausibly be related to the selective pressure of ICS therapy. This raises the intriguing possibility that steroid nonresponsiveness, as observed in some asthmatic patients despite adherence to therapy,52 might be due in part to the presence of airway bacteria with steroid-degrading capacity.

Given the existence of airway microbiota, it is conceivable that the coincident presence of multiple potentially pathogenic bacteria contributes to persistent airway inflammation in asthma. This possibility is supported by the finding of neutrophilia in bronchial airway samples of some patients with severe asthma despite regular use of high-dose ICSs.53 Such inflammation could represent an appropriate response by the host to inappropriate airway microbial colonization. For example, Sphingomonadaceae, members of which were significantly correlated with bronchial hyperresponsiveness, are characterized by the presence of cell-membrane glycosphingolipids, which are recognized by and can activate invariant natural killer T (iNKT) cells, resulting in induction of IL-4 and IL-13.54, 55 Glycolipid activation of iNKT cells can also induce airway hyperreactivity independent of conventional CD4+ T cells.56 Although the presence of iNKT cells in asthmatic airways remains controversial,57, 58 our findings suggest a possible role for Sphingomonadaceae in airway pathophysiology and the possibility that the observed variability of iNKT cell populations in asthmatic patients might depend in part on the relative abundance of these species.

We recognize that this pilot study is limited by the relatively small number of subjects, especially healthy control subjects, and by the absence of asthmatic patients not taking an ICS. Nonetheless, the data indicate that airway bacterial burden in subjects without airway disease is much lower than in asthmatic patients requiring ICS therapy. Because all asthmatic patients were required to be taking standardized ICSs for the parent study,23 we are unable to infer whether the higher bacterial burden among asthmatic patients is a function of the disease itself or ICS treatment, although the relationships of bacterial burden and the composition and diversity of the microbiota to bronchial hyperresponsiveness are robust. Determining the effects of ICS use on the airway microbiome will require further investigation and is an important question given the wide use of ICS therapies in airway diseases, which, among patients with chronic obstructive pulmonary disease, has been associated with an increased risk of pneumonia.59

The array-based analysis in this study identified a greater diversity of airway microbiota than has previously been described, although to our knowledge, there has been only one other study involving asthmatic patients that applied 16S rRNA phylogenetic analysis.7 In addition to clone library-sequencing validations, our results also agreed with those of Hilty et al7 in that we detected all bacterial phyla and genera identified in their study of the lower airways of 11 adult and 13 pediatric asthmatic patients on ICS therapy. The PhyloChip platform permits in-depth analysis of relatively large numbers of samples by using a standardized assay that applies several levels of stringent criteria for determining the microbiota profile in a given sample. It is impractical to perform sequence corroboration for every array-based positive hit. Moreover, the extent of microbial diversity detected is dependent on the community structure and the depth of sequencing performed.60 Our study included the largest number of clones sequenced to date for asthmatic airway samples and confirmed the presence of many array-detected phylotypes, indicating that complex microbiota do exist in this niche.

Finally, we observed that asthmatic patients exhibiting a significant decrease in bronchial hyperresponsiveness after clarithromycin treatment in the parent study possessed greater pretreatment airway bacterial diversity than nonresponders. Despite the small number of subjects in this exploratory analysis, the trend is interesting given that prior studies have also observed decreased bronchial responsiveness after prolonged macrolide therapy.4, 61 Furthermore, although the antimicrobial susceptibility of the entire airway microbiota is unknown, conceivably, a multitude of community members might be sensitive to macrolides, resulting in a reduction in bacterial diversity, burden, or both on treatment and manifesting clinically in a reduction in bronchial responsiveness. Thus although the efficacy of macrolides in airway disease is often attributed to their anti-inflammatory properties, their effectiveness could also reflect an extended spectrum of antibacterial activity against members of the airway microbiota, the composition of which, as in the setting of ICSs, might influence outcomes.

In summary, several features of the airway microbial community are significantly associated with the degree of bronchial hyperresponsiveness among patients with suboptimally controlled asthma. Despite limitations, these study findings are notable in providing the first evidence for the potential functional, physiologic, and clinical relevance of the airway microbiome in asthma. Our results suggest that variations in airway microbiome structure and function might exert distinct effects that contribute to asthma heterogeneity and provide novel targets and hypotheses for future studies on disease mechanisms.

Clinical implications

Relationships between airway microbiota and clinical features of asthma were investigated by using culture-independent approaches. Among adults with suboptimally controlled asthma, bronchial hyperresponsiveness correlated strongly with airway microbiota composition and diversity.

Back to Article Outline

 

We thank the Steering Committee of the Asthma Clinical Research Network, which permitted and indeed encouraged our use of clinical data from the MIA trial for our analysis of relationships with bronchial microbial community composition, and we also thank the Asthma Clinical Research Network investigators (see the Acknowledgments section) and coordinators who obtained and provided the protected bronchial brush samples from asthmatic patients. We also thank Yvette Piceno and Gary Andersen at the Lawrence Berkeley National Laboratory for their PhyloChip expertise and advice.

The Asthma Clinical Research Network investigators were as follows: Bill T. Ameredes, Eugene Bleecker, William J. Calhoun, Mario Castro, Reuben Cherniack, Vernon M. Chinchilli, Timothy J. Craig, Loren Denlinger, Emily DiMango, John Fahy, Elliot Israel, Nizar Jarjour, Monica Kraft, Stephen C. Lazarus, Robert F. Lemanske, Stephen Peters, Joe Ramsdell, Robert A. Smith, Christine A. Sorkness, Stanley J. Szefler, Michael J. Walter, Stephen Wasserman, and Michael E. Wechsler.

Back to Article Outline

Methods 

Subject characterization and study procedures 

For the MIA clinical trial, adult men and women 18 to 60 years of age meeting the diagnostic criteria for asthma were enrolled through 10 Asthma Clinical Research Network clinical research centers. Subjects had to have suboptimal asthma control, which was defined as persistent symptoms at least as severe as those of mild persistent asthma, as reflected by an Asthma Control Questionnaire score of greater than 1.25, after 4 weeks of standardized treatment with 88 μg of fluticasone proprionate administered twice daily through a metered-dose inhaler (GlaxoSmithKline, London, United Kingdom). Exclusion criteria included a smoking history of greater than 10 pack-years, unstable asthma, respiratory tract infection or exacerbation less than 6 weeks before enrollment, postbronchodilator FEV1 of less than 60% of predicted value, and other pulmonary disease. The institutional review board at each center approved all study protocols. After the run-in treatment period with inhaled fluticasone, sputum induction was performed for differential cell counts, followed by flexible fiberoptic bronchoscopy 3 to 10 days later after achievement of local anesthesia, as previously described.E1 This was performed before randomization to treatment with clarithromycin or placebo in the parent trial. During bronchoscopy, 3 bronchial brushings with triple-lumen protected specimen brushes were obtained for microbiome analysis, primarily from the upper or middle lobes of the lung, and immediately placed into 2 mL of RNALater. Samples were placed at 4°C for at least 24 hours before shipping on blue ice to the University of California, San Francisco. In addition, healthy nonsmoking adults were enrolled separately at the University of California, San Francisco, and ruled out for asthma and atopy based on the presence of normal spirometric results, normal bronchial reactivity to methacholine (PC20 >16 mg/mL), and negative responses to skin testing with 10 mixes of allergens common to northern California.

Sample preparations and array processing 

Sterilized forceps were used to transfer 3 protected specimen bronchial brushes to a single lysing matrix B tube containing 0.1-mm silica spheres (MP Biomedicals, Solon, Ohio). Original collection tubes were centrifuged at 12,000g, RNALater was removed, and Buffer RLT (650 μL, AllPrep kit; Qiagen, Valencia, Calif) was used to resuspend the pellet before transfer into the lysing matrix B tube containing the brushes. Bead beating at 5.5 m/s for 30 seconds was carried out with the FastPrep system (MP Biomedicals). Tubes were then centrifuged at 2,000g for 2 minutes, and the supernatant was transferred to a DNA spin column of the AllPrep Kit and further processed according to kit instructions.

Samples were screened by means of PCR for the presence of a 16S rRNA gene product with 100 ng of extracted DNA and 15 μmol each of the universal 16S rRNA primersE2 Bact-27F (5′-AGAGTTTGATCCTGGCTCAG-3′) and Bact-1492R (5′-GGTTACCTTGTTACGACTT-3′) in 50-μL volume reactions at an annealing temperature of 53°C for 40 cycles in an Eppendorf MasterCycler gradient PCR machine. In addition to negative controls, parallel poison control PCR reactions (with Pseudomonas aeruginosa DNA) were performed to rule out PCR inhibition for those samples with no evidence of 16S rRNA PCR product on a 1% ethidium bromide–stained agarose gel. Those with a PCR-positive screen for a 16S rRNA PCR product then proceeded to sample preparation for array hybridization and processing (Fig E1). For this, eight 25-μL PCR reactions, each containing 100 ng of extracted sample DNA and 8 μmol of each of Bact-27F and Bact-1492R primers, were performed in 30 cycles of amplification across a gradient of annealing temperatures (48°C to 58°C). 16S rRNA amplicons for each sample were then pooled, purified, and gel quantified before array hybridization (E-gel 2% and E-gel Low Range Quantitative Ladder, Invitrogen). A standard concentration of 250 ng (or 100 ng in 7 samples, as discussed in the Results section) of DNAse-fragmented and biotin-labeled 16S rRNA amplicons from each sample spiked with known concentrations of synthetic 16S rRNA gene fragments and non-16S rRNA gene fragments as internal controls for interarray normalization was hybridized to each array, followed by washing, staining, and scanning of arrays (GeneArray scanner; Affymetrix, Santa Clara, Calif) as previously described.E3., E4.

Q-PCR 

Triplicate 25-μL reactions were performed in 21 samples by using the QuantiTect SYBR Green PCR master mix kit (Qiagen), 10 ng of extracted DNA, and the 16S rRNA Q-PCR primersE5 338-F (5′-ACTCCTACGGGAGGCAGCAG-3′) and 518-R (5′- ATTACCGCGGCTGCTGG-3′) at a final concentration of 0.3 μmol/L under the following conditions: 95°C for 10 minutes, followed by 40 cycles of denaturation at 95°C for 30 seconds, annealing at 53°C for 30 seconds, extension at 72°C for 30 seconds, and final extension at 72°C for 10 minutes. Standard curves at serial log concentrations of 1 × 102 to 1 × 108 were generated from a PCR product obtained by using the universal 16S rRNA primers Bact-27F and Bact-1492R to calculate total 16S rRNA gene copy numbers. Regression coefficients (r2) for all standard curves were greater than 0.99. Total 16S rRNA copy numbers for each sample were normalized to mammalian β-actin copy numbers in each sample, which were determined by means of Q-PCR with β-actin primers, as previously described.E6

Microarray data analysis 

After scaling and normalization of the raw data, relatively conservative criteria were applied to determine the presence or absence of bacterial taxa by using previously described methods.E3., E4. Briefly, taxon presence was determined based on a positive-fraction value corresponding to 90% of the probe sets for an individual taxon. Each taxon probe set is comprised of a minimum of 11 probe pairs (a perfect match and a mismatch probe differing at a single central oligonucleotide) to distinguish taxa on the basis of at least 2 discriminatory loci on the 16S rRNA gene. Relative abundance estimates were based on averaged, log-transformed fluorescence intensities for each taxon. Sorenson distance measures were calculated by using the software package PC-ORD version 5. Three samples from asthmatic patients were identified as statistical outliers with a mean relative Sorensen distance of greater than 2 SDs from the pooled mean and were not included in the array-based analyses of community structure and diversity.

16S rRNA clone library-sequencing 

Clone libraries were constructed by using the same 16S rRNA amplicon pool that was hybridized to the array. 16S rRNA PCR products were gel purified, and 3′ adenine overhangs were reintroduced after amplification through Taq DNA polymerase for improved cloning efficiency. Purified amplicons were ligated to the pCR4-TOPO vector by using the TOPO TA Cloning Kit (Invitrogen). Ligation mixtures were electroporated into One Shot TOP10 Electrocomp Escherichia coli cells (Invitrogen) and plated on selective agar. Two hundred ninety-nine to 374 clones per library were sequenced on an ABI PRISM 3730 capillary DNA sequencer (Applied Biosystems, Foster City, Calif) at the Department of Energy’s Joint Genome Institute (Walnut Creek, Calif), according to standard Joint Genome Institute protocols (www.jgi.doe.gov). Bidirectional sequence reads were end paired, trimmed for vector and PCR primer sequence, and filtered for quality. This yielded 1,997 sequences for 6 asthmatic libraries, which were aligned by using the NAST algorithm for multiple sequence alignment, checked for chimeras with Bellepheron, and classified by means of G2-chip taxonomy, all available as Greengenes Web-based tools (http://greengenes.lbl.gov). Taxa identified by means of clone library sequencing and by means of corresponding array data in these subjects were compared to determine the number of phylotypes detected by using the array, by sequencing alone, or by both methods.

Back to Article Outline

Appendix E1 

Collaborators 

The following Asthma Clinical Research Network sites and investigators participated in the parent MIA trial, which obtained the clinical data analyzed in this study and additionally obtained the bronchial brushings analyzed in this study only:

Brigham & Women’s Hospital, Boston, Mass: E. Israel, M. E. Wechsler

National Jewish Health, Denver, Colo: R. J. Martin, R. M. Cherniack, S. J. Szefler, E. R. Sutherland

University of Wisconsin School of Medicine and Public Health, Madison, Wis: R. F. Lemanske, C. A. Sorkness, N. N. Jarjour, L. Denlinger.

University of California, San Francisco: H. A. Boushey, J. V. Fahy, S. C. Lazarus

Columbia University College of Physicians and Surgeons, New York, NY: E. DiMango

Duke University, Durham, NC: M. Kraft,

University of Texas Medical Branch, Galveston, Tex: W. J. Calhoun, B. T. Ameredes

Washington University, School of Medicine, St Louis, Mo: M. Castro, M. Walter

University of California, San Diego Medical Center, Calif: J. Ramsdell, S. I. Wasserman

Wake Forest University, School of Medicine, Winston-Salem, NC: E. Bleecker, D. Meyers, S. P. Peters, W. C. Moore, R. Pascual

Pennsylvania State College of Medicine, Hershey, Pa: V. M. Chinchilli, T. J. Craig, N. Icitovic, T. S. King

Back to Article Outline

Fig E1. 

Back to Article Outline

Supplementary data 

Tables E1 and E2.

Back to Article Outline

References 

  1. Bisgaard H, Hermansen MN, Buchvald F, Loland L, Halkjaer LB, Bonnelykke K, et al. Childhood asthma after bacterial colonization of the airway in neonates. N Engl J Med. 2007;357:1487–1495
  2. Martin RJ, Kraft M, Chu HW, Berns EA, Cassell GH. A link between chronic asthma and chronic infection. J Allergy Clin Immunol. 2001;107:595–601
  3. Black PN, Blasi F, Jenkins CR, Scicchitano R, Mills GD, Rubinfeld AR, et al. Trial of roxithromycin in subjects with asthma and serological evidence of infection with Chlamydia pneumoniae. Am J Respir Crit Care Med. 2001;164:536–541
  4. Kostadima E, Tsiodras S, Alexopoulos EI, Kaditis AG, Mavrou I, Georgatou N, et al. Clarithromycin reduces the severity of bronchial hyperresponsiveness in patients with asthma. Eur Respir J. 2004;23:714–717
  5. Kraft M, Cassell GH, Pak J, Martin RJ. Mycoplasma pneumoniae and Chlamydia pneumoniae in asthma: effect of clarithromycin. Chest. 2002;121:1782–1788
  6. Flanagan JL, Brodie EL, Weng L, Lynch SV, Garcia O, Brown R, et al. Loss of bacterial diversity during antibiotic treatment of intubated patients colonized with Pseudomonas aeruginosa. J Clin Microbiol. 2007;45:1954–1962
  7. Hilty M, Burke C, Pedro H, Cardenas P, Bush A, Bossley C, et al. Disordered microbial communities in asthmatic airways. PLoS One. 2010;5:e8578
  8. Klepac-Ceraj V, Lemon KP, Martin TR, Allgaier M, Kembel SW, Knapp AA, et al. Relationship between cystic fibrosis respiratory tract bacterial communities and age, genotype, antibiotics and Pseudomonas aeruginosa. Environ Microbiol. 2010;12:1293–1303
  9. Ledder RG, Gilbert P, Huws SA, Aarons L, Ashley MP, Hull PS, et al. Molecular analysis of the subgingival microbiota in health and disease. Appl Environ Microbiol. 2007;73:516–523
  10. Turnbaugh PJ, Hamady M, Yatsunenko T, Cantarel BL, Duncan A, Ley RE, et al. A core gut microbiome in obese and lean twins. Nature. 2009;457:480–484
  11. Wen L, Ley RE, Volchkov PY, Stranges PB, Avanesyan L, Stonebraker AC, et al. Innate immunity and intestinal microbiota in the development of Type 1 diabetes. Nature. 2008;455:1109–1113
  12. Brodie EL, DeSantis TZ, Parker JP, Zubietta IX, Piceno YM, Andersen GL. Urban aerosols harbor diverse and dynamic bacterial populations. Proc Natl Acad Sci U S A. 2007;104:299–304
  13. Green BD, Keller M. Capturing the uncultivated majority. Curr Opin Biotechnol. 2006;17:236–240
  14. Lin LH, Wang PL, Rumble D, Lippmann-Pipke J, Boice E, Pratt LM, et al. Long-term sustainability of a high-energy, low-diversity crustal biome. Science. 2006;314:479–482
  15. Brodie EL, Desantis TZ, Joyner DC, Baek SM, Larsen JT, Andersen GL, et al. Application of a high-density oligonucleotide microarray approach to study bacterial population dynamics during uranium reduction and reoxidation. Appl Environ Microbiol. 2006;72:6288–6298
  16. Sagaram US, DeAngelis KM, Trivedi P, Andersen GL, Lu SE, Wang N. Bacterial diversity analysis of Huanglongbing pathogen-infected citrus, using PhyloChip arrays and 16S rRNA gene clone library sequencing. Appl Environ Microbiol. 2009;75:1566–1574
  17. Schonmann S, Loy A, Wimmersberger C, Sobek J, Aquino C, Vandamme P, et al. 16S rRNA gene-based phylogenetic microarray for simultaneous identification of members of the genus Burkholderia. Environ Microbiol. 2009;11:779–800
  18. La Duc MT, Osman S, Vaishampayan P, Piceno Y, Andersen G, Spry JA, et al. Comprehensive census of bacteria in clean rooms by using DNA microarray and cloning methods. Appl Environ Microbiol. 2009;75:6559–6567
  19. Huang YJ, Kim E, Cox MJ, Brodie EL, Brown R, Wiener-Kronish JP, et al. A persistent and diverse airway microbiota present during chronic obstructive pulmonary disease exacerbations. OMICS. 2010;14:9–59
  20. Cox MJ, Huang YJ, Fujimura KE, Liu JT, McKean M, Boushey HA, et al. Lactobacillus casei abundance is associated with profound shifts in the infant gut microbiome. PLoS One. 2010;5:e8745
  21. DeSantis TZ, Brodie EL, Moberg JP, Zubieta IX, Piceno YM, Andersen GL. High-density universal 16S rRNA microarray analysis reveals broader diversity than typical clone library when sampling the environment. Microb Ecol. 2007;53:371–383
  22. Dethlefsen L, Huse S, Sogin ML, Relman DA. The pervasive effects of an antibiotic on the human gut microbiota, as revealed by deep 16S rRNA sequencing. PLoS Biol. 2008;6:e280
  23. Sutherland ER, King TS, Icitovic N, Ameredes BT, Bleecker E, Boushey HA, et al. A trial of clarithromycin for the treatment of suboptimally controlled asthma. J Allergy Clin Immunol. 2010;126:747–753
  24. Huang YJ, Lynch SV, Brodie EL, Wiener-Kronish JP, Sutherland E, Investigators A, et al. Detection of diverse bacterial populations in asthmatic airways by 16S rRNA PhyloChip [abstract]. Am J Respir Crit Care Med. 2008;177:A335
  25. Lynch SV, Huang YJ, Nelson CE, Brodie EL, DeSantis TZ, Baek MS, et al. Diverse bacterial communities in the airways of patients with sub-optimally controlled asthma [abstract]. Am J Respir Crit Care Med. 2010;181:A5093
  26. Juniper EF, O’Byrne PM, Guyatt GH, Ferrie PJ, King DR. Development and validation of a questionnaire to measure asthma control. Eur Respir J. 1999;14:902–907
  27. Lane DJ. 16S/23S rRNA sequencing. Chichester (United Kingdom): Wiley; 1991;
  28. Ramette A. Multivariate analyses in microbial ecology. FEMS Microbiol Ecol. 2007;62:142–160
  29. Storey JD, Tibshirani R. Statistical significance for genomewide studies. Proc Natl Acad Sci U S A. 2003;100:9440–9445
  30. Hill MO. Diversity and evenness: a unifying notation and its consequences. Ecology. 1973;54:428–432
  31. McCune B, Mefford MJ. PC-ORD v5. Multivariate analysis of ecological data. Gleneden Beach (OR): MJM Software Design. 2006;
  32. DeSantis TZ, Hugenholtz P, Keller K, Brodie EL, Larsen N, Piceno YM, et al. NAST: a multiple sequence alignment server for comparative analysis of 16S rRNA genes. Nucleic Acids Res. 2006;34:W394–W399
  33. Huber T, Faulkner G, Hugenholtz P. Bellerophon: a program to detect chimeric sequences in multiple sequence alignments. Bioinformatics. 2004;20:2317–2319
  34. DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, et al. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006;72:5069–5072
  35. Price MN, Dehal PS, Arkin AP. FastTree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol Biol Evol. 2009;26:1641–1650
  36. Letunic I, Bork P. Interactive tree of life (iTOL): an online tool for phylogenetic tree display and annotation. Bioinformatics. 2007;23:127–128
  37. Nelson KE, Weinstock GM, Highlander SK, Worley KC, Creasy HH, Wortman JR, et al. A catalog of reference genomes from the human microbiome. Science. 2010;328:994–999
  38. Chen T, Yu WH, Izard J, Baranova OV, Lakshmanan A, Dewhirst FE. The Human Oral Microbiome Database: a web accessible resource for investigating oral microbe taxonomic and genomic information. Oxford: Database; 2010;baq013
  39. Haldar P, Pavord ID, Shaw DE, Berry MA, Thomas M, Brightling CE, et al. Cluster analysis and clinical asthma phenotypes. Am J Respir Crit Care Med. 2008;178:218–224
  40. Moore WC, Meyers DA, Wenzel SE, Teague WG, Li H, Li X, et al. Identification of asthma phenotypes using cluster analysis in the severe asthma research program. Am J Respir Crit Care Med. 2010;181:315–323
  41. Woodruff PG, Modrek B, Choy DF, Jia G, Abbas AR, Ellwanger A, et al. T-helper type 2-driven inflammation defines major subphenotypes of asthma. Am J Respir Crit Care Med. 2009;180:388–395
  42. Willing B, Halfvarson J, Dicksved J, Rosenquist M, Jarnerot G, Engstrand L, et al. Twin studies reveal specific imbalances in the mucosa-associated microbiota of patients with ileal Crohn’s disease. Inflamm Bowel Dis. 2009;15:653–660
  43. Dethlefsen L, McFall-Ngai M, Relman DA. An ecological and evolutionary perspective on human-microbe mutualism and disease. Nature. 2007;449:811–818
  44. Lloyd CM, Saglani S. Asthma and allergy: the emerging epithelium. Nat Med. 2010;16:273–274
  45. Beaumont HJ, van Schooten B, Lens SI, Westerhoff HV, van Spanning RJ. Nitrosomonas europaea expresses a nitric oxide reductase during nitrification. J Bacteriol. 2004;186:4417–4421
  46. Payne DN, Adcock IM, Wilson NM, Oates T, Scallan M, Bush A. Relationship between exhaled nitric oxide and mucosal eosinophilic inflammation in children with difficult asthma, after treatment with oral prednisolone. Am J Respir Crit Care Med. 2001;164:1376–1381
  47. Jatakanon A, Lim S, Kharitonov SA, Chung KF, Barnes PJ. Correlation between exhaled nitric oxide, sputum eosinophils, and methacholine responsiveness in patients with mild asthma. Thorax. 1998;53:91–95
  48. Duncan SH, Richardson AJ, Kaul P, Holmes RP, Allison MJ, Stewart CS. Oxalobacter formigenes and its potential role in human health. Appl Environ Microbiol. 2002;68:3841–3847
  49. Coenye T, Goris J, Spilker T, Vandamme P, LiPuma JJ. Characterization of unusual bacteria isolated from respiratory secretions of cystic fibrosis patients and description of Inquilinus limosus gen. nov., sp. nov. J Clin Microbiol. 2002;40:2062–2069
  50. Gohler A, Xiong G, Paulsen S, Trentmann G, Maser E. Testosterone-inducible regulator is a kinase that drives steroid sensing and metabolism in Comamonas testosteroni. J Biol Chem. 2008;283:17380–17390
  51. Horinouchi M, Kurita T, Yamamoto T, Hatori E, Hayashi T, Kudo T. Steroid degradation gene cluster of Comamonas testosteroni consisting of 18 putative genes from meta-cleavage enzyme gene tesB to regulator gene tesR. Biochem Biophys Res Commun. 2004;324:597–604
  52. Campbell JD, Borish L, Haselkorn T, Rasouliyan L, Lee JH, Wenzel SE, et al. The response to combination therapy treatment regiments in severe/difficult-to-treat asthma. Eur Respir J. 2008;32:1237–1242
  53. Wenzel SE, Szefler SJ, Leung DY, Sloan SI, Rex MD, Martin RJ. Bronchoscopic evaluation of severe asthma. Persistent inflammation associated with high dose glucocorticoids. Am J Respir Crit Care Med. 1997;156:737–743
  54. Akbari O, Stock P, Meyer E, Kronenberg M, Sidobre S, Nakayama T, et al. Essential role of NKT cells producing IL-4 and IL-13 in the development of allergen-induced airway hyperreactivity. Nat Med. 2003;9:582–588
  55. Kinjo Y, Wu D, Kim G, Xing GW, Poles MA, Ho DD, et al. Recognition of bacterial glycosphingolipids by natural killer T cells. Nature. 2005;434:520–525
  56. Meyer EH, Goya S, Akbari O, Berry GJ, Savage PB, Kronenberg M, et al. Glycolipid activation of invariant T cell receptor+ NK T cells is sufficient to induce airway hyperreactivity independent of conventional CD4+ T cells. Proc Natl Acad Sci U S A. 2006;103:2782–2787
  57. Akbari O, Faul JL, Hoyte EG, Berry GJ, Wahlstrom J, Kronenberg M, et al. CD4+ invariant T-cell-receptor+ natural killer T cells in bronchial asthma. N Engl J Med. 2006;354:1117–1129
  58. Vijayanand P, Seumois G, Pickard C, Powell RM, Angco G, Sammut D, et al. Invariant natural killer T cells in asthma and chronic obstructive pulmonary disease. N Engl J Med. 2007;356:1410–1422
  59. Drummond MB, Dasenbrook EC, Pitz MW, Murphy DJ, Fan E. Inhaled corticosteroids in patients with stable chronic obstructive pulmonary disease: a systematic review and meta-analysis. JAMA. 2008;300:2407–2416
  60. Qin J, Li R, Raes J, Arumugam M, Burgdorf KS, Manichanh C, et al. A human gut microbial gene catalogue established by metagenomic sequencing. Nature. 2010;464:59–65
  61. Amayasu H, Yoshida S, Ebana S, Yamamoto Y, Nishikawa T, Shoji T, et al. Clarithromycin suppresses bronchial hyperresponsiveness associated with eosinophilic inflammation in patients with asthma. Ann Allergy Asthma Immunol. 2000;84:594–598

Back to Article Outline

References 

  1. Kraft M, Martin RJ, Lazarus SC, Fahy JV, Boushey HA, Lemanske RF, et al. Airway tissue mast cells in persistent asthma: predictor of treatment failure when patients discontinue inhaled corticosteroids. Chest. 2003;124:42–50
  2. Lane DJ. 16S/23S rRNA sequencing. Chichester (United Kingdom): Wiley; 1991;
  3. Flanagan JL, Brodie EL, Weng L, Lynch SV, Garcia O, Brown R, et al. Loss of bacterial diversity during antibiotic treatment of intubated patients colonized with Pseudomonas aeruginosa. J Clin Microbiol. 2007;45:1954–1962
  4. DeSantis TZ, Brodie EL, Moberg JP, Zubieta IX, Piceno YM, Andersen GL. High-density universal 16s rRNA microarray analysis reveals broader diversity than typical clone library when sampling the environment. Microb Ecol. 2007;53:371–383
  5. Einen J, Thorseth IH, Ovreas L. Enumeration of archaea and bacteria in seafloor basalt using real-time quantitative pcr and fluorescence microscopy. FEMS Microbiol Lett. 2008;282:182–187
  6. Matsui K, Tanaka N, Nishikawa A. Lipopolysaccharide of Haemophilus influenzae induces interleukin-5 mRNA expression in human peripheral blood mononuclear cells. J Interferon Cytokine Res. 2001;21:439–443

 Supported by the National Heart, Lung, and Blood Institute (NHLBI; U10 HL 074204) and by the Strategic Asthma Basic Research Center at the University of California, San Francisco, supported by the Sandler Family Foundation. Y.J.H. was funded by National Institutes of Health (NIH)/NHLBI grant T32 HL007185 and by a University of California Tobacco-related Disease Research Program award (17FT-0040). C.E.N. is funded by NSF 0709975 (to C.E.N. and J.M. Melack). S.V.L. is funded by NIH/National Institute of Allergy and Infectious Diseases (NIAID) grant U01 AI075410. E.L.B., T.Z.D., and J.B. are funded under the auspices of the University of California under contract number DOE DE-AC02-05CH11231.

 Disclosure of potential conflict of interest: T. Z. DeSantis is a part-time employee of PhyloTech, Inc. J. Bristow receives research support from the DOE Joint Genome Institute. J. P. Weiner-Kronish is a board member of the Foundation of Anesthesia Education and Research. E. R. Sutherland is an advisor and DSMB member for GlaxoSmithKline, is an advisor for Dey, is a DSMB member for Merck, and receives research support from the National Institutes of Health, Novartis, and Boehringer-Ingelheim. R. J. Martin receives research support from the National Heart, Lung, and Blood Institute of the National Institutes of Health. M. Castro is a consultant for NKT Therapeutics, Schering-Plough, Asthmatx, and Cephalon; is on the Advisory Board for Genentech; is on the speakers’ bureau for Astra-Zeneca, Boehringer-Ingelheim, Pfizer, Merck, and GlaxoSmithKline; has received grant support from Asthmatx, Amgen, Ception, Genentech, MedImmune, Merck, Novartis, the National Institutes of Health, GlaxoSmithKline, and the American Lung Association; and has received royalties from Elsevier. L. C. Denlinger receives research support from the National Institutes of Health (NIH)–National Heart, Lung, and Blood Institute (NHLBI). M. Kraft has received research support from GlaxoSmithKline, Merck, Asthmatx, GE Healthcare, Novartis, and Genentech. S. P. Peters receives grant support from the NIH-NHLBI. S. I. Wasserman has provided legal consultation services/expert witness testimony in cases related to mold toxicity and transfer factor and is president of the American Board of Allergy and Immunology. M. E. Wechsler receives research support from the NHLBI. H. A. Boushey is an ad-hoc consultant for Kalobios, is on the advisory committee for Pharmaxis, is on ad-hoc advisory committees for Glaxo-SmithKline and Merck, and receives research support from GlaxoSmithKline. S. V. Lynch receives research support from the National Institutes of Health. The rest of the authors have declared that they have no conflict of interest.

PII: S0091-6749(10)01760-4

doi:10.1016/j.jaci.2010.10.048

The Journal of Allergy and Clinical Immunology
Volume 127, Issue 2 , Pages 372-381.e3, February 2011