Cross-platform Data Analysis Reveals a Generic Gene Expression Signature for Microsatellite Instability in Colorectal Cancer

MSIsensor-RNA: Microsatellite Instability Detection for Bulk and Single-cell Gene Expression Data

Microsatellite instability (MSI) is an indispensable biomarker in cancer immunotherapy. Currently, MSI scoring methods by high-throughput omics methods have gained popularity and demonstrated better performance than the gold standard method for MSI detection. However, the MSI detection method on expression data, especially single-cell expression data, is still lacking, limiting the scope of clinical application and prohibiting the investigation of MSI at a single-cell level. Herein, we developed MSIsensor-RNA, an accurate, robust, adaptable, and standalone software to detect MSI status based on expression values of MSI-associated genes. We demonstrated the favorable performance and promise of MSIsensor-RNA in both bulk and single-cell gene expression data in multiplatform technologies including RNA sequencing (RNA-seq), microarray, and single-cell RNA-seq. MSIsensor-RNA is a versatile, efficient, and robust method for MSI status detection from both bulk and single-cell gene expression data in clinical studies and applications. MSIsensor-RNA is available at https://github.com/xjtu-omics/msisensor-rna.

Detecting microsatellite instability by length comparison of microsatellites in the 3' untranslated region with RNA-seq

Microsatellite instability (MSI), a phenomenon caused by deoxyribonucleic acid (DNA) mismatch repair system deficiencies, is an important biomarker in cancer research and clinical diagnostics. MSI detection often involves next-generation sequencing data, with many studies focusing on DNA. Here, we introduce a novel approach by measuring microsatellite lengths directly from ribonucleic acid sequencing (RNA-seq) data and comparing its distribution to detect MSI. Our findings reveal distinct instability patterns between MSI-high (MSI-H) and microsatellite stable samples, indicating the efficacy of RNA-based MSI detection. Additionally, microsatellites in the 3'-untranslated regions showed the greatest predictive value for MSI detection. Notably, this efficacy extends to detecting MSI-H samples even in tumors not commonly associated with MSI. Our approach highlights the utility of RNA-seq data in MSI detection, facilitating more precise diagnostics through the integration of various biological data.

MSIMEP: Predicting microsatellite instability from microarray DNA methylation tumor profiles

Deficiency in DNA MMR activity results in tumors with a hypermutator phenotype, termed microsatellite instability (MSI). Beyond its utility in Lynch syndrome screening algorithms, today MSI has gained importance as predictive biomarker for various anti-PD-1 therapies across many different tumor types. Over the past years, many computational methods have emerged to infer MSI using either DNA- or RNA-based approaches. Considering this together with the fact that MSI-high tumors frequently exhibit a hypermethylated phenotype, herein we developed and validated MSIMEP, a computational tool for predicting MSI status from microarray DNA methylation tumor profiles of colorectal cancer samples. We demonstrated that MSIMEP optimized and reduced models have high performance in predicting MSI in different colorectal cancer cohorts. Moreover, we tested its consistency in other tumor types with high prevalence of MSI such as gastric and endometrial cancers. Finally, we demonstrated better performance of both MSIMEP models vis-à-vis a MLH1 promoter methylation-based one in colorectal cancer.

Using empirical biological knowledge to infer regulatory networks from multi-omics data

Background

Integration of multi-omics data can provide a more complex view of the biological system consisting of different interconnected molecular components, the crucial aspect for developing novel personalised therapeutic strategies for complex diseases. Various tools have been developed to integrate multi-omics data. However, an efficient multi-omics framework for regulatory network inference at the genome level that incorporates prior knowledge is still to emerge.

Results

We present IntOMICS, an efficient integrative framework based on Bayesian networks. IntOMICS systematically analyses gene expression, DNA methylation, copy number variation and biological prior knowledge to infer regulatory networks. IntOMICS complements the missing biological prior knowledge by so-called empirical biological knowledge, estimated from the available experimental data. Regulatory networks derived from IntOMICS provide deeper insights into the complex flow of genetic information on top of the increasing accuracy trend compared to a published algorithm designed exclusively for gene expression data. The ability to capture relevant crosstalks between multi-omics modalities is verified using known associations in microsatellite stable/instable colon cancer samples. Additionally, IntOMICS performance is compared with two algorithms for multi-omics regulatory network inference that can also incorporate prior knowledge in the inference framework. IntOMICS is also applied to detect potential predictive biomarkers in microsatellite stable stage III colon cancer samples.

