Dataset Information

Identification of ferroptosis genes in immune infiltration and prognosis in thyroid papillary carcinoma using network analysis.

ABSTRACT:

Background

Papillary thyroid carcinoma (PTC) is the most common thyroid cancer. While many patients survive, a portion of PTC cases display high aggressiveness and even develop into refractory differentiated thyroid carcinoma. This may be alleviated by developing a novel model to predict the risk of recurrence. Ferroptosis is an iron-dependent form of regulated cell death (RCD) driven by lethal accumulation of lipid peroxides, is regulated by a set of genes and shows a variety of metabolic changes. To elucidate whether ferroptosis occurs in PTC, we analyse the gene expression profiles of the disease and established a new model for the correlation.

Methods

The thyroid carcinoma (THCA) datasets were downloaded from The Cancer Genome Atlas (TCGA), UCSC Xena and MisgDB, and included 502 tumour samples and 56 normal samples. A total of 60 ferroptosis related genes were summarised from MisgDB database. Gene set enrichment analysis (GSEA) and Gene set variation analysis (GSVA) were used to analyse pathways potentially involving PTC subtypes. Single sample GSEA (ssGSEA) algorithm was used to analyse the proportion of 28 types of immune cells in the tumour immune infiltration microenvironment in THCA and the hclust algorithm was used to conduct immune typing according to the proportion of immune cells. Spearman correlation analysis was performed on the ferroptosis gene expression and the correlation between immune infiltrating cells proportion. We established the WGCNA to identify genes modules that are highly correlated with the microenvironment of immune invasion. DEseq2 algorithm was further used for differential analysis of sequencing data to analyse the functions and pathways potentially involving hub genes. GO and KEGG enrichment analysis was performed using Clusterprofiler to explore the clinical efficacy of hub genes. Univariate Cox analysis was performed for hub genes combined with clinical prognostic data, and the results was included for lasso regression and constructed the risk regression model. ROC curve and survival curve were used for evaluating the model. Univariate Cox analysis and multivariate Cox analysis were performed in combination with the clinical data of THCA and the risk score value, the clinical efficacy of the model was further evaluated.

Results

We identify two subtypes in PTC based on the expression of ferroptosis related genes, with the proportion of cluster 1 significantly higher than cluster 2 in ferroptosis signature genes that are positively associated. The mutations of Braf and Nras are detected as the major mutations of cluster 1 and 2, respectively. Subsequent analyses of TME immune cells infiltration indicated cluster 1 is remarkably richer than cluster 2. The risk score of THCA is in good performance evaluated by ROC curve and survival curve, in conjunction with univariate Cox analysis and multivariate Cox analysis results based on the clinical data shows that the risk score of the proposed model could be used as an independent prognostic indicator to predict the prognosis of patients with papillary thyroid cancer.

Conclusions

Our study finds seven crucial genes, including Ac008063.2, Apoe, Bcl3, Acap3, Alox5ap, Atxn2l and B2m, and regulation of apoptosis by parathyroid hormone-related proteins significantly associated with ferroptosis and immune cells in PTC, and we construct the risk score model which can be used as an independent prognostic index to predict the prognosis of patients with PTC.

SUBMITTER: Lin R

