DNA methylation, transcriptome and genetic copy number signatures of diffuse cerebral WHO grade II/III gliomas resolve cancer heterogeneity and development

Pan-cancer analysis identifies olfactory receptor family 7 subfamily A member 5 as a potential biomarker for glioma

Background

Human olfactory receptors (ORs) account for approximately 60% of all human G protein-coupled receptors. The functions of ORs extend beyond olfactory perception and have garnered significant attention in tumor biology. However, a comprehensive pan-cancer analysis of ORs in human cancers is lacking.

Methods

Using data from public databases, such as HPA, TCGA, GEO, GTEx, TIMER2, TISDB, UALCAN, GEPIA2, and GSCA, this study investigated the role of olfactory receptor family 7 subfamily A member 5 (OR7A5) in various cancers. Functional analysis of OR7A5 in LGG and GBM was performed using the CGGA database. Molecular and cellular experiments were performed to validate the expression and biological function of OR7A5 in gliomas.

Results

The results revealed heightened OR7A5 expression in certain tumors, correlating with the expression levels of immune checkpoints and immune infiltration. In patients with gliomas, the expression levels of OR7A5 were closely associated with adverse prognosis, 1p/19p co-deletion status, and wild-type IDH status. Finally, in vitro experiments confirmed the inhibitory effect of OR7A5 knockdown on the proliferative capacity of glioma cells and on the expression levels of proteins related to lipid metabolism.

Conclusion

This study establishes OR7A5 as a novel biomarker, potentially offering a novel therapeutic target for gliomas.

Pediatric diffuse intrinsic pontine glioma radiotherapy response prediction: MRI morphology and T2 intensity-based quantitative analyses

OBJECTIVES

An easy-to-implement MRI model for predicting partial response (PR) postradiotherapy for diffuse intrinsic pontine glioma (DIPG) is lacking. Utilizing quantitative T2 signal intensity and introducing a visual evaluation method based on T2 signal intensity heterogeneity, and compared MRI radiomic models for predicting radiotherapy response in pediatric patients with DIPG.

METHODS

We retrospectively included patients with brainstem gliomas aged ≤ 18 years admitted between July 2011 and March 2023. Applying Response Assessment in Pediatric Neuro-Oncology criteria, we categorized patients into PR and non-PR groups. For qualitative analysis, tumor heterogeneity vision was classified into four grades based on T2-weighted images. Quantitative analysis included the relative T2 signal intensity ratio (rT2SR), extra pons volume ratio, and tumor ring-enhancement volume. Radiomic features were extracted from T2-weighted and T1-enhanced images of volumes of interest. Univariate analysis was used to identify independent variables related to PR. Multivariate logistic regression was performed using significant variables (p < 0.05) from univariate analysis.

RESULTS

Of 140 patients (training n = 109, and test n = 31), 64 (45.7%) achieved PR. The AUC of the predictive model with extrapontine volume ratio, rT2SRmax-min (rT2SRdif), and grade was 0.89. The AUCs of the T2-weighted and T1WI-enhanced models with radiomic signatures were 0.84 and 0.81, respectively. For the 31 DIPG test sets, the AUCs were 0.91, 0.83, and 0.81, for the models incorporating the quantitative features, radiomic model (T2-weighted images, and T1W1-enhanced images), respectively.

CONCLUSION

Combining T2-weighted quantification with qualitative and extrapontine volume ratios reliably predicted pediatric DIPG radiotherapy response.

CLINICAL RELEVANCE STATEMENT

Combining T2-weighted quantification with qualitative and extrapontine volume ratios can accurately predict diffuse intrinsic pontine glioma (DIPG) radiotherapy response, which may facilitate personalized treatment and prognostic assessment for patients with DIPG.

KEY POINTS

Early identification is crucial for radiotherapy response and risk stratification in diffuse intrinsic pontine glioma. The model using tumor heterogeneity and quantitative T2 signal metrics achieved an AUC of 0.91. Using a combination of parameters can effectively predict radiotherapy response in this population.

CD57 defines a novel cancer stem cell that drive invasion of diffuse pediatric-type high grade gliomas

BACKGROUND

Diffuse invasion remains a primary cause of treatment failure in pediatric high-grade glioma (pHGG). Identifying cellular driver(s) of pHGG invasion is needed for anti-invasion therapies.

METHODS

Ten highly invasive patient-derived orthotopic xenograft (PDOX) models of pHGG were subjected to isolation of matching pairs of invasive (HGGINV) and tumor core (HGGTC) cells.

RESULTS

pHGGINV cells were intrinsically more invasive than their matching pHGGTC cells. CSC profiling revealed co-positivity of CD133 and CD57 and identified CD57+CD133- cells as the most abundant CSCs in the invasive front. In addition to discovering a new order of self-renewal capacities, i.e., CD57+CD133- > CD57+CD133+ > CD57-CD133+ > CD57-CD133- cells, we showed that CSC hierarchy was impacted by their spatial locations, and the highest self-renewal capacities were found in CD57+CD133- cells in the HGGINV front (HGGINV/CD57+CD133- cells) mediated by NANOG and SHH over-expression. Direct implantation of CD57+ (CD57+/CD133- and CD57+/CD133+) cells into mouse brains reconstituted diffusely invasion, while depleting CD57+ cells (i.e., CD57-CD133+) abrogated pHGG invasion.

CONCLUSION

We revealed significantly increased invasive capacities in HGGINV cells, confirmed CD57 as a novel glioma stem cell marker, identified CD57+CD133- and CD57+CD133+ cells as a new cellular driver of pHGG invasion and suggested a new dual-mode hierarchy of HGG stem cells.

Integration of Multi-Omics Data for the Classification of Glioma Types and Identification of Novel Biomarkers

Glioma is currently one of the most prevalent types of primary brain cancer. Given its high level of heterogeneity along with the complex biological molecular markers, many efforts have been made to accurately classify the type of glioma in each patient, which, in turn, is critical to improve early diagnosis and increase survival. Nonetheless, as a result of the fast-growing technological advances in high-throughput sequencing and evolving molecular understanding of glioma biology, its classification has been recently subject to significant alterations. In this study, we integrate multiple glioma omics modalities (including mRNA, DNA methylation, and miRNA) from The Cancer Genome Atlas (TCGA), while using the revised glioma reclassified labels, with a supervised method based on sparse canonical correlation analysis (DIABLO) to discriminate between glioma types. We were able to find a set of highly correlated features distinguishing glioblastoma from lower-grade gliomas (LGGs) that were mainly associated with the disruption of receptor tyrosine kinases signaling pathways and extracellular matrix organization and remodeling. Concurrently, the discrimination of the LGG types was characterized primarily by features involved in ubiquitination and DNA transcription processes. Furthermore, we could identify several novel glioma biomarkers likely helpful in both diagnosis and prognosis of the patients, including the genes PPP1R8, GPBP1L1, KIAA1614, C14orf23, CCDC77, BVES, EXD3, CD300A, and HEPN1. Collectively, this comprehensive approach not only allowed a highly accurate discrimination of the different TCGA glioma patients but also presented a step forward in advancing our comprehension of the underlying molecular mechanisms driving glioma heterogeneity. Ultimately, our study also revealed novel candidate biomarkers that might constitute potential therapeutic targets, marking a significant stride toward personalized and more effective treatment strategies for patients with glioma.

Deciphering Glioblastoma: Fundamental and Novel Insights into the Biology and Therapeutic Strategies of Gliomas

Gliomas constitute a diverse and complex array of tumors within the central nervous system (CNS), characterized by a wide range of prognostic outcomes and responses to therapeutic interventions. This literature review endeavors to conduct a thorough investigation of gliomas, with a particular emphasis on glioblastoma (GBM), beginning with their classification and epidemiological characteristics, evaluating their relative importance within the CNS tumor spectrum. We examine the immunological context of gliomas, unveiling the intricate immune environment and its ramifications for disease progression and therapeutic strategies. Moreover, we accentuate critical developments in understanding tumor behavior, focusing on recent research breakthroughs in treatment responses and the elucidation of cellular signaling pathways. Analyzing the most novel transcriptomic studies, we investigate the variations in gene expression patterns in glioma cells, assessing the prognostic and therapeutic implications of these genetic alterations. Furthermore, the role of epigenetic modifications in the pathogenesis of gliomas is underscored, suggesting that such changes are fundamental to tumor evolution and possible therapeutic advancements. In the end, this comparative oncological analysis situates GBM within the wider context of neoplasms, delineating both distinct and shared characteristics with other types of tumors.

