Crosstalk of angiogenesis-related subtypes, establishment of a prognostic signature and immune infiltration characteristics in colorectal adenocarcinoma.

Colorectal adenocarcinoma (COAD) is one of the most common malignancies and angiogenesis is vital to the development of cancer. Here, we explored the roles of angiogenesis-related genes (ARGs) that affect the prognosis of COAD and constructed risk models to assess patient prognosis, immune characteristics, and treatment outcomes. We comprehensively characterized the transcriptional and genetic modifications of 48 ARGs in COAD and evaluated the expression patterns. We identified two ARG subgroups using the consensus clustering algorithm. Based on the differentially expressed genes (DEGs) of two ARG subtypes, we calculated risk score, namely ARG_scores, and calssified COAD patients into different risk groups. To investigate the expression of ARG_score-related genes, qRT-PCR was performed. Subsequently, we mapped the nomogram to visually and accurately describe the value of the application of ARG_score. Finally, the correlation between ARG_score and clinical features, immune infiltration along with drug sensitivity were explored. We identified two ARG related subgroups and there were great differences in overall survival (OS) and tumor microenvironment. Then, we created an ARG_score for predicting overall survival based on eight DEGs and confirmed its reliable predictive power in COAD patients, with higher ARG_score associated with worse prognosis. Furthermore, eight ARG_score-related genes expression was investigated by qRT-PCR. To make the ARG_score clinically feasible, we created a highly reliable nomogram. We also found a higher proportion of microsatellite instability-high (MSI-H) and higher tumor mutational burden (TMB) in the high-risk group. In addition, ARG_score was notably correlated with cancer stem cell indices and drug sensitivity. This scoring model has potential clinical application value in the prognosis, immune microenvironment and therapeutic drug sensitivity of COAD, which provides new insights for personalized treatment.

Cui G ，Liu J ，Wang M ，Shon K ，Wang C ，Wei F ，Sun Z ... -  《Frontiers in Immunology》

Identification of angiogenesis-related subtypes, the development of a prognosis model, and features of tumor microenvironment in colon cancer.

Angiogenesis is associated with tumor progression, prognosis, and treatment effect. However, the angiogenesis' underlying mechanisms in the tumor microenvironment (TME) still remain unclear. Understanding the dynamic interactions between angiogenesis and TME in colon adenocarcinoma (COAD) is necessary. We downloaded the transcriptome data and corresponding clinical data of colon cancer patients from The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) databases, respectively. We identified two distinct angiogenesis-related molecular subtypes (subtype A and subtype B) and assessed the clinical features, prognosis, and infiltrating immune cells of patients in the two subtypes. According to the prognostic differential genes, we defined two different gene clusters to further explore the correlation between angiogenesis and tumor heterogeneity. Then, we construct the prognostic risk scoring model angiogenesis-related gene (ARG-score) including seven genes (ARMCX2, latent transforming growth factor β binding protein 1, ADAM8, FABP4, CCL11, CXCL11, ITLN1) using Lasso-multivariate cox method. We analyzed the correlation between ARG-score and prognosis, clinicopathological features, TME, molecular feature, cancer stem cells (CSCs), and microsatellite instability (MSI) status. To assess the application value of ARG-score in clinical treatment, immunophenotype score was used to predict patients' immunotherapy response in colon cancer. We found the mutations of ARGs in TCGA-COAD dataset from genetic levels and discussed their expression patterns based on TCGA and GEO datasets. We observed important differences in clinicopathological features, prognosis, immune feature, molecular feature between the two molecular subgroups. Then, we established an ARG-score for predicting OS and validated its predictive capability. A high ARG-score characterized by higher transcription level of ARGs, suggested lower MSI-high (MSI-H), lower immune score, and worse clinical stage and survival outcome. Additionally, the ARG-score was remarkably related to the CSCs index and immunotherapy sensitivity. We found two new molecular subtypes and two gene clusters based on ARGs and established an ARG-score. Multilayered analysis revealed that ARGs were remarkably correlated to the heterogeneity of colon cancer patients and explained the process of tumorigenesis and progression better. The ARG-score can help us better assess patients' survival outcomes and provide guidance for individualized treatment.

Wang F ，Wang C ，Li B ，Wang G ，Meng Z ，Han J ，Guo G ，Yu B ，Wang G ... -  《-》

Transcriptomic correlates of cell cycle checkpoints with distinct prognosis, molecular characteristics, immunological regulation, and therapeutic response in colorectal adenocarcinoma.

Wang H ，Wang W ，Wang Z ，Li X ... -  《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》

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》