Identification of colon cancer subtypes based on multi-omics data-construction of methylation markers for immunotherapy

Background

Being the most widely used biomarker for immunotherapy, the microsatellite status has limitations in identifying all patients who benefit in clinical practice. It is essential to identify additional biomarkers to guide immunotherapy. Aberrant DNA methylation is consistently associated with changes in the anti-tumor immune response, which can promote tumor progression. This study aims to explore immunotherapy biomarkers for colon cancers from the perspective of DNA methylation.

Methods

The related data (RNA sequencing data and DNA methylation data) were obtained from The Cancer Genome Atlas (TCGA) and UCSC XENA database. Methylation-driven genes (MDGs) were identified through the Pearson correlation analysis. Unsupervised consensus clustering was conducted using these MDGs to identify distinct clusters of colon cancers. Subsequently, we evaluated the immune status and predicted the efficacy of immunotherapy by tumor immune dysfunction and exclusion (Tide) score. Finally, The Quantitative Differentially Methylated Regions (QDMR) software was used to identify the specific DNA methylation markers within particular clusters.

Results

A total of 282 MDGs were identified by integrating the DNA methylation and RNA-seq data. Consensus clustering using the K-means algorithm revealed that the optimal number of clusters was 4. It was revealed that the composition of the tumor immune microenvironment (TIME) in Cluster 1 was significantly different from others, and it exhibited a higher level of tumor mutation burdens (TMB) and stronger anti-tumor immune activity. Furthermore, we identified three specific hypermethylation genes that defined Cluster 1 (PCDH20, APCDD1, COCH). Receiver operating characteristic (ROC) curves demonstrated that these specific markers could effectively distinguish Cluster 1 from other clusters, with an AUC of 0.947 (95% CI 0.903-0.990). Finally, we selected clinical samples for immunohistochemical validation.

Conclusion

In conclusion, through the analysis of DNA methylation, consensus clustering of colon cancer could effectively identify the cluster that benefit from immunotherapy along with specific methylation biomarkers.

Prognostic Value of Tumor Size in Colon Cancer-Smaller is Better?

BACKGROUND

The prognostic value of tumor size in colon cancer remains controversial. This study aimed to reveal the correlation between tumor size and prognosis of colon cancer.

METHODS

A total of 491 patients with colon cancer were included in this study. The correlation of tumor size with prognosis, mismatch repair status, and other clinicopathological characteristics as well as tumor microenvironment was analyzed.

RESULTS

For stage IIA microsatellite stable (MSS) colon cancer, tumors sized <3.5 cm and ≥5 cm were associated with a poorer disease free survival (DFS) compared with tumors sized between 3.5 and 5 cm (P = .002). Small tumor size (HR = 5.098, P = .001) and large tumor size (HR = 2.749, P = .029) were found to be independent prognostic factors for stage IIA MSS colon cancer. Moreover, high expression of transgelin (TAGLN), a marker of cancer-associated fibroblasts (CAFs), was found to be an independent prognostic factor for poorer DFS (HR = 9.651, P = .009), which was also associated with smaller tumor size (P = .027).

CONCLUSION

Small (<3.5 cm) and large (≥5 cm) tumor sizes are associated with decreased DFS in stage IIA MSS colon cancer. Enrichment of TAGLN+ CAFs is associated with decreased DFS and small tumor size.

NAT10 mediated mRNA acetylation modification patterns associated with colon cancer progression and microsatellite status

N4-acetylcytidine (ac4C) is one type of RNA modification found in eukaryotes. RNA acetylation modifications are gradually expanding in oncology. However, the role of RNA acetylation modifications in colorectal cancer and its association with colorectal cancer microsatellite status remain unclear. Using public databases and in vitro experiments, we verified the expression and biological function of NAT10, as the key RNA acetylation modification enzyme, in colorectal cancer. The results showed that NAT10 was highly expressed in colorectal cancer, and significantly promoted colorectal cancer cell proliferation. NAT10 was also involved in several aspects of cell homoeostasis such as ion transport, calcium-dependent phospholipid binding, and RNA stability. NAT10 expression positively correlated with immune infiltration in colorectal cancer. We further constructed a risk regression model for mRNA acetylation in colorectal cancer using acetylation-related differential genes. We found that tumour immune infiltration, microsatellite instability (MSI) proportion, tumour immune mutation burden, and patient response to immunotherapy were positively correlated with risk scores. For the first time, our study showed that the level of mRNA acetylation modification level is elevated in colorectal cancer and positively correlates with immune infiltration and microsatellite status of patients. Based on our findings, NAT10 may be a new target for colorectal cancer treatment.