Temporal changes of gene expression in health, schizophrenia, bipolar disorder, and major depressive disorder

The molecular events underlying the development, manifestation, and course of schizophrenia, bipolar disorder, and major depressive disorder span from embryonic life to advanced age. However, little is known about the early dynamics of gene expression in these disorders due to their relatively late manifestation. To address this, we conducted a secondary analysis of post-mortem prefrontal cortex datasets using bioinformatics and machine learning techniques to identify differentially expressed gene modules associated with aging and the diseases, determine their time-perturbation points, and assess enrichment with expression quantitative trait loci (eQTL) genes. Our findings revealed early, mid, and late deregulation of expression of functional gene modules involved in neurodevelopment, plasticity, homeostasis, and immune response. This supports the hypothesis that multiple hits throughout life contribute to disease manifestation rather than a single early-life event. Moreover, the time-perturbed functional gene modules were associated with genetic loci affecting gene expression, highlighting the role of genetic factors in gene expression dynamics and the development of disease phenotypes. Our findings emphasize the importance of investigating time-dependent perturbations in gene expression before the age of onset in elucidating the molecular mechanisms of psychiatric disorders.

Improved prognostic stratification of patients with isocitrate dehydrogenase-mutant astrocytoma

Prognostic factors and standards of care for astrocytoma, isocitrate dehydrogenase (IDH)-mutant, CNS WHO grade 4, remain poorly defined. Here we sought to explore disease characteristics, prognostic markers, and outcome in patients with this newly defined tumor type. We determined molecular biomarkers and assembled clinical and outcome data in patients with IDH-mutant astrocytomas confirmed by central pathology review. Patients were identified in the German Glioma Network cohort study; additional cohorts of patients with CNS WHO grade 4 tumors were identified retrospectively at two sites. In total, 258 patients with IDH-mutant astrocytomas (114 CNS WHO grade 2, 73 CNS WHO grade 3, 71 CNS WHO grade 4) were studied. The median age at diagnosis was similar for all grades. Karnofsky performance status at diagnosis inversely correlated with CNS WHO grade (p < 0.001). Despite more intensive treatment upfront with higher grade, CNS WHO grade was strongly prognostic: median overall survival was not reached for grade 2 (median follow-up 10.4 years), 8.1 years (95% CI 5.4-10.8) for grade 3, and 4.7 years (95% CI 3.4-6.0) for grade 4. Among patients with CNS WHO grade 4 astrocytoma, median overall survival was 5.5 years (95% CI 4.3-6.7) without (n = 58) versus 1.8 years (95% CI 0-4.1) with (n = 12) homozygous CDKN2A deletion. Lower levels of global DNA methylation as detected by LINE-1 methylation analysis were strongly associated with CNS WHO grade 4 (p < 0.001) and poor outcome. MGMT promoter methylation status was not prognostic for overall survival. Histomolecular stratification based on CNS WHO grade, LINE-1 methylation level, and CDKN2A status revealed four subgroups of patients with significantly different outcomes. In conclusion, CNS WHO grade, global DNA methylation status, and CDKN2A homozygous deletion are prognostic in patients with IDH-mutant astrocytoma. Combination of these parameters allows for improved prediction of outcome. These data aid in designing upcoming trials using IDH inhibitors.

Genomic Profiling of Lower-Grade Gliomas Subtype with Distinct Molecular and Clinicopathologic Characteristics via Altered DNA-Damage Repair Features

Lower WHO grade II and III gliomas (LGGs) exhibit significant genetic and transcriptional heterogeneity, and the heterogeneity of DNA damage repair (DDR) and its relationship to tumor biology, transcriptome, and tumor microenvironment (TME) remains poorly understood. In this study, we conducted multi-omics data integration to investigate DDR alterations in LGG. Based on clinical parameters and molecular characteristics, LGG patients were categorized into distinct DDR subtypes, namely, DDR-activated and DDR-suppressed subtypes. We compared gene mutation, immune spectrum, and immune cell infiltration between the two subtypes. DDR scores were generated to classify LGG patients based on DDR subtype features, and the results were validated using a multi-layer data cohort. We found that DDR activation was associated with poorer overall survival and that clinicopathological features of advanced age and higher grade were more common in the DDR-activated subtype. DDR-suppressed subtypes exhibited more frequent mutations in IDH1. In addition, we observed significant upregulation of activated immune cells in the DDR-activated subgroup, which suggests that immune cell infiltration significantly influences tumor progression and immunotherapeutic responses. Furthermore, we constructed a DDR signature for LGG using six DDR genes, which allowed for the division of patients into low- and high-risk groups. Quantitative real-time PCR results showed that CDK1, CDK2, TYMS, SMC4, and WEE1 were significantly upregulated in LGG samples compared to normal brain tissue samples. Overall, our study sheds light on DDR heterogeneity in LGG and provides insight into the molecular pathways of DDR involved in LGG development.

Transcriptional states of CAR-T infusion relate to neurotoxicity – lessons from high-resolution single-cell SOM expression portraying

Anti-CD19 CAR-T cell immunotherapy is a hopeful treatment option for patients with B cell lymphomas, however it copes with partly severe adverse effects like neurotoxicity. Single-cell resolved molecular data sets in combination with clinical parametrization allow for comprehensive characterization of cellular subpopulations, their transcriptomic states, and their relation to the adverse effects. We here present a re-analysis of single-cell RNA sequencing data of 24 patients comprising more than 130,000 cells with focus on cellular states and their association to immune cell related neurotoxicity. For this, we developed a single-cell data portraying workflow to disentangle the transcriptional state space with single-cell resolution and its analysis in terms of modularly-composed cellular programs. We demonstrated capabilities of single-cell data portraying to disentangle transcriptional states using intuitive visualization, functional mining, molecular cell stratification, and variability analyses. Our analysis revealed that the T cell composition of the patient’s infusion product as well as the spectrum of their transcriptional states of cells derived from patients with low ICANS grade do not markedly differ from those of cells from high ICANS patients, while the relative abundancies, particularly that of cycling cells, of LAG3-mediated exhaustion and of CAR positive cells, vary. Our study provides molecular details of the transcriptomic landscape with possible impact to overcome neurotoxicity.

Epigenomic interplay in tumor heterogeneity: Potential of epidrugs as adjunct therapy

"Heterogeneity" in tumor mass has immense importance in cancer progression and therapy. The impact of tumor heterogeneity is an emerging field and not yet fully explored. Tumor heterogeneity is mainly considered as intra-tumor heterogeneity and inter-tumor heterogeneity based on their origin. Intra-tumor heterogeneity refers to the discrepancy within the same cancer mass while inter-tumor heterogeneity refers to the discrepancy between different patients having the same tumor type. Both of these heterogeneity types lead to variation in the histopathological as well as clinical properties of the cancer mass which drives disease resistance towards therapeutic approaches. Cancer stem cells (CSCs) act as pinnacle progenitors for heterogeneity development along with various other genetic and epigenetic parameters that are regulating this process. In recent times epigenetic factors are one of the most studied parameters that drive oxidative stress pathways essential during cancer progression. These epigenetic changes are modulated by various epidrugs and have an impact on tumor heterogeneity. The present review summarizes various aspects of epigenetic regulation in the tumor microenvironment, oxidative stress, and progression towards tumor heterogeneity that creates complications during cancer treatment. This review also explores the possible role of epidrugs in regulating tumor heterogeneity and personalized therapy against drug resistance.

Classifying Germinal Center Derived Lymphomas-Navigate a Complex Transcriptional Landscape

