Integrated profiling identifies ferredoxin 1 as an immune-related biomarker of malignant phenotype in glioma.

Background: Glioma, a highly resistant and recurrent type of central nervous system tumor, poses a significant challenge in terms of effective drug treatments and its associated mortality rates. Despite the discovery of Ferredoxin 1 (FDX1) as a crucial participant in cuproptosis, an innovative mechanism of cellular demise, its precise implications for glioma prognosis and tumor immune infiltration remain inadequately elucidated. Methods: To analyze pan-cancer data, we employed multiple public databases. Gene expression evaluation was performed using tissue microarray (TMA) and single-cell sequencing data. Furthermore, four different approaches were employed to assess the prognostic importance of FDX1 in glioma. We conducted the analysis of differential expression genes (DEGs) and Gene Set Enrichment Analysis (GSEA) to identify immune-related predictive signaling pathways. Somatic mutations were assessed using Tumor Mutation Burden (TMB) and waterfall plots. Immune cell infiltration was evaluated with five different algorithms. Furthermore, we performed in vitro investigations to evaluate the biological roles of FDX1 in glioma. Results: Glioma samples exhibited upregulation of FDX1, which in turn predicted poor prognosis and was positively associated with unfavorable clinicopathological characteristics. Notably, the top four enriched signaling pathways were immune-related, and the discovery revealed a connection between the expression of FDX1 and the frequency of mutations or the TMB. The FDX1_high group exhibited heightened infiltration of immune cells, and there existed a direct association between the expression of FDX1 and the regulation of immune checkpoint. In vitro experiments demonstrated that FDX1 knockdown reduced proliferation, migration, invasion and transition from G2 to M phase in glioma cells. Conclusion: In glioma, FDX1 demonstrated a positive association with the advancement of malignancy and changes in the infiltration of immune cells.

Identification of three tumor antigens and immune subtypes for mRNA vaccine development in diffuse glioma.

Rationale: Diffuse glioma patients have high mortality and recurrence despite multimodal therapies. This study aims to identify the potential tumor antigens for mRNA vaccines and subtypes suitable for the immunotherapy of patients with diffuse glioma. Methods: Gene expression profiles and corresponding clinical information were obtained from the Chinese Glioma Genome Atlas (CGGA) and the Cancer Genome Atlas (TCGA) databases. Genetic alterations were extracted from cBioPortal. Differential gene analysis, survival analysis, correlation analysis, consensus clustering analysis, and immune cell infiltration analysis were conducted based on the various databases. Finally, the hub genes, the modules related to tumor antigens, and the immune subtypes were identified using WGCNA method. Results: Three over-expressed, amplified, and mutated tumor antigens, including KDR, COL1A2, and SAMD9, were associated with clinical outcomes. The expression of the three genes had a positive correlation with the abundance of antigen-presenting cells (APCs) and APC marker expression. Subsequently, three immune subtypes (Ims1, Ims2, and Ims3) were distinguished in the TCGA cohort, which exhibited distinct molecular, cellular, and clinical characteristics consistent with the CGGA cohort. Diffuse gliomas with subtype Ims1 were more malignant with immunosuppressive phenotypes and more associated with poor prognosis than the other two subtypes. The three antigens and the immune checkpoints were differentially expressed among the three immune subtypes. Finally, functional enrichment analysis of the genes related to tumor antigens and immune subtypes suggested that they are enriched in many immune-associated processes. Conclusions: KDR, COL1A2, and SAMD9 are potential antigens for developing mRNA vaccines against diffuse glioma. The results suggest that immunotherapy targeting these three antigens is more suitable for patients with subtype Ims1. This study provides insights into immunotherapy for diffuse glioma.

METTL7B is a novel prognostic biomarker of lower-grade glioma based on pan-cancer analysis.

