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

Abstract Background Diffuse lower WHO grade II and III gliomas (LGG) are slowly progressing brain tumors, many of which eventually transform into a more aggressive type. LGG is characterized by widespread genetic and transcriptional heterogeneity, yet little is known about the heterogeneity of the DNA methylome, its function in tumor biology, coupling with the transcriptome and tumor microenvironment and its possible impact for tumor development. Methods We here present novel DNA methylation data of an LGG-cohort collected in the German Glioma Network containing about 85% isocitrate dehydrogenase (IDH) mutated tumors and performed a combined bioinformatics analysis using patient-matched genome and transcriptome data. Results Stratification of LGG based on gene expression and DNA-methylation provided four consensus subtypes. We characterized them in terms of genetic alterations, functional context, cellular composition, tumor microenvironment and their possible impact for treatment resistance and prognosis. Glioma with astrocytoma-resembling phenotypes constitute the largest fraction of nearly 60%. They revealed largest diversity and were divided into four expression and three methylation groups which only partly match each other thus reflecting largely decoupled expression and methylation patterns. We identified a novel G-protein coupled receptor and a cancer-related ‘keratinization’ methylation signature in in addition to the glioma-CpG island methylator phenotype (G-CIMP) signature. These different signatures overlap and combine in various ways giving rise to diverse methylation and expression patterns that shape the glioma phenotypes. The decrease of global methylation in astrocytoma-like LGG associates with higher WHO grade, age at diagnosis and inferior prognosis. We found analogies between astrocytoma-like LGG with grade IV IDH-wild type tumors regarding possible worsening of treatment resistance along a proneural-to-mesenchymal axis. Using gene signature-based inference we elucidated the impact of cellular composition of the tumors including immune cell bystanders such as macrophages. Conclusions Genomic, epigenomic and transcriptomic factors act in concert but partly also in a decoupled fashion what underpins the need for integrative, multidimensional stratification of LGG by combining these data on gene and cellular levels to delineate mechanisms of gene (de-)regulation and to enable better patient stratification and individualization of treatment