Identification of a Genome Instability-Associated LncRNA Signature for Prognosis Prediction in Colon Cancer

Identification of Genomic Instability-Associated LncRNAs as Potential Therapeutic Targets in Lung Adenocarcinoma

BACKGROUND/OBJECTIVES

Non-small cell lung cancer (NSCLC) is the most common type of cancer, with lung adenocarcinoma (LUAD) as the predominant subtype. Despite advancements in targeted therapies, many NSCLC patients still experience poor outcomes due to treatment resistance and disease progression. Genomic instability (GI), a hallmark of cancer, defined as the increased tendency of DNA mutations and alterations, is closely linked to cancer initiation, progression, and resistance to therapy. Emerging evidence suggests that long non-coding RNAs (lncRNAs)-molecules longer than 200 nucleotides that do not encode proteins but regulate gene expression-play critical roles in cancer biology and are associated with GI. However, the relationship between GI and lncRNA expression in LUAD remains poorly understood.

METHODS

In this study, we analyzed the transcript profiles of lncRNAs and mRNAs from LUAD samples in The Cancer Genome Atlas (TCGA) database and classified them based on their Homologous Recombination Deficiency (HRD) score. The HRD score is an unweighted sum of three independent DNA-based measures of genomic instability: loss of heterozygosity, telomeric allelic imbalance, and large-scale transitions. We then performed a differential gene expression analysis to identify lncRNAs and mRNAs that were either upregulated or downregulated in samples with high HRD scores compared to those with low HRD scores. Following this, we conducted a correlation analysis to assess the significance of the association between HRD scores and the expression of both lncRNAs and mRNAs.

RESULTS

We identified 30 differentially expressed lncRNAs and 200 mRNAs associated with genomic instability. Using an RNA interactome database from sequencing experiments, we found evidence of interactions between GI-associated lncRNAs (GI-lncRNAs) and GI-associated mRNAs (GI-mRNAs). Further investigation showed that some GI-lncRNAs play regulatory and functional roles in LUAD and other diseases. We also found that GI-lncRNAs have potential as prognostic biomarkers, particularly when integrated with HRD stratification. The expression of specific GI-lncRNAs was associated with primary therapy response and immune infiltration in LUAD. Additionally, we identified existing drugs that could modulate GI-lncRNAs, offering potential therapeutic strategies to address GI in LUAD.

CONCLUSIONS

Our findings suggest that GI-associated lncRNAs could serve as valuable biomarkers for LUAD prognosis and therapeutic response. Furthermore, modulating these lncRNAs presents potential treatment avenues to address genomic instability in LUAD.

Integrative Analysis of Whole-Genome and Transcriptomic Data Reveals Novel Variants in Differentially Expressed Long Noncoding RNAs Associated with Asthenozoospermia

Background/Objectives: Asthenozoospermia, characterized by reduced sperm motility, is a common cause of male infertility. Emerging evidence suggests that noncoding RNAs, particularly long noncoding RNAs (lncRNAs), play a critical role in the regulation of spermatogenesis and sperm function. Coding regions have a well-characterized role and established predictive value in asthenozoospermia. However, this study was designed to complement previous findings and provide a more holistic understanding of asthenozoospermia, this time focusing on noncoding regions. This study aimed to identify and prioritize variants in differentially expressed (DE) lncRNAs found exclusively in asthenozoospermic men, focusing on their impact on lncRNA structure and lncRNA-miRNA-mRNA interactions. Methods: Whole-genome sequencing (WGS) was performed on samples from asthenozoospermic and normozoospermic men. Additionally, an RNA-seq dataset from normozoospermic and asthenozoospermic individuals was analyzed to identify DE lncRNAs. Bioinformatics analyses were conducted to map unique variants on DE lncRNAs, followed by prioritization based on predicted functional impact. The structural impact of the variants and their effects on lncRNA-miRNA interactions were assessed using computational tools. Gene ontology (GO) and KEGG pathway analyses were employed to investigate the affected biological processes and pathways. Results: We identified 4173 unique variants mapped to 258 DE lncRNAs. After prioritization, 5 unique variants in 5 lncRNAs were found to affect lncRNA structure, while 20 variants in 17 lncRNAs were predicted to disrupt miRNA-lncRNA interactions. Enriched pathways included Wnt signaling, phosphatase binding, and cell proliferation, all previously implicated in reproductive health. Conclusions: This study identifies specific variants in DE lncRNAs that may play a role in asthenozoospermia. Given the limited research utilizing WGS to explore the role of noncoding RNAs in male infertility, our findings provide valuable insights and a foundation for future studies.

