Skip to main content

Oral-microbiome-derived signatures enable non-invasive diagnosis of laryngeal cancers

Abstract

Background

Recent studies have uncovered that the microbiota in patients with head and neck cancers is significantly altered and may drive cancer development. However, there is limited data to explore the unique microbiota of laryngeal squamous cell carcinoma (LSCC), and little is known regarding whether the oral microbiota can be utilized as an early diagnostic biomarker.

Methods

Using 16S rRNA gene sequencing, we characterized the microbiome of oral rinse and tissue samples from 77 patients with LSCC and 76 control patients with vocal polyps, and then performed bioinformatic analyses to identify taxonomic groups associated with clinicopathologic features.

Results

Multiple bacterial genera exhibited significant differences in relative abundance when stratifying by histologic and tissue type. By exploiting the distinct microbial abundance and identifying the tumor-associated microbiota taxa between patients of LSCC and vocal polyps, we developed a predictive classifier by using rinse microbiota as key features for the diagnosis of LSCC with 85.7% accuracy.

Conclusion

This is the first evidence of taxonomical features based on the oral rinse microbiome that could diagnose LSCC. Our results revealed the oral rinse microbiome is an understudied source of clinical variation and represents a potential non-evasive biomarker of LSCC.

Introduction

The human body harbors trillions of microorganisms with various functions that are closely related to our health. Being one of the largest habitats of microorganisms, there are more than 1000 different kinds of microorganisms in the oral cavity. Plenty of studies showed that altered oral microbial profile and their metabolites are associated with head and neck squamous cell carcinoma (HNSCC) [1, 2] and esophageal cancer [3], and they might also influence remote organs through the digestive tract [4], thus contributing to colorectal cancer [5], pancreatic cancer [6], and lung cancer [7]. Through extensive studies, some specific taxa were identified to be related to oral/oropharyngeal cancers, such as Fusobacterium, Peptostreptococcus, Prevotella, Veillonella, Capnocytophaga, etc. Several mechanisms of carcinogenic action of oral bacteria have been proposed [8,9,10]. The first is bacterial stimulation of chronic inflammation. Anaerobic species such as Fusobacterium, Prevotella, and Porphyromonas are responsible for periodontal diseases and lead to chronic inflammation that may facilitate cell proliferation, angiogenesis, and oncogene activation. The second mechanism involves antiapoptotic activity. Bacteria such as Porphyromonas gingivalis inhibit apoptosis by modulating several pathways [11]. The third mechanism is the direct or indirect production of carcinogenic substances. Several oral bacteria are capable of metabolizing alcohol to carcinogenic acetaldehyde [12]. Besides, genotoxic metabolites of certain bacteria could cause DNA damage or produce free radicals and affect reactive oxygen species (ROS) [13, 14].

Laryngeal squamous cell carcinoma (LSCC) is the second most prevalent in all respiratory system carcinomas after lung cancer [15]. LSCC is also common among HNSCC. Approximately 184,615 new cases and 99,840 deaths were estimated attributing to LSCC worldwide in 2020 [15]. However, due to the lack of early detection, most patients have developed advanced LSCC upon diagnosis [16, 17], which always leads to regional recurrence and distant metastase even after surgical resection. The pathogenesis of LSCC is not fully uncovered, and several risk factors are related to LSCC, such as smoking, alcohol, and human papillomavirus (HPV) infection.

Recently, several studies of head and neck cancer showed that alterations in microbiota may drive HNSCCs, and some potentially oncogenic bacteria are identified [9, 18, 19]. However, prior work investigating the microbiome of head and neck cancer mainly focused on oral/oropharyngeal cancer. There is very limited data to explore the unique microbiota of laryngeal cancer, in addition, whether the oral microbiome could serve as a biomarker for the non-invasive early detection of laryngeal cancer has not been studied. The oral microbiome from oral rinse has been considered as a good tool to study the microbiota of oral/oropharyngeal cancer due to the close anatomical sites [1, 20]. It is also a promising biomarker for the early detection of cancer with lower cost and noninvasive access. As the oral microbiota can migrate to the laryngopharyngeal region through inhalation and swallowing, this study aimed to identify alterations in the microbiome of LSCC patients from control-group patients with tissue samples and oral rinse samples, and explore whether oral microbiota can be applied as non-evasive diagnostic markers in LSCC patients.

Materials and methods

Clinical samples and database