Classification of lymphoid neoplasms is based mainly on histologic, immunologic, and (rarer) genetic features. It has been supplemented by gene expression profiling (GEP) in the last decade. Despite the considerable success, particularly in associating lymphoma subtypes with specific transcriptional programs and classifier signatures of up- or downregulated genes, competing molecular classifiers were often proposed in the literature by different groups for the same classification tasks to distinguish, e.g., BL versus DLBCL or different DLBCL subtypes. Moreover, rarer sub-entities such as MYC and BCL2 "double hit lymphomas" (DHL), IRF4-rearranged large cell lymphoma (IRF4-LCL), and Burkitt-like lymphomas with 11q aberration pattern (mnBLL-11q) attracted interest while their relatedness regarding the major classes is still unclear in many respects. We explored the transcriptional landscape of 873 lymphomas referring to a wide spectrum of subtypes by applying self-organizing maps (SOM) machine learning. The landscape reveals a continuum of transcriptional states activated in the different subtypes without clear-cut borderlines between them and preventing their unambiguous classification. These states show striking parallels with single cell gene expression of the active germinal center (GC), which is characterized by the cyclic progression of B-cells. The expression patterns along the GC trajectory are discriminative for distinguishing different lymphoma subtypes. We show that the rare subtypes take intermediate positions between BL, DLBCL, and FL as considered by the 5th edition of the WHO classification of haemato-lymphoid tumors in 2022. Classifier gene signatures extracted from these states as modules of coregulated genes are competitive with literature classifiers. They provide functional-defined classifiers with the option of consenting redundant classifiers from the literature. We discuss alternative classification schemes of different granularity and functional impact as possible avenues toward personalization and improved diagnostics of GC-derived lymphomas.

Integrated Multi-Omics Maps of Lower-Grade Gliomas

Multi-omics high-throughput technologies produce data sets which are not restricted to only one but consist of multiple omics modalities, often as patient-matched tumour specimens. The integrative analysis of these omics modalities is essential to obtain a holistic view on the otherwise fragmented information hidden in this data. We present an intuitive method enabling the combined analysis of multi-omics data based on self-organizing maps machine learning. It "portrays" the expression, methylation and copy number variations (CNV) landscapes of each tumour using the same gene-centred coordinate system. It enables the visual evaluation and direct comparison of the different omics layers on a personalized basis. We applied this combined molecular portrayal to lower grade gliomas, a heterogeneous brain tumour entity. It classifies into a series of molecular subtypes defined by genetic key lesions, which associate with large-scale effects on DNA methylation and gene expression, and in final consequence, drive with cell fate decisions towards oligodendroglioma-, astrocytoma- and glioblastoma-like cancer cell lineages with different prognoses. Consensus modes of concerted changes of expression, methylation and CNV are governed by the degree of co-regulation within and between the omics layers. The method is not restricted to the triple-omics data used here. The similarity landscapes reflect partly independent effects of genetic lesions and DNA methylation with consequences for cancer hallmark characteristics such as proliferation, inflammation and blocked differentiation in a subtype specific fashion. It can be extended to integrate other omics features such as genetic mutation, protein expression data as well as extracting prognostic markers.

A Novel Multi-Omics Analysis Model for Diagnosis and Survival Prediction of Lower-Grade Glioma Patients

Background

Lower-grade gliomas (LGGs) are characterized by remarkable genetic heterogeneity and different clinical outcomes. Classification of LGGs is improved by the development of molecular stratification markers including IDH mutation and 1p/19q chromosomal integrity, which are used as a hallmark of survival and therapy sensitivity of LGG patients. However, the reproducibility and sensitivity of the current classification remain ambiguous. This study aimed to construct more accurate risk-stratification approaches.

Methods

According to bioinformatics, the sequencing profiles of methylation and transcription and imaging data derived from LGG patients were analyzed and developed predictable risk score and radiomics score. Moreover, the performance of predictable models was further validated.

Results

In this study, we determined a cluster of 6 genes that were correlated with IDH mutation/1p19q co-deletion status. Risk score model was calculated based on 6 genes and showed gratifying sensitivity and specificity for survival prediction and therapy response of LGG patients. Furthermore, a radiomics risk score model was established to noninvasively assist judgment of risk score in pre-surgery. Taken together, a predictable nomogram that combined transcriptional signatures and clinical characteristics was established and validated to be preferable to the histopathological classification. Our novel multi-omics nomograms showed a satisfying performance. To establish a user-friendly application, the nomogram was further developed into a web-based platform: https://drw576223193.shinyapps.io/Nomo/, which could be used as a supporting method in addition to the current histopathological-based classification of gliomas.

Conclusions

Our novel multi-omics nomograms showed the satisfying performance of LGG patients and assisted clinicians to draw up individualized clinical management.

Differential expression of microRNAs targeting genes associated with the development of high-grade gliomas

Background

Highly malignant high-grade gliomas are tumors of the central nervous system (CNS). They are solid tumors arising from transformed cells of the brain and/or the spinal cord. In recent years, the expression of genes and regulating miRNAs in glial brain tumors has been actively studied. The present study is devoted to assessing the expression levels of miR-215-5p, miR-22-3p, miR-122-5p, miR-107, miR-324-5p, miR-34a-5p, miR-155. -5p, miR-21-5p, miR-497-5p, miR-330-3p, miR-146a-5p, miR-92a-1-5p, miR-326 and target genes EGFR, SMAD4, SMAD7, SMO, NOTCH1, NOTCH2, HIF1A, EGLIN1/3, KDM1B, KDM1A, MSI1, MSI2, TET1 in high-grade glioma tissues.

Results

As a result of the analysis of the levels of relative expression of the studied genes, there are significant changes (p < 0.05) in tumor tissue for genes: EGFR, SMAD4, SMAD7, SMO, HIF1A, EGLN1/3. We obtained data on a significant change (p < 0.05) in the levels of relative expression for microRNA: hsa-miR-215-5p, hsa-miR-22-3p, hsa-miR-122-5p, hsa-miR-107, hsa-miR-324-5p, hsa-miR-155-5p, hsa-miR-21-5p, hsa-miR-330-3p, hsa-miR-326. Data on the association of overall survival in patients with high-grade glioma and the level of relative expression of the EGFR and HIF1A genes were obtained. The obtained data demonstrate the association of overall survival of patients with high-grade glioma and the level of relative expression of EGFR, HIF1A and hsa-miR-22-3p, hsa-miR-107 and hsa-miR-330-3p.

Conclusions

The obtained data on the expression of genes and microRNAs expand the understanding of the biology of the development of high-grade glial tumors. These data demonstrate new potential therapeutic and prognostic goals in high-grade gliomas.

Treatment Response and Prognosis Evaluation in High-Grade Glioma: An Imaging Review Based on MRI

In recent years, the development of advanced magnetic resonance imaging (MRI) technology and machine learning (ML) have created new tools for evaluating treatment response and prognosis of patients with high-grade gliomas (HGG); however, patient prognosis has not improved significantly. This is mainly due to the heterogeneity between and within HGG tumors, resulting in standard treatment methods not benefitting all patients. Moreover, the survival of patients with HGG is not only related to tumor cells, but also to noncancer cells in the tumor microenvironment (TME). Therefore, during preoperative diagnosis and follow-up treatment of patients with HGG, noninvasive imaging markers are needed to characterize intratumoral heterogeneity, and then to evaluate treatment response and predict prognosis, timeously adjust treatment strategies, and achieve individualized diagnosis and treatment. In this review, we summarize the research progress of conventional MRI, advanced MRI technology, and ML in evaluation of treatment response and prognosis of patients with HGG. We further discuss the significance of the TME in the prognosis of HGG patients, associate imaging features with the TME, indirectly reflecting the heterogeneity within the tumor, and shifting treatment strategies from tumor cells alone to systemic therapy of the TME, which may be a major development direction in the future. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY STAGE: 4.

Grading of adult diffuse gliomas according to the 2021 WHO Classification of Tumors of the Central Nervous System

The grading of gliomas based on histological features has been a subject of debate for several decades. A consensus has not yet been reached because of technical limitations and inter-observer variations. While the traditional grading system has failed to stratify the risk of IDH-mutant astrocytoma, canonical histological and proliferative markers may be applicable to the risk stratification of IDH-wild-type astrocytoma. Numerous studies have examined molecular markers in order to obtain more clinically relevant information that will improve the risk stratification of gliomas. The CDKN2A/B homozygous deletion for IDH-mutant astrocytoma and the following three criteria for IDH-wild-type astrocytoma: the concurrent gain of whole chromosome 7 and loss of whole chromosome 10, TERT promoter mutations, and EGFR amplification, were identified as independent molecular markers of the worst clinical outcomes. Therefore, the 2021 World Health Organization (WHO) Classification of Tumors of the Central Nervous System adopted these molecular markers into the revised grading criteria of IDH-mutant and -wild-type astrocytoma, respectively, as a grading system within tumor types. Of note, several recent studies have shown that some low-grade IDH-wild-type astrocytoma lacking both the molecular glioblastoma signature and genetic alterations typical of pediatric-type gliomas may demonstrate a relatively indolent clinical course, suggesting the existence of lower-grade adult IDH-wild-type astrocytoma. In terms of oligodendroglioma, IDH-mutant, and 1p/19q codeleted, consistent makers that predict poor outcomes have not yet been identified, and, thus, the current criteria have remained unchanged. Molecular testing to fulfill the revised WHO criteria is, however, not always available worldwide, and in that case, an integrated diagnosis combining all available complementary information is highly recommended. This review discusses controversial issues surrounding legacy grading systems and newly identified potential genetic markers of adult diffuse gliomas and provides perspectives on future grading systems.

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⁶-methyladenosine (m⁶A), 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, m⁶A 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, m⁶A status, and stemness characteristics in LGG. Furthermore, these lncRNAs may be potential and alternative therapeutic targets.