Disulfideptosis-associated lncRNAs reveal features of prognostic, immune escape, tumor mutation, and tumor malignant progression in renal clear cell carcinoma

PURPOSE

Investigating the role of lncRNAs associated with the latest cell death mode (Disulfideptosis) in renal clear cell carcinoma, as well as their correlation with tumor prognosis, immune escape, immune checkpoints, tumor mutational burden, and malignant tumor progression. Searching for potential biomarkers and targets for renal clear cell carcinoma.

METHODS

Downloaded the expression profile data and clinical data of 533 cases of renal clear cell carcinoma from the TCGA database, and randomly divided them into a test set (267 cases) and a validation set (266 cases). Based on previous research, 13 genes associated with Disulfideptosis were obtained. Using R software, lncRNAs with a differential expression that is related to the prognosis of renal clear cell carcinoma and associated with Disulfideptosis were screened out. After univariate Cox regression analysis, Lasso regression analysis, and multivariate Cox regression analysis, lncRNAs with independent predictive ability were obtained. A predictive risk model was established based on the risk scores. Verification was carried out between the obtained high-risk and low-risk groups and their subgroups (including Age, Gender, tumor mutational burden (TMB), tumor grading, and staging). Subsequently, a nomogram was established, and a calibration curve was generated for verification. Performed GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) functional enrichment analyses. Downloaded the values of Tumor Immune Dysfunction and Exclusion (TIDE) for all samples and calculated the difference between the high and low-risk groups. Selected human renal tumor cell lines (786-O, OS-RC-2, A-498, ACHN) and human renal cortex proximal tubule epithelial cell line (HK-2). The RNA expression levels of the above lncRNAs in each cell line were analyzed using RT-qPCR (Real-time Quantitative PCR Detecting System). Used siRNA (small interfering RNA) to knock down FAM225B in 786-O and OS-RC-2 cell lines, and then performed in vitro cell experiments to validate the functional characteristics of FAM225B.

RESULTS

Our constructed predictive model includes 5 lncRNAs with an independent predictive ability (FAM225B, ZNF503-AS1, SPINT1-AS1, WWC2-AS2, LINC01338), which can effectively distinguish between patients in high and low-risk groups and their subgroups. The 1, 3, and 5-year AUC (Area Under the ROC Curve) values of the established nomogram are 0.756, 0.752, and 0.781, respectively. The 5-year AUC value is higher compared to other clinical characteristics (Age: 0.598, Gender: 0.488, Grade: 0.680, Stage: 0.717). After the knockdown of FAM225B, the proliferation, migration, and invasion abilities of renal cancer cell lines OS-RC-2 and 786-O all decreased.

CONCLUSION

We have constructed and validated a prognostic model based on Disulfideptosis-associated lncRNAs. This model can effectively predict the high or low risk of patient prognosis and can distinguish the tumor cell mutational burden and immune escape capabilities among high-risk and low-risk patients. This predictive model can serve as an independent prognostic factor for renal clear cell carcinoma, providing a new direction for personalized treatment of patients with renal clear cell carcinoma.

Multi-Omics Mining of lncRNAs with Biological and Clinical Relevance in Cancer

In this review, we provide a general overview of the current panorama of mining strategies for multi-omics data to investigate lncRNAs with an actual or potential role as biological markers in cancer. Several multi-omics studies focusing on lncRNAs have been performed in the past with varying scopes. Nevertheless, many questions remain regarding the pragmatic application of different molecular technologies and bioinformatics algorithms for mining multi-omics data. Here, we attempt to address some of the less discussed aspects of the practical applications using different study designs for incorporating bioinformatics and statistical analyses of multi-omics data. Finally, we discuss the potential improvements and new paradigms aimed at unraveling the role and utility of lncRNAs in cancer and their potential use as molecular markers for cancer diagnosis and outcome prediction.

A tumor mutational burden-derived immune computational framework selects sensitive immunotherapy/chemotherapy for lung adenocarcinoma populations with different prognoses

Background

