scholarly journals Identification of a novel autophagy signature for predicting survival in patients with lung adenocarcinoma

PeerJ ◽  
2021 ◽  
Vol 9 ◽  
pp. e11074
Author(s):  
Jin Duan ◽  
Youming Lei ◽  
Guoli Lv ◽  
Yinqiang Liu ◽  
Wei Zhao ◽  
...  

Background Lung adenocarcinoma (LUAD) is the most commonhistological lung cancer subtype, with an overall five-year survivalrate of only 17%. In this study, we aimed to identify autophagy-related genes (ARGs) and develop an LUAD prognostic signature. Methods In this study, we obtained ARGs from three databases and downloaded gene expression profiles from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database. We used TCGA-LUAD (n = 490) for a training and testing dataset, and GSE50081 (n = 127) as the external validation dataset.The least absolute shrinkage and selection operator (LASSO) Cox and multivariate Cox regression models were used to generate an autophagy-related signature. We performed gene set enrichment analysis (GSEA) and immune cell analysis between the high- and low-risk groups. A nomogram was built to guide the individual treatment for LUAD patients. Results We identified a total of 83 differentially expressed ARGs (DEARGs) from the TCGA-LUAD dataset, including 33 upregulated DEARGs and 50 downregulated DEARGs, both with thresholds of adjusted P < 0.05 and |Fold change| > 1.5. Using LASSO and multivariate Cox regression analyses, we identified 10 ARGs that we used to build a prognostic signature with areas under the curve (AUCs) of 0.705, 0.715, and 0.778 at 1, 3, and 5 years, respectively. Using the risk score formula, the LUAD patients were divided into low- or high-risk groups. Our GSEA results suggested that the low-risk group were enriched in metabolism and immune-related pathways, while the high-risk group was involved in tumorigenesis and tumor progression pathways. Immune cell analysis revealed that, when compared to the high-risk group, the low-risk group had a lower cell fraction of M0- and M1- macrophages, and higher CD4 and PD-L1 expression levels. Conclusion Our identified robust signature may provide novel insight into underlying autophagy mechanisms as well as therapeutic strategies for LUAD treatment.

2021 ◽  
Author(s):  
Meiling Jin ◽  
Diangeng Li

Abstract BackgroundPapillary renal cell carcinoma (PRCC) is a common renal cell carcinoma. Recent studies have reported that ferroptosis is involved in the occurrence and development of tumors. Long non-coding RNAs could be used as independent biomarkers for the diagnosis and prognosis of a variety of tumors, and many lncRNAs are related to the pathogenesis of PRCC. However, there are few studies on the ferroptosis-related lncRNAs of PRCC. This study aimed to establish ferroptosis-related lncRNAs prognostic signature in patients with PRCC.MethodsGene expression profile and clinical information of patients with PRCC were obtained from The Cancer Genome Atlas (TCGA) database. Lasso-Penalzed Cox regression and univariate Cox regression analysis were utilized for model construction. The Kaplan-Meier (K-M) and Receiver Operating Characteristic (ROC) curves were plotted to validate the predictive effect of the prognostic signature. Immune cell infiltration and immune function were compared between the high-risk and low-risk groups. Chemotherapy sensitivity analysis was also performed. ResultsWe constructed a prognostic signature consisted of 15 ferroptosis-related lncRNAs. The Kaplan-Meier curves validated the fine predictive accuracy of the prognostic signature (p < 0.001). The area under the curve (AUC) of the lncRNAs signature was 0.930, exhibiting robust prognostic capacity. The high-risk group had a greater degree of immune cell infiltration, such as B cells, T cell CD8, macrophages, NK cell, etc., compared with the low-risk group. There were significant differences in inflammation-promoting, parainflammation and Type I IFN reponse between the low-risk and high-risk groups. The expressions of immune checkpoints including CD80, IDO1, LAG3, etc. were significantly higher in high-risk group. Chemotherapy sensitivity analysis showed that MNX1-AS1, ZFAS1, MIR4435-2HG and ADAMTS9-AS1 were significantly correlated with the sensitivity of some chemotherapy drugs. ConclusionWe demonstrated that a ferroptosis-related lncRNAs prognostic signature could be a novel biomarker for PRCC. Our findings could provide a new insight for immune research and treatment strategies for patients with PRCC.


2020 ◽  
Vol 2020 ◽  
pp. 1-15
Author(s):  
Jianfeng Zheng ◽  
Benben Cao ◽  
Xia Zhang ◽  
Zheng Niu ◽  
Jinyi Tong