PROVIDER: S-EPMC8314640 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BACKGROUND:Papillary thyroid carcinoma (PTC) is the most common type of thyroid cancer (TC), accounting for more than 80% of all cases. Ferroptosis is a novel iron-dependent and Reactive oxygen species (ROS) reliant type of cell death which is distinct from the apoptosis, necroptosis and pyroptosis. Considerable studies have demonstrated that ferroptosis is involved in the biological process of various cancers. However, the role of ferroptosis in PTC remains unclear. This study aims at exploring the expression of ferroptosis-related genes (FRG) and their prognostic values in PTC. METHODS:A ferroptosis-related gene signature was constructed using lasso regression analysis through the PTC datasets of the Cancer Genome Atlas (TCGA). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed to investigate the bioinformatics functions of significantly different genes (SDG) of ferroptosis. Additionally, the correlations of ferroptosis and immune cells were assessed through the single-sample gene set enrichment analysis (ssGSEA) and CIBERSORT database. Finally, SDG were test in clinical PTC specimens and normal thyroid tissues. RESULTS:LASSO regression model was utilized to establish a novel FRG signature with 10 genes (ANGPTL7, CDKN2A, DPP4, DRD4, ISCU, PGD, SRXN1, TF, TFRC, TXNRD1) to predicts the prognosis of PTC, and the patients were separated into high-risk and low-risk groups by the risk score. The high-risk group had poorer survival than the low-risk group (p < 0.001). Receiver operating characteristic (ROC) curve analysis confirmed the signature's predictive capacity. Multivariate regression analysis identified the prognostic signature-based risk score was an independent prognostic indicator for PTC. The functional roles of the DEGs in the TGCA PTC cohort were explored using GO enrichment and KEGG pathway analyses. Immune related analysis demonstrated that the most types of immune cells and immunological function in the high-risk group were significant different with those in the low-risk group. Quantitative Real-Time Polymerase Chain Reaction (qRT-PCR) verified the SDG have differences in expression between tumor tissue and normal thyroid tissue. In addition, cell experiments were conducted to observe the changes in cell morphology and expression of signature's genes with the influence of ferroptosis induced by sorafenib. CONCLUSIONS:We identified differently expressed FRG that may involve in PTC. A ferroptosis-related gene signature has significant values in predicting the patients' prognoses and targeting ferroptosis may be an alternative for PTC's therapy.

Project description:Thyroid associated ophthalmopathy (TAO) is an orbital autoimmune inflammatory disease that is commonly associated with thyroid dysfunction. Although the etiology of TAO is unclear, ROS accumulation and oxidative stress have been closely linked to the pathogenesis of TAO. Ferroptosis is an iron-dependent programmed cell death characterized by intracellular labile iron levels, excessive accumulation of reactive oxygen species (ROS) and lipid peroxidation. Currently, there are few reports regarding the role of ferroptosis in TAO. This article aimed to identify ferroptosis-related genes (FRGs) with diagnostic and therapeutic potential in TAO and explore their relationship with immune cells and lncRNAs. GSE58331 was downloaded from Gene Expression Omnibus (GEO) database. A total of 162 DEGs were identified between 27 TAO samples and 22 health samples from GSE58331, among which six FRGs (CYBB, CTSB, SLC38A1, TLR4, PEX3, and ABCC1) were obtained. The AUC of SLC38A1, TLR4, PEX3 in lacrimal gland tissues was greater than 80 which suggested high diagnostic value in TAO. The result of immune cell infiltrate analysis indicated increased infiltration of monocytes (p < 0.001), macrophages M0(p = 0.039), mast cells activated (p = 0.008), and neutrophils (p = 0.045) in orbital tissues from TAO patients. Meanwhile, mast cells resting (p = 0.043) and macrophages M2 (p = 0.02) showed reduced infiltration in TAO samples. There were no gender differences in immune cell infiltration in the TAO patients. Two differentially expressed lncRNAs, LINC01140 and ZFHX4-AS1, in TAO groups were identified as ferroptosis-related lncRNAs. CYBB-LINC01140-TLR4, CYBB- LINC01140- SLC38A1, TLR4- LINC01140- SLC38A1, and CTSB- ZFHX4-AS1- CYBB may be potential RNA regulatory pathways in TAO. Targeted drugs and transcription factors for differential expressed FRGs were also screened out in our study. In vitro, experiments revealed that CTSB, PEX3, ABCC1 and ZFHX4-AS1(lncRNA) were differentially expressed in orbital fibroblasts (OFs) between TAO groups and healthy controls at the transcriptional level.