Lung adenocarcinoma (LUAD) kills millions of people every year. Recently, FDA and researchers proved the significance of high tumor mutational burden (TMB) in treating solid tumors. But no scholar has constructed a TMB-derived computing framework to select sensitive immunotherapy/chemotherapy for the LUAD population with different prognoses.

Methods

The datasets were collected from TCGA, GTEx, and GEO. We constructed the TMB-derived immune lncRNA prognostic index (TILPI) computing framework based on TMB-related genes identified by weighted gene co-expression network analysis (WGCNA), oncogenes, and immune-related genes. Furthermore, we mapped the immune landscape based on eight algorithms. We explored the immunotherapy sensitivity of different prognostic populations based on immunotherapy response, tumor immune dysfunction and exclusion (TIDE), and tumor inflammation signature (TIS) model. Furthermore, the molecular docking models were constructed for sensitive drugs identified by the pRRophetic package, oncopredict package, and connectivity map (CMap).

Results

The TILPI computing framework was based on the expression of TMB-derived immune lncRNA signature (TILncSig), which consisted of AC091057.1, AC112721.1, AC114763.1, AC129492.1, LINC00592, and TARID. TILPI divided all LUAD patients into two populations with different prognoses. The random grouping verification, survival analysis, 3D PCA, and ROC curve (AUC=0.74) firmly proved the reliability of TILPI. TILPI was associated with clinical characteristics, including smoking and pathological stage. Furthermore, we estimated three types of immune cells threatening the survival of patients based on multiple algorithms. They were macrophage M0, T cell CD4 Th2, and T cell CD4 memory activated. Nevertheless, five immune cells, including B cell, endothelial cell, eosinophil, mast cell, and T cell CD4 memory resting, prolonged the survival. In addition, the immunotherapy response and TIDE model proved the sensitivity of the low-TILPI population to immunotherapy. We also identified seven intersected drugs for the LUAD population with poor prognosis, which included docetaxel, gemcitabine, paclitaxel, palbociclib, pyrimethamine, thapsigargin, and vinorelbine. Their molecular docking models and best binding energy were also constructed and calculated.

Conclusions

We divided all LUAD patients into two populations with different prognoses. The good prognosis population was sensitive to immunotherapy, while the people with poor prognosis benefitted from 7 drugs.

Epigenetic reprogramming in cancer: From diagnosis to treatment

Disruption of the epigenetic program of gene expression is a hallmark of cancer that initiates and propagates tumorigenesis. Altered DNA methylation, histone modifications and ncRNAs expression are a feature of cancer cells. The dynamic epigenetic changes during oncogenic transformation are related to tumor heterogeneity, unlimited self-renewal and multi-lineage differentiation. This stem cell-like state or the aberrant reprogramming of cancer stem cells is the major challenge in treatment and drug resistance. Given the reversible nature of epigenetic modifications, the ability to restore the cancer epigenome through the inhibition of the epigenetic modifiers is a promising therapy for cancer treatment, either as a monotherapy or in combination with other anticancer therapies, including immunotherapies. Herein, we highlighted the main epigenetic alterations, their potential as a biomarker for early diagnosis and the epigenetic therapies approved for cancer treatment.

Long Non-coding RNAs and CRISPR-Cas Edition in Tumorigenesis

Long non-coding RNAs (lncRNAs) are one of the most abundant and heterogeneous transcripts with key roles in chromatin remodeling and gene regulation at the transcriptional and post-transcriptional levels. Due to their role in cell growth and differentiation, lncRNAs have emerged as an important biomarker in cancer diagnosis, prognosis, and targeted treatment. Recent studies have focused on elucidating lncRNA function during malignant transformation, tumor progression and drug resistance. The advent of the CRISPR system has made it possible to precisely edit complex genomic loci such as lncRNAs. Thus, we summarized the advances in CRISPR-Cas approaches for functional studies of lncRNAs including gene knockout, knockdown, overexpression and RNA targeting in tumorigenesis and drug resistance. Additionally, we highlighted the perspectives and potential applications of CRISPR approaches to treat cancer, as an emerging and promising target therapy.

Long non‑coding RNAs interact with RNA‑binding proteins to regulate genomic instability in cancer cells (Review)