Cervical cancer (CC) is a common gynecological malignancy for which prognostic and therapeutic biomarkers are urgently needed. The signature based on immune-related lncRNAs (IRLs) of CC has never been reported. This study is aimed at establishing an IRL signature for patients with CC. A cohort of 326 CC and 21 normal tissue samples with corresponding clinical information was included in this study. Twenty-eight IRLs were collected according to the Pearson correlation analysis between the immune score and lncRNA expression ( p < 0.01 ). Four IRLs (BZRAP1-AS1, EMX2OS, ZNF667-AS1, and CTC-429P9.1) with the most significant prognostic values ( p < 0.05 ) were identified which demonstrated an ability to stratify patients into the low-risk and high-risk groups by developing a risk score model. It was observed that patients in the low-risk group showed longer overall survival (OS) than those in the high-risk group in the training set, valid set, and total set. The area under the curve (AUC) of the receiver operating characteristic curve (ROC curve) for the four-IRL signature in predicting the one-, two-, and three-year survival rates was larger than 0.65. In addition, the low-risk and high-risk groups displayed different immune statuses in GSEA. These IRLs were also significantly correlated with immune cell infiltration. Our results showed that the IRL signature had a prognostic value for CC. Meanwhile, the specific mechanisms of the four IRLs in the development of CC were ascertained preliminarily.


2020 ◽  
Vol 13 (1) ◽  
Author(s):  
Yinglian Pan ◽  
Li Ping Jia ◽  
Yuzhu Liu ◽  
Yiyu Han ◽  
Qian Li ◽  
...  

Abstract Background In this study we aimed to identify a prognostic signature in BRCA1/2 mutations to predict disease progression and the efficiency of chemotherapy ovarian cancer (OV), the second most common cause of death from gynecologic cancer in women worldwide. Methods Univariate Cox proportional-hazards and multivariate Cox regression analyses were used to identifying prognostic factors from data obtained from The Cancer Genome Atlas (TCGA) database. The area under the curve of the receiver operating characteristic curve was assessed, and the sensitivity and specificity of the prediction model were determined. Results A signature consisting of two long noncoding RNAs(lncRNAs), Z98885.2 and AC011601.1, was selected as the basis for classifying patients into high and low-risk groups (median survival: 7.2 years vs. 2.3 years). The three-year overall survival (OS) rates for the high- and low-risk group were approximately 38 and 100%, respectively. Chemotherapy treatment survival rates indicated that the high-risk group had significantly lower OS rates with adjuvant chemotherapy than the low-risk group. The one-, three-, and five-year OS were 100, 40, and 15% respectively in the high-risk group. The survival rate of the high-risk group declined rapidly after 2 years of OV chemotherapy treatment. Multivariate Cox regression associated with other traditional clinical factors showed that the 2-lncRNA model could be used as an independent OV prognostic factor. Analyses of data from the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) indicated that these signatures are pivotal to cancer development. Conclusion In conclusion, Z98885.2 and AC011601.1 comprise a novel prognostic signature for OV patients with BRCA1/2 mutations, and can be used to predict prognosis and the efficiency of chemotherapy.


2021 ◽  
Vol 11 ◽  
Author(s):  
Xiaoxia Tong ◽  
Xiaofei Qu ◽  
Mengyun Wang

BackgroundCutaneous melanoma (CM) is one of the most aggressive cancers with highly metastatic ability. To make things worse, there are limited effective therapies to treat advanced CM. Our study aimed to investigate new biomarkers for CM prognosis and establish a novel risk score system in CM.MethodsGene expression data of CM from Gene Expression Omnibus (GEO) datasets were downloaded and analyzed to identify differentially expressed genes (DEGs). The overlapped DEGs were then verified for prognosis analysis by univariate and multivariate COX regression in The Cancer Genome Atlas (TCGA) datasets. Based on the gene signature of multiple survival associated DEGs, a risk score model was established, and its prognostic and predictive role was estimated through Kaplan-Meier (K-M) analysis and log-rank test. Furthermore, the correlations between prognosis related genes expression and immune infiltrates were analyzed via Tumor Immune Estimation Resource (TIMER) site.ResultsA total of 103 DEGs were obtained based on GEO cohorts, and four genes were verified in TCGA datasets. Subsequently, four genes (ADAMDEC1, GNLY, HSPA13, and TRIM29) model was developed by univariate and multivariate Cox regression analyses. The K-M plots showed that the high-risk group was associated with shortened survival than that in the low-risk group (P &lt; 0.0001). Multivariate analysis suggested that the model was an independent prognostic factor (high-risk vs. low-risk, HR= 2.06, P &lt; 0.001). Meanwhile, the high-risk group was prone to have larger breslow depth (P&lt; 0.001) and ulceration (P&lt; 0.001).ConclusionsThe four-gene risk score model functions well in predicting the prognosis and treatment response in CM and will be useful for guiding therapeutic strategies for CM patients. Additional clinical trials are needed to verify our findings.


