Dataset Information

Identification of Iron Metabolism-Related Genes as Prognostic Indicators for Lower-Grade Glioma.

ABSTRACT: Lower-grade glioma (LGG) is characterized by genetic and transcriptional heterogeneity, and a dismal prognosis. Iron metabolism is considered central for glioma tumorigenesis, tumor progression and tumor microenvironment, although key iron metabolism-related genes are unclear. Here we developed and validated an iron metabolism-related gene signature LGG prognosis. RNA-sequence and clinicopathological data from The Cancer Genome Atlas (TCGA) and the Chinese Glioma Genome Atlas (CGGA) were downloaded. Prognostic iron metabolism-related genes were screened and used to construct a risk-score model via differential gene expression analysis, univariate Cox analysis, and the Least Absolute Shrinkage and Selection Operator (LASSO)-regression algorithm. All LGG patients were stratified into high- and low-risk groups, based on the risk score. The prognostic significance of the risk-score model in the TCGA and CGGA cohorts was evaluated with Kaplan-Meier (KM) survival and receiver operating characteristic (ROC) curve analysis. Risk- score distributions in subgroups were stratified by age, gender, the World Health Organization (WHO) grade, isocitrate dehydrogenase 1 (IDH1) mutation status, the O⁶-methylguanine-DNA methyl-transferase (MGMT) promoter-methylation status, and the 1p/19q co-deletion status. Furthermore, a nomogram model with a risk score was developed, and its predictive performance was validated with the TCGA and CGGA cohorts. Additionally, the gene set enrichment analysis (GSEA) identified signaling pathways and pathological processes enriched in the high-risk group. Finally, immune infiltration and immune checkpoint analysis were utilized to investigate the tumor microenvironment characteristics related to the risk score. We identified a prognostic 15-gene iron metabolism-related signature and constructed a risk-score model. High risk scores were associated with an age of > 40, wild-type IDH1, a WHO grade of III, an unmethylated MGMT promoter, and 1p/19q non-codeletion. ROC analysis indicated that the risk-score model accurately predicted 1-, 3-, and 5-year overall survival rates of LGG patients in the both TCGA and CGGA cohorts. KM analysis showed that the high-risk group had a much lower overall survival than the low-risk group (P < 0.0001). The nomogram model showed a strong ability to predict the overall survival of LGG patients in the TCGA and CGGA cohorts. GSEA analysis indicated that inflammatory responses, tumor-associated pathways, and pathological processes were enriched in high-risk group. Moreover, a high risk score correlated with the infiltration immune cells (dendritic cells, macrophages, CD4+ T cells, and B cells) and expression of immune checkpoint (PD1, PDL1, TIM3, and CD48). Our prognostic model was based on iron metabolism-related genes in LGG, can potentially aid in LGG prognosis, and provides potential targets against gliomas.

SUBMITTER: Xu S

PROVIDER: S-EPMC8458946 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundThe thyroid cancer subtype that occurs more frequently is papillary thyroid carcinoma (PTC). Despite a good surgical outcome, treatment with traditional antitumor therapy does not offer ideal results for patients with radioiodine resistance, recurrence, and metastasis. The evidence for the connection between iron metabolism imbalance and cancer development and oncogenesis is growing. Nevertheless, the iron metabolism impact on PTC prognosis is still indefinite.MethodsHerein, we acquired the medical data and gene expression of individuals with PTC from The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) database. Typically, three predictive iron metabolism-related genes (IMRGs) were examined and employed to build a risk score (RS) model via the least absolute shrinkage and selection operator (LASSO) regression, univariate Cox, and differential gene expression analyses. Then we analyzed somatic mutation and immune cell infiltration among RS groups. We also validated the prognostic value of two IMRGs (SFXN3 and TFR2) by verifying their biological function through in vitro experiments.ResultsBased on RS, all patients with PTC were stratified into low- and high-risk groups, where Kaplan-Meier analysis indicated that disease-free survival (DFS) in the high-risk group was much lower than in the low-risk group (P < 0.0001). According to ROC analysis, the RS model successfully predicted the 1-, 3-, and 5-year DFS of individuals with PTC. Additionally, in the TCGA cohort, a nomogram model with RS was developed and exhibited a strong capability to anticipate PTC patients' DFS. In the high-risk group, the enriched pathological processes and signaling mechanisms were detected utilizing the gene set enrichment analysis (GSEA). Moreover, the high-risk group had a significantly higher level of BRAF mutation frequency, tumor mutation burden, and immune cell infiltration than the low-risk group. In vitro experiments found that silencing SFXN3 or TFR2 significantly reduced cell viability.ConclusionCollectively, our predictive model depended on IMRGs in PTC, which could be potentially utilized to predict the PTC patients' prognosis, schedule follow-up plans, and provide potential targets against PTC.

