Abstract
Background: The traditional TNM staging system is often insufficient to differentiate the survival discrepancies of colorectal cancer (CRC) patients at TNM stage I/II. Our study aimed to reclassify stage I/II CRC patients into several subgroups with different prognoses and explore their suitable therapeutic methods.
Methods: Single-cell RNA (scRNA) sequencing data, bulk RNA sequencing data, and clinicopathological information of CRC patients were enrolled from the TCGA and GEO databases. The tumor microenvironment of CRC tissues was accessed by the ESTIMATE algorithm. The prognostic genes were identified by Cox regression analysis. GO and KEGG analyses were conducted in the DAVID database. GSEA analysis was performed for annotation of the correlated gene sets.
Results: We successfully reclassified stage I/II CRC patients into two subgroups and discovered that patients in cluster-2 underwent worse overall survival than those in cluster-1. GSEA analysis showed that immune-associated gene sets were positively enriched in cluster-2. Besides, the differentially expressed genes (DEGs) between cluster-1 and cluster-2 patients also participated in immune-related biological processes and signaling pathways. Moreover, we found that more immune cells infiltrated the microenvironment of cluster-2 patients compared to that of cluster-1 patients, such as Tregs and tumor-associated macrophages. ScRNA sequencing analysis uncovered that most of the enriched immune-associated signaling in cluster-2 patients was mainly attributed to these upregulated immune cells whose infiltration levels were also high in CRC tissues rather than in normal tissues. In addition, we demonstrated that the expression of immune checkpoint genes was significantly higher in cluster-2 patients compared to cluster-1 patients. ScRNA sequencing analysis revealed that the infiltrated CD8+T cells in CRC were naïve T cells and can be activated into effector T cells after immune checkpoint blockade (ICB) treatment.
Conclusion: TNM stage I/II CRC patients can be divided into two subgroups, which have different overall survival rates, tumor microenvironment, and response to ICB therapy.
Introduction
Colorectal cancer (CRC) ranks the second leading cause of tumor-related mortality worldwide (Sung et al., 2021). Although the American Joint Committee on Cancer (AJCC) TNM staging system has been widely applied to predict prognosis and formulate therapeutic strategies for CRC patients, it is insufficient to differentiate the survival discrepancies of TNM stage I/II CRC patients. For instance, the prognosis of part CRC patients at stage II was worse than that at stage III (Hari et al., 2013).
Tumor cells live in a complex microenvironment that is composed of various stromal cells, immune cells, extracellular matrix molecules, and cytokines (Wu and Dai, 2017). Mounting studies have proved that the abnormal tumor microenvironment (TME) plays a critical role in the progression and treatment of cancer (Bruni et al., 2020). For example, inducible co-stimulator-activated CD4+ T cells are triggers of antitumor immunity in early-stage breast cancer (Zhou et al., 2021). Besides, Shi et al. (2022) systematically profiled a single-cell immune signature to assess anti-PD-1 immunotherapy efficacy of early-stage hepatocellular carcinoma. Although differential gene expression of tumor-infiltrating CD33 myeloid cells in advanced-versus early-stage CRC has been reported (Toor et al., 2021), it remains elusive whether there is a discrepant TME among CRC at TNM stage I/II.
The modality of immune checkpoint blockade (ICB) has revolutionized the treatment of advanced solid tumors over the last decade (Burtness et al., 2019; Mok et al., 2019). Recently, several ongoing clinical trials suggested that integrating ICB into the neoadjuvant treatment of early-stage triple-negative breast cancer and non-small cell lung cancer improved patients’ survival without adding substantial toxicity (Gobbini and Giaj Levra, 2018; Schmid et al., 2020). Whereas there is no literature about ICB treatment in TNM stage I/II CRC. Moreover, due to the heterogeneity of tumors, selecting TNM stage I/II CRC patients who are more suitable for ICB therapy can promote personalized therapy and avoid overtreatment.
In the present study, we discovered that TNM stage I/II CRC patients can be reclassified into two subgroups with different overall survival rates which was mainly attributed to the distinct immune microenvironment of tumors. Moreover, we revealed that this TNM stage I/II CRC patients with poor outcomes owned higher expression levels of immune checkpoint genes and were more suitable for ICB treatment.
Materials and methods
Bulk RNA sequencing and ScRNA sequencing analysis
The high-throughput bulk RNA sequencing data and clinicopathological characteristics of CRC patients were downloaded from the TCGA database deposited in the University of California, Santa Cruz (UCSC) Xena browser (https://xenabrowser.net/datapages/) and GEO database (GSE17536 and GSE39582) (Marisa et al., 2013; Chen et al., 2019). The transcription values of genes in the enrolled datasets had been transformed into a normalized count. The single-cell RNA (scRNA) sequencing data of CRC tissues were enrolled from two GEO datasets (GSE146771 and GSE122969) (Kurtulus et al., 2019; Zhang et al., 2020) and analyzed in the public Tumor Immune Single-cell Hub (TISCH) database (http://tisch.comp-genomics.org/home/).
Functional enrichment analysis
Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of genes was conducted in the DAVID database (https://david.ncifcrf.gov). The results of the functional enrichment analysis were visualized through an online tool, OmicShare (http://www.omicshare.com/tools). The GSEA analysis was implemented based on the MSigDB database (http://www.gsea-msigdb.org/gsea/msigdb/index.jsp). The enriched gene sets with a p-value < 0.05 and FDR value <0.25 were identified to be significant ones.
Evaluating the infiltration levels of immune and stromal cells
The Estimation of Stromal and Immune cells in Malignant Tumor tissues using Expression data (ESTIMATE) algorithm was used to assess the portion of immune and stromal cells (Yoshihara et al., 2013). The infiltration levels of specific immune cells were estimated by the CIBERSORT, CIBERSORT-ABS, EPIC, and XCELL algorithms. The immune infiltration analysis was performed with the online tool TIMER2 (http://timer.cistrome.org).
Statistical analysis
Statistical analysis was performed by using GraphPad Prism 8.0 (GraphPad Software, United States) and R software (R 4.1). The Kaplan-Meier (KM) curve with a log rank test was used to compare the significant difference in prognosis between the two groups. The prognostic genes were identified by the univariate Cox regression analysis. The statistical difference between the two groups was analyzed through the Wilcoxon test. p-value < 0.05 was considered statistically significant.
Results
Reclassifying TNM stage I/II CRC patients into two subgroups with different prognosis
We first conducted a univariate Cox regression analysis to screen the prognostic genes among TNM stage I/II CRC patients based on the TCGA CRC cohort and selected 46 prognostic genes with p< 0.01 (Figure 1A). Subsequently, the R package of “Consensus ClusterPlus” was used to test whether these prognostic genes could reclassify TNM stage I/II CRC patients into novel subclusters, and the result indicated that the optimal clustering was two (Figures 1B,C). Based on the unsupervised clustering, this TNM stage I/II CRC patients (N = 153) were well divided into two distinct clusters (Figure 1D). The KM curve analysis showed that stage I/II CRC patients in cluster-2 underwent worse OS than that in cluster-1 (HR = 5.50%, 95%CI: 2.21–13.70, p < 0.001) (Figure 1E). Intriguingly, although CRC patients in cluster-2 underwent better OS than patients at the TNM-IV stage (N = 64) (HR = 2.42%, 95%CI: 1.35–4.33, p = 0.0024) (Figure 1F), there was no significant difference in the OS between CRC patients in cluster-2 and patients at the TNM-III stage (N = 126) (HR = 1.04%, 95%CI: 0.58–1.87, p = 0.89) (Figure 1G).
FIGURE 1
Validating the survival discrepancies of TNM stage I/II CRC patients
Two GEO datasets (GSE39582 and GSE17536) were implemented to validate the survival discrepancies of TNM stage I/II CRC patients. Firstly, based on univariate Cox regression analysis, 27 prognostic genes (p< 0.001) and 103 prognostic genes (p< 0.01) were identified in GSE39582 and GSE17536, respectively (Figure 2A, Supplementary Figure S1A). The proportion of ambiguous clustering analysis suggested that the lowest clusters were two (Figures 2B,C, Supplementary Figure S1B,C). Based on the unsupervised clustering, TNM stage I/II CRC patients in both GSE39582 (N = 302) and GSE17536 (N = 81) were reclassified into two subgroups (Figure 2D, Supplementary Figure S1D). The KM curve analysis showed that no matter in GSE39582 or GSE17536, TNM stage I/II CRC patients in cluster-2 underwent worse OS than that in cluster-1 (Figure 2E, Supplementary Figure S1E). As expected, analysis of the GSE39582 dataset exhibited that although CRC patients in cluster-2 underwent better OS than patients at the TNM-IV stage (N = 59) (HR = 3.81%, 95%CI: 2.19–6.45, p < 0.001) (Figure 2F), there was no significant difference in the OS between CRC patients in cluster-2 and patients at the TNM-III stage (N = 208) (HR = 0.98%, 95%CI: 0.69–1.41, p = 0.94) (Figure 2G). Consistently, based on the GSE17536 dataset, CRC patients in cluster-2 underwent better OS than patients at the TNM-IV stage (N = 39) (HR = 3.30%, 95%CI: 1.83–5.94, p < 0.001) (Supplementary Figure S1F), but there was no significant difference in the OS between CRC patients in cluster-2 and patients at the TNM-III stage (N = 57) (HR = 0.87%, 95%CI: 0.45–1.70, p = 0.67) (Supplementary Figure S1G). Taken together, we demonstrated that TNM stage I/II CRC patients can be reclassified into two novel subgroups with distinct overall survival rates.
FIGURE 2
The different immune regulation systems between TNM stage I/II CRC patients in cluster-1 and cluster-2
Based on the TCGA CRC cohort, we carried out a GSEA analysis to explore the difference between TNM stage I/II CRC patients in cluster-1 and cluster-2, and the results showed that gene sets of HALLMARK_INFLAMMATORY_RESPONSE, HALLMARK_COMPLEMENT, HALLMARK_INTERFERON_GAMMA_RESPONSE, HALLMARK_INTERFERON_ALPHA_RESPONSE, HALLMARK_ALLOGRAFT_REJECTION, HALLMARK_IL6_JAK_STAT3_SIGNALING, HALLMARK_IL2_STAT5_SIGNALING, and HALLMARK_INFLAMMATORY_RESPONSE were positively enriched in CRC patients in cluster-2 (Figure 3A). Subsequently, 2,374 differentially expressed genes (DEGs) with p< 0.001 were identified between CRC patients in cluster-1 and cluster-2 (Figure 3B). Biological process analysis showed that these DEGs were mainly involved in immune regulation, such as innate immune response, T cell co-stimulation, T cell activation, immune effector process, lymphocyte migration, and antigen processing and presentation (Figure 3C). In addition, KEGG pathway analysis exhibited that most DEGs participated in immune-related signaling pathways, such as T cell receptor signaling pathway, natural killer cell-mediated cytotoxicity, PD-L1 expression and PD-1 checkpoint pathway in cancer, and chemokine signaling pathway (Figure 3D).
FIGURE 3
Validating the distinct immune regulation systems between TNM stage I/II CRC patients in cluster-1 and cluster-2
The GSE17536 and GSE39582 datasets were used to validate the distinct immune regulation between TNM stage I/II CRC patients in cluster-1 and cluster-2. GSEA analysis of GSE17536 showed that the immune-associated gene sets were positively enriched in CRC patients in cluster-2 (Figure 4A). Consistently, GSEA analysis of GSE39582 exhibited that the immune-associated gene sets were negatively correlated with CRC patients in cluster-1 (Supplementary Figure S2A). Subsequently, 1,466 and 3,419 DEGs with p< 0.001 between CRC patients in cluster-1 and cluster-2 were identified in GSE17536 and GSE39582, respectively (Figure 4B, Supplementary Figure S2B). Biological process analysis based on GSE17536 and GSE39582 both showed that part DEGs were also associated with immune regulation, such as positive regulation of T cell proliferation and positive regulation of IL-8 production (Figure 4C, Supplementary Figure S2C). Moreover, KEGG pathway analysis demonstrated that part DEGs in GSE17536 and GSE39582 both participated in immune-related signaling pathways, such as TNF signaling pathway, natural killer cell-mediated cytotoxicity, PD-L1 expression and PD-1 checkpoint pathway in cancer, and chemokine receptor interaction (Figure 4D, Supplementary Figure S2D). Therefore, the different immune regulations were partly responsible for the discrepant prognosis of TNM stage I/II CRC patients.
FIGURE 4
The discrepant tumor microenvironments between TNM Stage I/II CRC patients in cluster-1 and cluster-2
To investigate the tumor microenvironments between CRC patients in cluster-1 and cluster-2, we first calculated their microenvironment scores based on the TCGA CRC cohort. The results showed that the microenvironment scores were significantly upregulated in TNM stage I/II CRC patients in cluster-2 compared with that in cluster-1 (Figure 5A). Subsequently, we analyzed the immune scores and stromal scores. Intriguingly, we found the immune scores, rather than stromal scores, were significantly elevated in CRC patients in cluster-2 (Figures 5B,C). Next, we compared the infiltrated levels of immune cells between the two groups based on four algorithms. As shown in Figure 5D, most immune cells were significantly enriched in CRC patients in cluster-2, such as CD8+ T cells, Tregs, resting NK cells, tumor-associated macrophages (TAMs), and resting mast cells.
FIGURE 5
To validate our findings, the GSE17536 dataset was analyzed. As shown in Figures 5E–G, the microenvironment scores, immune scores, and stromal scores were significantly upregulated in TNM stage I/II CRC patients in cluster-2 compared with that in cluster-1. In line with the above results, the infiltrated levels of most immune cells were significantly higher in CRC patients in cluster-2 than those in cluster-1 (Figure 5H). Overall, our results indicated that the different prognoses of TNM stage I/II CRC patients in two subgroups were partly attributed to their distinct tumor microenvironment.
Exploring the dysregulated immune signaling and differentially infiltrated immune cells in CRC through ScRNA sequencing
A single-cell RNA (scRNA) sequencing data enrolled from GSE146771 were analyzed to further explore the dysregulated immune signaling and differentially infiltrated immune cells in CRC. Firstly, 10,468 cells collected from normal adjacent tissue (NATs), CRC tissues, and peripheral blood mononuclear cells (PBMC), were well divided into 15 cell subtypes and the corresponding marker genes were exhibited (Figures 6A,B). Intriguingly, these immune cells, such as Tregs, CD8+T cells, B cells, TAMs, and mast cells, whose infiltration levels were significantly elevated in TNM stage I/II CRC patients of cluster-2 were also enriched in CRC tissues compared to NATs (Figure 6C). Our previous GSEA analysis based on bulk RNA sequencing data exhibited that several immune-associated gene sets were positively enriched in cluster-2. Interestingly, scRNA sequencing analysis revealed that these gene sets were enriched in these immune cells whose infiltration levels were significantly elevated in CRC tissues and in CRC patients of cluster-2, including HALLMARK_INFLAMMATORY_RESPONSE, HALLMARK_COMPLEMENT, HALLMARK_INTERFERON_GAMMA_RESPONSE, HALLMARK_INTERFERON_ALPHA_RESPONSE, HALLMARK_ALLOGRAFT_REJECTION, HALLMARK_IL6_JAK_STAT3_SIGNALING, HALLMARK_IL2_STAT5_SIGNALING, and HALLMARK_INFLAMMATORY_RESPONSE (Figure 6D). Our scRNA sequencing analysis demonstrated that the dysregulated immune signaling was mainly enriched in differentially infiltrated immune cells in CRC patients of cluster-2.
FIGURE 6
TNM stage I/II CRC patients in cluster-2 were more suitable for ICB treatment
Subsequently, we examined the expression of immune checkpoint genes. Analysis of the TCGA CRC cohort and GSE17536 dataset revealed that most immune checkpoint genes were significantly upregulated in TNM stage I/II CRC patients of cluster-2 compared with that of cluster-1 (Figures 7A,B). In addition, the Chi-square test revealed that CRC patients in cluster-2 own more high or low microsatellite instability (MSI-H/L) status and less microsatellite stability (MSS) status compared to those in cluster-1 (p < 0.001) (Figure 7C). Given the high infiltrated levels of CD8+T cells in cluster-2 patients, based on scRNA sequencing data of the GSE122969 dataset, we simulated the changes of CD8+T cells in tumors before and after immune checkpoint blockade (ICB) treatment. Firstly, 5,457 immune cells achieved from the tumor-bearing (MC38 cells) mice before and after anti-PD-1/TIM3 treatment were well divided into nine cell subtypes (Figure 7D). The corresponding marker genes of cells were exhibited in Figure 7E. As shown in Figure 7F, the infiltrated CD8+T cells in CRC tissues were naïve CD8+ T cells (CD8Tn) and after ICB treatment, more central memory CD8+ T cells (CD8Tcm), effective CD8+ T cells (CD8Teff), and exhausted CD8+ T cells (CD8Tex) were enriched in tumor tissues. Taken together, our results uncovered that TNM stage I/II CRC patients in cluster-2 were more suitable for ICB treatment.
FIGURE 7
Discussion
For decades, the AJCC TNM classification system provides a guideline for surgical resection, adjuvant chemotherapy, as well as patient outcomes for a variety of cancers (Locker et al., 2006). However, unusual clinical outcomes are often observed in patients at TNM stage I/II (Nagtegaal et al., 2011; Cadiz et al., 2018). For instance, the prognosis of part CRC patients at stage II was worse than that at stage III (Hari et al., 2013). In this study, we verified that TNM stage I/II CRC patients can be well divided into two novel subgroups with distinct overall survival rates. Besides, there was even no difference in prognosis between CRC patients in cluster-2 and advanced CRC patients. Therefore, the therapeutic strategy for TNM stage I/II CRC patients in cluster-2 should be different from that in cluster-1.
Emerging evidence uncovered that the tumor microenvironment plays a critical role in tumor progression and that the pre-existing antitumor adaptive immune reaction is vital for patient survival (Galon et al., 2006). For example, tumor cells can enhance macrophage-mediated immunosuppression and subsequently suppress CD8+ T cytotoxic function to accelerate metastasis (Zhuang et al., 2020). Similarly, our results revealed that the DEGs between TNM stage I/II CRC patients in cluster-1 and cluster-2 mainly participated in immune-related biological processes and signaling pathways. Subsequently, we discovered that more immune cells infiltrated the tumor tissues of CRC patients in cluster-2 compared with that in cluster-1, such as Treg cells, mast cells, TAMs, CD8+ T cells, and B cells. Treg cells suppress abnormal/excessive immune responses to maintain immune homeostasis (Kumar et al., 2020). Treg cells are often involved in tumor development and progression by inhibiting antitumor immunity (Ohue and Nishikawa, 2019). TAMs are also critical regulators of tumors and are significantly associated with metastasis and drug resistance of cancer cells (Guan et al., 2021; Ma et al., 2021). Recently, the advances in macrophage-based cancer immunotherapy have attracted more and more attention (Baradaran et al., 2022). For example, Wang et al. have constructed an engineering endogenous TAM-targeted biomimetic system to reprogram tumor immunosuppressive microenvironment and enhance chemo-immunotherapy (Wang et al., 2021).
As we know, active CD8+ T cells bind and kill tumor cells by secreting granzymes, perforin, and cathepsin C (Basu et al., 2016). Interestingly, the infiltration levels of CD8+ T cells were also markedly upregulated in TNM stage I/II CRC patients of cluster-2 whose prognosis was poor. ScRNA sequencing technology provided a possibility to deeply analyze the subtypes of various cells which often changed the traditional opinions. For example, oncogenic and tumor-suppressing fibroblasts and macrophages were uncovered in the same tumor tissues (Sebastian et al., 2020; Liang et al., 2021). Based on scRNA sequencing, our study revealed that most infiltrated CD8+ T cells in tumor tissues were exhausted CD8+ T cells that have lost their cytotoxicity. Tumors can induce CD8+ T cell exhaustion and inhibit its activation via expressing immune escape factors, such as PD-L1 (Zhang et al., 2018). Recently, the antibody-based ICB treatment has been applied to improve CD8+ T cells’ priming ability and to establish a durable and efficient antitumor immunity (Borst et al., 2018). ICB does not act on the tumor cell itself but directs membrane ligands or receptors to enhance T cell response (van de Ven and Borst, 2015). We discovered that most immune checkpoint genes were upregulated in TNM stage I/II CRC patients in cluster-2, suggesting that these patients may be more suitable for ICB treatment. To validate our hypothesis, we then simulated the changes of CD8+ T cells in CRC tissues before and after ICB treatment based on scRNA-sequencing. We found that after ICB treatment, more activated CD8 + T cells (CD8Tcm and CD8Teff) infiltrated tumor tissues. Although our study provided theoretical support, whether ICB treatment could improve the prognosis of TNM stage I/II CRC patients in cluster-2 should be further investigated in clinical trials.
Indeed, there were several limitations in our study. First, although multiple independent datasets were enrolled to confirm the correctness of the data, it was better to personally detect these parameters. Second, it is necessary to consider the expense and the testing period about the classification when our findings were applied to clinical practice. Third, whether TNM stage I/II CRC patients in cluster-2 were more suitable for ICB treatment should be further validated in clinical.
In conclusion, based on bulk RNA sequencing and scRNA sequencing, we first reclassified CRC patients at TNM stage I/II into two novel subgroups with different overall survival rates, tumor microenvironment, and response to ICB treatment.
Statements
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 in the article/Supplementary Material.
Author contributions
XL, SW, and YP, conceptualization, supervision, funding acquisition, writing original draft preparation; XL and JQ, methodology, software, validation, formal analysis, investigation, data curation; JN and HS, validation, writing review and editing.
Funding
This work was supported by grants from Jiangsu Provincial Key Research and Development Plan (Grant No. BE2019614), Key Project of Science and Technology Development of Nanjing Medicine (ZKX21042), Elderly Health Research Project of Jiangsu Province (Grant No. LR2021017), Specialized Cohort Research Project of Nanjing Medical University (NMUC2020035).
Acknowledgments
We acknowledge the TCGA and GEO databases for providing platforms and contributors for uploading their meaningful datasets.
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/fgene.2022.948920/full#supplementary-material
References
1
BaradaranA.AsadzadehZ.HemmatN.BaghbanzadehA.ShadbadM. A.KhosraviN.et al (2022). The cross-talk between tumor-associated macrophages and tumor endothelium: Recent advances in macrophage-based cancer immunotherapy. Biomed. Pharmacother.146, 112588. 10.1016/j.biopha.2021.112588
2
BasuR.WhitlockB. M.HussonJ.Le Floc'hA.JinW.Oyler-YanivA.et al (2016). Cytotoxic T cells use mechanical force to potentiate target cell killing. Cell165 (1), 100–110. 10.1016/j.cell.2016.01.021
3
BorstJ.AhrendsT.BąbałaN.MeliefC. J. M.KastenmüllerW. (2018). CD4(+) T cell help in cancer immunology and immunotherapy. Nat. Rev. Immunol.18 (10), 635–647. 10.1038/s41577-018-0044-0
4
BruniD.AngellH. K.GalonJ. (2020). The immune contexture and Immunoscore in cancer prognosis and therapeutic efficacy. Nat. Rev. Cancer20 (11), 662–680. 10.1038/s41568-020-0285-7
5
BurtnessB.HarringtonK. J.GreilR.SoulièresD.TaharaM.de CastroG.Jr.et al (2019). Pembrolizumab alone or with chemotherapy versus cetuximab with chemotherapy for recurrent or metastatic squamous cell carcinoma of the head and neck (KEYNOTE-048): A randomised, open-label, phase 3 study. Lancet394 (10212), 1915–1928. 10.1016/s0140-6736(19)32591-7
6
CadizF.GormazJ. G.BurottoM. (2018). Breast cancer staging: Is TNM ready to evolve?J. Glob. Oncol.4, 1–3. 10.1200/jgo.17.00004
7
ChenM. S.LoY. H.ChenX.WilliamsC. S.DonnellyJ. M.CrissZ. K.et al (2019). Growth factor-independent 1 is a tumor suppressor gene in colorectal cancer. Mol. Cancer Res.17 (3), 697–708. 10.1158/1541-7786.Mcr-18-0666
8
GalonJ.CostesA.Sanchez-CaboF.KirilovskyA.MlecnikB.Lagorce-PagèsC.et al (2006). Type, density, and location of immune cells within human colorectal tumors predict clinical outcome. Science313 (5795), 1960–1964. 10.1126/science.1129139
9
GobbiniE.Giaj LevraM. (2018). Is there a room for immune checkpoint inhibitors in early stage non-small cell lung cancer?J. Thorac. Dis.10 (13), S1427–s1437. 10.21037/jtd.2018.01.81
10
GuanW.LiF.ZhaoZ.ZhangZ.HuJ.ZhangY. (2021). Tumor-associated macrophage promotes the survival of cancer cells upon docetaxel chemotherapy via the CSF1/CSF1R-CXCL12/CXCR4 Axis in castration-resistant prostate cancer. Genes (Basel)12 (5), 773. 10.3390/genes12050773
11
HariD. M.LeungA. M.LeeJ. H.SimM. S.VuongB.ChiuC. G.et al (2013). AJCC cancer staging manual 7th edition criteria for colon cancer: Do the complex modifications improve prognostic assessment?J. Am. Coll. Surg.217 (2), 181–190. 10.1016/j.jamcollsurg.2013.04.018
12
KumarP.SainiS.PrabhakarB. S. (2020). Cancer immunotherapy with check point inhibitor can cause autoimmune adverse events due to loss of Treg homeostasis. Semin. Cancer Biol.64, 29–35. 10.1016/j.semcancer.2019.01.006
13
KurtulusS.MadiA.EscobarG.KlapholzM.NymanJ.ChristianE.et al (2019). Checkpoint blockade immunotherapy induces dynamic changes in PD-1(-)CD8(+) tumor-infiltrating T cells. Immunity50 (1), 181–194. e186. 10.1016/j.immuni.2018.11.014
14
LiangL.YuJ.LiJ.LiN.LiuJ.XiuL.et al (2021). Integration of scRNA-seq and bulk RNA-seq to analyse the heterogeneity of ovarian cancer immune cells and establish a molecular risk model. Front. Oncol.11, 711020. 10.3389/fonc.2021.711020
15
LockerG. Y.HamiltonS.HarrisJ.JessupJ. M.KemenyN.MacdonaldJ. S.et al (2006). ASCO 2006 update of recommendations for the use of tumor markers in gastrointestinal cancer. J. Clin. Oncol.24 (33), 5313–5327. 10.1200/jco.2006.08.2644
16
MaJ.ShayitiF.MaJ.WeiM.HuaT.ZhangR.et al (2021). Tumor-associated macrophage-derived CCL5 promotes chemotherapy resistance and metastasis in prostatic cancer. Cell Biol. Int.45 (10), 2054–2062. 10.1002/cbin.11630
17
MarisaL.de ReynièsA.DuvalA.SelvesJ.GaubM. P.VescovoL.et al (2013). Gene expression classification of colon cancer into molecular subtypes: Characterization, validation, and prognostic value. PLoS Med.10 (5), e1001453. 10.1371/journal.pmed.1001453
18
MokT. S. K.WuY. L.KudabaI.KowalskiD. M.ChoB. C.TurnaH. Z.et al (2019). Pembrolizumab versus chemotherapy for previously untreated, PD-L1-expressing, locally advanced or metastatic non-small-cell lung cancer (KEYNOTE-042): A randomised, open-label, controlled, phase 3 trial. Lancet393 (10183), 1819–1830. 10.1016/s0140-6736(18)32409-7
19
NagtegaalI. D.QuirkeP.SchmollH. J. (2011). Has the new TNM classification for colorectal cancer improved care?Nat. Rev. Clin. Oncol.9 (2), 119–123. 10.1038/nrclinonc.2011.157
20
OhueY.NishikawaH. (2019). Regulatory T (Treg) cells in cancer: Can Treg cells be a new therapeutic target?Cancer Sci.110 (7), 2080–2089. 10.1111/cas.14069
21
SchmidP.CortesJ.PusztaiL.McArthurH.KümmelS.BerghJ.et al (2020). Pembrolizumab for early triple-negative breast cancer. N. Engl. J. Med.382 (9), 810–821. 10.1056/NEJMoa1910549
22
SebastianA.HumN. R.MartinK. A.GilmoreS. F.PeranI.ByersS. W.et al (2020). Single-cell transcriptomic analysis of tumor-derived fibroblasts and normal tissue-resident fibroblasts reveals fibroblast heterogeneity in breast cancer. Cancers (Basel)12 (5), E1307. 10.3390/cancers12051307
23
ShiJ.LiuJ.TuX.LiB.TongZ.WangT.et al (2022). Single-cell immune signature for detecting early-stage HCC and early assessing anti-PD-1 immunotherapy efficacy. J. Immunother. Cancer10 (1), e003133. 10.1136/jitc-2021-003133
24
SungH.FerlayJ.SiegelR. L.LaversanneM.SoerjomataramI.JemalA.et al (2021). Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. Ca. Cancer J. Clin.71 (3), 209–249. 10.3322/caac.21660
25
ToorS. M.TahaR. Z.Sasidharan NairV.SalehR.MurshedK.Abu NadaM.et al (2021). Differential gene expression of tumor-infiltrating CD33(+) myeloid cells in advanced- versus early-stage colorectal cancer. Cancer Immunol. Immunother.70 (3), 803–815. 10.1007/s00262-020-02727-0
26
van de VenK.BorstJ. (2015). Targeting the T-cell co-stimulatory CD27/CD70 pathway in cancer immunotherapy: Rationale and potential. Immunotherapy7 (6), 655–667. 10.2217/imt.15.32
27
WangY.YuJ.LuoZ.ShiQ.LiuG.WuF.et al (2021). Engineering endogenous tumor-associated macrophage-targeted biomimetic nano-RBC to reprogram tumor immunosuppressive microenvironment for enhanced chemo-immunotherapy. Adv. Mat.33 (39), e2103497. 10.1002/adma.202103497
28
WuT.DaiY. (2017). Tumor microenvironment and therapeutic response. Cancer Lett.387, 61–68. 10.1016/j.canlet.2016.01.043
29
YoshiharaK.ShahmoradgoliM.MartínezE.VegesnaR.KimH.Torres-GarciaW.et al (2013). Inferring tumour purity and stromal and immune cell admixture from expression data. Nat. Commun.4, 2612. 10.1038/ncomms3612
30
ZhangL.LiZ.SkrzypczynskaK. M.FangQ.ZhangW.O'BrienS. A.et al (2020). Single-cell analyses inform mechanisms of myeloid-targeted therapies in colon cancer. Cell181 (2), 442–459. e429. 10.1016/j.cell.2020.03.048
31
ZhangY.WuL.LiZ.ZhangW.LuoF.ChuY.et al (2018). Glycocalyx-mimicking nanoparticles improve anti-PD-L1 cancer immunotherapy through reversion of tumor-associated macrophages. Biomacromolecules19 (6), 2098–2108. 10.1021/acs.biomac.8b00305
32
ZhouW.YuM.PanH.QiuW.WangH.QianM.et al (2021). Microwave ablation induces Th1-type immune response with activation of ICOS pathway in early-stage breast cancer. J. Immunother. Cancer9 (4), e002343. 10.1136/jitc-2021-002343
33
ZhuangH.DaiX.ZhangX.MaoZ.HuangH. (2020). Sophoridine suppresses macrophage-mediated immunosuppression through TLR4/IRF3 pathway and subsequently upregulates CD8(+) T cytotoxic function against gastric cancer. Biomed. Pharmacother.121, 109636. 10.1016/j.biopha.2019.109636
Summary
Keywords
colorectal cancer, immune, prognosis, immune checkpoint blockade treatment, single-cell RNA sequencing, TCGA, GEO
Citation
Liu X, Qin J, Nie J, Sun H, Pan Y and Wang S (2022) Reclassifying TNM stage I/II colorectal cancer into two subgroups with different overall survival, tumor microenvironment, and response to immune checkpoint blockade treatment. Front. Genet. 13:948920. doi: 10.3389/fgene.2022.948920
Received
20 May 2022
Accepted
05 September 2022
Published
21 September 2022
Volume
13 - 2022
Edited by
Zheng Liu, National Cancer Center of China, China
Reviewed by
Hao Wang, Second Military Medical University, China
Cheng-Le Zhuang, Tongji University, China
Updates
Copyright
© 2022 Liu, Qin, Nie, Sun, Pan and Wang.
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: Yuqin Pan, panyuqin01@163.com; Shukui Wang, sk_wang@njmu.edu.cn
This article was submitted to Cancer Genetics and Oncogenomics, a section of the journal Frontiers in Genetics
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.