Conclusions

We provide IntOMICS, a framework for multi-omics data integration using a novel approach to biological knowledge discovery. IntOMICS is a powerful resource for exploratory systems biology and can provide valuable insights into the complex mechanisms of biological processes that have a vital role in personalised medicine.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04891-9.

Tumor microenvironment-aware, single-transcriptome prediction of microsatellite instability in colorectal cancer using meta-analysis

Detecting microsatellite instability (MSI) in colorectal cancers (CRCs) is essential because it is the determinant of treatment strategies, including immunotherapy and chemotherapy. Yet, no attempt has been made to exploit transcriptomic profile and tumor microenvironment (TME) of it to unveil MSI status in CRC. Hence, we developed a novel TME-aware, single-transcriptome predictor of MSI for CRC, called MAP (Microsatellite instability Absolute single sample Predictor). MAP was developed utilizing recursive feature elimination-random forest with 466 CRC samples from The Cancer Genome Atlas, and its performance was validated in independent cohorts, including 1118 samples. MAP showed robustness and predictive power in predicting MSI status in CRC. Additional advantages for MAP were demonstrated through comparative analysis with existing MSI classifier and other cancer types. Our novel approach will provide access to untouched vast amounts of publicly available transcriptomic data and widen the door for MSI CRC research and be useful for gaining insights to help with translational medicine.

Experimental and Meta-Analytic Validation of RNA Sequencing Signatures for Predicting Status of Microsatellite Instability

Microsatellite instability (MSI) is an important diagnostic and prognostic cancer biomarker. In colorectal, cervical, ovarian, and gastric cancers, it can guide the prescription of chemotherapy and immunotherapy. In laboratory diagnostics of susceptible tumors, MSI is routinely detected by the size of marker polymerase chain reaction products encompassing frequent microsatellite expansion regions. Alternatively, MSI status is screened indirectly by immunohistochemical interrogation of microsatellite binding proteins. RNA sequencing (RNAseq) profiling is an emerging source of data for a wide spectrum of cancer biomarkers. Recently, three RNAseq-based gene signatures were deduced for establishing MSI status in tumor samples. They had 25, 15, and 14 gene products with only one common gene. However, they were developed and tested on the incomplete literature of The Cancer Genome Atlas (TCGA) sampling and never validated experimentally on independent RNAseq samples. In this study, we, for the first time, systematically validated these three RNAseq MSI signatures on the literature colorectal cancer (CRC) (n = 619), endometrial carcinoma (n = 533), gastric cancer (n = 380), uterine carcinosarcoma (n = 55), and esophageal cancer (n = 83) samples and on the set of experimental CRC RNAseq samples (n = 23) for tumors with known MSI status. We found that all three signatures performed well with area under the curve (AUC) ranges of 0.94-1 for the experimental CRCs and 0.94-1 for the TCGA CRC, esophageal cancer, and uterine carcinosarcoma samples. However, for the TCGA endometrial carcinoma and gastric cancer samples, only two signatures were effective with AUC 0.91-0.97, whereas the third signature showed a significantly lower AUC of 0.69-0.88. Software for calculating these MSI signatures using RNAseq data is included.

Multi-omics characterization and validation of MSI-related molecular features across multiple malignancies

HEADINGS AIMS

To establish a microsatellite instability (MSI) predictive model in pan-cancer and compare the multi-omics characterization of MSI-related molecular features.

MATERIALS AND METHODS

We established a 15-gene signature for predicting MSI status and performed a systematic assessment of MSI-related molecular features including gene and miRNA expression, DNA methylation, and somatic mutation, in approximately 10,000 patients across 30 cancer types from The Cancer Genome Atlas, Gene Expression Omnibus database, and our institution. Then we identified common MSI-associated dysregulated molecular features across six cancers and explored their mutual interfering relationships and the drug sensitivity.