Impact of CDKN2A/B Homozygous Deletion on the Prognosis and Biology of IDH-Mutant Glioma

Although hotspot mutations in isocitrate dehydrogenase (IDH) genes are associated with favorable clinical outcomes in glioma, CDKN2A/B homozygous deletion has been identified as an independent predicator of poor prognosis. Accordingly, the 2021 edition of the World Health Organization (WHO) classification of tumors of the central nervous system (CNS) has adopted this molecular feature by upgrading IDH-mutant astrocytoma to CNS WHO grade IV, even in the absence of glioblastoma-specific histological features—necrosis and microvascular proliferation. This new entity of IDH-mutant astrocytoma not only signifies an exception to the generally favorable outcome of IDH-mutant glioma, but also brings into question whether, and, if so, how, CDKN2A/B homozygous deletion overrides the anti-tumor activity of IDH mutation by promoting the proliferation of stem/neural progenitor-like cells. Understanding the mechanism by which IDH mutation requires intact tumor-suppressor genes for conferring favorable outcome may improve therapeutics.

The Transcriptome and Methylome of the Developing and Aging Brain and Their Relations to Gliomas and Psychological Disorders

Mutually linked expression and methylation dynamics in the brain govern genome regulation over the whole lifetime with an impact on cognition, psychological disorders, and cancer. We performed a joint study of gene expression and DNA methylation of brain tissue originating from the human prefrontal cortex of individuals across the lifespan to describe changes in cellular programs and their regulation by epigenetic mechanisms. The analysis considers previous knowledge in terms of functional gene signatures and chromatin states derived from independent studies, aging profiles of a battery of chromatin modifying enzymes, and data of gliomas and neuropsychological disorders for a holistic view on the development and aging of the brain. Expression and methylation changes from babies to elderly adults decompose into different modes associated with the serial activation of (brain) developmental, learning, metabolic and inflammatory functions, where methylation in gene promoters mostly represses transcription. Expression of genes encoding methylome modifying enzymes is very diverse reflecting complex regulations during lifetime which also associates with the marked remodeling of chromatin between permissive and restrictive states. Data of brain cancer and psychotic disorders reveal footprints of pathophysiologies related to brain development and aging. Comparison of aging brains with gliomas supports the view that glioblastoma-like and astrocytoma-like tumors exhibit higher cellular plasticity activated in the developing healthy brain while oligodendrogliomas have a more stable differentiation hierarchy more resembling the aged brain. The balance and specific shifts between volatile and stable and between more irreversible and more plastic epigenomic networks govern the development and aging of healthy and diseased brain.

A Novel Prognostic Tool for Glioma Based on Enhancer RNA-Regulated Immune Genes

Background: Gliomas are the most malignant tumors of the nervous system. Even though their survival outcome is closely affected by immune-related genes (IRGs) in the tumor microenvironment (TME), the corresponding regulatory mechanism remains poorly characterized.

Methods: Specific enhancer RNAs (eRNAs) can be found in tumors, where they control downstream genes. The present study aimed to identify eRNA-regulated IRGs, evaluate their influence on the TME, and use them to construct a novel prognostic model for gliomas.

Results: Thirteen target genes (ADCYAP1R1, BMP2, BMPR1A, CD4, DDX17, ELN, FGF13, MAPT, PDIA2, PSMB8, PTPN6, SEMA6C, and SSTR5) were identified and integrated into a comprehensive risk signature, which distinguished two risk subclasses. Discrepancies between these subclasses were compared to explore potential mechanisms attributed to eRNA-regulated genes, including immune cell infiltration, clinicopathological features, survival outcomes, and chemotherapeutic drug sensitivity. Furthermore, the risk signature was used to construct a prognostic tool that was evaluated by calibration curve, clinical utility, Harrell’s concordance index (0.87; 95% CI: 0.84–0.90), and time-dependent receiver operator characteristic curves (AUCs: 0.93 and 0.89 at 3 and 5 years, respectively). The strong reliability and robustness of the established prognostic tool were validated in another independent cohort. Finally, potential subtypes were explored in patients with grade III tumors.

Conclusion: Overall, eRNAs were associated with immune-related dysfunctions in the TME. Targeting of IRGs regulated by eRNAs could improve immunotherapeutic/therapeutic outcomes.

Metabolic Signature-Based Subtypes May Pave Novel Ways for Low-Grade Glioma Prognosis and Therapy

Metabolic signatures are frequently observed in cancer and are starting to be recognized as important regulators for tumor progression and therapy. Because metabolism genes are involved in tumor initiation and progression, little is known about the metabolic genomic profiles in low-grade glioma (LGG). Here, we applied bioinformatics analysis to determine the metabolic characteristics of patients with LGG from the Cancer Genome Atlas (TCGA) and the Chinese Glioma Genome Atlas (CGGA). We also performed the ConsensusClusterPlus, the CIBERSORT algorithm, the Estimate software, the R package "GSVA," and TIDE to comprehensively describe and compare the characteristic difference between three metabolic subtypes. The R package WGCNA helped us to identify co-expression modules with associated metabolic subtypes. We found that LGG patients were classified into three subtypes based on 113 metabolic characteristics. MC1 patients had poor prognoses and MC3 patients obtained longer survival times. The different metabolic subtypes had different metabolic and immune characteristics, and may have different response patterns to immunotherapy. Based on the metabolic subtype, different patterns were exhibited that reflected the characteristics of each subtype. We also identified eight potential genetic markers associated with the characteristic index of metabolic subtypes. In conclusion, a comprehensive understanding of metabolism associated characteristics and classifications may improve clinical outcomes for LGG.

Distinct DNA Methylation Patterns of Subependymal Giant Cell Astrocytomas in Tuberous Sclerosis Complex

Tuberous sclerosis complex (TSC) is a monogenic disorder caused by mutations in either the TSC1 or TSC2 gene, two key regulators of the mechanistic target of the rapamycin complex pathway. Phenotypically, this leads to growth and formation of hamartomas in several organs, including the brain. Subependymal giant cell astrocytomas (SEGAs) are low-grade brain tumors commonly associated with TSC. Recently, gene expression studies provided evidence that the immune system, the MAPK pathway and extracellular matrix organization play an important role in SEGA development. However, the precise mechanisms behind the gene expression changes in SEGA are still largely unknown, providing a potential role for DNA methylation. We investigated the methylation profile of SEGAs using the Illumina Infinium HumanMethylation450 BeadChip (SEGAs n = 42, periventricular control n = 8). The SEGA methylation profile was enriched for the adaptive immune system, T cell activation, leukocyte mediated immunity, extracellular structure organization and the ERK1 & ERK2 cascade. More interestingly, we identified two subgroups in the SEGA methylation data and show that the differentially expressed genes between the two subgroups are related to the MAPK cascade and adaptive immune response. Overall, this study shows that the immune system, the MAPK pathway and extracellular matrix organization are also affected on DNA methylation level, suggesting that therapeutic intervention on DNA level could be useful for these specific pathways in SEGA. Moreover, we identified two subgroups in SEGA that seem to be driven by changes in the adaptive immune response and MAPK pathway and could potentially hold predictive information on target treatment response.

Transcriptome profile of the sinoatrial ring reveals conserved and novel genetic programs of the zebrafish pacemaker

BACKGROUND

Sinoatrial Node (SAN) is part of the cardiac conduction system, which controls the rhythmic contraction of the vertebrate heart. The SAN consists of a specialized pacemaker cell population that has the potential to generate electrical impulses. Although the SAN pacemaker has been extensively studied in mammalian and teleost models, including the zebrafish, their molecular nature remains inadequately comprehended.

RESULTS

