A Gene Signature Derived from the Loss of CDKN1A (p21) Is Associated with CMS4 Colorectal Cancer

SIMPLE SUMMARY: A gene signature derived from the loss of CDKN1A (p21) gene, obtained in HCT116 p21-/- colorectal cancer cells, is identified in a large cohort of primary colorectal (CRC) tumors and is associated with the Consensus Molecular Subtype (CMS) of colon cancer that has a worse relapse-free and overall survival, that is, CMS4 (also called mesenchymal subtype). The presented gene signature can help to uncover the early molecular mechanisms of epithelial–mesenchymal transition (EMT), which is known to be associated with high stemness and drug resistance. ABSTRACT: The epithelial–mesenchymal transition (EMT) is associated with tumor aggressiveness and increased invasion, migration, metastasis, angiogenesis, and drug resistance. Although the HCT116 p21-/- cell line is well known for its EMT-associated phenotype, with high Vimentin and low E-cadherin protein levels, the gene signature of this rather intermediate EMT-like cell line has not been determined so far. In this work, we present a robust molecular and bioinformatics analysis, to reveal the associated gene expression profile and its correlation with different types of colorectal cancer tumors. We compared the quantitative signature obtained with the NanoString platform with the expression profiles of colorectal cancer (CRC) Consensus Molecular Subtypes (CMS) as identified, and validated the results in a large independent cohort of human tumor samples. The expression signature derived from the p21-/- cells showed consistent and reliable numbers of upregulated and downregulated genes, as evaluated with two machine learning methods against the four CRC subtypes (i.e., CMS1, 2, 3, and 4). High concordance was found between the upregulated gene signature of HCT116 p21-/- cells and the signature of the CMS4 mesenchymal subtype. At the same time, the upregulated gene signature of the native HCT116 cells was similar to that of CMS1. Using a multivariate Cox regression model to analyze the survival data in the CRC tumor cohort, we selected genes that have a predictive risk power (with a significant gene risk incidence score). A set of genes of the mesenchymal signature was proven to be significantly associated with poor survival, specifically in the CMS4 CRC human cohort. We suggest that the gene signature of HCT116 p21-/- cells could be a suitable metric for mechanistic studies regarding the CMS4 signature and its functional consequences in CRC. Moreover, this model could help to discover the molecular mechanisms of intermediate EMT, which is known to be associated with extraordinarily high stemness and drug resistance.