Telomere-related gene risk model for prognosis prediction in colorectal cancer

Colorectal cancer (CRC) is the third-most prevalent cancer globally. The biological significance of telomeres in CRC carcinogenesis and progression is underscored by accumulating data. Nevertheless, not much is known about how telomere-related genes (TRGs) affect CRC prognosis. Therefore, the aim of this study was to investigate the role of TRGs in CRC prognosis.

We retrospectively obtained the expression profiles and clinical data of CRC patients from public databases. Utilizing least absolute shrinkage and selection operator (LASSO) regression analysis, we created a telomere-related risk model to predict survival outcomes, identifying ten telomere-related differentially expressed genes (TRDEGs). Based on TRDEGs, we stratified patients from The Cancer Genome Atlas (TCGA) into low- and high-risk subsets. Subsequently, we conducted comprehensive analyses, including survival assessment, immune cell infiltration, drug sensitivity, and prediction of molecular interactions using Kaplan-Meier curves, ESTIMATE, CIBERSORT, OncoPredict, and other approaches.

The model showed exceptional predictive accuracy for survival. Significant differences in survival were observed between the two groups of participants grouped according to the model (P<0.001), and this difference was further confirmed in the external validation set (GSE39582) (P=0.004). Additionally, compared to the low-risk group, the high-risk group exhibited significantly advanced tumor node metastasis (TNM) stages, lower proportions of activated CD4⁺ T cells, effector memory CD4⁺ T cells, and memory B cells, but increased ratios of M2 macrophages and regulatory T cells (Tregs), elevated tumor immune dysfunction and exclusion (TIDE) scores, and diminished sensitivity to dabrafenib, lapatinib, camptothecin, docetaxel, and telomerase inhibitor IX, reflecting the signature's capacity to distinguish clinical pathological characteristics, immune environment, and drug efficacy. Finally, we validated the expression of the ten TRDEGs (ACACB, TPX2, SRPX, PPARGC1A, CD36, MMP3, NAT2, MMP10, HIGD1A, and MMP1) through quantitative real-time polymerase chain reaction (qRT-PCR) and found that compared to normal cells, the expression levels of ACACB, HIGD1A, NAT2, PPARGC1A, and TPX2 in CRC cells were elevated, whereas those of CD36, SRPX, MMP1, MMP3, and MMP10 were reduced.

Overall, we constructed a telomere-related biomarker capable of predicting prognosis and treatment response in CRC individuals, offering potential guidance for drug therapy selection and prognosis prediction.