To characterize the molecular profile of the zebrafish sinoatrial ring (SAR) and elucidate the mechanism of pacemaker function, we utilized the transgenic line sqet33mi59BEt to isolate cells of the SAR of developing zebrafish embryos and profiled their transcriptome. Our analyses identified novel candidate genes and well-known conserved signaling pathways involved in pacemaker development. We show that, compared to the rest of the heart, the zebrafish SAR overexpresses several mammalian SAN pacemaker signature genes, which include hcn4 as well as those encoding calcium- and potassium-gated channels. Moreover, genes encoding components of the BMP and Wnt signaling pathways, as well as members of the Tbx family, which have previously been implicated in pacemaker development, were also overexpressed in the SAR. Among SAR-overexpressed genes, 24 had human homologues implicated in 104 different ClinVar phenotype entries related to various forms of congenital heart diseases, which suggest the relevance of our transcriptomics resource to studying human heart conditions. Finally, functional analyses of three SAR-overexpressed genes, pard6a, prom2, and atp1a1a.2, uncovered their novel role in heart development and physiology.

CONCLUSION

Our results established conserved aspects between zebrafish and mammalian pacemaker function and revealed novel factors implicated in maintaining cardiac rhythm. The transcriptome data generated in this study represents a unique and valuable resource for the study of pacemaker function and associated heart diseases.

Electrotherapies for Glioblastoma

Non-thermal, intermediate frequency (100-500 kHz) electrotherapies present a unique therapeutic strategy to treat malignant neoplasms. Here, pulsed electric fields (PEFs) which induce reversible or irreversible electroporation (IRE) and tumour-treating fields (TTFs) are reviewed highlighting the foundations, advances, and considerations of each method when applied to glioblastoma (GBM). Several biological aspects of GBM that contribute to treatment complexity (heterogeneity, recurrence, resistance, and blood-brain barrier(BBB)) and electrophysiological traits which are suggested to promote glioma progression are described. Particularly, the biological responses at the cellular and molecular level to specific parameters of the electrical stimuli are discussed offering ways to compare these parameters despite the lack of a universally adopted physical description. Reviewing the literature, a disconnect is found between electrotherapy techniques and how they target the biological complexities of GBM that make treatment difficult in the first place. An attempt is made to bridge the interdisciplinary gap by mapping biological characteristics to different methods of electrotherapy, suggesting important future research topics and directions in both understanding and treating GBM. To the authors' knowledge, this is the first paper that attempts an in-tandem assessment of the biological effects of different aspects of intermediate frequency electrotherapy methods, thus offering possible strategies toward GBM treatment.

Differential and Prognostic Significance of HOXB7 in Gliomas

Diffuse glioma is the most common primary tumor of the central nervous system. The prognosis of the individual tumor is heavily dependent on its grade and subtype. Homeobox B7 (HOXB7), a member of the homeobox family, is abnormally overexpressed in a variety of tumors. However, its function in glioma is unclear. In this study, HOXB7 mRNA and protein expression levels were analyzed in 401 gliomas from the CGGA RNA-seq database (325 cases) and our hospital (76 cases). HOXB7 expression, at both mRNA and protein levels, were upregulated in glioblastoma (GBM) and isocitrate dehydrogenase 1 (IDH1) wild-type glioma tissues. Kaplan–Meier with log-rank test showed that patients with high HOXB7 expression had a poor prognosis (p < 0.0001). Moreover, HOXB7 protein was deleted in 90.9% (20/22) of oligodendrogliomas and 13.0% (3/23) of astrocytomas. The sensitivity and specificity of HOXB7 protein deletion in oligodendroglioma were 90.9% (20/22) and 87.0% (20/23), respectively. To verify the reliability of using HOXB7 in differentiating oligodendroglioma, we used 1p/19q fluorescence in situ hybridization (FISH) testing as a positive control. The Cohen’s kappa coefficient of HOXB7 immunohistochemistry staining and 1p/19q FISH testing was 0.778 (95% CI: 0.594–0.962, p < 0.001). In conclusion, HOXB7 is an independent predictor of poor prognosis in all grade gliomas. Additionally, HOXB7 is also a highly sensitive and specific indicator to differentiate oligodendroglioma from astrocytoma.

Lower-Grade Gliomas: Predicting DNA Methylation Subtyping and its Consequences on Survival with MR Features

RATIONALE AND OBJECTIVES

To explore associations between MR imaging features, DNA methylation subtyping, and survival in lower-grade gliomas (LGG).

MATERIALS AND METHODS

The MR data from 170 patients generated with the Cancer Imaging Archive were reviewed. The correlation was evaluated by Fisher's Exact Test, Pearson Chi-Square and binary regression analysis. Survival analysis was conducted by using time-dependent ROC analysis and the Kaplan-Meier method (the worst prognosis subgroup).

RESULTS

Identified were 9 (5.3%) M1-subtype, 18 (10.6%) M2-subtype, 48 (28.2%) M3-subtype, 31 (18.2%) M4-subtype and 64 (37.6%) M5-subtype. Patients with M4-subtype had the shortest median OS (49.3 vs. 28.4) months(p < 0.05). The time-dependent ROC for the M4-subtype was 0.83 (95% confidence interval 0.72-0.95) for survival at 12 months, 0.82 (95% confidence interval 0.70-0.94) for survival at 24 months, and 0.74 (95% confidence interval 0.62-0.86) for survival at 36 months. After uni- and multivariate analysis, a nomogram was built based on proportion contrast-enhanced (CE) tumor, extranodular growth, volume_{_cutoff_median}, and location. For the prediction of M4-subtype, the nomogram showed good discrimination, with an area under the curve (AUC) of 0.886 (95% CI: 0.820-952) and was well calibrated. On multivariate logistic regression analysis, volume ≥60cm³ (OR: 0.200; p < 0.001; 95%CI: 0.048-0.834) was associated with M1-subtype (AUC: 0.690). Hemorrhage (OR: 5.443; p = 0.002; 95%CI: 1.844-16.069) and volume > median (OR: 3.256; p = 0.05; 95%CI: 0.992-10.686) were associated with M2-subtype (AUC: 0.733). Proportion CE tumor<=5% (OR: 3.968; P=0.002; 95%CI: 1.634-9.635) was associated with M3-subtype (AUC: 0.632). Poorly-defined (OR: 2.258; p = 0.05; 95%CI: 1.000-5.101) and volume > median (OR: 2.447; p = 0.01; 95%CI: 1.244-4.813) were associated with M5-subtype (AUC: 0.645). Decision curve analysis indicated predictions for all models were clinically useful.

CONCLUSION

This preliminary radiogenomics analysis of lower-grade gliomas demonstrated associations between MR features and DNA methylation subtyping. The shortest survival was observed in patients with M4-subtype. And we have constructed nomogram that enables more accurate predictions of M4-subtype.

EVA1C Is a Potential Prognostic Biomarker and Correlated With Immune Infiltration Levels in WHO Grade II/III Glioma

Background

Immunotherapy is an effective therapeutic approach for multiple human cancer types. However, the correlations between EVA1C and patients’ prognosis as well as immune infiltration remain obscure. Herein, we employed transcriptomic and clinical data extracted from two independent databases to systematically investigate the role of EVA1C in the oncological context.

Methods

The differential expression of EVA1C was analyzed via TCGA and Oncomine databases. We evaluated the influence of EVA1C on clinical prognosis using Kaplan-Meier plotter. We then used the expression profiler to calculate stromal score, immune score, and ESTIMATE score based on the ESTIMATE algorithm. The abundance of infiltrating immune cells was calculated via TIMER. The correlations between EVA1C expression and immune infiltration levels were analyzed in two independent cohorts.

Results

In patients with World Health Organization (WHO) grade II/III glioma, high EVA1C expression was associated with malignant clinicopathological features and poor overall survival in both cohorts. EVA1C expression was positively associated with immune infiltration levels of B cell, CD4+ T cell, neutrophil, macrophage, and dendritic cells (DCs). Besides, EVA1C expression strongly correlated with diverse immune marker sets. And the predictive power of EVA1C was better than that of other indicators in predicting high immune infiltration levels in glioma.

Conclusions

For the first time, we identified the overexpression of EVA1C in glioma, which was tightly correlated with the high infiltration levels of multiple immune cells as well as poor prognosis. Meanwhile, EVA1C might be a potential biomarker for predicting high immune infiltration in WHO grade II/III gliomas.

High-Resolution Cartography of the Transcriptome and Methylome Landscapes of Diffuse Gliomas