KEY FINDINGS

we demonstrated the model's high prediction performance and found the samples with high-MSI were mainly distributed in six cancers: BRCA, COAD, LUAD, LIHC, STAD, and UCEC. We found RPL22L1 was up-regulated in the high-MSI group of 5/6 cancer types. CYP27A1 and RAI2 were down-regulated in 4/6 cancer types. More than 20 miRNAs and 39 DMGs were found up-regulated in MSI-H at least three cancers. We discovered some drugs, including OSI-027 and AZD8055 had a higher sensitivity in the high MSI-score group. Functional enrichment analysis revealed the correlation between MSI score and APM score, HLA score, or glycolysis score. The complicated regulatory mechanism of tumor MSI status in multiple dimensions was explored by an integrated analysis of the correlations among MSI-related genes, miRNAs, methylation, and drug response data.

SIGNIFICANCE

Our pan-cancer study provides a valuable predictive model and a comprehensive atlas of tumor MSI, which may guide more precise and personalized therapeutic strategies for tumor patients.

Novel Epigenetic Eight-Gene Signature Predictive of Poor Prognosis and MSI-Like Phenotype in Human Metastatic Colorectal Carcinomas

Simple Summary

The global methylation profile of two human metastatic colorectal carcinoma subgroups with significantly different outcomes (primary-resistant versus drug-sensitive tumors) was analyzed and compared with the gene expression and methylation data from The Cancer Genome Atlas COlon ADenocarcinoma (TCGA COAD) metastatic colorectal carcinoma dataset with the aim to identify a prognostic signature of functionally methylated genes. A novel epigenetic eight-gene signature, with hypermethylation of the promoter regions, was identified and validated for its capacity to predict poor outcome, which had a CpG-island methylator phenotype (CIMP)-high status and microsatellite instability (MSI)-like phenotype.

Abstract

Epigenetics is involved in tumor progression and drug resistance in human colorectal carcinoma (CRC). This study addressed the hypothesis that the DNA methylation profiling may predict the clinical behavior of metastatic CRCs (mCRCs). The global methylation profile of two human mCRC subgroups with significantly different outcome was analyzed and compared with gene expression and methylation data from The Cancer Genome Atlas COlon ADenocarcinoma (TCGA COAD) and the NCBI GENE expression Omnibus repository (GEO) GSE48684 mCRCs datasets to identify a prognostic signature of functionally methylated genes. A novel epigenetic signature of eight hypermethylated genes was characterized that was able to identify mCRCs with poor prognosis, which had a CpG-island methylator phenotype (CIMP)-high and microsatellite instability (MSI)-like phenotype. Interestingly, methylation events were enriched in genes located on the q-arm of chromosomes 13 and 20, two chromosomal regions with gain/loss alterations associated with adenoma-to-carcinoma progression. Finally, the expression of the eight-genes signature and MSI-enriching genes was confirmed in oxaliplatin- and irinotecan-resistant CRC cell lines. These data reveal that the hypermethylation of specific genes may provide prognostic information that is able to identify a subgroup of mCRCs with poor prognosis.

PreMSIm: An R package for predicting microsatellite instability from the expression profiling of a gene panel in cancer

Microsatellite instability (MSI) is a genomic property of the cancers with defective DNA mismatch repair and is a useful marker for cancer diagnosis and treatment in diverse cancer types. In particular, MSI has been associated with the active immune checkpoint blockade therapy response in cancer. Most of computational methods for predicting MSI are based on DNA sequencing data and a few are based on mRNA expression data. Using the RNA-Seq pan-cancer datasets for three cancer cohorts (colon, gastric, and endometrial cancers) from The Cancer Genome Atlas (TCGA) program, we developed an algorithm (PreMSIm) for predicting MSI from the expression profiling of a 15-gene panel in cancer. We demonstrated that PreMSIm had high prediction performance in predicting MSI in most cases using both RNA-Seq and microarray gene expression datasets. Moreover, PreMSIm displayed superior or comparable performance versus other DNA or mRNA-based methods. We conclude that PreMSIm has the potential to provide an alternative approach for identifying MSI in cancer.