Project description:Adrenocortical carcinoma (ACC) is a rare malignant tumor with poor prognosis. Ferroptosis, a new form of cell death, differs from other forms of cell death and plays a vital role in tumor progress. Our study aimed to establish a ferroptosis-related signature with prognostic value in ACC. RNA-seq data and corresponding clinical characteristics for ACC were downloaded from TCGA and GEO databases. Genes included in ferroptosis risk signature were assessed by univariable and multivariable Cox regression analysis as well as lasso regression analysis. The prognostic value of the ferroptosis risk signature was assessed using K-M and ROC curves. Furthermore, we performed GSEA to discover the enriched gene sets in high-risk group. Additionally, TIMER website was applied to detect a possible connection between the signature and immune cells infiltration. ssGSEA was performed to evaluate scores of immune cells and immune-related pathways in two groups. A ferroptosis signature comprised of six genes (SLC7A11, TP53, HELLS, ACSL4, PCBP2, and HMGB1) was constructed to predict prognosis and reflect the immune infiltration in ACC. Patients in high-risk group were inclined to have worse prognosis. The ferroptosis model performed well in predicting prognosis and could be served as an independent indicator in ACC. GSEA revealed that gene sets correlated with biological processes including cell cycle, DNA replication, base excision repair, and P53 signaling pathway were highly enriched in high-risk group. In addition, we discovered that the expressional levels of hub genes were linked to six immune cells' infiltration in ACC tumor. ssGSEA revealed that contents of most immune cells significantly decreased in the high-risk group. In conclusion, the novel ferroptosis risk signature could be useful in predicting prognosis and reflecting immune infiltration in ACC. It also brings us new insights into the possible value of targeting ferroptosis during the therapy of ACC.

Project description:ObjectiveThe primary objective of this study is to conduct a comprehensive screening and analysis of differentially expressed genes related to disulfidoptosis (DEDRGs) in thyroid carcinoma (THCA). This entails delving into the intricate characterization of immune cell infiltration within the THCA context and subsequently formulating and validating a novel prognostic model.MethodTo achieve our objectives, we first delineated two distinct subtypes of disulfidoptosis-related genes (DRGs) via consensus clustering methodology. Subsequently, employing the limma R package, we identified the DEDRGs critical for our investigation. These DEDRGs underwent meticulous validation across various databases, alongside an in-depth analysis of gene regulation. Employing functional enrichment techniques, we explored the potential molecular mechanisms underlying disulfidoptosis in THCA. Furthermore, we scrutinized the immune landscape within the two identified subtypes utilizing CIBERSORT and ESTIMATE algorithms. The construction of the prognostic model for THCA entailed intricate methodologies including univariate, multivariate Cox regression, and LASSO regression algorithms. The validity and efficacy of our prognostic model were corroborated through Kaplan-Meier survival curves and ROC curves. Additionally, a nomogram was meticulously formulated to facilitate the prediction of patient prognosis. To fortify our findings, we conducted a comprehensive Bayesian co-localization analysis coupled with rigorous in vitro experimentation, aimed at unequivocally establishing the validity of the identified DEDRGs.ResultOur analyses unveiled Cluster C1, characterized by elevated expression levels of DEDRGs, as harboring a favorable prognosis accompanied by abundant immune cell infiltration. Correlation analyses underscored predominantly positive associations among the DEDRGs, further affirming their significance in THCA. Differential expression patterns of DEDRGs between tumor samples and normal tissues were evident across the GEPIA and HPA databases. Insights from the TIMER database underscored a robust correlation between DEDRGs and immune cell infiltration. KEGG analysis elucidated the enrichment of DEDRGs primarily in pivotal pathways including MAPK, PPAR signaling pathway, and Proteoglycans in cancer. Furthermore, analyses using CIBERSORT and ESTIMATE algorithms shed light on the crucial role played by DEDRGs in shaping the immune microenvironment. The prognostic model, anchored by five genes intricately associated with THCA prognosis, exhibited commendable predictive accuracy and was intricately linked to the tumor immune microenvironment. Notably, patients categorized with low-risk scores stood to potentially benefit more from immunotherapy. The validation of DEDRGs unequivocally underscores the protective role of INF2 in THCA.ConclusionIn summary, our study delineates two discernible subtypes intricately associated with DRGs, revealing profound disparities in immune infiltration and survival prognosis within the THCA milieu. The implications of our findings extend to potential treatment strategies for THCA patients, which could entail targeted interventions directed towards DEDRGs and prognostic genes, thereby influencing disulfidptosis and the immune microenvironment. Moreover, the robust predictive capability demonstrated by our prognostic model, based on the five genes (ANGPTL7, FIRRE, ODAPH, PROKR1, SFRP5), underscores its potential clinical utility in guiding personalized therapeutic approaches for THCA patients.