Project description:Previous studies have shown that the prognosis of patients with lower-grade glioma (LGG) is closely related to the infiltration of immune cells and the expression of long non-coding RNAs (lncRNAs). In this paper, we applied single-sample gene set enrichment analysis (ssGSEA) algorithm to evaluate the expression level of immune genes from tumor tissues in The Cancer Genome Atlas (TCGA) database, and divided patients into the high immune group and the low immune group, which were separately analyzed for differential expression. Venn analysis was taken to select 36 immune-related lncRNAs. To construct a prognostic model of LGG based on immune-related lncRNAs, we divided patients into a training set and a verification set at a ratio of 2:1. Univariate Cox regression and the Least Absolute Shrinkage and Selection Operator (LASSO) regression were performed to select 11 immune-related lncRNAs associated with the prognosis of LGG, and based on these selected lncRNAs, the risk scoring model was constructed. Through Kaplan-Meier analysis, the overall survival (OS) of patients in the high-risk group was significantly lower than that of the low-risk group. Then, established a nomogram including age, gender, neoplasm histologic grade, and risk score. Meanwhile, the predictive performance of the model was evaluated by calculating the C-index, drawing the calibration chart, the clinical decision curve as well as the Receiver Operating Characteristic (ROC) curve. Similar results were obtained by utilizing the validation data to verify the above consequences. Based on the TIMER database, the correlation analysis showed that the 11 immune-related lncRNAs risk score of LGG were in connection with the infiltration of the subtypes of immune cells. Subsequently, we performed enrichment analysis, whose results showed that these immune-related lncRNAs played important roles in the progress of LGG. In conclusion, these 11 immune-related lncRNAs have the potential to predict the prognosis of patients with LGG, which may play a key role in the development of LGG.

Project description:Energy metabolic processes play important roles for tumor malignancy, indicating that related protein-coding genes and regulatory upstream genes (such as long noncoding RNAs (lncRNAs)) may represent potential biomarkers for prognostic prediction. This study will develop a new energy metabolism-related lncRNA-mRNA prognostic signature for lower-grade glioma (LGG) patients. A GSE4290 dataset obtained from Gene Expression Omnibus was used for screening the differentially expressed genes (DEGs) and lncRNAs (DELs). The Cancer Genome Atlas (TCGA) dataset was used as the prognosis training set, while the Chinese Glioma Genome Atlas (CGGA) was for the validation set. Energy metabolism-related genes were collected from the Molecular Signatures Database (MsigDB), and a coexpression network was established between energy metabolism-related DEGs and DELs to identify energy metabolism-related DELs. Least absolute shrinkage and selection operator (LASSO) analysis was performed to filter the prognostic signature which underwent survival analysis and nomogram construction. A total of 1613 DEGs and 37 DELs were identified between LGG and normal brain tissues. One hundred and ten DEGs were overlapped with energy metabolism-related genes. Twenty-seven DELs could coexpress with 67 metabolism-related DEGs. LASSO regression analysis showed that 9 genes in the coexpression network were the optimal signature and used to construct the risk score. Kaplan-Meier curve analysis showed that patients with a high risk score had significantly worse OS than those with a low risk score (TCGA: HR = 3.192, 95%CI = 2.182-4.670; CGGA: HR = 1.922, 95%CI = 1.431-2.583). The predictive accuracy of the risk score was also high according to the AUC of the ROC curve (TCGA: 0.827; CGGA: 0.806). Multivariate Cox regression analyses revealed age, IDH1 mutation, and risk score as independent prognostic factors, and thus, a prognostic nomogram was established based on these three variables. The excellent prognostic performance of the nomogram was confirmed by calibration and discrimination analyses. In conclusion, our findings provided a new biomarker for the stratification of LGG patients with poor prognosis.

Project description:BackgroundTo identify prognostic hub genes which associated with tumor microenvironment (TME) in lower grade glioma (LGG) of central nervous system.MethodsWe downloaded LGG patients gene transcriptome profiles of the central nervous system in The Cancer Genome Atlas (TCGA) database. Clinical characteristics and survival data through the Genomic Data Commons (GDC) tool were extracted. We used limma package for normalization processing. Scores of immune, stromal and ESTIMATE were calculated using ESTIMATE algorithm. Then, box plots were applied to explore the association between immune scores, stromal scores, ESTIMATE scores and histological type, tumor grade. Kaplan-Meier (K-M) analysis was utilized to explore the prognostic value of scores. Furthermore, heatmaps and volcano plots were applied for visualizing expression of differential expressed-gene screening and cluster analysis. Venn plots were constructed to screen the intersected differentially expressed genes (DEGs). In addition, enrichment of functions and signaling pathways and Gene Set Enrichment Analysis (GESA) of the DEGs were performed. Then we used protein-protein interaction (PPI) network and Cytoscape software to identify hub genes. We evaluated the prognostic value of hub genes and risk score (RS) calculated based on multivariate cox regression analysis. Finally, relationships of hub genes with the TME of LGG patients were evaluated based on tumor immune estimation resource (TIMER) database.ResultsGene expression profiles and clinical data of 514 LGG samples were extracted and the results revealed that higher scores were significantly related with histological types and higher tumor grade (P<0.0001, respectively). Besides, higher scores were associated with worse survival outcomes in immune scores (P=0.0167), stromal scores (P=0.0035) and ESTIMATE scores (P=0.0190). Then, 785 up-regulated intersected genes and 357 down-regulated intersected genes were revealed. Functional enrichment analysis revealed that intersected genes were associated with immune response, inflammatory response, plasma membrane and receptor activity. After PPI network construction and cytoHubba analysis, 25 tumor immune-related hub genes were identified and enriched pathways were identified by GSEA. Besides, receiver operating characteristic (ROC) curves showed significantly predictive accuracy [area under curve (AUC) =0.771] of RS. Furthermore, significant prognostic values of hub genes were observed, and the relationships between hub genes and LGG TME were demonstrated.ConclusionsWe identified 25 TME-related genes which significantly associated with overall survival in patients with central nervous system LGG from TCGA database.

Dataset Information

Identification of Iron Metabolism-Related Genes as Prognostic Indicators for Lower-Grade Glioma.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets