Correlations Between Tumor Mutation Burden and Immunocyte Infiltration and Their Prognostic Value in Colon Cancer

Identification and validation of a novel six-gene signature based on mucinous adenocarcinoma-related gene molecular typing in colorectal cancer

BACKGROUND AND OBJECTIVES

Colorectal mucinous adenocarcinoma (MAC) is a particular pathological type that has yet to be thoroughly studied. This study aims to investigate the characteristics of colorectal MAC-related genes in colorectal cancer (CRC), explore the role of MAC-related genes in accurately classifying CRC, and further construct a prognostic signature.

METHODS

CRC samples were collected from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO). MAC-related differentially expressed genes (DEGs) were analyzed in TCGA samples. Based on colorectal MAC-related genes, TCGA CRC samples were molecularly typed by the non-negative matrix factorization (NMF). According to the molecular subtype characteristics, the RiskScore signature was constructed through univariate Cox, the least absolute shrinkage and selection operator (LASSO), and multivariate Cox regression analyses. Clinical significance in CRC of the RiskScore signature was analyzed. A nomogram was further built based on the RiskScore signature.

RESULTS

From the colorectal MAC-related genes, three distinct molecular subtypes were identified. A RiskScore signature composed of six CRC subtype-related genes (CALB1, MMP1, HOXC6, ZIC2, SFTA2, and HYAL1) was constructed. Patients with high-RiskScores had the worse prognoses. RiskScores led to differences in gene mutation characteristics, antitumor drug sensitivity, and tumor microenvironment of CRC. A nomogram based on the signature was developed to predict the one-, three-, and five-year survival of CRC patients.

CONCLUSION

MAC-related genes were able to classify CRC. A RiskScore signature based on the colorectal MAC-related molecular subtype was constructed, which had important clinical significance for guiding the accurate stratification of CRC patients.

PD-1 and PD-L1 inhibitors in cold colorectal cancer: challenges and strategies

Colorectal cancer (CRC) is the second most common cause of cancer mortality, with mismatch repair proficient (pMMR) and/or microsatellite stable (MSS) CRC making up more than 80% of metastatic CRC. Programmed death-ligand 1 (PD-L1) and programmed death 1 (PD-1) immune checkpoint inhibitors (ICIs) are approved as monotherapy in many cancers including a subset of advanced or metastatic colorectal cancer (CRC) with deficiency in mismatch repair (dMMR) and/or high microsatellite instability (MSI-H). However, proficient mismatch repair and microsatellite stable (pMMR/MSS) cold CRCs have not shown clinical response to ICIs alone. To potentiate the anti-tumor response of PD-L1/PD-1 inhibitors in patients with MSS cold cancer, combination strategies currently being investigated include dual ICI, and PD-L1/PD-1 inhibitors in combination with chemotherapy, radiotherapy, vascular endothelial growth factor (VEGF) /VEGF receptor (VEGFR) inhibitors, mitogen-activated protein kinase (MEK) inhibitors, and signal transducer and activation of transcription 3 (STAT3) inhibitors. This paper will review the mechanisms of PD-1/PD-L1 ICI resistance in pMMR/MSS CRC and potential combination strategies to overcome this resistance, summarize the published clinical experience with different combination therapies, and make recommendations for future avenues of research.

Mitochondrial-related genes markers that predict survival in patients with head and neck squamous cell carcinoma affect immunomodulation through hypoxia, glycolysis, and angiogenesis pathways

Mitochondria play a crucial role in the occurrence and development of tumors. We used mitochondria-related genes for consistent clustering to identify three stable molecular subtypes of head and neck squamous cell carcinoma (HNSCC) with different prognoses, mutations, and immune characteristics. Significant differences were observed in clinical characteristics, immune microenvironment, immune cell infiltration, and immune cell scores. TP53 was the most significantly mutated; cell cycle-related pathways and tumorigenesis-related pathways were activated in different subtypes. Risk modeling was conducted using a multifactor stepwise regression method, and nine genes were identified as mitochondria-related genes affecting prognosis (DKK1, EFNB2, ITGA5, AREG, EPHX3, CHGB, P4HA1, CCND1, and JCHAIN). Risk score calculations revealed significant differences in prognosis, immune cell scores, immune cell infiltration, and responses to conventional chemotherapy drugs. Glycolysis, angiogenesis, hypoxia, and tumor-related pathways were positively correlated with the RiskScore. Clinical samples were subjected to qPCR to validate the results. In this work, we constructed a prognostic model based on the mitochondrial correlation score, which well reflects the risk and positive factors for the prognosis of patients with HNSCC. This model can be used to guide individualized adjuvant and immunotherapy in patients with HNSCC.

