Developing a prognosis and chemotherapy evaluating model for colon adenocarcinoma based on mitotic catastrophe-related genes.

Liu Y ，Zhao Y ，Zhang S ，Rong S ，He S ，Hua L ，Wang X ，Chen H ... -  《Scientific Reports》

Construction of a novel choline metabolism-related signature to predict prognosis, immune landscape, and chemotherapy response in colon adenocarcinoma.

Colon adenocarcinoma (COAD) is a common digestive system malignancy with high mortality and poor prognosis. Accumulating evidence indicates that choline metabolism is closely related to tumorigenesis and development. However, the efficacy of choline metabolism-related signature in predicting patient prognosis, immune microenvironment and chemotherapy response has not been fully clarified. Choline metabolism-related differentially expressed genes (DEGs) between normal and COAD tissues were screened using datasets from The Cancer Genome Atlas (TCGA), Kyoto Encyclopedia of Genes and Genomes (KEGG), AmiGO2 and Reactome Pathway databases. Two choline metabolism-related genes (CHKB and PEMT) were identified by univariate and multivariate Cox regression analyses. TCGA-COAD was the training cohort, and GSE17536 was the validation cohort. Patients in the high- and low-risk groups were distinguished according to the optimal cutoff value of the risk score. A nomogram was used to assess the prognostic accuracy of the choline metabolism-related signature. Calibration curves, decision curve analysis (DCA), and clinical impact curve (CIC) were used to improve the clinical applicability of the prognostic signature. Gene Ontology (GO) and KEGG pathway enrichment analyses of DEGs in the high- and low-risk groups were performed. KEGG cluster analysis was conducted by the KOBAS-i database. The distribution and expression of CHKB and PEMT in various types of immune cells were analyzed based on single-cell RNA sequencing (scRNA-seq). The CIBERSORT and ESTIMATE algorithms evaluated tumor immune cell infiltration in the high- and low-risk groups. Evaluation of the half maximal inhibitory concentration (IC50) of common chemotherapeutic drugs based on the choline metabolism-related signature was performed. Small molecule compounds were predicted using the Connectivity Map (CMap) database. Molecular docking is used to simulate the binding conformation of small molecule compounds and key targets. By immunohistochemistry (IHC), Western blot, quantitative reverse transcription-polymerase chain reaction (qRT-PCR) experiments, the expression levels of CHKB and PEMT in human, mouse, and cell lines were detected. We constructed and validated a choline metabolism-related signature containing two genes (CHKB and PEMT). The overall survival (OS) of patients in the high-risk group was significantly worse than that of patients in the low-risk group. The nomogram could effectively and accurately predict the OS of COAD patients at 1, 3, and 5 years. The DCA curve and CIC demonstrate the clinical utility of the nomogram. scRNA-seq showed that CHKB was mainly distributed in endothelial cells, while PEMT was mainly distributed in CD4+ T cells and CD8+ T cells. In addition, multiple types of immune cells expressing CHKB and PEMT differed significantly. There were significant differences in the immune microenvironment, immune checkpoint expression and chemotherapy response between the two risk groups. In addition, we screened five potential small molecule drugs that targeted treatment for COAD. Finally, the results of IHC, Western blot, and qRT-PCR consistently showed that the expression of CHKB in human, mouse, and cell lines was elevated in normal samples, while PMET showed the opposite trend. In conclusion, we constructed a choline metabolism-related signature in COAD and revealed its potential application value in predicting the prognosis, immune microenvironment, and chemotherapy response of patients, which may lay an important theoretical basis for future personalized precision therapy.

Liu C ，Liu D ，Wang F ，Liu Y ，Xie J ，Xie J ，Xie Y ... -  《Frontiers in Immunology》

An Intratumor Heterogeneity-Related Signature for Predicting Prognosis, Immune Landscape, and Chemotherapy Response in Colon Adenocarcinoma.

Colon adenocarcinoma (COAD) is a frequent malignancy of the digestive system with a poor prognosis and high mortality rate worldwide. Intratumor heterogeneity (ITH) is associated with tumor progression, poor prognosis, immunosuppression, and therapy resistance. However, the relationship between ITH and prognosis, the immune microenvironment, and the chemotherapy response in COAD patients remains unknown, and this knowledge is urgently needed. We obtained clinical information and gene expression data for COAD patients from The Cancer Genome Atlas (TCGA) database. The DEPTH2 algorithm was utilized to evaluate the ITH score. X-tile software was used to determine the optimal cutoff value of the ITH score. The COAD patients were divided into high- and low-ITH groups based on the cutoff value. We analyzed prognosis, tumor mutation burden (TMB), gene mutations, and immune checkpoint expression between the high- and low-ITH groups. Differentially expressed genes (DEGs) in the high- and low-ITH groups were subjected to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses. We performed univariate Cox regression and least absolute shrinkage and selection operator (LASSO) regression analyses to screen the prognosis-related genes for the construction of an ITH-related prognostic signature. The nomogram was used to predict the overall survival (OS) of COAD patients. The protein-protein interaction (PPI) network was constructed by using the GeneMANIA database. Principal component analysis (PCA) and single-sample gene set enrichment analysis (ssGSEA) were employed to explore the differences in biological pathway activation status between the high- and low-risk groups. The proportion and type of tumor-infiltrating immune cells were evaluated by the CIBERSORT and ESTIMATE algorithms. Additionally, we assessed the chemotherapy response and predicted small-molecule drugs for treatment. Finally, the expression of the prognosis-related genes was validated by using the UALCAN database and Human Protein Atlas (HPA) database. The OS of the high-ITH group was worse than that of the low-ITH group. A positive correlation between ITH and TMB was identified. In subgroups stratified by age, gender, and tumor stage, the OS of the low-ITH group remained better than that of the high-ITH group. There were dramatic differences in the mutated genes, single nucleotide variant classes, variant types, immune checkpoints and cooccurring and mutually exclusive mutations of the DEGs between the high- and low-ITH groups. Based on the DEGs between the high- and low-ITH groups, we constructed a five-gene signature consisting of CEACAM5, ENO2, GABBR1, MC1R, and SLC44A4. The COAD patients were divided into high- and low-risk groups according to the median risk score. The OS of the high-risk group was worse than that of the low-risk group. The nomogram was used to accurately predict the 1-, 3- and 5-year OS of COAD patients and showed good calibration and moderate discrimination ability. The stromal score, immune score, and ESTIMATE score of the high-risk group were significantly higher than those of the low-risk group, whereas tumor purity showed the opposite trend. The patients classified by the risk score had distinguishable sensitivity to chemotherapeutic drugs. Finally, two public databases confirmed that CEACAM5 and SLC44A4 were upregulated in normal tissues compared with COAD tissues, and ENO2, GABBR1, and MC1R were upregulated in COAD tissues compared with normal tissues. Overall, we identified an ITH-related prognostic signature for COAD that was closely related to the tumor microenvironment and chemotherapy response. This signature may help clinicians make more personalized and precise treatment decisions for COAD patients.

