Search

Abstract Background The fusiform aneurysm is a nonsaccular dilatation affecting the entire vessel wall over a short distance. Although PDGFRB somatic variants have been identified in fusiform intracranial aneurysms, the molecular and cellular mechanisms driving fusiform intracranial aneurysms due to PDGFRB somatic variants remain poorly understood. Methods In this study, single-cell sequencing and immunofluorescence were employed to investigate the phenotypic changes in smooth muscle cells within fusiform intracranial aneurysms. Whole-exome sequencing revealed the presence of PDGFRB gene mutations in fusiform intracranial aneurysms. Subsequent immunoprecipitation experiments further explored the functional alterations of these mutated PDGFRB proteins. For the common c.1684 mutation site of PDGFRβ, we established mutant smooth muscle cell lines and zebrafish models. These models allowed us to simulate the effects of PDGFRB mutations. We explored the major downstream cellular pathways affected by PDGFRBY562D mutations and evaluated the potential therapeutic effects of Ruxolitinib. Results Single-cell sequencing of two fusiform intracranial aneurysms sample revealed downregulated smooth muscle cell markers and overexpression of inflammation-related markers in vascular smooth muscle cells, which was validated by immunofluorescence staining, indicating smooth muscle cell phenotype modulation is involved in fusiform aneurysm. Whole-exome sequencing was performed on seven intracranial aneurysms (six fusiform and one saccular) and PDGFRB somatic mutations were detected in four fusiform aneurysms. Laser microdissection and Sanger sequencing results indicated that the PDGFRB mutations were present in smooth muscle layer. For the c.1684 (chr5: 149505131) site mutation reported many times, further cell experiments showed that PDGFRBY562D mutations promoted inflammatory-related vascular smooth muscle cell phenotype and JAK-STAT pathway played a crucial role in the process. Notably, transfection of PDGFRBY562D in zebrafish embryos resulted in cerebral vascular anomalies. Ruxolitinib, the JAK inhibitor, could reversed the smooth muscle cells phenotype modulation in vitro and inhibit the vascular anomalies in zebrafish induced by PDGFRB mutation. Conclusion Our findings suggested that PDGFRB somatic variants played a role in regulating smooth muscle cells phenotype modulation in fusiform aneurysms and offered a potential therapeutic option for fusiform aneurysms.

Abstract Background The study aimed to investigate the association between nonalcoholic fatty liver disease (NAFLD) and ischemic stroke events after revascularization in patients with Moyamoya disease (MMD). Methods This study prospectively enrolled 275 MMD patients from September 2020 to December 2021. Patients with alcoholism and other liver diseases were excluded. NAFLD was confirmed by CT imaging or abdominal ultrasonography. Stroke events and modified Rankin Scale (mRS) scores at the latest follow-up were compared between the two groups. Results A total of 275 patients were enrolled in the study, among which 65 were diagnosed with NAFLD. Univariate logistic regression analysis showed that NAFLD (P = 0.029) was related to stroke events. Multivariate logistic regression analysis showed that NAFLD is a predictor of postoperative stroke in MMD patients (OR = 27.145, 95% CI = 2.031–362.81, P = 0.013). Kaplan-Meier analysis showed that compared with MMD patients with NAFLD, patients in the control group had a longer stroke-free time (P = 0.004). Univariate Cox analysis showed that NAFLD (P = 0.016) was associated with ischemic stroke during follow-up in patients with MMD. Multivariate Cox analysis showed that NAFLD was an independent risk factor for stroke in patients with MMD (HR = 10.815, 95% CI = 1.259–92.881, P = 0.030). Furthermore, fewer patients in the NAFLD group had good neurologic status (mRS score ≤ 2) than the control group (P = 0.005). Conclusion NAFLD was an independent risk factor for stroke in patients with MMD after revascularization and worse neurological function outcomes.