Molecular mechanisms of lower-grade (II-III) diffuse gliomas (LGG) are still poorly understood, mainly because of their heterogeneity. They split into astrocytoma- (IDH-A) and oligodendroglioma-like (IDH-O) tumors both carrying mutations(s) at the isocitrate dehydrogenase (IDH) gene and into IDH wild type (IDH-wt) gliomas of glioblastoma resemblance. We generated detailed maps of the transcriptomes and DNA methylomes, revealing that cell functions divided into three major archetypic hallmarks: (i) increased proliferation in IDH-wt and, to a lesser degree, IDH-O; (ii) increased inflammation in IDH-A and IDH-wt; and (iii) the loss of synaptic transmission in all subtypes. Immunogenic properties of IDH-A are diverse, partly resembling signatures observed in grade IV mesenchymal glioblastomas or in grade I pilocytic astrocytomas. We analyzed details of coregulation between gene expression and DNA methylation and of the immunogenic micro-environment presumably driving tumor development and treatment resistance. Our transcriptome and methylome maps support personalized, case-by-case views to decipher the heterogeneity of glioma states in terms of data portraits. Thereby, molecular cartography provides a graphical coordinate system that links gene-level information with glioma subtypes, their phenotypes, and clinical context.

Online database for brain cancer-implicated genes: exploring the subtype-specific mechanisms of brain cancer

Background

Brain cancer is one of the eight most common cancers occurring in people aged 40+ and is the fifth-leading cause of cancer-related deaths for males aged 40–59. Accurate subtype identification is crucial for precise therapeutic treatment, which largely depends on understanding the biological pathways and regulatory mechanisms associated with different brain cancer subtypes. Unfortunately, the subtype-implicated genes that have been identified are scattered in thousands of published studies. So, systematic literature curation and cross-validation could provide a solid base for comparative genetic studies about major subtypes.

Results

Here, we constructed a literature-based brain cancer gene database (BCGene). In the current release, we have a collection of 1421 unique human genes gathered through an extensive manual examination of over 6000 PubMed abstracts. We comprehensively annotated those curated genes to facilitate biological pathway identification, cancer genomic comparison, and differential expression analysis in various anatomical brain regions. By curating cancer subtypes from the literature, our database provides a basis for exploring the common and unique genetic mechanisms among 40 brain cancer subtypes. By further prioritizing the relative importance of those curated genes in the development of brain cancer, we identified 33 top-ranked genes with evidence mentioned only once in the literature, which were significantly associated with survival rates in a combined dataset of 2997 brain cancer cases.

Conclusion

BCGene provides a useful tool for exploring the genetic mechanisms of and gene priorities in brain cancer. BCGene is freely available to academic users at http://soft.bioinfo-minzhao.org/bcgene/.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07793-x.

Integrative Analysis of DNA Methylation and Transcriptome Identifies a Predictive Epigenetic Signature Associated With Immune Infiltration in Gliomas

Glioma is the most common primary brain tumor with poor prognosis and high mortality. The purpose of this study was to use the epigenetic signature to predict prognosis and evaluate the degree of immune infiltration in gliomas. We integrated gene expression profiles and DNA methylation data of lower-grade glioma and glioblastoma to explore epigenetic differences and associated differences in biological function. Cox regression and lasso analysis were used to develop an epigenetic signature based on eight DNA methylation sites to predict prognosis of glioma patients. Kaplan–Meier analysis showed that the overall survival time of high- and low-risk groups was significantly separated, and ROC analysis verified that the model had great predictive ability. In addition, we constructed a nomogram based on age, sex, 1p/19q status, glioma type, and risk score. The epigenetic signature was obviously associated with tumor purity, immune checkpoints, and tumor-immune infiltrating cells (CD8+ T cells, gamma delta T cells, M0 macrophages, M1 macrophages, M2 macrophages, activated NK cells, monocytes, and activated mast cells) and thus, it may find application as a guide for the evaluation of immune infiltration or in treatment decisions in immunotherapy.

Deciphering the Transcriptomic Heterogeneity of Duodenal Coeliac Disease Biopsies

Coeliac disease (CD) is a clinically heterogeneous autoimmune disease with variable presentation and progression triggered by gluten intake. Molecular or genetic factors contribute to disease heterogeneity, but the reasons for different outcomes are poorly understood. Transcriptome studies of tissue biopsies from CD patients are scarce. Here, we present a high-resolution analysis of the transcriptomes extracted from duodenal biopsies of 24 children and adolescents with active CD and 21 individuals without CD but with intestinal afflictions as controls. The transcriptomes of CD patients divide into three groups-a mixed group presenting the control cases, and CD-low and CD-high groups referring to lower and higher levels of CD severity. Persistence of symptoms was weakly associated with subgroup, but the highest marsh stages were present in subgroup CD-high, together with the highest cell cycle rates as an indicator of virtually complete villous atrophy. Considerable variation in inflammation-level between subgroups was further deciphered into immune cell types using cell type de-convolution. Self-organizing maps portrayal was applied to provide high-resolution landscapes of the CD-transcriptome. We find asymmetric patterns of miRNA and long non-coding RNA and discuss the effect of epigenetic regulation. Expression of genes involved in interferon gamma signaling represent suitable markers to distinguish CD from non-CD cases. Multiple pathways overlay in CD biopsies in different ways, giving rise to heterogeneous transcriptional patterns, which potentially provide information about etiology and the course of the disease.

A comprehensive model including preoperative peripheral blood inflammatory markers for prediction of the prognosis of diffuse spinal cord astrocytoma following surgery

PURPOSE

Due to the rarity of diffuse spinal cord astrocytoma, an effective model is still lacking to stratify their prognosis. Here, we aimed to establish a prognostic model through comprehensively evaluating clinicopathological features and preoperative peripheral blood inflammatory markers in 89 cases.

METHODS

We performed univariate and multivariate Cox regression to identify prognosis factors. The Kaplan-Meier curves and ROC curves were employed to compare the prognostic value of selected factors.

RESULTS

In addition to clinicopathological factors, we revealed the preoperative peripheral blood leukocyte count, neutrophils-to-lymphocytes ratio (NLR), and platelet-to-lymphocyte ratio (PLR) were also significantly correlated with overall survival of spinal cord astrocytoma in univariate Cox regression, and NLR was still significant in multivariate Cox analysis. Further, we demonstrated that NLR ≤ 3.65 and preoperative McCormick score (MMS) ≤ 3 were independently correlated with better survival of WHO grade IV tumors. Meanwhile, Ki-67 < 10% and resection extent ≥ 90% were independent prognostic factors in WHO grade II/III tumors. Finally, we developed a prognostic model that had better predictive efficiencies than WHO grade and histological grade for 1-year (AUC = 76.6), 2- year (AUC = 80.9), and 3-year (AUC = 80.3) survival. This model could classify tumors into 4 classifications with increasingly poor prognosis: 1, WHO grade II/III, with Ki-67 < 10% and resection extent ≥ 90%; 2, WHO grade II/III, Ki-67 ≥ 10% or resection < 90%; 3, WHO grade IV, NLR ≤ 3.65 and MMS ≤ 3; 4, WHO grade IV, with NRL > 3.65 or MMS = 4.

CONCLUSION

We successfully constructed a comprehensive prognostic model including preoperative peripheral blood inflammatory markers, which can stratify diffuse spinal cord astrocytoma into 4 subgroups.

DZIP3 is a key factor to stratify IDH1 wild-type lower-grade gliomas

BACKGROUND

Malignant glioma is the most common form of primary malignant brain cancer. Heterogeneity is the hallmark of glioma. DAZ-interacting zinc finger 3 (DZIP3), acts as an RNA-binding RING-type ubiquitin ligase; however, its function in glioma is yet unclear.

RESULTS

The DZIP3 expression was related to the World Health Organization (WHO) grade and isocitrate dehydrogenase 1(IDH1) status, as well as the clinical outcome. Malignant cases exhibit lower DZIP3 expression. DZIP3 was an independent predictive factor of good prognosis in all grade and lower grade gliomas (p < 0.0001). Gene enrichment analysis and immunohistochemistry indicated that DZIP3 affected the biological behavior of glioma through the angiogenesis pathway. Moreover, based on DZIP3 expression, IDH1 wild-type lower-grade gliomas could be divided into two groups with different survival time.

CONCLUSION

