A mitophagy-related gene signature associated with prognosis and immune microenvironment in colorectal cancer.

Colorectal cancer (CRC) is a heterogeneous disease and one of the most prevalent malignancies worldwide. Previous research has demonstrated that mitophagy is crucial to developing colorectal cancer. This study aims to examine the association between mitophagy-related genes and the prognosis of CRC patients. Gene expression profiles and clinical information of CRC patients were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Univariate Cox regression and the least absolute shrinkage and selection operator (LASSO) regression analysis were applied to establish a prognostic signature using mitophagy related genes. Kaplan-Meier and receiver operating characteristic (ROC) curves were used to analyze patient survival and predictive accuracy. Meanwhile, we also used the Genomics of Drug Sensitivity in Cancer (GDSC) database and Tumor Immune Dysfunction and Exclusion (TIDE) algorithm to estimate the sensitivity of chemotherapy, targeted therapy and immunotherapy. ATG14 overexpression plasmid was used to regulate the ATG14 expression level in HCT116 and SW480 cell lines, and cell counting kit-8, colony formation and transwell migration assay were performed to validate the function of ATG14 in CRC cells. A total of 22 mitophagy-driven genes connected with CRC survival were identified, and then a novel prognostic signature was established based on 10 of them (AMBRA1, ATG14, MAP1LC3A, MAP1LC3B, OPTN, VDAC1, ATG5, CSNK2A2, MFN1, TOMM22). Patients were divided into high-risk and low-risk groups based on the median risk score, and the survival of patients in the high-risk group was significantly shorter in both the training cohort and two independent cohorts. ROC curve showed that the area under the curves (AUC) of 1-, 3- and 5-year survival were 0.66, 0.66 and 0.64, respectively. Multivariate Cox regression analysis confirmed the independent prognostic value of the signature. Then we constructed a Nomogram combining the risk score, age and M stage, which had a concordance index of survival prediction of 0.77 (95% CI 0.71-0.83) and more robust predictive accuracy. Results showed that CD8+ T cells, regulatory T cells and activated NK cells were significantly more enriched in the high-risk group. Furthermore, patients in the high-risk group are more sensitive to targeted therapy or chemotherapy, including bosutinib, elesclomol, lenalidomide, midostaurin, pazopanib and sunitinib, while the low-risk group is more likely to benefit from immunotherapy. Finally, in vitro study confirmed the oncogenic significance of ATG14 in both HCT116 and SW480 cells, whose overexpression increased CRC cell proliferation, colony formation, and migration. In conclusion, we developed a novel mitophagy-related gene signature that can be utilized not only as an independent predictive biomarker but also as a tool for tailoring personalizing treatment for CRC patients, and we confirmed ATG14 as a novel oncogene in CRC.

Zhang C ，Zeng C ，Xiong S ，Zhao Z ，Wu G ... -  《Scientific Reports》

Development of a novel hypoxia-immune-related LncRNA risk signature for predicting the prognosis and immunotherapy response of colorectal cancer.

Colorectal cancer (CRC) is one of the most common digestive system tumors worldwide. Hypoxia and immunity are closely related in CRC; however, the role of hypoxia-immune-related lncRNAs in CRC prognosis is unknown. Data used in the current study were sourced from the Gene Expression Omnibus and The Cancer Genome Atlas (TCGA) databases. CRC patients were divided into low- and high-hypoxia groups using the single-sample gene set enrichment analysis (ssGSEA) algorithm and into low- and high-immune groups using the Estimation of STromal and Immune cells in MAlignant Tumours using Expression data (ESTIMATE) algorithm. Differentially expressed lncRNAs (DElncRNAs) between low- and high-hypoxia groups, low- and high-immune groups, and tumor and control samples were identified using the limma package. Hypoxia-immune-related lncRNAs were obtained by intersecting these DElncRNAs. A hypoxia-immune-related lncRNA risk signature was developed using univariate Cox regression and least absolute shrinkage and selection operator (LASSO) analyses. The tumor microenvironments in the low- and high-risk groups were evaluated using ssGSEA, ESTIMATE, and the expression of immune checkpoints. The therapeutic response in the two groups was assessed using TIDE, IPS, and IC50. A ceRNA network based on signature lncRNAs was constructed. Finally, we used RT-qPCR to verify the expression of hypoxia-immune-related lncRNA signatures in normal and cancer tissues. Using differential expression analysis, and univariate Cox and LASSO regression analyses, ZNF667-AS1, LINC01354, LINC00996, DANCR, CECR7, and LINC01116 were selected to construct a hypoxia-immune-related lncRNA signature. The performance of the risk signature in predicting CRC prognosis was validated in internal and external datasets, as evidenced by receiver operating characteristic curves. In addition, we observed significant differences in the tumor microenvironment and immunotherapy response between low- and high-risk groups and constructed a CECR7-miRNA-mRNA regulatory network in CRC. Furthermore, RT-qPCR results confirmed that the expression patterns of the six lncRNA signatures were consistent with those in TCGA-CRC cohort. Our study identified six hypoxia-immune-related lncRNAs for predicting CRC survival and sensitivity to immunotherapy. These findings may enrich our understanding of CRC and help improve CRC treatment. However, large-scale long-term follow-up studies are required for verification.