Identification of a 3-gene signature for predicting the prognosis of stage II colon cancer based on microsatellite status

Background

Though colon cancer (CC) is one of the most malignant tumors across the world, CC patients with microsatellite instability-high (MSI-H) in stage II seem to have a better prognosis. However, the molecular mechanisms underlying the phenomena haven't been elucidated yet.

Methods

This study enrolled 322 CCs with known microsatellite status from GSE143985, GSE39582 and GSE92921 in the Gene Expression Omnibus (GEO) database. Robust rank aggregation (RRA) analysis, univariate Cox regression analysis and multivariate Cox stepwise regression analysis were performed to identify genes and construct risk score signature. Kaplan-Meier and receiver operating characteristic (ROC) curves analyses were used to evaluate the prognostic value of the signature. The potential mechanisms underlying this signature were assessed in the Metascape database, gene set enrichment analysis (GSEA) and immune infiltration analysis.

Results

RRA analysis identified 40 differently expressed genes (DEGs). A 3-gene risk score signature (MKQ signature) associated with disease-free survival (DFS) was generated. DFS was significantly longer in CC patients with lower than higher scores (P=0.0046). The areas under curves (AUCs) of the time-dependent ROC curves of MKQ signature at 1-, 3- and 5-year DFS were 1, 0.963 and 0.961 respectively. Recurrence-free survival (RFS) was significantly longer in patients in GSE39582 with lower than higher risk scores (P=0.032). The AUCs for 1-, 3- and 5-year RFS in GSE39582 were 0.63, 0.618 and 0.583, respectively, validating the value of the MKQ signature. Functional annotation and GSEA revealed that the MKQ signature was associated with multiple immune-related pathways. Immune cell infiltration was found to differ in patients differing in the MKQ signature.

Conclusions

Gene expression and microsatellite status identified a 3-gene signature (MKQ signature) that could facilitate risk-stratified management in patients with stage II CC. Dysregulation of MSMB, KRT23, and QPRT can serve as prognostic markers in stage II CC.

Construction of a novel methylation-related prognostic model for colorectal cancer based on microsatellite status

The present study aimed to construct a novel methylation-related prognostic model based on microsatellite status that may enhance the prognosis of colorectal cancer (CRC) from methylation and microsatellite status perspective. DNA methylation and mRNA expression data with clinical information were downloaded from The Cancer Genome Atlas (TCGA) data set. The samples were divided into microsatellite stability and microsatellite instability group, and CIBERSORT was used to assess the immune cell infiltration characteristics. After identifying the differentially methylated genes and differentially expression genes using R packages, the methylation-driven genes were further identified. Prognostic genes that were used to establish the methylation-related risk score model were generated by the univariate and multivariate Cox regression model. Finally, we established and evaluated the methylation-related prognostic model for CRC patients. A total of 69 MDGs were obtained and three of these genes (MIOX, TH, DKFZP434K028) were selected to construct the prognostic model. Patients in the low-risk score group had a conspicuously better overall survival than those in the high-risk score group (p < .0001). The area under the receiver operating characteristic curve for this model was 0.689 at 3 years, 0.674 at 4 years, and 0.658 at 5 years. The Wilcoxon test showed that higher risk score was associated with higher T stage (p = .01), N stages (p = .0028), metastasis (p = .013), and advanced pathological stage (p = .0013). However, the more instability of microsatellite status, the lower risk score of CRC patients (p = .0048). Our constructed methylation-related prognostic model based on microsatellite status presents potential significance in assessing recurrence risk stratification, tumor staging, and immunotherapy for CRC patients.

Quantitative Prediction of Microsatellite Instability in Colorectal Cancer With Preoperative PET/CT-Based Radiomics

OBJECTIVES

Microsatellite instability (MSI) status is an important hallmark for prognosis prediction and treatment recommendation of colorectal cancer (CRC). To address issues due to the invasiveness of clinical preoperative evaluation of microsatellite status, we investigated the value of preoperative 18F-FDG PET/CT radiomics with machine learning for predicting the microsatellite status of colorectal cancer patients.

METHODS