2021 ◽  
Author(s):  
Yali Zhong ◽  
Xiaobin Luo ◽  
Fubing Yang ◽  
Xinling Song

Abstract Object: Immune related genes play an important role in the process of tumor genesis and development. Therefore, we aim to find the Immune genes which are related to the prognosis of glioma patients, and to explore the infiltration of Immune cells in glioma microenvironment. Methods We downloaded the data of the glioma samples from the CGGA database, and performed batch correction to screen the primary glioma samples for subsequent analysis. Then the ESTIMATE algorithm was used to deal with the Stromal scores and Immune scores of the primary glioma samples, and the difference was analyzed. Then the common Immune related genes (IRGs) were obtained by intersecting with the Immune genes in the ImmPort database. Moreover, we used common IRGs to construct protein-protein interaction (PPI) networks, from which we screened the top 30 genes with high connectivity, and Lasso regression was used to screen the IRGs. Lastly, we obtained the combined genes, which were overlapped both in the top 30 high-connection genes and Lasso regression genes. The final genes were used to construct COX risk prediction models. The accuracy of the model were verified by the TCGA glioma data, and the model genes were analyzed for Immune-related pathways, as well as the Hallmark and KEGG enrichment. Additionally, we used CIBERSOFT algorithm to estimate the Immune cell content of the samples, and analyzed the differences, correlations and survival of the Immune cells in high and low risk groups. Results Firstly, a total of 117 IRGs were obtained from the gene sets, which were overlapped in the data of Stromal score, Immune score and ImmPort database. Secondly, the top 30 genes were selected after the PPI network, and another 26 genes were screened out after the Lasso regression algorithm. And then, six coexist IRGs were obtained from the intersecting sets. Furthermore, the COX risk prediction model was constructed and tested, showing that the overall survival rate of the high-risk group was about 50% of that of the low-risk group. We observed that the high-risk group were enriched in Immune response and Immune process. Most importantly, in KEGG pathways, the high-risk groups were mainly enriched in p53 signaling pathway, JAK-STAT signaling pathway, pathways in cancer and cell cycle. By estimating the Immune cell contents, we also found that the Immune cell Plasma cells, T cells CD8, T cells CD4 naïve, T cells regulatory (Tregs), Macrophages M0 and Neutrophils were higher in high-risk groups, when compared to the low-risk group, with significant difference. Finally, the correlation analysis showed that the degree of Immune infiltration in high-risk groups was related to T cells regulatory (Tregs), Macrophages M0 and Neutrophils. Conclusion A COX risk prediction model of 6 genes was successfully constructed, which was enriched in Immune-related pathways. Meanwhile, survival analysis and TCGA data validation revealed significant differences in the model genes in the overall survival of the glioma patients, and the degree of Immune infiltration in the model was associated with T cells regulatory (Tregs), Macrophages M0 and Neutrophils.


2020 ◽  
Author(s):  
Kui Wu ◽  
Yongjie Shui ◽  
Wenzheng Sun ◽  
Sheng Lin ◽  
Haowen Pang

Abstract Objective This study aimed to develop and validate the combination of radiomic features and clinical characteristics that can predict patient survival in HCC with PVTT treated with SBRT. Materials and Methods The prediction model was developed in a primary cohort of 70 patients with HCC and PVTT treated with SBRT, using data acquired between December 2015 and June 2017. The radiomic features were extracted from computed tomography (CT) scans. A least absolute shrinkage and selection operator regression model was used to build the radiomic feature. Multivariate Cox-regression hazard models were created for analyzing survival outcomes and the radiomic features and clinical characteristics were presented with a nomogram. The area under the curve (AUC) of the receiver operating characteristic curve was used to evaluate the model. Participants were divided into a high-risk group and a low-risk group based on the radiomic features. Results A total of seven radiomic features and five clinical characteristics were extracted for survival analysis. A combination of the radiomic features and clinical characteristics resulted in better performance for the estimation of overall survival (OS) [AUC = 0.859 (CI: 0.770–0.948)] than that with clinical characteristics alone [AUC = 0.761 (CI: 0.641–0.881)]. These patients were divided into high-risk and low-risk groups according to the radiomic features. Conclusion This study demonstrated that a nomogram of combined radiomic features and clinical characteristics can be conveniently used to facilitate individualized preoperative prediction of OS in patients with HCC with PVTT before SBRT.