Construction and Validation of a Novel Prognosis Model in Colon Cancer Based on Cuproptosis-Related Long Non-Coding RNAs

Colon cancer (CC) is one of the most common (6%) malignancies and leading cause of cancer-associated death (more than 0.5 million) worldwide, which demands reliable prognostic biomarkers. Cuproptosis is a novel modality of regulated cell death triggered by the accumulation of intracellular copper. LncRNAs have been reported as prognostic signatures in different types of tumors. However, the correlation between cuproptosis-related lncRNAs (CRLs) and CC remains unclear. Data of CC patients were downloaded from public databases. The prognosis-associated CRLs were identified by co-expression analysis and univariate Cox. Least absolute shrinkage and selection operator were utilized to construct the CRLs-based prognostic signature in silico for CC patients. CRLs level was validated in human CC cell lines and patient tissues. ROC curve and Kaplan-Meier curve results revealed that high CRLs-risk score was associated with poor prognosis in CC patients. Moreover, the nomogram revealed that this model possessed a steady prognostic prediction capability with C-index as 0.68. More importantly, CC patients with high CRLs-risk score were more sensitive to eight targeted therapy drugs. The prognostic prediction power of the CRLs-risk score was further confirmed by cell lines, tissues and two independent CC cohorts. This study constructed a novel ten-CRLs-based prognosis model for CC patients. The CRLs-risk score is expected to serve as a promising prognostic biomarker and predict targeted therapy response in CC patients.

Predictive Models for Colon Adenocarcinoma Diagnosis, Prognosis, and Immune Microenvironment Based on 2 Hypoxia-Related Genes: KDM3A and ENO3

Background: Hypoxia is known to play a critical role in tumor occurrence, progression, prognosis, and therapy resistance. However, few studies have investigated hypoxia markers for diagnosing and predicting prognosis in colon adenocarcinoma (COAD). This study aims to identify a hypoxia genes-based biomarker for predicting COAD patients' prognosis and response to immunotherapy on an individual basis. Methods: Hypoxia-related genes were extracted from the Molecular Signatures Database. Gene expression, clinical data, and mutation data of COAD were collected retrospectively from the Cancer Genome Atlas, the Gene Expression Omnibus, and the International Cancer Genome Consortium databases. Univariate and multivariate cox regression, and the least absolute shrinkage and selection operator method were used to select the genes most associated with the prognosis of COAD patients. Kaplan-Meier survival analysis, receiver operating characteristic curves, calibration curves, and decision curve analyses were performed to validate the efficacy of the signature in predicting the prognosis of COAD patients. EdU incorporation assays, cell survival assays, western blot assays, and trans-well invasion assays were performed to further confirm the function of the screened genes in tumorigenesis. Results: ENO3 and KDM3A were identified as key genes for constructing prognostic and diagnostic signatures, which were found to be independent risk factors for predicting the prognosis and diagnosis of COAD patients. Using these signatures, COAD patients could be stratified into high-risk and low-risk groups, with the latter exhibiting better overall survival outcomes. Moreover, the high-risk group displayed elevated levels of immune checkpoint genes and tumor mutation burden, indicating that these patients may benefit from immune checkpoint inhibitor therapy. Conclusion: The signature developed in this study demonstrates excellent efficacy in prognosticating the outcomes of COAD patients. Moreover, it can serve as a valuable tool for clinicians to identify COAD patients who are suitable for ICI therapy.

A Novel m7G-Related Gene Signature Predicts the Prognosis of Colon Cancer