Methyltransferase-like 7B (METTL7B) is a member of the methyltransferase-like protein family that plays an important role in the development and progression of tumors. However, its prognostic value and the correlation of METTL7B expression and tumor immunity in some cancers remain unclear. By analyzing online data, we found that METTL7B is abnormally overexpressed in multiple human tumors and plays an important role in the overall survival (OS) of patients with 8 cancer types and disease-free survival (DFS) of patients with 5 cancer types. Remarkably, METTL7B expression was positively correlated with the OS and DFS of patients with lower-grade glioma (LGG). In addition, a positive correlation between METTL7B expression and immune cell infiltration in LGG was observed. Moreover, we identified a strong correlation between METTL7B expression and immune checkpoint gene expression in kidney chromophobe (KICH), LGG and pheochromocytoma and paraganglioma (PCPG). Furthermore, METTL7B was involved in the extracellular matrix (ECM) and immune-related pathways in LGGs. Finally, in vitro experiments showed that knockdown of METTL7B inhibited the growth, migration, invasion and the epithelial-mesenchymal transition (EMT) of LGG cells. METTL7B expression potentially represents a novel prognostic biomarker due to its significant association with immune cell infiltration in LGG.

Development of an Immune-Related LncRNA Prognostic Signature for Glioma.

INTRODUCTION: Glioma is the most common primary cancer of the central nervous system with dismal prognosis. Long noncoding RNAs (lncRNAs) have been discovered to play key roles in tumorigenesis in various cancers, including glioma. Because of the relevance between immune infiltrating and clinical outcome of glioma, identifying immune-related lncRNAs is urgent for better personalized management. MATERIALS AND METHODS: Single-sample gene set enrichment analysis (ssGSEA) was applied to estimate immune infiltration, and glioma samples were divided into high immune cell infiltration group and low immune cell infiltration group. After screening differentially expressed lncRNAs in two immune groups, least absolute shrinkage and selection operator (LASSO) Cox regression analysis was performed to construct an immune-related prognostic signature. Additionally, we explored the correlation between immune infiltration and the prognostic signature. RESULTS: A total of 653 samples were appropriate for further analyses, and 10 lncRNAs were identified as immune-related lncRNAs in glioma. After univariate Cox regression and LASSO Cox regression analysis, six lncRNAs were identified to construct a prognostic signature for glioma, which could be taken as independent prognostic factors in both univariate and multivariate Cox regression analyses. Moreover, risk score was significantly correlated with all the 29 immune-related checkpoint expression (p < 0.05) in ssGSEA except neutrophils (p = 0.43). CONCLUSION: The study constructed an immune-related prognostic signature for glioma, which contributed to improve clinical outcome prediction and guide immunotherapy.

Identification of collagen genes related to immune infiltration and epithelial-mesenchymal transition in glioma.

BACKGROUND: Gliomas account for the majority of fatal primary brain tumors, and there is much room for research in the underlying pathogenesis, the multistep progression of glioma, and how to improve survival. In our study, we aimed to identify potential biomarkers or therapeutic targets of glioma and study the mechanism underlying the tumor progression. METHODS: We downloaded the microarray datasets (GSE43378 and GSE7696) from the Gene Expression Omnibus (GEO) database. Then, we used weighted gene co-expression network analysis (WGCNA) to screen potential biomarkers or therapeutic targets related to the tumor progression. ESTIMATE (Estimation of STromal and Immune cells in MAlignant Tumors using Expression data) algorithm and TIMER (Tumor Immune Estimation Resource) database were used to analyze the correlation between the selected genes and the tumor microenvironment. Real-time reverse transcription polymerase chain reaction was used to measure the selected gene. Transwell and wound healing assays were used to measure the cell migration and invasion capacity. Western blotting was used to test the expression of epithelial-mesenchymal transition (EMT) related markers. RESULTS: We identified specific module genes that were positively correlated with the WHO grade but negatively correlated with OS of glioma. Importantly, we identified that 6 collagen genes (COL1A1, COL1A2, COL3A1, COL4A1, COL4A2, and COL5A2) could regulate the immunosuppressive microenvironment of glioma. Moreover, we found that these collagen genes were significantly involved in the EMT process of glioma. Finally, taking COL3A1 as a further research object, the results showed that knockdown of COL3A1 significantly inhibited the migration, invasion, and EMT process of SHG44 and A172 cells. CONCLUSIONS: In summary, our study demonstrated that collagen genes play an important role in regulating the immunosuppressive microenvironment and EMT process of glioma and could serve as potential therapeutic targets for glioma management.

