The Interferon Gamma-Related Long Noncoding RNA Signature Predicts Prognosis and Indicates Immune Microenvironment Infiltration in Colon Adenocarcinoma.

Colon adenocarcinoma (COAD) is one of the most common clinically malignant tumours of the digestive system, with high incidence and mortality and poor prognosis. Interferon-gamma (IFN-γ) and long noncoding RNAs (lncRNAs) have prognostic values and were closely associated with immune microenvironment in COAD. Thus, identifying IFN-γ-related lncRNAs may be valuable in predicting the survival of patients with COAD. In this study, we identified IFN-γ-related lncRNAs and divided COAD patients from the Cancer Genome Atlas (TCGA) database into training and validation sets. Pearson's correlation analysis and least absolute shrinkage and selection operator (LASSO) Cox regression were performed to select IFN-γ-related lncRNA-associated prognoses. Thirteen lncRNAs (AC025165.8, AC091633.3, FENDRR, LINC00882, LINC01828, LINC01829, MYOSLID, RP11-154H23.4, RP11-20J15.3, RP11-324L17.1, RP11-342A23.2, RP11-805I24.3, SERTAD4-AS1) were identified to construct an IFN-γ-related lncRNA prognostic signature in TCGA training (n =213) and validation (n =213) cohorts. COAD patient risk scores were calculated and classified into high- and low-risk groups based on the median value of the risk scores in each dataset. We compared the overall survival (OS) of patients stratified by age, gender, and stage. The OS in the high-risk group was significantly shorter than that in the low-risk group. In addition, the clinical nomogram incorporating the prognostic signature and clinical features showed a high concordance index of 0.78 and accurately predicted 1-, 3-, and 5-year survival times among COAD patients in the high- and low-risk groups. Based on the risk model, the high- and low-risk groups exhibited distinct differences in the immune system by gene set enrichment analysis (GSEA) functional annotation, and differentially expressed genes (DEGs) between the high- and low-risk groups were subjected to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. We investigated the expression of multiple immune checkpoint genes in the high- and low-risk groups and plotted Kaplan-Meier survival curves, indicating that immune checkpoint genes, such as LAG3 and PD. L1, STING and TIM 3, were also expressed differently between the two risk groups. Subsequently, there were dramatic differences in mutated genes, SNV (single nucleotide variants) classes, variant types and variant allele frequencies between low- and high-risk patients with COAD. Patients stratified by risk scores had different sensitivities to common chemotherapeutic agents. Finally, we used quantitative real-time polymerase chain reaction (qRT-PCR) assays to demonstrate that three lncRNAs were significantly differentially expressed in COAD tissues and adjacent normal tissues. Considered together, a thirteen-lncRNA prognostic signature has great potential to be a prognostic biomarker and could play an essential role in the immune microenvironment of COAD.

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

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》

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》

Construction and validation of a cuproptosis-related lncRNA signature as a novel and robust prognostic model for colon adenocarcinoma.

Cuproptosis, a newly identified form of programmed cell death, is thought to play a role in tumorigenesis. Long non-coding RNAs (lncRNAs) are reported to be associated with tumor progression and prognosis in colon adenocarcinoma (COAD). However, the role and prognostic value of cuproptosis-related lncRNAs in COAD remains unknown. This study is devoted to constructing and validating a cuproptosis-related lncRNA signature that can predict COAD patient outcomes using bioinformatics methods. The COAD mRNA and lncRNA expression profiles and corresponding clinical data were downloaded from The Cancer Genome Atlas (TCGA) database and 2,567 cuproptosis-related lncRNAs were obtained. A 10 cuproptosis-related-lncRNA prognostic signature was then constructed using the least absolute shrinkage and selection operator (LASSO) algorithm and Cox regression model and patients were divided into high- and low-risk groups. Kaplan-Meier analysis, receiver operating characteristic (ROC) curve, and a nomogram were employed to evaluate the predictive power of the signature. The immune characteristics and drug sensitivity were also investigated based on the signature. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) was performed to verify the risk model. In vitro experiments were conducted to validate the expression of the ten lncRNAs during cuproptosis. The high-risk group was associated with shorter overall survival (OS) time in COAD patients (p<0.001). Multivariate Cox regression indicated that a high-risk score was an independent risk factor for poor prognosis (p<0.001). ROC curve analysis was performed to confirm the validity of the signature (area under the curve (AUC) at 3 years: 0.879). Gene Ontology (GO) enrichment analysis revealed that the signature was highly correlated with the immune response in biological processes. The immune function, the score of the immune cells, and the expression of immune checkpoints were significantly different between the two risk groups. Three drugs, LAQ824, FH535, YM155, were found to be more sensitive in the high-risk group. Finally, the expression levels of the ten lncRNAs comprising the signature were tested by qRT-PCR. A ten-cuproptosis-related lncRNA signature was constructed that provided promising insights into personalized prognosis and drug selection among COAD patients.

Xu M ，Mu J ，Wang J ，Zhou Q ，Wang J ... -  《Frontiers in Oncology》