Patients were enrolled in this study at Peking Union Medical College Hospital (PUMCH) from 2020 to 2022. All the patients were divided into two groups: (1) tumor group, 77 patients who had pathologically confirmed, previously untreated LSCC undergoing radical resection; (2) control group, which was composed of 76 patients who had vocal cord polyps matched by age and gender. The exclusion criteria were as follows: (i) antibiotics therapy in a month; (ii) infection with HBV, HCV, syphilis, or HIV; (iii) a history of malignant tumors, chemotherapy, or radiotherapy. TNM stages of all participants were identified according to National Comprehensive Cancer Network (NCCN) Guidelines in 2021. Clinical information such as age, sex, tobacco, and alcohol consumption was obtained from the medical history of the participants. “Smokers” were defined as patients who had smoked > 100 cigarettes during their lifetime; otherwise, they were defined as “never smokers”. “Drinkers” were defined as those who drank alcohol at least once a week for a year or longer; whereas “never drinkers” had drunk less. Before sampling, participants were banned from dieting, smoking, and oral hygiene prophylaxis for at least 2 h. Participants rinsed the mouth vigorously with 10 ml sterile saline for 30 s in the operating room to sterilely collect oral rinses. Then, cancerous tissues or polyp tissues from the central area of the lesions were obtained in the surgery. Tissue samples were placed into sterile 2 ml Eppendorf tubes (Axygen, USA) and then frozen at − 80 degrees Celsius before further processing.

16S rRNA sequencing data analysis

Total genomic DNA was extracted from the specimens using the QIAampFast DNA Stool Mini Kit (Qiagen, Hilden, Germany). The V3-V4 region of the bacterial 16S rRNA gene was amplified with PCR using the forward primer 341F (5′-barcode-CCTAYGGGRBGCASCAG-3′) and reverse primer 806R (5′-GGACTACNNGGGTATCTAAT-3′). PCR was performed in a mixture of 30 μl as previously described. Pooled PCR products were extracted using 2% agarose gels and then purified using a GeneJET Gel Extraction Kit (Thermo Fisher Scientific). NEB Next® Ultra™ DNA Library Prep Kit (New England Biolabs) was applied for generation of the sequencing libraries, and the index codes were added. The qualified library was sequenced on an Illumina NovaSeq platform, and 250 bp paired-end reads were generated according to the standard protocols.

Bioinformatics analysis

Paired-end reads were merged into tags with FLASH (version 1.2.7) [21]. Raw tags were filtered using the command split_libraries_fastq.py in QIIME (version 1.8) [22]. Chimeric tags were removed with USEARCH (version 6.1) [23] after quality control. The satisfied reads were clustered into operational taxonomic units (OTUs) at 97% similarity using QIIME. Taxonomic information on OTUs was obtained with the Ribosomal Database Project classifier [24] against the Greengenes database (release 13.8) [25]. Alpha diversity indices and β-diversity distance matrices measuring the pairwise difference among samples were calculated with QIIME. Differences in the Chao1 index, Shannon index, Simpson index, and observed OTUs were detected using the Wilcoxon rank sum test. Principal coordinates analysis (PCoA) was conducted using R software (version 4.0.2) [26], PERMANOVA was performed with the vegan [27] package in R to clarify differences in microbial communities between groups. The Wilcoxon rank sum test was applied to identify differentially abundant taxa between groups, and a false discovery rate < 0.05 was considered to be statistically significant.

Using genus profiles, samples were randomly divided into a training set and a test set by using the “sample” function (80% of samples as training dataset and 20% of samples as testing dataset). Only genera that were present in at least 10% of subjects were considered in the analyses. The number of selected biomarkers were determined after tenfold cross-validation that was repeated for five times based on the training set profile. Then, the final random forest classification models were constructed, and the mean decrease accuracy of selected markers were used to evaluate the importance of biomarkers. The area under the curve (AUC) was used to measure the performance of the models when applied to the training and test sets using the ROCR (version 1.0–11) [28] package in R, and receiver-operating characteristic (ROC) curves were calculated using the pROC (version 1.17.0.1) [29] package in R. The code for bioinformatics analysis can be requested at https://github.com/Jora1991/Oral-microbiome-for-laryngeal-cancers.

Results

Diversity dysbiosis of the tumor and oral microbiome in LSCC