Colon cancer (CC), one of the most common malignancies worldwide, lacks an effective prognostic prediction biomarker. N7-methylguanosine (m7G) methylation is a common RNA modification type and has been proven to influence tumorigenesis. However, the correlation between m7G-related genes and CC remains unclear. The gene expression levels and clinical information of CC patients were downloaded from public databases. Twenty-nine m7G-related genes were obtained from the published literature. Via unsupervised clustering based on the expression levels of m7G-related genes, CC patients were divided into three m7G clusters. Based on differentially expressed genes (DEGs) from the above three groups, CC patients were further divided into three gene clusters. The m7G score, a prognostic model, was established using principal component analysis (PCA) based on 15 prognosis-associated m7G genes. KM curve analysis demonstrated that the overall survival rate was remarkably higher in the high-m7G score group, which was much more significant in advanced CC patients as confirmed by subgroup analysis. Correlation analysis indicated that the m7G score was associated with tumor mutational burden (TMB), PD-L1 expression, immune infiltration, and drug sensitivity. The expression level of prognosis-related m7G genes was further confirmed in human CC cell lines and samples. This study established an m7G gene-based prognostic model (m7G score), which demonstrated the important roles of m7G-related genes during CC initiation and progression. The m7G score could be a practical biomarker to predict immunotherapy response and prognosis in CC patients.

Prognostic Impact of TP53 Mutations and Tumor Mutational Load in Colorectal Cancer

The DNA damage response (DDR) is critical for maintaining genome stability, and abnormal DDR—resulting from mutations in DNA damage-sensing and repair proteins—is a hallmark of cancer. Here, we aimed to investigate the predictive power of DDR gene mutations and the tumor mutational load (TML) for survival outcomes in a cohort of 22 rectal cancer patients who received pre-operative neoadjuvant therapy. Univariate analysis revealed that TML-high and TP53 mutations were significantly associated with worse overall survival (OS) with TML-high retaining significance in multivariate analyses. Kaplan–Meier survival analyses further showed TML-high was associated with worse disease-free (p = 0.036) and OS (p = 0.024) results in our patient cohort. A total of 53 somatic mutations were identified in 22 samples with eight (36%) containing mutations in DDR genes, including ATM, ATR, CHEK2, MRE11A, RAD50, NBN, ERCC2 and TP53. TP53 was the most frequently mutated gene, and TP53 mutations were significantly associated with worse OS (p = 0.023) in Kaplan–Meier survival analyses. Thus, our data indicate that TML and TP53 mutations have prognostic value for rectal cancer patients and may be important independent biomarkers for patient management. This suggests that prognostic determination for rectal cancer patients receiving pre-operative neoadjuvant therapy should include consideration of the initial TML and tumor genetic status.

Genomic signature of MTOR could be an immunogenicity marker in human colorectal cancer

Background

The mTOR signaling pathway plays an important role in cancer. As a master regulator, the status of MTOR affects pathway activity and the efficacy of mTOR inhibitor therapy. However, little research has been performed to explore MTOR in colorectal cancer (CRC).

Methods

In this study, gene expression and clinical data were analyzed using The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) databases. Signaling pathways related to MTOR in CRC were identified by Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and gene set enrichment analysis (GSEA). Somatic mutation data were downloaded from TCGA and analyzed using the maftools R package. Tumor Immune Estimation Resource (TIMER) and CIBERSORT were used to analyze correlations between MTOR and tumor-infiltrating immune cells (TIICs). Finally, we detected MTOR mutations in a CRC cohort from our database using whole-exome sequencing.

Results

We found that MTOR was overexpressed in Asian CRC patients and associated with a poor prognosis. Enrichment analysis showed that MTOR was involved in metabolism, cell adhesion, and translation pathways in CRC. High MTOR expression was correlated with high tumor mutation burden (TMB) and several TIICs. Finally, we found that the mTOR signaling pathway was activated in CRC lines characterized by microsatellite instability (MSI), and the frequency of MTOR mutations was higher in MSI-high (MSI-H) patients than in microsatellite stable (MSS) patients.

Conclusions

MTOR may represent a comprehensive indicator of prognosis and immunological status in CRC. The genomic signatures of MTOR may provide guidance for exploring the role of mTOR inhibitors in CRC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-09901-w.

Establishment of a circular RNA regulatory stemness-related gene pair signature for predicting prognosis and therapeutic response in colorectal cancer

Background

Colorectal cancer (CRC) is a common malignant tumor of the digestive tract with a poor prognosis. Cancer stem cells (CSCs) affect disease outcomes and treatment responses in CRC. We developed a circular RNA (circRNA) regulatory stemness-related gene pair (CRSRGP) signature to predict CRC patient prognosis and treatment effects.

Methods

