Yuan Liu
Zhengrui Li
Yanxu Qi4,5
Ling Zhang- 1Department of Oral and Maxillofacial-Head and Neck Oncology, Shanghai Ninth People’s Hospital, Shanghai Jiao Tong University School of Medicine, College of Stomatology, Shanghai Jiao Tong University, Shanghai, China
- 2National Center for Stomatology & National Clinical Research Center for Oral Diseases, Shanghai, China
- 3Shanghai Key Laboratory of Stomatology, Shanghai, China
- 4Department of Oral and Maxillofacial Surgery, School and Hospital of Stomatology, Cheeloo College of Medicine, Shandong University, Shandong, China
- 5Shandong Key Laboratory of Oral Tissue Regeneration & Shandong Engineering Laboratory for Dental Materials and Oral Tissue Regeneration, Shandong University, Shandong, China
- 6Department of Oral and Maxillofacial Surgery, College and Hospital of Stomatology, Guangxi Medical University, Nanning, Guangxi, China
The relationship between oral squamous cell carcinoma (OSCC) development and the microbiome has attracted increasing attention. The depth of invasion (DOI) is an important indicator of tumor progression, staging and prognosis, and the change in the oral microbiome based on the DOI is unclear. This report describes the use of metagenomic analyses to investigate the relationship between the oral microbiome and the DOI. Forty patients in different DOI categories were recruited; 10 healthy people served as the control group. Swab samples collected from the participants were subjected to metagenomic analyses, and the oral microbial communities and their functions were investigated. The abundances of Fusobacterium nucleatum, Capnocytophaga sputigena, Porphyromonas endodontalis, and Gemella haemolysans were significantly increased in the patients compared with the controls. The abundances of some bacteria exhibited a stage-related trend. The abundances of P. endodontalis, Gemella morbillorum and G. haemolysans increased with increasing DOI. In contrast, the abundances of Prevotella melaninogenica, Haemophilus parainfluenzae and Neisseria flavescens decreased with increasing DOI. Based on receiver operating characteristic (ROC) curve analysis, eight species were found to have predictive value: Rothia mucilaginosa, P. melaninogenica, H. parainfluenzae, and N. flavescens in the healthy control group and P. endodontalis, G. morbillorum, G. haemolysans and Fusobacterium periodonticum in the high DOI group. In the functional analysis, several metabolic pathways were decreased, whereas flagellar assembly and bacterial chemotaxis showed an increasing trend as the disease progressed. Biofilm formation, flagella, lipopolysaccharide (LPS) and other virulence factors exhibited staging-related changes. These pathogenic pathways and factors had a clear correlation with specific pathogens. In particular, when OSCC progressed to the late stage, microbial diversity and functional potential changed greatly.
Introduction
In 2017, the 8th edition of the Union for International Cancer Control and American Joint Committee on Cancer staging manual incorporated tumor depth of invasion (DOI) in the T category of the TNM (tumor, node, metastasis) staging system for oral squamous cell carcinoma (OSCC) (Amin et al., 2017; Brierley et al., 2017). The DOI is measured from the level of the basement membrane of the closest adjacent normal mucosa, and a “plumb line” is dropped from this plane to the deepest point of tumor invasion (Lydiatt et al., 2017). According to the guidelines, the DOI serves as an important reference for OSCC staging. Compared with the thickness of the tumor, the DOI better reflects the degree of tumor invasion and infiltration. A deep DOI has a statistically significant association with more advanced disease, including a higher T category and N category and extracapsular spread, which is well established as an independent predictor of recurrence and survival (Tan et al., 2012; Ganly et al., 2013; Liu et al., 2016).
Given the important relationship between bacteria and human diseases, an increasing number of studies have linked bacteria and tumors. Recent research has shown that some microbes promote carcinogenesis, such as Helicobacter pylori, which promotes the development of cancer through epithelial injury and inflammation (Smoot, 1997). In addition, the International Agency for Research on Cancer has noted the following as class I carcinogen pathogens: H. pylori, hepatitis B and C viruses, human papillomavirus and Epstein Barr virus (IARC Working Group on the Evaluation of Carcinogenic Risks to Humans, 2012). Other research has indicated that cancers can be promoted by dysbiotic microbiomes (Schwabe and Jobin, 2013); for instance, disorders of intestinal microorganisms can promote colorectal cancer. The main mechanisms include chronic inflammation, immune regulation and microbial metabolites (Schwabe and Jobin, 2013; Meng et al., 2018). The oral cavity, as part of the digestive tract, is also one of the most bacteria-rich parts of the human body. The composition and function of bacteria in OSCC deserve attention because they may play an important role in the development of cancer.
Human bacteria isolated using culture-based methods usually represent nearly 70% of oral microbiota taxa (Dewhirst et al., 2010). Culture-independent analysis using next-generation sequencing fills the gap and is crucial for defining and understanding the microbiome and its key roles in human disease. Studies have used 16S rDNA to analyze microbial diversity in the context of OSCC. 16S amplicons have obvious advantages in the study of microbial diversity, but they are also flawed compared to metagenomic sequence, such as chimeras during PCR amplification and differential resolution across taxa, raising less convincing conclusions (Chen et al., 2021). In addition, although studies have analyzed the microbial diversity associated with different OSCC clinical stages (Yang et al., 2018), this grouping method includes some redundant factors, such as regional lymph node status and distant metastasis. There is insufficient evidence to show that these factors are clearly related to microorganisms. Exploring the microbiome under different DOI backgrounds can link it to the patient’s tumor stage and prognosis.
We determined the DOI through postoperative pathology and then sorted the patients into three groups. Consistent with the guidelines, we chose 5 mm and 10 mm as the categorical limits (Lydiatt et al., 2017). Thus, the DOI of the DOI-low group was less than or equal to 5 mm, that of the DOI-medium group was greater than 5 mm and less than or equal to 10 mm, and that of the DOI-high group was greater than 10 mm. According to the guidelines, the DOI is an important factor in determining tumor staging; thus, the DOI can be viewed as a snapshot of the dynamic changes in tumor progression. We used swabs to sample from the OSCC area for metagenomic analysis. We described the bacterial diversity and functional transitions associated with the DOI. First, we attempted to clarify whether there was a significant change in the overall microbiome or a significant change in the abundance of key bacteria during tumor progression. Second, we explored changes in microbiome function, related metabolic pathways and virulence factors during tumor progression through metagenomics technology. These results provide a basis for subsequent in-depth studies.
Materials and Methods
Study Population and Covariate Assessment
This study was approved by the Ethics Committee of Shanghai Ninth People’s Hospital. Our study has been registered in the Chinese clinical trial registry. The participants were drawn from hospital patients and their family members. The participant covariate information, such as drinking status, smoking status, body mass index (BMI) and sampling sites, was obtained from medical histories and questionnaires. The patient group inclusion criteria were as follows: pathological diagnosis of OSCC and absence of other mucosal lesions in the oral cavity. The patient group exclusion criteria were as follows: antibiotics taken for one month, treatment with radiotherapy and chemotherapy or a history of malignant tumors. Forty patients with OSCC were recruited from August 2020 to February 2021. Additionally, a total of 10 healthy control subjects were recruited. The inclusion criteria were as follows: no oral disease, such as periodontal disease and oral mucosal disease, found through physical examination. The exclusion criteria were as follows: antibiotics taken for one month, treatment with radiotherapy and chemotherapy or a history of malignant tumors. We used disposable sterile nylon flocking swabs for sampling and stored the samples in a prepared oral swab preservation solution (mainly Tris, EDTA and antiseptic) to prevent DNA degradation. The patient sample came from the surface of the tumor. For patients with tongue cancer, the samples were taken at the edge of the tongue. The buccal sampling was from the middle of the buccal mucosa, not the lips or pharynx. Sampling at the bottom of the mouth was performed in front of the sublingual caruncle. For gum cancer, we only selected those patients with cancer in the posterior gum. The palate sampling site was the soft palate. We guaranteed the consistency of the sampling process. All the samples were kept on ice and transported to the laboratory within 2 h after collection and then stored at −80°C in the laboratory until subsequent use.
DNA Extraction, Metagenomic Library Preparation and Sequencing
Total DNA was extracted from the samples using protease K cleavage combined with the phenol chloroform extraction method (Sambrook and Russell, 2006). DNA purity was verified by running the samples on 1.2% agarose gels. The DNA concentration was determined using a Qubit Fluorometer (Thermo Fisher, Waltham, MA, United States). Extracted DNA was sheared on a Covaris M220 (Covaris, Woburn, MA, United States) programmed to generate 300-bp fragments. The sequencing libraries were constructed with a NEBNext® Ultra™ DNA Library Prep Kit for Illumina® (NEB, Ipswich, MA, United States). The products were purified using AgarosAgencourt AMPure XP (Beckman, Brea, CA, United States) and quantified using the GenNext™ NGS Library Quantification Kit (Toyobo, Japan). The libraries were sequenced using an Illumina NovaSeq 6000 and 150-bp paired-end technology.
Bioinformatics Analysis
The raw FASTQ files were demultiplexed based on the index. The raw, paired-end reads were trimmed and quality controlled using Trimmomatic to remove low-quality reads and filter out reads of less than 50 bp (window size of 50 bp; if the average quality was lower than 20, the base after the window was cropped) (Bolger et al., 2014). Cutadapt software was used to remove reads containing N bases (N parameter setting 10) and reads whose overlap with the adapter exceeded 15 bp (Martin, 2011). Default parameters were used with the BWA software to map the human reference genome GRCh38/hg38, and host gene sequences and highly similar contaminated reads were removed (Li and Durbin, 2010). The optimized sequences were assembled using Megahit under default settings (Li et al., 2015). Putative genes were predicted using Prodigal in metagenome mode (-p meta) (Hyatt et al., 2012). We used CD-HIT software (parameters: 95% identity, 90% coverage) to cluster the predicted gene sequences of all the samples (Fu et al., 2012). The longest gene of each category was selected as the representative sequence to construct a non-redundant gene set. Using the genomeCoverageBed method of Bedtools (Quinlan and Hall, 2010), the gene abundance in each sample was calculated and normalized to 100%. Then, by using DIAMOND software (Buchfink et al., 2015), we mapped the gene set with the NCBI NR database using BLASTP with an e-value of 1e-5 and selected the best hit as the species information of the sequence. The taxonomic information included kingdom, phylum, class, order, family, genus, and species.
The β diversity was calculated using Bray–Curtis dissimilarity (R version 3.3.2, vegan package 2.4–4). The P-values were adjusted based on the false discovery rate (FDR) using the Benjamini–Hochberg method (Benjamini and Hochberg, 1995). Permutational multivariate analysis of variance (PERMANOVA) using the “adonis” function in the R Vegan package was performed to assess the effects of phenotype on gene/taxa profiles. To understand the functions of the differentially expressed genes, a functional analysis was performed based on the VFDB database (Chen et al., 2005) and KEGG database (Kanehisa and Goto, 2000), and the mapping ORF was selected using BLAST (BLAST Version 2.2.28 +)1 with an e-value of 1e-5 to calculate the relative abundance of different functions.
For statistical analysis, the abundances of each species and gene were determined to be significantly elevated or depleted in each of the stages using Wilcoxon and Kruskal–Wallis non-parametric tests and by using the ggpubr package in R. Correlation test used Spearman rank correlation analysis. Comparisons between the four groups for participant characteristics were performed with ANOVA or Pearson’s chi-square test or non-parametric test in GraphPad 8.0 (GraphPad Software, San Diego, CA, United States). A p < 0.05 was considered statistically significant.
H&E Staining
After being fixed with 4% paraformaldehyde, the cancer tissue was dehydrated with a concentration gradient from 70 to 100% ethanol, and the dehydration time of each concentration was 15 min. Then, xylene was used to replace the ethanol. The tissue blocks were put into a mold containing melted paraffin wax for embedding. After solidification, the sample was placed at 4°C overnight for sectioning. The section thickness was 5 μm. Then, the sections were deparaffinized in xylene and covered with an alcohol gradient solution for hydration (Fischer et al., 2008). Next, we used hematoxylin for dyeing for 5 min. After this step, the sections were treated with 1% hydrochloric acid alcohol for 2 s and rinsed with water for 5 min. Then, eosin staining was performed for 1 min. After that, the sections were dehydrated in a gradient alcohol and xylene solution. Finally, the sections were sealed with neutral gum for observation.
Results
Characteristics of the Participants and Microbiome Composition Between the Patients and Controls
To reduce system errors, the proportion of samples from each site in each group was similar. In addition, to reduce system errors, we limited the covariate factors that might affect the results. The cancer patients and controls had similar proportions of sex, age, alcohol consumption, smoking status, and BMI. A total of 50 samples were sequenced (Table 1). We used pathological sections of OSCC with H&E staining to demonstrate the DOI (Supplementary Figure 1). The “plumb line” represented the DOI.
Table 1. Clinical characteristics of the subjects.
After quality filtering, assembly, gene prediction and gene set construction, 44,456,926 genes from 50 samples were identified for subsequent analysis (Supplementary Table 1 and Supplementary Figure 2). In the β diversity comparison (comparing the difference between community and community from one angle), the Bray–Curtis analysis showed significant differences in the microbiomes of the patients and healthy controls (PERMANOVA, F = 6.747, R2 = 0.1232, p < 0.001). Using a cluster dendrogram to display and compare the similarity between samples, the similarity of the control group samples was closer than that of the DOI groups (Supplementary Figure 3). A principal coordinate analysis (PCoA) ordination separated the patients and controls. The x-axis (PC1) of microbiota composition explained 22.31% of the variation, and the y-axis (PC2) explained 10.8% of the variation (Figure 1A). The composition of the OSCC microbiome was different from that of the control at the β diversity-based gene level.
Figure 1. Microbial diversity between the patient group and the healthy control group. (A) PCoA was done on Bray-Curtis distances displaying β-diversity of the microbiome between the patients and controls. The x- and y-axes represent the first and second principal coordinates, respectively, and the proportion of variance. (B,C) At the phylum and class levels, the bacteria changed significantly. (D) Comparing the two groups, the cumulative graph shows 24 species with significant changes. Because the results of metagenomics included bacteria, fungi, and archaea, the total abundance of bacteria was less than 100%. The designation s_unannotated indicates that the species information cannot be matched through the NCBI database (using Wilcoxon non-parametric tests, *p < 0.05, **p < 0.01, ***p < 0.001).
We compared the oral microbiome profiles of the OSCC and control groups. The average abundance of phyla in the patient group was lower than that in the control group (total average abundance in the patients was 501,400 vs. 654,092 in the controls). Based on abundance, the phyla Firmicutes, Bacteroidetes, Proteobacteria, Actinobacteria, Fusobacteria, and Spirochetes were the main bacterial groups, which is consistent with previous reports (Cui et al., 2019). These taxa accounted for 97.31% (patients) and 97.88% (controls) of the entire microbiome (Supplementary Table 2). Among them, Spirochetes was significantly higher in the patient group (Figure 1B). Although the abundance of Bacteroidetes and Fusobacteria increased in the OSCC group, the difference was not significant. At the class level, Bacilli, Betaproteobacteria, Actinobacteria and Negativicutes showed clear decreases in the OSCC group (Figure 1C).
At the species level, several species exhibited significantly different abundances in the patient group (Figure 1D and Table 2). Among them, the colorectal cancer-related pathogen Fusobacterium nucleatum (F. nucleatum) exhibited a significant increase in abundance in the patients.
Table 2. Significant changes in species between the two groups.
Stage-Related Microbial Diversity Between the Patients and Healthy Controls
We conducted a stratified analysis between the DOI-high, DOI-medium, DOI-low, and healthy control groups. Compared with the healthy control group, we found microbiome shifts in the DOI-low, DOI-medium and DOI-high groups that were distinct across the stages. Bray–Curtis dissimilarity was used to describe the β diversity. The Adonis test showed that the comparison between the three DOI groups and the control group was significant (PERMANOVA, p < 0.001), but there was no significant difference between the three DOI groups. The PCoA plot showed the result (Figure 2A).
Figure 2. Microbial diversity between the four groups. (A) PCoA was done on Bray-Curtis distances displaying β-diversity of the microbiome between the three DOI groups and the control. (B) Among the four groups, the main phylum changed obviously. Because the results of metagenomics included bacteria, fungi, and archaea, the total abundance of bacteria was less than 100%. (C) Phylum distribution of the number of species that were either elevated or depleted in each of the three stages compared with the healthy controls. (D) The ROC curve shows the significance of bacteria in predicting the grouping of the control and DOI-high groups. (E) The PCoA diagrams were done on the Bray-Curtis matrix based on the genes of annotated species which differed significantly between the four groups. The mapping method was to pick out the genes of the corresponding species, and then the vegdist function (method = Bray) in vegan package (version 2.5.5) was used to calculate the distance based on the abundance of the genes. The pcoa function in ape package (version 5.3) was used for calculation, and then the ggplot function (version 3.3.5) of the ggplot2 package was used for PCoA visualization.
At the phylum level, the abundance of Proteobacteria decreased as the disease progressed (p = 0.046). The abundances of Firmicutes (p = 0.022) and Actinobacteria (p < 0.001) decreased and that of Spirochetes (p = 0.036) increased in the three DOI groups. The three phyla showed similar changes in species taxa, and the Spirochetes increased significantly in the DOI group (Figures 2B,C). In order to better reveal the effect of DOI and sampling site on bacteria, we drew stacked bar graph and boxplot to show the relative abundance of bacteria. Sampling from tongue had no obvious effect on bacterial composition at the phylum taxa, except that Actinobacteria differed significantly (Supplementary Figure 4).
At the class level, similar to the comparison between the control and patient groups, Bacilli (p = 0.007), Betaproteobacteria (p = 0.016) and Actinobacteria (p < 0.001) showed a significant decrease in all three DOI groups, particularly in the DOI-medium and DOI-high groups. The classes Flavobacteriia (p = 0.0145) and Spirochaetia (p = 0.0342) were significantly higher in the three DOI groups (Supplementary Table 3).
We analyzed the changes in the abundance of 290 species, and the abundance of a variety of bacteria showed obvious stage-related changes (Figure 3). We sorted the bacteria in species taxa in the control group and the DOI-high group according to their abundance, and retained the bacteria with significant differences (p < 0.05). Interestingly, the top four abundances in the DOI-high group and the control group changed significantly. There were Rothia mucilaginosa, Prevotella melaninogenica, Haemophilus parainfluenzae and Neisseria flavescens in the healthy control group and Porphyromonas endodontalis, Gemella morbillorum, Gemella haemolysans and Fusobacterium periodonticum in the DOI-high group. These bacterial abundances differed significantly between the control and DOI-high groups. We used receiver operating characteristic (ROC) curves to analyze whether the differences had discriminative significance. The ROC curves for the training set showed a remarkable ability to discriminate the healthy control group by specificity and sensitivity and a reasonable performance in the DOI-high group. The area under the curve (AUC) was 0.992 (95% CI 0.997–1.000) and 0.832 (95% CI 0.723–0.942) (Figure 2D).
Figure 3. In total, 290 differentially abundant species are shown in a phylogenetic tree, grouped in the phyla Proteobacteria, Bacteroidetes, Fusobacteria, Firmicutes and Actinobacteria. In the outer circles, species are marked for relative abundance. The darker the color is, the higher the abundance. Columns were used to display the abundance changes in nine species between groups. (Dunn’s test to compare differences between groups, *p < 0.05, **p < 0.01, ***p < 0.001).
Notably, two species with increased abundance in the patient group, Prevotella intermedia (p = 0.166) and F. nucleatum (p = 0.053), increased with an increasing DOI, but the difference was not significant; however, the high p value may be due to the large standard deviation.
PCoA plots were generated, and ANOSIM was used to test the differences between the groups. The following species differences were identified: Streptococcus pseudopneumoniae (R = 0.142, p = 0.001); R. mucilaginosa (R = 0.197, p = 0.001), Mycobacteroides abscessus (R = 0.201, p = 0.001), N. flavescens (R = 0.138, p = 0.003), P. intermedia (R = 0.106, p = 0.008), G. morbillorum (R = 0.231, p = 0.001), G. haemolysans (R = 0.072, p = 0.044), and F. nucleatum (R = 0.086, p = 0.03). ANOSIM analysis indicated that there were differences among the four groups in β diversity (Figure 2E), especially between the control and DOI-high groups (the yellow and blue ellipse).
Functional Analysis of the Microbiome
We explored the functional features of the oral microbiome across the four groups in our study by annotating the gene catalog based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) modules. In the four-group comparison, as the disease progressed, the expression of some functional genes decreased simultaneously in level 2 (Figure 4A), such as ABC transporters and purine metabolism, particularly in the DOI-medium and DOI-high groups. Some KEGG pathway relative abundances increased as cancer progressed, such as bacterial chemotaxis, bacterial invasion of epithelia and flagellar assembly (Supplementary Figure 5).
Figure 4. Microbiome functional changes between groups. (A) The KEGG database counted the abundance of each knumber, found the knumber contained in the corresponding pathway, summed the abundances of all the knumbers contained in the pathway, and then normalized it to 100%. After that, finding the corresponding level 2 of the pathway, and calculating their separately relative abundance. KEGG pathways that significantly decreased as the disease progressed. (B,C) The VFDB database using mapped ORF in gene set calculated the abundance from the aspects of gene and VF number. It also added the abundance of genes or the same VF number, and then normalized it to 100%. Analysis of significant differences in virulence-related genes and virulence factors among the four groups. Some virulence factors had increased expression with tumor progression. (using Kruskal-Wallis non-parametric tests, *p < 0.05, **p < 0.01, ***p < 0.001).
The abundance of virulence factors was analyzed using the VFDB database. Compared with that in the control group, the abundance of related virulence factors increased in the patient group, mainly including biofilm-related virulence factors such as alginate and inflammation and immune-related virulence factors such as lipopolysaccharide (LPS), hemolysin, and flagella. Between the four groups, these virulence factors showed a trend of increasing abundance as the disease progressed. The analysis of virulence factor genes showed that the higher the DOI was, the higher the abundance of some virulence factor genes, such as hlyb, adeG, flhF, lpxD, and ompA (Figures 4B,C).
The above results showed that a variety of bacteria changed significantly according to the staging. We further carried out Spearman correlation analysis of bacteria and microbiome functions among the four groups. First, we analyzed the correlation between bacteria and KEGG pathways. Bacteria whose abundance decreased as the disease progressed, such as N. flavescens, H. parainfluenzae, R. mucilaginosa, and P. melaninogenica, were positively correlated with metabolism-related pathways, such as DNA replication, mismatch repair, and amino acid metabolism, which explained the overall decrease in metabolism-related potential. F. nucleatum, Capnocytophaga sputigena, P. endodontalis, and G. haemolysans showed obvious correlations with bacterial chemotaxis and flagellar assembly (Figure 5A). These four pathogens were positively correlated with multiple virulence genes, such as algR and ompA. Similarly, analysis of virulence factors showed that they were positively correlated with multiple virulence factors (Figures 5B,C).
Figure 5. Spearman correlation analysis between bacteria and function among the four groups. We selected bacteria that significantly increased or decreased during the progression of OSCC and analyzed the correlation between KEGG pathways, virulence genes, and virulence factors. (A) Correlation analysis between bacteria and KEGG pathways. F. nucleatum, C. sputigena, P. endodontalis, and G. haemolysans showed a significant positive correlation with bacterial chemotaxis and flagellar assembly. (B) Correlation analysis between bacteria and virulence genes. The pathogens F. nucleatum showed a significant positive correlation with the algR, flhF, ompA, lpxD, and ybtP genes. (C) Correlation analysis between bacteria and virulence factors. P. endodontalis and Pse5Ac7Ac, O-antigen showed a significant positive correlation with Pse5Ac7Ac and O-antigen. (using Kruskal-Wallis non-parametric tests,*p < 0.05, **p < 0.01).
Discussion
This study sheds light on the link between the microbiome and the DOI and illustrates the dynamic microbial changes during OSCC progression. The comparison between patients and healthy controls is partially consistent with previous studies, but the sequencing of the metagenomics provides additional information.
Regarding the changes in the microbiome at different tumor stages, we believe that during the progression of the cancer, the tumor microenvironment changes, and the abundance of highly abundant bacteria in the healthy state is reduced due to changes in the living environment, which is reflected in the decreasing expression of metabolism or other functions. For example, the phylum Actinobacteria, which acts as a human commensal in the oropharynx, gastrointestinal tract, and female genital tract, decreased in patients, particularly in the DOI-high group. Studies show that Actinomyces spp. exert a protective effect through the secretion of protease inhibitors that inhibit tumorigenesis (Hozumi et al., 1972). The progression of the tumor led to a low abundance of this taxon, and its tumor-suppressor effect was weakened. Some opportunistic oral pathogens are well suited for survival in the context of cancer, and their abundances increase. The increase in F. nucleatum in patients was statistically significant, and this finding is consistent with previous reports (Zhang et al., 2019). This result is important in the context of what is known regarding F. nucleatum and colorectal cancer (CRC), which is associated with several mechanisms that promote cancer progression, such as adhesion to CDH1 and TLR4 (Rubinstein et al., 2013; Yang Y. et al., 2017). Recent studies have pointed out that F. nucleatum abundance correlates with high glucose metabolism in patients with CRC, which supports carcinogenesis by increasing CRC cell glucose metabolism by elevating ENO1-IT, which acts as a modular guider for KAT7 histone acetyltransferase (Hong et al., 2021). The species F. nucleatum has attracted attention in OSCC; F. nucleatum can enhance the expression of STAT3, JAK1, and MYC under co-culture conditions to promote cell proliferation (Harrandah et al., 2020). In addition, F. nucleatum contributes to oral cancer cell proliferation via the Ku70/p53 pathway (Geng et al., 2020). A study using co-culture of F. nucleatum and human immortalized oral epithelial cells (HOIECs) observed the differential expression of 353 mRNAs, and the expression of multiple oncogenes, such as CREM, CREB1 and NCOA, was upregulated. An analysis of oral cancer samples in a database indicated the same trend (Zhang et al., 2021). In addition, F. nucleatum bound and activated the cell inhibitory receptor CEACAM1 on CEACAM1 + TILs and CEACAM1 + tumor cells, indicating its potential importance in modulating antitumor immunity, which helped the tumor evade immune cell attack through an additive mechanism (Gur et al., 2019). Despite the lack of relevant mechanistic research, a P. intermedia infection can increase cancer risk (Mai et al., 2016). Our research also revealed some trends that contrast with the results from previous research. The abundance of F. periodonticum, a pathogen, gradually decreased as the disease progressed (Yang et al., 2018).
We think the four bacteria found at high abundance in the control group are of great significance. As the disease progressed, their abundance decreased, and significant differences occurred in the DOI-high stage. The ROC analysis showed that eight species had diagnostic and predictive significance. Given this transitional state, we believe that the four species in the control group may play a central role in the community or represent relatively conserved bacterial species. During the development of the disease, the abundances of these species changed, and these changes were clear in the DOI-high group. In the DOI-high group, some bacteria with significantly increased abundance can interact with the host and play a pathogenic role. The association of members of Porphyromonas with colorectal cancers has garnered attention (Ahn et al., 2013; Koliarakis et al., 2019). Our results showed that P. endodontalis, a gram-negative organism considered to be a pulpal pathogen, was an important species in the DOI-high group. Its LPS can induce alveolar bone resorption via the Wnt5a/NF-κB pathway (Mirucki et al., 2014). The species P. endodontalis can reactivate latent Epstein–Barr virus, which is highly correlated with nasopharyngeal carcinoma (Makino et al., 2018). The species G. morbillorum is an aerobic gram-positive coccus and is deemed to be a normal inhabitant of the oral cavity (Chotai et al., 2012). However, when appearing in the intestine, it may be associated with CRC (Kwong et al., 2018). G. haemolysans is an opportunistic pathogen reported to infect immunocompromised patients or cause a poor dental state (Fangous et al., 2016). It is also involved in various infections, mainly endocarditis and eye infections. These three species occupied an important position in the microbiome composition of late-stage OSCC, which was significantly different from that of the control group. However, it is not yet clear whether they produce virulence factors that promote tumor progression.
When analyzing functional changes, the KEGG pathways within the four groups changed, and a considerable number of metabolic pathways declined as the disease progressed. This decline may be related to the bacteria with high abundance and significant changes between groups. Under healthy conditions, various bacteria coexist in a specific way. Some bacteria can produce or break down glycoproteins and then form biofilms, which become the basis for bacterial survival and metabolism (Marsh and Zaura, 2017). Through the KEGG pathways, we found that sugar metabolism, nucleotide metabolism, and amino acid-related metabolism in OSCC patients showed a decreasing trend. These declining functional potential were consistent with the decline in the overall abundance of bacteria at the phylum level. Correlation analysis between bacteria and KEGG confirmed these hypotheses. Metabolic pathways were significantly positively correlated with changes in bacterial abundance, such as in H. parainfluenzae and R. mucilaginosa, whose abundance significantly decreased in the late stage. Some pathogenic bacteria had a strong correlation with pathogenic pathways, such as bacterial chemotaxis and flagellar assembly. Bacterial chemotaxis allows bacteria to grow in suitable locations and is closely related to infectious diseases. For H. pylori, bacterial chemotaxis helps this bacterium overcome the harsh stomach environment and guide it toward its preferred niche, and it plays a role in modulating host immune responses (Johnson and Ottemann, 2018). At present, no studies have confirmed the role of chemotaxis in the occurrence of OSCC.
Some virulence factors and genes were highly expressed. F. nucleatum, C. sputigena, P. endodontalis, and G. haemolysans were positively correlated with multiple virulence genes, for example, P. endodontalis with the hlyB, algR, flhF, ompA, adeG, flmH, lpxD, ybtP genes and F. nucleatum with the algR, flhF, ompA, lpxD, and ybtP genes. Reportedly, to avoid complement-mediated attack, the product of the gene OmpA binds to the C4b-binding protein, which allows microorganisms to survive (Prasadarao et al., 2002). HlyB is a hemolysin-related gene (Zaitseva et al., 2005), and hemolysin can cause inflammation and can promote the secretion of IFN-γ from NK cells, thereby regulating the immune response (Guan et al., 2021). AdeG can encode the AdeFGH efflux pump (He et al., 2015), and algR encodes alginate (Kong et al., 2015), both of which play an important role in the formation of biofilms (Song et al., 2003). Biofilms contribute to the persistence of bacteria in the oral cavity, acting as an adhesin, preventing the bacteria from being expelled and making it more difficult for phagocytes to ingest and kill the bacteria (Mosaddad et al., 2019). LpxD is related to LPS synthesis (Luke et al., 2010; Lin et al., 2012), and LPS plays an important role in inflammation and the regulation of immunity (Kanevskiy et al., 2019). LPS potentially enhances prostate cancer metastasis through NF-κB activation (Jain et al., 2019), and it has also been found to promote breast and lung cancers (Liu et al., 2017; Yang N. et al., 2017). FlmH is a gene related to flagellar synthesis (Zhang et al., 2020). The virulence factors type IV pili and Pse5ac7ac all show increased expression and are closely related to the formation of flagella. Flagella are necessary for motility, adhesion and evasion of the host immune system. For example, H. pylori colonizes the gastric mucosal epithelium using flagella and invades the mucosal epithelium, which is an important basis for malignant transformation (Waskito et al., 2018; Chidwick and Fascione, 2020). There is still a lack of research on virulence factors and carcinogenic mechanisms in OSCC. Although bacteria and virulence factors show obvious correlations, further basic research is needed for verification.
Conclusion
Our research shows that the microbiome of patients with oral cancer changes as the disease progresses, particularly when the tumor progresses to a more advanced stage. Functional analysis shows that various metabolic potential of bacteria exhibit a decreasing trend in patients and are related to disease progression, but virulence factor-related potential show an increasing trend with disease progression. Some pathogenic pathways, virulence genes and factors were obviously correlated with specific bacteria.
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: https://www.ncbi.nlm.nih.gov/, PRJNA735337.
Ethics Statement
The studies involving human participants were reviewed and approved by the Ethics Committee of Shanghai Ninth People’s Hospital. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
YL and ZL contributed to the study conception and drafted and critically revised the manuscript. YQ, XW, and LZ contributed to the data analysis and critically revised the manuscript. All authors gave final approval and agreed to be accountable for all aspects of the work.
Funding
The work was supported by the National Natural Science Foundation of China (No. 81771127).
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.
Publisher’s Note
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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.795777/full#supplementary-material
Footnotes
References
Ahn, J., Sinha, R., Pei, Z., Dominianni, C., Wu, J., Shi, J., et al. (2013). Human gut microbiome and risk for colorectal cancer. J. Natl. Cancer Inst. 105, 1907–1911. doi: 10.1093/jnci/djt300
Amin, M. B., Edge, S. B., Greene, F. L., Byrd, D. R., Brookland, R. K., Washington, M. K., et al. (2017). AJCC Cancer Staging Manual. Berlin: Springer.
Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. Series B Methodol. 57, 289–300.
Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120. doi: 10.1093/bioinformatics/btu170
Brierley, J. D., Gospodarowicz, M. K., and Wittekind, C. (2017). TNM Classification of Malignant Tumours. Hoboken: Wiley.
Buchfink, B., Xie, C., and Huson, D. H. (2015). Fast and sensitive protein alignment using DIAMOND. Nat. Methods 12, 59–60. doi: 10.1038/nmeth.3176
Chen, C., Zhou, Y., Fu, H., Xiong, X., Fang, S., Jiang, H., et al. (2021). Expanded catalog of microbial genes and metagenome-assembled genomes from the pig gut microbiome. Nat. Commun. 12:1106. doi: 10.1038/s41467-021-21295-0
Chen, L., Yang, J., Yu, J., Yao, Z., Sun, L., Shen, Y., et al. (2005). VFDB: a reference database for bacterial virulence factors. Nucleic Acids Res. 33, D325–D328. doi: 10.1093/nar/gki008
Chidwick, H. S., and Fascione, M. A. (2020). Mechanistic and structural studies into the biosynthesis of the bacterial sugar pseudaminic acid (Pse5Ac7Ac). Org. Biomol. Chem. 18, 799–809. doi: 10.1039/c9ob02433f
Chotai, S., Moon, H. J., Kim, J. H., Kim, J. H., Chung, H. S., Park, Y. K., et al. (2012). Brain abscess caused by Gemella morbillorum: case report and review of the literature. Turk. Neurosurg. 22, 374–377. doi: 10.5137/1019-5149.jtn.3634-10.0
Cui, J., Cui, H., Yang, M., Du, S., Li, J., Li, Y., et al. (2019). Tongue coating microbiome as a potential biomarker for gastritis including precancerous cascade. Protein Cell 10, 496–509. doi: 10.1007/s13238-018-0596-6
Dewhirst, F. E., Chen, T., Izard, J., Paster, B. J., Tanner, A. C. R., Yu, W.-H., et al. (2010). The human oral microbiome. J. Bacteriol. 192, 5002–5017. doi: 10.1128/JB.00542-10
Fangous, M. S., Hémon, F., Graf, P., Samier-Guérin, A., Alavi, Z., Le Bars, H., et al. (2016). Bone infections caused by Gemella haemolysans. Med. Mal. Infect. 46, 449–452. doi: 10.1016/j.medmal.2016.06.005
Fischer, A. H., Jacobson, K. A., Rose, J., and Zeller, R. (2008). Hematoxylin and eosin staining of tissue and cell sections. CSH Protoc. 2008:pdb.prot4986. doi: 10.1101/pdb.prot4986
Fu, L., Niu, B., Zhu, Z., Wu, S., and Li, W. (2012). CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinformatics 28, 3150–3152. doi: 10.1093/bioinformatics/bts565
Ganly, I., Goldstein, D., Carlson, D. L., Patel, S. G., O’Sullivan, B., Lee, N., et al. (2013). Long-term regional control and survival in patients with “low-risk,” early stage oral tongue cancer managed by partial glossectomy and neck dissection without postoperative radiation: the importance of tumor thickness. Cancer 119, 1168–1176. doi: 10.1002/cncr.27872
Geng, F., Zhang, Y., Lu, Z., Zhang, S., and Pan, Y. (2020). Fusobacterium nucleatum caused DNA damage and promoted cell proliferation by the Ku70/p53 pathway in oral cancer cells. DNA Cell Biol. 39, 144–151. doi: 10.1089/dna.2019.5064
Guan, Z., Liu, Y., Liu, C., Wang, H., Feng, J., and Yang, G. (2021). Staphylococcus aureus β-Hemolysin Up-Regulates the Expression of IFN-γ by Human CD56(bright) NK Cells. Front. Cell Infect. Microbiol. 11:658141. doi: 10.3389/fcimb.2021.658141
Gur, C., Maalouf, N., Shhadeh, A., Berhani, O., Singer, B. B., Bachrach, G., et al. (2019). Fusobacterium nucleatum supresses anti-tumor immunity by activating CEACAM1. Oncoimmunology 8:e1581531. doi: 10.1080/2162402x.2019.1581531
Harrandah, A. M., Chukkapalli, S. S., Bhattacharyya, I., Progulske-Fox, A., and Chan, E. K. L. (2020). Fusobacteria modulate oral carcinogenesis and promote cancer progression. J. Oral Microbiol. 13:1849493. doi: 10.1080/20002297.2020.1849493
He, X., Lu, F., Yuan, F., Jiang, D., Zhao, P., Zhu, J., et al. (2015). Biofilm formation caused by clinical Acinetobacter baumannii isolates is associated with overexpression of the AdeFGH Efflux Pump. Antimicrob. Agents Chemother. 59, 4817–4825. doi: 10.1128/aac.00877-15
Hong, J., Guo, F., Lu, S. Y., Shen, C., Ma, D., Zhang, X., et al. (2021). F. nucleatum targets lncRNA ENO1-IT1 to promote glycolysis and oncogenesis in colorectal cancer. Gut 70, 2123–2137. doi: 10.1136/gutjnl-2020-322780
Hozumi, M., Ogawa, M., Sugimura, T., Takeuchi, T., and Umezawa, H. (1972). Inhibition of tumorigenesis in mouse skin by leupeptin, a protease inhibitor from Actinomycetes. Cancer Res. 32, 1725–1728.
Hyatt, D., LoCascio, P. F., Hauser, L. J., and Uberbacher, E. C. (2012). Gene and translation initiation site prediction in metagenomic sequences. Bioinformatics 28, 2223–2230. doi: 10.1093/bioinformatics/bts429
IARC Working Group on the Evaluation of Carcinogenic Risks to Humans (2012). Biological agents. Volume 100 B. A review of human carcinogens. IARC Monogr. Eval. Carcinog Risks Hum. 100, 1–441.
Jain, S., Dash, P., Minz, A. P., Satpathi, S., Samal, A. G., Behera, P. K., et al. (2019). Lipopolysaccharide (LPS) enhances prostate cancer metastasis potentially through NF-κB activation and recurrent dexamethasone administration fails to suppress it in vivo. Prostate 79, 168–182. doi: 10.1002/pros.23722
Johnson, K. S., and Ottemann, K. M. (2018). Colonization, localization, and inflammation: the roles of H. pylori chemotaxis in vivo. Curr. Opin. Microbiol. 41, 51–57. doi: 10.1016/j.mib.2017.11.019
Kanehisa, M., and Goto, S. (2000). KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 28, 27–30. doi: 10.1093/nar/28.1.27
Kanevskiy, L. M., Erokhina, S. A., Streltsova, M. A., Ziganshin, R. H., Telford, W. G., Sapozhnikov, A. M., et al. (2019). The role of O-Antigen in LPS-induced activation of human NK cells. J. Immunol. Res. 2019:3062754. doi: 10.1155/2019/3062754
Koliarakis, I., Messaritakis, I., Nikolouzakis, T. K., Hamilos, G., Souglakos, J., and Tsiaoussis, J. (2019). Oral bacteria and intestinal dysbiosis in colorectal cancer. Int. J. Mol. Sci. 20:4146. doi: 10.3390/ijms20174146
Kong, W., Zhao, J., Kang, H., Zhu, M., Zhou, T., Deng, X., et al. (2015). ChIP-seq reveals the global regulator AlgR mediating cyclic di-GMP synthesis in Pseudomonas aeruginosa. Nucleic Acids Res. 43, 8268–8282. doi: 10.1093/nar/gkv747
Kwong, T. N. Y., Wang, X., Nakatsu, G., Chow, T. C., Tipoe, T., Dai, R. Z. W., et al. (2018). Association between bacteremia from specific microbes and subsequent diagnosis of colorectal cancer. Gastroenterology 155, 383–390.e8. doi: 10.1053/j.gastro.2018.04.028
Li, D., Liu, C. M., Luo, R., Sadakane, K., and Lam, T. W. (2015). MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics 31, 1674–1676. doi: 10.1093/bioinformatics/btv033
Li, H., and Durbin, R. (2010). Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics 26, 589–595. doi: 10.1093/bioinformatics/btp698
Lin, L., Tan, B., Pantapalangkoor, P., Ho, T., Baquir, B., Tomaras, A., et al. (2012). Inhibition of LpxC protects mice from resistant Acinetobacter baumannii by modulating inflammation and enhancing phagocytosis. mBio 3, e00312–12. doi: 10.1128/mBio.00312-12
Liu, K. Y., Durham, J. S., Wu, J., Anderson, D. W., Prisman, E., and Poh, C. F. (2016). Nodal disease burden for early-stage oral cancer. JAMA Otolaryngol. Head Neck Surg. 142, 1111–1119. doi: 10.1001/jamaoto.2016.2241
Liu, X., Zhao, W., Wang, W., Lin, S., and Yang, L. (2017). Puerarin suppresses LPS-induced breast cancer cell migration, invasion and adhesion by blockage NF-κB and Erk pathway. Biomed. Pharmacother. 92, 429–436. doi: 10.1016/j.biopha.2017.05.102
Luke, N. R., Sauberan, S. L., Russo, T. A., Beanan, J. M., Olson, R., Loehfelm, T. W., et al. (2010). Identification and characterization of a glycosyltransferase involved in Acinetobacter baumannii lipopolysaccharide core biosynthesis. Infect. Immun. 78, 2017–2023. doi: 10.1128/iai.00016-10
Lydiatt, W. M., Patel, S. G., O’Sullivan, B., Brandwein, M. S., Ridge, J. A., Migliacci, J. C., et al. (2017). Head and Neck cancers-major changes in the American Joint Committee on cancer eighth edition cancer staging manual. CA Cancer J. Clin. 67, 122–137. doi: 10.3322/caac.21389
Mai, X., Genco, R. J., LaMonte, M. J., Hovey, K. M., Freudenheim, J. L., Andrews, C. A., et al. (2016). Periodontal pathogens and risk of incident cancer in postmenopausal females: the buffalo osteoperio study. J. Periodontol. 87, 257–267. doi: 10.1902/jop.2015.150433
Makino, K., Takeichi, O., Imai, K., Inoue, H., Hatori, K., Himi, K., et al. (2018). Porphyromonas endodontalis reactivates latent Epstein-Barr virus. Int. Endod. J. 51, 1410–1419. doi: 10.1111/iej.12959
Marsh, P. D., and Zaura, E. (2017). Dental biofilm: ecological interactions in health and disease. J. Clin. Periodontol. 44, S12–S22. doi: 10.1111/jcpe.12679
Martin, M. (2011). Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 17:3. doi: 10.14806/ej.17.1.200
Meng, C., Bai, C., Brown, T. D., Hood, L. E., and Tian, Q. (2018). Human gut microbiota and gastrointestinal cancer. Genom. Proteom. Bioinform. 16, 33–49. doi: 10.1016/j.gpb.2017.06.002
Mirucki, C. S., Abedi, M., Jiang, J., Zhu, Q., Wang, Y. H., Safavi, K. E., et al. (2014). Biologic activity of porphyromonas endodontalis complex lipids. J. Endod. 40, 1342–1348. doi: 10.1016/j.joen.2014.02.017
Mosaddad, S. A., Tahmasebi, E., Yazdanian, A., Rezvani, M. B., Seifalian, A., Yazdanian, M., et al. (2019). Oral microbial biofilms: an update. Eur. J. Clin. Microbiol. Infect. Dis. 38, 2005–2019. doi: 10.1007/s10096-019-03641-9
Prasadarao, N. V., Blom, A. M., Villoutreix, B. O., and Linsangan, L. C. (2002). A novel interaction of outer membrane protein A with C4b binding protein mediates serum resistance of Escherichia coli K1. J. Immunol. 169, 6352–6360. doi: 10.4049/jimmunol.169.11.6352
Quinlan, A. R., and Hall, I. M. (2010). BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics 26, 841–842. doi: 10.1093/bioinformatics/btq033
Rubinstein, M. R., Wang, X., Liu, W., Hao, Y., Cai, G., and Han, Y. W. (2013). Fusobacterium nucleatum promotes colorectal carcinogenesis by modulating E-cadherin/β-catenin signaling via its FadA adhesin. Cell Host Microbe 14, 195–206. doi: 10.1016/j.chom.2013.07.012
Sambrook, J., and Russell, D. W. (2006). Purification of nucleic acids by extraction with phenol: chloroform. CSH Protoc. 2006:pdb.prot4455.
Schwabe, R. F., and Jobin, C. (2013). The microbiome and cancer. Nat. Rev. Cancer 13, 800–812. doi: 10.1038/nrc3610
Smoot, D. T. (1997). How does Helicobacter pylori cause mucosal damage? Direct mechanisms. Gastroenterology 113, S31–S34.
Song, Z., Wu, H., Ciofu, O., Kong, K. F., Høiby, N., Rygaard, J., et al. (2003). Pseudomonas aeruginosa alginate is refractory to Th1 immune response and impedes host immune clearance in a mouse model of acute lung infection. J. Med. Microbiol. 52, 731–740. doi: 10.1099/jmm.0.05122-0
Tan, W. J., Chia, C. S., Tan, H. K., Soo, K. C., and Iyer, N. G. (2012). Prognostic significance of invasion depth in oral tongue squamous cell carcinoma. ORL J. Otorhinolaryngol. Relat. Spec. 74, 264–270. doi: 10.1159/000343796
Waskito, L. A., Salama, N. R., and Yamaoka, Y. (2018). Pathogenesis of Helicobacter pylori infection. Helicobacter 23:e12516. doi: 10.1111/hel.12516
Yang, C. Y., Yeh, Y. M., Yu, H. Y., Chin, C. Y., Hsu, C. W., Liu, H., et al. (2018). Oral microbiota community dynamics associated with oral squamous cell carcinoma staging. Front. Microbiol. 9:862. doi: 10.3389/fmicb.2018.00862
Yang, N., Liang, Y., Yang, P., and Ji, F. (2017). Propofol suppresses LPS-induced nuclear accumulation of HIF-1α and tumor aggressiveness in non-small cell lung cancer. Oncol. Rep. 37, 2611–2619. doi: 10.3892/or.2017.5514
Yang, Y., Weng, W., Peng, J., Hong, L., Yang, L., Toiyama, Y., et al. (2017). Fusobacterium nucleatum increases proliferation of colorectal cancer cells and tumor development in mice by activating toll-like receptor 4 signaling to nuclear factor-kappaB, and Up-regulating expression of MicroRNA-21. Gastroenterology 152, 851–866.e24. doi: 10.1053/j.gastro.2016.11.018
Zaitseva, J., Jenewein, S., Oswald, C., Jumpertz, T., Holland, I. B., and Schmitt, L. (2005). A molecular understanding of the catalytic cycle of the nucleotide-binding domain of the ABC transporter HlyB. Biochem. Soc. Trans. 33, 990–995. doi: 10.1042/bst20050990
Zhang, K., He, J., Cantalano, C., Guo, Y., Liu, J., and Li, C. (2020). FlhF regulates the number and configuration of periplasmic flagella in Borrelia burgdorferi. Mol. Microbiol. 113, 1122–1139. doi: 10.1111/mmi.14482
Zhang, L., Liu, Y., Zheng, H. J., and Zhang, C. P. (2019). The oral microbiota may have influence on oral cancer. Front. Cell Infect. Microbiol. 9:476. doi: 10.3389/fcimb.2019.00476
Keywords: metagenomic, oral squamous cell carcinoma, depth of invasion, oral microbiome, cancer progression
Citation: Liu Y, Li Z, Qi Y, Wen X and Zhang L (2022) Metagenomic Analysis Reveals a Changing Microbiome Associated With the Depth of Invasion of Oral Squamous Cell Carcinoma. Front. Microbiol. 13:795777. doi: 10.3389/fmicb.2022.795777
Received: 10 December 2021; Accepted: 17 January 2022;
Published: 09 February 2022.
Edited by:
Jozsef Soki, University of Szeged, HungaryReviewed by:
Pourya Gholizadeh, Tabriz University of Medical Sciences, IranDaniel R. Utter, Harvard University, United States
Copyright © 2022 Liu, Li, Qi, Wen and Zhang. 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: Ling Zhang, dG9wZ3VuMTEyOEAxNjMuY29t
†These authors have contributed equally to this work