IRF2-ferroptosis related gene is associated with prognosis and EMT in gliomas

Ferroptosis is a new type of programmed cell death that has excellent anti-tumor potential in different tumors. However, the research on ferroptosis in glioma is still incomplete. In this study, we aimed to revealed the relationship between ferroptosis-related genes (FRGs) and glioma. We collected gene expression profiles of glioma patients from the TCGA and CGGA databases. All glioma samples were classified into five subtypes using the R software ConsensusClusterPlus. Subsequently, we performed single sample gene set enrichment analysis (ssGSEA) to explore the correlation between different subtypes and immune status and ferroptosis. Then co-expression modules were constructed via weighted gene co-expression network analysis (WGCNA). A Gene Ontology (GO) analysis was conducted to analyze the potential biological functions of the genes in the modules. Finally, we identified 10 hub genes using the PPI network. The in vitro experiments were used to validate our predictions. We found that the expression level of IRF2 is positively correlated with the grade of glioma. The overexpression of IRF2 could protect glioma cells from ferroptosis and enhance the invasive and migratory abilities. Silence of IRF2 had the opposite effect. In conclusion, we demonstrated a novel ferroptosis-related signature for predicting prognosis, and IRF2 could be a potential biomarker for diagnosis and treatment in glioma.     

Ferroptosis-related gene signature correlates with the tumor immune features and predicts the prognosis of glioma patients

Background: Glioma is a malignant intracranial tumor and the most fatal cancer. The role of ferroptosis in the clinical progression of gliomas is unclear.Method: Univariate and least absolute shrinkage and selection operator (Lasso) Cox regression methods were used to develop a ferroptosis-related signature (FRSig) using a cohort of glioma patients from the Chinese Glioma Genome Atlas (CGGA), and was validated using an independent cohort of glioma patients from The Cancer Genome Atlas (TCGA). A single-sample gene set enrichment analysis (ssGSEA) was used to calculate levels of the immune infiltration. Multivariate Cox regression was used to determine the independent prognostic role of clinicopathological factors and to establish a nomogram model for clinical application.Results: We analyzed the correlations between the clinicopathological features and ferroptosis-related gene (FRG) expression and established an FRSig to calculate the risk score for individual glioma patients. Patients were stratified into two subgroups with distinct clinical outcomes. Immune cell infiltration in the glioma microenvironment and immune-related indexes were identified that significantly correlated with the FRSig, the tumor mutation burden (TMB), copy number alteration (CNA), and immune checkpoint expression was also significantly positively correlated with the FRSig score. Ultimately, an FRSig-based nomogram model was constructed using the independent prognostic factors age, World Health Organization (WHO) grade, and FRSig score.Conclusion: We established the FRSig to assess the prognosis of glioma patients. The FRSig also represented the glioma microenvironment status. Our FRSig will contribute to improve patient management and individualized therapy by offering a molecular biomarker signature for precise treatment.     

Multi-omics Analysis of Ferroptosis Regulation Patterns and Characterization of Tumor Microenvironment in Patients with Oral Squamous Cell Carcinoma

Ferroptosis is a newly recognized mechanism of regulated cell death. It was reported to be highly associated with immune therapy and chemotherapy. However, its mechanism of regulation in the tumor microenvironment (TME) and influence on oral squamous cell carcinoma (OSCC) therapy are unknown. We identified a ferroptosis-specific gene-expression signature, an FPscore, developed by a principal component analysis (PCA) algorithm to evaluate the ferroptosis regulation patterns of individual tumor. Multi-omics analysis of ferroptosis regulation patterns was conducted. Three distinct ferroptosis regulation subtypes, which linked to outcomes and the clinical relevance of each patient, were established. A high FPscore of patients with OSCC was associated with a favorable prognosis, a ferroptosis-related immune-activation phenotype, potential sensitivities to the chemotherapy and immunotherapy. Importantly, a high FPscore correlated with a low gene copy number burden and high immune checkpoint expressions. We validated the prognostic value of the FPscore using independent immunotherapy and pan-cancer cohorts. Comprehensive evaluation of individual tumors with distinct ferroptosis regulation patterns provides new mechanistic insights, which may be clinically relevant for the application of combination therapies in OSCC.     

