- Open Access
Enrichment of lung microbiome with supraglottic taxa is associated with increased pulmonary inflammation
Microbiome volume 1, Article number: 19 (2013)
The lung microbiome of healthy individuals frequently harbors oral organisms. Despite evidence that microaspiration is commonly associated with smoking-related lung diseases, the effects of lung microbiome enrichment with upper airway taxa on inflammation has not been studied. We hypothesize that the presence of oral microorganisms in the lung microbiome is associated with enhanced pulmonary inflammation. To test this, we sampled bronchoalveolar lavage (BAL) from the lower airways of 29 asymptomatic subjects (nine never-smokers, 14 former-smokers, and six current-smokers). We quantified, amplified, and sequenced 16S rRNA genes from BAL samples by qPCR and 454 sequencing. Pulmonary inflammation was assessed by exhaled nitric oxide (eNO), BAL lymphocytes, and neutrophils.
BAL had lower total 16S than supraglottic samples and higher than saline background. Bacterial communities in the lower airway clustered in two distinct groups that we designated as pneumotypes. The rRNA gene concentration and microbial community of the first pneumotype was similar to that of the saline background. The second pneumotype had higher rRNA gene concentration and higher relative abundance of supraglottic-characteristic taxa (SCT), such as Veillonella and Prevotella, and we called it pneumotypeSCT. Smoking had no effect on pneumotype allocation, α, or β diversity. PneumotypeSCT was associated with higher BAL lymphocyte-count (P= 0.007), BAL neutrophil-count (P= 0.034), and eNO (P= 0.022).
A pneumotype with high relative abundance of supraglottic-characteristic taxa is associated with enhanced subclinical lung inflammation.
The gut microbiome modulates host mucosal immune response [1, 2]. In contrast, despite emerging data on airway microbiota, little is known about the role of the lung microbiome in modulating pulmonary mucosal immune response. Along the human airways, structures above the vocal cords are exposed to high bacterial burden producing contamination of lower airway samples with oropharyngeal secretions [3, 4]. Classical culture techniques fail to fully describe microbial communities of the lower airways due to difficulties in growing fastidious bacteria, thus identification using marker gene sequencing is a more promising approach [5–11]. Such studies identified bacterial rRNA genes of oral cavity anaerobes such as Prevotella and Veillonella in the lower airways of normal individuals [6, 8, 12–16].
Microaspiration of small volumes of oropharyngeal secretions occurs in healthy people . The risk for microaspiration is increased in smoking-related lung disease due to reduced coordination of breathing with swallowing and gastro-esophageal reflux [18, 19]. Both microaspiration and impaired mechanical clearance in smokers may lead to increased lower airway colonization with oral microbiota. Prior studies suggest that in moderate to severe chronic obstructive pulmonary disease (COPD), the lung microbiome is enriched with bacteria or bacterial products common to the oral cavity [7, 9, 20, 21]. However, these studies have focused on advanced COPD patients, in whom frequent antibiotic and corticosteroid use may affect the bacterial communities of the lower airways. Studies of early disease and asymptomatic cases would avoid these potential confounding effects.
Increased toll-like receptor signaling has been associated with pulmonary inflammation in advanced COPD, offering a mechanisms by which microbial inhabitants in the lung might be relevant for the development of smoking-related lung injury [22, 23]. Furthermore, randomized clinical trials provide indications that antibiotics, especially macrolides, may reduce COPD exacerbations [24, 25]. However, it is unclear whether the beneficial effect is due to the antibiotic or to the anti-inflammatory properties of these drugs. An in-depth understanding of the lung microbiome and its association with mucosal inflammatory response is needed to understand potential mechanisms of lung mucosal immune regulation. Here, we hypothesized that the enrichment of the lung microbiome of asymptomatic subjects with supraglottic-characteristic taxa is associated with lung inflammation. To avoid potential confounders, we selected asymptomatic smokers and never-smokers with preserved lung function and no recent exposure to antibiotics or immune modulators, to assess lung microbiome characteristics.
We enrolled 29 asymptomatic subjects(nine never-smokers, 14 former-smokers, and six current-smokers) for research bronchoscopy. All participants signed informed consent to take part in this study and the research protocol was approved by the human subjects review committees of New York University and by Bellevue Hospital Center (New York, NY) institutional review boards. For all participants, exclusion criteria were: FEV1<70%; recent treatment with antibiotics or steroids in the prior 3 months; cardiovascular disease (abnormal EKG, known or suspected coronary artery disease, or congestive heart failure); diabetes mellitus; renal or liver disease; lung cancer; heavy alcohol use (>6 beers daily). During the week prior to bronchoscopy, all participants underwent research pulmonary function testing for physiologic phenotyping, which included spirometry, lung volumes, diffusion, impulse oscillometry (IOS), and exhaled nitric oxide (eNO) as a non-invasive measure of airway inflammation . In eight participants, eNO could not be reliably measured due to high levels of ambient NO that was associated with inability to identify a clear plateau on the exhaled NO vs. time tracing.
For all 29 participants, we used a nasal approach to perform research bronchoscopy. Since we previously observed low bacterial loads in nasal brushing , we used an optimized nasal bronchoscopy in an attempt to minimize the potential for contamination of the bronchoscope channel with upper airway microbiota. Given the concern for potential carry-over of upper airway microbiome to the lower airways, we avoid suctioning until the scope was in position for sampling. If during bronchoscopy, visualization was obscured, we instilled saline through the bronchoscope channel to clear the optics. In a subgroup of 15 participants, we used two bronchoscopes to evaluate potential carry-over. In this subgroup, we obtained a background sample from sterile saline and sterile normal saline passed through the bronchoscope channel. We then passed the first bronchoscope until vocal cords were visualized to obtain a supraglottic sample. The second bronchoscope was passed to obtain BAL in two different segments. The initial BAL was from a segment of the left lung (lingula) and the second BAL was from a segment of the right lung (right middle lobe, Additional file 1: Figure S1). We used three aliquots of 50 mL of saline for each BAL. On average, 50% of the instilled volume was recovered. The BAL samples were immediately placed on ice and processed within 30 min of acquisition. Total cell count and BAL cell differentials (500 cells counted) were performed to assess counts of macrophages, lymphocytes, and neutrophils. BAL fluid aliquots were frozen at −80°C for DNA isolation.
Bacterial DNA quantification and 454 pyrosequencing
DNA was isolated from acellular BAL fluid after centrifugation at 500 g for 10 min at 4°C. Three methods were used to ensure uniform bacterial lysis: a freeze-thaw cycle, lysozyme, and a 56°C heat step at the beginning of the DNA isolation process. DNA was then extracted with an ion exchange column (Qiagen). Total bacterial and human DNA levels were determined by quantitative PCR (qPCR). We performed high-throughput sequencing of the bacterial 16S rRNA gene using the Roche 454 amplicon sequencing protocol. Each sample was individually barcoded using a unique 12nt sequence during amplification of the V1-V2 regions of the 16S gene . Further details described in Additional file 2 on line supplement.
Upstream informatics analysis of the 16S sequences
The obtained 16S rRNA sequences were analyzed using the Quantitative Insights into Microbial Ecology (QIIME) pipeline for analysis of community sequence data . Processing consisted of the following steps: (1) demultiplexing and filtering of short (<150 nt) and low quality reads; (2) de novo clustering of the sequences into operational taxonomic units (OTUs) with the UCLUST program using a 97% similarity threshold ; (3) taxonomical assignment of each OTU by running the RDP Classifier  at 80% bootstrap confidence on a selected representative sequence from each OTU; (4) alignment of representative sequences using PyNAST  with the Greengenes core-set alignment template; (5) phylogenetic tree reconstruction from the representative sequences for each OTUs using the FASSTTREE program ; and (6) UniFrac distance calculations . For each sample, the proportion of reads at the OTU or genus levels was used as a measure of the relative abundance of each type of bacteria in a specimen. The exact commands run are available in the Supplement.
Univariate analysis of association between microbiome and immunological phenotype
The absolute OTU sequence counts were normalized to obtain the relative abundances of the microbiota in each sample. These relative abundances at 97% OTU similarity and each of the five higher taxonomic levels (phylum, class, order, family, genus) were tested for univariate associations with clinical variables. To decrease the number of features, we only focused on major taxa and OTUs, defined as those having mean relative abundance >1% in at least one sample.
Since the distributions of microbiome data are non-normal, and no distribution-specific tests are available, we used non-parametric tests of association. For association with discrete factors, we used either the Mann–Whitney test (in the case of two categories) or the Kruskal Wallis ANOVA (in case of >2 categories). For tests of association with continuous variables, we used non-parametric Spearman correlation tests. False discovery rate (FDR) was used to control for multiple testing . Weighted UniFrac was used to measure β diversity of bacterial communities and to perform principal coordinate analysis (PCoA) . We used the ade4 package in R to PCoA on weighted Unifrac distances . To avoid negative eigenvalues in the analysis, we used the Cailliez method to convert the weighted Unifrac distance matrix into a closest corresponding matrix with Euclidean properties, which was further used for PCoA . Hierarchical clustering on the relative abundance profiles was used to establish deep branching, which was interpreted as evidence for distinct pneumotypes.
Classification of streptococcus OTUs
For classification of Streptococcus OTUs, 16S rRNA sequences were aligned to the Greengenes database using online Blast tool (http://greengenes.lbl.gov/cgi-bin/nph-blast_interface.cgi). Top hits were retained. A perfect match was used for final classification if found, otherwise top hit is reported. In case two species we found more than once, both are reported.
Table 1 shows demographics, clinical characteristics, pulmonary function, and BAL cell counts for the 29 participants of this cohort. Enrolled subjects had normal FEV1 (median (IQR) = 94(86–104)% predicted) and FVC (97(85–104)% predicted), and normal or mild obstruction as quantified by FEV1/VC (73.4 (70.0-83.4)). Total airway resistance as assessed by IOS also was in the normal range (R5 = 3.0(2.5-3.4) cmH2O/L/s). Diffusion was normal (91(80–105)% predicted). Exhaled NO ranged from normal to elevated (range, 6.9-46.3 ppb). BAL cell differentials were normal with major predominance of alveolar macrophages and no neutrophilia or eosinophilia. In total, these data show that our cohort of asymptomatic subjects had no major pulmonary function abnormalities or obvious clinical evidence of disease.
Bacterial 16S rRNA gene qPCR and sequencing
We next examined bacterial load in the upper airways, lower airways, and background. Supraglottic samples had the highest 16S rRNA gene concentration (3,237,500 (276,625-6,262,500) copies/μL, Figure 1). 16S concentrations in BAL were approximately 50-fold lower (68,875 (29,475-183,312) copies/μL, P<0.001). Despite extensive overlap in bacterial load between background and BAL samples, background samples had lower 16S rRNA gene concentration (41,195 (18,035-79,875) copies/μL, P= 0.018). After trimming and quality control filtering, we obtained 641,847 rDNA sequences.
454 sequencing yield no significant difference in the number of reads per sample between rRNA obtained from background, BAL, and supraglottic specimens (median = 6,116, 5,532, and 7,961 high quality reads per sample, respectively, P= 0.12). The number of OTUs observed at 97% homology was significantly lower in background samples (93 (79–112)) compared with BAL (136 (98–158), P<0.001) and supraglottic samples (112 (102–140), P= 0.004).
To ensure that research bronchoscopy did not systematically carry-over supraglottic secretions into the lower airway, we compared the matched sequencing samples from 15 subjects (three never-smokers and 12 smokers) from the supraglottic area, with the first BAL obtained from a segment in the lingula and the second BAL obtained from a right middle lobe segment. UniFrac distances were calculated between supraglottic and both first and second BAL (Figure 2). There were no significantly different UniFrac distances to supraglottic microbiome between first and second BAL, providing evidence that systematic carry-over did not occur. Analysis using unweighted UniFrac analysis yielded similar results. Figure 3 shows data from the 15 patients with complete supraglottic, first BAL, and second BAL data. The relative abundance of highly represented taxa in the supraglottic (Prevotella and Veillonella) and saline background (Staphylococcus and Propionibacterium) are shown in the bar graph with 16S rRNA copies/mL and Unifrac distance shown below each sample. There was no consistent ‘dilution’ of bacterial load and/or relative abundances of supraglottic taxa in subsequent BAL samples. As an example, in case S10, the first BAL had lower relative abundance of Prevotella and Veillonella and lower total bacterial load by qPCR (6,902 copies/mL) than the supraglottic sample (bacterial load = 34,250,000 copies/mL). Compared to first BAL, the second BAL had increased relative abundance of Prevotella and Veillonella (total bacterial load 103,300 copies/mL). Since the bronchoscope channel was not in contact with the upper airway between the first and second BAL, the increased relative abundance of Prevotella/Veillonella and total bacterial load between the first and second BAL in this example could not have been due to carry-over of supraglottic secretions contaminating the bronchoscope. From these observations, we conclude that enrichment of lower airway microbiome with supraglottic-characteristic taxa (SCT) most likely results from microaspiration, rather than from bronchoscopic carry-over.
Distinct pneumotypes in the lung microbiome
We next asked whether distinct microbiomes could be identified in the lower airways. We used unsupervised hierarchical clustering to assess the structure of the microbiome in the background, supraglottic, and BAL samples. Background sterile saline contained high relative abundances of Propionibacterium, Staphylococcus, Corynebacterium, Stenotrophomonas, and Pseudomonas (Figure 4, panel A). Supraglottic samples contained high relative abundances of Prevotella, Veillonella, Streptococcus, Fusobacterium, and Porphyromonas. PCoA of weighted UniFrac distances confirms that the supraglottic and background samples are distinct (panel B). Similar results were found when unweighted UniFrac was used (data not shown).
Unsupervised hierarchical clustering of BAL samples exhibited deep branches into two distinct clusters that we called ‘pneumotypes’ (Figure 4, panel C). The first pneumotype is similar to background samples with high relative abundance of Propionibacterium (panel C). This similarity to background microbiome makes it difficult to determine presence of unique lung taxa and therefore we labeled it undetermined pneumotype (pneumotypeUN). This pneumotype was present in 12/20 (60%) smokers and 5/9 (55%) never-smokers. PCoA of weighted UniFrac distances shows that BAL samples characterized as pneumotypeUN were in similar spatial location as background samples (panels B and D, Additional file 3: Figure S2). Similar results were found when unweighted UniFrac was used.
The second pneumotype had high relative abundance of supraglottic-characteristic taxa (SCT), such as Prevotella and Veillonella (Figure 4, panels A and C). We named this group pneumotypeSCT. PCoA of weighted UniFrac distances indicated that BAL samples characterized as pneumotypeSCT were in similar spatial locations as supraglottic samples (panels B and D, Additional file 3: Figure S2). PneumotypeSCT had four-fold higher total 16S rRNA gene concentrations compared with pneumotypeUN (183,312 (72,343-320,442)vs. 41,300 (22,375-96,000) copies/mL, P= 0.003). The 16S concentration in pneumotypeUN was not significantly different from background.
Compared with the never-smoker group, the lung microbiome of smokers had similar distribution along PCoA as well as α and β diversity (Additional file 4: Figure S3). Smoking did not have significant impact on pneumotype allocation.
Since 16S rRNA from the genus Streptococcus was found in background, BAL, and supraglottic samples, we evaluated the relative abundances of each Streptococcus OTU in all samples (Figure 5). Classification of Streptococcus OTU by Greengenes Blast showed that background samples were enriched with Streptococcus thermophilus (relative abundance 0.023 (0.010-0.036)), a known environmental contaminant . Similarly, pneumotypeUN had higher relative abundance of Streptococcus thermophilus than pneumotypeSCT (0.014 (0.005-0.044)vs. 0.003 (0.002-0.017), P= 0.04). Conversely, supraglottic samples were enriched with Streptococcus mitis (relative abundance 0.0520 (range, 0.0084-0.1625)), a well-known member of the oral microbiota and an opportunistic pathogen. PneumotypeSCT also was enriched with Streptococcus mitis (0.0057 (0.0024-0.0104)vs. 0.0441 (0.0266-0.0617), for pneumotypeUN and pneumotypeSCT respectively, P<0.0001).
The presence of bacterial DNA in the sterile saline used to perform the BAL led to the expected presence of background in each BAL sample. We therefore tested whether removal of sequences commonly found in the saline altered pneumotype definition. Subtraction of these ‘saline OTUs’ (≥1% relative abundance in saline) from BAL samples did not alter sample allocation in the unsupervised hierarchical clustering (data not shown). However, subtraction revealed high relative abundances of Sphingomonas, Tropheryma, Acidovorax, and Asticcacaulis in some pneumotypeUN samples, but not in pneumotypeSCT samples. UniFrac distances between pneumotypeUN and saline (0.72 ± 0.02) were lower than the distance between pneumotypeSCT and saline (0.78 ± 0.04, P= 0.002), also demonstrating a lower signal/noise ratio for samples in pneumotype 1.
Pneumotype and lung inflammation
Since lymphocytes and neutrophils are recruited into the alveolar space during inflammation, we used BAL cell counts as measures of lung inflammation. Compared with pneumotypeUN, pneumotypeSCT was associated with higher numbers of lymphocytes (11,130 (4,640-16,069)vs. 23,315 (10,256-47,706) cells/mL for pneumotypeUN and pneumotypeSCT respectively, P= 0.007) and neutrophils in BAL (3,333 (1,411-4,327)vs. 4,537 (3,211-8,259) cells/mL, P= 0.034, Figure 6A and B). We also evaluated eNO obtained prior to the bronchoscopy as an independent measure of lung inflammation (Figure 6C). Compared with pneumotypeUN, pneumotypeSCT was associated with higher levels of eNO (14.5 (9.0-23.3)vs. 27.5 (15.8-41.5) ppb, P= 0.022). Since smoking is known to decrease eNO we investigated whether the association between pneumotypeSCT and increased eNO was independent of smoking. First, eNO levels in never-smokers and smokers were similar (17.8 (10.9-36.0)vs. 20.0 (13.0-32.6), P= 0.71). Further, using a linear regression model that considered eNO as dependent variable and smoking and pneumotype as dichotomized predictors, we found that pneumotypeSCT was independently associated with eNO (B = 12.3 ppb, P= 0.028) while smoking was not significantly associated with levels of eNO (B = 0.5 ppb, P= 0.9). We then evaluated whether relative abundances of specific supraglottic-characteristic taxa were associated with lung inflammation. As shown in Figure 6D, E, and F, the relative abundance of Veillonella was positively correlated with both BAL inflammatory cells and eNO. Similarly, relative abundance of Prevotella correlated with BAL inflammatory cells and eNO (data not shown). These data indicate that a pneumotype defined by supraglottic-characteristic taxa is associated with inflammatory cells in the lung and with a non-invasive marker of airway inflammation.
This paper described distinct microbial communities that defined two pneumotypes in asymptomatic subjects. PneumotypeUN had low bacterial load and a taxonomic distribution similar (although not identical) to the background microbiome found in sterile saline. PneumotypeSCT had a higher 16S rRNA gene load and was enriched with supraglottic taxa such as Prevotella and Veillonella. PneumotypeSCT subjects had subclinical lung inflammation with increased BAL inflammatory cells and eNO. These data suggest an association between a distinct human microbiome and inflammation in the lung.
Supraglottic secretions have a 50-fold higher 16S rRNA gene concentration than BAL; therefore, carry-over of oropharyngeal secretions to the lower airways is a major challenge for bronchoscopy. Prior studies suggest that the lower airway microbiome have similar microbial composition to the upper airway, raising concerns for carry-over of oropharyngeal secretions to the lower airways [5, 6, 9, 10]. The route used for bronchoscopic approach might have significant impact since the microbiome of the oral and nasal cavity are different [40, 41]. We used the nasal route because of evidence that nasal brushes have low 16S rRNA gene concentrations . Our nasal bronchoscopic technique may have reduced contamination with oral secretions, enabling the uncovering of pneumotypeUN with low abundance of supraglottic-characteristic OTUs.
Three lines of evidence suggest that pneumotypeSCT was not due to carry-over of supraglottic taxa during bronchoscopy. First, had there been carry-over, the microbiome of first BAL would more closely resemble the supraglottic microbiome than the second BAL. In fact, UniFrac distances between supraglottic samples and the first BAL and second BAL were not different. Second, it is likely that the increased lymphocytes and neutrophils in the alveolar space of pneumotypeSCT subjects were almost certainly present prior to bronchoscopy since cell recruitment into the alveolar space probably requires more time than the duration of the bronchoscopy procedure. Third, exhaled NO was increased in pneumotypeSCT patients before the research bronchoscopy, demonstrating that this inflammatory phenotype existed prior to bronchoscopy.
In our study population, we found two distinct lung microbiomes. PneumotypeUN had low 16S rRNA concentration, little evidence of supraglottic-characteristic OTUs and a taxonomic composition similar to the normal saline used in the research bronchoscopy. These results differ from recently published series in which oral microbiota was present in most of the BAL samples [6, 11]. Technical differences between the studies may account for this discrepancy, since we aimed to minimize any potential carry-over with our bronchoscopic technique. The presence of microbial DNA in sterile saline made it difficult to ascertain which bacterial OTUs, if any, were present in the lung of pneumotypeUN subjects prior to the bronchoscopy. To reduce background noise, we subtracted saline-characteristic OTUs from the BAL. This revealed the presence of Sphingomonas, Tropheryma, Acidovorax, and Asticcacaulis in some pneumotypeUN but not in pneumotypeSCT samples. Our finding of Tropheryma in 2/16 patients in pneumotypeUN (relative abundance of 0.02 in the never-smoker and 0.11 in the smoker) is consistent with a recent report describing Tropheryma in the lower airways of 12% to 40% of normal subjects . Studies with larger cohorts using low contamination bronchoscopic techniques, deeper sequencing, and bioinformatic source tracking approaches are needed to better describe pneumotype(s) in subjects with low bacterial rRNA gene concentrations [42, 43].
Our clustering analysis allowed us to define pneumotypeSCT, which was present in 8/20 smokers (40%) and 4/9 never-smokers (45%). Prior investigations have attributed the observed increased abundance of oral-characteristic taxa in the lower airways to microaspiration or carry-over. This study extends prior investigations by showing that systematic carry-over does not occur and that the increased abundance of supraglottic-characteristic flora in pneumotypeSCT was associated with increased BAL neutrophils, lymphocytes, and eNO. Since pneumotypeSCT also has high relative abundance of Prevotella/Veillonella and high inflammatory cells in the lung, there may be a common link between microaspiration of Prevotella/Veillonella and inflammation. Importantly, the BAL cell differentials were normal in both never-smokers and smokers, suggesting that differences in lung inflammation between pneumotypeUN and pneumotypeSCT represent subclinical variation.
Data suggest that smoking alters the upper airway microbiome leading to enrichment with Veillonella[11, 41]. Persistent microorganisms contribute to pulmonary inflammation in current and former smokers with advanced COPD . Similar to recent studies , we found no significant difference in α or β diversity between never-smokers and smokers. It is possible that increasing the sequencing depth would reveal smoking or early COPD-specific microbes present at lower relative abundance. It also is possible that investigation of larger cohorts would reveal more substantial microbiome differences between these never-smokers and smokers. However, perturbation of the lung microbiome may occur only in more advanced disease due to progressive microaspiration and impaired bacterial clearance. Alteration of the lung microbiome also may follow treatments for COPD with antibiotics and/or inhaled steroids. The early nature of disease in our smoking cohort or the absence of recent drug treatment may account for the differences between our observations and the prior reports [7, 8, 21].
The association between specific pneumotype and airway inflammation was evaluated utilizing eNO based on prior observations describing increased eNO in pulmonary infections and following lipopolysaccharide administration [45–47]. Our results extend these observations by demonstrating that elevated eNO may occur in association with a specific pneumotype even in the absence of overt pneumonia. Increased eNO levels were associated with both pneumotypeSCT and enrichment with supraglottic taxa such as Veillonella suggesting that the presence of subclinical inflammatory changes may occur in response to pathogen-associated molecules. Although smoking may influence eNO levels, in our cohort, the association of increased eNO with pneumotypeSCT was independent of smoking history [48, 49]. Further investigation will be needed to determine whether a change in pneumotype leads to increase levels of eNO. The enhanced subclinical inflammation observed in pneumotypeSCT may warrant investigation into whether this pneumotype is associated with early physiologic markers of airway dysfunction in smokers at risk for COPD.
This study has several limitations. DNA molecular techniques cannot determine bacterial viability. Further, the association between inflammation and supraglottic-characteristic taxa in the lower airways does not imply that Prevotella and Veillonella in the lower airways cause inflammation. Additionally, the relatively small sample size, and the depth of sequencing may have prevented us from observing subtle but potentially important differences between smokers and never-smokers. Finally, our cross-sectional design does not allow us to measure the temporal stability of the lower airway microbiome. Evaluation of the impact of pneumotype resilience on lung inflammation is important to understand the potential role of the lung microbiome in the development of lung injury.
We observed that a pneumotype enriched with supraglottic-characteristic bacteria Prevotella and Veillonella was associated with higher inflammatory markers, consistent with the hypothesis that oral bacteria in the lung produce subclinical pulmonary inflammation. In this relatively healthy cohort, the association between microbiome and lung inflammation exists in asymptomatic subjects, both never-smokers and smokers. In advancing COPD, mucocilliary dysfunction worsens and the risk for microaspiration increases [50, 51]. This could introduce high relative abundance of oral flora into the lower airways leading to pneumotypeSCT and its associated increased inflammation. In addition, it is plausible that the progressive immune dysfunction of COPD changes the lung microbiome [52, 53]. Further understanding of the association between pneumotypes and pulmonary inflammation will be required to relate the above observations to the development of lung diseases and the potential impact of antibiotics.
Availability of supporting data
The dataset(s) supporting the results of this article available in the dbGaP repository (phs000633.v1.p1).
Chronic obstructive pulmonary disease
Exhaled nitric oxide
False discovery rate
Forced expiratory volume in 1 s
Functional residual capacity
Forced vital capacity
Operational taxonomic units
Principal coordinate analysis
Quantitative Insights into Microbial Ecology
Total lung capacity.
Ivanov II, Frutos Rde L, Manel N, Yoshinaga K, Rifkin DB, Sartor RB, Finlay BB, Littman DR: Specific microbiota direct the differentiation of IL-17-producing T-helper cells in the mucosa of the small intestine. Cell Host Microbe. 2008, 4: 337-349. 10.1016/j.chom.2008.09.009.
Chung H, Pamp SJ, Hill JA, Surana NK, Edelman SM, Troy EB, Reading NC, Villablanca EJ, Wang S, Mora JR, Umesaki Y, Mathis D, Benoist C, Relman DA, Kasper DL: Gut immune maturation depends on colonization with a host-specific microbiota. Cell. 2012, 149: 1578-1593. 10.1016/j.cell.2012.04.037.
Rogers GB, Daniels TW, Tuck A, Carroll MP, Connett GJ, David GJ, Bruce KD: Studying bacteria in respiratory specimens by using conventional and molecular microbiological approaches. BMC Pulm Med. 2009, 9: 14-10.1186/1471-2466-9-14.
Cabello H, Torres A, Celis R, El-Ebiary M, Puig de la Bellacasa J, Xaubet A, Gonzalez J, Agusti C, Soler N: Bacterial colonization of distal airways in healthy subjects and chronic lung disease: a bronchoscopic study. Eur Respir J. 1997, 10: 1137-1144. 10.1183/09031936.97.10051137.
Charlson ES, Bittinger K, Chen J, Diamond JM, Li H, Collman RG, Bushman FD: Assessing bacterial populations in the lung by replicate analysis of samples from the upper and lower respiratory tracts. PLoS One. 2012, 7: e42786-10.1371/journal.pone.0042786.
Charlson ES, Bittinger K, Haas AR, Fitzgerald AS, Frank I, Yadav A, Bushman FD, Collman RG: Topographical continuity of bacterial populations in the healthy human respiratory tract. Am J Respir Crit Care Med. 2011, 184: 957-963. 10.1164/rccm.201104-0655OC.
Erb-Downward JR, Thompson DL, Han MK, Freeman CM, McCloskey L, Schmidt LA, Young VB, Toews GB, Curtis JL, Sundaram B, Martinez FJ, Huffnagle GB: Analysis of the lung microbiome in the “healthy” smoker and in COPD. PLoS One. 2011, 6: e16384-10.1371/journal.pone.0016384.
Hilty M, Burke C, Pedro H, Cardenas P, Bush A, Bossley C, Davies J, Ervine A, Poulter L, Pachter L, Moffatt MF, Cookson WO: Disordered microbial communities in asthmatic airways. PLoS One. 2010, 5: e8578-10.1371/journal.pone.0008578.
Pragman AA, Kim HB, Reilly CS, Wendt C, Isaacson RE: The lung microbiome in moderate and severe chronic obstructive pulmonary disease. PLoS One. 2012, 7: e47305-10.1371/journal.pone.0047305.
Lozupone C, Cota-Gomez A, Palmer BE, Linderman DJ, Charlson ES, Sodergren E, Mitreva M, Abubucker S, Martin J, Yao G, Campbell TB, Flores SC, Ackerman G, Stombaugh J, Ursell L, Beck JM, Curtis JL, Young VB, Lynch SV, Huang L, Weinstock GM, Knox KS, Twigg H, Morris A, Ghedin E, Bushman FD, Collman RG, Knight R, Fontenot AP, Lung HIV Microbiome Project: Widespread colonization of the lung by Tropheryma whipplei in HIV infection. Am J Respir Crit Care Med. 2013, 187: 1110-1117. 10.1164/rccm.201211-2145OC.
Morris A, Beck JM, Schloss PD, Campbell TB, Crothers K, Curtis JL, Flores SC, Fontenot AP, Ghedin E, Huang L, Jablonski K, Kleerup E, Lynch SV, Sodergren E, Twigg H, Young VB, Bassis CM, Venkataraman A, Schmidt TM, Weinstock GM, Lung HIV Microbiome Project: Comparison of the respiratory microbiome in healthy nonsmokers and smokers. Am J Respir Crit Care Med. 2013, 187: 1067-1075. 10.1164/rccm.201210-1913OC.
Monso E, Rosell A, Bonet G, Manterola J, Cardona PJ, Ruiz J, Morera J: Risk factors for lower airway bacterial colonization in chronic bronchitis. Eur Respir J. 1999, 13: 338-342. 10.1034/j.1399-3003.1999.13b20.x.
Monso E, Ruiz J, Rosell A, Manterola J, Fiz J, Morera J, Ausina V: Bacterial infection in chronic obstructive pulmonary disease. A study of stable and exacerbated outpatients using the protected specimen brush. Am J Respir Crit Care Med. 1995, 152: 1316-1320. 10.1164/ajrccm.152.4.7551388.
Zalacain R, Sobradillo V, Amilibia J, Barron J, Achotegui V, Pijoan JI, Llorente JL: Predisposing factors to bacterial colonization in chronic obstructive pulmonary disease. Eur Respir J. 1999, 13: 343-348. 10.1034/j.1399-3003.1999.13b21.x.
Huang YJ, Lynch SV: The emerging relationship between the airway microbiota and chronic respiratory disease: clinical implications. Expert Rev Respir Med. 2011, 5: 809-821. 10.1586/ers.11.76.
Huang YJ, Nelson CE, Brodie EL, Desantis TZ, Baek MS, Liu J, Woyke T, Allgaier M, Bristow J, Wiener-Kronish JP, Sutherland ER, King TS, Icitovic N, Martin RJ, Calhoun WJ, Castro M, Denlinger LC, Dimango E, Kraft M, Peters SP, Wasserman SI, Wechsler ME, Boushey HA, Lynch SV, National Heart, Lung, and Blood Institute’s Asthma Clinical Research Network: Airway microbiota and bronchial hyperresponsiveness in patients with suboptimally controlled asthma. J Allergy Clin Immunol. 2011, 127: 372-381. 10.1016/j.jaci.2010.10.048. e371-373
Gleeson K, Eggli DF, Maxwell SL: Quantitative aspiration during sleep in normal subjects. Chest. 1997, 111: 1266-1272. 10.1378/chest.111.5.1266.
Cvejic L, Harding R, Churchward T, Turton A, Finlay P, Massey D, Bardin PG, Guy P: Laryngeal penetration and aspiration in individuals with stable COPD. Respirology. 2011, 16: 269-275. 10.1111/j.1440-1843.2010.01875.x.
Rascon-Aguilar IE, Pamer M, Wludyka P, Cury J, Coultas D, Lambiase LR, Nahman NS, Vega KJ: Role of gastroesophageal reflux symptoms in exacerbations of COPD. Ches. 2006, 130: 1096-1101. 10.1378/chest.130.4.1096.
Sze MA, Dimitriu PA, Hayashi S, Elliott WM, McDonough JE, Gosselink JV, Cooper J, Sin DD, Mohn WW, Hogg JC: The lung tissue microbiome in chronic obstructive pulmonary disease. Am J Respir Crit Care Med. 2012, 185: 1073-1080. 10.1164/rccm.201111-2075OC.
Huang YJ, Kim E, Cox MJ, Brodie EL, Brown R, Wiener-Kronish JP, Lynch SV: A persistent and diverse airway microbiota present during chronic obstructive pulmonary disease exacerbations. OMICS. 2010, 14: 9-59. 10.1089/omi.2009.0100.
Baines KJ, Simpson JL, Gibson PG: Innate immune responses are increased in chronic obstructive pulmonary disease. PLoS One. 2011, 6: e18426-10.1371/journal.pone.0018426.
Pace E, Ferraro M, Siena L, Melis M, Montalbano AM, Johnson M, Bonsignore MR, Bonsignore G, Gjomarkaj M: Cigarette smoke increases toll-like receptor 4 and modifies lipopolysaccharide-mediated responses in airway epithelial cells. Immunology. 2008, 124: 401-411. 10.1111/j.1365-2567.2007.02788.x.
Seemungal TA, Wilkinson TM, Hurst JR, Perera WR, Sapsford RJ, Wedzicha JA: Long-term erythromycin therapy is associated with decreased chronic obstructive pulmonary disease exacerbations. Am J Respir Crit Care Med. 2008, 178: 1139-1147. 10.1164/rccm.200801-145OC.
Albert RK, Connett J, Bailey WC, Casaburi R, Cooper JA, Criner GJ, Curtis JL, Dransfield MT, Han MK, Lazarus SC, Make B, Marchetti N, Martinez FJ, Madinger NE, McEvoy C, Niewoehner DE, Porsasz J, Price CS, Reilly J, Scanlon PD, Sciurba FC, Scharf SM, Washko GR, Woodruff PG, Anthonisen NR, COPD Clinical Research Network: Azithromycin for prevention of exacerbations of COPD. N Engl J Med. 2011, 365: 689-698. 10.1056/NEJMoa1104623.
Documents ATS: ATS/ERS recommendations for standardized procedures for the online and offline measurement of exhaled lower respiratory nitric oxide and nasal nitric oxide, 2005. Am J Respir Crit Care Med. 2005, 171: 912-930.
Segal LN, Methe BA, Nolan A, Hoshino Y, Rom WN, Dawson R, Bateman E, Weiden MD: HIV-1 and bacterial pneumonia in the era of antiretroviral therapy. Proc Am Thorac Soc. 2011, 8: 282-287. 10.1513/pats.201006-044WR.
Fierer N, Hamady M, Lauber CL, Knight R: The influence of sex, handedness, and washing on the diversity of hand surface bacteria. Proc Natl Acad Sci USA. 2008, 105: 17994-17999. 10.1073/pnas.0807920105.
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Pena AG, Goodrich JK, Gordon JI, Huttley GA, Kelley ST, Knights D, Koenig JE, Ley RE, Lozupone CA, McDonald D, Muegge BD, Pirrung M, Reeder J, Sevinsky JR, Turnbaugh PJ, Walters WA, Widmann J, Yatsunenko T, Zaneveld J, Knight R: QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010, 7: 335-336. 10.1038/nmeth.f.303.
Edgar RC: Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010, 26: 2460-2461. 10.1093/bioinformatics/btq461.
Wang Q, Garrity GM, Tiedje JM, Cole JR: Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007, 73: 5261-5267. 10.1128/AEM.00062-07.
Caporaso JG, Bittinger K, Bushman FD, DeSantis TZ, Andersen GL, Knight R: PyNAST: a flexible tool for aligning sequences to a template alignment. Bioinformatics. 2010, 26: 266-267. 10.1093/bioinformatics/btp636.
Price MN, Dehal PS, Arkin AP: FastTree 2–approximately maximum-likelihood trees for large alignments. PLoS One. 2010, 5: e9490-10.1371/journal.pone.0009490.
Lozupone C, Lladser ME, Knights D, Stombaugh J, Knight R: UniFrac: an effective distance metric for microbial community comparison. ISME J. 2011, 5: 169-172. 10.1038/ismej.2010.133.
Reiner A, Yekutieli D, Benjamini Y: Identifying differentially expressed genes using false discovery rate controlling procedures. Bioinformatics. 2003, 19: 368-375. 10.1093/bioinformatics/btf877.
Lozupone C, Knight R: UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol. 2005, 71: 8228-8235. 10.1128/AEM.71.12.8228-8235.2005.
Dray SaD AB: The ade4 package: implementing the duality diagram for ecologists. J Stat Softw. 2007, 22: 1-20.
Cailliez F: The analytical solution of the additive constant problem. Psychometrika. 1983, 48: 305-310. 10.1007/BF02294026.
DeSantis TZ, Hugenholtz P, Keller K, Brodie EL, Larsen N, Piceno YM, Phan R, Andersen GL: NAST: a multiple sequence alignment server for comparative analysis of 16S rRNA genes. Nucleic Acids Res. 2006, 34: W394-399. 10.1093/nar/gkl244.
Lemon KP, Klepac-Ceraj V, Schiffer HK, Brodie EL, Lynch SV, Kolter R: Comparative analyses of the bacterial microbiota of the human nostril and oropharynx. M Bio. 2010, 1: e00129-10.
Charlson ES, Chen J, Custers-Allen R, Bittinger K, Li H, Sinha R, Hwang J, Bushman FD, Collman RG: Disordered microbial communities in the upper respiratory tract of cigarette smokers. PLoS One. 2010, 5: e15216-10.1371/journal.pone.0015216.
Knights D, Kuczynski J, Charlson ES, Zaneveld J, Mozer MC, Collman RG, Bushman FD, Knight R, Kelley ST: Bayesian community-wide culture-independent microbial source tracking. Nat Methods. 2011, 8: 761-763. 10.1038/nmeth.1650.
Chen J, Bittinger K, Charlson ES, Hoffmann C, Lewis J, Wu GD, Collman RG, Bushman FD, Li H: Associating microbiome composition with environmental covariates using generalized UniFrac distances. Bioinformatics. 2012, 28: 2106-2113. 10.1093/bioinformatics/bts342.
Cosio MG, Saetta M, Agusti A: Immunologic aspects of chronic obstructive pulmonary disease. N Engl J Med. 2009, 360: 2445-2454. 10.1056/NEJMra0804752.
Adrie C, Monchi M, Dinh-Xuan AT, Dall'Ava-Santucci J, Dhainaut JF, Pinsky MR: Exhaled and nasal nitric oxide as a marker of pneumonia in ventilated patients. Am J Respir Crit Care Med. 2001, 163: 1143-1149. 10.1164/ajrccm.163.5.9906049.
Dressel H, de la Motte D, Reichert J, Ochmann U, Petru R, Angerer P, Holz O, Nowak D, Jorres RA: Exhaled nitric oxide: independent effects of atopy, smoking, respiratory tract infection, gender and height. Respir Med. 2008, 102: 962-969. 10.1016/j.rmed.2008.02.012.
Soop A, Sollevi A, Weitzberg E, Lundberg JO, Palm J, Albert J: Exhaled NO and plasma cGMP increase after endotoxin infusion in healthy volunteers. Eur Respir J. 2003, 21: 594-599. 10.1183/09031936.03.00008603.
Kharitonov SA, Robbins RA, Yates D, Keatings V, Barnes PJ: Acute and chronic effects of cigarette smoking on exhaled nitric oxide. Am J Respir Crit Care Med. 1995, 152: 609-612. 10.1164/ajrccm.152.2.7543345.
Malinovschi A, Janson C, Holmkvist T, Norback D, Merilainen P, Hogman M: Effect of smoking on exhaled nitric oxide and flow-independent nitric oxide exchange parameters. Eur Respir J. 2006, 28: 339-345. 10.1183/09031936.06.00113705.
Clunes LA, Davies CM, Coakley RD, Aleksandrov AA, Henderson AG, Zeman KL, Worthington EN, Gentzsch M, Kreda SM, Cholon D, Bennett WD, Riordan JR, Boucher RC, Tarran R: Cigarette smoke exposure induces CFTR internalization and insolubility, leading to airway surface liquid dehydration. FASEB J. 2012, 26: 533-545. 10.1096/fj.11-192377.
Simet SM, Sisson JH, Pavlik JA, Devasure JM, Boyer C, Liu X, Kawasaki S, Sharp JG, Rennard SI, Wyatt TA: Long-term cigarette smoke exposure in a mouse model of ciliated epithelial cell function. Am J Respir Cell Mol Biol. 2010, 43: 635-640. 10.1165/rcmb.2009-0297OC.
Hodge S, Hodge G, Scicchitano R, Reynolds PN, Holmes M: Alveolar macrophages from subjects with chronic obstructive pulmonary disease are deficient in their ability to phagocytose apoptotic airway epithelial cells. Immunol Cell Biol. 2003, 81: 289-296. 10.1046/j.1440-1711.2003.t01-1-01170.x.
Taylor AE, Finney-Hayward TK, Quint JK, Thomas CM, Tudhope SJ, Wedzicha JA, Barnes PJ, Donnelly LE: Defective macrophage phagocytosis of bacteria in COPD. Eur Respir J. 2010, 35: 1039-1047. 10.1183/09031936.00036709.
This research was supported in part by grants: K24 AI080298; CTSI Grant #UL1 TR000038; EDRN 5U01CA086137-13; Diane Belfer Program for Human Microbial Ecology; Michael Saperstein Scholarship Fund; R01DK090989; UH2 AR57506.
The authors declare that they have no competing interests.
Conception and design: LNS, MJB, MDW. Acquisition of data: LNS, RK, KIB, RMG. Analysis and interpretation of data: LNS, AVA, JCC, BW, HC, MJB, MDW. Drafting or revising of article: LNS, AVA, JCC, KIB, RMG, WNR, MJB, MDW. Final approval of the manuscript: LNS, AVA, JCC, MJB, MDW. All authors read and approved the final manuscript.
An erratum to this article is available at http://dx.doi.org/10.1186/2049-2618-2-21.
Electronic supplementary material
Additional file 1: Figure S1: Sampling scheme of background, supraglottic and BAL. Thirty-four background samples were obtained from either sterile saline used for BAL (n = 6) or saline passed through bronchoscope’s suctioning channel prior to procedure (n = 28). All bronchoscopies were performed nasally. In the 15 of the 29 studied subjects where supraglottic samples were obtained, we used a separate bronchoscope (Bronchoscope 1). This bronchoscope was passed without suctioning until vocal cords were visualized, at which time sample was obtained and then scope was withdrawn. The supraglottic sample was obtained by flushing 10 cc of normal saline through this bronchoscope. For all BAL samples, a separate bronchoscope (Bronchoscope 2) was passed without suctioning until wedged, at which time BAL was obtained. In the subset of 15 of the 29 studied subjects, BAL was differentially obtained from lingula and right middle lobe to evaluate for carry-over of supraglotic microorganisms. For rest of analysis, BAL obtained in the lingula and right middle lobe was pooled in all 29 subjects. (PDF 420 KB)
Additional file 2:Online supplement: Enrichment of lung microbiome with supraglottic taxa is associated with increased pulmonary inflammation.(PDF 146 KB)
Additional file 3: Figure S2: Evaluation of BAL microbiome compared with background and supraglottic microbiome in Principal Coordinate Analysis (PCoA). PCoA (x axis PC1 = 26.8% vs. y axis PC2 = 8.56%) based on weighted UniFrac distances for microbiome of background (black), BAL (never-smokers in green, smokers in blue), and supraglottic area (red). BAL samples had the highest variability as expressed by their distribution along PC1. Some BAL samples overlapped with background microbiome whereas others overlapped with supraglottic microbiome. (A) PCoA weighted by relative abundances of Propionibacterium (black boxes) showed higher relative abundances for this taxa in BAL samples that overlapped with background samples. (B,C) PCoA weighted for relative abundances of Prevotella and Veillonella (black boxes) showed higher relative abundances for these taxa among BAL samples that overlapped with supraglottic samples. (PDF 433 KB)
Additional file 4: Figure S3: Comparison between the lung microbiome of never-smokers and asymptomatic smokers. () PCoA based on weighted UniFrac distances for never smokers and smokers. PC1, PC2, and PC3 represent 43.7% of the variability on the data. Data shows complete overlapping of circle of inertia between smokers and never smokers. (B) α-diversity, expressed as richness, was similar in never-smokers and smokers. (C) β–diversity, based on weighted UniFrac distance for pairwise comparisons, among and between never-smoker and smoker subjects also was not significantly different (mean±SEM) between the groups. (PDF 278 KB)
About this article
Cite this article
Segal, L.N., Alekseyenko, A.V., Clemente, J.C. et al. Enrichment of lung microbiome with supraglottic taxa is associated with increased pulmonary inflammation. Microbiome 1, 19 (2013). https://doi.org/10.1186/2049-2618-1-19