Abstract Background Moyamoya disease (MMD) is a rare and complex cerebrovascular disorder characterized by the progressive narrowing of the internal carotid arteries and the formation of compensatory collateral vessels. The etiology of MMD remains enigmatic, making diagnosis and management challenging. The MOYAOMICS project was initiated to investigate the molecular underpinnings of MMD and explore potential diagnostic and therapeutic strategies. Methods The MOYAOMICS project employs a multidisciplinary approach, integrating various omics technologies, including genomics, transcriptomics, proteomics, and metabolomics, to comprehensively examine the molecular signatures associated with MMD pathogenesis. Additionally, we will investigate the potential influence of gut microbiota and brain-gut peptides on MMD development, assessing their suitability as targets for therapeutic strategies and dietary interventions. Radiomics, a specialized field in medical imaging, is utilized to analyze neuroimaging data for early detection and characterization of MMD-related brain changes. Deep learning algorithms are employed to differentiate MMD from other conditions, automating the diagnostic process. We also employ single-cellomics and mass cytometry to precisely study cellular heterogeneity in peripheral blood samples from MMD patients. Conclusions The MOYAOMICS project represents a significant step toward comprehending MMD’s molecular underpinnings. This multidisciplinary approach has the potential to revolutionize early diagnosis, patient stratification, and the development of targeted therapies for MMD. The identification of blood-based biomarkers and the integration of multiple omics data are critical for improving the clinical management of MMD and enhancing patient outcomes for this complex disease.

Abstract At present, there is limited research on the mechanisms underlying moyamoya disease (MMD). Herein, we aimed to determine the role of glutamine in MMD pathogenesis, and 360 adult patients were prospectively enrolled. Human brain microvascular endothelial cells (HBMECs) were subjected to Integrin Subunit Beta 4 (ITGB4) overexpression or knockdown and atorvastatin. We assessed factors associated with various signaling pathways in the context of the endothelial‐to‐mesenchymal transition (EndMT), and the expression level of related proteins was validated in the superficial temporal arteries of patients. We found glutamine levels were positively associated with a greater risk of stroke (OR = 1.599, p = 0.022). After treatment with glutamine, HBMECs exhibited enhanced proliferation, migration, and EndMT, all reversed by ITGB4 knockdown. In ITGB4‐transfected HBMECs, the MAPK–ERK–TGF–β/BMP pathway was activated, with Smad4 knockdown reversing the EndMT. Furthermore, atorvastatin suppressed the EndMT by inhibiting Smad1/5 phosphorylation and promoting Smad4 ubiquitination in ITGB4‐transfected HBMECs. We also found the protein level of ITGB4 was upregulated in the superficial temporal arteries of patients with MMD. In conclusion, our study suggests that glutamine may be an independent risk factor for hemorrhage or infarction in patients with MMD and targeting ITGB4 could potentially be therapeutic approaches for MMD.

Abstract Background Moyamoya disease (MMD) stands as a prominent cause of stroke among children and adolescents in East Asian populations. Although a growing body of evidence suggests that dysregulated inflammation and autoimmune responses might contribute to the development of MMD, a comprehensive and detailed understanding of the alterations in circulating immune cells associated with MMD remains elusive. Methods In this study, we employed a combination of single‐cell RNA sequencing (scRNA‐seq), mass cytometry and RNA‐sequencing techniques to compare immune cell profiles in peripheral blood samples obtained from patients with MMD and age‐matched healthy controls. Results Our investigation unveiled immune dysfunction in MMD patients, primarily characterized by perturbations in T‐cell (TC) subpopulations, including a reduction in effector TCs and an increase in regulatory TCs (Tregs). Additionally, we observed diminished natural killer cells and dendritic cells alongside heightened B cells and monocytes in MMD patients. Notably, within the MMD group, there was an augmented proportion of fragile Tregs, whereas the stable Treg fraction decreased. MMD was also linked to heightened immune activation, as evidenced by elevated expression levels of HLA‐DR and p‐STAT3. Conclusions Our findings offer a comprehensive view of the circulating immune cell landscape in MMD patients. Immune dysregulation in patients with MMD was characterized by alterations in T‐cell populations, including a decrease in effector T‐cells and an increase in regulatory T‐cells (Tregs), suggest a potential role for disrupted circulating immunity in the aetiology of MMD.