An inflammation-associated ferroptosis signature optimizes the diagnosis,prognosis evaluation and immunotherapy options in hepatocellular carcinoma

Inflammation and ferroptosis crosstalk complexly with immune microenvironment of hepatocellular carcinoma (HCC), thus affecting the efficacy of immunotherapy. Herein, our aim was to identify the inflammation-associated ferroptosis (IAF) biomarkers for contributing HCC. A total of 224 intersecting DEGs identified from different inflammation- and ferroptosis-subtypes were set as IAF genes. Seven of them including ADH4, APOA5, CFHR3, CXCL8, FTCD, G6PD and PON1 were used for construction of a risk model which classified HCC patients into two groups (high and low risk). HCC patients in the high-risk group exhibited shorter survival rate and higher immune score, and were predicted to have higher respond rate in immune checkpoint inhibition (ICI) therapy. Levels of the seven genes were significantly changed in HCC tissues in comparison to adjacent tissues. After inserting the gene expression into the risk model, we found that the risk model exhibited the higher diagnostic value for distinguish HCC tissues compared each single gene. Furthermore, HCC tissues from our research group with high-risk score exhibited more cases of microsatellite instability (MSI), heavier tumour mutational burden (TMB), higher expression level of PDL1 and cells with CD8. Knockdown of APOA5 reduced HCC cell proliferation combining with elevating inflammation and ferroptosis levels. In conclusion, we considered APOA5 maybe a novel target for suppressing HCC via simultaneously elevating inflammation and ferroptosis levels, and signature constructed by seven IAF genes including ADH4, APOA5, CFHR3, CXCL8, FTCD, G6PD and PON1 can act as a biomarker for optimising the diagnosis, prognosis evaluation and immunotherapy options in HCC patients.     

急性髓系白血病铁死亡相关基因筛选及预后模型的建立

基于急性髓系白血病(Acute Myeloid Leukemia,AML)临床大数据及多组学数据库探讨铁死亡相关基因在AML中的作用,并建立铁死亡基因表达相关预后模型。整合TCGA数据库中151例AML患者和GTEx数据库中337例正常人外周血的临床和转录组数据。将Wilcoxon检验和单因素Cox分析结果取交集,筛选出预后相关差异表达基因(Differential Expression Genes, DEGs),使用Lasso回归建立基因标志物预后模型,利用受试者工作特征曲线(Receiver Operating Characteristic Curve,ROC曲线)评价预测价值,Kaplan-Meier法进行生存分析,对AML患者临床数据进行单因素和多因素Cox回归分析,使用差异基因表达分析等方法比较高、低风险患者间的组学差异,最后,利用BeatAML数据库对基因标志物进行验证。将差异基因表达分析和单因素分析结果取交集,得到13个预后相关DEGs。构建了8个基因标志物的预后评分模型,并将患者分为高、低风险两组;ROC曲线分析证实了模型良好的预测性能;生存分析提示高、低风险组患者的生存率具有显著差异;单因素分析显示年龄和风险评分与患者整体生存显著相关,多因素分析显示,年龄和风险评分是独立预后指标。在2个风险组之间筛选出384个DEGs,GO富集分析结果显示,富集的基因大多与中性粒细胞和白细胞的趋化与迁移等免疫相关分子和通路显著相关,KEGG富集通路主要与TNF信号通路、细胞因子与细胞因子受体相互作用相关。BeatAML数据库验证结果显示,5个基因与预后显著相关。铁死亡相关基因在AML中显著表达,且高风险患者预后较差,该研究对AML铁死亡相关潜在生物标志物的发现和应用奠定了一定的基础。     

DNA Methylation Patterns Can Estimate Nonequivalent Outcomes of Breast Cancer with the Same Receptor Subtypes

