Abstract
Background:
The functional role of lung microbiota has attracted an accumulating attention recently, but the profile and functional role of the lung microbiota in patients with lung cancer remained largely unknown.
Methods:
To evaluate the association of the microbiota with lung cancer, we performed comparative analysis of the lung microbiota using 16S rRNA amplicon sequencing approach in the paired bronchoalveolar lavage fluid (BALF) samples (paired samples from cancerous lung and the contralateral non-cancerous lung) from 50 cancer patients with unilateral lobar masses.
Results:
We found that the relative abundance of phylum Tenericutes, its class Mollicutes, its order Entomoplasmatales, its family Spiroplasmataceae, and its genus Spiroplasma was significantly increased in cancerous lung, but the relative abundance of phylum Bacteroidetes, its class Bacteroidia, and its order Bacteroidales was significantly decreased in cancerous lung. In addition, the relative abundance of family Leuconostocaceae and its genus Weissella was significantly increased in cancerous lung.
Conclusion:
Our findings provide insights into a change of lung microbiota community associated with the development of lung cancer.
Introduction
Lung cancer is the leading cause of cancer-related morbidity and mortality worldwide (). Multiple studies suggest that several factors, such as smoking (), inflammatory response (), and infections with several microbiota (–), are associated with lung cancer development.
Since mucosal tissue harbors the largest surface area in the human body, the lung is exposed to many airborne microbiota and environmental factors through inhalation. Thus, the relationship between the lung microbiota and lung disease is of particular interest. Healthy human lungs were traditionally deemed to be sterile owing to the inability of culturing bacteria from lower airway samples through routine microbiological methods (). However, since the first culture-independent study found microbiota in asthmatic airways (), series of studies using diverse molecular approaches have indicated that the healthy lung is also home to a complex diversity of microbiota (, ).
The functional role of respiratory microbiota has recently attracted an accumulating attention (, ) and additionally several studies found that variations of local microbial community were associated with the exacerbation of several respiratory disorders, including asthma, chronic obstructive pulmonary disease, and cystic fibrosis (, –).
However, the profile and functional role of the lung microbiota in patients with lung cancer remained largely unknown. To date, three studies used high-throughput sequencing for characterizing the lung microbiota in lung cancer patients, using bronchoalveolar lavage fluid (BALF) samples (), protected specimen brush (), and tissue samples (), respectively. The two studies using BALF and protected specimen brush enrolled fewer patients (n < 30) (, ), and another study collected lung biopsy in non-cancerous tissue. Moreover, in healthy lungs, spatial variation of microbiota within an individual was significantly less than variation across individuals () and thus comparative analysis of the lung microbiota between lung cancer patients and healthy control subjects could produce inaccurate results regarding microbiota related with lung cancer.
To evaluate the association of the microbiota with lung cancer, we performed comparative analysis of the lung microbiota in the paired BALF samples (one from the cancerous lung lobe, the other from the contralateral corresponding non-cancerous lung lobe) from patients with lung cancer.
Materials and Methods
Subject Recruitment and Sample Collection
This study received ethical approval from the Institutional Review Board of the Peking University School of Oncology, China, and informed consent was obtained from all subjects. Fifty cancer patients with unilateral lobar masses were selected from patients consented to bronchoscopy examination at Beijing Cancer Hospital, Beijing, China. Demographic and clinical data, including age, gender, smoking, drinking, family history, and final diagnosis were obtained from each participant. Transbronchoendoscope, which avoided contamination of the upper respiratory tract or oral microbiota, was performed to obtain paired BALF samples in lung cancer patients (one from the cancerous lung, the other from the contralateral non-cancerous lung). All samples were immediately frozen and maintained at −80°C until further DNA extraction.
DNA Extraction, 16S rRNA Amplification and Sequencing
Bacterial DNA was isolated from airway samples using PowerSoil DNA Isolation Kit (MoBio Laboratories, Carlsbad, CA, United States) following the manufacturer’s instructions. PCR was used with unique primers to amplify the V3–V4 region of the 16S rRNA to obtain an amplicon library from BALF samples, and amplicons were sequenced using the Miseq platform, as previously described ().
Sequence Data Analysis
Clean data were trimmed from raw data by trimmomatic 0.33. Operational taxonomic units (OTUs) were classified based on 97% similarity after chimeric sequences were removed using UPARSE (version 2.7.1). The phylogenetic affiliation of each 16S rRNA gene sequence was analyzed by RDP Classifier1 against the Silva (SSU128) 16S rRNA database using a confidence threshold of 70%. α-diversity (within-subject diversity) was assessed on the basis of the non-parametric Shannon Wiener diversity index and Simpson diversity index. The relationship between overall microbiota composition (β-diversity) and lung cancer was assessed by analysis of weighted and unweighted UniFrac distances calculated in QIIME (). To visualize separation of samples based on pairwise distances, principal coordinate analysis (PCoA) plots were generated using the first three principal coordinates. The QIIME pipeline was also used to generate PCoA plots to visualize the unweighted UniFrac dissimilarity. We also broadly compared bacterial taxa (phylum to genus) between cases and controls. We limited our analysis of bacterial phyla to those with mean relative abundance ≥0.01%. For lower-level taxa (class to genus), we limited analysis to those with mean relative abundance ≥0.0001%. Paired t-test and Wilcoxon signed-ranks test were used with false discovery rate (FDR) to adjust P-values for multiple testing to detect taxa with differential abundance among groups. α-diversity difference between cancer and contralateral normal lung was examined in conditional logistic analysis. Bar plots and PCoA plots were all generated in R2.
Accession Number
Raw sequencing data were submitted to the Sequence Read Archive (BioProject ID: PRJNA552813).
Results
Characteristics of the Study Participants
This study prospectively enrolled 50 patients, including 46 primary lung cancer patients and 4 patients with metastatic cancer in lung. The clinical characteristics of patients are listed in Table 1. In total, 34 males and 16 females were enrolled and additionally 29 smokers and 14 drinkers were included. The average age was 60 years old, and 9 patients have a family history of cancer. There were 14 cases (28%) diagnosed with lung squamous cell carcinoma, 12 cases (24%) with lung adenocarcinoma, 7 cases (14%) with non-specified non-small cell lung cancer, 13 cases (26%) with small cell lung cancer, 2 cases with colorectal cancer, 1 case with breast cancer, and 1 case with renal clear cell cancer. Among all patients, 5 patients (10%) were diagnosed in stages IA–IIIA and 44 patients (90%) had more advanced cancer stages (IIIB, IV) (Table 1). In 46 lung cancer patients, none previously had received chemotherapy, radiation therapy, or other treatments.
TABLE 1
| Variable | Value |
| Age | |
| Mean ± SD | 60 ± 9.6 |
| Gender | |
| Male | 34 |
| Female | 16 |
| Smoking | |
| Yes | 29 |
| No | 21 |
| Drinking | |
| Yes | 14 |
| No | 36 |
| Histology type | |
| Adenocarcinoma | 12 |
| Squamous carcinoma | 14 |
| Small cell lung cancer | 13 |
| NOS | 7 |
| Lung metastatic cancer | 4 |
| AJCC 8th stage | |
| IA | 1 |
| IIB | 2 |
| IIIA | 2 |
| IIIB | 16 |
| IV | 28 |
| T stage | |
| T1 | 2 |
| T2 | 12 |
| T3 | 18 |
| T4 | 16 |
| Lymph node involvement | |
| N0 | 2 |
| N1 | 8 |
| N2 | 15 |
| N3 | 23 |
| Family history | |
| Yes | 9 |
| No | 41 |
Characteristics of subjects.
NOS: non-specified non-small cell lung cancer. AJCC: American Joint Committee on Cancer.
Taxonomic Profiles of the Cancer and Normal Lung Microbiota
Bacterial OTUs were classified into 46 and 43 phyla, 123 and 112 classes, 177 and 158 orders, 342 and 313 families, and 690 and 631 genera in cancer lung and normal lung samples, respectively (Supplementary Tables 2–6). Four phyla (Firmicutes, Proteobacteria, Fusobacteria, and Bacteroidetes) were detected in all cancer lung and normal lung samples. Seven genera (Prevotella_7, Gemella, Porphyromonas, Haemophilus, Fusobacterium, Veillonella, and Streptococcus) and 10 genera (Neisseria, Porphyromonas, Fusobacterium, Campylobacter, Alloprevotella, Prevotella_6, Streptococcus, Veillonella, Haemophilus, and Prevotella_7) were detected in all cancer lung and normal lung samples, respectively (Supplementary Tables 1–5). The three dominant phyla, classes, orders, and families were the same among cancer lung and normal lung samples. The two dominant genera were the same, but the third dominant genera were Alloprevotella and Prevotella among cancer lung and normal lung samples, respectively (Figure 1).
FIGURE 1
Overall Microbiota Diversity in Relation to Cancer and Normal Lung
Cancer lung was not significantly different from normal lung in α-diversity (within-subject diversity), as measured by the Shannon (P = 0.871) and Simpson diversity index (P = 0.627), or overall lung microbiota composition (β-diversity, between-subject diversity), as measured by unweighted and weighted UniFrac distances (P = 0.952) (Figure 2).
FIGURE 2
Comparison of Microbiota From Cancer and Normal Lung
We found that greater abundance of phylum Tenericutes (P = 0.007) was associated with increased risk of lung cancer, as well as its class Mollicutes (P = 0.007), its order Entomoplasmatales (P = 0.003), its family Spiroplasmataceae (P = 0.003), and its genus Spiroplasma (P = 0.003). The abundance of phylum Bacteroidetes was associated with reduced risk of lung cancer (P = 0.009), as well as its class Bacteroidia (P = 0.01), and its order Bacteroidales (P = 0.01). Family Leuconostocaceae (P = 0.009) and its genus Weissella (P = 0.009) were positively associated with lung cancer risk, although its phylum Firmicutes, class Bacilli, and order Lactobacillales were not significantly associated with lung cancer risk (Figure 3 and Supplementary Table 6).
FIGURE 3
Discussion
Despite accumulating evidence linking microbiota to cancer, the role of lung microbiota in lung cancer has not been fully elucidated. In this study, we firstly characterized the taxonomic profiles of the microbiota in BALF from bilateral lungs (cancer lung and normal lung) within the same patients with lung cancer and found that increased relative abundance of phylum Tenericutes, its class Mollicutes, its order Entomoplasmatales, its family Spiroplasmataceae, and its genus Spiroplasma could be associated with lung cancer risk, but decreased relative abundance of phylum Bacteroidetes, its class Bacteroidia, and its order Bacteroidales could be related with risk of lung cancer. In addition, increased relative abundance of family Leuconostocaceae and its genus Weissella could be associated with lung cancer risk.
To date, three studies used high-throughput sequencing for characterizing the lung microbiota in lung cancer patients, using BALF samples (20 lung cancer patients) (), protected specimen brush samples (24 lung cancer patients) (), and tissue samples (31 lung cancer patients) (), respectively. Since these studies enrolled less lung patients and used different sampling methods, their findings were quite different.
We found that three dominant phyla were Firmicutes, Proteobacteria, and Bacteroidetes in normal and cancer lung, similar to results of the above three studies (). However, three dominant genera in normal and cancer lung in our study were similar to Lee et al.’s and Liu et al.’s results (, ) but different from Yu et al.’s findings (), indicating that sampling methods could significantly affect microbial detection since our study’s sampling method was similar to Lee et al.’s and Liu et al.’s sampling methods but different from Yu et al.’s sampling method. Since lung biopsy in tumor tissue could be an invasive approach and moreover lung biopsy in non-cancerous tissue and healthy human subjects is not ethical, it is more convenient to detect lung microbiota using BALF or protected specimen brush.
We found that neither α-diversity nor β-diversity was significantly different between cancer and health lung, similar to Liu et al.’s findings (). However, Lee et al.’s () and Yu et al.’s () studies reported significantly higher and lower α-diversity in cancer lung than normal lung, respectively. Similar to our study, Liu et al.’s study compared α-diversity and β-diversity between cancerous lesion and contralateral non-cancerous site within the same patient with lung cancer, but Lee et al.’s () study compared patients with lung cancer and patients with benign diseases and Yu et al.’s () study obtained lung tissue samples not BALF, which indicated that significant variation across different individuals and sampling methods could contribute to controversial findings regarding diversity of lung microbiota.
Among the above three studies (–), each study found specific microbiota associated with lung cancer, but none of these microbiota was repeatedly found to be associated with lung cancer in all these three studies, which could be attributed to small sample size and different sampling methods. We found that increased relative abundance of not only genus Spiroplasma but also its family, order, class, and phylum was significantly associated with the lung cancer risk, which could indicate their reliable associations, although these three previous studies did not find the significant association of these microbial taxa with lung cancer (–). Moreover, we found decreased relative abundance of phylum Bacteroidetes, its class Bacteroidia, and its order Bacteroidales in cancer lung, which was similar to the findings in Vogtmann’s study using oral wash (). In addition, we found increased relative abundance of family Leuconostocaceae and its genus Weissella in cancer lung. Previous studies reported that infections with Weissella spp. primarily occurred in immunocompromised status and/or those with other underlying medical conditions (), which could account for increased abundance of the genus Weissella in patients with lung cancer.
To date, our study had the largest sample of patients with lung cancer, included paired samples derived from cancer and normal lung within the same patients, and used a transbronchoendoscope, which prevented contamination of the upper respiratory tract or oral microbiota to obtain paired BALF. Although our study has the above advantages, it also has some limitations. First, it has a cross-sectional study design and therefore cannot identify a causal relationship between microbiota and lung cancer. Second, although our study had the largest sample of patients with lung cancer until now, the sample size was moderate and thus our study did not evaluate the influence of age, gender, smoking status, antibiotic use, and lung function on lung microbiota. Further large-scale studies are needed to investigate the role of microbiota in patients with lung cancer. Third, although we compared the lung microbiota in the paired BALF samples (paired samples from cancerous lung and the contralateral non-cancerous lung) to reduce variations of lung microbiota across individuals, the variation across individuals still makes it hard to infer whether differences within one individual really reflect causal relationship.
Conclusion
In conclusion, our study found that there were differences in the microbial communities of cancer and normal lung within the same patient with lung cancer, and the genera Spiroplasma and Weissella were significantly enriched in cancer lung. These findings indicate that within individuals with unilateral lung cancers, there are differences in the microbiota between the cancerous and non-cancerous sides, but further large-scale studies and function of microbiota in animal model systems are essential to investigate potential use for prevention, diagnosis, and treatment strategies targeting the lung microbiota.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: Sequence Read Archive (BioProject ID: PRJNA552813) (https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA552813).
Ethics statement
The studies involving human participants were reviewed and approved by the Institutional Review Board of the Peking University School of Oncology. The patients/participants provided their written informed consent to participate in this study.
Author contributions
ZW and CZ designed the research. MZ, TA, and CZ conducted the experiments. CZ and MZ analyzed the data. CZ, ZW, and MZ wrote the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the Beijing Natural Science Foundation (Grant No. 7182029 to CZ and 7162038 to MZ), Capital’s Funds for Health Improvement and Research (Grant No. 2020-4-1028 to CZ), the Beijing Municipal Science & Technology Commission (Grant No. Z181100001718104 to ZW), the Beijing Hospitals Authority Youth Programme (Grant No. QML20191105 to CZ), the Beijing Municipal Administration of Hospitals’ Youth Program (Grant No. QML20161101 to MZ), the Peking University Medicine Fund of Fostering Young Scholars’ Scientific & Technological Innovation (Grant No. BMU2020PYB017 to CZ), and the Science Foundation of Peking University Cancer Hospital (Grant No. 18-02 to ZW and 2020-4 to MZ).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2020.01584/full#supplementary-material
References
1.
BrayFFerlayJSoerjomataramISiegelRLTorreLAJemalA.Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries.CA Cancer J Clin. (2018) 68:394–424. 10.3322/caac.21492
2.
HechtSS.Lung carcinogenesis by tobacco smoke.Int J Cancer. (2012) 131:2724–32. 10.1002/ijc.27816
3.
GomesMTeixeiraALCoelhoAAraujoAMedeirosR.The role of inflammation in lung cancer.Adv Exp Med Biol. (2014) 816:1–23. 10.1007/978-3-0348-0837-8_1
4.
ZhanPSuoLJQianQShenXKQiuLXYuLKet alChlamydia pneumoniae infection and lung cancer risk: a meta-analysis.Eur J Cancer. (2011) 47:742–7. 10.1016/j.ejca.2010.11.003
5.
LinTYHuangWYLinJCLinCLSungFCKaoCHet alIncreased lung cancer risk among patients with pneumococcal pneumonia: a nationwide population-based cohort study.Lung. (2014) 192:159–65. 10.1007/s00408-013-9523-z
6.
FabbrizziAAmedeiALavoriniFRendaTFontanaG.The lung microbiome: clinical and therapeutic implications.Intern Emerg Med. (2019) 14:1241–50. 10.1007/s11739-019-02208-y
7.
O’DwyerDNDicksonRPMooreBB.The lung microbiome, immunity, and the pathogenesis of chronic lung disease.J Immunol. (2016) 196:4839–47. 10.4049/jimmunol.1600279
8.
HiltyMBurkeCPedroHCardenasPBushABossleyCet alDisordered microbial communities in asthmatic airways.PLoS One. (2010) 5:e8578. 10.1371/journal.pone.0008578
9.
DicksonRPErb-DownwardJRMartinezFJHuffnagleGB.The Microbiome and the respiratory tract.Ann Rev Physiol. (2016) 78:481–504. 10.1146/annurev-physiol-021115-105238
10.
DicksonRPHuffnagleGB.The lung microbiome: new principles for respiratory bacteriology in health and disease.PLoS Pathog. (2015) 11:e1004923. 10.1371/journal.ppat.1004923
11.
SzeMADimitriuPASuzukiMMcDonoughJECampbellJDBrothersJFet alHost response to the lung microbiome in chronic obstructive pulmonary disease.Am J Respirat Crit Care Med. (2015) 192:438–45. 10.1164/rccm.201502-0223OC
12.
GatesKMartinezF.The human microbiome in the lung: are infections contributing to lung health and disease?Chronic Obstruct Pulmon Dis. (2016) 3:466–72. 10.15326/jcopdf.3.1.2015.0174
13.
SzeMADimitriuPAHayashiSElliottWMMcDonoughJEGosselinkJVet alThe lung tissue microbiome in chronic obstructive pulmonary disease.Am J Respirat Crit Care Med. (2012) 185:1073–80. 10.1164/rccm.201111-2075OC
14.
MarslandBJGollwitzerES.Host-microorganism interactions in lung diseases.Nat Rev Immunol. (2014) 14:827–35. 10.1038/nri3769
15.
MolyneauxPLWillis-OwenSAGCoxMJJamesPCowmanSLoebingerMet alHost-microbial interactions in idiopathic pulmonary fibrosis.Am J Respirat Crit Care Med. (2017) 195:1640–50. 10.1164/rccm.201607-1408OC
16.
LeeSHSungJYYongDChunJKimSYSongJHet alCharacterization of microbiome in bronchoalveolar lavage fluid of patients with lung cancer comparing with benign mass like lesions.Lung Cancer. (2016) 102:89–95. 10.1016/j.lungcan.2016.10.016
17.
LiuHXTaoLLZhangJZhuYGZhengYLiuDet alDifference of lower airway microbiome in bilateral protected specimen brush between lung cancer patients with unilateral lobar masses and control subjects.Int J Cancer. (2018) 142:769–78. 10.1002/ijc.31098
18.
YuGGailMHConsonniDCarugnoMHumphrysMPesatoriACet alCharacterizing human lung tissue microbiota and its relationship to epidemiological and clinical features.Genome Biol. (2016) 17:163. 10.1186/s13059-016-1021-1
19.
DicksonRPErb-DownwardJRFreemanCMMcCloskeyLBeckJMHuffnagleGBet alSpatial variation in the healthy human lung microbiome and the adapted island model of lung biogeography.Ann Am Thorac Soc. (2015) 12:821–30. 10.1513/AnnalsATS.201501-029OC
20.
ZhangCLiuYGaoWPanYGaoYShenJet alThe direct and indirect association of cervical microbiota with the risk of cervical intraepithelial neoplasia.Cancer Med. (2018) 7:2172–9. 10.1002/cam4.1471
21.
LozuponeCKnightR.UniFrac: a new phylogenetic method for comparing microbial communities.Appl Environ Microbiol. (2005) 71:8228–35. 10.1128/AEM.71.12.8228-8235.2005
22.
VogtmannEGoedertJJ.Epidemiologic studies of the human microbiome and cancer.Br J Cancer. (2016) 114:237–42. 10.1038/bjc.2015.465
23.
FuscoVQueroGMChoGSKabischJMeskeDNeveHet alThe genus Weissella: taxonomy, ecology and biotechnological potential.Front Microbiol. (2015) 6:155. 10.3389/fmicb.2015.00155
Summary
Keywords
microbiota, 16S rRNA, lung cancer, bronchoalveolar lavage fluid, cancerous lung and the contralateral non-cancerous lung
Citation
Zhuo M, An T, Zhang C and Wang Z (2020) Characterization of Microbiota in Cancerous Lung and the Contralateral Non-Cancerous Lung Within Lung Cancer Patients. Front. Oncol. 10:1584. doi: 10.3389/fonc.2020.01584
Received
01 June 2020
Accepted
22 July 2020
Published
24 August 2020
Volume
10 - 2020
Edited by
Noemi Reguart, Hospital ClÃnic de Barcelona, Spain
Reviewed by
Amedeo Amedei, University of Florence, Italy; Desiree Hao, University of Calgary, Canada
Updates
Copyright
© 2020 Zhuo, An, Zhang and Wang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Chaoting Zhang, chaotingzhang@bjmu.edu.cnZiping Wang, wangzp2007@126.com
†These authors have contributed equally to this work
This article was submitted to Thoracic Oncology, a section of the journal Frontiers in Oncology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.