The circRNA, miRNA, and mRNA expression profiles and clinical information of CRC patients were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. CRSRGPs were established based on stemness-related genes in the competing endogenous RNA (ceRNA) network. A CRSRGP signature was generated using the least absolute shrinkage and selection operator (Lasso) and Cox regression analysis of TCGA training set. The prognosis was predicted by generating a nomogram integrating the CRSRGP signature and clinicopathologic features. The model was validated in an external validation set (GSE17536). The antitumor drug sensitivity and immunotherapy responses of CRC patients in the high-risk group (HRG) and low-risk group (LRG) were evaluated by the pRRophetic algorithm and immune checkpoint analysis.

Results

We established an 18-CRSRGP signature to predict the prognosis and treatment responses of CRC patients. In the training and external validation sets, risk scores were used to categorize CRC patients into the HRG and LRG. The Kaplan–Meier analysis showed a poor prognosis for patients in the HRG and that subgroups with different clinical characteristics had significantly different prognoses. A multivariate Cox analysis revealed that the CRSRGP signature was an independent prognostic factor. The nomogram integrating clinical features and the CRSRGP signature efficiently predicted CRC patient prognosis, outperformed the current TNM staging system, and had improved practical clinical value. Anticancer drug sensitivity predictions revealed that the tumors of patients in the HRG were more sensitive to pazopanib, sunitinib, gemcitabine, lapatinib, and cyclopamine. Analysis of immune checkpoint markers demonstrated that patients in the HRG were more likely to benefit from immunotherapy.

Conclusion

An efficient, reliable tool for evaluating CRC patient prognosis and treatment response was established based on the 18-CRSRGP signature and nomogram.

The long non-coding RNA T cell leukemia homeobox 1 neighbor enhances signal transducer and activator of transcription 5A phosphorylation to promote colon cancer cell invasion, migration, and metastasis

Colon cancer is among the most prevalent gastrointestinal tumor types. The long noncoding RNA (lncRNA) T cell leukemia homeobox 1 neighbor (TLX1NB) is up-regulated in colorectal cancer (CRC). However, the functional role of this lncRNA in colon cancer remains unknown. In our study, we investigated the clinical significance of TLX1NB in colon cancer through bioinformatics analysis and explored its role in migration, invasion and metastasis of colon cancer cell with a series of experiments. Firstly, TLX1NB was up-regulated in colon cancer tissues and increased TLX1NB expression was significantly associated with advanced N stages. In wound healing assays and transwell assays, TLX1NB overexpression promoted HCT116 cell migration and invasion while TLX1NB knockdown inhibited SW620 cell migration and invasion. In vivo, TLX1NB knockdown suppressed pulmonary metastasis of SW620 cell and vimentin expression but increased E-cadherin expression. Then, TLX1NB overexpression enhanced signal transducer and activator of transcription 5A (STAT5A) phosphorylation and TLX1NB knockdown suppressed STAT5A phosphorylation. Moreover, the inhibition of STAT5A phosphorylation reversed TLX1NB overexpression-associated increase in HCT116 cell migratory and invasive activity. In conclusion, TLX1NB enhances STAT5A phosphorylation to promote colon cancer cell invasion, migration, and metastasis.

KRAS: A Druggable Target in Colon Cancer Patients

Mutations in KRAS are among the most frequent aberrations in cancer, including colon cancer. KRAS direct targeting is daunting due to KRAS protein resistance to small molecule inhibition. Moreover, its elevated affinity to cellular guanosine triphosphate (GTP) has made the design of specific drugs challenging. Indeed, KRAS was considered ‘undruggable’. KRASG12C is the most commonly mutated variant of KRAS in non-small cell lung cancer. Currently, the achievements obtained with covalent inhibitors of this variant have given the possibility to assess the best therapeutic approach to KRAS-driven tumors. Mutation-related biochemical assets and the tissue of origin are expected to influence responses to treatment. Further attempts to obtain mutant-specific KRAS (KRASG12C) switch-II covalent inhibitors are ongoing and the results are promising. Drugs targeted to block KRAS effector pathways could be combined with direct KRAS inhibitors, immunotherapy or T cell-targeting approaches in KRAS-mutant tumors. The development of valuable combination regimens will be essential against potential mechanisms of resistance that may arise during treatment.

A Novel Glycolysis and Hypoxia Combined Gene Signature Predicts the Prognosis and Affects Immune Infiltration of Patients with Colon Cancer