A total of 153 patients were enrolled in this study at Peking Union Medical College Hospital (PUMCH) from 2020 to 2022. Most patients were between 56 and 70 years old (66.0%), and average age of patients was 60.4. The male/female ratio was 8.6:1 (137 males:16 females). According to cancer statistics, laryngeal cancer occurs more commonly in men than in women, most frequently diagnosed among people aged 55–74 [15, 30]. The male/female ratio and age range of LSCC patients in this study approximately corresponded with the epidemiology of laryngeal cancer in literature. The baseline characteristics of the entire cohort are presented in Table 1. As is shown, no significant differences were found between the LSCC group and the control group regarding age or gender. Smokers and drinkers were more common in the LSCC group, which was inevitable. Community alpha diversity was assessed with observed taxonomic units (OTUs), Simpson, Shannon, and Chao1 indices (Fig. 1A). Statistical P values of alpha diversity indices among different groups are presented in Additional file 1: Table S1. The oral rinse samples exhibited more elevated microbial diversity indices compared with the tissue samples in both the tumor group and control group (P value < 0.05 for the Simpson index, Shannon index, and observed OTUs). The Simpson index and Shannon index were significantly higher in the tumor tissues than in the control tissues. Although the oral rinse samples did not show significant differences in α-diversity between the tumor and control group, it is worth noting that the in-group variation of oral rinse samples was lower than that of tissue samples. Then, principal coordinate analysis (PCoA) was performed based on the Bray–Curtis distance matrix from the genus profile to investigate β-diversity. As shown in Fig. 1B, the oral rinse samples were clustered closely with each other and separated from the tissue samples, in which a clear distinction between tumor tissue and control tissue samples was found based on permutational multivariate analysis of variance (PERMANOVA). In summary, oral rinse samples presented lower within-group variation compared to laryngeal tissue samples.

Table 1 Demographics and clinical characteristics of the study cohort
Fig. 1
figure 1

Microbiota diversity in LSCC (N = 77) and control samples (N = 76). A. Comparison of microbial α-diversity between groups. P values were calculated from the Wilcoxon rank sum test. P < 0.01 was labeled with “+”, and “*” represented P < 0.05. B Principal coordinate analysis (PCoA) plots. The P value was derived from permutational multivariate analysis of variance (PERMANOVA)

Microbial composition alterations of LSCC tumor and oral rinse samples

Subsequently, the relative abundance of the top 20 taxa was compared between specimens from LSCC patients and specimens from control patients (Fig. 2). At the genus level, tumor tissues contained a much higher percentage of Fusobacterium (0.088 ± 0.115 (mean ± SD) vs. 0.022 ± 0.028, tumor vs. control tissues), Pseudomonas (0.083 ± 0.160 vs. 0.055 ± 0.101), and Acinetobacter (0.055 ± 0.127 vs. 0.003 ± 0.013) compared to control tissues, whereas Ralstonia (0.016 ± 0.026 vs. 0.233 ± 0.241) were prominent constituents of most control tissues, but were scarce in tumor tissues. Oral rinse samples displayed similar microbiome composition and were characterized by the same general in both the tumor group and control group (Prevotella, Neisseria, Streptococcus, Haemophilus, and Alloprevotella), which also implied lower within-group variation compared to laryngeal tissue samples. As expected, significant differences in the microbiota composition at the genus level were discovered among different groups of samples (Fig. 3A). Compared with the control tissues, the tumor tissues exhibited a marked increase in the abundance of several genera, such as Flavobacterium (P = 0.008, 0.003 ± 0.007 vs. 0.000 ± 0.000, tumor vs. control tissues), Klebsiella (P = 0.002, 0.017 ± 0.038 vs. 0.000 ± 0.001), Mycoplasma (P = 0.001, 0.002 ± 0.011 vs. 0.000 ± 0.000), etc., while the abundances of Ralstonia (P = 0.000, 0.012 ± 0.026 vs. 0.233 ± 0.241, tumor vs. control tissues), Actinobacillus (P = 0.000, 0.000 ± 0.000 vs. 0.002 ± 0.005), Streptococcus (P = 0.000, 0.024 ± 0.039 vs. 0.151 ± 0.238), Lactobacillus (P = 0.000, 0.002 ± 0.006 vs. 0.010 ± 0.035), Rothia (P = 0.002, 0.003 ± 0.006 vs. 0.015 ± 0.041), and several other genera were significantly lower in the tumor tissues. In the oral rinse samples, Saccharopolyspora (P = 0.018, 0.001 ± 0.003 vs. 0.000 ± 0.002, tumor vs. control rinses) and Actinobacillus (P = 0.008, 0.019 ± 0.050 vs. 0.005 ± 0.010) were more abundant from tumor patients than those from control patients. The heatmap is colored according to the relative abundances (Fig. 3B) showing the genera differences among the four sample groups (the tumor tissue samples, the control tissue samples, the oral rinse samples from LSCC patients, and the oral rinse samples from control patients). We found that the tissue samples were separated from the oral rinse samples and that the tumor tissue samples were distinct from the control tissue samples. The genera Veillonella, Fusobacterium, Ralstonia, etc., contributed to this separation. Additionally, 39 genera were observed only in the differential analysis of the tissue samples (Leptotrichia, Veillonella, Solobacterium, etc.; shown in Fig. 3B). Three unique genera, Alistipes, Selenomonas, and Saccharopolysora were observed only in the differential analysis of the oral rinse samples. Four genera, L.NK4A136, Alicycliphilus, Actinobacillus, and Moraxella were detected in the differential analysis of both the oral rinse and tissue samples. More details about the percentage values and statistical P values of genera among different groups are presented in Additional file 1: Table S2.

