Epigenetic landscape reorganization and reactivation of embryonic development genes are associated with malignancy in IDH-mutant astrocytoma.

Accurate grading of IDH-mutant gliomas defines patient prognosis and guides the treatment path. Histological grading is however difficult and, apart from CDKN2A/B homozygous deletions in IDH-mutant astrocytomas, there are no other objective molecular markers used for grading. Experimental Design: RNA-sequencing was conducted on primary IDH-mutant astrocytomas (n=138) included in the prospective CATNON trial, which was performed to assess the prognostic effect of adjuvant and concurrent temozolomide. We integrated the RNA sequencing data with matched DNA-methylation and NGS data. We also used multi-omics data from IDH-mutant astrocytomas included in the TCGA dataset and validated results on matched primary and recurrent samples from the GLASS-NL study. We used the DNA-methylation profiles to generate a Continuous Grading Coefficient (CGC) that is based on classification scores derived from a CNS-tumor classifier. We found that the CGC was an independent predictor of survival outperforming current WHO-CNS5 and methylation-based classification. Our RNA-sequencing analysis revealed four distinct transcription clusters that were associated with i) an upregulation of cell cycling genes; ii) a downregulation of glial differentiation genes; iii) an upregulation of embryonic development genes (e.g. HOX, PAX and TBX) and iv) an upregulation of extracellular matrix genes. The upregulation of embryonic development genes was associated with a specific increase of CpG island methylation near these genes.

Transcriptome analysis reveals tumor microenvironment changes in glioblastoma.

A better understanding of transcriptional evolution of IDH-wild-type glioblastoma may be crucial for treatment optimization. Here, we perform RNA sequencing (RNA-seq) (n = 322 test, n = 245 validation) on paired primary-recurrent glioblastoma resections of patients treated with the current standard of care. Transcriptional subtypes form an interconnected continuum in a two-dimensional space. Recurrent tumors show preferential mesenchymal progression. Over time, hallmark glioblastoma genes are not significantly altered. Instead, tumor purity decreases over time and is accompanied by co-increases in neuron and oligodendrocyte marker genes and, independently, tumor-associated macrophages. A decrease is observed in endothelial marker genes. These composition changes are confirmed by single-cell RNA-seq and immunohistochemistry. An extracellular matrix-associated gene set increases at recurrence and bulk, single-cell RNA, and immunohistochemistry indicate it is expressed mainly by pericytes. This signature is associated with significantly worse survival at recurrence. Our data demonstrate that glioblastomas evolve mainly by microenvironment (re-)organization rather than molecular evolution of tumor cells.

The EGFRvIII transcriptome in glioblastoma: A meta-omics analysis.

BACKGROUND: EGFR is among the genes most frequently altered in glioblastoma, with exons 2-7 deletions (EGFRvIII) being among its most common genomic mutations. There are conflicting reports about its prognostic role and it remains unclear whether and how it differs in signaling compared with wildtype EGFR. METHODS: To better understand the oncogenic role of EGFRvIII, we leveraged 4 large datasets into 1 large glioblastoma transcriptome dataset (n = 741) alongside 81 whole-genome samples from 2 datasets. RESULTS: The EGFRvIII/EGFR expression ratios differ strongly between tumors and range from 1% to 95%. Interestingly, the slope of relative EGFRvIII expression is near-linear, which argues against a more positive selection pressure than EGFR wildtype. An absence of selection pressure is also suggested by the similar survival between EGFRvIII-positive and -negative glioblastoma patients. EGFRvIII levels are inversely correlated with pan-EGFR (all wildtype and mutant variants) expression, which indicates that EGFRvIII has a higher potency in downstream pathway activation. EGFRvIII-positive glioblastomas have a lower CDK4 or MDM2 amplification incidence than EGFRvIII-negative (P = .007), which may point toward crosstalk between these pathways. EGFRvIII-expressing tumors have an upregulation of "classical" subtype genes compared to those with EGFR-amplification only (P = 3.873e-6). Genomic breakpoints of the EGFRvIII deletions have a preference toward the 3'-end of the large intron-1. These preferred breakpoints preserve a cryptic exon resulting in a novel EGFRvIII variant and preserve an intronic enhancer. CONCLUSIONS: These data provide deeper insights into the complex EGFRvIII biology and provide new insights for targeting EGFRvIII mutated tumors.