Breast cancer has various molecular subtypes and displays high heterogeneity. Aberrant DNA methylation is involved in tumor origin, development and progression. Moreover, distinct DNA methylation patterns are associated with specific breast cancer subtypes. We explored DNA methylation patterns in association with gene expression to assess their impact on the prognosis of breast cancer based on Infinium 450K arrays (training set) from The Cancer Genome Atlas (TCGA). The DNA methylation patterns of 12 featured genes that had a high correlation with gene expression were identified through univariate and multivariable Cox proportional hazards models and used to define the methylation risk score (MRS). An improved ability to distinguish the power of the DNA methylation pattern from the 12 featured genes (p = 0.00103) was observed compared with the average methylation levels (p = 0.956) or gene expression (p = 0.909). Furthermore, MRS provided a good prognostic value for breast cancers even when the patients had the same receptor status. We found that ER-, PR- or Her2- samples with high-MRS had the worst 5-year survival rate and overall survival time. An independent test set including 28 patients with death as an outcome was used to test the validity of the MRS of the 12 featured genes; this analysis obtained a prognostic value equivalent to the training set. The predict power was validated through two independent datasets from the GEO database. The DNA methylation pattern is a powerful predictor of breast cancer survival, and can predict outcomes of the same breast cancer molecular subtypes.     

Identification of methylation markers for the prediction of nodal metastasis in oral and oropharyngeal squamous cell carcinoma

Hypermethylation is an important mechanism for the dynamic regulation of gene expression, necessary for metastasizing tumour cells. Our aim is to identify methylation tumour markers that have a predictive value for the presence of regional lymph node metastases in patients with oral and oropharyngeal squamous cell carcinoma (OOSCC). Significantly differentially expressed genes were retrieved from four reported microarray expression profiles comparing pN0 and pN+ head-neck tumours, and one expression array identifying functionally hypermethylated genes. Additional metastasis-associated genes were included from the literature. Thus genes were selected that influence the development of nodal metastases and might be regulated by methylation. Methylation-specific PCR (MSP) primers were designed and tested on 8 head-neck squamous cell carcinoma cell lines and technically validated on 10 formalin-fixed paraffin-embedded (FFPE) OOSCC cases. Predictive value was assessed in a clinical series of 70 FFPE OOSCC with pathologically determined nodal status. Five out of 28 methylation markers (OCLN, CDKN2A, MGMT, MLH1 and DAPK1) were frequently differentially methylated in OOSCC. Of these, MGMT methylation was associated with pN0 status (P = 0.02) and with lower immunoexpression (P = 0.02). DAPK1 methylation was associated with pN+ status (P = 0.008) but did not associate with protein expression. In conclusion, out of 28 candidate genes, two (7%) showed a predictive value for the pN status. Both genes, DAPK1 and MGMT, have predictive value for nodal metastasis in a clinical group of OOSCC. Therefore DNA methylation markers are capable of contributing to diagnosis and treatment selection in OOSCC. To efficiently identify additional new methylation markers, genome-wide methods are needed.     

Identification of metabolism-associated molecular subtype in ovarian cancer

Ovarian cancer (OC) is the most lethal gynaecological cancer with genomic complexity and extensive heterogeneity. This study aimed to characterize the molecular features of OC based on the gene expression profile of 2752 previously characterized metabolism-relevant genes and provide new strategies to improve the clinical status of patients with OC. Finally, three molecular subtypes (C1, C2 and C3) were identified. The C2 subtype displayed the worst prognosis, upregulated immune-cell infiltration status and expression level of immune checkpoint genes, lower burden of copy number gains and losses and suboptimal response to targeted drug bevacizumab. The C1 subtype showed downregulated immune-cell infiltration status and expression level of immune checkpoint genes, the lowest incidence of BRCA mutation and optimal response to targeted drug bevacizumab. The C3 subtype had an intermediate immune status, the highest incidence of BRCA mutation and a secondary optimal response to bevacizumab. Gene signatures of C1 and C2 subtypes with an opposite expression level were mainly enriched in proteolysis and immune-related biological process. The C3 subtype was mainly enriched in the T cell-related biological process. The prognostic and immune status of subtypes were validated in the Gene Expression Omnibus (GEO) dataset, which was predicted with a 45-gene classifier. These findings might improve the understanding of the diversity and therapeutic strategies for OC.     

