- 1Department of Medical Oncology, St. Vincent’s Hospital, College of Medicine, The Catholic University of Korea, Seoul, Republic of Korea
- 2Division of Oncology, Department of Medicine, Stanford University School of Medicine, Palo Alto, CA, United States
- 3Department of Electrical Engineering, Stanford University, Palo Alto, CA, United States
Background: Colon microbiome composition contributes to the pathogenesis of colorectal cancer (CRC) and prognosis. We analyzed 16S rRNA sequencing data from tumor samples of patients with metastatic CRC and determined the clinical implications.
Materials and methods: We enrolled 133 patients with metastatic CRC at St. Vincent Hospital in Korea. The V3-V4 regions of the 16S rRNA gene from the tumor DNA were amplified, sequenced on an Illumina MiSeq, and analyzed using the DADA2 package.
Results: After excluding samples that retained <5% of the total reads after merging, 120 samples were analyzed. The median age of patients was 63 years (range, 34–82 years), and 76 patients (63.3%) were male. The primary cancer sites were the right colon (27.5%), left colon (30.8%), and rectum (41.7%). All subjects received 5-fluouracil-based systemic chemotherapy. After removing genera with <1% of the total reads in each patient, 523 genera were identified. Rectal origin, high CEA level (≥10 ng/mL), and presence of lung metastasis showed higher richness. Survival analysis revealed that the presence of Prevotella (p = 0.052), Fusobacterium (p = 0.002), Selenomonas (p<0.001), Fretibacterium (p = 0.001), Porphyromonas (p = 0.007), Peptostreptococcus (p = 0.002), and Leptotrichia (p = 0.003) were associated with short overall survival (OS, <24 months), while the presence of Sphingomonas was associated with long OS (p = 0.070). From the multivariate analysis, the presence of Selenomonas (hazard ratio [HR], 6.35; 95% confidence interval [CI], 2.38–16.97; p<0.001) was associated with poor prognosis along with high CEA level.
Conclusion: Tumor microbiome features may be useful prognostic biomarkers for metastatic CRC.
1 Introduction
The human microbiome plays a fundamental role in maintaining physiological homeostasis. Microbial dysbiosis contributes to the pathophysiology of human diseases, including colorectal cancer (CRC), which has a close physical proximity with the diverse microbial populations in the colon (1, 2). The tumor-associated microbiome plays an important role in colorectal cancer biology given its direct contact with the local tumor microenvironment (TME). DNA sequencing analysis conducted directly on CRC tumors provides insights into the cancer-associated microbiota that is in direct contact with the local TME of CRCs. As an example, many studies investigating the tumor microbiome of CRC have demonstrated an enrichment of Fusobacterium nucleatum (F. nucleatum) in colon carcinogenesis (3–5).
Other studies have delineated how specific microbiome may impact CRC pathophysiology. Preclinical studies have shown that Bacteroides fragilis and Escherichia coli secrete enterotoxins that induce pro-inflammatory immune responses and DNA damage. The direct adhesion of F. nucleatum to colon cells has been shown to exert similar effects. There is mounting evidence supporting the hypothesis that these effects promote colon carcinogenesis (6–10). Identification of specific microbiome alterations may provide clinically relevant information when screening for early-stage CRC (11–16).
Specific microbiome features are also associated with CRC prognosis. For example, the presence of F. nucleatum in CRC tissues is correlated with advanced stage, microsatellite instability, and poor prognosis in patients with CRC who have undergone surgical resection (17, 18). Microbiome composition also influences the efficacy of cancer immunotherapy (19–21) and chemotherapy via various mechanisms (22–26).
Most microbiome studies on CRC have been conducted using tumor samples from surgical resections, which are typically performed with curative intent in lower-stage (i.e., I–III) patients. Therefore, little is known about the clinical implications of the microbiome composition in patients with stage IV metastatic CRC receiving systemic treatment. In this study, we conducted microbiome 16S rRNA sequencing using tumor tissues from patients with initially stage IV CRC and determined the correlation between the tissue microbiome and clinical outcomes.
2 Materials and methods
2.1 Study population
All patients had a histological diagnosis of colorectal adenocarcinoma, initially as a stage IV disease, and received 5-fluouracil (5-FU)-based systemic chemotherapy and/or biological agents (cetuximab or bevacizumab) after tissue sampling. Treatment response was determined using the treatment evaluation criteria for solid tumors (v. 1.1) (27).
This study was approved by the Institutional Review boards (IRB) of the St. Vincent’s Hospital (No. VC19TESI0114), and Stanford Hospital (IRB number 55131). All the patients provided written informed consent to participate in this study. This study was conducted in accordance with the principles of the Declaration of Helsinki.
2.2 Sequencing analysis
The original sample set comprised 133 CRC tissue samples biopsied from patients with stage IV disease. Two scrolls of 7 to 8 µm sections were sliced from the formalin-fixed, paraffin-embedded (FFPE) block and stored at –20°C until their processing. Total genomic DNA was extracted using Maxwell® 16 FFPE Tissue LEV DNA Purification Kit from Promega (Madison, Wisconsin, USA). ZymoBIOMIC Microbial Community DNA Standard (Zymo Research, Catalog No. D6300) was used as a positive control, and DNA-free water was used as a negative control for 16S sequencing and data analysis. For 16S rRNA gene amplicon library preparation, we used 16S metagenomic sequencing library preparation according to the manufacturer’s protocol (Illumina, San Diego, California, USA). Briefly, 200 ng of mucosal DNA was amplified using primers targeting the V3-V4 variable region of the 16S rRNA gene.
16S forward primer:
5’-TCGTCGGCAGCGTCAGATGTGTATAAGAGACAGCCTACG GGNGGCWGCAG-3’;
16S reverse primer:
5’-GTCTCGTGGGCTCG GAGATGTGTAT AAGAGACAGGACTACHVGGGTATCTAATCC-3”.
The DNA sample was adjusted to a final concentration of 12 pM, subjected to PCR amplification and sequence library processing, and sequenced on an Illumina MiSeq platform (Illumina Technologies, USA).
2.3 16S rRNA data analysis
Cutadapt software (version 1.14,-m30) was used to remove the primer and adapter sequences. After filtering low-quality reads and merging, 13 samples had low number of reads and were excluded from the study. This step left 120 samples for further downstream analysis.
The Divisive Amplicon Denoising Algorithm 2 (DADA2) program was used to analyze the data and generate amplicon sequence variant (ASV) tables (Supplementary Figure S1). Then, we converted these tables into phyloseq objects using the “phyloseq” package (v. 1.24.2). Features with ambiguous genus annotations were discarded, and ASVs found in at least 1% of each sample were selected. Microbiomes detected in the negative control were subtracted from the in silico dataset.
Two alpha indices were chosen to depict taxon diversity between the groups: richness determined by number of Observed-ASVs per sample (observed ASV) and evenness determined by diversity of ASVs per sample (Shannon index). The beta diversity to depict intergroup dissimilarities was measured using UniFrac distances and visualized using principal component analysis (PCA). Linear discriminant analysis with effect size (LEfSe) algorithm was used to define differential microbiome patterns (28). The relative microbiome abundance between the groups was visualized using a heat map.
2.4 Statistical analysis
We performed a chi-squared test, Wilcoxon signed-rank test, and Spearman’s correlation to evaluate the results among the different pipelines. The Mann–Whitney U test or one-way ANOVA test was used for alpha diversity comparison. Permutational multivariate analysis of variance (PERMANOVA) was used to compare the beta diversity.
Progression-free survival (PFS) was measured from the date of first-line systemic chemotherapy to the date of progression or censored at the last follow-up date. Overall survival (OS) was measured from the date of first-line systemic chemotherapy until death from any cause, or the last censored date during follow-up. PFS and OS were calculated using the Kaplan–Meier method, and differences in survival between groups were compared using the log-rank test. Variables with p value <0.07 were included in the Cox regression multivariable model, and adjusted hazard ratios (HR) with 95% confidence intervals (CI) were calculated. A p value <0.05 was considered statistically significant. Statistical analyses were performed using R software (version 3.6.1).
3 Results
3.1 Patient characteristics
The baseline characteristics of the 120 enrolled patients are summarized in Table 1. The median age was 63 years (range, 34–82 years), and 76 patients (63.3%) were male. All patients were of Korean ethnicity. The distribution of the primary tumor sites was 27.5% in the right colon, 30.8% in the left colon (excluding the rectum), and 41.7% in the rectum. Twenty-one (18.5%) patients showed an involvement of three metastatic organs or more, which was defined as a high tumor burden.
During the median follow-up period of 17.9 months (1.0–59.0), 51 (42.5%) patients showed responses (complete or partial responses), and the median PFS was 8.9 months (1.0–52.0) after first-line treatment. Twenty-nine (24.2%) patients survived more than 24 months, and the median OS was 18.1 months (3.6–59.4).
3.2 Microbiome composition
In total, 13.3 million reads were generated from 120 CRC tissues. The median read length was approximately 392 bp (range, 283–590 bp). A total of 43,551 ASVs were assigned based on 97% similarity; however, 30,693 (70.5%) remained unclassified at the genus level. We removed unclassified ASV, and filtered out genera that showed an abundance of less than 1% in each patient to remove possible artifacts of sequencing errors or clustering algorithms. Finally, 523 genera were identified and analyzed. Pseudomonas was the most abundant, followed by Clostridium, Bacteroides, Rombroutsia, Flectobacillus, and Staphylococcus, in descending order (Supplementary Figure S2).
3.3 Association of alpha diversity and clinical parameters
Alpha diversity describes the microbial diversity in terms of richness (observed ASVs) and evenness (Shannon index). Alpha diversity analysis revealed that cancer tissues originating from the rectum showed higher microbiome richness than those from left- or right-sided CRC (p = 0.018), but there was no difference in microbiome diversity (Figure 1A). A trend towards higher microbiome richness (p = 0.057), but lower diversity (p = 0.047) was observed in CRC tissues with high serum CEA level (≥ 10 ng/mL, Figure 1B). CRC tissues with peritoneal invasion (cT4) showed lower microbiome diversity than those without peritoneal invasion (p = 0.016, Figure 1C). Patients with lung metastasis showed higher microbiome richness than those without lung metastasis (p = 0.026; Figure 1D). There were no significant differences in alpha diversity according to clinical outcomes, including response rate, progression, and survival.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1310054/fonc-13-1310054-HTML-r1/image_m/fonc-13-1310054-g001.jpg)
Figure 1 Microbial diversity. Alpha-diversity index according to primary site (A), CEA level (B), clinical T stage (C), and presence of lung metastasis (D).
3.4 Association of beta diversity and clinical parameters
Beta diversity differed according to primary site (p = 0.03, Figure 2A), serum CEA level (p = 0.015, Figure 2B), or clinical T stage (p = 0.08, Figure 2C). Patients who survived longer than 24 months (long-term survivors) showed distinct beta diversity compared to those who died within 24 months (short-term survivors, p = 0.021, Figure 2D). Short-term survival was associated with Bacteroides, Prevotella, Fusobacterium, Selenomonas, Fretibacterium, and Porphyromonas. In contrast, Sphingomonas, Corynebacterium, Curvibacter, Enhydrobacter, Paracoccus, Rothia, Caulobacter, and Dermacoccus were associated with long-term survival (Figure 3). The mean relative abundances of Bacteroides, Selenomonas, Fretibacterium, and Peptostreptococcus were higher in short-term survivors, while Rothia, Curvibacter, Corynebacterium, and Sphingomonas were more abundant in long-term survivors (Table 2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1310054/fonc-13-1310054-HTML-r1/image_m/fonc-13-1310054-g002.jpg)
Figure 2 Principle co-ordinates analysis based on the Bray-Curtis distance according to primary site (A), CEA level (B), clinical T stage (C), and overall survival (D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1310054/fonc-13-1310054-HTML-r1/image_m/fonc-13-1310054-g003.jpg)
Figure 3 Survival analysis. Differential microbiome abundance between short term (<24 months) and long term (≥24 months) survivors analyzed by linear discriminate analysis with effect size measurement (A). The listed genera were significantly gathered for their respective groups (p<0.05). Heatmap represents the relative abundance of each genera between two groups (B).
3.5 Prognostic implications of cancer tissue microbiome
The survival analysis revealed that the presence of Prevotella (p = 0.050), Fusobacterium (p = 0.002), Selenomonas (p <0.001), Fretibacterium (p = 0.001), Porphyromonas (p = 0.007), Peptostreptococcus (p = 0.002), and Leptotrichia (p = 0.003) were associated with poor OS, while the presence of Sphingomonas was associated with good OS (p = 0.070) (Figure 4). In a multivariate analysis, the presence of Selenomonas (HR, 6.35; 95% CI, 2.38–16.97; p <0.001) remained a significant poor prognostic factor, along with high serum CEA level (Table 3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1310054/fonc-13-1310054-HTML-r1/image_m/fonc-13-1310054-g004.jpg)
Figure 4 Kaplan-Meier analysis of overall survival according to the presence or absence of specific microbes in the colorectal tissue. Prevotella (A), Fusobacterium (B), Selenomonas (C), Fretibacterium (D), Porphyromonas (E), Peptostreptococcus (F), Leptotrichia (G), and Sphingomonas (H).
4 Discussion
We characterized the composition of the tissue microbiome in patients with metastatic CRC and identified potential microbiome genera that may contribute to the prediction of survival outcomes in this population. Some clinical parameters are associated with microbiome diversity. The rectal origin showed higher microbiome richness and distinct beta diversity, in accordance with previous reports (29, 30). Patients with high CEA levels showed higher microbiome richness, but lower diversity. Peritoneal invasion (cT4) was associated with lower diversity. These are clinical parameters of poor prognosis in CRC; however, we did not find any significant differences in alpha diversity according to the clinical outcomes. Microbiome feature was associated with survival outcomes. The presence of Prevotella, Fusobacterium, Selenomonas, Fretibacterium, Porphyromonas, Peptostreptococcus and Leptotrichia was associated with shorter OS, while Sphingomonas was associated with longer OS.
One strength of our study is that we performed microbiome analysis using tissue samples from patients with metastatic CRC instead of stool samples. The colon tissue microbiome directly reflects features of the TME, and thus is less affected by extrinsic factors such as diet, drug exposure, gut preparation, or sample storage processes than stool samples (31). Studies comparing local tumor tissues with stool microbiome patterns have yielded discordant results (11, 29, 30). Thus, direct determination of the tissue-associated microbiome from primary CRC tumors may prove useful in future studies to confirm its prognostic implications.
Interestingly, the microbiomes related to poor prognosis in our study belonged to the oral cavity. Microbiomes living in the oral cavity can translocate to the lower gut and participate in pathogenesis (32). Significant enrichment of orally originated microbiomes has been observed in CRC tissues compared to normal tissues (12, 16, 29, 33). Fusobacteria, particularly F. nucleatum, is an established part of the oral microbiome and has the potential to contribute to CRC carcinogenesis. Its presence in CRC is associated with prognosis after surgical resection (8, 10, 17). These bacteria induce EMT, dysregulate the E-cadherin/β-catenin complex, and secrete matrix metalloproteinases or chemokines, which enhance tumor aggressiveness and metastasis (34–37). Few studies reported a prognostic role of F. nucleatum in the metastatic setting (38, 39). In the present study, Fusobacterium was associated with poor survival outcomes. At the species level, F. nucleatum was detected in 15 (12.5%) patients, and its presence was significantly associated with shorter OS (p <0.001; data not shown).
The prognostic significance of other bacteria in CRC has not been widely studied (16). We also identified other bacterial genera that may be contributing factors. For example, Porphyromonas promotes cell proliferation or inflammatory responses via interaction with TLR2 or TLR4 in preclinical models, and is associated with the prognosis of lung or esophageal cancer (40, 41). Peptostreptococcus involved in CRC carcinogenesis via direct interaction with integrin α2/β1 activates PI3K-ALK-NF-ĸB signaling and immune modulation (42). Leptotrichia and Selenomonas are associated with high-grade CRC; however, their prognostic significance has not yet been reported (30). We also found that some prognostic parameters correlated with each microbiome. Fusobacterium was detected more frequently in patients with high CEA level (p <0.001) and high tumor burden (p = 0.009). Patients with a high tumor burden showed more frequent Leptotrichia (p = 0.022) and Fretibacterium (p = 0.048). In this study, Sphingomonas was associated with favorable prognosis, which has been reported in breast, lung, and gastric cancer (43–45).
Most patients with CRC received 5-FU based chemotherapy. F. nucleatum may cause 5-FU chemoresistance via mechanisms involving aberrant EMT, autophagy, apoptosis, or tumor cell cycle regulation in CRC (26, 34, 46). Porphyromonas enhances the viability of cancer cells exposed to chemotherapeutic agents, including 5-FU, by modifying apoptosis and cell cycle signaling in esophageal cancer (40). Several studies have reported that short-chain fatty acids (SCFA), including those produced by butyrate-producing bacteria, may be associated with 5-FU chemo-sensitivity in CRC. Yuan et al. reported that antibiotic administration decreased 5-FU efficacy and was associated with a decrease of SCFA-producing bacteria, such as Blautia, Roseburia, Anaerotruncus, Oscillibacter, and Bacteroidetes, which could induce inflammation and abrogate antitumor activity. Another study reported that Roseburia, Dorea, and Anaerostipes in stool samples of patients with locally advanced rectal cancer were overrepresented in 5-FU chemoradiation responders (47). In the present study, these microbiomes were not associated with 5-FU efficacy or clinical outcomes. Porphyromonas or Prevotella, which might belong to the butyrate-producing microbiome (48–50), are associated with poor prognosis.
Our study has some limitations. First, it was conducted at a single institution, and the sample size was relatively small. Due to its retrospective nature, some important information, including prior antibiotic or probiotic exposure, was missing. Second, we used FFPE samples, which may present sequencing quality or contamination issues during processing. To minimize these issues, we filtered out reads of low quality, and ASVs with ambiguous genera annotations or those found in less than 1% of the samples were also discarded. Some microbiome was found in a small number of patients, so results should be interpreted in caution. Third, we investigated genera instead of species because of the low sensitivity of the 16S rRNA sequencing method for classifying species (51). Shotgun metagenomic sequencing was not performed because of high costs and excessive contamination with human DNA in the tissue samples (3, 4). Some studies have suggested that finer taxonomic resolution could impede the clinical availability of microbiome analyses for CRC (37).
In conclusion, we performed comprehensive tumor microbiome sequencing in patients with metastatic CRC and identified distinct tumor microbiome features associated with survival outcomes in this population. The presence of Selenomonas, along with the typical oral carcinogenic microbiome in CRC tissues, may be a useful prognostic biomarker. To confirm their clinical value, these microbiomes should be validated from RNA-ISH method and in stool samples from a large cohort.
Data availability statement
The data presented in the study are deposited in the NIH's Short Read Archive repository, accession number PRJNA808420: https://www.ncbi.nlm.nih.gov/sra/PRJNA808420.
Ethics statement
The studies involving humans were approved by the institutional review boards (IRB) of St. Vincent’s Hospital (No. VC19TESI0114) and Stanford Hospital (IRB number 55131). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
HA: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Resources, Writing – original draft. MP: Formal analysis, Writing – original draft. HL: Formal analysis, Software, Supervision, Validation, Writing – review & editing. BL: Data curation, Methodology, Supervision, Writing – review & editing. DP: Formal analysis, Writing – review & editing. AA: Project administration, Resources, Writing – review & editing. AH: Data curation, Methodology, Supervision, Writing – review & editing. GS: Data curation, Methodology, Supervision, Writing – review & editing. HJ: Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was funded by the Clayville Foundation, Fund for Innovation in Cancer Informatics, Catholic Medical Center Research Foundation in the program year of 2023, and a National Research Foundation of Korea (NRF) grant funded by the Korean government (MSIT) (No. RS-2023-00249672).
Acknowledgments
We would like to thank Editage (www.editage.co.kr) for English language editing.
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.1310054/full#supplementary-material.
References
1. Selber-Hnatiw S, Rukundo B, Ahmadi M, Akoubi H, Al-Bizri H, Aliu AF, et al. Human gut microbiota: toward an ecology of disease. Front Microbiol (2017) 8:1265. doi: 10.3389/fmicb.2017.01265
2. Garrett WS. Cancer and the microbiota. Sci (New York NY) (2015) 348(6230):80–6. doi: 10.1126/science.aaa4972
3. Kostic AD, Gevers D, Pedamallu CS, Michaud M, Duke F, Earl AM, et al. Genomic analysis identifies association of Fusobacterium with colorectal carcinoma. Genome Res (2012) 22(2):292–8. doi: 10.1101/gr.126573.111
4. Castellarin M, Warren RL, Freeman JD, Dreolini L, Krzywinski M, Strauss J, et al. Fusobacterium nucleatum infection is prevalent in human colorectal carcinoma. Genome Res (2012) 22(2):299–306. doi: 10.1101/gr.126516.111
5. Villéger R, Lopès A, Veziant J, Gagnière J, Barnich N, Billard E, et al. Microbial markers in colorectal cancer detection and/or prognosis. World J Gastroenterol (2018) 24(22):2327–47. doi: 10.3748/wjg.v24.i22.2327
6. Wu S, Rhee KJ, Albesiano E, Rabizadeh S, Wu X, Yen HR, et al. A human colonic commensal promotes colon tumorigenesis via activation of T helper type 17 T cell responses. Nat Med (2009) 15(9):1016–22. doi: 10.1038/nm.2015
7. Goodwin AC, Destefano Shields CE, Wu S, Huso DL, Wu X, Murray-Stewart TR, et al. Polyamine catabolism contributes to enterotoxigenic Bacteroides fragilis-induced colon tumorigenesis. Proc Natl Acad Sci USA (2011) 108(37):15354–9. doi: 10.1073/pnas.1010203108
8. Rubinstein MR, Wang X, Liu W, Hao Y, Cai G, Han YW. Fusobacterium nucleatum promotes colorectal carcinogenesis by modulating E-cadherin/β-catenin signaling via its FadA adhesin. Cell Host Microbe (2013) 14(2):195–206. doi: 10.1016/j.chom.2013.07.012
9. Arthur JC, Perez-Chanona E, Mühlbauer M, Tomkovich S, Uronis JM, Fan TJ, et al. Intestinal inflammation targets cancer-inducing activity of the microbiota. Sci (New York NY) (2012) 338(6103):120–3. doi: 10.1126/science.1224820
10. Kostic AD, Chun E, Robertson L, Glickman JN, Gallini CA, Michaud M, et al. Fusobacterium nucleatum potentiates intestinal tumorigenesis and modulates the tumor-immune microenvironment. Cell Host Microbe (2013) 14(2):207–15. doi: 10.1016/j.chom.2013.07.007
11. Zeller G, Tap J, Voigt AY, Sunagawa S, Kultima JR, Costea PI, et al. Potential of fecal microbiota for early-stage detection of colorectal cancer. Mol Syst Biol (2014) 10:766. doi: 10.15252/msb.20145645
12. Thomas AM, Manghi P, Asnicar F, Pasolli E, Armanini F, Zolfo M, et al. Metagenomic analysis of colorectal cancer datasets identifies cross-cohort microbial diagnostic signatures and a link with choline degradation. Nat Med (2019) 25(4):667–78. doi: 10.1038/s41591-019-0405-7
13. Zhu J, Hua RX, Jiang J, Zhao LQ, Sun X, Luan J, et al. Association studies of ERCC1 polymorphisms with lung cancer susceptibility: a systematic review and meta-analysis. PloS One (2014) 9(5):e97616. doi: 10.1371/journal.pone.0097616
14. Amitay EL, Krilaviciute A, Brenner H. Systematic review: Gut microbiota in fecal samples and detection of colorectal neoplasms. Gut Microbes (2018) 9(4):293–307. doi: 10.1080/19490976.2018.1445957
15. Ai L, Tian H, Chen Z, Chen H, Xu J, Fang JY. Systematic evaluation of supervised classifiers for fecal microbiota-based prediction of colorectal cancer. Oncotarget (2017) 8(6):9546–56. doi: 10.18632/oncotarget.14488
16. Xu S, Yin W, Zhang Y, Lv Q, Yang Y, He J. Foes or friends? Bacteria enriched in the tumor microenvironment of colorectal cancer. Cancers (2020) 12(2):372. doi: 10.3390/cancers12020372
17. Mima K, Nishihara R, Qian ZR, Cao Y, Sukawa Y, Nowak JA, et al. Fusobacterium nucleatum in colorectal carcinoma tissue and patient prognosis. Gut (2016) 65(12):1973–80. doi: 10.1136/gutjnl-2015-310101
18. Viljoen KS, Dakshinamurthy A, Goldberg P, Blackburn JM. Quantitative profiling of colorectal cancer-associated bacteria reveals associations between fusobacterium spp., enterotoxigenic Bacteroides fragilis (ETBF) and clinicopathological features of colorectal cancer. PloS One (2015) 10(3):e0119462. doi: 10.1371/journal.pone.0119462
19. Routy B, Le Chatelier E, Derosa L, Duong CPM, Alou MT, Daillère R, et al. Gut microbiome influences efficacy of PD-1-based immunotherapy against epithelial tumors. Sci (New York NY) (2018) 359(6371):91–7. doi: 10.1126/science.aan3706
20. Gopalakrishnan V, Helmink BA, Spencer CN, Reuben A, Wargo JA. 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
21. Sivan A, Corrales L, Hubert N, Williams JB, Aquino-Michaels K, Earley ZM, et al. Commensal Bifidobacterium promotes antitumor immunity and facilitates anti-PD-L1 efficacy. Sci (New York NY) (2015) 350(6264):1084–9. doi: 10.1126/science.aac4255
22. Alexander JL, Wilson ID, Teare J, Marchesi JR, Nicholson JK, Kinross JM. Gut microbiota modulation of chemotherapy efficacy and toxicity. Nat Rev Gastroenterol Hepatol (2017) 14(6):356–65. doi: 10.1038/nrgastro.2017.20
23. Geller LT, Barzily-Rokni M, Danino T, Jonas OH, Shental N, Nejman D, et al. Potential role of intratumor bacteria in mediating tumor resistance to the chemotherapeutic drug gemcitabine. Sci (New York NY) (2017) 357(6356):1156–60. doi: 10.1126/science.aah5043
24. Wallace BD, Roberts AB, Pollet RM, Ingle JD, Biernat KA, Pellock SJ, et al. Structure and inhibition of microbiome β-glucuronidases essential to the alleviation of cancer drug toxicity. Chem Biol (2015) 22(9):1238–49. doi: 10.1016/j.chembiol.2015.08.005
25. Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet CC, Al-Ghalith GA, et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol (2019) 37(8):852–7. doi: 10.1038/s41587-019-0209-9
26. Yu T, Guo F, Yu Y, Sun T, Ma D, Han J, et 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
27. Eisenhauer EA, Therasse P, Bogaerts J, Schwartz LH, Sargent D, Ford R, et al. New response evaluation criteria in solid tumours: revised RECIST guideline (version 1.1). Eur J Cancer (Oxford Engl 1990) (2009) 45(2):228–47. doi: 10.1016/j.ejca.2008.10.026
28. Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, et al. Metagenomic biomarker discovery and explanation. Genome Biol (2011) 12(6):1–18. doi: 10.1186/gb-2011-12-6-r60
29. Flemer B, Lynch DB, Brown JM, Jeffery IB, Ryan FJ, Claesson MJ, et al. Tumour-associated and non-tumour-associated microbiota in colorectal cancer. Gut (2017) 66(4):633–43. doi: 10.1136/gutjnl-2015-309595
30. Zwinsová B, Petrov VA, Hrivňáková M, Smatana S, Micenková L, Kazdová N, et al. Colorectal tumour mucosa microbiome is enriched in oral pathogens and defines three subtypes that correlate with markers of tumour progression. Cancers (2021) 13(19):4799. doi: 10.3390/cancers13194799
31. Yachida S, Mizutani S, Shiroma H, Shiba S, Nakajima T, Sakamoto T, et al. Metagenomic and metabolomic analyses reveal distinct stage-specific phenotypes of the gut microbiota in colorectal cancer. Nat Med (2019) 25(6):968–76. doi: 10.1038/s41591-019-0458-7
32. Komiya Y, Shimomura Y, Higurashi T, Sugi Y, Arimoto J, Umezawa S, et al. Patients with colorectal cancer have identical strains of Fusobacterium nucleatum in their colorectal cancer and oral cavity. Gut (2019) 68(7):1335–7. doi: 10.1136/gutjnl-2018-316661
33. Zhao L, Cho WC, Nicolls MR. Colorectal cancer-associated microbiome patterns and signatures. Front Genet (2021) 12:787176. doi: 10.3389/fgene.2021.787176
34. Ma CT, Luo HS, Gao F, Tang QC, Chen W. Fusobacterium nucleatum promotes the progression of colorectal cancer by interacting with E-cadherin. Oncol Lett (2018) 16(2):2606–12. doi: 10.3892/ol.2018.8947
35. Rubinstein MR, Baik JE, Lagana SM, Han RP, Raab WJ, Sahoo D, et al. Fusobacterium nucleatum promotes colorectal cancer by inducing Wnt/β-catenin modulator Annexin A1. EMBO Rep (2019) 20(4):e47638. doi: 10.15252/embr.201847638
36. Uitto VJ, Baillie D, Wu Q, Gendron R, Grenier D, Putnins EE, et al. Fusobacterium nucleatum increases collagenase 3 production and migration of epithelial cells. Infect Immun (2005) 73(2):1171–9. doi: 10.1128/IAI.73.2.1171-1179.2005
37. Bersanelli M, Santoni M, Ticinesi A, Buti S. The urinary microbiome and anticancer immunotherapy: the potentially hidden role of unculturable microbes. Target Oncol (2019) 14(3):247–52. doi: 10.1007/s11523-019-00643-7
38. Lee DW, Han SW, Kang JK, Bae JM, Kim HP, Won JK, et al. Association between fusobacterium nucleatum, pathway mutation, and patient prognosis in colorectal cancer. Ann Surg Oncol (2018) 25(11):3389–95. doi: 10.1245/s10434-018-6681-5
39. Yan X, Liu L, Li H, Qin H, Sun Z. Clinical significance of Fusobacterium nucleatum, epithelial-mesenchymal transition, and cancer stem cell markers in stage III/IV colorectal cancer patients. OncoTargets Ther (2017) 10:5031–46. doi: 10.2147/OTT.S145949
40. Gao S, Liu Y, Duan X, Liu K, Mohammed M, Gu Z, et al. Porphyromonas gingivalis infection exacerbates oesophageal cancer and promotes resistance to neoadjuvant chemotherapy. Br J Cancer (2021) 125(3):433–44. doi: 10.1038/s41416-021-01419-5
41. Liu Y, Yuan X, Chen K, Zhou F, Yang H, Yang H, et al. Clinical significance and prognostic value of Porphyromonas gingivalis infection in lung cancer. Transl Oncol (2021) 14(1):100972. doi: 10.1016/j.tranon.2020.100972
42. Long X, Wong CC, Tong L, Chu ESH, Ho Szeto C, Go MYY, et al. Peptostreptococcus anaerobius promotes colorectal carcinogenesis and modulates tumour immunity. Nat Microbiol (2019) 4(12):2319–30. doi: 10.1038/s41564-019-0541-3
43. Dong Z, Chen B, Pan H, Wang D, Liu M, Yang Y, et al. Detection of microbial 16S rRNA gene in the serum of patients with gastric cancer. Front Oncol (2019) 9:608. doi: 10.3389/fonc.2019.00608
44. Huang YF, Chen YJ, Fan TC, Chang NC, Chen YJ, Midha MK, et al. Analysis of microbial sequences in plasma cell-free DNA for early-onset breast cancer patients and healthy females. BMC Med Genomics (2018) 11(Suppl 1):16. doi: 10.1186/s12920-018-0329-y
45. Peters BA, Hayes RB, Goparaju C, Reid C, Pass HI, Ahn J. The microbiome in lung cancer tissue and recurrence-free survival. Cancer Epidemiol Biomarkers Prev (2019) 28(4):731–40. doi: 10.1158/1055-9965.EPI-18-0966
46. Ternes D, Karta J, Tsenkova M, Wilmes P, Haan S, Letellier E. Microbiome in colorectal cancer: how to get from meta-omics to mechanism? Trends Microbiol (2020) 28(5):401–23. doi: 10.1016/j.tim.2020.01.001
47. Yi Y, Shen L, Shi W, Xia F, Zhang H, Wang Y, et al. Gut microbiome components predict response to neoadjuvant chemoradiotherapy in patients with locally advanced rectal cancer: A prospective, longitudinal study. Clin Cancer Res an Off J Am Assoc Cancer Res (2021) 27(5):1329–40. doi: 10.1158/1078-0432.CCR-20-3445
48. Wilson MJ, Hall V, Brazier J, Lewis MAO. Evaluation of a phenotypic scheme for identification of the 'butyrate-producing' Peptostreptococcus species. J Med Microbiol (2000) 49(8):747–51. doi: 10.1099/0022-1317-49-8-747
49. Okumura S, Konishi Y, Narukawa M, Sugiura Y, Yoshimoto S, Arai Y, et al. Gut bacteria identified in colorectal cancer patients promote tumourigenesis via butyrate secretion. Nat Commun (2021) 12(1):1–14. doi: 10.1038/s41467-021-25965-x
50. Vital M, Howe AC, Tiedje JM. Revealing the bacterial butyrate synthesis pathways by analyzing (meta)genomic data. mBio (2014) 5(2):e00889. doi: 10.1128/mBio.00889-14
Keywords: tissue microbiome, metastasis, colon, rectum, adenocarcinoma
Citation: An HJ, Partha MA, Lee H, Lau BT, Pavlichin DS, Almeda A, Hooker AC, Shin G and Ji HP (2024) Tumor-associated microbiome features of metastatic colorectal cancer and clinical implications. Front. Oncol. 13:1310054. doi: 10.3389/fonc.2023.1310054
Received: 09 October 2023; Accepted: 20 December 2023;
Published: 18 January 2024.
Edited by:
Antonio Mario Scanu, University of Sassari, ItalyReviewed by:
Rachel Violet Purcell, University of Otago, New ZealandRomain Villéger, University of Poitiers, France
Copyright © 2024 An, Partha, Lee, Lau, Pavlichin, Almeda, Hooker, Shin and Ji. 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: Hanlee P. Ji, genomics_ji@stanford.edu