Molecular Subtypes Based on the Stemness Index Predict Prognosis in Glioma Patients.

Glioma is the common histological subtype of malignancy in the central nervous system, with high morbidity and mortality. Glioma cancer stem cells (CSCs) play essential roles in tumor recurrence and treatment resistance. Thus, exploring the stem cell-related genes and subtypes in glioma is important. In this study, we collected the RNA-sequencing (RNA-seq) data and clinical information of glioma patients from The Cancer Genome Atlas (TCGA) and Chinese Glioma Genome Atlas (CGGA) databases. With the differentially expressed genes (DEGs) and weighted gene correlation network analysis (WGCNA), we identified 86 mRNA expression-based stemness index (mRNAsi)-related genes in 583 samples from TCGA RNA-seq dataset. Furthermore, these samples from TCGA database could be divided into two significantly different subtypes with different prognoses based on the mRNAsi corresponding gene, which could also be validated in the CGGA database. The clinical characteristics and immune cell infiltrate distribution of the two stemness subtypes are different. Then, functional enrichment analyses were performed to identify the different gene ontology (GO) terms and pathways in the two different subtypes. Moreover, we constructed a stemness subtype-related risk score model and nomogram to predict the prognosis of glioma patients. Finally, we selected one gene (ETV2) from the risk score model for experimental validation. The results showed that ETV2 can contribute to the invasion, migration, and epithelial-mesenchymal transition (EMT) process of glioma. In conclusion, we identified two distinct molecular subtypes and potential therapeutic targets of glioma, which could provide new insights for the development of precision diagnosis and prognostic prediction for glioma patients.

Multi-Omics Analysis Based on Genomic Instability for Prognostic Prediction in Lower-Grade Glioma.

Background: Lower-grade gliomas (LGGs) are a heterogeneous set of gliomas. One of the primary sources of glioma heterogeneity is genomic instability, a novel characteristic of cancer. It has been reported that long noncoding RNAs (lncRNAs) play an essential role in regulating genomic stability. However, the potential relationship between genomic instability and lncRNA in LGGs and its prognostic value is unclear. Methods: In this study, the LGG samples from The Cancer Genome Atlas (TCGA) were divided into two clusters by integrating the lncRNA expression profile and somatic mutation data using hierarchical clustering. Then, with the differentially expressed lncRNAs between these two clusters, we identified genomic instability-related lncRNAs (GInLncRNAs) in the LGG samples and analyzed their potential function and pathway by co-expression network. Cox and least absolute shrinkage and selection operator (LASSO) regression analyses were conducted to establish a GInLncRNA prognostic signature (GInLncSig), which was assessed by internal and external verification, correlation analysis with somatic mutation, independent prognostic analysis, clinical stratification analysis, and model comparisons. We also established a nomogram to predict the prognosis more accurately. Finally, we performed multi-omics-based analyses to explore the relationship between risk scores and multi-omics data, including immune characteristics, N 6-methyladenosine (m6A), stemness index, drug sensitivity, and gene set enrichment analysis (GSEA). Results: We identified 52 GInLncRNAs and screened five from them to construct the GInLncSig model (CRNDE, AC025171.5, AL390755.1, AL049749.1, and TGFB2-AS1), which could independently and accurately predict the outcome of patients with LGG. The GInLncSig model was significantly associated with somatic mutation and outperformed other published signatures. GSEA revealed that metabolic pathways, immune pathways, and cancer pathways were enriched in the high-risk group. Multi-omics-based analyses revealed that T-cell functions, m6A statuses, and stemness characteristics were significantly disparate between two risk subgroups, and immune checkpoints such as PD-L1, PDCD1LG2, and HAVCR2 were significantly highly expressed in the high-risk group. The expression of GInLncSig prognostic genes dramatically correlated with the sensitivity of tumor cells to chemotherapy drugs. Conclusion: A novel signature composed of five GInLncRNAs can be utilized to predict prognosis and impact the immune status, m6A status, and stemness characteristics in LGG. Furthermore, these lncRNAs may be potential and alternative therapeutic targets.