基于SNP共表达网络肝癌分子分型及预后分析

基于SNP突变数据与mRNA表达谱关联分析,构建一种肝癌分子分型方法并对比不同分型预后的差异,并对不同分型肝癌的发生发展机制进一步研究。首先通过TCGA数据库收集359例肝细胞癌患者的SNP突变数据和mRNA表达数据,采用Wilcoxon秩和检验,筛选突变后差异表达基因,并通过生物信息学工具String和Cytoscape构建差异表达基因的蛋白互作网络,筛选连接度最高的10个Hub基因。利用Consensus Cluster Plus软件包,基于Hub基因mRNA表达水平构建NMF分子分型模型,再结合生存数据评估各分型患者的预后。最后利用加权基因共表达网络分析(WGCNA),识别与肝癌分子分型相关的模块,并针对关键模块的基因进行通路富集,从而对不同分型肝癌的基因表达谱进行比较。结果：NMF模型将肝癌分为高危、低危2个分型,其中CDKN2A和FOXO1基因对分型贡献度高。生存分析显示低危组患者的生存情况显著优于高危组,高危组富集多个与肿瘤细胞侵蚀、转移、复发过程相关的信号通路,低危组则与细胞周期和胰液分泌相关。本研究在无先验性信息的前提下,基于突变后显著差异表达的Hub基因表达水平构建的...     

EYA4 serves as a prognostic biomarker in hepatocellular carcinoma and suppresses tumour angiogenesis and metastasis

Eye absent homolog 4 (EYA4) has been demonstrated to be down‐regulated in hepatocellular carcinoma (HCC), but its biological function and the mechanism in HCC angiogenesis and metastasis remain largely unknown. Herein, we showed that EYA4 expression was frequently low in HCC tissue samples compared with matched adjacent non‐tumourous tissues. In the analysis of 302 HCC specimens, we revealed that decreased expression of EYA4 correlated with tumour differentiation status. Univariate and multivariate analyses identified EYA4 as an independent risk factor for recurrence‐free survival (RFS) and overall survival (OS) among the 302 patients. Functional assays showed that forced expression of EYA4 suppressed HCC cell migration, invasion and capillary tube formation of endothelial cells in vitro, as well as in vivo tumour angiogenesis and metastasis in a mouse model. Furthermore, mechanism study exhibited that EYA4 could inhibit HCC angiogenesis and metastasis by inhibiting c‐JUN/VEGFA pathway. Together, we provide proof that EYA4 is a novel tumour suppressor in HCC and a new prognostic biomarker and therapeutic target in HCC.     

Loss‐of‐function of miR‐142 by hypermethylation promotes TGF‐β‐mediated tumour growth and metastasis in hepatocellular carcinoma

Objectives

Hypermethylation‐induced epigenetic silencing of tumour suppressor genes (TSGs) are frequent events during carcinogenesis. MicroRNA‐142 (miR‐142) is found to be dysregulated in cancer patients to participate into tumour growth, metastasis and angiogenesis. However, the tumour suppressive role of miR‐142 and the status of methylation are not fully understood in hepatocellular carcinoma (HCC).

Methods

Hepatocellular carcinoma tissues and corresponding non‐neoplastic tissues were collected. The expression and function of miR‐142 and TGF‐β in two HCC cell lines were determined. The miRNA‐mRNA network of miR‐142 was analysed in HCC cell lines.

Results

We found that the miR‐142 expression was reduced in tumour tissues and two HCC cell lines HepG2 and SMMC7721, which correlated to higher TNM stage, metastasis and differentiation. Moreover, miR‐142 was identified to directly target and inhibit transforming growth factor β (TGF‐β), leading to decreased cell vitality, proliferation, EMT and the ability of pro‐angiogenesis in TGF‐β‐dependent manner. Interestingly, the status of methylation of miR‐142 was analysed and the results found the hypermethylated miR‐142 in tumour patients and cell lines. The treatment of methylation inhibitor 5‐Aza could restore the expression of miR‐142 to suppress the TGF‐β expression, which impaired TGF‐β‐induced tumour growth.