A total of 173 patients that underwent 18F-FDG PET/CT scans before operations were retrospectively analyzed in this study. The microsatellite status for each patient was identified as microsatellite instability-high (MSI-H) or microsatellite stable (MSS), according to the test for mismatch repair gene proteins with immunohistochemical staining methods. There were 2,492 radiomic features in total extracted from 18F-FDG PET/CT imaging. Then, radiomic features were selected through multivariate random forest selection and univariate relevancy tests after handling the imbalanced dataset through the random under-sampling method. Based on the selected features, we constructed a BalancedBagging model based on Adaboost classifiers to identify the MSI status in patients with CRC. The model performance was evaluated by the area under the curve (AUC), sensitivity, specificity, and accuracy on the validation dataset.

RESULTS

The ensemble model was constructed based on two radiomic features and achieved an 82.8% AUC for predicting the MSI status of colorectal cancer patients. The sensitivity, specificity, and accuracy were 83.3, 76.3, and 76.8%, respectively. The significant correlation of the selected two radiomic features with multiple effective clinical features was identified (p < 0.05).

CONCLUSION

18F-FDG PET/CT radiomics analysis with the machine learning model provided a quantitative, efficient, and non-invasive mechanism for identifying the microsatellite status of colorectal cancer patients, which optimized the treatment decision support.

Relationship between microsatellite status and immune microenvironment of colorectal cancer and its application to diagnosis and treatment

Due to advances in understanding the immune microenvironment of colorectal cancer (CRC), microsatellite classification (dMMR/MSI-H and pMMR/MSS) has become a key biomarker for the diagnosis and treatment of CRC patients and therefore has important clinical value. Microsatellite status is associated with a variety of clinicopathological features and affects drug resistance and the prognosis of patients. CRC patients with different microsatellite statuses have different compositions and distributions of immune cells and cytokines within their tumor microenvironments (TMEs). Therefore, there is great interest in reversing or reshaping CRC TMEs to transform immune tolerant "cold" tumors into immune sensitive "hot" tumors. This requires a thorough understanding of differences in the immune microenvironments of MSI-H and MSS type tumors. This review focuses on the relationship between CRC microsatellite status and the immune microenvironment. It focuses on how this relationship has value for clinical application in diagnosis and treatment, as well as exploring the limitations of its current application.

Synchronous Colorectal Cancer: Improving Accuracy of Detection and Analyzing Molecular Heterogeneity-The Main Keys for Optimal Approach

BACKGROUND

In patients with synchronous colorectal cancer (SCRC), understanding the underlying molecular behavior of such cases is mandatory for designing individualized therapy. The aim of this paper is to highlight the importance of transdisciplinary evaluation of the pre- and post-operative assessment of patients with SCRCs, from imaging to molecular investigations.

METHODS

Six patients with SCRCs presented with two carcinomas each. In addition to the microsatellite status (MSS), the epithelial mesenchymal transition was checked in each tumor using the biomarkers β-catenin and E-cadherin, same as KRAS and BRAF mutations.

RESULTS

In two of the patients, the second tumor was missed at endoscopy, but diagnosed by a subsequent computed-tomography-scan (CT-scan). From the six patients, a total of 11 adenocarcinomas (ADKs) and one squamous cell carcinoma (SCC) were analyzed. All the examined carcinomas were BRAF-wildtype microsatellite stable tumors with an epithelial histological subtype. In two of the six cases, KRAS gene status showed discordance between the two synchronous tumors, with mutations in the index tumors and wildtype status in the companion ones.

CONCLUSIONS

Preoperative CT-scans can be useful for detection of synchronous tumors which may be missed by colonoscopy. Where synchronous tumors are identified, therapy should be based on the molecular profile of the indexed tumors.

Prognostic relevance of programmed death-ligand 1 expression and microsatellite status in small bowel adenocarcinoma