Abstract Introduction Increasing evidence has revealed the key activity of protein disulfide isomerase A4 (PDIA4) in the endoplasmic reticulum stress (ERS) response. However, the role of PDIA4 in regulating glioblastoma (GBM)-specific pro-angiogenesis is still unknown. Methods The expression and prognostic role of PDIA4 were analyzed using a bioinformatics approach and were validated in 32 clinical samples and follow-up data. RNA-sequencing was used to search for PDIA4-associated biological processes in GBM cells, and proteomic mass spectrum (MS) analysis was used to screen for potential PDIA4 substrates. Western blotting, real-time quantitative polymerase chain reaction (RT-qPCR), and enzyme-linked immunosorbent assays (ELISA) were used to measure the levels of the involved factors. Cell migration and tube formation assays determined the pro-angiogenesis activity of PDIA4 in vitro. An intracranial U87 xenograft GBM animal model was constructed to evaluate the pro-angiogenesis role of PDIA4 in vivo. Results Aberrant overexpression of PDIA4 was associated with a poor prognosis in patients with GBM, although PDIA4 could also functionally regulate intrinsic GBM secretion of vascular endothelial growth factor-A (VEGF-A) through its active domains of Cys-X-X-Cys (CXXC) oxidoreductase. Functionally, PDIA4 exhibits pro-angiogenesis activity both in vitro and in vivo, and can be upregulated by ERS through transcriptional regulation of X-box binding protein 1 (XBP1). The XBP1/PDIA4/VEGFA axis partially supports the mechanism underlying GBM cell survival under ER stress. Further, GBM cells with higher expression of PDIA4 showed resistance to antiangiogenic therapy in vivo. Conclusions Our findings revealed the pro-angiogenesis role of PDIA4 in GBM progression and its potential impact on GBM survival under a harsh microenvironment. Targeting PDIA4 might help to improve the efficacy of antiangiogenic therapy in patients with GBM.

BackgroundBrain metastases (BMs) are the most common central nervous system (CNS) malignant tumors, with rapid disease progression and extremely poor prognosis. The heterogeneity between primary lung cancers and BMs leads to the divergent efficacy of the adjuvant therapy response to primary tumors and BMs. However, the extent of heterogeneity between primary lung cancers and BMs, and the evolutionary process remains little known.MethodsTo deeply insight into the extent of inter-tumor heterogeneity at a single-patient level and the process of these evolutions, we retrospectively analyzed a total of 26 tumor samples from 10 patients with matched primary lung cancers and BMs. One patient underwent four times brain metastatic lesion surgery with diverse locations and one operation for the primary lesion. The genomic and immune heterogeneity between primary lung cancers and BMs were evaluated by utilizing whole-exome sequencing (WESeq) and immunohistochemical analysis.ResultsIn addition to inheriting genomic phenotype and molecular phenotype from the primary lung cancers, massive unique genomic phenotype and molecular phenotype were also observed in BMs, which revealed unimaginable complexity of tumor evolution and extensive heterogeneity among lesions at a single-patient level. By analysis of a multi-metastases case (Case 3) of cancer cells’ subclonal composition, we found similar multiple subclonal clusters in the four spatial and temporal isolated brain metastatic focus, with the characteristics of polyclonal dissemination. Our study also verified that the expression level of immune checkpoints-related molecule Programmed Death-Ligand 1 (PD-L1) (P = 0.0002) and the density of tumor-infiltrating lymphocytes (TILs) (P = 0.0248) in BMs were significantly lower than that in paired primary lung cancers. Additionally, tumor microvascular density (MVD) also differed between primary tumors and paired BMs, indicating that temporal and spatial diversity profoundly contributes to the evolution of BMs heterogeneity.ConclusionOur study revealed the significance of temporal and spatial factors to the evolution of tumor heterogeneity by multi-dimensional analysis of matched primary lung cancers and BMs, which also provided novel insight for formulating individualized treatment strategies for BMs.