Genomic instability, a feature of most cancers, contributes to malignant cell transformation and cancer progression due to the accumulation of genetic alterations. Genomic instability is reflected at numerous levels, from single nucleotide to the chromosome levels. However, the exact molecular mechanisms and regulators of genomic instability in cancer remain unclear. Growing evidence indicates that the binding of long non-coding RNAs (lncRNAs) to protein chaperones confers a variety of regulatory functions, including managing of genomic instability. The aim of the present review was to examine the roles of mitosis, telomeres, DNA repair, and epigenetics in genomic instability, and the mechanisms by which lncRNAs regulate them by binding proteins in cancer cells. This review contributes to our understanding of the role of lncRNAs and genomic instability in cancer and can potentially provide entry points and molecular targets for cancer therapies.

Comprehensive Analyses of Mutation-Derived Long-Chain Noncoding RNA Signatures of Genome Instability in Kidney Renal Papillary Cell Carcinoma

Background: The role of long-chain noncoding RNA (lncRNA) in genomic instability has been demonstrated to be increasingly importance. Therefore, in this study, lncRNAs associated with genomic instability were identified and kidney renal papillary cell carcinoma (KIRP)-associated predictive features were analysed to classify high-risk patients and improve individualised treatment.

Methods: The training (n = 142) and test (n = 144) sets were created using raw RNA-seq and patient’s clinical data of KIRP obtained from The Cancer Genome Atlas (TCGA).There are 27 long-chain noncoding RNAs (lncRNAs) that are connected with genomic instability, these lncRNAs were identified using the ‘limma’ R package based on the numbers of somatic mutations and lncRNA expression profiles acquired from KIRP TCGA cohort. Furthermore, Cox regression analysis was carried out to develop a genome instability-derived lncRNA-based gene signature (GILncSig), whose prognostic value was confirmed in the test cohort as well as across the entire KIRP TCGA dataset.

Results: A GILncSig derived from three lncRNAs (BOLA3-AS1, AC004870, and LINC00839), which were related with poor KIRP survival, was identified, which was split up into high- and low-risk groups. Additionally, the GILncSig was found to be an independent prognostic predictive index in KIRP using univariate and multivariate Cox analysis. Furthermore, the prognostic significance and characteristics of GilncSig were confirmed in the training test and TCGA sets. GilncSig also showed better predictive performance than other prognostic lncRNA features.

Conclusion: The function of lncRNAs in genomic instability and the genetic diversity of KIRP were elucidated in this work. Moreover, three lncRNAs were screened for prediction of the outcome of KIRP survival and novel insights into identifying cancer biomarkers related to genomic instability were discussed.

Identification of LncRNA Prognostic Signature Associated With Genomic Instability in Pancreatic Adenocarcinoma

Background

Genomic instability (GI) is a critical feature of cancer which plays a key role in the occurrence and development of pancreatic adenocarcinoma (PAAD). Long non-coding RNA (LncRNA) is an emerging prognostic biomarker because it is involved in regulating GI. Recently, researchers used such GI-related LncRNAs (GILncRNAs) to establish a prognostic signature for patients with cancer and helped in predicting the overall prognosis of the patients. However, it is evident that patients with PAAD still lack such prognostic signature constructed with GILncRNA.

Methods

The present study screened GILncRNAs from 83 patients with PAAD. Prognosis-related GILncRNAs were identified by univariate Cox regression analysis. The correlation coefficients of these GILncRNAs were obtained by multivariate Cox regression analysis and used to construct a signature. The signature in the present study was then assessed through survival analysis, mutation correlation analysis, independent prognostic analysis, and clinical stratification analysis in the training set and validated in the testing as well as all TCGA set. The current study performed external clinical relevance validation of the signature and validated the effect of AC108134.2 in GILncSig on PAAD using in vitro experiments. Finally, the function of GILncRNA signature (GILncSig) dependent on Gene Ontology enrichment analysis was explored and chemotherapeutic drug sensitivity analysis was also performed.

Results

Results of the present study found that a total of 409 GILncRNAs were identified, 5 of which constituted the prognostic risk signature in this study, namely, AC095057.3, AC108134.2, AC124798.1, AL606834.1, and AC104695.4. It was found that the signature of the present study was better than others in predicting the overall survival and applied to patients with PAAD of all ages, genders, and tumor grades. Further, it was noted that the signature of the current study in the GSE102238, was correlated with tumor length, and tumor stage of patients with PAAD. In vitro, functional experiments were used in the present study to validate that AC108134.2 is associated with PAAD genomic instability and progression. Notably, results of the pRRophetic analysis in the current study showed that the high-risk group possessed reverse characteristics and was sensitive to chemotherapy.