Purpose

We aimed to characterize the expression patterns of glycolysis and hypoxia genes in colon cancers as well as their value in prognosis and immune microenvironment.

Methods

The expression profiles were acquired from the Cancer Genome Atlas database. Enrichment of hypoxia and glycolysis gene sets in colon cancer was identified by gene set enrichment analysis. Then, a prognostic signature was built up after Cox regression analyses, and overall survival analysis validated the predictive ability. Immune status and infiltration in cancer tissues were explored using the single sample gene set enrichment analysis and CIBERSORT algorithm. A nomogram model integrating clinical variables and the gene signature was established and assessed.

Results

Altogether, 378 cancer and 39 control cases were enrolled. Three glycolysis gene sets and two hypoxia gene sets were enriched in colon cancer (P < 0.05). Five independent genes (ENO3, GPC1, P4HA1, SPAG4, and STC2) were significantly correlated with prognosis of colon cancer patients. Patients with higher risks had significantly better prognosis than those with lower risks (P = 0.002 and AUC = 0.750), which was also observed in the elderly, female and stage I–II subgroups (P < 0.05). In high-risk cases, proportion of NK cells resting increased (P < 0.05) while that of dendritic cells activated (P < 0.05), dendritic cells resting (P < 0.01) and monocytes (P < 0.01) decreased. Besides, expressions of 22 checkpoint genes were found abnormal in groups with different risks (P < 0.05). The predictive nomogram presented satisfactory performance with C-index of 0.771 (0.712–0.830). The area under ROC curve was 0.796 and 0.803 for 3- and 5-year survival prediction, respectively.

Conclusion

A glycolysis and hypoxia combined gene signature was a promising method to evaluate the prognosis and immune infiltration of colon cancer patients, which may provide a new tool for cancer management.

Association between Altered Oncogenic Signaling Pathways and Overall Survival of Patients with Metastatic Colorectal Cancer

Systemic characterization of genomic alterations into signaling pathways helps to understand the molecular pathogenies of colorectal cancer; however, their clinical implications remain unclear. Here, 128 patients with metastatic colorectal cancer (mCRC) receiving targeted next generation sequencing were retrospectively enrolled to analyze the impact of altered oncogenic pathways on clinical outcome. The datasets from Memorial Sloan Kettering Cancer Center were used for validation. In 123 patients with non-MSI-high tumor, the most common mutated gene was TP53 (84.6%), followed by APC (78.0%), KRAS (49.6%), and SMAD4 (22.8%). When mutated genes were allocated into signaling pathways defined as The Cancer Genome Atlas Pan-Cancer Analysis Project, alterations of cell cycle, Wnt, p53, RTK-RAS, PI3K, TGF-β, Notch, and Myc pathways were identified in 88%, 87%, 85%, 75%, 28%, 26%, 17%, and 10% of mCRC tissues, respectively. The survival analyses revealed that Myc and TGF-β pathway alterations were associated with a shorter overall survival (OS) (hazard ratio [HR]: 2.412; 95% confidence interval [CI]: 1.139–5.109; p = 0.018 and HR: 2.754; 95% CI: 1.044–7.265; p = 0.033, respectively). The negative prognostic impact of altered TGF-β pathway was maintained in patients receiving an anti-EGFR antibody. The OS of patients with mCRC carrying MYC and BRAF mutation was shorter than those with either MYC or BRAF mutation (HR: 4.981, 95% CI: 0.296–83.92; p = 0.02). These findings have clinical implications, such as prognosis prediction, treatment guidance, and molecular-targeted therapy development.

TCGA database analysis of the tumor mutation burden and its clinical significance in colon cancer

Background

Colon cancer is one of the most common malignant tumors, with high rates of incidence and death. The tumor mutational burden (TMB), which is characterized by microsatellite instability, has been becoming a powerful predictor which can show tumor behavior and response to immunotherapy.

Methods

In this study, we analyzed 437 mutation data of colon cancer samples obtained from The Cancer Genome Atlas (TCGA) and divided patients into low- and high-TMB groups according to the TMB value. Then we identified differentially-expressed genes (DEGs), conducted immune cell infiltration and survival analyses between groups.

Results