2021 ◽  
Vol 2021 ◽  
pp. 1-20
Author(s):  
Zhicheng Zhuang ◽  
Huajun Cai ◽  
Hexin Lin ◽  
Bingjie Guan ◽  
Yong Wu ◽  
...  

Background. Pyroptosis has been confirmed as a type of inflammatory programmed cell death in recent years. However, the prognostic role of pyroptosis in colon cancer (CC) remains unclear. Methods. Dataset TCGA-COAD which came from the TCGA portal was taken as the training cohort. GSE17538 from the GEO database was treated as validation cohorts. Differential expression genes (DEGs) between normal and tumor tissues were confirmed. Patients were classified into two subgroups according to the expression characteristics of pyroptosis-related DEGs. The LASSO regression analysis was used to build the best prognostic signature, and its reliability was validated using Kaplan–Meier, ROC, PCA, and t-SNE analyses. And a nomogram based on the multivariate Cox analysis was developed. The enrichment analysis was performed in the GO and KEGG to investigate the potential mechanism. In addition, we explored the difference in the abundance of infiltrating immune cells and immune microenvironment between high- and low-risk groups. And we also predicted the association of common immune checkpoints with risk scores. Finally, we verified the expression of the pyroptosis-related hub gene at the protein level by immunohistochemistry. Results. A total of 23 pyroptosis-related DEGs were identified in the TCGA cohort. Patients were classified into two molecular clusters (MC) based on DEGs. Kaplan–Meier survival analysis indicated that patients with MC1 represented significantly poorer OS than patients with MC2. 13 overall survival- (OS-) related DEGs in MCs were used to construct the prognostic signature. Patients in the high-risk group exhibited poorer OS compared to those in the low-risk group. Combined with the clinical features, the risk score was found to be an independent prognostic factor of CC patients. The above results are verified in the external dataset GSE17538. A nomogram was established and showed excellent performance. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses indicated that the varied prognostic performance between high- and low-risk groups may be related to the immune response mediated by local inflammation. Further analysis showed that the high-risk group has stronger immune cell infiltration and lower tumor purity than the low-risk group. Through the correlation between risk score and immune checkpoint expression, T-cell immunoglobulin and mucin domain-containing protein 3 (TIM-3) was predicted as a potential therapeutic target for the high-risk group. Conclusion. The 13-gene signature was associated with OS, immune cells, tumor purity, and immune checkpoints in CC patients, and it could provide the basis for immunotherapy and predicting prognosis and help clinicians make decisions for individualized treatment.


2021 ◽  
Author(s):  
Fang Wen ◽  
Xiaoxue Chen ◽  
Wenjie Huang ◽  
Shuai Ruan ◽  
Suping Gu ◽  
...  

Abstract Background: The diagnosis rate and mortality of gastric cancer (GC) are among the highest in the global, so it is of great significance to predict the survival time of GC patients. Ferroptosis and iron-metabolism make a critical impact on tumor development and are closely linked to the treatment of cancer and the prognosis of patients. However, the predictive value of the genes involved in ferroptosis and iron-metabolism in GC and their effects on immune microenvironment remain to be further clarified.Methods: In this study, the RNA sequence information and general clinical indicators of GC patients were acquired from the public databases. We first systematically screen out 134 DEGs and 13 PRGs related to ferroptosis and iron-metabolism. Then, we identified six PRDEGs (GLS2, MTF1, SLC1A5, SP1, NOX4, and ZFP36) based on the LASSO-penalized Cox regression analysis. The 6-gene prognostic risk model was established in the TCGA cohort and the GC patients were separated into the high- and the low-risk groups through the risk score median value. GEO cohort was used for verification. The expression of PRDEGs was verified by quantitative QPCR.Results: Our study demonstrated that patients in the low-risk group had a higher survival probability compared with those in high-risk group. In addition, univariate and multivariate Cox regression analyses confirmed that the risk score was an independent prediction parameter. The ROC curve analysis and nomogram manifested that the risk model had the high predictive ability and was more sensitive than general clinical features. Furthermore, compared with the high-risk group, the low-risk group had higher TMB and a longer 5-year survival period. In the immune microenvironment of GC, there were also differences in immune function and highly infiltrated immune cells between the two risk groups.Conclusions: The prognostic risk model based on the six genes associated with ferroptosis and iron-metabolism has a good performance for predicting the prognosis of patients with GC. The treatment of cancer by inducing tumor ferroptosis or mediating tumor iron-metabolism, especially combined with immunotherapy, provides a new possibility for individualized treatment of GC patients.


