scholarly journals Development and validation of a novel 3-gene prognostic model for pancreatic adenocarcinoma based on ferroptosis-related genes

2022 ◽  
Vol 22 (1) ◽  
Author(s):  
Jihua Yang ◽  
XiaoHong Wei ◽  
Fang Hu ◽  
Wei Dong ◽  
Liao Sun

Abstract Background Molecular markers play an important role in predicting clinical outcomes in pancreatic adenocarcinoma (PAAD) patients. Analysis of the ferroptosis-related genes may provide novel potential targets for the prognosis and treatment of PAAD. Methods RNA-sequence and clinical data of PAAD was downloaded from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) public databases. The PAAD samples were clustered by a non-negative matrix factorization (NMF) algorithm. The differentially expressed genes (DEGs) between different subtypes were used by “limma_3.42.2” package. The R software package clusterProfiler was used for functional enrichment analysis. Then, a multivariate Cox proportional and LASSO regression were used to develop a ferroptosis-related gene signature for pancreatic adenocarcinoma. A nomogram and corrected curves were constructed. Finally, the expression and function of these signature genes were explored by qRT-PCR, immunohistochemistry (IHC) and proliferation, migration and invasion assays. Results The 173 samples were divided into 3 categories (C1, C2, and C3) and a 3-gene signature model (ALOX5, ALOX12, and CISD1) was constructed. The prognostic model showed good independent prognostic ability in PAAD. In the GSE62452 external validation set, the molecular model also showed good risk prediction. KM-curve analysis showed that there were significant differences between the high and low-risk groups, samples with a high-risk score had a worse prognosis. The predictive efficiency of the 3-gene signature-based nomogram was significantly better than that of traditional clinical features. For comparison with other models, that our model, with a reasonable number of genes, yields a more effective result. The results obtained with qPCR and IHC assays showed that ALOX5 was highly expressed, whether ALOX12 and CISD1 were expressed at low levels in tissue samples. Finally, function assays results suggested that ALOX5 may be an oncogene and ALOX12 and CISD1 may be tumor suppressor genes. Conclusions We present a novel prognostic molecular model for PAAD based on ferroptosis-related genes, which serves as a potentially effective tool for prognostic differentiation in pancreatic cancer patients.

2021 ◽  
Author(s):  
Jihua Yang ◽  
Xiaohong Wei ◽  
Fang Hu ◽  
Dong Wei ◽  
Liao Sun

Abstract Background: Molecular markers play an important role in predicting clinical outcomes in pancreatic adenocarcinoma (PAAD) patients. Analysis of the ferroptosis-related genes may provide novel potential targets for the prognosis and treatment of PAAD. Methods: RNA-sequence data of PAAD was downloaded from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) public databases. The PAAD samples were clustered by a non-negative matrix factorization (NMF) algorithm. The R software package clusterProfiler was used for functional enrichment analysis. By LASSO regression, we established a multi-gene prognostic model, and the model was validated by qPCR and immunohistochemistry and T3M4 cell line. Results: Three molecular subtypes were identified, and a 3-gene signature model (ALOX5, ALOX12, and CISD1) was constructed. The prognostic model showed good independent prognostic ability in PAAD. In the GSE62452 external validation set, the molecular model also showed good risk prediction. The predictive efficiency of the 3-gene signature-based nomogram was significantly better than that of traditional clinical features. Finally, qPCR and immunohistochemical staining and cell function results suggested that ALOX5 may be an oncogene and ALOX12 and CISD1 may be tumor suppressor genes. Conclusions: We present a novel prognostic molecular model for PAAD based on ferroptosis-related genes, which serves as a potentially effective tool for prognostic differentiation in pancreatic cancer patients.


2020 ◽  
Author(s):  
Chuxiang Lei ◽  
Wenlin Chen ◽  
Yuekun Wang ◽  
Binghao Zhao ◽  
Penghao Liu ◽  
...  