The higher TMB of the patients with colon cancer predicts a poorer prognosis. Functional analysis was performed to assess the prognostic value of the top 30 core genes. The CIBER-SORT algorithm was used to investigate the correlation between the immune cells and TMB subtypes. An immune prognosis model was constructed to screen out immune genes related to prognosis, and the tumor immunity assessment resource (TIMER) was then used to determine the correlation between gene expression and the abundance of tumor-infiltrating immune cell subsets in colon cancer. We observed that APC, TP53, TTN, KRAS, MUC16, SYNE1, PIK3CA have higher somatic mutations. DEGs enrichment analysis showed that they are involved in the regulation of neuroactive ligand-receptor interaction, the Cyclic adenosine monophosphate (cAMP) signaling pathway, the calcium signaling pathway, and pantothenate and Coenzyme A (CoA) biosynthesis. The difference in the abundance of various white blood cell subtypes showed that Cluster of Differentiation 8 (CD8) T cells (P=0.008), activated CD4 memory T cells (P=0.019), M1 macrophages (P=0.002), follicular helper T cells (P=0.034), activated Natural killer (NK cell) cells (P=0.017) increased remarkably, while M0 macrophages significantly reduced (P=0.025). The two immune model genes showed that secretin (SCT) was negatively correlated with survival, while Guanylate cyclase activator 2A (GUCA2A) was positively correlated.

Conclusions

This study conducted a systematically comprehensive analysis of the prediction and clinical significance of TMB in colon cancer in identification, monitoring, and prognosis of colon cancer, and providing reference information for immunotherapy.

Multiomics analysis of tumor mutational burden across cancer types

Whether tumor mutational burden (TMB) is related to improved survival outcomes or the promotion of immunotherapy in various malignant tumors remains controversial, and we lack a comprehensive understanding of TMB across cancers. Based on the data obtained from The Cancer Genome Atlas (TCGA), we conducted a multiomics analysis of TMB across 21 cancer types to identify characteristics related to TMB and determine the mechanism as it relates to prognosis, gene expression, gene mutation and signaling pathways. In our study, TMB was found to have a significant relationship with prognosis for 21 tumors, and the relationship was different in different tumors. TMB may also be related to different outcomes for patients with different tumor subtypes. TMB was confirmed to be correlated with clinical information, such as age and sex. Mutations in GATA3 and MAP3K1 in beast invasive carcinoma (BRCA), TCF7L2 in colon adenocarcinoma (COAD), NFE2L2 in esophageal carcinoma (ESCA), CIC and IDH1 in brain lower grade glioma (LGG), CDH1 in stomach adenocarcinoma (STAD), and TP53 in uterine corpus endometrial carcinoma (UCEC) were demonstrated to be correlated with lower TMB. Moreover, we identified differentially expressed genes (DEGs) and differentially methylated regions (DMRs) according to different TMB levels in 21 cancers. We also investigated the correlation between enrichment of signaling pathways, immune cell infiltration and TMB. In conclusion, we identified multiomic characteristics related to the TMB in 21 tumors, providing support for a comprehensive understanding of the role of TMB in different tumors.

Tumor Mutational Burden Predicting the Efficacy of Immune Checkpoint Inhibitors in Colorectal Cancer: A Systematic Review and Meta-Analysis

Objectives

For colorectal cancer patients, traditional biomarker deficient mismatch repair/microsatellite instability (dMMR/MSI) is an accurate predictor of immune checkpoint inhibitors (ICIs). Recent years, researchers considered tumor mutation burden (TMB) as another predictive biomarker which means the number of nonsynonymous mutations in cancer cells. Several studies have proven that TMB can evaluate the efficacy of ICI therapy in diverse types of cancer, especially in non-small cell lung cancer and melanoma. However, studies on the association between TMB and the response to ICI therapy in colorectal cancer alone are still lacking. In this study, we aim to verify the effect of TMB as a biomarker in predicting the efficacy of ICIs in colorectal cancer.

Methods

We searched the PubMed and Ovid MEDLINE databases up to May 1, 2021 and screened studies for eligibility. Thirteen studies published from 2015 to 2021 with 5062 patients were included finally. We extracted and calculated hazard ratios (HRs) and odds ratios (ORs) of overall survival (OS) and objective response rates (ORRs) and their 95% confidence intervals (95% CIs). Pooled HR and OR were evaluated to compare OS and ORR between TMB-high and TMB-low groups in colorectal cancer patients. Meanwhile, we assessed heterogeneity with the I² statistic and p-values and performed publication bias assessments, sensitivity analyses, and subgroup analyses to search the cause of heterogeneity.

