Abstract
In this study, a next-generation sequencing strategy on 16S ribosomal RNA (16S rRNA) gene was employed to analyze 70 oral samples from 32 patients with throat cancer, nine patients with vocal cord polyp, and 29 healthy individuals (normal controls). Using this strategy, we demonstrated, for the first time, that the salivary microbiota of cancer patients were significantly different from those of patients with a polyp and healthy individuals. We observed that the beta diversity of the cancer group was divergent from both the normal and polyp groups, while alpha-diversity indices such as the Chao1 estimator (P = 8.1e-05), Simpson (P = 0.0045), and Shannon (P = 0.0071) were significantly reduced in cancer patients compared with patients containing a polyp and normal healthy individuals. Linear discriminant analysis (LDA) and Kruskal–Wallis test analyses and real-time quantitative polymerase chain reaction (qPCR) verification test revealed that the genera Aggregatibacter, Pseudomonas, Bacteroides, and Ruminiclostridium were significantly enriched in the throat cancer group compared with the vocal cord polyp and normal control groups (score value >2). Finally, diagnostic models based on putatively important constituent bacteria were constructed with 87.5% accuracy [area under the curve (AUC) = 0.875, 95% confidence interval (CI): 0.695–1]. In summary, in this study we characterized, for the first time, the oral microbiota of throat cancer patients without smoking history. We speculate that these results will help in the pathogenic mechanism and early diagnosis of throat cancer.
Introduction
Throat cancer is a type of head and neck carcinoma. This cancer type can be subcategorized into cancer of the oropharynx, hypopharynx, nasopharynx, and larynx depending on where the cancer initially occurred (Edwards and Mendes, 2011). Hypopharyngeal cancer has the highest mortality rate among all the head and neck carcinoma types. The latter cancer has one of the worst prognoses of any malignant head and neck cancer, and the 5-year survival rate is ~15–65% (Edwards and Mendes, 2011; Krstevska et al., 2012; Garden et al., 2016). The 5-year survival rate for early-stage (T1–T2) patients is ~60%, while an overall survival rate of only 17–32% is observed for late-stage patients (Zeng et al., 2010). Furthermore, hypopharyngeal cancer has the highest rate of metastases, the earliest incidence of nodal metastases, and the highest rate of distant metastases (Elias et al., 1995; Johansen et al., 2000). Unfortunately, throat cancers, especially hypopharyngeal cancer, are difficult to detect at the early stages, and by the time the disease is diagnosed, it is usually at an advanced stage (Elias et al., 1995; Johansen et al., 2000; Garden, 2001). Thus, research into the mechanisms that underpin throat cancer and strategies that facilitate the detection of specific biomarkers for early prevention and diagnosis is urgent.
Use of tobacco and (or) alcohol are the highest risk factors for throat cancer, with ~85% of patients having a history of tobacco usage (Anantharaman et al., 2011). Patients with laryngeal cancer and a history of tobacco usage have a 20% higher mortality incidence than laryngeal cancer patients with no history of smoking (Maier et al., 1992). Additionally, infection with human papilloma virus (HPV) is another factor that contributes to throat carcinogenesis (Fakhry et al., 2006). However, the association between HPV and throat cancer has been less well-studied (Mineta et al., 1998; Joo et al., 2014).
The relationship between microbiota composition and cancer occurrence has been extensively studied in recent years. Studies have revealed that oral, intestinal, gastric, and throat mucosal layers are colonized by human commensal bacterial populations, which play important roles in human health (Shin et al., 2017; Dejea et al., 2018; Gao L. et al., 2018). Variations in microbiota diversity and abundance have been shown to be associated with cancer (Tözün and Vardareli, 2016; Yang et al., 2018). Helicobacter pylori infection is an important pathogenic risk factor contributing to the development of gastric cancer (Correa et al., 1975; Lee et al., 2010; Maldonado-Contreras et al., 2011). Helicobacter pylori infection can lead to gastric dysbiosis, which induces cancer development (Bruno et al., 2018). In addition, changes in the gut microbiota are also associated with the development of colorectal carcinogenesis (Nakatsu et al., 2015). The relative pathogenic microbiome has been identified and validated (Amitay et al., 2017; Wong et al., 2017; Chung et al., 2018), and these microbial biomarkers can be used as an aid for early diagnosis of colorectal cancer. Interestingly, the oral microbiota may be diagnostic indicators of oral cancer (Mager et al., 2005). The saliva bacteria Porphyromonas gingivalis, Fusobacterium nucleatum, and Prevotella intermedia have been reported to significantly associate with oral squamous cell carcinoma (Mager et al., 2005; Gholizadeh et al., 2016). Besides, Gong et al. (2017) observed that an imbalance in the throat microbiota profile can increase the risk of laryngeal cancer (Gong et al., 2017). Hayes et al. (2018) demonstrated that the relative abundance of oral bacteria including Corynebacterium and Kingella was negatively associated with the occurrence of laryngeal cancer (Hayes et al., 2018). Although little has been studied, previous studies indicated that the oral bacterium may play an essential role in throat carcinogenesis. Further studies need to be done to confirm the relationship between the human microbiome and throat cancer.
In this study, we investigated the link between the oral microbiome and throat cancer while also attempting to identify bacterial taxa that might represent risk factors in relation to throat carcinoma. We hope that these results provide a platform to further investigate the pathogenic mechanism of oral bacteria for throat cancer.
Patients and Methods
Patients
The study was approved by ethics committee of Beijing Tongren Hospital, Capital Medical University (approval number TRECKY2018-007). All participants upon recruitment signed informed consent. A total of 70 subjects were enrolled from April 17, 2017, to January 12, 2018, including 32 throat cancer patients, nine vocal cord polyp patients, and 29 healthy controls (sample information is in Supplementary Table 1). The individuals in the throat cancer and vocal cord polyp groups were between 45 and 71 years of age; the control group was assembled with individuals of similar age. The majority of the cases were male due to disease morbidity (Gale et al., 2010). All recruited participants had no history of tobacco or alcohol use. Subjects who used antibiotics in the preceding 6 months were excluded from the study.
Saliva specimens were collected from subjects in Beijing Tongren Hospital by using the Oragene DNA Sample Collection Kit (OG-500), which permits saliva specimens to be preserved for years at room temperature without DNA degradation. A volume of 2 ml of saliva was collected from each subject and mixed well with the Oragene solution. All samples were stored at −80°C for further analysis.
DNA Extraction and Sequencing
Oral bacterial genomic DNA was isolated using the QIAamp DNA Mini Kit (Qiagen, Hilden, Germany), following the manufacturer's instructions according to DNA purification protocol for blood/body fluids. Quantification of isolated DNA was performed using a Qubit Fluorometer (Invitrogen, Life Technologies, Grand Island, NY, USA). The V3 sequence of the 16S ribosomal RNA (16S rRNA) gene was amplified using V3 region-specific primers (forward primer; CCTACGGGNGGCWGCAG: reverse primer; ATTCCGCGGCTGGCA) (Lane, 1991). The V3 region sequences were sequenced on the Illumina Miseq system.
Sequence Analysis
Filtering steps were performed to remove low-quality reads. The filtered 16S rRNA gene sequencing data were analyzed by Quantitative Insights Into Microbial Ecology (QIIME) (V1.9) (Caporaso et al., 2010). The reads were clustered into operational taxonomic units (OTUs) using UCLUST with 97% similarity. The sequences were taxonomically assigned based on the Greengenes Database (DeSantis et al., 2006). Alpha-diversity analyses were performed using alpha diversity. The Shannon index, Simpson diversity index, the Abundance-based Coverage Estimator (ACE), Chao1 index, and Good's coverage index were analyzed. Differences in alpha diversity (single index) were compared using the t-test controlled with 103 Monte Carlo permutations. Principal coordinate analysis (PCoA) was calculated based on unweighted and weighted UniFrac distance matrices to observe the differences between individuals or groups; PC1 represents the principal coordinate component that can explain the variation in data as much as possible, and PC2 is the principal coordinate component that explains the largest proportion in the remaining degree of variation. The most important differences in relative abundances among the groups were analyzed using linear discriminant analysis (LDA) effect size (LEfSe) with the throat cancer-related bacteria. Bacterial taxa with LDA scores >4 and a p-value <0.05 were considered strikingly enriched. Taxa with LDA scores >2 and a p-value <0.05 were regarded to be significantly different.
Real-Time Quantitative Polymerase Chain Reaction
The abundance of Aggregatibacter, Pseudomonas, Bacteroides, Faecalibacterium, and Ruminiclostridium was investigated in the healthy sample, throat cancer, and polyp groups by real-time quantitative polymerase chain reaction (qPCR). A total of 10 saliva samples were randomly selected from each group, and real-time qPCR was performed independently on an ABI 7500 Genetic Analyzer (Thermo Fisher Scientific Inc., Carlsbad, CA, USA).
The PCR profile was as follows: 50°C for 2 min, 95°C for 2 min, 40 cycles of 95°C for 15 s, and 60°C for 1 min. A total of 50 ng of extracted DNA and 0.5 μl of each primer (10 pmol) were added to 10 μl of PowerUP SYBR Green Master Mix (Thermo), and the volume was adjusted to 20 μl following the addition of distilled water. All assays were performed in triplicate for each sample. A negative control consisting of DNA-free water was used for each assay. The real-time qPCR data analysis was performed with the ABI 7500 Software with manually set threshold. The Ct value was <35 for targeted samples and >38 for negative control. Finally, ΔCt was calculated according to the method ΔCt = Cttarget-Ctcontrol and relative abundances = POWER(2, –ΔCt). Primers used for the quantification of Pseudomonas, Faecalibacterium, and Bacteroides and the universal primer in the study were referred from previous publications. The primers of Pseudomonas (forward primer: 5′-CTACGGGAGGCAGCAGTGG-3′, reverse primer: 5′-TCGGTAACGTCAAAACAGCAAAGT-3′) (Purohit et al., 2003), Faecalibacterium (forward primer: 5′-CTAACTACGTGCCAGCAGCC-3′, reverse primer: 5′-GCCTTCGCCACTGGTGTTCC-3′) (Duan et al., 2016), Bacteroides (forward primer: 5′-GAGAGGAAGGTCCCCCAC-3′, reverse primer: 5′-CGCTACTTGGCTGGTTCAG-3′) (Guo et al., 2010), and the universal primer (forward primer: 5′-ACTCCTACGGGAGGCAGCAG-3′, reverse primer: 5′-ATTACCGCGGCTGCTGG-3′) (Fierer et al., 2005) were used, respectively. The sequence encoding 16S rRNA was retrieved from the U.S. National Center for Biotechnology Information (NCBI) GeneBank database to design primers for specific detection of Ruminiclostridium (accession nos. NR_026104, NR_026102, NR_112037, and NR_117165, NR_117164) (forward primer: 5′-GGTGAGTAACGCGTGGGTAA-3′, reverse primer: 5′-AACTAGCTAATCGGACGCGG-3′) and Aggregatibacter (accession nos. CP012067 and CP012959) (forward primer: 5′-ACGGGTGAGTAATGCTTGG-3′, reverse primer: 5′-GAGATCGTCGGCTTGGTAGG-3′). The conserved sequences were selected using ClustalW program, and the primers were designed from this highly conserved region.
Results
Oral Microbiota Profiles Generated by 16S rRNA Sequencing
In this study, 32 patients with throat cancer (17 patients with hypopharyngeal cancer and 15 patients with laryngeal cancer), nine patients with a vocal cord polyp, and 29 healthy individuals (controls) were enrolled (sample information is in Supplementary Table 1). Subjects in different groups were similar with respect to age, race, and body mass index (BMI). The subjects had no history of tobacco or alcohol use. The saliva microbiota sequences were analyzed by barcoded 16S rRNA next-generation sequencing. After filtering poor-quality reads, we collected an average of 138,153 clean reads and 223 OTUs per sample (Supplementary Table 1). We identified 504 OTUs in the throat cancer group, 475 OTUs in the vocal cord polyp group, and 488 OTUs in the normal control group (shown in Figure 1A). Additionally, 468 OTUs were shared among these groups, and 16, 1, and 2 unique OTUs were observed for the throat cancer, vocal cord polyp, and normal control groups, respectively (Supplementary Table 2).
Figure 1
The total oral microbiota OTU composition was found similar with that of a previous study in the Human Oral Microbiome Database (HOMD) (Chen et al., 2010), and the dominant phyla and genera are shown in Figures 1B,C. The class, order, and family profiles for all samples are shown in Supplementary Figure 1. The most prevalent phyla among the throat cancer, vocal cord polyp, and normal healthy groups were Proteobacteria, Bacteroidetes, Firmicutes, Actinobacteria, and Fusobacteria in the subjects analyzed, while Streptococcus, Klebsiella, Neisseria, Prevotella, Haemophilus, Capnocytophaga, Veillonella, Porphyromonas, and Rothia were among the most prevalent genera.
Variations in Oral Microbiome Composition Among Different Groups
Alpha diversity and beta diversity were evaluated to compare diversity in the oral microbiome among the different groups. The samples from the throat cancer patients exhibited remarkably reduced microbial diversity indices compared with the healthy control cases (P = 8.1e-05, 0.0045, and 0.0071 for the Chao1, Simpson, and Shannon indices, respectively) (Figures 2A–C). Beta diversity was calculated by PCoA. Following analysis of unweighted distance matrices (P = 0.016), we observed that the microbiota profile of patients with throat carcinoma was significantly different from that of the other groups (Figure 2D).
Figure 2
LEfSe analysis (Segata et al., 2011) was also performed to compare the salivary microbial community structure of patients with throat cancer to that of vocal cord polyp patients or normal control individuals. In total, 76 taxa differed between the throat cancer group and the healthy control group (LDA scores were >2); 15 taxa were significantly enriched in the throat cancer group (Figures 3A,B). One class (Alphaproteobacteria), two orders (Rickettsiales and Burkholderiales), four families (Bacteroidaceae, Bacteroidales Incertae_Sedis, Burkholderiaceae, and Pseudomonadaceae), and eight genera (including PrevotellaceaeNK3B31group, Bacteroides, Phocaeicola, Lautropia, Ruminiclostridium6, Faecalibacterium, Aggregatibacter, and Pseudomonas) were significantly enriched in the throat cancer group compared with the normal control group. Compared to the vocal cord polyp group, the Pseudomonadaceae family and genera Pseudomonas, Aggregatibacter, and Enterobacter were significantly enriched in the throat cancer group (Figures 3C,D).
Figure 3
Validation and Identification of Diagnostic Markers for Throat Cancer
In order to identify potentially useful markers for the diagnosis of throat cancer, an LEfSe analysis was performed to compare patients with throat cancer to vocal cord polyp patients and normal healthy individuals. A total of 74 taxa were identified with LDA scores >2. Two families (Pseudomonadaceae and Bacteroidaceae) and five genera (Aggregatibacter, Pseudomonas, Bacteroides, Faecalibacterium, and Ruminiclostridium) were significantly enriched in the throat cancer group (Figures 4A,B).
Figure 4
A non-parametric Kruskal–Wallis test was performed to confirm differences in abundances among the throat cancer, vocal cord polyp, and normal healthy groups for each bacterial taxon. The analysis confirmed that Pseudomonas, Aggregatibacter, Bacteroides, Faecalibacterium, and Ruminiclostridium6 were prominently ample in the throat cancer group (P < 0.05) (Table 1). Real-time PCR was performed to further analyze the relative abundance of Pseudomonas, Aggregatibacter, Bacteroides, Faecalibacterium, and Ruminiclostridium6 among the different patient groups. The relative abundance of Pseudomonas, Aggregatibacter, Bacteroides, and Ruminiclostridium6 was significantly high among the throat cancer group (P < 0.05) (Figure 5A). Receiver operating characteristic (ROC) curve analysis was subsequently performed with the putative throat cancer markers (Pseudomonas, Aggregatibacter, Bacteroides, and Ruminiclostridium6). As revealed by the area under the curve [AUC = 0.875, 95% confidence interval (CI): 0.695–1], this model exhibited a robust and statistically significant diagnostic accuracy (Figure 5B).
Table 1
| Genera | p-value | Relative abundance mean | ||
|---|---|---|---|---|
| Health | Cancer | Polyp | ||
| Corynebacterium | 0.007 | 0.002 | 0.002 | 0.002 |
| Corynebacterium1 | 0.041 | 0 | 0 | 0 |
| Rothia | 0 | 0.055 | 0.021 | 0.034 |
| Bacteroides | 0.028 | 0 | 0.002 | 0.001 |
| Phocaeicola | 0.001 | 0 | 0 | 0 |
| Porphyromonas | 0.003 | 0.051 | 0.024 | 0.039 |
| Prevotella1 | 0.003 | 0.002 | 0.001 | 0.009 |
| Prevotella2 | 0.001 | 0.006 | 0.003 | 0.005 |
| Prevotella7 | 0.041 | 0.078 | 0.055 | 0.057 |
| PrevotellaceaeNK3B31group | 0.022 | 0 | 0 | 0 |
| RikenellaceaeRC9gutgroup | 0.022 | 0.001 | 0 | 0.002 |
| Bergeyella | 0.004 | 0.004 | 0.002 | 0.002 |
| DefluviitaleaceaeUCG011 | 0.027 | 0 | 0 | 0.001 |
| Parvimonas | 0.012 | 0.003 | 0.001 | 0.003 |
| Anaerovorax | 0 | 0.001 | 0 | 0 |
| FamilyXIIIUCG001 | 0.008 | 0 | 0 | 0 |
| Eubacteriumnodatumgroup | 0.048 | 0.003 | 0.002 | 0.002 |
| Eubacteriumsaphenumgroup | 0.012 | 0 | 0 | 0.001 |
| Johnsonella | 0.018 | 0.001 | 0.001 | 0.001 |
| Lachnoclostridium5 | 0.01 | 0.006 | 0.004 | 0.004 |
| Roseburia | 0.005 | 0 | 0 | 0 |
| Peptococcus | 0.022 | 0 | 0 | 0.001 |
| Filifactor | 0 | 0.002 | 0.001 | 0.006 |
| Peptostreptococcus | 0.005 | 0.006 | 0.002 | 0.004 |
| Eubacteriumyuriigroup | 0 | 0.001 | 0.001 | 0.001 |
| Faecalibacterium | 0.005 | 0 | 0 | 0 |
| Ruminiclostridium6 | 0.022 | 0 | 0 | 0 |
| RuminococcaceaeUCG014 | 0.024 | 0.003 | 0.002 | 0.005 |
| Bulleidia | 0.045 | 0 | 0 | 0 |
| Veillonella | 0.026 | 0.044 | 0.033 | 0.029 |
| Fusobacterium | 0.001 | 0.023 | 0.012 | 0.029 |
| Leptotrichia | 0.023 | 0.012 | 0.008 | 0.014 |
| Lautropia | 0.007 | 0.005 | 0.005 | 0.008 |
| Neisseria | 0.016 | 0.142 | 0.114 | 0.063 |
| Helicobacter | 0.036 | 0 | 0 | 0 |
| Actinobacillus | 0 | 0.009 | 0.003 | 0.004 |
| Aggregatibacter | 0.001 | 0.012 | 0.016 | 0.009 |
| Haemophilus | 0.002 | 0.081 | 0.043 | 0.048 |
| Moraxella | 0.026 | 0.004 | 0.001 | 0 |
| Pseudomonas | 0.001 | 0 | 0.016 | 0 |
| CandidatusSaccharimonas | 0.007 | 0.002 | 0.001 | 0.002 |
| Treponema2 | 0.017 | 0.005 | 0.004 | 0.01 |
| Acholeplasma | 0.016 | 0.001 | 0.001 | 0.001 |
The difference in abundance for each bacterial taxon in samples from the throat cancer patients, vocal cord polyp patients, and normal control group.
Figure 5
Discussion
Relevance between oral microbiome and throat cancer has been reported in a previous study (Hayes et al., 2018). The relative abundance of Corynebactrium and Kingella was associated to the risk of laryngeal tumor with smoking usage. As commonly known, smoking was the major risk factor for the incidence of throat cancer (Lubin et al., 2010; Anantharaman et al., 2011; Amma et al., 2013). It has been reported that smoking could influence the oral microbiota colonization and lead to upper respiratory disease (Charlson et al., 2010; Yu et al., 2017). Smoking may affect the correlation of the oral bacteria with laryngeal cancer, and the correlation needs to be studied independently.
In the current study, we identified the oral bacterial markers for diagnosis of throat cancer, using a cohort of subjects with either hypopharyngeal carcinoma or laryngeal cancer, vocal cord polyp, and normal healthy individuals. All the subjects had no history of smoking. This result revealed, for the first time, that the oral bacterial features of throat cancer patients were different from those of vocal cord polyp patients and normal healthy individuals. The PCoA analysis of unweighted distance matrices revealed that the microbiota profile for the throat cancer group was significantly different from that of the healthy control group. In addition, we also observed significant reductions in microbiota diversity and richness in the throat tumor samples. Statistical analysis revealed that the Chao1 (P = 8.1e-05), Simpson (P = 0.0045), and Shannon (P = 0.0071) indices were significantly reduced in the saliva samples of tumor patients. Previous studies have shown that reductions in the microbial diversity of the human microbiome are linked to cancer occurrence (Schmidt et al., 2014; Gao J. et al., 2018). The associated reduction in bacterial diversity reflected the altering of the abundance in certain bacteria and their metabolic pathways. The “oncogenic bacteria” have been reported to promote tumorigenesis through the indirect inflammatory effects and/or the direct epigenetic mechanism (Sauid and Lois, 2015; Atarashi et al., 2017; Chung et al., 2018). The oral “oncogenic bacteria” associated with throat cancer may migrate with the saliva flow and colonize in the throat by swallowing, which will induce cancer indirectly or directly as a consequence.
Through the LDA analysis on throat cancer vs. normal group, throat cancer vs. vocal polyp, and throat cancer vs. the other two matching groups, bacterial taxa, which might cause throat cancer, were identified. We found that the genera PrevotellaceaeNK3B31group, Bacteroides, Phocaeicola, Lautropia, Ruminiclostridium, Faecalibacterium, Aggregatibacter, and Pseudomonas were significantly enriched in the throat cancer group. These taxa associated only with throat cancer, and they may play important roles in the occurrence of cancer. Compared to the vocal cord polyp group, Pseudomonas, Aggregatibacter, and Enterobacter were significantly enriched in the throat cancer group. In the throat cancer groups, Pseudomonadaceae, Bacteroidaceae, Aggregatibacter, Pseudomonas, Bacteroides, Faecalibacterium, and Ruminiclostridium had the highest LDA scores. Increase in the abundance of Pseudomonas, Aggregatibacter, Bacteroides, and Ruminiclostridium in the throat cancer group (P < 0.05) was also observed by using the Kruskal–Wallis test and real-time qPCR. However, the LDA analysis had shown some anomalies for several microbial groups. Some bacterial taxa associated with several groups. For example, Rickettsiales and Lautropia were related with throat cancer, and Spirochaetales was health-associated according to the LDA analysis on throat cancer vs. health group, while they appeared as associated with vocal polyps on the basis of LDA analysis on the three groups. These results indicated that these bacteria may not play important roles in the associated groups.
The ROC models were constructed using the aforementioned obviously different taxa with random combination (data not shown). The ROC model with the taxa Aggregatibacter, Pseudomonas, Bacteroides, and Ruminiclostridium performed better than with the other taxa in the ROC model analysis. It proved to be robust with statistically significant diagnostic accuracy and had an acquired accuracy of 87.5%. These results indicated that Pseudomonas, Aggregatibacter, Bacteroides, and Ruminiclostridium may have a significant correlation with throat cancer. Aggregatibacter, Pseudomonas, and Bacteroides were reported to associate with aggressive periodontitis and oral infection (Leys et al., 2002; Souto et al., 2014; Cordero and Varela-Calviño, 2018). Periodontitis and oral infection are risk factors for head and neck cancer (Xian-Tao et al., 2013a,b, 2015; Moritani et al., 2015). Pseudomonas, a gram-negative bacterium with carbon degradation activity, was the most common pathogen for hospital infection (Hota et al., 2009; Decker and Palmore, 2014; Hung et al., 2016). It can invade the epithelium and induce immune responses through activation of protein kinase C alpha (PKCα), c-Jun N-terminal kinase (JNK), extracellular-regulated protein kinases (ERK1/2), nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB), and/or glutamic acid–leucine–arginine–positive CXC chemokines (Gregson et al., 2013; Chiang-Wen et al., 2018; Curran et al., 2018). Aggregatibacter was associated with aggressive periodontitis (Ando et al., 2010), which can induce inflammatory response through cytolethal distending toxin (CDT), leukotoxin, and lipopolysaccharide (LPS) (Herbert et al., 2016). Bacteroides fragilis has been verified to correlate with diarrheal disease and inflammatory bowel and can induce colon carcinogenesis through indirect inflammatory effects and direct epigenetic pathways (Toprak et al., 2006; Chung et al., 2018). It can cleave catenins and lead activation of oncogene c-myc (Hardy et al., 2000). As a conclusion, Aggregatibacter, Pseudomonas, and Bacteroides may have individual or combined inducing effects on tumorigenesis through indirect inflammatory and (or) direct epigenetic pathway. Further study is required to confirm the relevant mechanism between tumorigenesis and infection of these bacteria.
In summary, we performed a comprehensive comparison of the oral bacterial taxa among throat cancer patients, vocal cord polyp patients, and healthy individuals. The study revealed different oral microbiota profiles among the analyzed groups. The analysis also detected reduced diversity and richness in the oral bacterial community in throat cancer patients, indicating that oral bacteria may play an important role in the occurrence of throat cancer. A total of four genera, Pseudomonas, Aggregatibacter, Bacteroides, and Ruminiclostridium, were identified significantly associated with throat cancer. However, the pathogenesis that cause this phenomenon still needs to be clarified, and large cohort studies are required to verify the results from this study.
Statements
Ethics statement
The study was approved by ethics committee of Beijing Tongren Hospital, Capital Medical University (approval number TRECKY2018-007). All participants upon recruitment signed informed consent.
Author contributions
LW, ZH, WG, and TS conceived and designed the study and contributed to the writing of the manuscript. GY and WG helped in the data acquisition. LW, GY, YG, YZ, MZ, YL, and PS analyzed and interpreted the data. All the authors read and approved the final manuscript.
Funding
This project was supported by the National Natural Science Foundation of China (no. 81670946) and Beijing Municipal Administration of Hospitals Incubating Program (PX2017032) Capital's Funds for Health Improvement and Research (CFH2018-1-2052).
Acknowledgments
We thank Yunfeng Duan, CAS Key Laboratory of Pathogenic Microbiology and Immunology, Institute of Microbiology, Chinese Academy of Sciences, for the help in the bioinformatics analysis.
Conflict of interest
LW, YG, YZ, MZ, YL, PS, and TS were employed by Beijing Cheer Land Biotechnology Co., Ltd., CL Investment Group. The remaining 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/fcimb.2019.00205/full#supplementary-material
Table S1Sample information in the study.
Table S2Unique OTUs in throat cancer group, vocal cord polyp group, and normal control group.
Supplementary Figure 1The oral microbiota profiles of samples from throat cancer patients, vocal cord polyp patients, and healthy individuals. (A) The relative abundance of the dominant classes in the throat cancer, vocal cord polyp, and normal control groups. (B) The relative abundance of the dominant orders in the throat cancer, vocal cord polyp, and normal control groups. (C) The relative abundance of the dominant families in the throat cancer, vocal cord polyp, and normal control groups.
References
1
AmitayE.WernerS.VitalM.PieperD.HöflerD.GierseI.et al. (2017). Fusobacterium and colorectal cancer: causal factor or passenger? Results from a large colorectal cancer screening study. Carcinogenesis38, 781–788. 10.1093/carcin/bgx053
2
AmmaJ. P.AthiraN.SuminoriA.PalethG.ChihayaK.ChrisB. (2013). Associations of tobacco use and alcohol drinking with laryngeal and hypopharyngeal cancer risks among men in Karunagappally, Kerala, India -Karunagappally cohort study. PLoS ONE8:e73716. 10.1371/journal.pone.0073716
3
AnantharamanD.MarronM.LagiouP.SamoliE.AhrensW.PohlabelnH.et al. (2011). Population attributable risk of tobacco and alcohol for upper aerodigestive tract cancer. Oral Oncol.47, 725–731. 10.1016/j.oraloncology.2011.05.004
4
AndoE. S.De-GennaroL. A.FaveriM.FeresM.DirienzoJ. M.MayerM. P. A. (2010). Immune response to cytolethal distending toxin of Aggregatibacter actinomycetemcomitans in periodontitis patients. J. Periodontol. Res.45, 471–480. 10.1111/j.1600-0765.2009.01260.x
5
AtarashiK.SudaW.LuoC.KawaguchiT.MotooI.NarushimaS.et al. (2017). Ectopic colonization of oral bacteria in the intestine drives T1 cell induction and inflammation. Science358, 359–365. 10.1126/science.aan4526
6
BrunoG.RoccoG.ZaccariP.PorowskaB.MascellinoM.SeveriC. (2018). Infection and gastric dysbiosis: can probiotics administration be useful to treat this condition?Can. J. Infect Dis. Med. Microbiol.2018:6237239. 10.1155/2018/6237239
7
CaporasoJ.KuczynskiJ.StombaughJ.BittingerK.BushmanF.CostelloE.et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods7, 335–336. 10.1038/nmeth.f.303
8
CharlsonE. S.ChenJ.Custers-AllenR.BittingerK.LiH.SinhaR.et al. (2010). Disordered microbial communities in the upper respiratory tract of cigarette smokers. PLoS ONE5:e15216. 10.1371/journal.pone.0015216
9
ChenT.YuW. H.IzardJ.BaranovaO. V.LakshmananA.DewhirstF. E. (2010). The human oral microbiome database: a web accessible resource for investigating oral microbe taxonomic and genomic information. Database J. Biol. Databases and Curation. 10.1093/database/baq013
10
Chiang-WenL.Cheng-HsunW.Yao-ChangC.Yuh-LienC.Kuo-TingC.Chu-ChunC.et al. (2018). Carbon monoxide releasing molecule-2 attenuates Pseudomonas aeruginosa-induced ROS-dependent ICAM-1 expression in human pulmonary alveolar epithelial cells. Redox Boil.18, 93–103. 10.1016/j.redox.2018.07.001
11
ChungL.Thiele OrbergE.GeisA.ChanJ.FuK.DeStefano ShieldsC.et al. (2018). Bacteroides fragilis toxin coordinates a pro-carcinogenic inflammatory cascade via targeting of colonic epithelial cells. Cell. Host Microbe23, 203–214.e205. 10.1016/j.chom.2018.01.007
12
CorderoO. J.Varela-CalviñoR. (2018). Oral hygiene might prevent cancer. Heliyon4:e00879. 10.1016/j.heliyon.2018.e00879
13
CorreaP.HaenszelW.CuelloC.TannenbaumS.ArcherM. (1975). A model for gastric cancer epidemiology. Lancet2, 58–60. 10.1016/S0140-6736(75)90498-5
14
CurranC. S.BoligT.Torabi-PariziP. (2018). Mechanisms and targeted therapies for Pseudomonas aeruginosa lung infection. Am. J. Respir. Crit. Care Med.197, 708–727. 10.1164/rccm.201705-1043SO
15
DeckerB. K.PalmoreT. N. (2014). Hospital water and opportunities for infection prevention. Curr. Infect. Dis. Rep.16:432. 10.1007/s11908-014-0432-y
16
DejeaC.FathiP.CraigJ.BoleijA.TaddeseR.GeisA.et al. (2018). Patients with familial adenomatous polyposis harbor colonic biofilms containing tumorigenic bacteria. Science359, 592–597. 10.1126/science.aah3648
17
DeSantisT.HugenholtzP.LarsenN.RojasM.BrodieE.KellerK.et al. (2006). Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl. Environ. Microbiol.72, 5069–5072. 10.1128/AEM.03006-05
18
DuanC.CuiY.ZhaoY.ZhaiJ.ZhangB.ZhangK.et al. (2016). Evaluation of faecalibacterium 16s rDNA genetic markers for accurate identification of swine faecal waste by quantitative PCR. J. Environ. Manag.181, 193–200. 10.1016/j.jenvman.2016.06.022
19
EdwardsA. A.MendesR. L. (2011). Head and neck cancer, in Specialist Training in Oncology, eds AjithkumarT. V.HatcherH. M. (London: Mosby Ltd), 77–78. 10.1016/B978-0-7234-3458-0.00013-0
20
EliasM. M.HilgersF. J.KeusR. B.GregorR. T.HartA. A.BalmA. J. (1995). Carcinoma of the pyriform sinus: a retrospective analysis of treatment results over a 20-year period. Clin. Otolaryngol. Allied Sci.20, 249–253. 10.1111/j.1365-2273.1995.tb01860.x
21
FakhryC.D'SouzaG.SugarE.WeberK.GoshuE.MinkoffH.et al. (2006). Relationship between prevalent oral and cervical human papillomavirus infections in human immunodeficiency virus-positive and -negative women. J. Clin. Microbiol.44, 4479–4485. 10.1128/JCM.01321-06
22
FiererN.JacksonJ. A.VilgalysR.JacksonR. B. (2005). Assessment of soil microbial community structure by use of taxon-specific quantitative PCR assays. Appl. Environ. Microbiol.71, 4117–4120. 10.1128/AEM.71.7.4117-4120.2005
23
GaleN.MichaelsL.LuzarB.PoljakM.ZidarN.FischingerJ.et al. (2010). Current review on squamous intraepithelial lesions of the larynx. Histopathology54, 639–656. 10.1111/j.1365-2559.2008.03111.x
24
GaoJ.ZhangY.GerhardM.Mejias-LuqueR.ZhangL.ViethM.et al. (2018). Association between gut microbiota and -related gastric lesions in a high-risk population of gastric cancer. Front. Cell. Infect. Microbiol.8:202. 10.3389/fcimb.2018.00202
25
GaoL.XuT.HuangG.JiangS.GuY.ChenF. (2018). Oral microbiomes: more and more importance in oral cavity and whole body. Protein Cell9, 488–500. 10.1007/s13238-018-0548-1
26
GardenA. S. (2001). Organ preservation for carcinoma of the larynx and hypopharynx. Hematol. Oncol. Clin. North Am.15, 243–260. 10.1016/S0889-8588(05)70210-0
27
GardenA. S.BeadleB. M.MorrisonW. H. (2016). Chapter 35—larynx and hypopharynx cancer, in Clinical Radiation oncology, 4th Edn., eds GundersonL. L.TepperJ. E. (Philadelphia, PA: Elsevier), 649–672.e644. 10.1016/B978-0-323-24098-7.00035-6
28
GholizadehP.EslamiH.YousefiM.AsgharzadehM.AghazadehM.KafilH. S. (2016). Role of oral microbiome on oral cancers, a review. Biomed. Pharmacother.84, 552–558. 10.1016/j.biopha.2016.09.082
29
GongH.ShiY.XiaoX.CaoP.WuC.TaoL.et al. (2017). Alterations of microbiota structure in the larynx relevant to laryngeal carcinoma. Sci. Rep.7:5507. 10.1038/s41598-017-05576-7
30
GregsonA. L.WangX.WeigtS. S.PalchevskiyV.LynchJ. P.RossD. J.et al. (2013). Interaction between Pseudomonas and CXC chemokines increases risk of bronchiolitis obliterans syndrome and death in lung transplantation. Am. J. Respir. Crit. Care Med.187, 518–526. 10.1164/rccm.201207-1228OC
31
GuoX.XiaX.TangR.ZhouJ.ZhaoH.WangK. (2010). Development of a real-time PCR method for Firmicutes and Bacteroidetes in faeces and its application to quantify intestinal population of obese and lean pigs. Lett. Appl. Microbiol.47, 367–373. 10.1111/j.1472-765X.2008.02408.x
32
HardyR. G.MeltzerS. J.JankowskiJ. A. (2000). ABC of colorectal cancer. Molecular basis for risk factors. BMJ321, 886–889. 10.1136/bmj.321.7265.886
33
HayesR. B.AhnJ.FanX.PetersB. A.MaY.YangL.et al. (2018). Association of oral microbiome with risk for incident head and neck squamous cell cancer. JAMA Oncol.4, 358–365. 10.1001/jamaoncol.2017.4777
34
HerbertB. A.NovinceC. M.KirkwoodK. L. (2016). Aggregatibacter actinomycetemcomitans, a potent immunoregulator of the periodontal host defense system and alveolar bone homeostasis. Mol. Oral Microbiol.31, 207–227. 10.1111/omi.12119
35
HotaS.HirjiZ.StocktonK.LemieuxC.DedierH.WolfaardtG.et al. (2009). Outbreak of multidrug-resistant Pseudomonas aeruginosa colonization and infection secondary to imperfect intensive care unit room design. Infect Control Hosp Epidemiol.30, 25–33. 10.1086/592700
36
HungC. S.ZingarelliS.NadeauL. J.BiffingerJ. C.DrakeC. A.CrouchA. L.et al. (2016). Carbon catabolite repression and impranil polyurethane degradation in Pseudomonas protegens strain Pf-5. Appl. Environ. Microbiol.82, 6080–6090. 10.1128/AEM.01448-16
37
JohansenL. V.GrauC.OvergaardJ. (2000). Hypopharyngeal squamous cell carcinoma: treatment results in 138 consecutively admitted patients. Acta Radiol. Ther. Phys. Biol.39, 529–536. 10.1080/028418600750013465
38
JooY.YooI.LeeY.ChoK.ParkJ.NamI.et al. (2014). Association between the standardized uptake value and high-risk HPV in hypopharyngeal squamous cell carcinoma. Acta Otolaryngol.134, 1062–1070. 10.3109/00016489.2014.905701
39
KrstevskaV.StojkovskiI.Zafirova-IvanovskaB. (2012). Concurrent radiochemotherapy in locally–regionally advanced oropharyngeal squamous cell carcinoma: analysis of treatment results and prognostic factors. Radiat. Oncol.7:78. 10.1186/1748-717X-7-78
40
LaneD. J. (1991). 16S/23S rRNA sequencing, in Nucleic Acidtechniques in Bacterial Systematics, eds StackebrandtE.GoodfellowM. (New York, NY: John Willey and Sons), 115–175.
41
LeeI.KimJ.ChoiY.PillingerM.KimS.BlaserM.et al. (2010). Helicobacter pylori CagA phosphorylation status determines the gp130-activated SHP2/ERK and JAK/STAT signal transduction pathways in gastric epithelial cells. J. Biol. Chem.285, 16042–16050. 10.1074/jbc.M110.111054
42
LeysE. J.LyonsS. R.MoeschbergerM. L.RumpfR. W.GriffenA. L. (2002). Association of Bacteroides forsythus and a novel Bacteroides phylotype with periodontitis. J. Clin. Microbiol.40, 821–825. 10.1128/JCM.40.3.821-825.2002
43
LubinJ. H.GaudetM. M.OlshanA. F.KelseyK.BoffettaP.BrennanP.et al. (2010). Body mass index, cigarette smoking, and alcohol consumption and cancers of the oral cavity, pharynx, and larynx: modeling odds ratios in pooled case-control data. Am. J. Epidemiol.171, 1250–1261. 10.1093/aje/kwq088
44
MagerD. L.HaffajeeA. D.DevlinP. M.NorrisC. M.PosnerM. R.GoodsonJ. M. (2005). The salivary microbiota as a diagnostic indicator of oral cancer: a descriptive, non-randomized study of cancer-free and oral squamous cell carcinoma subjects. J. Transl. Med.3:27. 10.1186/1479-5876-3-27
45
MaierH.DietzA.GewelkeU.HellerW. D.WeidauerH. (1992). Tobacco and alcohol in the risk of head and neck cancer. Clin. Investig.70, 320–327. 10.1007/BF00184668
46
Maldonado-ContrerasA.GoldfarbK.Godoy-VitorinoF.KaraozU.ContrerasM.BlaserM.et al. (2011). Structure of the human gastric bacterial community in relation to Helicobacter pylori status. ISME J.5, 574–579. 10.1038/ismej.2010.149
47
MinetaH.OginoT.AmanoH. M.OhkawaY.ArakiK.TakebayashiS.et al. (1998). Human papilloma virus (HPV) type 16 and 18 detected in head and neck squamous cell carcinoma. Anticancer Res.18, 4765–4768.
48
MoritaniK.TakeshitaT.ShibataY.NinomiyaT.KiyoharaY.YamashitaY. (2015). Acetaldehyde production by major oral microbes. Oral Dis.21, 748–754. 10.1111/odi.12341
49
NakatsuG.LiX.ZhouH.ShengJ.WongS.WuW.et al. (2015). Gut mucosal microbiome across stages of colorectal carcinogenesis. Nat. Commun.6:8727. 10.1038/ncomms9727
50
PurohitH. J.RajeD. V.KapleyA. (2003). Identification of signature and primers specific to genus Pseudomonas using mismatched patterns of 16s rDNA sequences. BMC Bioinformatics4:19. 10.1186/1471-2105-4-19
51
SauidI.LoisN. (2015). Helicobacter pylori and gastric cancer: a state of the art review. Gastroenterol. Hepatol. Bed Bench8, S6–S14.
52
SchmidtB.KuczynskiJ.BhattacharyaA.HueyB.CorbyP.QueirozE.et al. (2014). Changes in abundance of oral microbiota associated with oral cancer. PLoS ONE9:e98741. 10.1371/journal.pone.0098741
53
SegataN.IzardJ.WaldronL.GeversD.MiropolskyL.GarrettW.et al. (2011). Metagenomic biomarker discovery and explanation. Genome Biol.12:R60. 10.1186/gb-2011-12-6-r60
54
ShinJ.LuoT.KamarajanP.FennoJ.RickardA.KapilaY. (2017). Microbial communities associated with primary and metastatic head and neck squamous cell carcinoma—a high fusobacterial and low streptococcal signature. Sci. Rep.7:9934. 10.1038/s41598-017-09786-x
55
SoutoR.Silva-BoghossianC. M.ColomboA. P. V. (2014). Prevalence of Pseudomonas aeruginosa and Acinetobacter spp. in subgingival biofilm and saliva of subjects with chronic periodontal infection. Braz. J. Microbiol. 45, 495–501. 10.1590/S1517-83822014000200017
56
ToprakN. U.YagciA.GulluogluB. M.AkinM. L.DemirkalemP.CelenkT.et al. (2006). A possible role of Bacteroides fragilis enterotoxin in the aetiology of colorectal cancer. Clin. Microbiol. Infect.12, 782–786. 10.1111/j.1469-0691.2006.01494.x
57
TözünN.VardareliE. (2016). Gut microbiome and gastrointestinal cancer: les liaisons dangereuses. J. Clin. Gastroenterol.50, S191–S196. 10.1097/MCG.0000000000000714
58
WongS.ZhaoL.ZhangX.NakatsuG.HanJ.XuW.et al. (2017). Gavage of fecal samples from patients with colorectal cancer promotes intestinal carcinogenesis in germ-free and conventional mice. Gastroenterology153, 1621–1633.e1626. 10.1053/j.gastro.2017.08.022
59
Xian-TaoZ.Ai-PingD.ChengL.Ling-YunX.Yu-MingN.Wei-DongL.et al. (2013a). Periodontal disease and risk of head and neck cancer: a meta-analysis of observational studies. PLoS One8:e79017. 10.1371/journal.pone.0079017
60
Xian-TaoZ.WeiL.WeiH.QuanW.YiG.Wei-DongL.et al. (2013b). Tooth loss and head and neck cancer: a meta-analysis of observational studies. PLoS One8:e79074. 10.1371/journal.pone.0079074
61
Xian-TaoZ.Wei-DongL.ChaoZ.JingL.Shi-YiC.WeiH. (2015). Meta-analysis on the association between toothbrushing and head and neck cancer. Oral Oncol.51, 446–451. 10.1016/j.oraloncology.2015.02.095
62
YangS.HuangH.FanW.JongY.ChenM.HuangC.et al. (2018). Compositional and functional variations of oral microbiota associated with the mutational changes in oral cancer. Oral Oncol.77, 1–8. 10.1016/j.oraloncology.2017.12.005
63
YuG.PhillipsS.GailM. H.GoedertJ. J.HumphrysM. S.RavelJ.et al. (2017). The effect of cigarette smoking on the oral and nasal microbiota. Microbiome5:3. 10.1186/s40168-016-0226-6
64
ZengW.QiuG.ZhouH.PengJ.ChenM.TanS. N.et al. (2010). Community structure and dynamics of the free and attached microorganisms during moderately thermophilic bioleaching of chalcopyrite concentrate. Bioresour. Technol.101, 7079–7086. 10.1016/j.biortech.2010.04.003
Summary
Keywords
throat cancer, early diagnosis, oral health, microbiome, 16S rRNA gene, next-generation sequencing
Citation
Wang L, Yin G, Guo Y, Zhao Y, Zhao M, Lai Y, Sui P, Shi T, Guo W and Huang Z (2019) Variations in Oral Microbiota Composition Are Associated With a Risk of Throat Cancer. Front. Cell. Infect. Microbiol. 9:205. doi: 10.3389/fcimb.2019.00205
Received
16 November 2018
Accepted
27 May 2019
Published
03 July 2019
Volume
9 - 2019
Edited by
Ulvi Kahraman Gürsoy, University of Turku, Finland
Reviewed by
J. Christopher Fenno, University of Michigan, United States; Gilad Bachrach, Hadassah Medical Center, Israel
Updates
Copyright
© 2019 Wang, Yin, Guo, Zhao, Zhao, Lai, Sui, Shi, Guo and Huang.
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: Taiping Shi tpshi@126.comWei Guo guowei_1983@163.comZhigang Huang huangzhigang1963@sohu.com
This article was submitted to Microbiome in Health and Disease, a section of the journal Frontiers in Cellular and Infection Microbiology
†These authors have contributed equally to this work
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.