Conclusions

In conclusion, it was evident that the GILncSig used in the present study has good prognostic performance. Therefore, the signature may become a potential sensitive biological indicator of PAAD chemotherapy, which may help in clinical decision-making and management of patients with cancer.

Identification of a 5-Nutrient Stress-Sensitive Gene Signature to Predict Survival for Colorectal Cancer

Background

The high heterogeneity and the complexity of the tumor microenvironment of colorectal cancer (CRC) have enhanced the difficulty of prognosis prediction based on conventional clinical indicators. Recent studies revealed that tumor cells could overcome various nutritional deficiencies by gene regulation and metabolic remodeling. However, whether differentially expressed genes (DEGs) in CRC cells under kinds of nutrient deficiency could be used to predict prognosis remained unveiled.

Methods

Three datasets (GSE70976, GSE13548, and GSE116087), in which colon cancer cells were, respectively, cultured in serum-free, glucose-free, or glutamine-free medium, were included to delineate the profiles of gene expression by nutrient stress. DEGs were figured out in three datasets, and gene functional analysis was performed. Survival analyses and Cox proportional hazards model were then used to identify nutrient stress sensitive genes in CRC datasets (GSE39582 and TCGA COAD). Then, a 5-gene signature was constructed and the risk scores were also calculated. Survival analyses, cox analyses, and nomogram were applied to predict the prognosis of patients with colorectal cancer. The effectiveness of the risk model was also tested.

Results

A total of 48 genes were found to be dysregulated in serum, glucose, or glutamine-deprived CRC cells, which were mainly enriched in cell cycle and endoplasmic reticulum stress pathways. After further analyses, 5 genes, MCM5, MCM6, CDCA2, GINS2, and SPC25, were identified to be differentially expressed in CRC and be related to prognosis of in CRC datasets. We used the above nutrient stress-sensitive genes to construct a risk scoring model. CRC samples in the datasets were divided into low-risk and high-risk groups. Data showed that higher risk scores were associated with better outcomes and risk scores decreased significantly with tumor procession. Moreover, the risk score could be used to predict the probability of survival based on nomogram.

Conclusions

The 5-nutrient stress-sensitive gene signature could act as an independent biomarker for survival prediction of CRC patients.

Selection of lncRNAs That Influence the Prognosis of Osteosarcoma Based on Copy Number Variation Data

Osteosarcoma is the most common primary malignancy in the musculoskeletal system. It is reported that copy number variation- (CNV-) derived lncRNAs contribute to the progression of osteosarcoma. However, whether CNV-derived lncRNAs affect the prognosis of osteosarcoma remains unclear. Here, we obtained osteosarcoma-related CNV data and gene expression profiles from The Cancer Genome Atlas (TCGA) database. CNV landscape analysis indicated that copy number amplification of lncRNAs was more frequent than deletion in osteosarcoma samples. Thirty-four CNV-lncRNAs with DNA-CNV frequencies greater than 30% and their corresponding 294 mRNAs were identified. Gene Ontology (GO) and Kyoto Encyclopedia of Gene and Genome (KEGG) pathway enrichment analyses revealed that these mRNAs were mainly enriched in olfaction, olfactory receptor activity, and olfactory transduction processes. Furthermore, we predicted that a total of 23 genes were cis-regulated by 16 CNV-lncRNAs, while 30 transcription factors (TFs) were trans-regulated by 5 CNV-lncRNAs. Through $t$ -tests, univariate Cox regression analysis, and the least absolute shrinkage and selection operator (LASSO), we constructed a CNV-related risk model including 3 lncRNAs (AC129492.1, PSMB1, and AC037459.4). The Kaplan-Meier (K-M) curves indicated that patients with high-risk scores showed poor prognoses. The areas under the receiver operating characteristic (ROC) curves (AUC) for predicting 3-, 5-, and 7-year overall survival (OS) were greater than 0.7, showing a satisfactory predictive efficiency. Gene set enrichment analysis (GSEA) revealed that the prognostic signature was intimately linked to skeletal system development, immune regulation, and inflammatory response. Collectively, our study developed a novel 3-CNV-lncRNA prognostic signature that would provide theoretical guidance for the clinical prognostic management of osteosarcoma.