Luan L ，Dai Y ，Shen T ，Yang C ，Chen Z ，Liu S ，Jia J ，Li Z ，Fang S ，Qiu H ，Cheng X ，Yang Z ... -  《Frontiers in Immunology》

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》

Identification of prognostic immune-related gene signature associated with tumor microenvironment of colorectal cancer.

The tumor microenvironment (TME) has significantly correlation with tumor occurrence and prognosis. Our study aimed to identify the prognostic immune-related genes (IRGs)in the tumor microenvironment of colorectal cancer (CRC). Transcriptome and clinical data of CRC cases were downloaded from TCGA and GEO databases. Stromal score, immune score, and tumor purity were calculated by the ESTIMATE algorithm. Based on the scores, we divided CRC patients from the TCGA database into low and high groups, and the differentially expressed genes (DEGs) were identified. Immune-related genes (IRGs) were selected by venn plots. To explore underlying pathways, protein-protein interaction (PPI) networks and functional enrichment analysis were used. After utilizing LASSO Cox regression analysis, we finally established a multi-IRGs signature for predicting the prognosis of CRC patients. A nomogram consists of the thirteen-IRGs signature and clinical parameters was developed to predict the overall survival (OS). We investigated the association between prognostic validated IRGs and immune infiltrates by TIMER database. Gene expression profiles and clinical information of 1635 CRC patients were collected from the TCGA and GEO databases. Higher stromal score, immune score and lower tumor purity were observed positive correlation with tumor stage and poor OS. Based on stromal score, immune score and tumor purity, 1517 DEGs, 1296 DEGs, and 1892 DEGs were identified respectively. The 948 IRGs were screened by venn plots. A thirteen-IRGs signature was constructed for predicting survival of CRC patients. Nomogram with a C-index of 0.769 (95%CI, 0.717-0.821) was developed to predict survival of CRC patients by integrating clinical parameters and thirteen-IRGs signature. The AUC for 1-, 3-, and 5-year OS were 0.789, 0.783 and 0.790, respectively. Results from TIMER database revealed that CD1B, GPX3 and IDO1 were significantly related with immune infiltrates. In this study, we established a novel thirteen immune-related genes signature that may serve as a validated prognostic predictor for CRC patients, thus will be conducive to individualized treatment decisions.

Wang Y ，Li W ，Jin X ，Jiang X ，Guo S ，Xu F ，Su X ，Wang G ，Zhao Z ，Gu X ... -  《BMC CANCER》

Integrative analysis of the prognostic value and immune microenvironment of mitophagy-related signature for multiple myeloma.

Multiple myeloma (MM) is a fatal malignant tumor in hematology. Mitophagy plays vital roles in the pathogenesis and drug sensitivity of MM. We acquired transcriptomic expression data and clinical index of MM patients from NCI public database, and 36 genes involved in mitophagy from the gene set enrichment analysis (GSEA) database. Least absolute shrinkage and selection operator (LASSO) Cox regression analysis was conducted to construct a risk score prognostic model. Kaplan-Meier survival analysis and receiver operation characteristic curves (ROC) were conducted to identify the efficiency of prognosis and diagnosis. ESTIMATE algorithm and immune-related single-sample gene set enrichment analysis (ssGSEA) was performed to uncover the level of immune infiltration. QRT-PCR was performed to verify gene expression in clinical samples of MM patients. The sensitivity to chemotherapy drugs was evaluated upon the database of the genomics of drug sensitivity in cancer (GDSC). Fifty mitophagy-related genes were differently expressed in two independent cohorts. Ten out of these genes were identified to be related to MM overall survival (OS) rate. A prognostic risk signature model was built upon on these genes: VDAC1, PINK1, VPS13C, ATG13, and HUWE1, which predicted the survival of MM accurately and stably both in training and validation cohorts. MM patients suffered more adverse prognosis showed more higher risk core. In addition, the risk score was considered as an independent prognostic element for OS of MM patients by multivariate cox regression analysis. Functional pathway enrichment analysis of differentially expressed genes (DEGs) based on risk score showed terms of cell cycle, immune response, mTOR pathway, and MYC targets were obviously enriched. Furthermore, MM patients with higher risk score were observed lower immune scores and lower immune infiltration levels. The results of qRT-PCR verified VDAC1, PINK1, and HUWE1 were dysregulated in new diagnosed MM patients. Finally, further analysis indicated MM patients showed more susceptive to bortezomib, lenalidomide and rapamycin in high-risk group. Our research provided a neoteric prognostic model of MM based on mitophagy genes. The immune infiltration level based on risk score paved a better understanding of the participation of mitophagy in MM.

Jia Y ，Liu R ，Shi L ，Feng Y ，Zhang L ，Guo N ，He A ，Kong G ... -  《BMC CANCER》