Cancer-specific functional profiling in microsatellite-unstable (MSI) colon and endometrial cancers using combined differentially expressed genes and biclustering analysis

Abstract

Microsatellite-unstable (MSI) cancers have distinct genetic and clinical features from microsatellite-stable cancers, but the molecular functional differences between MSI cancers originating from different tissues or organs have not been well studied because the application of usual differentially expressed gene (DEG) analysis is error-prone, producing too many noncancer-specific normally functioning genes. To maximize therapeutic efficacy, biomarkers reflecting cancer-specific differences between MSI cancers of different tissue origins should be identified. To identify functional differences between MSI colon and endometrial cancers, we combined DEG analysis and biclustering instead of DEG analysis alone and refined functionally relevant biclusters reflecting genuine functional differences between the 2 tumors. Specifically, using The Cancer Genome Atlas and genome-tissue expression as data sources, gene ontology (GO) enrichment tests were performed after routinely identifying DEGs between the 2 tumors with the exclusion of DEGs identified in their normal counterparts. Cancer-specific biclusters and associated enriched GO terms were obtained by biclustering with enrichment tests for the preferences for cancer type (either colon or endometrium) and GO enrichment tests for each cancer-specific bicluster, respectively. A novel childness score was developed to select functionally relevant biclusters among cancer-specific biclusters based on the extent to which the enriched GO terms of the biclusters tended to be child terms of the enriched GO terms in DEGs. The selected biclusters were tested using survival analysis to validate their clinical significance. We performed multiple sequential analyses to produce functionally relevant biclusters from the RNA sequencing data of MSI colon and endometrial cancer samples and their normal counterparts. We identified 3066 cancer-specific DEGs. Biclustering analysis revealed 153 biclusters and 41 cancer-specific biclusters were selected using Fisher exact test. A mean childness score over 0.6 was applied as the threshold and yielded 8 functionally relevant biclusters from cancer-specific biclusters. Functional differences appear to include gland cavitation and the TGF-β receptor, G protein, and cytokine pathways. In the survival analysis, 6 of the 8 functionally relevant biclusters were statistically significant. By attenuating noise and applying a synergistic contribution of DEG results, we refined candidate biomarkers to complement tissue-specific features of MSI tumors.