Background: Small bowel adenocarcinoma (SBA) is a dreadful disease. Patient prognosis is limited due to late presentation and ineffective chemotherapy. PD-1/PD-L1 checkpoint immunotherapy is regarded as a promising approach in several cancer entities. The association of PD-1/PD-L1 expression and its impact on patient prognosis with SBA is unclear. Material and methods: Seventy-five consecutive patients who underwent surgery for SBA were retrospectively analyzed and stained for PD-L1 expression in the tumour or the stroma. Analysis of mismatch repair genes was performed to determine microsatellite status. Kaplan-Meier estimate was used to analyze patient survival. Univariate and multivariable Cox regression-analyses were used to assess the impact of PD-L1 expression and microsatellite status on patient survival.Results: PD-L1 was weakly upregulated within the tumour or the stroma and associated with prolonged survival (p = .0071 and p = .0472, respectively). Fifty-one tumours (68%) revealed microsatellite stability (MSS) and 24 tumours (32%) were microsatellite instable (MSI) without correlating with patient survival (p = .611). Neither PD-L1 expression in the tumour nor in the stroma was identified as an independent risk factor influencing survival (p = .572 and p = .3055).Conclusion: Although PD-L1 expression is associated with prolonged survival, it was not identified as an independent prognostic marker. Microsatellite status did not influence long-term survival.

Impact of microsatellite status on chemotherapy for colorectal cancer patients with KRAS or BRAF mutation

KRAS and BRAF mutations are frequently detected in cases of colorectal cancer (CRC). The microsatellite status of patients with CRC and mutated KRAS/BRAF is important when determining cancer therapy. In the present study, the microsatellite status and genetic polymorphisms of KRAS (codons 12 and 13) and BRAF (V600E) were characterized in CRC tissue. The mismatch repair activity and oncogenic potential of KRAS were assessed by immunoblots from two KRAS-mutated CRC cell lines, SW480 and HCT116, with different microsatellite statuses, following treatment with 5-fluorouracil (5-FU) and oxaliplatin. Of all the 205 patients with CRC enrolled in the present study, 31.2% (64 of 205) had a KRAS or BRAF mutation, and 79.7% (51 of 64) of these patients with a KRAS/BRAF mutation exhibited microsatellite stability (MSS), indicating that microsatellite status is correlated with KRAS/BRAF mutation (P=0.027). A higher proportion (39.0%, 41 of 105) of elderly patients (≥62.6 years) had mutated KRAS or BRAF than younger patients (<62.6 years; 23.0%, 23 of 100; P=0.013). In the subgroup of 154 patients with MSS, patients without the KRAS or BRAF mutation (n=110) had longer disease-specific survival rates (58.8±9.4%) than patients with KRAS or BRAF mutations (n=44; 50.6±11.0%; P=0.043). Cytoplasmic KRAS levels decreased whereas nuclear MutS protein homolog 2 (MSH2) levels increased slightly in CRC HCT116 cells that were microsatellite instable, following treatment with 76.9 µM 5-FU for 2 days. In microsatellite stable SW480 cells, MSH2 levels markedly increased in the nucleus following 150 µM oxaliplatin treatment for 3 days. However, no significant change was observed regarding KRAS distribution in these cells. The results of the present study suggest that it is important to identify patients with CRC who may benefit from adjuvant chemotherapy with 5-FU or oxaliplatin, particularly CRC patients with MSS and mutated KRAS or BRAF, who have poorer overall survival rates than patients with microsatellite instability. Knowledge of the microsatellite status of patients and whether they harbor KRAS or BRAF mutations may enable more effective therapeutic strategies to be developed. Further prospective studies are required to validate the findings of the current study.

Gastric cancer in young vs old Romanian patients: immunoprofile with emphasis on maspin and mena protein reactivity

Increasing number of early-onset gastric carcinomas (GCs) and controversial results regarding the differences among young and older patients with this type of cancer are the reasons why correlation of clinicopathological factors with molecular markers is necessary. The aim of our study was to compare the demographic, clinical and immunohistochemical (IHC) aspects in Romanian patients with GC diagnosed below and above 45 years old. In 191 samples provided from patients with GC, the clinicopathological parameters were correlated with a panel of 15 antibodies: E-cadherin, HER-2, VEGF, CD31, CD105, COX-2, maspin, bax, bcl-2, p53, Ki67, MLH-1, MSH-2, mena protein and vimentin. Compared to the conventional cases, GCs diagnosed below 45 years old were more frequently located at the gastroesophageal junction and presented a higher percentage of lymph node metastases. The diffuse type E-cadherin/mena/p53/Ki67/bax-negative cases that displayed nuclear maspin positivity were also more frequently in younger patients. The intestinal type early-onset GCs were the most angiogenic ones, the apoptotic rate being lower than in the intestinal type GCs of the aged. Compared to the conventional cases, in the early-onset GCs the nuclear maspin-mediated antiproliferative activity is more intense in diffuse type while the mena-dependent tumor cell proliferation is more characteristic for intestinal type GCs.