Abstract Background. Glioblastoma (GBM) is the most common primary malignant intracranial tumor and is closely related to metabolic alterations. However, few accepted prognostic models are currently available, especially models based on metabolic genes. Methods . Transcriptome data were obtained for all patients diagnosed with GBM from the Gene Expression Omnibus (GEO) (training cohort, n=369) and The Cancer Genome Atlas (TCGA) (validation cohort, n=152) with the following variables: age at diagnosis, sex, follow-up and overall survival (OS). Metabolic genes according to Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were filtered, and a Lasso regression model was constructed. Survival was assessed by univariate or multivariate Cox proportional hazards regression and Kaplan-Meier analysis, and we also conducted an independent external validation to examine the model. Results. There were 341 metabolic genes that showed significant differences between normal brain tissues and GBM tissues in both the training and validation cohorts, among which 56 genes were significantly correlated with the OS of patients. Lasso regression revealed that the metabolic prognostic model was composed of 18 genes, including COX10 , COMT , and GPX2 , with protective effects, as well as OCRL and RRM2 , with unfavorable effects. Patients classified as high-risk by the risk score from this model had markedly shorter OS than low-risk patients ( P <0.0001), and this significant result was also observed in the independent external validation cohort ( P <0.001). Conclusions . The prognosis of GBM was dramatically related to metabolic pathways, and our metabolic prognostic model had high accuracy and application value in predicting the OS of GBM patients. Background. Glioblastoma (GBM) is the most common primary malignant intracranial tumor and is closely related to metabolic alterations. However, few accepted prognostic models are currently available, especially models based on metabolic genes. Methods . Transcriptome data were obtained for all patients diagnosed with GBM from the Gene Expression Omnibus (GEO) (training cohort, n=369) and The Cancer Genome Atlas (TCGA) (validation cohort, n=152) with the following variables: age at diagnosis, sex, follow-up and overall survival (OS). Metabolic genes according to Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were filtered, and a Lasso regression model was constructed. Survival was assessed by univariate or multivariate Cox proportional hazards regression and Kaplan-Meier analysis, and we also conducted an independent external validation to examine the model. Results. There were 341 metabolic genes that showed significant differences between normal brain tissues and GBM tissues in both the training and validation cohorts, among which 56 genes were significantly correlated with the OS of patients. Lasso regression revealed that the metabolic prognostic model was composed of 18 genes, including COX10 , COMT , and GPX2 , with protective effects, as well as OCRL and RRM2 , with unfavorable effects. Patients classified as high-risk by the risk score from this model had markedly shorter OS than low-risk patients ( P <0.0001), and this significant result was also observed in the independent external validation cohort ( P <0.001).Conclusions . The prognosis of GBM was dramatically related to metabolic pathways, and our metabolic prognostic model had high accuracy and application value in predicting the OS of GBM patients.


2022 ◽  
Vol 2022 ◽  
pp. 1-16
Author(s):  
Jin Zhou ◽  
Zheming Liu ◽  
Huibo Zhang ◽  
Tianyu Lei ◽  
Jiahui Liu ◽  
...  

Purpose. Recent researches showed the vital role of BACH1 in promoting the metastasis of lung cancer. We aimed to explore the value of BACH1 in predicting the overall survival (OS) of early-stage (stages I-II) lung adenocarcinoma. Patients and Methods. Lung adenocarcinoma cases were screened from the Cancer Genome Atlas (TCGA) database. Functional enrichment analysis was performed to obtain the biological mechanisms of BACH1. Gene set enrichment analysis (GSEA) was performed to identify the difference of biological pathways between high- and low-BACH1 groups. Univariate and multivariate COX regression analysis had been used to screen prognostic factors, which were used to establish the BACH1 expression-based prognostic model in the TCGA dataset. The C-index and time-dependent AUC curve were used to evaluate predictive power of the model. External validation of prognostic value was performed in two independent datasets from Gene Expression Omnibus (GEO). Decision analysis curve was finally used to evaluate clinical usefulness of the BACH1-based model beyond pathologic stage alone. Results. BACH1 was an independent prognostic factor for lung adenocarcinoma. High-expression BACH1 cases had worse OS. BACH1-based prognostic model showed an ideal C-index and t -AUC and validated by two GEO datasets, independently. More importantly, the BACH1-based model indicated positive clinical applicability by DCA curves. Conclusion. Our research confirmed that BACH1 was an important predictor of prognosis in early-stage lung adenocarcinoma. The higher the expression of BACH1, the worse OS of the patients.


