Construction of a ferroptosis-related five-lncRNA signature for predicting prognosis and immune response in thyroid carcinoma.

Thyroid carcinoma (THCA) is the most common endocrine-related malignant tumor. Despite the good prognosis, some THCA patients may deteriorate into more aggressive diseases, leading to poor survival. This may be alleviated by developing a novel model to predict the risk of THCA, including recurrence and survival. Ferroptosis is an iron-dependent, oxidative, non-apoptotic form of cell death initially described in mammalian cells, and plays an important role in various cancers. To explore the potential prognostic value of ferroptosis in THCA, ferroptosis-related long non-coding RNAs (FRLs) were used to construct model for risk prediction of THCA. RNA-sequencing data of THCA patients and ferroptosis-related genes were downloaded from The Cancer Genome Atlas (TCGA) and FerrDb, respectively. A total of 502 patients with complete data were randomly separated into a training cohort and a validation cohort at the ratio of 2:1. The Pearson correlation coefficients were calculated to determine the correlation between ferroptosis-related genes (FRGs) and the corresponding lncRNAs, and those meeting the screening conditions were defined as FRLs. Gene Expression Omnibus (GEO) database and qRT-PCR were used to verify the expression level of FRLs in THCA tissues. Univariate and multivariate cox regression analysis were performed to construct a FRLs signature based on lowest Akaike information criterion (AIC) value in the training cohort, then further tested in the validation cohort and the entire cohort. Gene set enrichment analysis (GSEA) and functional enrichment analysis were used to analyze the biological functions and signal pathways related to differentially expressed genes between the high-risk and low-risk groups. Finally, the relative abundance of different tumor-infiltrating immune cells were calculated by CIBERSORT algorithm. The patients were divided into high-risk group and low-risk group based on a 5-FRLs signature (AC055720.2, DPP4-DT, AC012038.2, LINC02454 and LINC00900) in training cohort, validation cohort and entire cohort. Through Kaplan-Meier analysis and area under ROC curve (AUC) value, patients in the high-risk group exhibited worse prognosis than patients in the low-risk group. GEO database and qRT-PCR confirmed that LINC02454 and LINC00900 were up-regulated in THCA. Univariate and multivariate cox regression analyses showed that the risk score was an independent prognostic indicator. GSEA and functional enrichment analysis confirmed that immune-related pathways against cancer were significantly activated in the low-risk THCA patients. Further analysis showed that the immune cells such as plasma cells, T cells CD8 and macrophages M1, and the expression of immune checkpoint molecules, including PD-1, PD-L1, CTLA4, and LAG3, were remarkably higher in the low-risk group. Our study used the TCGA THCA dataset to construct a novel FRLs prognostic model which could precisely predict the prognosis of THCA patients. These FRLs potentially mediate anti-tumor immunity and serve as therapeutic targets for THCA, which provided the novel insight into treatment of THCA.

Qin Y ，Zhang D ，Zhang H ，Hou L ，Wang Z ，Yang L ，Zhang M ，Zhao G ，Yao Q ，Ling R ，Zhang J ... -  《Cancer Cell International》

Development and Validation of a Novel Ferroptosis-Related LncRNA Signature for Predicting Prognosis and the Immune Landscape Features in Uveal Melanoma.

Ferroptosis is a newly iron-dependent mode of programmed cell death that is involved in a variety of malignancies. But no research has shown a link between ferroptosis-related long non-coding RNAs (FRLs) and uveal melanoma (UM). We aimed to develop a predictive model for UM and explore its potential function in relation to immune cell infiltration. Identification of FRLs was performed using the Cancer Genome Atlas (TCGA) and FerrDb databases. To develop a prognostic FRLs signature, univariate Cox regression and least absolute shrinkage and selection operator (LASSO) were used in training cohort. Kaplan-Meier (K-M) and receiver operating characteristic (ROC) curve analyses were used to assess the reliability of the risk model. The immunological functions of FRLs signature were determined using gene set enrichment analysis (GSEA). Immunological cell infiltration and immune treatment were studied using the ESTIMATE, CIBERSORT, and ssGSEA algorithms. Finally, in vitro assays were carried out to confirm the biological roles of FRLs with known primer sequences (LINC00963, PPP1R14B.AS1, and ZNF667.AS1). A five-genes novel FRLs signature was identified. The mean risk score generated by this signature was used to create two risk groups. The high-risk score UM patients had a lower overall survival rate. The area under the curve (AUC) of ROC and K-M analysis further validated the strong prediction capacity of the prognostic signature. Immune cells such as memory CD8 T cells, M1 macrophages, monocytes, and B cells showed a substantial difference between the two groups. GSEA enrichment results showed that the FRLs signature was linked to certain immune pathways. Moreover, UM patients with high-risk scores were highly susceptible to several chemotherapy drugs, such as cisplatin, imatinib, bortezomib, and pazopanib. Finally, the experimental validation confirmed that knockdown of three identified lncRNA (LINC00963, PPP1R14B.AS1, and ZNF667.AS1) suppressed the invasive ability of tumor cells in vitro. The five-FRLs (AC104129.1, AC136475.3, LINC00963, PPP1R14B.AS1, and ZNF667.AS1) signature has effects on clinical survival prediction and selection of immunotherapies for UM patients.