Conclusion

These findings implicated that miR‐142 was a tumour suppressor gene in HCC and often hyermethylated to increase TGF‐β‐induced development of hepatocellular carcinoma.

     

Mechanisms of ferroptosis with immune infiltration and inflammatory response in rotator cuff injury

The processes driving ferroptosis and rotator cuff (RC) inflammation are yet unknown. The mechanism of ferroptosis and inflammation involved in the development of RC tears was investigated. The Gene Expression Omnibus database was used to obtain the microarray data relevant to the RC tears for further investigation. In this study, we created an RC tears rat model for in vivo experimental validation. For the additional function enrichment analysis, 10 hub ferroptosis-related genes were chosen to construct the correlation regulation network. In RC tears, it was discovered that genes related to hub ferroptosis and hub inflammatory response were strongly correlated. The outcomes of in vivo tests showed that RC tears were related to Cd68-Cxcl13, Acsl4-Sat1, Acsl3-Eno3, Acsl3-Ccr7, and Ccr7-Eno3 pairings in regulating ferroptosis and inflammatory response. Thus, our results show an association between ferroptosis and inflammation, providing a new avenue to explore the clinical treatment of RC tears.     

Construction of a prognostic model for non-small-cell lung cancer based on ferroptosis-related genes

We wished to construct a prognostic model based on ferroptosis-related genes and to simultaneously evaluate the performance of the prognostic model and analyze differences between high-risk and low-risk groups at all levels. The gene-expression profiles and relevant clinical data of patients with non-small-cell lung cancer (NSCLC) were downloaded from public databases. Differentially expressed genes (DEGs) were obtained by analyzing differences between cancer tissues and paracancerous tissues, and common genes between DEGs and ferroptosis-related genes were identified as candidate ferroptosis-related genes. Next, a risk-score model was constructed using univariate Cox analysis and least absolute shrinkage and selection operator (Lasso) analysis. According to the median risk score, samples were divided into high-risk and low-risk groups, and a series of bioinformatics analyses were conducted to verify the predictive ability of the model. Single-sample gene set enrichment analysis (ssGSEA) was used to investigate differences in immune status between high-risk and low-risk groups, and differences in gene mutations between the two groups were investigated. A risk-score model was constructed based on 21 ferroptosis-related genes. A Kaplan–Meier curve and receiver operating characteristic curve showed that the model had good prediction ability. Univariate and multivariate Cox analyses revealed that ferroptosis-related genes associated with the prognosis may be used as independent prognostic factors for the overall survival time of NSCLC patients. The pathways enriched with DEGs in low-risk and high-risk groups were analyzed, and the enriched pathways were correlated significantly with immunosuppressive status.     

Identification of immune infiltration landscape on prognosis and therapy of the ferroptosis-related genes signature in breast cancer

Ferroptosis is a unique iron-dependent cell death mechanism characterized by the generation of lipid reactive oxygen species (ROS) in cancer cells, which leads to mitochondrial metabolic dysregulation. However, how could the tumor immune microenvironment (TIME) modulates ferroptosis remains unclear. Thus, by integrating multiple algorithms, we revealed the novel functional and immune patterns of the ferroptosis-related genes (FRGs) in breast cancer. Five prognostic FRGs were finally selected for the prognostic signature and four of which were identified as the independent biomarkers for immunotherapies. The consensus cluster analysis illustrated the FRGs were characterized by the metabolism dysfunction and immune infiltration cells, meanwhile, these FRGs have the same stem cell characteristics and response efficacy to the immunotherapies. In conclusion, a comprehensive analysis of the FRGs in breast cancer was conducted to develop a prognostic gene signature. Functional and immunological evidence of vulnerabilities in the interaction between ferroptosis and the TIME was also revealed. Further data and research are required.