2021 ◽  
Vol 12 ◽  
Author(s):  
Siqin Zhang ◽  
Shaoyong Chen ◽  
Yuchen Wang ◽  
Yuxiang Zhan ◽  
Jiarui Li ◽  
...  

Although immunotherapy has recently demonstrated a substantial promise in treating advanced thyroid carcinoma (THCA), it is not appropriate for all THCA patients. As a result, this study aims to identify biomarkers for predicting immunotherapy efficacy and prognosis in THCA patients based on a constructed prognostic model. The transcriptomic and corresponding clinical data of THCA patients were obtained from the Cancer Genome Atlas (TCGA) database. We identified differentially expressed genes (DEGs) between THCA and normal samples and performed an intersection analysis of DEGs with immune-related genes (IRGs) downloaded from the ImmPort database. Functional enrichment analysis was performed on the chosen immune-related DEGs. Subsequently, Cox and LASSO regression analyses were conducted to obtain three hub immune-related DEGs, including PPBP, SEMA6B, and GCGR. Following that, a prognostic risk model was established and validated based on PPBP, SEMA6B, and GCGR genes to predict immunotherapy efficacy and THCA prognosis. Finally, we investigated the association between the constructed risk model and tumor mutational burden (TMB), abundance of tumor-infiltrating immune cells (TICs) as well as immunotherapeutic targets (PDL-1, PD-1, and CTLA4) in THCA. THCA patients in the high-risk score (RS) group showed higher TMB levels and worse prognosis than the low RS group. Patients in the high-RS group had higher proportions of monocytes, M2 macrophages, and activated dendritic cells, whereas those in the low-RS group exhibited higher numbers of M1 macrophages and dendritic resting cells. Our data implied that the constructed THCA prognostic model was sound and we concluded that the THCA patients having high TMB and low PD-L1 expression levels might respond poorly to immunotherapy. Taken together, we constructed a novel prognostic model for THCA patients to predict their prognosis and immunotherapy efficacy, providing a viable option for the future management of THCA patients in the clinic.


2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Jiahua Liu ◽  
Chunhui Jiang ◽  
Chunjie Xu ◽  
Dongyang Wang ◽  
Yuguang Shen ◽  
...  

AbstractThe overall survival of metastatic colon adenocarcinoma (COAD) remains poor, so it is important to explore the mechanisms of metastasis and invasion. This study aimed to identify invasion-related genetic markers for prognosis prediction in patients with COAD. Three molecular subtypes (C1, C2, and C3) were obtained based on 97 metastasis-related genes in 365 COAD samples from The Cancer Genome Atlas (TCGA). A total of 983 differentially expressed genes (DEGs) were identified among the different subtypes by using the limma package. A 6-gene signature (ITLN1, HOXD9, TSPAN11, GPRC5B, TIMP1, and CXCL13) was constructed via Lasso-Cox analysis. The signature showed strong robustness and could be used in the training, testing, and external validation (GSE17537) cohorts with stable predictive efficiency. Compared with other published signatures, our model showed better performance in predicting outcomes. Pan-cancer expression analysis results showed that ITLN1, TSPAN11, CXCL13, and GPRC5B were downregulated and TIMP1 was upregulated in most tumor samples, including COAD, which was consistent with the results of the TCGA and GEO cohorts. Western blot analysis and immunohistochemistry were performed to validate protein expression. Tumor immune infiltration analysis results showed that TSPAN11, GPRC5B, TIMP1, and CXCL13 protein levels were significantly positively correlated with CD4+ T cells, macrophages, neutrophils, and dendritic cells. Further, the TIMP1 and CXCL13 proteins were significantly related to the tumor immune infiltration of CD8+ T cells. We recommend using our signature as a molecular prognostic classifier to assess the prognostic risk of patients with COAD.


PeerJ ◽  
2019 ◽  
Vol 7 ◽  
pp. e8128 ◽  
Author(s):  
Cheng Yue ◽  
Hongtao Ma ◽  
Yubai Zhou