Results

The TMB-high patient group had a longer OS than the TMB-low patient group (HR = 0.68, 95% CI: 0.51, 0.92, p = 0.013) among colorectal cancer patients receiving ICIs. In addition, the TMB-high patient group was superior in terms of ORR (OR = 19.25, 95% CI: 10.06, 36.82, p < 0.001) compared to the TMB-low patient group.

Conclusions

In conclusion, this meta-analysis revealed that TMB can be used as a potential predictive biomarker of colorectal cancer patients receiving ICI therapy. Nevertheless, this finding is not stable enough. Therefore, many more randomized controlled trials are needed to prove that TMB is reliable enough to be used clinically to predict the efficacy of immunotherapy in colorectal cancer. And the most relevant biomarker remains to be determined when TMB high overlaps with other biomarkers like MSI and TILs.

Establishment of a novel ferroptosis-related lncRNA pair prognostic model in colon adenocarcinoma

Long non-coding RNAs (lncRNAs) have been reported to be prognostic factors for cancer. Ferroptosis is an iron-dependent process of programmed cell death. Here, we established a ferroptosis-related lncRNA (frlncRNA) pair signature and revealed its prognostic value in colon adenocarcinoma (COAD) by analyzing the data from The Cancer Genome Atlas (TCGA). FrlncRNAs were identified based on co-expression analysis using the Pearson correlation. Differentially expressed frlncRNAs (DEfrlncRNAs) were recognized and paired, followed by prognostic assessment using univariate Cox regression analysis. The least absolute shrinkage and selection operator (LASSO) penalized Cox analysis was used to determine and construct a risk score prognostic model, by which the receiver operating characteristic (ROC) curves for predicting the overall survival (OS) were conducted. Following the evaluation of whether it was an independent prognostic factor, correlations between the risk score model and clinicopathological characteristics, hypoxia- and immune-related factors, and somatic variants were investigated. In total, 148 DEfrlncRNA pairs were identified, 25 of which were involved in a risk score prognostic signature. The area under ROC curves (AUCs) representing the predictive effect for 1-, 3-, and 5-year survival rates were 0.860, 0.885, and 0.934, respectively. The risk score model was confirmed as an independent prognostic factor and was significantly superior to the clinicopathological characteristics. Correlation analyses showed disparities in clinicopathological characteristics, hypoxia- and immune-related factors, and somatic variants, as well as specific signaling pathways between high- and low-risk groups. The novel risk score prognostic model constructed by pairing DEfrlncRNAs showed promising clinical prediction value in COAD.

The Adaptive Immune Landscape of the Colorectal Adenoma-Carcinoma Sequence

BACKGROUND

The tumor immune microenvironment exerts a pivotal influence in tumor initiation and progression. The aim of this study was to analyze the immune context of sporadic and familial adenomatous polyposis (FAP) lesions along the colorectal adenoma-carcinoma sequence (ACS).

METHODS

We analyzed immune cell counts (CD3+, CD4+, CD8+, Foxp3+, and CD57+), tumor mutation burden (TMB), MHC-I expression and PD-L1 expression of 59 FAP and 74 sporadic colorectal lesions, encompassing adenomas with low-grade dysplasia (LGD) (30 FAP; 30 sporadic), adenomas with high-grade dysplasia (22 FAP; 30 sporadic), and invasive adenocarcinomas (7 FAP; 14 sporadic).

RESULTS

The sporadic colorectal ACS was characterized by (1) a stepwise decrease in immune cell counts, (2) an increase in TMB and MHC-I expression, and (3) a lower PD-L1 expression. In FAP lesions, we observed the same patterns, except for an increase in TMB along the ACS. FAP LGD lesions harbored lower Foxp3+ T cell counts than sporadic LGD lesions. A decrease in PD-L1 expression occurred earlier in FAP lesions compared to sporadic ones.

CONCLUSIONS

The colorectal ACS is characterized by a progressive loss of adaptive immune infiltrate and by the establishment of a progressively immune cold microenvironment. These changes do not appear to be related with the loss of immunogenicity of tumor cells, or to the onset of an immunosuppressive tumor microenvironment.