Dataset Information

Analysis of core genes for colorectal cancer prognosis based on immune and stromal scores

ABSTRACT:

Background

Colorectal cancer (CRC) is one of the most common malignancies.An early diagnosis and an accurate prognosis are major focuses of CRC research. Tumor microenvironment cells and the extent of infiltrating immune and stromal cells contribute significantly to the tumor prognosis.

Methods

Immune and stromal scores were calculated based on the ESTIMATE algorithm using the sample expression profile of the The Cancer Genome Atlas (TCGA) database. GSE102479 was used as the validation database. Differentially expressed genes whose expression was significantly associated with the prognosis of CRC patients were identified based on the immune matrix score. Survival analysis was conducted on the union of the differentially expressed genes. A protein–protein interaction (PPI) network was constructed using the STRING database to identify the closely connected modules. To conduct functional enrichment analysis of the relevant genes, GO and KEGG pathway analyses were performed with Cluster Profiler. Pivot analysis of the ncRNAs and TFs was performed by using the RAID2.0 database and TRRUST v2 database. TF-mRNA regulatory relationships were analyzed in the TRRUST V2 database. Hubgene targeting relationships were screened in the TargetScan, miRTarBase and miRDB databases. The SNV data of the hub genes were analyzed by using the R maftools package. A ROC curve was drawn based on the TCGA database. The proportion of immune cells was estimated using CIBERSORT and the LM22 feature matrix.

Results

The results showed that the matrix score was significantly correlated with colorectal cancer stage T. A total of 789 differentially expressed genes and 121 survival-related prognostic genes were identified. The PPI network showed that 22 core genes were related to the CRC prognosis. Furthermore, four ncRNAs that regulated the core prognosis genes, 11 TFs with regulatory effects on the core prognosis genes, and two drugs, quercetin and pseudoephedrine, that have regulatory effects on colorectal cancer were also identified.

Conclusions

We obtained a list of tumor microenvironment-related genes for CRC patients. These genes could be useful for determining the prognosis of CRC patients. To confirm the function of these genes, additional experiments are necessary.

SUBMITTER: Zhu Y

PROVIDER: S-EPMC8607933 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Background: Colorectal cancer (CRC) is one gastrointestinal malignancy, accounting for 10% of cancer diagnoses and cancer-related deaths worldwide each year. Therefore, it is urgent to identify genes involved in CRC predicting the prognosis. Methods: CRC's data were acquired from the Gene Expression Omnibus (GEO) database (GSE39582 and GSE41258 datasets) and The Cancer Genome Atlas (TCGA) database. The differentially expressed necroptosis-related genes (DENRGs) were sorted out between tumor and normal tissues. Univariate Cox regression analysis and least absolute shrinkage and selectionator operator (LASSO) analysis were applied to selected DENRGs concerning patients' overall survival and to construct a prognostic biomarker. The effectiveness of this biomarker was assessed by the Kaplan-Meier curve and the receiver operating characteristic (ROC) analysis. The GSE39582 dataset was utilized as external validation for the prognostic signature. Moreover, using univariate and multivariate Cox regression analyses, independent prognostic factors were identified to construct a prognostic nomogram. Next, signaling pathways regulated by the signature were explored through the gene set enrichment analysis (GSEA). The single sample gene set enrichment analysis (ssGSEA) algorithm and tumor immune dysfunction and exclusion (TIDE) were used to explore immune correlation in the two groups, high-risk and low-risk ones. Finally, prognostic genes' expression was examined in the GSE41258 dataset. Results: In total, 27 DENRGs were filtered, and a necroptosis-related prognostic signature based on 6 DENRGs was constructed, which may better understand the overall survival (OS) of CRC. The Kaplan-Meier curve manifested the effectiveness of the prognostic signature, and the ROC curve showed the same result. In addition, univariate and multivariate Cox regression analyses revealed that age, pathology T, and risk score were independent prognostic factors, and a nomogram was established. Furthermore, the prognostic signature was most significantly associated with the apoptosis pathway. Meanwhile, 24 immune cells represented significant differences between two groups, like the activated B cell. Furthermore, 32 immune checkpoints, TIDE scores, PD-L1 scores, and T-cell exclusion scores were significantly different between the two groups. Finally, a 6-gene prognostic signature represented different expression levels between tumor and normal samples significantly in the GSE41258 dataset. Conclusion: Our study established a signature including 6 genes and a prognostic nomogram that could significantly assess the prognosis of patients with CRC.

Project description:The tumor environment is of vital importance for the incidence and development of colorectal cancer. Increasing evidence in recent years has elaborated the vital role of the tumor environment in cancer subtype classification and patient prognosis, but a comprehensive understanding of the colorectal tumor environment that is purely dependent on the stromal compartment is lacking. To decipher the tumor environment in colorectal cancer and explore the role of its immune context in cancer classification, we performed a gene expression microarray on the stromal compartment of colorectal cancer and adjacent normal tissues. Through the integrated analysis of our data with public gene expression microarray data of stromal and epithelial colorectal cancer tissues processed through laser capture microdissection, we identified four highly connected gene modules representing the biological features of four tissue compartments by applying a weighted gene coexpression network analysis algorithm and classified colorectal cancers into three immune subtypes by adopting a nearest template prediction algorithm. A systematic analysis of the four identified modules mainly reflected the close interplay between the biological changes of intrinsic and extrinsic characteristics at the initiation of colorectal cancer. Colorectal cancers were stratified into three immune subtypes based on gene templates identified from representative gene modules of the stromal compartment: active immune, active stroma, and mixed type. These immune subtypes differed by the immune cell infiltration pattern, expression of immune checkpoint inhibitors, mutation landscape, extent of mutation burden, extent of copy number burden, prognosis and chemotherapeutic sensitivity. Further analysis indicated that activation of the NF-kB signaling pathway was the major mechanism causing the no immune infiltration milieu in the active stroma subtype and that inhibitors of the NF-kB signaling pathway could be candidate drugs for treating patients with an active stroma. Overall, these results suggest that characterizing colorectal cancer by the tumor environment is of vital importance in predicting patients' clinical outcomes and helping guide precision and personalized treatment.