Abstract
Microbial communities significantly inhabit the human body. Evidence shows the interaction between the human microbiome and host cells plays a central role in multiple physiological processes and organ microenvironments. However, the majority of related studies focus on gut microbiota or specific tissues/organs, and the component signature of intratumor microbiota across various cancer types remains unclear. Here, we systematically analyzed the correlation between intratumor microbial signature with survival outcomes, genomic features, and immune profiles across 32 cancer types based on the public databases of Bacteria in Cancer (BIC) and The Cancer Genome Atlas (TCGA). Results showed the relative abundance of microbial taxa in tumors compared to normal tissues was observed as particularly noticeable. Survival analysis found that specific candidate microbial taxa were correlated with prognosis across various cancers. Then, a microbial-based scoring system (MS), which was composed of 64 candidate prognostic microbes, was established. Further analyses showed significant differences in survival status, genomic function, and immune profiles among the distinct MS subgroups. Taken together, this study reveals the diversity and complexity of microbiomes in tumors. Classifying cancer into different subtypes based on intratumor microbial signatures might reasonably reflect genomic characteristics, immune features, and survival status.
Introduction
Various and complex microorganisms inhabit the human body, composing what is called the human microbiota (1, 2). Evidence shows the interaction between the human microbiome and host cells plays a central role in multiple physiological processes and organ microenvironments (3, 4). For instance, the gut microbiota regulates host metabolism and the immune system through putative-specific microbes, metabolites, and toxins (5–7). Certain bacterial components of the human microbiota can drive tumorigenesis and development in various cancers. Numerous evidence found that the bacterium Helicobacter pylori contributes to atrophic gastritis, peptic ulcers, gastric cancers via Wnt/b-catenin pathway, and chronic inflammatory response (8–10). Fusobacterium nucleatum has been found highly abundant in colorectal cancer (CRC) tumors and metastasis tissues than in matched normal tissues, which was associated with poorer prognosis for CRC patients (11–13).
Microbial dysbiosis contributes to tumor susceptibility through complex mechanisms, including inducing tumorigenesis and progression through inflammation, remodeling immune and stromal cells in the tumor microenvironment, and interfering with anticancer drug pharmacodynamics (14–17). An analysis of 16 rRNAs found in stool discovered that the structure and function of gut microbiota in patients with lung cancer were unbalanced, and the imbalance between firmicutes and Bacteroides contributed to tumorigenesis and progression of lung cancer (18, 19). Moreover, the diversity of gut microbiome has been found to positively correlate with the efficacy of immunotherapy in various cancer types (20, 21). Fusobacterium nucleatum is an oral anaerobe that has been found to be prevalent in colorectal cancer and breast cancer, which promoted tumor growth and metastatic progression by attaching tumor-displayed Gal-GalNAc via Fap2 (13, 22, 23).
Although the majority of related studies focus on gut microbiota (24–26), several studies have recently characterized the existence, metabolic activity, and functional importance of intra-tumoral microbiota in various cancers (23, 27–30). Exploring the alteration in the microbial community derived from human tissues and organs will help us better understand the occurrence, progression, and therapeutic approaches for tumors. Bacteria in Cancer (BIC, http://bic.jhlab.tw/ ) reveals a collection of curated, decontaminated tissue-resident microbiota of 32 cancer types based on samples from the TCGA program (31). The microbial signatures of tumor and normal samples from 32 cancer types can be estimated from BIC at different taxonomic levels, which provides an excellent and powerful resource for studying the abundance and alternation of microbial components in various cancers.
The objective of the current study was to investigate the microbiota profile across 32 cancer types. First, we identified the relative abundance of microbiota in tumor tissues compared to normal tissues at various taxonomic levels. Next, we estimated the prognostic value of various microbial compositions. Then, three microbiome-based clusters were determined by NMF clustering analysis based on candidate prognostic microbiota. A microbial-based scoring system (MS) was established by applying the least absolute shrinkage and selection operator (LASSO) regression algorithm, a machine learning-based approach that selected features that were predictive of survival outcomes. Furthermore, the correlation between MS and patients’ survival outcomes, genomic features, and immune profiles was further investigated by combining the BIC microbial profile with the genomic and clinical data from the TCGA cohort.
Materials and methods
Data acquisition from BIC and TCGA
In this study, the microbiota profiles of samples from 10,362 cases (including 9,687 tumor tissues and 675 adjacent normal tissues across 32 cancer types) at the phylum, class, order, family, and genus levels were obtained from the BIC database (http://bic.jhlab.tw/ ); the corresponding genomic and clinical data for patients were obtained from the TCGA dataset (https://portal.gdc.cancer.gov/repository ). Figure 1 displays an overview of the study design, and the information on included tumor types is summarized in Table S1.
Figure 1
Microbiota abundance analysis of tumor and normal tissues
First, we calculated the relative abundance of microbiota between the tumor and normal groups at the phylum, class, order, family, and genus taxonomic levels. Then, we applied Partial least squares discrimination analysis (PLS-DA) to visualize the microbiota profile landscapes between tumor and normal samples by using the package “mixOmics”. Next, we performed a statistical analysis of metagenomic profiles (STAMP) to investigate the overall differences in microbiota profiles between the tumor and normal groups. Moreover, Linear discriminant analysis effect size (LEfSe) analysis was used to identify specific enriched microbial biomarkers for each group. Linear discriminant analysis (LDA) was further applied to evaluate the microbial effects for different groups.
Analysis of global microbiota profiles among 31 cancer types
We first identified the cross-tumor abundant microbiota taxa at the phylum, class, order, family, and genus levels, respectively. Then, we calculated the relative abundance of microbiota across 31 tumor types (except for GBM, which only had normal tissues in the BIC program) at the phylum, class, order, family, and genus taxonomic levels.
Non-negative matrix factorization clustering analysis
Patients with survival data and follow time ≥ 30 days were chosen for survival analyses. First, we investigated the prognostic significance of microbiota by performing a univariate Cox proportional hazards model. Then, NMF was applied to identify distinct microbiome-based clusters based on the abundance of candidate prognostic microbes. The optimal number of clusters and their stability were determined by the consensus clustering algorithm. The R package “NMF” was used to perform the consensus clustering (32).
Generation of a novel microbiome-based classification
To quantify the microbiome features of individual patients, we explored a novel microbial-based scoring system (MS) to investigate the microbiota features of the individual patients. Specifically, candidate prognostic microbiomes were chosen from the Lasso regression algorithm to construct the microbial scoring system. The MS was calculated by the corresponding coefficients of selected microbiota signatures:
MS = Σi Coefficient (microbiota)* Abundance (microbiota)
Where i represents the selected microbial signatures.
Gene set variation analysis and gene set enrichment analysis
GSVA enrichment analysis was performed to investigate the variation in biological processes between different MS subgroups by using “GSVA” R packages (33). The gene sets of “c2.all.v2022.1.Hs.symbols” and “h.all.v2022.1.Hs.symbols” were downloaded from the MSigDB database. Adjusted P value <0.05 was considered statistically significant. GSEA was used to explore the signaling enrichment between different MS subgroups by applying the “Clusterprofile” R package (34). The FDR q < 0.25 and P < 0.05 were considered statistically significant.
Estimation of TME cell infiltration
The ssGSEA algorithm was used to quantify the relative abundance of 29 immune cell types in the TME (35, 36). The relative abundance of each immune cell type in each sample was represented by the enrichment scores that were calculated by the ssGSEA analysis. The CIBERSORT algorithm was applied to analyze the compositions of 22 types of tumor-infiltrating immune cells among different MS subgroups (37).
Significantly mutated genes and tumor mutation burden in different MS subgroups
Using the R package maftools (38), the overall mutation landscape was summarized in patients with high and low MS subgroups in the TCGA cohort. Then, TMB scores based on the TGCA somatic mutation data were calculated to evaluate the mutation status between different MS subgroups.
Statistical analysis
Student’s t-tests were applied to analyze normally distributed variables and the Wilcoxon rank-sum test was performed to evaluate non-normally distributed variables. One-way ANOVA and Kruskal-Wallis tests were used to conduct difference comparisons of more than two groups. Kaplan-Meier survival analysis and Cox proportional hazards model were chosen to investigate the prognostic significance of microbiota and microbiota-based subtypes by applying the survival and survminer packages. All statistical P values were two sided, with p < 0.05 being statistically significant. All data processing was done in R 4.0.5 software.
Results
Differential microbiota signatures in tumor and normal tissues
Overall, we collected and integrated the microbiota profile and clinical characteristics of 10,362 samples from 32 cancer types. A total of 47, 56, 127, 303, and 1,607 microbial taxa were obtained for each sample at the phylum, class, order, family, and genus levels, respectively. First, the average relative abundance of differential microbiota between tumor and normal tissues was calculated (Figures 2A–E), which shows the phylum, class, order, family, and genus levels, respectively. Then, the PLS-DA plot exhibited the microbiota profile landscapes between the tumor and normal samples (Figure 2F). We next explored the differential microbial compositions for each group. Overall, 8, 10, 17, 25, and 66 differential microbial components were found between tumor and normal tissues at the phylum, class, order, family, and genus levels, respectively (Figures 2G).
Figure 2
At the phylum level, Chordata, Firmicutes, Fusobacteria, Bacteroidetes, and Planctomycetes were the main bacteria groups in tumor samples, while Proteobacteria and Actinobacteria were the main groups in normal tissues. At the class level, the microbial composition of Actinobacteria, Alphaproteobacteria, and Betaproteobacteria were significantly abundant in normal tissues, whereas Fusobacteriia, Clostridia, Bacilli, and Bacteroidia were increased in tumor samples. Other differential microbial composition signatures between tumor and normal tissues are shown in Figure 2G and Table S2. The top 30 differential microbial taxa between tumor and normal tissues are shown in Figure 2H.
LEfSe analysis helps to identify tumor- and normal-enriched microbiota
The LDA score of specific microbial taxa in the tumor and normal groups showed the compositional abundances of p_Bacteroidetes, c_Bacteroidetes, o_Nostocalles, f_Prevotellaceae, and g_Ruminococcuss were higher in tumor tissues, while the compositional abundances of p_Actinobacteria, c_Actinobacteria, o_Actinomycetales, f_Moraxellaceae, and g_Acinetobacter were enriched in normal tissues (Figure 3A and Table S3). Cladogram further showed the distinct tumor- and normal-enriched microbiota (Figure 3B). Next, we analyzed the cross-tissue abundant microbiome in 32 cancer types at the phylum, class, order, family, and genus levels, respectively. The average relative abundance of microbiota among 31 cancer types (except for GBM, which only had normal tissues in the BIC program) was further explored, which shown in Figures 3C–G. Moreover, the top 10 most differential microbial taxa across 31 cancer types were shown in Figures 3H–L at the phylum, class, order, family, and genus levels, respectively.
Figure 3
Exploring prognostic microbiota and different clusters mediated by microbiota
We first investigated specific candidate microbial taxa that correlated with survival outcomes by integrating the microbial abundance profile from BIC and survival information from TCGA. The results showed that a total of 182 microbiota were significantly associated with overall survival (OS) (Table S4) and 112 microbiomes were significantly associated with disease-specific survival (DSS) (Table S7). Then, Consensus Clustering analysis of the NMF algorithm was applied to classify patients with qualitatively different clusters based on the abundance of 182 prognostic microbiota (Figure 4A). Three distinct microbiome-based clusters were eventually identified, including 4,592 cases in cluster 1, 1,223 cases in cluster 2, and 3,460 cases in cluster 3. A tSNE plot further showed the performance of microbiome-based clusters to distinguish tumor samples in the TCGA cohort (Figure 4B). A Kaplan-Meier plot revealed a particularly prominent survival advantage in cluster 2, whereas the worst prognosis was found in cluster 3 (Figures 4C, D). The proportion of distinct clusters in 31 tumor types are shown in Figure 4J.
Figure 4
Construction of microbial signature and microbiome-based scoring system
The Lasso regression algorithm based on the 182 prognostic microbiota was performed to find candidate microbial signatures (Figures 4E–G). A total of 64 selected microbial taxa were identified from the Lasso regression algorithm (Table S5). Then, a novel microbiome-based scoring system was constructed to quantify the microbial profiles of individual patients, which we termed as microbial score (MS). Consistent with the NMF clustering analysis, two distinct MS subgroups were found and we named these two subgroups MS-low and -high. Further survival analysis indicated significant prognostic differences between the low- and high-MS subgroups (Figures 4H, I). The performance of 64 selected microbial signatures to classify different MS groups was further analyzed by using PLS-LA (Figure 4K). To better illustrate the association between the established MS with prognosis, an alluvial diagram was applied to visualize the attribute changes of individual patients (Figure 5A). LEfSe analysis was further performed to identify specific enriched microbial biomarkers for high- and low-MS groups (Figures 5B, C, and Table S6).
Figure 5
TME cell infiltration characteristics in distinct MS subgroups
The CIBERSORT algorithm was used to show the differences in the compositions of tumor microenvironment (TME) immune cell types between distinct MS subgroups. As shown in Figure 5D, remarkable differences in immune cell infiltration were observed between the high and low-MS groups, which suggested that intratumor microbiota plays an inevitable role in tumor microenvironment immune profiles. Furthermore, the GSVA algorithm showed significant differences in KEGG pathways, cancer hallmarks, and immune profiles among the distinct MS subgroups (Figures 5E–G). The GSEA algorithm showed nucleotide excision repair, cell cycle, DNA replication, homologous recombination, cytosolic DNA sensing pathway, pyrimidine metabolism, proteasome, spliceosome, and P53 signaling pathway were significantly enriched in the MS-high subgroup (Figures 6A–I).
Figure 6
Landscape of genomic variation and expression of different MS subgroups
We further analyzed the distribution differences of somatic mutation and TMB between low- and high-MS subgroups in the TCGA cohort. The top 20 genes of mutation frequency of the low- and high-MS subtypes are shown in Figure 7A, B, respectively. The top differential mutated genes between the low- and high-MS subtypes are shown Figure 7C. Moreover, a remarkable diversity of tumor mutation burden (TMB) was found between the distinct subgroups (Figure 7D). The abovementioned results indicated the potentially complex interaction between genomic variation and microbial components in cancers, which might be novel regulators for tumorigenesis and progression.
Figure 7
Discussion
Increasing evidence suggests microbiome plays an important role in carcinogenesis and progression of various cancer types (24, 39, 40). However, the pan-cancer microbial heterogeneity and their functions are least understood. Here, we conducted a comprehensive analysis of the intratumor microbiota across 32 cancer types, which could provide a better understanding of microbiota dysbiosis and establish a new foundation for studying host-microbiota interactions in tumorigenesis and malignancy in cancer.
We first characterized the global microbial composition at the phylum, class, order, family, and genus levels of tumor and adjacent normal tissues across 32 cancer types. Consistent with previous studies, the relative abundance of Firmicutes and Fusobacteria significantly enriched tumor tissues, whereas that of Proteobacteria, Actinobacteria, Alphaproteobacteria, and Betaproteobacteria were remarkably predominant in adjacent normal tissues (41–47). The LEfSe analysis further supported that microbial taxa, including p_Firmicutes/c_Clostridia, p_Bacteroidetes/c_Bacteroidetes, and p_Fusobacteria/c_Fusobacteria/o_Fusobacteriales, were tumor-enriched microbes, while p_Actinobacteria/c_Actinobacteria, p_Spirochaetes/c_Spirochaetes, and p_Proteobacteria/c_Alphaproteobacteria were normal-enriched microbes, suggesting that distinct microbial component in tumor and normal samples, and specific microbes might play an essential role in tumorigenesis and development of cancer.
Previous studies reported that certain intratumor microbiota of human tumors could also be a potential biomarker for survival outcomes and chemotherapy/immunotherapy (4, 23, 30, 46, 48, 49). Here, we identified specific microbial components that were associated with prognosis in patients from the TCGA cohort. Consistent with other studies, the relative abundance of Actinobacteria, Firmicutes, Fusobacteria, and Proteobacteria were found to be associated with prognosis in cancer. We further established a scoring system (MS) based on candidate microbial signatures that was somehow predictive of survival outcomes, molecular alternation, and immune profiles.
GSVA showed differential genomic function between the high- and low-MS groups. The high-MS group was remarkedly enriched in cell cycle and mismatch repair signaling pathways, while selenoamino acid metabolism and primary bile acid biosynthesis signaling pathways were found significantly enriched in the low-MS group. It is well established that cell cycle and mismatch repair are involved in cancer cell proliferation, sphere-forming capacity, metastasis, and chemotherapeutic/immunotherapeutic sensitiveness (50–54). The interaction between bile acid and microbiota was found to play an essential role in gastrointestinal inflammation and carcinogenesis (55). Considering that, specific microbiomes and their relative biological functions may be involved in the development and malignancy of tumors.
It has been well known that microbiota mediates the host immune system by a complicit mechanism (56, 57). In accordance with previous studies, we demonstrated that microbiota displayed an important role in tumor immune profiles. Furthermore, our study showed that the TMB encounters significant changes between the distinct MS groups. Furthermore, significant gene mutation diversity was observed between the high- and low-MS groups, which indicated the intratumor microbial component might also exert an effect on the immune profile and genomic heterogeneity of tumors.
Compared with recent studies on microbiota alternation in tumors, our research performed a more comprehensive investigation of microbial characteristics in 32 cancer types. However, several limitations need to be clarified in this study. First, the study was analyzed based on the TCGA dataset, so external verification should be performed based on clinical samples to eliminate the false correlation drawing from bioinformatics data in the future. Second, the current study is more of a observational research paper that focuses on the pan-cancer level. The alteration in the microbial community in specific cancer subtypes should be further investigated. Besides, the biological functions and underlying mechanisms of selected microbiota in this study are warranted for further experimental validation.
Conclusion
In this study, we conducted a comprehensive analysis of the intratumor microbiota in 32 cancer types. Significant differences in the microbial components were found between the tumor and adjacent normal tissues. Several candidate microbial biomarkers were further identified and correlated with tumor prognosis. The potential functions of these microbes in tumors merit further study. Furthermore, we established a microbial-based scoring system, which was significantly related to genomic characteristics, immune features, and survival status of patients in the TCGA cohort. We expect that our research will facilitate a better understanding of the intratumor microbiome and provide a new perspective on the role of the microbiome in tumors.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Ethics statement
The study was approved by the Ethics Institutional Review Board of Sir Run Run Shaw hospital Hospital of Zhejiang University, School of Medicine.
Author contributions
XY designed, performed the study, and wrote the manuscript. GF, YH, and HA contributed to data and statistical analysis. XY and ZJ supervised the study and revised the manuscript. All authors contributed to the article and approved the submitted version.
Acknowledgments
We sincerely acknowledge the TCGA and BIC datasets for providing their platforms and uploading their meaningful data.
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/fonc.2023.1127225/full#supplementary-material
References
1
CostelloEKLauberCLHamadyMFiererNGordonJIKnightRet al. Bacterial community variation in human body habitats across space and time. Science (2009) 326(5960):1694–7. doi: 10.1126/science.1177486
2
SenderRFuchsSMiloR. Revised estimates for the number of human and bacteria cells in the body. PloS Biol (2016) 14(8):e1002533. doi: 10.1371/journal.pbio.1002533
3
MatsonVChervinCSGajewskiTF. Cancer and the microbiome-influence of the commensal microbiota on cancer, immune responses, and immunotherapy. Gastroenterology (2021) 160(2):600–13. doi: 10.1053/j.gastro.2020.11.041
4
LauHCHSungJJYuJ. Gut microbiota: impacts on gastrointestinal cancer immunotherapy. Gut Microbes (2021) 13(1):1–21. doi: 10.1080/19490976.2020.1869504
5
ZhengDLiwinskiTElinavE. Interaction between microbiota and immunity in health and disease. Cell Res (2020) 30(6):492–506. doi: 10.1038/s41422-020-0332-7
6
MartinAMSunEWRogersGBKeatingDJ. The influence of the gut microbiome on host metabolism through the regulation of gut hormone release. Front Physiol (2019) 10:428. doi: 10.3389/fphys.2019.00428
7
BelkaidYHandTW. Role of the microbiota in immunity and inflammation. Cell (2014) 157(1):121–41. doi: 10.1016/j.cell.2014.03.011
8
MossSF. The clinical evidence linking helicobacter pylori to gastric cancer. Cell Mol Gastroenterol Hepatol (2017) 3(2):183–91. doi: 10.1016/j.jcmgh.2016.12.001
9
AmievaMPeekRMJr.Pathobiology of helicobacter pylori-induced gastric cancer. Gastroenterology (2016) 150(1):64–78. doi: 10.1053/j.gastro.2015.09.004
10
CoverTLBlaserMJ. Helicobacter pylori in health and disease. Gastroenterology (2009) 136(6):1863–73. doi: 10.1053/j.gastro.2009.01.073
11
BullmanSPedamalluCSSicinskaEClancyTEZhangXCaiDet al. Analysis of fusobacterium persistence and antibiotic response in colorectal cancer. Science (2017) 358(6369):1443–8. doi: 10.1126/science.aal5240
12
CastellarinMWarrenRLFreemanJDDreoliniLKrzywinskiMStraussJet al. Fusobacterium nucleatum infection is prevalent in human colorectal carcinoma. Genome Res (2012) 22(2):299–306. doi: 10.1101/gr.126516.111
13
MimaKNishiharaRQianZRCaoYSukawaYNowakJAet al. Fusobacterium nucleatum in colorectal carcinoma tissue and patient prognosis. Gut (2016) 65(12):1973–80. doi: 10.1136/gutjnl-2015-310101
14
KnippelRJDrewesJLSearsCL. The cancer microbiome: Recent highlights and knowledge gaps. Cancer Discovery (2021) 11(10):2378–95. doi: 10.1158/2159-8290.CD-21-0324
15
YuTGuoFYuYSunTMaDHanJet al. Fusobacterium nucleatum promotes chemoresistance to colorectal cancer by modulating autophagy. Cell (2017) 170(3):548–63 e16. doi: 10.1016/j.cell.2017.07.008
16
WangHRongXZhaoGZhouYXiaoYMaDet al. The microbial metabolite trimethylamine n-oxide promotes antitumor immunity in triple-negative breast cancer. Cell Metab (2022) 34(4):581–94 e8. doi: 10.1016/j.cmet.2022.02.010
17
FuAYaoBDongTChenYYaoJLiuYet al. Tumor-resident intracellular microbiota promotes metastatic colonization in breast cancer. Cell (2022) 185(8):1356–72 e26. doi: 10.1016/j.cell.2022.02.027
18
ZhengYFangZXueYZhangJZhuJGaoRet al. Specific gut microbiome signature predicts the early-stage lung cancer. Gut Microbes (2020) 11(4):1030–42. doi: 10.1080/19490976.2020.1737487
19
ZhuangHChengLWangYZhangYKZhaoMFLiangGDet al. Dysbiosis of the gut microbiome in lung cancer. Front Cell Infect Microbiol (2019) 9:112. doi: 10.3389/fcimb.2019.00112
20
JinYDongHXiaLYangYZhuYShenYet al. The diversity of gut microbiome is associated with favorable responses to anti-programmed death 1 immunotherapy in Chinese patients with NSCLC. J Thorac Oncol (2019) 14(8):1378–89. doi: 10.1016/j.jtho.2019.04.007
21
PinatoDJHowlettSOttavianiDUrusHPatelAMineoTet al. Association of prior antibiotic treatment with survival and response to immune checkpoint inhibitor therapy in patients with cancer. JAMA Oncol (2019) 5(12):1774–8. doi: 10.1001/jamaoncol.2019.2785
22
AbedJEmgardJEZamirGFarojaMAlmogyGGrenovAet al. Fap2 mediates fusobacterium nucleatum colorectal adenocarcinoma enrichment by binding to tumor-expressed gal-GalNAc. Cell Host Microbe (2016) 20(2):215–25. doi: 10.1016/j.chom.2016.07.006
23
ParhiLAlon-MaimonTSolANejmanDShhadehAFainsod-LeviTet al. Breast cancer colonization by fusobacterium nucleatum accelerates tumor growth and metastatic progression. Nat Commun (2020) 11(1):3259. doi: 10.1038/s41467-020-16967-2
24
HelminkBAKhanMAWHermannAGopalakrishnanVWargoJA. The microbiome, cancer, and cancer therapy. Nat Med (2019) 25(3):377–88. doi: 10.1038/s41591-019-0377-7
25
ZitvogelLMaYRaoultDKroemerGGajewskiTF. The microbiome in cancer immunotherapy: Diagnostic tools and therapeutic strategies. Science (2018) 359(6382):1366–70. doi: 10.1126/science.aar6918
26
GopalakrishnanVHelminkBASpencerCNReubenAWargoJA. The influence of the gut microbiome on cancer, immunity, and cancer immunotherapy. Cancer Cell (2018) 33(4):570–80. doi: 10.1016/j.ccell.2018.03.015
27
PooreGDKopylovaEZhuQCarpenterCFraraccioSWandroSet al. Microbiome analyses of blood and tissues suggest cancer diagnostic approach. Nature (2020) 579(7800):567–74. doi: 10.1038/s41586-020-2095-1
28
NejmanDLivyatanIFuksGGavertNZwangYGellerLTet al. The human tumor microbiome is composed of tumor type-specific intracellular bacteria. Science (2020) 368(6494):973–80. doi: 10.1126/science.aay9189
29
JinCLagoudasGKZhaoCBullmanSBhutkarAHuBet al. Commensal microbiota promote lung cancer development via gammadelta T cells. Cell (2019) 176(5):998–1013 e16.
30
PushalkarSHundeyinMDaleyDZambirinisCPKurzEMishraAet al. The pancreatic cancer microbiome promotes oncogenesis by induction of innate and adaptive immune suppression. Cancer Discovery (2018) 8(4):403–16. doi: 10.1158/2159-8290.CD-17-1134
31
ChenKPHsuCLOyangYJHuangHCJuanHF. BIC: a database for the transcriptional landscape of bacteria in cancer. Nucleic Acids Res (2022). doi: 10.1093/nar/gkac891
32
GaujouxRSeoigheC. A flexible r package for nonnegative matrix factorization. BMC Bioinf (2010) 11:367. doi: 10.1186/1471-2105-11-367
33
HanzelmannSCasteloRGuinneyJ. GSVA: Gene set variation analysis for microarray and RNA-seq data. BMC Bioinf (2013) 14:7. doi: 10.1186/1471-2105-14-7
34
SubramanianATamayoPMoothaVKMukherjeeSEbertBLGilletteMAet al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U.S.A. (2005) 102(43):15545–50. doi: 10.1073/pnas.0506580102
35
CharoentongPFinotelloFAngelovaMMayerCEfremovaMRiederDet al. Pan-cancer immunogenomic analyses reveal genotype-immunophenotype relationships and predictors of response to checkpoint blockade. Cell Rep (2017) 18(1):248–62. doi: 10.1016/j.celrep.2016.12.019
36
BarbieDATamayoPBoehmJSKimSYMoodySEDunnIFet al. Systematic RNA interference reveals that oncogenic KRAS-driven cancers require TBK1. Nature (2009) 462(7269):108–12. doi: 10.1038/nature08460
37
NewmanAMLiuCLGreenMRGentlesAJFengWXuYet al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods (2015) 12(5):453–7. doi: 10.1038/nmeth.3337
38
MayakondaALinDCAssenovYPlassCKoefflerH P. Maftools: efficient and comprehensive analysis of somatic variants in cancer. Genome Res (2018) 28(11):1747–56. doi: 10.1101/gr.239244.118
39
Sepich-PooreGDZitvogelLStraussmanRHastyJWargoJAKnightRet al. The microbiome and human cancer. Science (2021) 371(6536). doi: 10.1126/science.abc4552
40
CullinNAzevedo AntunesCStraussmanRStraussmanRStein-ThoeringerCKElinavEet al. Microbiome and cancer. Cancer Cell (2021) 39(10):1317–41. doi: 10.1016/j.ccell.2021.08.006
41
FlemerBLynchDBBrownJMJefferyIBRyanFJClaessonMJet al. Tumour-associated and non-tumour-associated microbiota in colorectal cancer. Gut (2017) 66(4):633–43. doi: 10.1136/gutjnl-2015-309595
42
LamontRJFitzsimondsZRWangHGaoS. Role of porphyromonas gingivalis in oral and orodigestive squamous cell carcinoma. Periodontol 2000 (2022) 89(1):154–65. doi: 10.1111/prd.12425
43
HozakaYOiHSatakeSUchinoYGotoYIdichiTet al. Are intratumoral microbiota involved in the progression of intraductal papillary mucinous neoplasms of the pancreas? Surgery (2022). doi: 10.1016/j.surg.2022.10.003
44
ZhaoQYangTYanYZhangYLiZWangYet al. Alterations of oral microbiota in Chinese patients with esophageal cancer. Front Cell Infect Microbiol (2020) 10:541144. doi: 10.3389/fcimb.2020.541144
45
GreathouseKLWhiteJRVargasAJBliskovskyVVBeckJAVon MuhlinenNet al. Interaction between the microbiome and TP53 in human lung cancer. Genome Biol (2018) 19(1):123. doi: 10.1186/s13059-018-1501-6
46
YamamuraKBabaYNakagawaSMimaKMiyakeKNakamuraKet al. Human microbiome fusobacterium nucleatum in esophageal cancer tissue is associated with prognosis. Clin Cancer Res (2016) 22(22):5574–81. doi: 10.1158/1078-0432.CCR-16-1786
47
FengQLiangSJiaHStadlmayrATangLLanZet al. Gut microbiome development along the colorectal adenoma-carcinoma sequence. Nat Commun (2015) 6:6528. doi: 10.1038/ncomms7528
48
TzengASangwanNJiaMLiuCCKeslarKSDowns-KellyEet al. Human breast microbiome correlates with prognostic features and immunological signatures in breast cancer. Genome Med (2021) 13(1):60. doi: 10.1186/s13073-021-00874-2
49
PengZChengSKouYWangZJinRHuHet al. The gut microbiome is associated with clinical response to anti-PD-1/PD-L1 immunotherapy in gastrointestinal cancer. Cancer Immunol Res (2020) 8(10):1251–61. doi: 10.1158/2326-6066.CIR-19-1014
50
LiuJPengYWeiW. Cell cycle on the crossroad of tumorigenesis and cancer therapy. Trends Cell Biol (2022) 32(1):30–44. doi: 10.1016/j.tcb.2021.07.001
51
ZhaoPLiLJiangXLiQ. Mismatch repair deficiency/microsatellite instability-high as a predictor for anti-PD-1/PD-L1 immunotherapy efficacy. J Hematol Oncol (2019) 12(1):54. doi: 10.1186/s13045-019-0738-1
52
BarettiMLeDT. DNA Mismatch repair in cancer. Pharmacol Ther (2018) 189:45–62. doi: 10.1016/j.pharmthera.2018.04.004
53
BolandCRGoelA. Microsatellite instability in colorectal cancer. Gastroenterology (2010) 138(6):2073–87 e3. doi: 10.1053/j.gastro.2009.12.064
54
EvanGIVousdenKH. Proliferation, cell cycle and apoptosis in cancer. Nature (2001) 411(6835):342–8. doi: 10.1038/35077213
55
JiaWXieGJiaW. Bile acid-microbiota crosstalk in gastrointestinal inflammation and carcinogenesis. Nat Rev Gastroenterol Hepatol (2018) 15(2):111–28. doi: 10.1038/nrgastro.2017.119
56
ZhouCBZhouYLFangJY. Gut microbiota in cancer immune response and immunotherapy. Trends Cancer (2021) 7(7):647–60. doi: 10.1016/j.trecan.2021.01.010
57
Wong-RolleAWeiHKZhaoCJinC. Unexpected guests in the tumor microenvironment: microbiome in cancer. Protein Cell (2021) 12(5):426–35. doi: 10.1007/s13238-020-00813-8
Summary
Keywords
microbiota signature, survival outcomes, genomic features, immune profiles, bacteria in cancer, the Cancer Genome Atlas
Citation
Yang X, An H, He Y, Fu G and Jiang Z (2023) Comprehensive analysis of microbiota signature across 32 cancer types. Front. Oncol. 13:1127225. doi: 10.3389/fonc.2023.1127225
Received
19 December 2022
Accepted
17 February 2023
Published
08 March 2023
Volume
13 - 2023
Edited by
Tezer Kutluk, Hacettepe University, Türkiye
Reviewed by
Li Zhang, University of Minnesota Twin Cities, United States; Fuyou Zhou, Anyang Tumor Hospital, China
Updates
Copyright
© 2023 Yang, An, He, Fu and Jiang.
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: Zhinong Jiang, 3200039@zju.edu.cn
This article was submitted to Cancer Epidemiology and Prevention, a section of the journal Frontiers in Oncology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.