Fig. 2
figure 2

Dominant genera in LSCC (N = 77) and control samples (N = 76). The relative abundance of the top 20 dominant genera in different groups were presented. Samples of each group were clustered and ordered based on the abundance profile of top 20 genera using “hclust” function in R

Fig. 3
figure 3

Differential microbial composition in LSCC tumor tissues and oral rinse samples. A Log2-fold change (Log2FC) of the genera was significantly different between tumor tissue and control tissue samples (above panel); Log2FC of the genera was significantly different between oral rinse samples from tumor patients and oral rinse samples from control patients (lower panel). Log2FC > 0 indicated that the genus was more abundant in tumor patients than in control patients, while log2FC < 0 indicated that the genus was less abundant in tumor patients than in control patients. B Colored heatmap according to the abundances of differential genera in A

Classification model for LSCC diagnosis based on oral rinse microbiota

Our results suggested that the abundance of several taxa was significantly different in the oral rinse samples between the LSCC group and the control group. Therefore, the oral microbiota has the potential to be a predictive marker indicating a risk for LSCC. To test this hypothesis, microbial biomarkers were selected first at the genus level using tenfold cross-validation that was repeated five times based on a training set of 121 samples (Table 1). Then, two random forest (RF) models were trained from the selected biomarkers utilizing the training sets from tissue samples (N = 61) and oral rinse samples (N = 60). The performances of the constructed models were evaluated by the area under the curve (AUC) of receiver operating characteristic (ROC) curves. A total of 15 genera, including Candidatus Saccharimonas, Ralstonia, Moraxella, Solobacterium, Neisseria, Streptococcus, Veillonella, and 8 others, were identified as microbial biomarkers that could discriminate LSCC patients from control patients in both tissue samples (Fig. 4A) and oral rinse samples (Fig. 4B), indicating the potential use of microbial features in the identification of LSCC. Specifically, the model for tissue samples exhibited excellent performance based on the training set (AUC = 98.1%) and the test set (N = 16, AUC = 83.3%) (Fig. 4A). Interestingly, compared with the previous model using tissue samples, the RF model for oral rinse samples showed a higher AUC based on the test set (N = 16, AUC = 85.7%) accompanied by a lower AUC based on the training set (AUC = 86.4%) (Fig. 4B), supporting the potential diagnostic value of our oral microbiota-based classifier for LSCC. Moreover, we also examined the potential ability of the 15 microbial biomarkers to discriminate LSCC based on PCoA and PERMANOVA. Significant separations (PERMANOVA, P < 0.05) between the LSCC and control groups were observed.

Fig. 4
figure 4

Identification of microbial-based markers for clinical diagnosis by random forest models. A Mean decrease accuracy (MDA) of variables from the RF model using the profile of 15 selected biomarkers in tissue samples (left). ROC curves are based on the training set and test set from tissue samples (middle). The AUC is displayed in the lower right-hand corner. PCoA based on the profile of 15 selected genera in tissue samples (right). The P value is derived from PERMANOVA. B MDA of 15 selected biomarkers (left), ROC curves, and AUC (middle) based on the RF model from oral rinse samples. PCoA and PERMANOVA based on the profile of 15 selected genera in oral rinse samples (right)

Discussion