Project description:Hepatocellular carcinoma (HCC) stands as the prevailing manifestation of primary liver cancer and continues to pose a formidable challenge to human well-being and longevity, owing to its elevated incidence and mortality rates. Nevertheless, the quest for reliable predictive biomarkers for HCC remains ongoing. Recent research has demonstrated a close correlation between ferroptosis and disulfidptosis, two cellular processes, and cancer prognosis, suggesting their potential as predictive factors for HCC. In this study, we employed a combination of bioinformatics algorithms and machine learning techniques, leveraging RNA sequencing data, mutation profiles, and clinical data from HCC samples in The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO), and the International Cancer Genome Consortium (ICGC) databases, to develop a risk prognosis model based on genes associated with ferroptosis and disulfidptosis. We conducted an unsupervised clustering analysis, calculating a risk score (RS) to predict the prognosis of HCC using these genes. Clustering analysis revealed two distinct HCC clusters, each characterized by significantly different prognostic and immune features. The median RS stratified HCC samples in the TCGA, GEO, and ICGC cohorts into high-and low-risk groups. Importantly, RS emerged as an independent prognostic factor in all three cohorts, with the high-risk group demonstrating poorer prognosis and a more active immunosuppressive microenvironment. Additionally, the high-risk group exhibited higher expression levels of tumor mutation burden (TMB), immune checkpoints (ICs), and human leukocyte antigen (HLA), suggesting a heightened responsiveness to immunotherapy. A cancer stem cell infiltration analysis revealed a higher similarity between tumor cells and stem cells in the high-risk group. Furthermore, drug sensitivity analysis highlighted significant differences in response to antitumor drugs between the two risk groups. In summary, our risk prognostic model, constructed based on ferroptosis-related genes associated with disulfidptosis, effectively predicts HCC prognosis. These findings hold potential implications for patient stratification and clinical decision-making, offering valuable theoretical insights in this field.

Project description:Background: One particular type of cellular death that is known as ferroptosis is caused by the excessive lipid peroxidation. It is a regulated form of cell death that can affect the response of the tumor cells. Currently, it is not known if the presence of this condition can affect the prognosis of patients with multiple myeloma (MM). Methods: In this study, we studied the expression differences and prognostic value of ferroptosis-related genes (FRGs) in MM, and established a ferroptosis risk scoring model. In order to improve the prediction accuracy and clinical applicability, a nomogram was also established. Through gene enrichment analysis, pathways closely related to high-risk groups were identified. We then explored the differences in risk stratification in drug sensitivity and immune patterns, and evaluated their value in prognostic prediction and treatment response. Lastly, we gathered MM cell lines and samples from patients to confirm the expression of marker FRGs using quantitative real-time PCR (qRT-PCR). Results: The ability to predict the survival of MM patients is a challenging issue. Through the use of a risk model derived from ferroptosis, we were able to develop a more accurate prediction of the disease's prognosis. They were then validated by a statistical analysis, which showed that the model is an independent factor in the prognosis of MM. Patients of high ferroptosis risk scores had a much worse chance of survival than those in the low-risk groups. The calibration and power of the nomogram were also strong. We noted that the link between the ferroptosis risk score and the clinical treatment was suggested by the FRG's significant correlation with the immune checkpoint genes and the medication sensitivity. We validated the predictive model using qRT-PCR. Conclusion: We demonstrated the association between FRGs and MM, and developed a new risk model for prognosis in MM patients. Our study sheds light on the potential clinical relevance of ferroptosis in MM and highlights its potential as a therapeutic target for patients with this disease.