In conclusion, the loss of DZIP3 may be involved in the mechanism of angiogenesis in the invasive biological process of glioma. These findings laid an understanding of DZIP3-specific clinical features in glioma.

METHODS

A total of 325 glioma patients from the Chinese Glioma Genome Atlas (CGGA) RNA-seq cohort comprised the training cohort, while 265 patients from the GSE 16011 array cohort formed the validation cohort. The mRNA expression of DZIP3 and clinical characteristics was assessed. DZIP3 protein expression and microvessel density (MVD) were evaluated by immunohistochemistry (IHC).

The Human Blood Transcriptome in a Large Population Cohort and Its Relation to Aging and Health

Background: The blood transcriptome is expected to provide a detailed picture of an organism's physiological state with potential outcomes for applications in medical diagnostics and molecular and epidemiological research. We here present the analysis of blood specimens of 3,388 adult individuals, together with phenotype characteristics such as disease history, medication status, lifestyle factors, and body mass index (BMI). The size and heterogeneity of this data challenges analytics in terms of dimension reduction, knowledge mining, feature extraction, and data integration.

Methods: Self-organizing maps (SOM)-machine learning was applied to study transcriptional states on a population-wide scale. This method permits a detailed description and visualization of the molecular heterogeneity of transcriptomes and of their association with different phenotypic features.

Results: The diversity of transcriptomes is described by personalized SOM-portraits, which specify the samples in terms of modules of co-expressed genes of different functional context. We identified two major blood transcriptome types where type 1 was found more in men, the elderly, and overweight people and it upregulated genes associated with inflammation and increased heme metabolism, while type 2 was predominantly found in women, younger, and normal weight participants and it was associated with activated immune responses, transcriptional, ribosomal, mitochondrial, and telomere-maintenance cell-functions. We find a striking overlap of signatures shared by multiple diseases, aging, and obesity driven by an underlying common pattern, which was associated with the immune response and the increase of inflammatory processes.

Conclusions: Machine learning applications for large and heterogeneous omics data provide a holistic view on the diversity of the human blood transcriptome. It provides a tool for comparative analyses of transcriptional signatures and of associated phenotypes in population studies and medical applications.

oposSOM-Browser: an interactive tool to explore omics data landscapes in health science

BACKGROUND

oposSOM is a comprehensive, machine learning based open-source data analysis software combining functionalities such as diversity analyses, biomarker selection, function mining, and visualization.

RESULTS

These functionalities are now available as interactive web-browser application for a broader user audience interested in extracting detailed information from high-throughput omics data sets pre-processed by oposSOM. It enables interactive browsing of single-gene and gene set profiles, of molecular 'portrait landscapes', of associated phenotype diversity, and signalling pathway activation patterns.

CONCLUSION

The oposSOM-Browser makes available interactive data browsing for five transcriptome data sets of cancer (melanomas, B-cell lymphomas, gliomas) and of peripheral blood (sepsis and healthy individuals) at www.izbi.uni-leipzig.de/opossom-browser .

Development of MR-based preoperative nomograms predicting DNA copy number subtype in lower grade gliomas with prognostic implication

OBJECTIVES

We aimed to determine the value of MR-based preoperative nomograms in predicting DNA copy number (CN) subtype in lower grade glioma (LGG) patients.

METHODS

The overall survival (OS) data were analyzed. MRI data of 170 subjects were retrospectively analyzed. The correlation was explored by univariate and multivariate regression analysis.

RESULTS

CN2 subtype was associated with shortest median OS (CN2 subtype vs. others: 46.8 vs. 221.7 months, p < 0.05). The time-dependent receiver operating characteristic for the CN2 subtype was 0.80 (95% CI: 0.74-0.85) for survival at 1 year, 0.80 (95% CI: 0.75-0.85) for survival at 2 years, and 0.77 (95% CI: 0.73-0.83) for survival at 3 years. On multivariate analysis, hemorrhage (OR: 0.118; p < 0.001; 95% CI: 0.037-0.376), poorly defined margin (OR: 4.592; p < 0.001; 95% CI: 1.965-10.730), extranodular growth (OR: 0.247; p = 0.006; 95% CI: 0.091-0.671), and volume ≥ 60 cm³ (OR: 4.734.256; p < 0.001; 95% CI: 2.051-10.924) were associated with CN1 subtype (AUC: 0.781). Proportion CE tumor (OR: 5.905; p = 0.007; 95% CI: 1.622-21.493), extranodular growth (OR: 9.047; p = 0.001; 95% CI: 2.349-34.846), width ≥ median (OR: 0.231; p = 0.049; 95% CI: 0.054-0.998), and depth ≥ median (OR: 0.192; p = 0.023; 95% CI: 0.046-0.799) were associated with CN2 subtype (AUC: 0.854). Necrosis/cystic (OR: 6.128; p = 0.007; 95% CI: 1.635-22.968), hemorrhage (OR: 5.752; p = 0.002; 95% CI: 1.953-16.942), poorly defined margin (OR: 0.164; p < 0.001; 95% CI: 0.063-0.427), and volume ≥ median (OR: 4.422; p < 0.001; 95% CI: 1.925-10.160) were associated with CN3 subtype (AUC: 0.808). All three nomograms showed good discrimination and calibration. Decision curve analysis supported that all nomograms were clinically useful. The average accuracy of the tenfold cross-validation was 0.680 (CN1), 0.794 (CN2), and 0.894 (CN3), respectively.

CONCLUSIONS

The shortest OS was observed in patients with CN2 subtype. This preliminary radiogenomics analysis revealed that the MR-based preoperative nomograms provide individualized prediction of DNA copy number subtype in LGG patients.

KEY POINTS

• This preliminary radiogenomics analysis of LGG revealed that the MR-based preoperative nomograms provide individualized prediction of DNA copy number subtype in LGG patients. • The AUC for the ROC curve was 0.781 for CN1 subtype, 0.854 for CN2 subtype, and 0.808 for CN3 subtype. Decision curve analysis supported that all nomograms were clinically useful. • The sensitivity was 0.779 (CN1), 0.731 (CN2), and 0.851 (CN3), respectively. The specificity was 0.664 (CN1), 0.872 (CN2), and 0.625 (CN3), respectively. And the accuracy was 0.717 (CN1), 0.849 (CN2), and 0.692 (CN3), respectively.

Zinc as a plausible epigenetic modulator of glioblastoma multiforme

Glioblastoma Multiforme (GBM) is an aggressive brain tumor (WHO grade 4 astrocytoma) with unknown causes and is associated with a reduced life expectancy. The available treatment options namely radiotherapy, surgery and chemotherapy have failed to improve life expectancy. Out of the various therapeutic approaches, epigenetic therapy is one of the most studied. Epigenetic therapy is involved in the effective treatment of GBM by inhibiting DNA methyltransferase, histone deacetylation and non-coding RNA. It also promotes the expression of the tumor suppressor gene and is involved in the suppression of the oncogene. Various targets are being studied to implement proper epigenetic regulation to control GBM effectively. Zinc is one of the micronutrients which is considered to maintain epigenetic regulation by promoting the proper DNA folding, protecting genetic material from the oxidative damage and controlling the enzyme activation involved in the epigenetic regulation. Here, we are discussing the importance of zinc in regulating the epigenetic modifications and assessing its role in glioblastoma research. The discussion also highlights the importance of artificial intelligence using epigenetics for envisaging the glioma progression, diagnosis and its management.

Identification of purity and prognosis-related gene signature by network analysis and survival analysis in brain lower grade glioma

Tumour microenvironment of brain lower grade glioma (LGG) consists of non‐tumour cells including stromal cells and immune cells mainly. These non‐tumour cells dilute the purity of LGG and play pivotal roles in tumour growth and development, thereby affecting patient prognosis. Tumour purity is also associated with molecular subtypes of LGG. In this study, we discovered the most relevant module to purity by weighted gene co‐expression network analysis (WGCNA) and afterwards performed consensus network analysis and survival analysis to filter 61 significant genes related to both purity and prognosis. In turn, we built a simplified model based on the calculation of purity score, and consensus measurement of purity estimation (CPE), with a satisfactory predictive performance by random forest regression. HLA‐E, MSN, GNG‐5, MYL12A, ITGB4, PDPN, AGTRAP, S100A4, PLSCR1, VAMP5 were selected as the most relevant genes correlating to both purity and prognosis. The risk score model based on the 10 genes could moderately predict patients’ overall survival. These 10 genes, respectively, were positively correlated positively to immunosuppressive cells like macrophage M2, but negatively correlated to patient prognosis, which may explain partially the poor prognosis with low‐purity group.