Ma X ，Yu S ，Zhao B ，Bai W ，Cui Y ，Ni J ，Lyu Q ，Zhao J ... -  《Frontiers in Immunology》

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

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. 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. 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. 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.

Lin R ，Fogarty CE ，Ma B ，Li H ，Ni G ，Liu X ，Yuan J ，Wang T ... -  《BMC GENOMICS》

Analysis of Ferroptosis-Related LncRNAs Signatures Associated with Tumor Immune Infiltration and Experimental Validation in Clear Cell Renal Cell Carcinoma.

Clear cell renal cell carcinoma (ccRCC) is the most aggressive subtype of renal cell carcinoma. Ferroptosis is an iron-dependent programmed cell death. Long non-coding RNAs (lncRNAs) emerge as a critical role in regulating cancer progression. This study aimed to identify molecular regulation of ferroptosis-related lncRNAs (FRLs) in ccRCC. The prognostic value of FRLs was investigated in ccRCC samples downloaded from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) dataset. The FRLs were screened out by Pearson correlation test. The 465 FRLs confirmed as potential prognostic factors through univariate Cox regression analysis were entered into Lasso and multivariate Cox regression to build a FRLs prognostic signature. A risk score based on the prognostic model divided ccRCC patients into low- and high-risk groups. A prognostic nomogram, derived from the prognostic signature and integrating clinical characteristics, was constructed. Gene set enrichment analysis (GSEA) revealed the immune- and tumor-associated pathways. Two distinct clusters were identified with different immune signatures through consensus clustering analysis. The prognostic value of some hub FRLs was externally validated via three GEO datasets (GSE46699, GSE53757 and GSE66272) and online databases. Finally, the three FRLs (LINC00460, LINC00941 and LINC02027) were verified through in vitro experiments. The FRLs prognostic signature, including 7 independent prognostic lncRNAs, exhibited good accuracy in predicting overall survival (OS) of ccRCC patients. This signature was correlated with immune infiltration and immune checkpoint blockade (ICB). We correlated two distinct clusters with immune infiltration signature of ccRCC. The worse prognosis of cluster 2 was probably mediated by immune evasion. We also found that the expression levels of LINC00460 and LINC00941 in ccRCC cell lines were higher than those in HK-2 cells, but LINC02027 showed the inverse trend. Collectively, our study demonstrated a FRLs prognostic signature which had great clinical value in prognosis assessment.

Chen X ，Tu J ，Ma L ，Huang Y ，Yang C ，Yuan X ... -  《International Journal of General Medicine》

Construction of a ferroptosis-related signature based on seven lncRNAs for prognosis and immune landscape in clear cell renal cell carcinoma.

Recent studies have demonstrated that long non-coding RNAs (lncRNAs) are involved in regulating tumor cell ferroptosis. However, prognostic signatures based on ferroptosis-related lncRNAs (FRLs) and their relationship to the immune microenvironment have not been comprehensively explored in clear cell renal cell carcinoma (ccRCC). In the present study, the expression profiles of ccRCC were acquired from The Cancer Genome Atlas (TCGA) database; 459 patient specimens and 69 adjacent normal tissues were randomly separated into training or validation cohorts at a 7:3 ratio. We identified 7 FRLs that constitute a prognostic signature according to the differential analysis, correlation analysis, univariate regression, and least absolute shrinkage and selection operator (LASSO) Cox analysis. To identify the independence of risk score as a prognostic factor, univariate and multivariate regression analyses were also performed. Furthermore, CIBERSORT was conducted to analyze the immune infiltration of patients in the high-risk and low-risk groups. Subsequently, the differential expression of immune checkpoint and m6A genes was analyzed in the two risk groups. A 7-FRLs prognostic signature of ccRCC was developed to distinguish patients into high-risk and low-risk groups with significant survival differences. This signature has great prognostic performance, with the area under the curve (AUC) for 1, 3, and 5 years of 0.713, 0.700, 0.726 in the training set and 0.727, 0.667, and 0.736 in the testing set, respectively. Moreover, this signature was significantly associated with immune infiltration. Correlation analysis showed that risk score was positively correlated with regulatory T cells (Tregs), activated CD4 memory T cells, CD8 T cells and follicular helper T cells, whereas it was inversely correlated with monocytes and M2 macrophages. In addition, the expression of fourteen immune checkpoint genes and nine m6A-related genes varied significantly between the two risk groups. We established a novel FRLs-based prognostic signature for patients with ccRCC, containing seven lncRNAs with precise predictive performance. The FRLs prognostic signature may play a significant role in antitumor immunity and provide a promising idea for individualized targeted therapy for patients with ccRCC.

Wei SY ，Feng B ，Bi M ，Guo HY ，Ning SW ，Cui R ... -  《BMC Medical Genomics》