ObjectiveMitochondrial dysfunction and oxidative stress are known to involved in tumor occurrence and progression. This study aimed to explore the molecular subtypes of lower-grade gliomas (LGGs) based on oxidative stress-related and mitochondrial-related genes (OMRGs) and construct a prognostic model for predicting prognosis and therapeutic response in LGG patients.MethodsA total of 223 OMRGs were identified by the overlap of oxidative stress-related genes (ORGs) and mitochondrial-related genes (MRGs). Using consensus clustering analysis, we identified molecular subtypes of LGG samples from TCGA database and confirmed the differentially expressed genes (DEGs) between clusters. We constructed a risk score model using LASSO regression and analyzed the immune-related profiles and drug sensitivity of different risk groups. The prognostic role of the risk score was confirmed using Cox regression and Kaplan-Meier curves, and a nomogram model was constructed to predict OS rates. We validated the prognostic role of OMRG-related risk score in three external datasets. Quantitative real-time PCR (qRT-PCR) and immunohistochemistry (IHC) staining confirmed the expression of selected genes. Furthermore, wound healing and transwell assays were performed to confirm the gene function in glioma.ResultsWe identified two OMRG-related clusters and cluster 1 was significantly associated with poor outcomes (P

ObjectiveMethionine sulfoxide (MetO) has been identified as a risk factor for vascular diseases and was considered as an important indicator of oxidative stress. However, the effects of MetO and its association with moyamoya disease (MMD) remained unclear. Therefore, we performed this study to evaluate the association between serum MetO levels and the risk of MMD and its subtypes.MethodsWe eventually included consecutive 353 MMD patients and 88 healthy controls (HCs) with complete data from September 2020 to December 2021 in our analyzes. Serum levels of MetO were quantified using liquid chromatography-mass spectrometry (LC–MS) analysis. We evaluated the role of MetO in MMD using logistic regression models and confirmed by receiver-operating characteristic (ROC) curves and area under curve (AUC) values.ResultsWe found that the levels of MetO were significantly higher in MMD and its subtypes than in HCs (p

ObjectiveIncreasing studies have indicated that senescence was associated with tumorigenesis and progression. Lower-grade glioma (LGG) presented a less invasive nature, however, its treatment efficacy and prognosis prediction remained challenging due to the intrinsic heterogeneity. Therefore, we established a senescence-related signature and investigated its prognostic role in LGGs.MethodsThe gene expression data and clinicopathologic features were from The Cancer Genome Atlas (TCGA) database. The experimentally validated senescence genes (SnGs) from the CellAge database were obtained. Then LASSO regression has been performed to build a prognostic model. Cox regression and Kaplan-Meier survival curves were performed to investigate the prognostic value of the SnG-risk score. A nomogram model has been constructed for outcome prediction. Immunological analyses were further performed. Data from the Chinese Glioma Genome Atlas (CGGA), Repository of Molecular Brain Neoplasia Data (REMBRANDT), and GSE16011 were used for validation.ResultsThe 6-SnG signature has been established. The results showed SnG-risk score could be considered as an independent predictor for LGG patients (HR=2.763, 95%CI=1.660-4.599, P

Abstract Background Glioma is the most common intracranial tumor. The inflammatory response actively participates in the malignancy of gliomas. There is still limited knowledge about the biological function of immune‐related genes (IRGs) and their potential involvement in the malignancy of gliomas. Methods We screened differentially expressed and survival‐associated IRGs, and explored their potential molecular characteristics. Then we developed a prognostic index derived from seven hub IRGs. A prognostic nomogram was built to indicate the prognostic value of the prognostic index and seven IRGs. We characterized the immune infiltration landscape to analyze tumor‐immune interactions. The real‐time quantitative polymerase chain reaction assay was performed to validate bioinformatics results. Results The differentially expressed IRGs are involved in cell chemotaxis, cytokine activity, and the chemokine‐mediated signaling pathway. The prognostic index derived from seven IRGs had clinical prognostic value in glioma, and positively correlated with the malignant clinicopathological characteristics. A nomogram further indicated that the prognostic index and seven hub IRGs had clinical prognostic value for gliomas. We revealed that the prognostic index could reflect the state of the glioma immune microenvironment. Conclusion This study demonstrates the importance of an IRG‐based prognostic index as a potential biomarker for predicting malignancy in gliomas.

Objective: Glioma is the most prevalent and fatal intracranial malignant tumor. Extracellular matrix protein 2 (ECM2) has rarely been studied in gliomas. Therefore, we explored the role of ECM2 in lower-grade gliomas (LGGs). Methods: The RNA-seq and clinicopathology data were obtained from the TCGA database. The immunohistochemical (IHC) staining was used to verify the expression of ECM2. Functional enrichment analyses, immune-related analyses, drug sensitivity, and mutation profile analyses were further conducted. Cox regression and Kaplan–Meier curves were utilized for survival analyses, while four external datasets were used to validate the prognostic role of ECM2. Furthermore, qRT-PCR, CCK-8, wound healing, and transwell assays were performed to confirm the function of ECM2 in gliomas. Results: The study found a significant upregulation of ECM2 expression with increasing glioma grades and a significant association between ECM2 expression and tumor immune infiltration. Cox regression verified the prognostic role of ECM2 in LGG patients (HR = 1.656, 95%CI = 1.055–2.600, p = 0.028). High ECM2 expression was significantly associated with poor outcome (p < 0.001). Four external datasets validated its prognostic value. After the knockdown of ECM2, the functional experiments showed a significant decrease in proliferation, migration, and invasion in glioma cell lines. Conclusion: The study suggested the potential of ECM2 as a novel immune-associated prognostic biomarker and therapeutic target for glioma patients.

Immune-related gene pairs (IRGPs) have been associated with prognosis in various cancer types, but few studies have examined their prognostic capabilities in glioma patients. Here, we gathered the gene expression and clinical profile data of primary lower-grade glioma (LGG) patients from The Cancer Genome Atlas (TCGA), the Chinese Glioma Genome Atlas (CGGA, containing CGGAseq1 and CGGAseq2), the Gene Expression Omnibus (GEO: GSE16011), and Rembrandt datasets. In the TCGA dataset, univariate Cox regression was performed to detect overall survival (OS)-related IRGs, Lasso regression, and multivariate Cox regression were used to screen robust prognosis-related IRGs, and 19 IRGs were selected for the construction of an IRGP prognostic signature. All patients were allotted to high- and low-risk subgroups based on the TCGA dataset median value risk score. Validation analysis indicated that the IRGP signature returned a stable prognostic value among all datasets. Univariate and multivariate Cox regression analyses indicated that the IRG -signature could efficiently predict the prognosis of primary LGG patients. The IRGP-signature-based nomogram model was built, revealing the reliable ability of the IRGP signature to predict clinical prognosis. The single-sample gene set enrichment analysis (ssGSEA) suggested that high-risk samples contained higher numbers of immune cells but featured lower tumor purity than low-risk samples. Finally, we verified the prognostic ability of the IRGP signature using experiments performed in LGG cells. These results indicated that the IRGP signature could be regarded as a stable prognostic assessment predictor for identifying high-risk primary LGG patients.

Lysine acetylation modification, which has key roles in cellular homeostasis as well as cancer malignancy, is dynamically regulated by lysine acetylation regulators (LARs). In our study, we found that most of 33 evaluated LARs were differentially expressed among 1,125 gliomas grouped by different clinicopathological characteristics. Consensus clustering was applied to 33 LARs, resulting in three glioma subtypes (LA1, 2, and 3). The LA3 subgroup was associated with the poorest clinical outcome, higher WHO grade, fewer isocitrate dehydrogenase mutations, and lower frequency of 1p/19q codeletion. Furthermore, gene set enrichment analysis indicated that eight tumor hallmarks were highly enriched in the LA3 subgroup. These results suggested that LARs are significantly related to glioma malignancy. We then designed a LAR-signature based on 14 overall survival (overall survival)-related LARs, and showed that the LAR-signature possesses strong and independent prognostic value for glioma patients in both training and validation datasets. Moreover, by interrogating single nucleotide polymorphism and copy number variation (CNV) data in The Cancer Genome Atlas dataset, we found that higher score of our risk signature is correlated with the hypermutation status of gliomas and that HDAC1(1p) was one of the oncogenes lost in 1p/19q codeletion events, while SIRT2(19q) and EP300(22q) may act as tumor suppressors in gliomas with 19q or 22q deletions, respectively. In conclusion, LARs are critical for the malignant development of gliomas, and our results are useful for prognostic stratification and development of novel assessment strategies for the prognosis of glioma patients.

RNA binding proteins (RBPs) have been reported to be involved in cancer malignancy but related functions in glioma have been less studied. Herein, we screened 14 prognostic RBP genes and constructed a risk signature to predict the prognosis of glioma patients. Univariate Cox regression was used to identify overall survival (OS)-related RBP genes. Prognostic RBP genes were screened and used to establish the RBP-signature using the least absolute shrinkage and selection operator (Lasso) method in The Cancer Genome Atlas (TCGA) cohort. The 14 RBP genes signature showed robust and stable prognostic value in the TCGA training (n = 562) cohort and in three independent validation cohorts (Chinese Glioma Genome Atlas [CGGA]seq1, CGGAseq2, and GSE16011 datasets comprising 303, 619, and 250 glioma patients, respectively). Risk scores were calculated for each patient and high-risk gliomas were defined by the median risk score in each cohort. Survival analysis in subgroups of glioma patients showed that the RBP-signature retained its prognostic value in low-grade gliomas (LGGs) and glioblastomas (GBM)s. Univariate and multivariate Cox regression analysis in each dataset and the meta cohort revealed that the RBP-signature stratification could efficiently recognize high-risk gliomas [Hazard Ratio (HR):3.662, 95% confidence interval (CI): 3.187–4.208, p < 0.001] and was an independent prognostic factor for OS (HR:1.594, 95% CI: 1.244–2.043, p < 0.001). Biological process and KEGG pathway analysis revealed the RBP gene signature was associated with immune cell activation, the p53 signaling pathway, and the PI3K-Akt signaling pathway and so on. Moreover, a nomogram model was constructed for clinical application of the RBP-signature, which showed stable predictive ability. In summary, the RBP-signature could be a robust indicator for prognostic evaluation and identifying high-risk glioma patients.

Background: Glioma is the most common primary intracranial tumor, accounting for the vast majority of intracranial malignant tumors. Aberrant expression of RNA:5-methylcytosine(m5C) methyltransferases have recently been the focus of research relating to the occurrence and progression of tumors. However, the prognostic value of RNA:m5C methyltransferases in glioma remains unclear. This study investigated RNA: m5C methyltransferase expression and defined its clinicopathological signature and prognostic value in gliomas.Methods: We obtained the RNA-sequence and Clinicopathological data of RNA:m5C methyltransferases underlying gliomas from the Chinese Glioma Genome Atlas (CGGA) and The Cancer Genome Atlas (TCGA) datasets. We analyzed the expression of RNA:m5C methyltransferase genes in gliomas with different clinicopathological characteristics and identified different subtypes using Consensus clustering analysis. Gene Ontology (GO) and Gene Set Enrichment Analysis (GSEA) was used to annotate the function of these genes. Univariate Cox regression and the least absolute shrinkage and selection operator (LASSO) Cox regression algorithm analyses were performed to construct the risk signature. Kaplan-Meier method and Receiver operating characteristic (ROC) curves were used to assess the overall survival of glioma patients. Additionally, Cox proportional regression model analysis was developed to address the connections between the risk scores and clinical factors.Results: We revealed the differential expression of RNA:m5C methyltransferase genes in gliomas with different clinicopathological features. Consensus clustering of RNA:m5C methyltransferases identified three clusters of gliomas with different prognostic and clinicopathological features. Meanwhile, functional annotations demonstrated that RNA:m5C methyltransferases were significantly associated with the malignant progression of gliomas. Thereafter, five RNA:m5C methyltransferase genes were screened to construct a risk signature that can be used to predict not only overall survival but also clinicopathological features in gliomas. ROC curves revealed the significant prognostic ability of this signature. In addition, Multivariate Cox regression analyses indicated that the risk score was an independent prognostic factor for glioma outcome.Conclusion: We demonstrated the prognostic role of RNA:m5C methyltransferases in the initiation and progression of glioma. We have expanded on the understanding of the molecular mechanism involved, and provided a unique approach to predictive biomarkers and targeted therapy for gliomas.

The prognostic value of N6-methylandenosine-related long non-coding RNAs (m6A-related lncRNAs) was investigated in 646 lower-grade glioma (LGG) samples from The Cancer Genome Atlas (TCGA) and the Chinese Glioma Genome Atlas (CGGA) datasets. We implemented Pearson correlation analysis to explore the m6A-related lncRNAs, and then univariate Cox regression analysis was performed to screen their prognostic roles in LGG patients. Twenty-four prognostic m6A-related lncRNAs were identified as prognostic lncRNAs and they were inputted in a least absolute shrinkage and selection operator (LASSO) Cox regression to establish a m6A-related lncRNA prognostic signature (m6A-LPS, including 9 m6A-related prognostic lncRNAs) in the TCGA dataset. Corresponding risk scores of patients were calculated and divided LGG patients into low- and high-risk subgroups by the median value of risk scores in each dataset. The m6A-LPS was validated in the CGGA dataset and it showed a robust prognostic ability in the stratification analysis. Principal component analysis showed that the low- and high-risk subgroups had distinct m6A status. Enrichment analysis indicated that malignancy-associated biological processes, pathways and hallmarks were more common in the high-risk subgroup. Moreover, we constructed a nomogram (based on m6A-LPS, age and World Health Organization grade) that had a strong ability to forecast the overall survival (OS) of the LGG patients in both datasets. We also establish a competing endogenous RNA (ceRNA) network based on seven of the twenty-four m6A-related lncRNAs. Besides, we also detected five m6A-related lncRNA expression levels in 22 clinical samples using quantitative real-time polymerase chain reaction assay.

Background: Epithelial-mesenchymal transition (EMT) is regulated by induction factors, transcription factor families and an array of signaling pathways genes, and has been implicated in the invasion and progression of gliomas.Methods: We obtained the Clinicopathological data sets from Chinese Glioma Genome Atlas (CGGA). The “limma” package was used to analyze the expression of EMT-related genes in gliomas with different pathological characteristics. We used the “ConsensusClusterPlus” package to divide gliomas into two groups to study their correlation with glioma malignancy. The least absolute shrinkage and selection operator (LASSO) Cox regression was applied to select seven prognosis-associated genes to build the risk signature, and the coefficients obtained from the LASSO algorithm were used to calculate the risk score which we applied to determine the prognostic value of the risk signature. Univariate and multivariate Cox regression analyses were used to determine whether the risk signature is an independent prognostic indicator.Results: We analyzed the differentially expressed 22 common epithelial-mesenchymal transition-associated genes in 508 gliomas graded by different clinicopathological features. Two glioma subgroups (EM1/2) were identified by consistent clustering of the proteins, of which the EM1 subgroup had a better prognosis than the EM2 subgroup, and the EM2 group was associated with cancer migration and proliferation. Significant enrichment analysis revealed that EMT-related transcriptional regulators and signaling pathways genes were highly related to glioma malignancies. Seven EMT-related genes were used to derive risk scores, which served as independent prognostic markers and prediction factors for the clinicopathological features of glioma. And we found the overall survival (OS) was significantly different between the low- and high-risk groups, the ROC curve indicated that the risk score can predict survival rates for glioma patients.Conclusion: EMT-related induction factors, transcriptional regulators and signaling pathways genes are important players in the malignant progression of glioma and may help in decision making regarding the choice of prognosis assessment and provide us clues to understand EMT epigenetic modification in glioma.