Identification of Prognostic Signatures of Alternative Splicing in Glioma

Alternative splicing (AS) is a ubiquitous mechanism in which pre-mRNA can be spliced into divergent variants and involved in carcinogenesis and progression in several cancers. In the present study, we systematically profiled prognostic AS signatures involving both low grade glioma (LGG) and glioblastoma (GBM) and investigated the association of AS signatures with tumor grade and IDH1 status in glioma. Percent spliced in (PSI) values and corresponding clinical data were obtained from TCGA SpliceSeq and TCGA data portal, respectively. Prognostic AS signatures were identified using univariate and stepwise multivariate Cox regression. Heatmap analysis was performed based on prognostic AS signatures. A prognostic signature was established with 69 and 88 AS events, including specific splicing events of MUTYH, STEAP3, and CTNNB1, in LGG and GBM cohorts, respectively. The area under the curve (AUC) of the prediction model was 0.968 at 2000 days of overall survival (OS) in the LGG cohort and 0.966 at 450 days of OS in the GBM cohort. In addition, these prognostic AS signatures could complement current molecular classification, such as IDH1 mutation, 1p/19q codeletion, and ATRX loss, of glioma and further identify potential subgroups of glioma with the same molecular features. In conclusion, our study systematically profiled prognostic AS events involving both low grade glioma and glioblastoma for the first time, which also shed light on the crosstalk between AS signatures and molecular features of glioma.

The epigenetics of brain tumors and its modulation during radiation: A review

The brain tumor is the abnormal growth of heterogeneous cells around the central nervous system and spinal cord. Most clinically prominent brain tumors affecting both adult and pediatric are glioblastoma, medulloblastoma, and ependymoma and they are classified according to their origin of tissue. Chemotherapy, radiotherapy, and surgery are important treatments available to date. However, these treatments fail due to multiple reasons, including chemoresistance and radiation resistance of cancer cells. Thus, there is a need of new therapeutic designs to target cell signaling and molecular events which are responsible for this resistance. Recently epigenetic changes received increased attention because it helps in understanding chromatin-mediated disease mechanism. The epigenetic modification alters chromatin structure that affects the docking site of many drugs which cause chemo-resistance of cancer therapy. This review centers the mechanism of how epigenetic changes affect the transcription repression and activation of various genes including Polycomb gene, V-Myc avian myelocytomatosis viral oncogene (MYCN). This review also put forth the pathway of radiation-induced reactive oxygen species generation and its role in epigenetic changes in the cellular level and its impact on tissue physiology. Additionally, there is a strong relationship between the behavior of an individual and environment-induced epigenetic regulation of gene expression. The review also discusses Transcriptome heterogeneity and role of tumor microenvironment in glioblastoma. Overall, this review emphasis important and novel epigenetic targets that could be of therapeutic benefit, which helps in overcoming the unsolved chromatin alteration in brain cancer.

Prognostic factor identification by analysis of the gene expression and DNA methylation data in glioma

Objective: This study was aimed to identify prognostic factors in glioma by analysis of the gene expression and DNA methylation data. Methods: The RNAseq and DNA methylation data associated with glioma were downloaded from GEO and TCGA databases to analyze the differentially expressed genes (DEGs) and methylated genes between tumor and normal tissues. Function and pathway analyses, co-expression network and survival analysis were performed based on these DEGs. The intersection genes of DEGs and differentially methylated genes were obtained followed by function analysis. Results: Total 2190 DEGs were identified between tumor and normal tissues, which were significantly enriched in neuron differentiation associated functions, as well as ribosome pathway. There were 6186 methylation sites (2834 up-regulated and 3352 down-regulated) with significant differences in tumor vs. normal. In the constructed co-expression network, DPP6, MAPK10 and RPL3 were hub genes. Survival analysis of 20 DEGs obtained 18 prognostic genes, among which 9 were differentially methylated, such as LHFPL tetraspan subfamily member 3 (LHFPL3), cadherin 20 (CDH20), complexin 2 (CPLX2), and tenascin R (TNR). The intersection of DEGs and differentially methylated genes (632 genes) were significantly enriched in functions of neuron differentiation. Conclusion: DPP6, MAPK10 and RPL3 may play important roles in tumorigenesis of glioma. Additionally, methylation of LHFPL3, CDH20, CPLX2, and TNR may serve as prognostic factors of glioma.

Radiomic profiles in diffuse glioma reveal distinct subtypes with prognostic value

PURPOSE

To evaluate a radiomic approach for the stratification of diffuse gliomas with distinct prognosis and provide additional resolution of their clinicopathological and molecular characteristics.

METHODS

For this retrospective study, a total of 704 radiomic features were extracted from the multi-channel MRI data of 166 diffuse gliomas. Survival-associated radiomic features were identified and submitted to distinguish glioma subtypes using consensus clustering. Multi-layered molecular data were used to observe the different clinical and molecular characteristics between radiomic subtypes. The relative profiles of an array of immune cell infiltrations were measured gene set variation analysis approach to explore differences in tumor immune microenvironment.

RESULTS

A total of 6 categories, including 318 radiomic features were significantly correlated with the overall survival of glioma patients. Two subgroups with distinct prognosis were separated by consensus clustering of radiomic features that significantly associated with survival. Histological stage and molecular factors, including IDH status and MGMT promoter methylation status were significant differences between the two subtypes. Furthermore, gene functional enrichment analysis and immune infiltration pattern analysis also hinted that the inferior prognosis subtype may more response to immunotherapy.

CONCLUSION

A radiomic model derived from multi-parameter MRI of the gliomas was successful in the risk stratification of diffuse glioma patients. These data suggested that radiomics provided an alternative approach for survival estimation and may improve clinical decision-making.

Comprehensive profiling identifies a novel signature with robust predictive value and reveals the potential drug resistance mechanism in glioma

Background

Gliomas are the most common and malignant brain tumors. The standard therapy is surgery combined with radiotherapy, chemotherapy, and/or other comprehensive methods. However, the emergence of chemoresistance is the main obstacle in treatment and its mechanism is still unclear.

Methods

We firstly developed a multi-gene signature by integrated analysis of cancer stem cell and drug resistance related genes. The Chinese Glioma Genome Atlas (CGGA, 325 samples) and The Cancer Genome Atlas (TCGA, 699 samples) datasets were then employed to verify the efficacy of the risk signature and investigate its significance in glioma prognosis. GraphPad Prism, SPSS and R language were used for statistical analysis and graphical work.

Results

This signature could distinguish the prognosis of patients, and patients with high risk score exhibited short survival time. The Cox regression and Nomogram model indicated the independent prognostic performance and high prognostic accuracy of the signature for survival. Combined with a well-known chemotherapy impact factor-MGMT promoter methylation status, this risk signature could further subdivide patients with distinct survival. Functional analysis of associated genes revealed signature-related biological process of cell proliferation, immune response and cell stemness. These mechanisms were confirmed in patient samples.

Conclusions

The signature was an independent and powerful prognostic biomarker in glioma, which would improve risk stratification and provide a more accurate assessment of personalized treatment.

Epigenetic Reprogramming for Targeting IDH-Mutant Malignant Gliomas

Targeting the epigenome has been considered a compelling treatment modality for several cancers, including gliomas. Nearly 80% of the lower-grade gliomas and secondary glioblastomas harbor recurrent mutations in isocitrate dehydrogenase (IDH). Mutant IDH generates high levels of 2-hydroxyglutarate (2-HG) that inhibit various components of the epigenetic machinery, including histone and DNA demethylases. The encouraging results from current epigenetic therapies in hematological malignancies have reinvigorated the interest in solid tumors and gliomas, both preclinically and clinically. Here, we summarize the recent advancements in epigenetic therapy for lower-grade gliomas and discuss the challenges associated with current treatment options. A particular focus is placed on therapeutic mechanisms underlying favorable outcome with epigenetic-based drugs in basic and translational research of gliomas. This review also highlights emerging bridges to combination treatment with respect to epigenetic drugs. Given that epigenetic therapies, particularly DNA methylation inhibitors, increase tumor immunogenicity and antitumor immune responses, appropriate drug combinations with immune checkpoint inhibitors may lead to improvement of treatment effectiveness of immunotherapy, ultimately leading to tumor cell eradication.