Genome Instability and Long Noncoding RNA Reveal Biomarkers for Immunotherapy and Prognosis and Novel Competing Endogenous RNA Mechanism in Colon Adenocarcinoma

Background: Long noncoding RNAs (lncRNAs) crucially modulate DNA damage responses/repair in cancer cells. However, the underlying regulatory role of genome integrity and its clinical value in colon adenocarcinoma (COAD) remains unclear. This study links genome instability to lncRNA using computational biology techniques, in attempt to propose novel biomarkers of immunotherapy outcome, and investigated a potential competing endogenous RNA (ceRNA) as a molecular regulatory mechanism.

Methods: TCGA-COAD patients were divided into genome unstable (GU)-like and genome stable (GS)-like clusters via hierarchical clustering to predict immunotherapy outcomes. Multivariate Cox model was established to predict the overall survival rate in COAD patients. Additionally, SVM and LASSO algorithms were applied to obtain hub lncRNAs. A novel genome instability-related ceRNA network was predicted with the Starbase 2.0 database. To better understand how these genes fundamentally interact during tumor progression and development, the mutation analysis and single-gene analysis for each gene was performed.

Results: In contrast to those in the GS-like cluster, GU-like-cluster patients demonstrated a higher tumor mutational burden (TMB)/microsatellite instability (MSI), DNA polymerase epsilon (POLE) mutation rate, and immune checkpoint expression, all indicate a greater predictive power for response rate for immunotherapy. The novel prognostic signature demonstrated an outstanding predictive performance (AUC > 0.70). The genes in the genome insatiability-related ceRNA network (including four axes: AL161772.1-has-miR-671-5p (hsa-miR-181d-5p, has-miR-106a-5p)-NINL, AL161772.1-has-miR-106a-5p-TNFSF11, AC124067.4-hsa-miR-92b-3p (hsa-miR-589-5p)-PHYHIPL, and BOLA3-AS1-has-miR-130b-3p-SALL4) were identified as critical regulators of tumor microenvironment infiltration, cancer stemness, and drug resistance. qPCR was performed to validate the expression patterns of these genes. Furthermore, the MSI-high proportion was greater in patients with mutated type than in those with the wild type according to all four target genes, indicating that these four genes modulate genomic integrity and could serve as novel immunotherapy biomarkers.

Conclusion: We demonstrated that genome instability-related lncRNA is a novel biomarker for immunotherapy outcomes and prognosis. A novel ceRNA network that modulates genomic integrity, including four lncRNA-miRNA-mRNA axes, was proposed.

A Novel Cancer Stemness-Related Signature for Predicting Prognosis in Patients with Colon Adenocarcinoma

Objective

To explore the cancer stemness features and develop a novel cancer stemness-related prognostic signature for colon adenocarcinoma (COAD).

Methods

We downloaded the mRNA expression data and clinical data of COAD from TCGA database and GEO database. Stemness index, mRNAsi, was utilized to investigate cancer stemness features. Weighted gene coexpression network analysis (WGCNA) was used to identify cancer stemness-related genes. Univariate and multivariate Cox regression analyses were applied to construct a prognostic risk cancer stemness-related signature. We then performed internal and external validation. The relationship between cancer stemness and COAD immune microenvironment was investigated.

Results

COAD patients with higher mRNAsi score or EREG-mRNAsi score have significant longer overall survival (OS). We identified 483 differently expressed genes (DEGs) between the high and low mRNAsi score groups. We developed a cancer stemness-related signature using fifteen genes (including RAB31, COL6A3, COL5A2, CCDC80, ADAM12, VGLL3, ECM2, POSTN, DPYSL3, PCDH7, CRISPLD2, COLEC12, NRP2, ISLR, and CCDC8) for prognosis prediction of COAD. Low-risk score was associated with significantly preferable OS in comparison with high-risk score, and the area under the ROC curve (AUC) for OS prediction was 0.705. The prognostic signature was an independent predictor for OS of COAD. Macrophages, mast cells, and T helper cells were the vital infiltration immune cells, and APC costimulation and type II IFN response were the vital immune pathways in COAD.

Conclusions

We developed and validated a novel cancer stemness-related prognostic signature for COAD, which would contribute to understanding of molecular mechanism in COAD.