Background Lung cancer has the highest morbidity and mortality worldwide, and lung adenocarcinoma (LADC) is the most common pathological subtype. Accumulating evidence suggests the tumor microenvironment (TME) is correlated with the tumor progress and the patient’s outcome. As the major components of TME, the tumor-infiltrated immune cells and stromal cells have attracted more and more attention. In this study, differentially expressed immune and stromal signature genes were used to construct a TME-related prognostic model for predicting the outcomes of LADC patients. Methods The expression profiles of LADC samples with clinical information were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO). The differentially expressed genes (DEGs) related to the TME of LADC were identified using TCGA dataset by Wilcoxon rank sum test. The prognostic effects of TME-related DEGs were analyzed using univariate Cox regression. Then, the least absolute shrinkage and selection operator (LASSO) regression was performed to reduce the overfit and the number of genes for further analysis. Next, the prognostic model was constructed by step multivariate Cox regression and risk score of each sample was calculated. Then, survival and Receiver Operating Characteristic (ROC) analyses were conducted to validate the model using TCGA and GEO datasets, respectively. The Kyoto Encyclopedia of Genes and Genomes analysis of gene signature was performed using Gene Set Enrichment Analysis (GSEA). Finally, the overall immune status, tumor purity and the expression profiles of HLA genes of high- and low-risk samples was further analyzed to reveal the potential mechanisms of prognostic effects of the model. Results A total of 93 TME-related DEGs were identified, of which 23 DEGs were up-regulated and 70 DEGs were down-regulated. The univariate cox analysis indicated that 23 DEGs has the prognostic effects, the hazard ratio ranged from 0.65 to 1.25 (p < 0.05). Then, seven genes were screened out from the 23 DEGs by LASSO regression method and were further analyzed by step multivariate Cox regression. Finally, a three-gene (ADAM12, Bruton Tyrosine Kinase (BTK), ERG) signature was constructed, and ADAM12, BTK can be used as independent prognostic factors. The three-gene signature well stratified the LADC patients in both training (TCGA) and testing (GEO) datasets as high-risk and low-risk groups, the 3-year area under curve (AUC) of ROC curves of three GEO sets were 0.718 (GSE3141), 0.646 (GSE30219) and 0.643 (GSE50081). The GSEA analysis indicated that highly expressed ADAM12, BTK, ERG mainly correlated with the activation of pathways involving in focal adhesion, immune regulation. The immune analysis indicated that the low-risk group has more immune activities and higher expression of HLA genes than that of the high-risk group. In sum, we identified and constructed a three TME-related DEGs signature, which could be used to predict the prognosis of LADC patients.


2020 ◽  
Vol 2020 ◽  
pp. 1-12 ◽  
Author(s):  
Xiaolin Yu ◽  
Xiaomei Zhang ◽  
Yanxia Zhang

Lung adenocarcinoma (LUAD) is a common subtype of lung cancer with a depressing survival rate. The reprogramming of tumor metabolism was identified as a new hallmark of cancer in tumor microenvironment (TME), and we made a comprehensive exploration to reveal the prognostic role of the metabolic-related genes. Transcriptome profiling data of LUAD were, respectively, downloaded from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database. Based on the extracted metabolic-related genes, a novel 5-gene metabolic prognostic signature (including GNPNAT1, LPGAT1, TYMS, LDHA, and PTGES) was constructed by univariate Cox regression and least absolute shrinkage and selection operator (LASSO) regression. This signature confirmed its robustness and accuracy by external validation in multiple databases. It could be an independent risk factor for LUAD, and the nomograms possessed moderately accurate performance with the C-index of 0.755 (95% confidence interval: 0.706–0.804) and 0.691 (95% confidence interval: 0.636–0.746) in training set and testing set. This signature could reveal the metabolic features according to the results of gene set enrichment analysis (GSEA) and meanwhile monitor the status of TME through ESTIMATE scores and the infiltration levels of immune cells. In conclusion, this gene signature is a cost-effective tool which could indicate the status of TME to provide more clues in the exploration of new diagnostic and therapeutic strategy.