While the role of oral microbiota in oral/oropharyngeal cancer has received much attention, few data in the literature investigated the association between microbial alteration and laryngeal cancer [31,32,33]. None of the studies included both oral rinse/saliva samples and tissue samples to study the alterations of oral microbiota and laryngeal microbial composition in LSCC patients. The present study included both sample types to investigate whether oral rinse samples can replace tissue samples in microbiome studies and be applied as early diagnostic markers in LSCC patients. Our study showed that although oral rinse samples exhibited remarkably elevated microbial diversity indices compared with the tissue samples, and the oral rinse samples did not show significant differences in α-diversity or β-diversity between the tumor and control group, oral rinse samples presented lower within-group variation compared with laryngeal tissue samples, which indicated their potential character as a stable diagnostic biomarker. Compared with the control tissues, the abundance of several genera was significantly higher in the laryngeal tissue of the LSCC group, including some notorious genera such as Fusobacterium, which is in accordance with previous studies [31, 34]. Fusobacterium is an invasive anaerobe mainly associated with periodontitis, which may lead to chronic inflammation and contribute to carcinogenesis. Recently, more and more studies have found an association between the abundance of fusobacterium and cancer, including oral/ oropharyngeal cancer [35], colorectal cancer [36, 37], and so on. In the oral rinse samples, taxa differences between the tumor group and the control group were not as large as those based on the tissue samples, however, the RF model using microbial biomarkers from oral rinse samples exhibited excellent performance based on the training set and the test set, indicating the diagnostic value of our oral microbiota-based classifier for LSCC.

Of note, no previous study had proven that a microbiota-based model could reliably detect laryngeal cancer through oral rinse samples with high sensitivity and specificity. We identified taxonomic features that yielded a model and developed a classifier, which was able to distinguish cancer status. Our results demonstrate that oral microbiome shifts might serve as a marker in advance of the clinical phenotypes detectable by endoscope or CT scans. As this study describes cohorts of patients with diagnosed laryngeal cancer, future work should examine if this signature can be applied to assess cancer onset in individuals with precancerous lesions such as vocal leukoplakia. This is a good clinical endpoint as the model enables non-invasive LSCC detection with a simple taxa quantification test on oral rinse samples. Our findings provide a strong rationale for developing microbiome-based liquid biopsy technology to prioritize at-risk individuals for clinical attention.

Although this is not the pilot study to reveal the specific microbial composition of LSCC patients, our study is the first one to include multiple sample types and investigate whether oral rinse samples can replace tissue samples in microbiome studies. Several unavoidable limitations inherent in this study should be acknowledged. First, this was a retrospective study conducted in a single institution and therefore involved possible selection bias. Second, although age and gender are paired in two groups, the influence of smoking and alcohol consumption on microbiota cannot be eliminated in this study. Tabbco and alcohol consumption could alter oral microbial composition [9, 38], and might lead to the natural selection of microbiota capable of a high rate of carcinogen metabolism, which may synergize with the primary risk factors such as alcohol abuse and smoking in cancer pathogenesis [39,40,41]. In the present study, there are no genera differences in oral rinses between smokers and non-smokers of control group, and only one genus (Ralstonia) was found in microbial differential analysis in oral rinses between drinkers and non-drinkers of control group. In brief, we believe that smoking and drinking showed limited influence on microbiota from oral rinses in this study, so the classification model for LSCC diagnosis based on microbiota from oral rinses is solid. Further studies conducted with more patients will be performed to allow sub-group analysis.

Conclusion

Through the analysis of alpha and beta diversity, we found that oral rinse samples presented lower within-group variation compared with laryngeal tissue samples, which indicated their potential character as a stable diagnostic biomarker. Then, our results proved that dysbiosis of the oral microbiome is a key hallmark of LSCC, and the oral microbiota-based model exhibited excellent performance in the discrimination of LSCC samples. To our knowledge, this is the first study to provide a strategy for the non-invasive prediction of LSCC by employing an oral microbiota-based model, which represents a promising early diagnostic biomarker of LSCC patients for future clinical use.

Availability of data and materials

The datasets used for the current study are available upon reasonable request from the corresponding author.

Abbreviations

LSCC:

Laryngeal squamous cell carcinoma

HNSCC:

Head and neck squamous cell carcinoma

OTU:

Observed taxonomic units

RF:

Random forest

AUC:

Area under the curve

ROC:

Receiver operating characteristic curves

PCoA:

Principal coordinate analysis

