RESUMO
Background: Colon cancer (CC) ranks as one of the leading causes of cancer-related mortality globally. Single-cell transcriptome sequencing (scRNA-seq) offers precise gene expression data for distinct cell types. This study aimed to utilize scRNA-seq and bulk transcriptome sequencing (bulk RNA-seq) data from CC samples to develop a novel prognostic model. Methods: scRNA-seq data was downloaded from the GSE161277 database. R packages including "Seurat", "Harmony", and "singleR" were employed to categorize eight major cell types within normal and tumor tissues. By comparing tumor and normal samples, differentially expressed genes (DEGs) across these major cell types were identified. Gene Ontology (GO) enrichment analyses of DEGs for each cell type were conducted using "Metascape". DEGs-based signature construction involved Cox regression and least absolute shrinkage operator (LASSO) analyses, performed on The Cancer Genome Atlas (TCGA) training cohort. Validation occurred in the GSE39582 and GSE33382 datasets. The expression pattern of prognostic genes was verified using spatial transcriptome sequencing (ST-seq) data. Ultimately, an established prognostic nomogram based on the gene signature and age was established and calibrated. Sensitivity to chemotherapeutic drugs was predicted with the "oncoPredict" R package. Results: Using scRNA-Seq data, we examined 33,213 cells, categorizing them into eight cell types within normal and tumor samples. GO enrichment analysis revealed various cancer-related pathways across DEGs in these cell types. Among the 55 DEGs identified via univariate Cox regression, four independent prognostic genes emerged: PTPN6, CXCL13, SPINK4, and NPDC1. Expression validation through ST-seq confirmed PTPN6 and CXCL13 predominance in immune cells, while SPINK4 and NPDC1 were relatively epithelial cell-specific. Creating a four-gene prognostic signature, Kaplan-Meier survival analyses emphasized higher risk scores correlating with unfavorable prognoses, confirmed across training and validation cohorts. The risk score emerged as an independent prognostic factor, supported by a reliable nomogram. Intriguingly, drug sensitivity analysis unveiled contrasting anti-cancer drug responses in the two risk groups, suggesting significant clinical implications. Conclusion: We developed a novel prognostic four-gene risk model, and these genes may act as potential therapeutic targets for CC.
RESUMO
Background: The extracellular matrix (ECM) plays a vital role in progression, expansion, and prognosis of malignancies. In this study, we aimed to explore a novel ECM-based prognostic model for patients with colon cancer (CC). Methods: ECM-related genes were obtained from Molecular Signatures database. Differential expression analysis was performed using the CC dataset from The Cancer Genome Atlas (TCGA) database. Four ECM-related genes related to overall survival were identified using the Cox regression and LASSO analysis. Then an ECM-related signature was developed and verified in three independent CC cohorts (GSE33882, GSE39582 and GSE29621) from the Gene Expression Omnibus (GEO). A prognostic nomogram was developed incorporating the ECM-related gene signature with clinical risk factors. CIBERSORT was used to explore the immune cell infiltration level. Human Protein Atlas (HPA) database was utilized to validate the expression levels of identified prognostic ECM genes. Results: Four ECM-related genes (CXCL13, CXCL14, SFRP5 and THBS4) were identified to develop an ECM-based gene signature and demarcated CC patients into the high- and low-risk groups. In training and validation datasets, patients in the low-risk group had better overall survival outcomes than those in the high-risk group (log-rank P<0.001). In addition, ECM-related signature was significantly associated with consensus molecular subtype 4 (CMS4) as well as other known clinical risk factors such as a higher Tumor, Nodal Involvement, Metastasis (TNM) stage. Moreover, the risk score derived from the ECM-based gene signature could be utilized as an independent prognostic factor for CC patients. A nomogram including the ECM-related gene signature, age and stage was developed to serve clinical practice. CIBERSORT analysis showed immune cell infiltration was different between high- and low-risk groups. The immunohistochemical results derived from HPA indicated differential expression of prognosis-related ECM genes in CC and normal tissues. Conclusions: In the present study, a novel risk model based on ECM-signature could effectively reflect individual risk classification and provide potential therapeutic targets for CC patients. Moreover, the prognostic nomogram may help predict individualized survival.
RESUMO
Background: Gastric cancer (GC) is a highly heterogeneous tumor with different responses to immunotherapy. Identifying immune subtypes and landscape of GC could improve immunotherapeutic strategies. Methods: Based on the abundance of tumor-infiltrating immune cells in GC patients from The Cancer Genome Atlas, we used unsupervised consensus clustering algorithm to identify robust clusters of patients, and assessed their reproducibility in an independent cohort from Gene Expression Omnibus. We further confirmed the feasibility of our immune subtypes in five independent pan-cancer cohorts. Finally, functional enrichment analyses were provided, and a deep learning model studying the pathological images was constructed to identify the immune subtypes. Results: We identified and validated three reproducible immune subtypes presented with diverse components of tumor-infiltrating immune cells, molecular features, and clinical characteristics. An immune-inflamed subtype 3, with better prognosis and the highest immune score, had the highest abundance of CD8+ T cells, CD4+ T-activated cells, follicular helper T cells, M1 macrophages, and NK cells among three subtypes. By contrast, an immune-excluded subtype 1, with the worst prognosis and the highest stromal score, demonstrated the highest infiltration of CD4+ T resting cells, regulatory T cells, B cells, and dendritic cells, while an immune-desert subtype 2, with an intermediate prognosis and the lowest immune score, demonstrated the highest infiltration of M2 macrophages and mast cells, and the lowest infiltration of M1 macrophages. Besides, higher proportion of EVB and MSI of TCGA molecular subtyping, over expression of CTLA4, PD1, PDL1, and TP53, and low expression of JAK1 were observed in immune subtype 3, which consisted with the results from Gene Set Enrichment Analysis. These subtypes may suggest different immunotherapy strategies. Finally, deep learning can predict the immune subtypes well. Conclusion: This study offers a conceptual frame to better understand the tumor immune microenvironment of GC. Future work is required to estimate its reference value for the design of immune-related studies and immunotherapy selection.
