Characterizing human lung tissue microbiota and its relationship to epidemiological and clinical features
© The Author(s). 2016
Received: 10 May 2016
Accepted: 7 July 2016
Published: 28 July 2016
The human lung tissue microbiota remains largely uncharacterized, although a number of studies based on airway samples suggest the existence of a viable human lung microbiota. Here we characterized the taxonomic and derived functional profiles of lung microbiota in 165 non-malignant lung tissue samples from cancer patients.
We show that the lung microbiota is distinct from the microbial communities in oral, nasal, stool, skin, and vagina, with Proteobacteria as the dominant phylum (60 %). Microbiota taxonomic alpha diversity increases with environmental exposures, such as air particulates, residence in low to high population density areas, and pack-years of tobacco smoking and decreases in subjects with history of chronic bronchitis. Genus Thermus is more abundant in tissue from advanced stage (IIIB, IV) patients, while Legionella is higher in patients who develop metastases. Moreover, the non-malignant lung tissues have higher microbiota alpha diversity than the paired tumors.
Our results provide insights into the human lung microbiota composition and function and their link to human lifestyle and clinical outcomes. Studies among subjects without lung cancer are needed to confirm our findings.
The human body harbors extraordinarily diverse communities of microbes (microbiota) that are increasingly thought to be crucial for human health. Recent studies have revealed intriguing correlations between specific patterns of human microbiota and various diseases, including autoimmune disorders, diabetes, obesity, and even psychiatric conditions [1–6].
The healthy human lung was traditionally considered sterile. However, since the first culture-independent report of microbiota in asthmatic airways , more than 30 studies using diverse molecular techniques have suggested that the healthy human lung is also home to bacteria (reviewed in [8, 9]). Because lung biopsy collection is not ethical in healthy human subjects, the study of lung microbiota has been mostly based on bronchoalveolar lavage (BAL), bronchoscopic brushing, or sputum samples. Reliance on these samples to determine lung microbiota is problematic due to contamination by the upper respiratory tract or oral microbiota . To date, only four studies on human lung tissue microbiota have been published. These studies have considerable limitations, including small sample size (n < 33) and use of samples mostly from patients with severe lung diseases, such as chronic obstructive pulmonary disease (COPD) or cystic fibrosis [10–13]. Therefore, the characteristics of lung tissue microbiota remain largely unknown.
In this study, we characterized the taxonomic and derived functional profiles of the microbiota in non-malignant lung tissue samples from 165 lung cancer patients and compared them with previously published profiles from other body sites (including oral cavity, nasal cavity, gut, skin, and vagina from the Human Microbiome Project ). We also evaluated the associations between features of non-malignant lung microbiota and epidemiological and clinical characteristics. Finally, we compared non-malignant with tumor lung microbiota.
Characteristics of the study participants
A total of 165 non-malignant lung tissue samples which generated at least 1000 sequence reads per sample (mean ± standard deviation (sd), 4091 ± 4167) were included in the analysis. Additional file 1: Table S1 describes the study population. The participants were mainly males (83 %) with a median age of 66.6 years; 53 % lived in the urban area of Milan (see Additional file 2: Figure S1 for a map of studied residential areas); most were smokers (51 % current and 43 % former) with a median of 45.3 pack-years and 43.5 years of smoking; 10 to 25 % self-reported a history of bronchitis, emphysema, and pneumonia; based on spirometry, 45 % subjects had a history of COPD. Most had tumor in the upper or lower lobe and 3 % had tumor in the principal bronchus; 38 % had squamous cell carcinoma and 59 % had adenocarcinoma; 92 % were diagnosed in stages IA, IB, IIA, IIB, and IIIA and only 8 % had more advanced cancer stages (IIIB, IV), as expected, since patients with later cancer stages are usually treated with systemic therapy instead of surgery. No patients had received chemotherapy, radiation therapy, or other treatments for lung cancer before surgery. The patients survived a median of 201.9 weeks after lung cancer diagnosis.
Taxonomic and functional profiles of the non-malignant lung tissue microbiota
We also examined the NIAID (National Institute of Allergy and Infectious Diseases) class A–C pathogen genera and opportunistic “pathogens” as defined by the PATRIC database . The potentially pathogenic genera Staphylococcus, Streptococcus, and Burkholderia were observed with low average relative abundances of around 2 %. Although Pseudomonas, which was frequently present in the lung tissue, is not included in the NIAID pathogen list, it could be pathogenic in immunosuppressed subjects. Other than these, “pathogens” were rarely observed in the non-malignant lung tissue (Additional file 1: Table S2).
Predicted functions based on Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt) analysis of 16S rRNA gene taxonomic data are shown in Fig. 1c for the most abundant modules. The functional profiles exhibited less variation across individuals than found for taxonomic profiles (compare Fig. 1a and 1c).
Analysis of negative controls
The number of reads in all negative controls (44–351 reads) was much lower than the number of reads in lung tissue samples (1551–12,340 reads) and the positive samples (2703–58,201) (Additional file 1: Table S3).
We found five operational taxonomic units (OTUs) that were shared across all negative controls. None of these five OTUs were present in the lung tissue samples.
At the phylum level, there was a strong difference in microbial profiles between the lung tissue samples and the 20 extraction negative controls (P = 0.001, performed on the Euclidean distance of phylum-level profiles by non-parametric permutation multivariate analysis of variance (MANOVA), Adonis test with 1000 permutations). The lung tissue samples also differed from the four PCR negative controls, although the P value was not statistically significant (P = 0.1) because of the small sample size (n = 4) of these negative controls.
There were 173 OTUs shared between the negative controls and the lung tissue samples. When we excluded the shared 173 OTUs from the original analysis, all results regarding lung tissue remained virtually unchanged.
Comparison of microbiotas from non-malignant lung tissue and other body sites
Figure 2c depicts the clustering of average phyla relative abundances by body site. The lung microbiota is distinct from that of other body sites in having a higher relative abundance of Proteobacteria, Thermi, and Cyanobacteria. Interestingly, Thermi had an average relative abundance of 8.8 % in lung but was rare in other body sites, with a highest average of 0.05 % in left antecubital fossa (Additional file 2: Figure S2). The Taq DNA polymerase used in this study was produced from Escherichia coli, not from Thermus aquaticus; therefore, it is unlikely that the Thermi observed in our samples were due to contamination. To confirm this, we re-sequenced five samples that originally included Thermi and five samples that originally had no Thermi. Five out of five positive samples remained positive in the re-sequencing data and five out of five negative samples remained negative in the re-sequencing data (Additional file 1: Table S4). These replication data argue strongly against laboratory contamination during PCR amplification or sequencing as the source of Thermus (Thermi) in the lung specimens, in which Thermus (Thermi) was one of the most abundant genera. In the repeated sequencing, the copy numbers of Thermi in the five positive cases were lower than in the original samples. This may be due to batch effects in the PCR amplification or because the DNA amount remaining for the replication assay was lower than the amount used for the original analysis.
Analogous clustering of predicted functions (Fig. 2d) shows that the lung microbiota had high relative abundance of the KEGG modules amino acid metabolism, xenobiotic biodegradation and metabolism, and lipid metabolism. Functional profiles and taxonomic profiles were correlated (Additional file 2: Figure S3).
Demographic and clinical associations of non-malignant lung tissue microbiota
Analysis of beta diversity by residential area or PM10 MANOVA adjusted for history of bronchitis and tumor stage) showed a statistically significant association based on unweighted UniFrac distance (P = 0.009 and 0.006, respectively) but not on weighted UniFrac distance (P = 0.14 and 0.16, respectively). However, when both residential area and PM10 were included in the same model, PM10 (P = 0.003), but not residential area (P = 0.17), remained statistically significant.
We observed no statistically significant differences among microbiota from various anatomical locations in the lung (Additional file 1: Table S5). Samples from the principal bronchus had non-significantly higher taxonomic alpha diversity than samples from lung lobes (observed species (mean ± sd), 116.6 ± 53.6 in principal bronchus versus 84.0 ± 24.0 in lobes), but this test was based on only five principal bronchus samples. Alpha diversity was similar in different lobes of the lung (Additional file 1: Table S5).
We observed a significant positive association of pack-years of cigarette smoking with taxonomic alpha diversity (Shannon index, Ptrend = 0.04 and observed species Ptrend = 0.04). No other significant association for lung microbiota measures was observed by smoking status, years of smoking, or cigarettes per day (Additional file 1: Table S6).
Compared with patients with no history of previous lung diseases, patients with a history of emphysema, COPD, and pneumonia had similar levels of taxonomic alpha diversity (Additional file 1: Table S7), but patients with a history of bronchitis had significantly decreased alpha diversity (observed species, P = 0.05; PD_whole_tree, P = 0.02). No difference in beta diversity and relative abundance of any taxa was found between patients with and without previous lung diseases (data not shown). The dominant phylum, Proteobacteria, did not differ by previous disease status, including COPD (mean ± sd, 0.60 ± 0.22 and 0.62 ± 0.25 for patients with and without COPD, respectively). We also examined the association of spirometry-based lung function measures, including forced vital capacity, forced expiratory volume in 1 s, peak expiratory flow, and the ratio of forced vital capacity and forced expiratory volume in 1 s, with lung tissue microbial features, including alpha and beta diversity and taxa relative abundance, and found no association (data not shown).
Microbiota from non-malignant lung tissue had significantly increased PD_whole_tree, but not Shannon’s index in late stages IIIB and IV (Fig. 3c). Beta diversity (adjusted MANOVA analysis) was also significantly associated with tumor stage (P = 0.001 for both unweighted and weighted UniFrac). The genus Thermus (Thermi) had significantly higher abundance in stages IIIB and IV (Fig. 3c). With respect to predicted function, microbiota significantly differed by tumor stage in predicted KEGG modules and pathways (Additional file 2: Figure S4). Specifically, compared with patients in stages IA to IIIA, microbiota in stage IV patients had increased relative abundance for the excretory system module and the amino acid metabolism, aldosterone regulated sodium reabsorption, and amoebiasis pathways. Moreover, microbiota in patients with both stage IIIB and IV had reduced abundance for signal transduction (Additional file 2: Figure S4).
We observed no difference in taxonomic alpha diversity or beta diversity by metastasis status after diagnosis. However, the non-malignant samples from patients who developed metastases, compared with those without metastases, had significantly increased relative abundance of Legionella (Proteobacteria) (mean ± sd, 0.003 ± 0.008 versus 0.001 ± 0.004, P(Bonferroni) = 0.01).
Difference between lung tumor and non-malignant tissue microbiotas
In the largest study of human non-malignant lung tissue to date, we describe the taxonomic and functional profiles of lung tissue microbiota. The lung tissue microbiota was clearly distinct from the microbiotas reported at other body sites (oral cavity, nasal cavity, gut, skin, and vagina). Moreover, it showed increased alpha diversity with environmental exposures such as air particulates, residence in high population density areas, and pack-years of tobacco smoking. Microbiota also varied by clinical endpoints, with increased alpha diversity in non-malignant lung tissue from advanced stages of cancer and decreased alpha diversity in lung affected by chronic bronchitis. Microbiota alpha diversity also significantly differed between non-malignant and tumor lung tissue.
Most previous studies on the airway microbiota were based on BAL, bronchoscopic brushing, or sputum samples [7, 10, 17–26]. A common concern of these samples is that they may be contaminated by the upper respiratory tract and oral microbiota . In our study, samples were surgically resected from lung tissue distant from the tumors and with no evidence of tumor nuclei. We found that five Proteobacteria genera had high relative abundance in lung tissue with and without COPD or other lung diseases and each genus was found in 80 % of samples. In contrast, previous studies of BAL showed high abundance of genera Prevotella (Bacteroidetes), Streptococcus (Firmicutes) and Veillonella (Firmicutes) (Additional file 1: Table S8), which are commonly found in the oral cavity (Additional file 1: Table S9). Reassuringly, the only four previous, very small studies of lung tissue also indicated members of the phylum Proteobacteria as predominant (Additional file 1: Table S6) [10–13]. Notably, we found that Thermi, although present in only 27 % of subjects, had high relative abundance in samples from these subjects (mean ± sd, 32 ± 20 %). To our knowledge, only one BAL-based study reported the presence of high relative abundance of Thermi (Thermus, ~96 %) in one healthy subject . Moreover, our data clearly show that the lung microbiota is distinct from digestive tract microbiota of healthy subjects, even at high taxonomic (phylum) and functional (KEGG module) level (Fig. 2a, b). Taken together these data suggest that the lung microbiota is unique.
A concern is that our lung tissue assays might be contaminated during DNA extraction, PCR amplification, or sequencing. Analyses of 24 negative controls and re-amplification and re-sequencing of DNA from ten previously analyzed lung specimens argue strongly against distortion of our results by contamination. Although Thermi has high resistance to environmental hazards  (Thermus acquaticus can survive at 50 to 80 °C and is the source of thermostable Taq DNA polymerase, which is commonly used for DNA amplification (https://en.wikipedia.org/wiki/Thermus_aquaticus)), we found no evidence for Thermi contamination. Repeat amplification and sequencing of five initially Thermi-negative DNA samples yielded five negative results and five initially Thermi-positive DNA samples yielded five repeat positive results. More generally, when all the 173 OTUs that were common to 24 negative control samples and to any lung cancer sample were excluded, the lung cancer analyses were virtually unchanged.
We found a positive association between microbiota alpha diversity in the lung and subjects’ residential area (from low to high density population), with evidence that this association reflected exposure to air pollution (PM10). We did not find significant differences in the lung tissue microbiota by smoking status, smoking intensity, or lifetime smoking duration, probably because there was no great variability across study subjects since most of them were heavy smokers. However, we did find significantly higher alpha diversity with increased pack-years of cigarette smoking. Together, these findings suggest that the lung microbiota may be altered by cumulative exposure to tobacco smoke and other air pollutants or life style conditions linked to high density population. Previously, analysis of 16 sputum samples revealed higher alpha diversity in samples from women in China who used smoky coal for cooking and heating compared with those using smokeless coal. Also, increased diversity and altered abundances of certain taxa in smokers’ versus nonsmokers’ subgingival samples have been shown . Chronic inhalation of dust or tobacco-related particles could allow increased diversity by impeding the dispersion and clearance of microbes from the bronchopulmonary system. Alternatively, particulates in the air could function as vectors for inhalation of microbes, as suggested by a study in which the dust from households with a dog or a cat had higher microbial diversity compared with the dust from households with no furred pets .
Long lasting and repetitive irritation of inhaled substances such as tobacco smoke, dust and silica may promote the development of bronchitis. Early clinical features of bronchitis include hyper-secretion of mucus and hypertrophy of sub-mucosal glands, eventually leading to chronic airway obstruction and possibly secondary growth of specific bacteria. This could be related to the lower alpha diversity we identified in subjects affected by chronic bronchitis. Clearly, larger clinical studies and investigations in model systems are warranted to further understand the viability of the microbiota, its role in chronic diseases, and its potential use for prevention and treatment strategies.
Despite painstaking exclusion of samples adjacent to tumors, we observed alteration of the lung microbiota in non-malignant tissue samples from subjects with advanced tumor stages (IIIB and IV). Specifically, these samples had higher phylogenetic diversity, high relative abundance of Thermus (Thermi), and increased/decreased abundance of several functional modules compared with samples from patients with earlier stages of lung cancer. Moreover, subjects who developed metastases had high relative abundance of Legionella (Proteobacteria). These data suggest that Thermus and Legionella might play a role in tumor progression, partially through the different microbiota functions, e.g., reduced signal transduction, increased excretory systems, amino acid metabolism, aldosterone-regulated sodium reabsorption, or amoebiasis pathways. Alternatively, tumor progression could affect the microenvironment and microbiota of a larger surrounding area. Given our collection methods and careful histological review, it is unlikely that the non-malignant tissue samples from advanced-stage subjects were contaminated by their tumor microbiota. In fact, the tumor microbiota showed low phylogenetic diversity that was unlike the corresponding non-malignant tissue. Moreover, while the microbiota differed in the tumor samples between subjects with adenocarcinoma and squamous cell carcinoma, the corresponding microbiota in these subjects’ non-malignant samples did not differ. This suggests that the microbiota from non-malignant lung tissue samples are different from that of the tumor lung tissue samples, as has been previously shown for tumor/non-malignant samples from colorectal cancer patients . It will be important to explore whether the microbiota in non-malignant lung tissue from advanced disease stages, or in tumor tissue, plays a role in tumor progression or is just a passive byproduct of tumor progression.
This study includes noteworthy strengths and limitations. It is the largest study of the non-malignant lung tissue microbiota to date. Moreover, we used uniform surgical procedures performed under sterile conditions for obtaining the surgical samples, which were frozen immediately. Also, we examined a comprehensive list of detailed and validated epidemiological and clinical variables in relation to features of the lung microbiota. Furthermore, we performed rigorous analysis to take into account correlations across subjects and paired samples and sequenced negative controls to exclude the possibility of contamination. One important caveat is that we compared our lung data with that of healthy individuals enrolled in the HMP. The two populations are very different, in age and health status, with young and extremely healthy individuals in the HMP study and older (average age 67 years) lung cancer patients in this study. Further, DNA extraction techniques, 16S rRNA gene primers, and sequencing platforms were all different and have been previously shown to potentially introduce some biases . However, we minimized the effects of these differences by restricting our comparisons to the highest and least variable taxonomic level (phylum level) and functional entity (KEGG module). In addition, a meta-analysis of microbiota studies revealed that differences in microbial populations across body sites are significantly larger than those driven by the experimental protocols, age, geography, and other population characteristics . Another limitation is that we could not study the effect of antibiotic use on the lung microbiota because most patients were treated with antibiotics at the time of surgery. In addition, as in most microbiota studies, we do not know whether the DNA we studied corresponded to living or dead bacteria and further studies are needed to address this issue. Furthermore, because we used a stringent threshold (1000 reads) to obtain reliable estimates of microbiota relative abundance, we had to exclude ~23 and ~45 % of samples from non-malignant and tumor sites, respectively. If we had used a less stringent threshold, e.g., 500 reads (which is still larger than the number of reads identified in most negative controls), we would have excluded only 13 and 25 % of samples, respectively. We opted to use a more stringent threshold since this is the first large study of microbiota in human lung tissue and it is important to report data with enough reads to characterize the lung bacterial community accurately. Finally, although we found no association between microbiota features and COPD or spirometry-based lung function, we cannot exclude that the lung cancer patients had abnormalities in their lungs that could affect our results. Our findings in lung cancer patients, although based on non-malignant tissues, may not be completely applicable to healthy subjects.
In the largest study of non-malignant lung tissue to date, we show that the lung microbiota has distinct features that differ from those of the oral cavity and other body sites and is dominated by Proteobacteria (60 %). The lung microbiota is affected by exposure to air pollution and tobacco smoking and is different in subjects with chronic bronchitis or advanced tumors. The genus Thermus is more abundant in tissue from advanced stage patients, while Legionella is higher in patients who develop metastases. Further studies in lung tissue and animal model systems are necessary to investigate the role of microbiota in the development of lung diseases and whether it can be exploited for treatment purposes.
Subject characteristics and epidemiological and clinical data collection
The study is nested in the Environment and Genetics in Lung cancer Etiology study (EAGLE), which was described in detail previously . In brief, EAGLE is an integrated population-based study of lung cancer with the aim to capture the major risk factors and genetic basis of lung cancer. The study was conducted in the Lombardy region of Italy and the catchment area included five cities (Milan, Monza, Brescia, Pavia, and Varese) and their surrounding municipalities (see Additional file 2: Figure S1 for a map). The lung cancer cases were enrolled from 13 hospitals that covered approximately 80 % of lung cancer cases from the catchment area. The epidemiological data were collected by a computer assisted personal interview and a self-administered questionnaire at the time of lung cancer diagnosis. The clinical data were collected by physicians from the clinical charts and hospital discharge records. The average annual atmospheric concentration (μg/m3) of particulate matter of 10 micrometers in diameter (PM10) in the subjects’ cities of residence during the year of their enrollment was estimated by combining land-use regression data with aerosol optical depth data from the MODIS (Moderate Resolution Imaging Spectroradiometer) instrument onboard the National Aeronautics and Space Administration’s Terra satellite .
Lung tissue samples were snap-frozen in liquid nitrogen within 20 minutes of surgical resection. Surgeons and pathologists were together in the surgery room at the time of resection and sample collection to ensure correct sampling of tissue from the tumor, the area adjacent to the tumor, and an additional area distant from the tumor (~1–5 cm), without adversely affecting the participant. The precise site of tissue sampling was indicated on a lung drawing and the pathologists classified the samples as tumor, adjacent lung tissue, and distant non-involved lung tissue. For the current study, we used lung tissue sampled from an area distant from the tumor (defined here as “non-malignant lung tissue”) to reduce the potential for local cancer field effects. For each subject, usually more than one non-malignant lung tissue sample was collected and at least one sample was examined by a pathologist to confirm the absence of tumor nuclei. All the tools and materials in contact with the lung tissues were sterile. Based on sample availability, we selected 233 non-malignant and 56 tumor samples. Results are based on 165 non-malignant and 31 tumor samples after quality control-determined exclusions.
16S rRNA gene sequence analysis
Fresh frozen lung tissue samples remained frozen while approximately 30 mg was subsampled for DNA extraction into pre-chilled 2.0 ml microcentrifuge tubes. Lysates for DNA extraction were generated by incubating 30 mg of tissue in 1 ml of 0.2 mg/ml Proteinase K (Ambion) in DNA lysis buffer (10 mM Tris-Cl (pH 8.0), 0.1 M EDTA (pH 8.0), and 0.5 % (w/v) SDS) for 24 h at 56 °C with shaking at 850 rpm in Thermomixer R (Eppendorf). DNA was extracted from the generated lysate using the QIAamp DNA Blood Maxi Kit (Qiagen) according to the manufacturer’s recommendation. The V3–V4 regions of the 16S rRNA gene were amplified and sequenced on an Illumina MiSeq instrument using the 300 paired-end protocol at the Institute of Genome Sciences, Genomic Resource Center, University of Maryland School of Medicine as described previously . We included two positive controls (one fecal and one vaginal sample) to examine the performance of the sequence run, 20 negative controls to examine the potential contamination by DNA extraction and PCR reagents, and four negative controls to exclude contamination during the PCR amplification process. We showed that our results are not affected by potential contamination.
Sequence reads were processed to remove low quality reads, specifically reads with average quality less than 20 over a 30-bp window based on the Phred algorithm. These were trimmed before the first base of the window and re-evaluated for length. Also removed were paired reads that had at least one of the reads with length less than 75 % of its original length, reads with less than 60 % similarity to Greengenes reference version 13_8 , and chimera reads (identified using UCHIME ). The remaining reads were clustered into OTUs at 97 % identity using the command pick_open_reference_otus.py in the software package Quantitative Insights into Microbial Ecology (QIIME 1.8.0) . The default parameters were used except method of usearch61 and percent_subsample of 0.1. OTUs with only one read or in only one sample were excluded.
Taxonomic alpha diversity was estimated as the number of 97 % identical OTUs (Observed_species), Shannon’s Index (using information of the relative abundance of observed species)  and phylogenetic diversity whole tree (PD_whole_tree, using information on both the relative abundance and phylogenetic tree of observed species)  by averaging over 20 rarefied tables (1000 reads/sample). Taxonomic beta diversity was measured as unweighted (presence/absence of observed species) and weighted UniFrac distance (also using information on the relative abundance of observed species) based on the OTU table . Relative abundance of taxa was calculated from unrarefied OTU table.
We downloaded the V3–V5 16S rRNA gene sequence data from the Human Microbiome Project (HMP, phase 1; http://hmpdacc.org/) and processed for comparison . The HMP data include 138–1623 samples with at least 1000 sequence reads per sample from each studied site (including oral cavity, nasal cavity, skin, stool, and vagina). Euclidean distance and Bray–Curtis distance were calculated based on phylum-level relative abundance for comparison between lung tissue and other body sites.
Functional prediction from 16S rRNA gene sequence data
PICRUSt 1.0.0 was used to predict the function of the microbiota from the 16S rRNA gene sequence taxonomic data for both the HMP dataset and this study using the KEGG database as reference [42, 43]. PICRUSt requires the use of Greengenes reference version 13_5 . Therefore, we reprocessed the sequence data in QIIME as previously but using Greengenes reference version gg_13_5 .
Euclidean distance and Bray–Curtis distance were calculated using the rarified KEGG orthologs (KO) table (430,000 predicted reads per sample) for comparison among lung samples and the KEGG module table for comparison with other body sites. Relative abundance of modules/pathways was calculated from unrarefied KO table.
All statistical analyses were performed in the R software (R Foundation for Statistical Computing, Vienna, Austria; http://www.r-project.org/). In boxplots, the black central lines represent the median and box edges the first and third quartiles. Wilcoxon rank-sum and Kruskal–Wallis tests were used for differences between categories, and Spearman correlation test was used for association of continuous variables. For the variables that showed significant associations with microbiome features, we used multiple linear regression models with microbiota measurements as the dependent variable to test the association while adjusting for other covariates (residential area, history of bronchitis, and tumor stage). For comparisons between non-malignant and tumor samples, the Wilcoxon signed-rank test was used for paired samples. For comparisons that included both paired and unpaired samples, individuals were stratified into three categories, those with tumor tissue only, those with non-malignant tissue only, and those with both types of tissue. A bootstrap was performed by resampling individuals with replacement within these strata to estimate the variances of mean differences. Jackknife analyses likewise removed individuals one at a time from these strata to account for correlations among means. Only the taxa (phylum, class, order, family, or genus) with relative abundance greater than 0.001 in at least 10 % of the samples were included in the analysis. Unless otherwise indicated, P values were Bonferroni-adjusted for multiple comparisons (R command P.adjust).
MANOVA in R (Adonis method in R Package Vegan ) was used to examine the association between beta diversity and individual epidemiological and clinical variables, adjusting for the variables showing significant associations with microbiota (including subjects’ residential area, history of bronchitis, and tumor stage). For some comparisons of beta-diversity, the jackknife procedure was used to compute variances that allowed for correlations within and between subjects.
BAL, bronchoalveolar lavage; COPD, chronic obstructive pulmonary disease; EAGLE, Environment and Genetics in Lung cancer Etiology study; HMP, Human Microbiome Project; KO, KEGG orthologs; MANOVA, Non-parametric permutation multivariate analysis of variance; NIAID, National Institute of Allergy and Infectious Diseases; OTU, operational taxonomic unit; PICRUSt, Phylogenetic Investigation of Communities by Reconstruction of Unobserved States; PM, particulate matter; sd, standard deviation
This work was supported by the Intramural Research Program of the Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, DHHS, Bethesda, MD, USA. We are deeply indebted to the EAGLE participants and the study collaborators (listed in http://dceg.cancer.gov/eagle). The principal investigators had full access to all the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis. We thank B. Ma from Ravel lab for de-noising and delivering the sequence data.
This study was supported by the Intramural Research Program of NIH, NCI, Division of Cancer Epidemiology and Genetics.
Availability of data and materials
The sequence data from this study have been submitted to NCBI BioProject (http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/bioproject) under accession number PRJNA303190.
GY and MTL designed the study. JR and MH performed the laboratory experiments. GY analyzed the data. MG provided statistical advice and supervised the statistical analyses. DC, MC, and AP collected and analyzed pollution data. GY, MG, and MTL wrote the manuscript. All authors contributed to the data interpretation and manuscript revision. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Ethics approval and consent to participate
The EAGLE study was approved by the institutional review board of each participating hospital and university in Italy and by the National Cancer Institute, Bethesda, MD, USA (IRB n. 01-C-N211). All participants provided written informed consent. All experimental methods used for this study comply with the Helsinki Declaration.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Round JL, Mazmanian SK. The gut microbiota shapes intestinal immune responses during health and disease. Nat Rev Immunol. 2009;9:313–23.View ArticlePubMedPubMed CentralGoogle Scholar
- Scher JU, Abramson SB. The microbiome and rheumatoid arthritis. Nat Rev Rheumatol. 2011;7:569–78.PubMedPubMed CentralGoogle Scholar
- Qin J, Li Y, Cai Z, Li S, Zhu J, Zhang F, Liang S, Zhang W, Guan Y, Shen D, et al. A metagenome-wide association study of gut microbiota in type 2 diabetes. Nature. 2012;490:55–60.View ArticlePubMedGoogle Scholar
- Ley RE, Turnbaugh PJ, Klein S, Gordon JI. Microbial ecology: human gut microbes associated with obesity. Nature. 2006;444:1022–3.View ArticlePubMedGoogle Scholar
- Turnbaugh PJ, Hamady M, Yatsunenko T, Cantarel BL, Duncan A, Ley RE, Sogin ML, Jones WJ, Roe BA, Affourtit JP, et al. A core gut microbiome in obese and lean twins. Nature. 2009;457:480–4.View ArticlePubMedGoogle Scholar
- Foster J, Neufeld KA. Gut-brain axis: how the microbiome influences anxiety and depression. Int J Neuropsychoph. 2014;17:27.Google Scholar
- Hilty M, Burke C, Pedro H, Cardenas P, Bush A, Bossley C, Davies J, Ervine A, Poulter L, Pachter L, et al. Disordered microbial communities in asthmatic airways. PLoS One. 2010;5:e8578.View ArticlePubMedPubMed CentralGoogle Scholar
- Berger G, Wunderink RG. Lung microbiota: genuine or artifact? Isr Med Assoc J. 2013;15:731–3.PubMedGoogle Scholar
- Dickson RP, Huffnagle GB. The lung microbiome: new principles for respiratory bacteriology in health and disease. PLoS Pathog. 2015;11:e1004923.View ArticlePubMedPubMed CentralGoogle Scholar
- Erb-Downward JR, Thompson DL, Han MK, Freeman CM, McCloskey L, Schmidt LA, Young VB, Toews GB, Curtis JL, Sundaram B, et al. Analysis of the lung microbiome in the “healthy” smoker and in COPD. PLoS One. 2011;6:e16384.View ArticlePubMedPubMed CentralGoogle Scholar
- Sze MA, Dimitriu PA, Suzuki M, McDonough JE, Campbell JD, Brothers JF, Erb-Downward JR, Huffnagle GB, Hayashi S, Elliott WM, et al. Host response to the lung microbiome in chronic obstructive pulmonary disease. Am J Respir Crit Care Med. 2015;192:438–45.View ArticlePubMedGoogle Scholar
- 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–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Willner D, Haynes MR, Furlan M, Schmieder R, Lim YW, Rainey PB, Rohwer F, Conrad D. Spatial distribution of microbial communities in the cystic fibrosis lung. ISME J. 2012;6:471–4.View ArticlePubMedGoogle Scholar
- Human Microbiome Project Consortium. Structure, function and diversity of the healthy human microbiome. Nature. 2012;486:207–14.View ArticleGoogle Scholar
- Wattam AR, Abraham D, Dalay O, Disz TL, Driscoll T, Gabbard JL, Gillespie JJ, Gough R, Hix D, Kenyon R, et al. PATRIC, the bacterial bioinformatics database and analysis resource. Nucleic Acids Res. 2014;42:D581–91.View ArticlePubMedGoogle Scholar
- Cesaroni G, Porta D, Badaloni C, Stafoggia M, Eeftens M, Meliefste K, Forastiere F. Nitrogen dioxide levels estimated from land use regression models several years apart and association with mortality in a large cohort study. Environ Health. 2012;11:48.View ArticlePubMedPubMed CentralGoogle Scholar
- Bassis CM, Erb-Downward JR, Dickson RP, Freeman CM, Schmidt TM, Young VB, Beck JM, Curtis JL, Huffnagle GB. Analysis of the upper respiratory tract microbiotas as the source of the lung and gastric microbiotas in healthy individuals. MBio. 2015;6:e00037.View ArticlePubMedPubMed CentralGoogle Scholar
- Borewicz K, Pragman AA, Kim HB, Hertz M, Wendt C, Isaacson RE. Longitudinal analysis of the lung microbiome in lung transplantation. FEMS Microbiol Lett. 2013;339:57–65.View ArticlePubMedGoogle Scholar
- 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–63.View ArticlePubMedPubMed CentralGoogle Scholar
- Dickson RP, Erb-Downward JR, Freeman CM, McCloskey L, Beck JM, Huffnagle GB, Curtis JL. Spatial variation in the healthy human lung microbiome and the adapted island model of lung biogeography. Ann Am Thorac Soc. 2015;12:821–30.View ArticlePubMedPubMed CentralGoogle Scholar
- Dickson RP, Martinez FJ, Huffnagle GB. The role of the microbiome in exacerbations of chronic lung diseases. Lancet. 2014;384:691–702.View ArticlePubMedPubMed CentralGoogle Scholar
- Marsland BJ, Yadava K, Nicod LP. The airway microbiome and disease. Chest. 2013;144:632–7.View ArticlePubMedGoogle Scholar
- Morris A, Beck JM, Schloss PD, Campbell TB, Crothers K, Curtis JL, Flores SC, Fontenot AP, Ghedin E, Huang L, et al. Comparison of the respiratory microbiome in healthy nonsmokers and smokers. Am J Respir Crit Care Med. 2013;187:1067–75.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Segal LN, Alekseyenko AV, Clemente JC, Kulkarni R, Wu B, Chen H, Berger KI, Goldring RM, Rom WN, Blaser MJ, Weiden MD. Enrichment of lung microbiome with supraglottic taxa is associated with increased pulmonary inflammation. Microbiome. 2013;1:19.View ArticlePubMedPubMed CentralGoogle Scholar
- Willner DL, Hugenholtz P, Yerkovich ST, Tan ME, Daly JN, Lachner N, Hopkins PM, Chambers DC. Reestablishment of recipient-associated microbiota in the lung allograft is linked to reduced risk of bronchiolitis obliterans syndrome. Am J Respir Crit Care Med. 2013;187:640–7.View ArticlePubMedGoogle Scholar
- Griffiths E, Gupta RS. Identification of signature proteins that are distinctive of the Deinococcus-Thermus phylum. Int Microbiol. 2007;10:201–8.PubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Fujimura KE, Johnson CC, Ownby DR, Cox MJ, Brodie EL, Havstad SL, Zoratti EM, Woodcroft KJ, Bobbitt KR, Wegienka G, et al. Man’s best friend? The effect of pet ownership on house dust microbial communities. J Allergy Clin Immunol. 2010; 126:410–412, 412 e411–413.Google Scholar
- Burns MB, Lynch J, Starr TK, Knights D, Blekhman R. Virulence genes are a signature of the microbiome in the colorectal tumor microenvironment. Genome Med. 2015;7:55.View ArticlePubMedPubMed CentralGoogle Scholar
- Lazarevic V, Gaia N, Girard M, Francois P, Schrenzel J. Comparison of DNA extraction methods in analysis of salivary bacterial communities. PLoS One. 2013;8:e67699.View ArticlePubMedPubMed CentralGoogle Scholar
- Lozupone C, Stombaugh J, Gonzalez A, Ackermann G, Wendel D, Vazquez-Baeza Y, Jansson JK, Gordon JI, Knight R. Meta-analyses of studies of the human microbiota. Genome Res. 2013;23:1704–14.View ArticlePubMedPubMed CentralGoogle Scholar
- Landi MT, Consonni D, Rotunno M, Bergen AW, Goldstein AM, Lubin JH, Goldin L, Alavanja M, Morgan G, Subar AF, et al. Environment And Genetics in Lung cancer Etiology (EAGLE) study: an integrative population-based case-control study of lung cancer. BMC Public Health. 2008;8:203.View ArticlePubMedPubMed CentralGoogle Scholar
- Nordio F, Kloog I, Coull BA, Chudnovsky A, Grillo P, Bertazzi PA, Baccarelli AA, Schwartz J. Estimating spatio-temporal resolved PM10 aerosol mass concentrations using MODIS satellite data and land use regression over Lombardy, Italy. Atmos Environ. 2013;74:227–36.View ArticleGoogle Scholar
- Fadrosh DW, Ma B, Gajer P, Sengamalay N, Ott S, Brotman RM, Ravel J. An improved dual-indexing approach for multiplexed 16S rRNA gene sequencing on the Illumina MiSeq platform. Microbiome. 2014;2:6.View ArticlePubMedPubMed CentralGoogle Scholar
- DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, Huber T, Dalevi D, Hu P, Andersen GL. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006;72:5069–72.View ArticlePubMedPubMed CentralGoogle Scholar
- Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27:2194–200.View ArticlePubMedPubMed CentralGoogle Scholar
- Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Pena AG, Goodrich JK, Gordon JI, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Shannon CE. The mathematical theory of communication. 1963. MD Comput. 1997;14:306–17.PubMedGoogle Scholar
- Faith DP, Baker AM. Phylogenetic diversity (PD) and biodiversity conservation: some bioinformatics challenges. Evol Bioinform. 2006;2:121–8.Google Scholar
- Lozupone C, Lladser ME, Knights D, Stombaugh J, Knight R. UniFrac: an effective distance metric for microbial community comparison. ISME J. 2011;5:169–72.View ArticlePubMedGoogle Scholar
- Langille MGI, Zaneveld J, Caporaso JG, McDonald D, Knights D, Reyes JA, Clemente JC, Burkepile DE, Thurber RLV, Knight R, et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotechnol. 2013;31:814–21.View ArticlePubMedPubMed CentralGoogle Scholar
- Kanehisa M, Goto S, Sato Y, Kawashima M, Furumichi M, Tanabe M. Data, information, knowledge and principle: back to metabolism in KEGG. Nucleic Acids Res. 2014;42:D199–205.View ArticlePubMedGoogle Scholar
- Dixon P. VEGAN, a package of R functions for community ecology. J Veg Sci. 2003;14:927–30.View ArticleGoogle Scholar