References

  1. Ganly I, Yang L, Giese RA, Hao Y, Nossa CW, Morris LGT, Rosenthal M, Migliacci J, Kelly D, Tseng W, Hu J, Li H, Brown S, Pei Z. Periodontal pathogens are a risk factor of oral cavity squamous cell carcinoma, independent of tobacco and alcohol and human papillomavirus. Int J Cancer. 2019;145(3):775–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Wolf A, Moissl-Eichinger C, Perras A, Koskinen K, Tomazic PV, Thurnher D. The salivary microbiome as an indicator of carcinogenesis in patients with oropharyngeal squamous cell carcinoma: a pilot study. Sci Rep. 2017;7(1):5867.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Peters BA, Wu J, Pei Z, Yang L, Purdue MP, Freedman ND, Jacobs EJ, Gapstur SM, Hayes RB, Ahn J. Oral microbiome composition reflects prospective risk for esophageal cancers. Cancer Res. 2017;77(23):6777–87.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Karpinski TM. Role of Oral Microbiota in Cancer Development. Microorganisms. 2019;7(1):20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Kostic AD, Chun E, Robertson L, Glickman JN, Gallini CA, Michaud M, Clancy TE, Chung DC, Lochhead P, Hold GL, El-Omar EM, Brenner D, Fuchs CS, Meyerson M, Garrett WS. Fusobacterium nucleatum potentiates intestinal tumorigenesis and modulates the tumor-immune microenvironment. Cell Host Microbe. 2013;14(2):207–15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Farrell JJ, Zhang L, Zhou H, Chia D, Elashoff D, Akin D, Paster BJ, Joshipura K, Wong DT. Variations of oral microbiota are associated with pancreatic diseases including pancreatic cancer. Gut. 2012;61(4):582–8.

    Article  CAS  PubMed  Google Scholar 

  7. Yan X, Yang M, Liu J, Gao R, Hu J, Li J, Zhang L, Shi Y, Guo H, Cheng J, Razi M, Pang S, Yu X, Hu S. Discovery and validation of potential bacterial biomarkers for lung cancer. Am J Cancer Res. 2015;5(10):3111–22.

    CAS  PubMed  PubMed Central  Google Scholar 

  8. Zhang Y, Wang X, Li H, Ni C, Du Z, Yan F. Human oral microbiota and its modulation for oral health. Biomed Pharmacother. 2018;99:883–93.

    Article  PubMed  Google Scholar 

  9. Hayes RB, Ahn J, Fan X, Peters BA, Ma Y, Yang L, Agalliu I, Burk RD, Ganly I, Purdue MP, Freedman ND, Gapstur SM, Pei Z. Association of oral microbiome with risk for incident head and neck squamous cell cancer. JAMA Oncol. 2018;4(3):358–65.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Bhatt AP, Redinbo MR, Bultman SJ. The role of the microbiome in cancer development and therapy. CA Cancer J Clin. 2017;67(4):326–44.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Gholizadeh P, Eslami H, Yousefi M, Asgharzadeh M, Aghazadeh M, Kafil HS. Role of oral microbiome on oral cancers, a review. Biomed Pharmacother. 2016;84:552–8.

    Article  CAS  PubMed  Google Scholar 

  12. Moritani K, Takeshita T, Shibata Y, Ninomiya T, Kiyohara Y, Yamashita Y. Acetaldehyde production by major oral microbes. Oral Dis. 2015;21(6):748–54.

    Article  CAS  PubMed  Google Scholar 

  13. Huycke MM, Abrams V, Moore DR. Enterococcus faecalis produces extracellular superoxide and hydrogen peroxide that damages colonic epithelial cell DNA. Carcinogenesis. 2002;23(3):529–36.

    Article  CAS  PubMed  Google Scholar 

  14. Goodwin AC, Destefano Shields CE, Wu S, Huso DL, Wu X, Murray-Stewart TR, Hacker-Prietz A, Rabizadeh S, Woster PM, Sears CL, Casero RA Jr. Polyamine catabolism contributes to enterotoxigenic Bacteroides fragilis-induced colon tumorigenesis. Proc Natl Acad Sci U S A. 2011;108(37):15354–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, Bray F. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin. 2021;71(3):209–49.

    Article  PubMed  Google Scholar 

  16. Groome PA, O’Sullivan B, Irish JC, Rothwell DM, Schulze K, Warde PR, Schneider KM, Mackenzie RG, Hodson DI, Hammond JA, Gulavita SP, Eapen LJ, Dixon PF, Bissett RJ, Mackillop WJ. Management and outcome differences in supraglottic cancer between Ontario, Canada, and the Surveillance, Epidemiology, and End Results areas of the United States. J Clin Oncol. 2003;21(3):496–505.

    Article  PubMed  Google Scholar 

  17. Marur S, Forastiere AA. Head and neck cancer: changing epidemiology, diagnosis, and treatment. Mayo Clin Proc. 2008;83(4):489–501.

    Article  PubMed  Google Scholar 

  18. Wang H, Funchain P, Bebek G, Altemus J, Zhang H, Niazi F, Peterson C, Lee WT, Burkey BB, Eng C. Microbiomic differences in tumor and paired-normal tissue in head and neck squamous cell carcinomas. Genome Med. 2017;9(1):14.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Frank DN, Qiu Y, Cao Y, Zhang S, Lu L, Kofonow JM, Robertson CE, Liu Y, Wang H, Levens CL, Kuhn KA, Song J, Ramakrishnan VR, Lu SL. A dysbiotic microbiome promotes head and neck squamous cell carcinoma. Oncogene. 2022;41(9):1269–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Lim Y, Fukuma N, Totsika M, Kenny L, Morrison M, Punyadeera C. The performance of an oral microbiome biomarker panel in predicting oral cavity and oropharyngeal cancers. Front Cell Infect Microbiol. 2018;8:267.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Magoc T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. 2011;27(21):2957–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet CC, Al-Ghalith GA, Alexander H, Alm EJ, Arumugam M, Asnicar F, Bai Y, Bisanz JE, Bittinger K, Brejnrod A, Brislawn CJ, Brown CT, Callahan BJ, Caraballo-Rodriguez AM, Chase J, Cope EK, Da Silva R, Diener C, Dorrestein PC, Douglas GM, Durall DM, Duvallet C, Edwardson CF, Ernst M, Estaki M, Fouquier J, Gauglitz JM, Gibbons SM, Gibson DL, Gonzalez A, Gorlick K, Guo J, Hillmann B, Holmes S, Holste H, Huttenhower C, Huttley GA, Janssen S, Jarmusch AK, Jiang L, Kaehler BD, Kang KB, Keefe CR, Keim P, Kelley ST, Knights D, Koester I, Kosciolek T, Kreps J, Langille MGI, Lee J, Ley R, Liu YX, Loftfield E, Lozupone C, Maher M, Marotz C, Martin BD, McDonald D, McIver LJ, Melnik AV, Metcalf JL, Morgan SC, Morton JT, Naimey AT, Navas-Molina JA, Nothias LF, Orchanian SB, Pearson T, Peoples SL, Petras D, Preuss ML, Pruesse E, Rasmussen LB, Rivers A, Robeson MS 2nd, Rosenthal P, Segata N, Shaffer M, Shiffer A, Sinha R, Song SJ, Spear JR, Swafford AD, Thompson LR, Torres PJ, Trinh P, Tripathi A, Turnbaugh PJ, Ul-Hasan S, van der Hooft JJJ, Vargas F, Vazquez-Baeza Y, Vogtmann E, von Hippel M, Walters W, Wan Y, Wang M, Warren J, Weber KC, Williamson CHD, Willis AD, Xu ZZ, Zaneveld JR, Zhang Y, Zhu Q, Knight R, Caporaso JG. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol. 2019;37(8):852–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26(19):2460–1.

    Article  CAS  PubMed  Google Scholar 

  24. Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73(16):5261–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. 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(7):5069–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. R.C.J.M.c. Team, R: A language and environment for statistical computing, 1 (2014).

  27. Dixon P. VEGAN, a package of R functions for community ecology. J Veg Sci. 2003;14(6):927–30.

    Article  Google Scholar 

  28. Sing T, Sander O, Beerenwinkel N, Lengauer T. ROCR: visualizing classifier performance in R. Bioinformatics. 2005;21(20):3940–1.

    Article  CAS  PubMed  Google Scholar 

  29. Robin X, Turck N, Hainard A, Tiberti N, Lisacek F, Sanchez JC, Muller M. pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinformatics. 2011;12:77.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Steuer CE, El-Deiry M, Parks JR, Higgins KA, Saba NF. An update on larynx cancer. CA Cancer J Clin. 2017;67(1):31–50.

    Article  PubMed  Google Scholar 

  31. Gong HL, Shi Y, Zhou L, Wu CP, Cao PY, Tao L, Xu C, Hou DS, Wang YZ. The composition of microbiome in larynx and the throat biodiversity between laryngeal squamous cell carcinoma patients and control population. PLoS ONE. 2013;8(6): e66476.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Gong H, Shi Y, Zhou X, Wu C, Cao P, Xu C, Hou D, Wang Y, Zhou L. Microbiota in the throat and risk factors for laryngeal carcinoma. Appl Environ Microbiol. 2014;80(23):7356–63.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Wang L, Yin G, Guo Y, Zhao Y, Zhao M, Lai Y, Sui P, Shi T, Guo W, Huang Z. Variations in oral microbiota composition are associated with a risk of throat cancer. Front Cell Infect Microbiol. 2019;9:205.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Gong H, Shi Y, Xiao X, Cao P, Wu C, Tao L, Hou D, Wang Y, Zhou L. Alterations of microbiota structure in the larynx relevant to laryngeal carcinoma. Sci Rep. 2017;7(1):5507.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Li Y, Tan X, Zhao X, Xu Z, Dai W, Duan W, Huang S, Zhang E, Liu J, Zhang S, Yin R, Shi X, Lu Z, Pan Y. Composition and function of oral microbiota between gingival squamous cell carcinoma and periodontitis. Oral Oncol. 2020;107: 104710.

    Article  CAS  PubMed  Google Scholar 

  36. Castellarin M, Warren RL, Freeman JD, Dreolini L, Krzywinski M, Strauss J, Barnes R, Watson P, Allen-Vercoe E, Moore RA, Holt RA. Fusobacterium nucleatum infection is prevalent in human colorectal carcinoma. Genome Res. 2012;22(2):299–306.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Repass J, Maherali N, Owen K, Reproducibility Project: Cancer Biology. Registered report: Fusobacterium nucleatum infection is prevalent in human colorectal carcinoma. Elife. 2016;5:299.

    Article  Google Scholar 

  38. Wu J, Peters BA, Dominianni C, Zhang Y, Pei Z, Yang L, Ma Y, Purdue MP, Jacobs EJ, Gapstur SM, Li H, Alekseyenko AV, Hayes RB, Ahn J. Cigarette smoking and the oral microbiome in a large study of American adults. ISME J. 2016;10(10):2435–46.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Galvao-Moreira LV, da Cruz MC. Oral microbiome, periodontitis and risk of head and neck cancer. Oral Oncol. 2016;53:17–9.

    Article  CAS  PubMed  Google Scholar 

  40. Chen J, Domingue JC, Sears CL. Microbiota dysbiosis in select human cancers: evidence of association and causality. Semin Immunol. 2017;32:25–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Salaspuro V, Salaspuro M. Synergistic effect of alcohol drinking and smoking on in vivo acetaldehyde concentration in saliva. Int J Cancer. 2004;111(4):480–3.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This study is supported by National High Level Hospital Clinical Research Funding (2022-PUMCH-B-094), CAMS Innovation Fund for Medical Sciences (CIFMS) (2021-I2M-1-023/2020-I2M-2-009), and Science & Technology Fundamental Resources Investigation Program (Grant No. 2022FY100800).

Author information

Authors and Affiliations

Authors

Contributions

PZ and XC provided the overall study design and co-wrote the paper. SY designed the study, collected samples, and drafted the paper. JC was responsible for the data analysis with software and cowriting the manuscript. YF, XX, and FY collected the clinical data and co-wrote the paper. YZ and GS were responsible for the optimization of the model and for cowriting the paper. All authors have approved the final version of the manuscript.

Corresponding authors

Correspondence to Peng Zhang or Xingming Chen.

Ethics declarations

Ethics approval and consent to participate

The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013). The current study was approved by the Ethics Committee of Peking Union Medical College Hospital (approval ID number: ZS-3148), Peking Union Medical College, and the Chinese Academy of Medical Sciences, Beijing, China. Written informed consent was obtained from all patients before enrollment.

Consent for publication

Not applicable.

Competing interests

All authors declare no conflicts of interest.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Statistical P values of alpha diversity indices among different groups. Table S2. Percentage values and statistical P values of genera among different groups.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yu, S., Chen, J., Zhao, Y. et al. Oral-microbiome-derived signatures enable non-invasive diagnosis of laryngeal cancers. J Transl Med 21, 438 (2023). https://doi.org/10.1186/s12967-023-04285-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12967-023-04285-2

Keywords