2022 ◽  
Author(s):  
Cong Zhang ◽  
Cailing Zeng ◽  
Shaoquan Xiong ◽  
Zewei Zhao ◽  
Guoyu Wu

Abstract Background: Colorectal cancer (CRC) is a heterogeneous disease and one of the most common malignancies in the world. Previous studies have found that mitophagy plays an important role in the progression of colorectal cancer. This study is aimed to investigate the relationship between mitophagy-related genes and the prognosis of patients with CRC.Methods: 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 used to establish the prognostic signature composed of mitophagy related genes. Kaplan-Meier curve and receiver operating characteristic (ROC) curve were used to analyze patient survival and verify the predictive accuracy of the signature, respectively. Construction of a nomogram prognostic prediction model was based on risk scores and clinicopathological parameters. Using 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. Results: A total of 44 mitophagy-driven genes connected with CRC survival were identified, and prognostic signature was established based on the expression of 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 than that of the low-risk group among the TCGA cohort (median OS 67.3 months vs not reached, p=0.00059) and two independent cohorts from GEO (median OS in GSE17536: 54.0 months vs not reached, p=0.0082; in GSE245: 7.7 months vs not reached, p=0.025). 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 sensitivity and specificity. Results showed that CD8+ T cells, regulatory T cells and activated NK cells were significantly more abundant in the high-risk group. Furthermore, patients in the high-risk group were more sensitive to potential targeted therapies, including Motesanib, ATRA, Olaparib, Selumetinib, AZD8055 and immunotherapy. Conclusion: In conclusion, we constructed and validated a novel mitophagy related gene signature that can be used as an independent prognostic biomarker for CRC, and may lead to better stratification and selection of precise treatment for CRC patients.


2020 ◽  
Vol 10 ◽  
Author(s):  
Qiongxuan Fang ◽  
Hongsong Chen

BackgroundHepatocellular carcinoma (HCC) is the seventh most common malignancy and the second most common cause of cancer-related deaths. Autophagy plays a crucial role in the development and progression of HCC.MethodsUnivariate and Lasso Cox regression analyses were performed to determine a gene model that was optimal for overall survival (OS) prediction. Patients in the GSE14520 and GSE54236 datasets of the Cancer Genome Atlas (TCGA) were divided into the high-risk and low-risk groups according to established ATG models. Univariate and multivariate Cox regression analyses were used to identify risk factors for OS for the purpose of constructing nomograms. Calibration and receiver operating characteristic (ROC) curves were used to evaluate model performance. Real-time PCR was used to validate the effects of the presence or absence of an autophagy inhibitor on gene expression in HepG2 and Huh7 cell lines.ResultsOS in the high-risk group was significantly shorter than that in the low-risk group. Gene set enrichment analysis (GSEA) indicated that the association between the low-risk group and autophagy- as well as immune-related pathways was significant. ULK2, PPP3CC, and NAFTC1 may play vital roles in preventing HCC progression. Furthermore, tumor environment analysis via ESTIMATION indicated that the low-risk group was associated with high immune and stromal scores. Based on EPIC prediction, CD8+ T and B cell fractions in the TCGA and GSE54236 datasets were significantly higher in the low-risk group than those in the high-risk group. Finally, based on the results of univariate and multivariate analyses three variables were selected for nomogram development. The calibration plots showed good agreement between nomogram prediction and actual observations. Inhibition of autophagy resulted in the overexpression of genes constituting the gene model in HepG2 and Huh7 cells.ConclusionsThe current study determined the role played by autophagy-related genes (ATGs) in the progression of HCC and constructed a novel nomogram that predicts OS in HCC patients, through a combined analysis of TCGA and gene expression omnibus (GEO) databases.


Sign in / Sign up

Export Citation Format

Share Document