Development and Validation of a Tumor Mutation Burden-Related Immune Prognostic Model for Lower-Grade Glioma.

Tumor mutation burden (TMB) is a useful biomarker to predict prognosis and the efficacy of immune checkpoint inhibitors (ICIs). In this study, we aimed to explore the prognostic value of TMB and the potential association between TMB and immune infiltration in lower-grade gliomas (LGGs). Somatic mutation and RNA-sequencing (RNA-seq) data were downloaded from the Cancer Genome Atlas (TCGA) database. TMB was calculated and patients were divided into high- and low-TMB groups. After performing differential analysis between high- and low-risk groups, we identified six hub TMB and immune-related genes that were correlated with overall survival in LGGs. Then, Gene Set Enrichment Analysis was performed to screen significantly enriched GO terms between the two groups. Moreover, an immune-related risk score system was developed by LASSO Cox analysis based on the six hub genes and was validated with the Chinese Glioma Genome Atlas dataset. Using the TIMER database, we further systematically analyzed the relationships between mutants of the six hub genes and immune infiltration levels, as well as the relationships between the immune-related risk score system and the immune microenvironment in LGGs. The results showed that TMB was negatively correlated with OS and high TMB might inhibit immune infiltration in LGGs. Furthermore, the risk score system could effectively stratify patients into low- and high-risk groups in both the training and validation datasets. Multivariate Cox analysis demonstrated that TMB was not an independent prognostic factor, but the risk score was. Higher infiltration of immune cells (B cells, CD4+ T cells, CD8+ T cells, neutrophils, macrophages, and dendritic cells) and higher levels of immune checkpoints (PD-1, CTLA-4, LAG-3, and TIM-3) were found in patients in the high-risk group. Finally, a novel nomogram model was constructed and evaluated to estimate the overall survival of LGG patients. In summary, our study provided new insights into immune infiltration in the tumor microenvironment and immunotherapies for LGGs.

Emerging insights into molecular mechanisms underlying pyroptosis and functions of inflammasomes in diseases.

Pyroptosis is a form of necrotic and inflammatory programmed cell death, which could be characterized by cell swelling, pore formation on plasma membranes, and release of proinflammatory cytokines (IL-1ß and IL-18). The process of pyroptosis presents as dual effects: protecting multicellular organisms from microbial infection and endogenous dangers; leading to pathological inflammation if overactivated. Two pathways have been found to trigger pyroptosis: caspase-1 mediated inflammasome pathway with the involvement of NLRP1-, NLRP3-, NLRC4-, AIM2-, pyrin-inflammasome (canonical inflammasome pathway) and caspase-4/5/11-mediated inflammasome pathway (noncanonical inflammasome pathway). Gasdermin D (GSDMD) has been proved to be a substrate of inflammatory caspases (caspase-1/4/5/11), and the cleaved N-terminal domain of GSDMD oligomerizes to form cytotoxic pores on the plasma membrane. Here, we mainly reviewed the up to date mechanisms of pyroptosis, and began with the inflammasomes as the activator of caspase-1/caspase-11, 4, and 5. We further discussed these inflammasomes functions in diseases, including infectious diseases, sepsis, inflammatory autoimmune diseases, and neuroinflammatory diseases.