2020 ◽  
Vol 21 (23) ◽  
pp. 9169
Author(s):  
Mingjun Zheng ◽  
Heather Mullikin ◽  
Anna Hester ◽  
Bastian Czogalla ◽  
Helene Heidegger ◽  
...  

(1) Background: Biomarkers might play a significant role in predicting the clinical outcomes of patients with ovarian cancer. By analyzing lipid metabolism genes, future perspectives may be uncovered; (2) Methods: RNA-seq data for serous ovarian cancer were downloaded from The Cancer Genome Atlas and Gene Expression Omnibus databases. The non-negative matrix factorization package in programming language R was used to classify molecular subtypes of lipid metabolism genes and the limma package in R was performed for functional enrichment analysis. Through lasso regression, we constructed a multi-gene prognosis model; (3) Results: Two molecular subtypes were obtained and an 11-gene signature was constructed (PI3, RGS, ADORA3, CH25H, CCDC80, PTGER3, MATK, KLRB1, CCL19, CXCL9 and CXCL10). Our prognostic model shows a good independent prognostic ability in ovarian cancer. In a nomogram, the predictive efficiency was notably superior to that of traditional clinical features. Related to known models in ovarian cancer with a comparable amount of genes, ours has the highest concordance index; (4) Conclusions: We propose an 11-gene signature prognosis prediction model based on lipid metabolism genes in serous ovarian cancer.


2021 ◽  
Vol 12 ◽  
Author(s):  
Quanxiao Li ◽  
Limin Jin ◽  
Meng Jin

Hepatocellular carcinoma (HCC) is the most common form of liver cancer with limited therapeutic options and low survival rate. The hypoxic microenvironment plays a vital role in progression, metabolism, and prognosis of malignancies. Therefore, this study aims to develop and validate a hypoxia gene signature for risk stratification and prognosis prediction of HCC patients. The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) databases were used as a training cohort, and one Gene Expression Omnibus database (GSE14520) was served as an external validation cohort. Our results showed that eight hypoxia-related genes (HRGs) were identified by the least absolute shrinkage and selection operator analysis to develop the hypoxia gene signature and demarcated HCC patients into the high- and low-risk groups. In TCGA, ICGC, and GSE14520 datasets, patients in the high-risk group had worse overall survival outcomes than those in the low-risk group (all log-rank P &lt; 0.001). Besides, the risk score derived from the hypoxia gene signature could serve as an independent prognostic factor for HCC patients in the three independent datasets. Finally, a nomogram including the gene signature and tumor-node-metastasis stage was constructed to serve clinical practice. In the present study, a novel hypoxia signature risk model could reflect individual risk classification and provide therapeutic targets for patients with HCC. The prognostic nomogram may help predict individualized survival.


Author(s):  
Jindong Xie ◽  
Yutian Zou ◽  
Feng Ye ◽  
Wanzhen Zhao ◽  
Xinhua Xie ◽  
...  

Regarded as the most invasive subtype, triple-negative breast cancer (TNBC) lacks the expression of estrogen receptors (ERs), progesterone receptors (PRs), and human epidermal growth factor receptor 2 (HER2) proteins. Platelets have recently been shown to be associated with metastasis of malignant tumors. Nevertheless, the status of platelet-related genes in TNBC and their correlation with patient prognosis remain unknown. In this study, the expression and variation levels of platelet-related genes were identified and patients with TNBC were divided into three subtypes. We collected cohorts from The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) databases. By applying the least absolute shrinkage and selection operator (LASSO) Cox regression method, we constructed a seven-gene signature which classified the two cohorts of patients with TNBC into low- or high-risk groups. Patients in the high-risk group were more likely to have lower survival rates than those in the low-risk group. The risk score, incorporated with the clinical features, was confirmed as an independent factor for predicting the overall survival (OS) time. Functional enrichment analyses revealed the involvement of a variety of vital biological processes and classical cancer-related pathways that could be important to the ultimate prognosis of TNBC. We then built a nomogram that performed well. Moreover, we tested the model in other cohorts and obtained positive outcomes. In conclusion, platelet-related genes were closely related to TNBC, and this novel signature could serve as a tool for the assessment of clinical prognosis.


Sign in / Sign up

Export Citation Format

Share Document