Liu C ，Liu D ，Wang F ，Xie J ，Liu Y ，Wang H ，Rong J ，Xie J ，Wang J ，Zeng R ，Zhou F ，Xie Y ... -  《Frontiers in Medicine》

Identification of a glycolysis- and lactate-related gene signature for predicting prognosis, immune microenvironment, and drug candidates in colon adenocarcinoma.

Liu C ，Liu D ，Wang F ，Xie J ，Liu Y ，Wang H ，Rong J ，Xie J ，Wang J ，Zeng R ，Zhou F ，Peng J ，Xie Y ... -  《Frontiers in Cell and Developmental Biology》

Identification of cuproptosis-related subtypes, construction of a prognosis model, and tumor microenvironment landscape in gastric cancer.

Cuproptosis is a novel identified regulated cell death (RCD), which is correlated with the development, treatment response and prognosis of cancer. However, the potential role of cuproptosis-related genes (CRGs) in the tumor microenvironment (TME) of gastric cancer (GC) remains unknown. Transcriptome profiling, somatic mutation, somatic copy number alteration and clinical data of GC samples were downloaded from the Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) database to describe the alterations of CRGs from genetic and transcriptional fields. Differential, survival and univariate cox regression analyses of CRGs were carried out to investigate the role of CRGs in GC. Cuproptosis molecular subtypes were identified by using consensus unsupervised clustering analysis based on the expression profiles of CRGs, and further analyzed by GO and KEGG gene set variation analyses (GSVA). Genes in distinct molecular subtypes were also analyzed by GO and KEGG gene enrichment analyses (GSEA). Differentially expressed genes (DEGs) were screened out from distinct molecular subtypes and further analyzed by GO enrichment analysis and univariate cox regression analysis. Consensus clustering analysis of prognostic DEGs was performed to identify genomic subtypes. Next, patients were randomly categorized into the training and testing group at a ratio of 1:1. CRG Risk scoring system was constructed through logistic least absolute shrinkage and selection operator (LASSO) cox regression analysis, univariate and multivariate cox analyses in the training group and validated in the testing and combined groups. Real-time quantitative polymerase chain reaction (RT-qPCR) was used to evaluate the expression of key Risk scoring genes. Sensitivity and specificity of Risk scoring system were examined by using receiver operating characteristic (ROC) curves. pRRophetic package in R was used to investigate the therapeutic effects of drugs in high- and low- risk score group. Finally, the nomogram scoring system was developed to predict patients' survival through incorporating the clinicopathological features and CRG Risk score. Most CRGs were up-regulated in tumor tissues and showed a relatively high mutation frequency. Survival and univariate cox regression analysis revealed that LIAS and FDX1 were significantly associated with GC patients' survival. After consensus unsupervised clustering analysis, GC patients were classified into two cuproptosis molecular subtypes, which were significantly associated with clinical features (gender, age, grade and TNM stage), prognosis, metabolic related pathways and immune cell infiltration in TME of GC. GO enrichment analyses of 84 DEGs, obtained from distinct molecular subtypes, revealed that DEGs primarily enriched in the regulation of metabolism and intracellular/extracellular structure in GC. Univariate cox regression analysis of 84 DEGs further screened out 32 prognostic DEGs. According to the expression profiles of 32 prognostic DEGs, patients were re-classified into two gene subtypes, which were significantly associated with patients' age, grade, T and N stage, and survival of patients. Nest, the Risk score system was constructed with moderate sensitivity and specificity. A high CRG Risk score, characterized by decreased microsatellite instability-high (MSI-H), tumor mutation burden (TMB) and cancer stem cell (CSC) index, and high stromal and immune score in TME, indicated poor survival. Four of five key Risk scoring genes expression were dysregulated in tumor compared with normal samples. Moreover, CRG Risk score was greatly related with sensitivity of multiple drugs. Finally, we established a highly accurate nomogram for promoting the clinical applicability of the CRG Risk scoring system. Our comprehensive analysis of CRGs in GC demonstrated their potential roles in TME, clinicopathological features, and prognosis. These findings may improve our understanding of CRGs in GC and provide new perceptions for doctors to predict prognosis and develop more effective and personalized therapy strategies.

Wang J ，Qin D ，Tao Z ，Wang B ，Xie Y ，Wang Y ，Li B ，Cao J ，Qiao X ，Zhong S ，Hu X ... -  《Frontiers in Immunology》