Identification and validation of the clinical prediction model and biomarkers based on chromatin regulators in colon cancer by integrated analysis of bulk- and single-cell RNA sequencing data.

BACKGROUNDChromatin regulators (CRs) are implicated in the development of cancer, but a comprehensive investigation of their role in colon adenocarcinoma (COAD) is inadequate. The purpose of this study is to find CRs that can provide recommendations for clinical diagnosis and treatment, and to explore the reasons why they serve as critical CRs.

METHODSWe obtained data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Weighted Gene Co-Expression Network Analysis (WGCNA) screened tumor-associated CRs. LASSO-Cox regression was used to construct the model and to screen key CRs together with support vector machine (SVM), the univariate Cox regression. We used single-cell data to explore the expression of CRs in cells and their communication. Immune infiltration, immune checkpoints, mutation, methylation, and drug sensitivity analyses were performed. Gene expression was verified by quantitative real-time reverse transcription-polymerase chain reaction (qRT-PCR). Pan-cancer analysis was used to explore the importance of hub CRs.

RESULTSWe finally obtained 32 tumor-associated CRs. The prognostic model was constructed based on RCOR2, PPARGC1A, PKM, RAC3, PHF19, MYBBP1A, ORC1, and EYA2 by the LASSO-Cox regression. Single-cell data revealed that the model was immune-related. Combined with immune infiltration analysis, immune checkpoint analysis, and tumor immune dysfunction and exclusion (TIDE) analysis, the low-score risk group had more immune cell infiltration and better immune response. Mutation and methylation analysis showed that multiple CRs may be mutated and methylated in colon cancer. Drug sensitivity analysis revealed that the low-risk group may be more sensitive to several drugs and PKM was associated with multiple drugs. Combined with machine learning, PKM is perhaps the most critical gene in CRs. Pan-cancer analysis showed that PKM plays a role in the prognosis of cancers.

CONCLUSIONSWe developed a prognostic model for COAD based on CRs. Increased expression of the core gene PKM is linked with a poor prognosis in several malignancies.