DNA methylation-based classification of central nervous system tumours

Accurate pathological diagnosis is crucial for optimal management of patients with cancer. For the approximately 100 known tumour types of the central nervous system, standardization of the diagnostic process has been shown to be particularly challenging-with substantial inter-observer variability in the histopathological diagnosis of many tumour types. Here we present a comprehensive approach for the DNA methylation-based classification of central nervous system tumours across all entities and age groups, and demonstrate its application in a routine diagnostic setting. We show that the availability of this method may have a substantial impact on diagnostic precision compared to standard methods, resulting in a change of diagnosis in up to 12% of prospective cases. For broader accessibility, we have designed a free online classifier tool, the use of which does not require any additional onsite data processing. Our results provide a blueprint for the generation of machine-learning-based tumour classifiers across other cancer entities, with the potential to fundamentally transform tumour pathology.

DNA methylation-based classification and grading system for meningioma: a multicentre, retrospective analysis

BACKGROUND

The WHO classification of brain tumours describes 15 subtypes of meningioma. Nine of these subtypes are allotted to WHO grade I, and three each to grade II and grade III. Grading is based solely on histology, with an absence of molecular markers. Although the existing classification and grading approach is of prognostic value, it harbours shortcomings such as ill-defined parameters for subtypes and grading criteria prone to arbitrary judgment. In this study, we aimed for a comprehensive characterisation of the entire molecular genetic landscape of meningioma to identify biologically and clinically relevant subgroups.

METHODS

In this multicentre, retrospective analysis, we investigated genome-wide DNA methylation patterns of meningiomas from ten European academic neuro-oncology centres to identify distinct methylation classes of meningiomas. The methylation classes were further characterised by DNA copy number analysis, mutational profiling, and RNA sequencing. Methylation classes were analysed for progression-free survival outcomes by the Kaplan-Meier method. The DNA methylation-based and WHO classification schema were compared using the Brier prediction score, analysed in an independent cohort with WHO grading, progression-free survival, and disease-specific survival data available, collected at the Medical University Vienna (Vienna, Austria), assessing methylation patterns with an alternative methylation chip.

FINDINGS

We retrospectively collected 497 meningiomas along with 309 samples of other extra-axial skull tumours that might histologically mimic meningioma variants. Unsupervised clustering of DNA methylation data clearly segregated all meningiomas from other skull tumours. We generated genome-wide DNA methylation profiles from all 497 meningioma samples. DNA methylation profiling distinguished six distinct clinically relevant methylation classes associated with typical mutational, cytogenetic, and gene expression patterns. Compared with WHO grading, classification by individual and combined methylation classes more accurately identifies patients at high risk of disease progression in tumours with WHO grade I histology, and patients at lower risk of recurrence among WHO grade II tumours (p=0·0096) from the Brier prediction test). We validated this finding in our independent cohort of 140 patients with meningioma.

INTERPRETATION

DNA methylation-based meningioma classification captures clinically more homogenous groups and has a higher power for predicting tumour recurrence and prognosis than the WHO classification. The approach presented here is potentially very useful for stratifying meningioma patients to observation-only or adjuvant treatment groups. We consider methylation-based tumour classification highly relevant for the future diagnosis and treatment of meningioma.

FUNDING

German Cancer Aid, Else Kröner-Fresenius Foundation, and DKFZ/Heidelberg Institute of Personalized Oncology/Precision Oncology Program.

Novel, improved grading system(s) for IDH-mutant astrocytic gliomas

According to the 2016 World Health Organization Classification of Tumors of the Central Nervous System (2016 CNS WHO), IDH-mutant astrocytic gliomas comprised WHO grade II diffuse astrocytoma, IDH-mutant (AII_IDHmut), WHO grade III anaplastic astrocytoma, IDH-mutant (AAIII_IDHmut), and WHO grade IV glioblastoma, IDH-mutant (GBM_IDHmut). Notably, IDH gene status has been made the major criterion for classification while the manner of grading has remained unchanged: it is based on histological criteria that arose from studies which antedated knowledge of the importance of IDH status in diffuse astrocytic tumor prognostic assessment. Several studies have now demonstrated that the anticipated differences in survival between the newly defined AII_IDHmut and AAIII_IDHmut have lost their significance. In contrast, GBM_IDHmut still exhibits a significantly worse outcome than its lower grade IDH-mutant counterparts. To address the problem of establishing prognostically significant grading for IDH-mutant astrocytic gliomas in the IDH era, we undertook a comprehensive study that included assessment of histological and genetic approaches to prognosis in these tumors. A discovery cohort of 211 IDH-mutant astrocytic gliomas with an extended observation was subjected to histological review, image analysis, and DNA methylation studies. Tumor group-specific methylation profiles and copy number variation (CNV) profiles were established for all gliomas. Algorithms for automated CNV analysis were developed. All tumors exhibiting 1p/19q codeletion were excluded from the series. We developed algorithms for grading, based on molecular, morphological and clinical data. Performance of these algorithms was compared with that of WHO grading. Three independent cohorts of 108, 154 and 224 IDH-mutant astrocytic gliomas were used to validate this approach. In the discovery cohort several molecular and clinical parameters were of prognostic relevance. Most relevant for overall survival (OS) was CDKN2A/B homozygous deletion. Other parameters with major influence were necrosis and the total number of CNV. Proliferation as assessed by mitotic count, which is a key parameter in 2016 CNS WHO grading, was of only minor influence. Employing the parameters most relevant for OS in our discovery set, we developed two models for grading these tumors. These models performed significantly better than WHO grading in both the discovery and the validation sets. Our novel algorithms for grading IDH-mutant astrocytic gliomas overcome the challenges caused by introduction of IDH status into the WHO classification of diffuse astrocytic tumors. We propose that these revised approaches be used for grading of these tumors and incorporated into future WHO criteria.

Cross-sectional and longitudinal changes in DNA methylation with age: an epigenome-wide analysis revealing over 60 novel age-associated CpG sites

Understanding the role of epigenetic modifications, e.g. DNA methylation, in the process of aging requires the characterization of methylation patterns in large cohorts. We analysed >480 000 CpG sites using Infinium HumanMethylation450 BeadChip (Illumina) in whole blood DNA of 965 participants of a population-based cohort study aged between 50 and 75 years. In an exploratory analysis in 400 individuals, 200 CpG sites with the highest Spearman correlation coefficients for the association between methylation and age were identified. Of these 200 CpGs, 162 were significantly associated with age, which was verified in an independent cohort of 498 individuals using mixed linear regression models adjusted for gender, smoking behaviour, age-related diseases and random batch effect and corrected for multiple testing by Bonferroni. In another independent cohort of 67 individuals without history of major age-related diseases and with a follow-up of 8 years, we observed a gain in methylation at 96% (52%, significant) of the positively age-associated CpGs and a loss at all (89%, significant) of the negatively age-associated CpGs in each individual while getting 8 years older. A regression model for age prediction based on 17 CpGs as predicting variables explained 71% of the variance in age with an average accuracy of 2.6 years. In comparison with cord blood samples obtained from the Ulm Birth Cohort Study, we observed a more than 2-fold change in mean methylation levels from birth to older age at 86 CpGs. We were able to identify 65 novel CpG sites with significant association of methylation with age.

Next-generation personalised medicine for high-risk paediatric cancer patients - The INFORM pilot study

The 'Individualized Therapy for Relapsed Malignancies in Childhood' (INFORM) precision medicine study is a nationwide German program for children with high-risk relapsed/refractory malignancies, which aims to identify therapeutic targets on an individualised basis. In a pilot phase, reported here, we developed the logistical and analytical pipelines necessary for rapid and comprehensive molecular profiling in a clinical setting. Fifty-seven patients from 20 centers were prospectively recruited. Malignancies investigated included sarcomas (n = 25), brain tumours (n = 23), and others (n = 9). Whole-exome, low-coverage whole-genome, and RNA sequencing were complemented with methylation and expression microarray analyses. Alterations were assessed for potential targetability according to a customised prioritisation algorithm and subsequently discussed in an interdisciplinary molecular tumour board. Next-generation sequencing data were generated for 52 patients, with the full analysis possible in 46 of 52. Turnaround time from sample receipt until first report averaged 28 d. Twenty-six patients (50%) harbored a potentially druggable alteration with a prioritisation score of 'intermediate' or higher (level 4 of 7). Common targets included receptor tyrosine kinases, phosphoinositide 3-kinase-mammalian target of rapamycin pathway, mitogen-activated protein kinase pathway, and cell cycle control. Ten patients received a targeted therapy based on these findings, with responses observed in some previously treatment-refractory tumours. Comparative primary relapse analysis revealed substantial tumour evolution as well as one case of unsuspected secondary malignancy, highlighting the importance of re-biopsy at relapse. This study demonstrates the feasibility of comprehensive, real-time molecular profiling for high-risk paediatric cancer patients. This extended proof-of-concept, with examples of treatment consequences, expands upon previous personalised oncology endeavors, and presents a model with considerable interest and practical relevance in the burgeoning era of personalised medicine.

Adult IDH wild type astrocytomas biologically and clinically resolve into other tumor entities

IDH wild type (IDHwt) anaplastic astrocytomas WHO grade III (AA III) are associated with poor outcome. To address the possibilities of molecular subsets among astrocytoma or of diagnostic reclassification, we analyzed a series of 160 adult IDHwt tumors comprising 120 AA III and 40 diffuse astrocytomas WHO grade II (A II) for molecular hallmark alterations and established methylation and copy number profiles. Based on molecular profiles and hallmark alterations the tumors could be grouped into four major sets. 124/160 (78 %) tumors were diagnosed as the molecular equivalent of conventional glioblastoma (GBM), and 15/160 (9 %) as GBM-H3F3A mutated (GBM-H3). 13/160 (8 %) exhibited a distinct methylation profile that was most similar to GBM-H3-K27, however, lacked the H3F3A mutation. This group was enriched for tumors of infratentorial and midline localization and showed a trend towards a more favorable prognosis. All but one of the 120 IDHwt AA III could be assigned to these three groups. 7 tumors recruited from the 40 A II, comprised a variety of molecular signatures and all but one were reclassified into distinct WHO entities of lower grades. Interestingly, TERT mutations were exclusively restricted to the molecular GBM (78 %) and associated with poor clinical outcome. However, the GBM-H3 group lacking TERT mutations appeared to fare even worse. Our data demonstrate that most of the tumors diagnosed as IDHwt astrocytomas can be allocated to other tumor entities on a molecular basis. The diagnosis of IDHwt diffuse astrocytoma or anaplastic astrocytoma should be used with caution.

Next-generation sequencing in routine brain tumor diagnostics enables an integrated diagnosis and identifies actionable targets

With the number of prognostic and predictive genetic markers in neuro-oncology steadily growing, the need for comprehensive molecular analysis of neuropathology samples has vastly increased. We therefore developed a customized enrichment/hybrid-capture-based next-generation sequencing (NGS) gene panel comprising the entire coding and selected intronic and promoter regions of 130 genes recurrently altered in brain tumors, allowing for the detection of single nucleotide variations, fusions, and copy number aberrations. Optimization of probe design, library generation and sequencing conditions on 150 samples resulted in a 5-workday routine workflow from the formalin-fixed paraffin-embedded sample to neuropathological report. This protocol was applied to 79 retrospective cases with established molecular aberrations for validation and 71 prospective cases for discovery of potential therapeutic targets. Concordance of NGS compared to established, single biomarker methods was 98.0 %, with discrepancies resulting from one case where a TERT promoter mutation was not called by NGS and three ATRX mutations not being detected by Sanger sequencing. Importantly, in samples with low tumor cell content, NGS was able to identify mutant alleles that were not detectable by traditional methods. Information derived from NGS data identified potential targets for experimental therapy in 37/47 (79 %) glioblastomas, 9/10 (90 %) pilocytic astrocytomas, and 5/14 (36 %) medulloblastomas in the prospective target discovery cohort. In conclusion, we present the settings for high-throughput, adaptive next-generation sequencing in routine neuropathology diagnostics. Such an approach will likely become highly valuable in the near future for treatment decision making, as more therapeutic targets emerge and genetic information enters the classification of brain tumors.

Integrated Molecular-Morphologic Meningioma Classification: A Multicenter Retrospective Analysis, Retrospectively and Prospectively Validated

PURPOSE

Meningiomas are the most frequent primary intracranial tumors. Patient outcome varies widely from benign to highly aggressive, ultimately fatal courses. Reliable identification of risk of progression for individual patients is of pivotal importance. However, only biomarkers for highly aggressive tumors are established (CDKN2A/B and TERT), whereas no molecularly based stratification exists for the broad spectrum of patients with low- and intermediate-risk meningioma.

METHODS

DNA methylation data and copy-number information were generated for 3,031 meningiomas (2,868 patients), and mutation data for 858 samples. DNA methylation subgroups, copy-number variations (CNVs), mutations, and WHO grading were analyzed. Prediction power for outcome was assessed in a retrospective cohort of 514 patients, validated on a retrospective cohort of 184, and on a prospective cohort of 287 multicenter cases.

RESULTS

Both CNV- and methylation family-based subgrouping independently resulted in increased prediction accuracy of risk of recurrence compared with the WHO classification (c-indexes WHO 2016, CNV, and methylation family 0.699, 0.706, and 0.721, respectively). Merging all risk stratification approaches into an integrated molecular-morphologic score resulted in further substantial increase in accuracy (c-index 0.744). This integrated score consistently provided superior accuracy in all three cohorts, significantly outperforming WHO grading (c-index difference P = .005). Besides the overall stratification advantage, the integrated score separates more precisely for risk of progression at the diagnostically challenging interface of WHO grade 1 and grade 2 tumors (hazard ratio 4.34 [2.48-7.57] and 3.34 [1.28-8.72] retrospective and prospective validation cohorts, respectively).

CONCLUSION

Merging these layers of histologic and molecular data into an integrated, three-tiered score significantly improves the precision in meningioma stratification. Implementation into diagnostic routine informs clinical decision making for patients with meningioma on the basis of robust outcome prediction.

Methylation-based classification of benign and malignant peripheral nerve sheath tumors

The vast majority of peripheral nerve sheath tumors derive from the Schwann cell lineage and comprise diverse histological entities ranging from benign schwannomas and neurofibromas to high-grade malignant peripheral nerve sheath tumors (MPNST), each with several variants. There is increasing evidence for methylation profiling being able to delineate biologically relevant tumor groups even within the same cellular lineage. Therefore, we used DNA methylation arrays for methylome- and chromosomal profile-based characterization of 171 peripheral nerve sheath tumors. We analyzed 28 conventional high-grade MPNST, three malignant Triton tumors, six low-grade MPNST, four epithelioid MPNST, 33 neurofibromas (15 dermal, 8 intraneural, 10 plexiform), six atypical neurofibromas, 43 schwannomas (including 5 NF2 and 5 schwannomatosis associated cases), 11 cellular schwannomas, 10 melanotic schwannomas, 7 neurofibroma/schwannoma hybrid tumors, 10 nerve sheath myxomas and 10 ganglioneuromas. Schwannomas formed different epigenomic subgroups including a vestibular schwannoma subgroup. Cellular schwannomas were not distinct from conventional schwannomas. Nerve sheath myxomas and neurofibroma/schwannoma hybrid tumors were most similar to schwannomas. Dermal, intraneural and plexiform neurofibromas as well as ganglioneuromas all showed distinct methylation profiles. Atypical neurofibromas and low-grade MPNST were indistinguishable with a common methylation profile and frequent losses of CDKN2A. Epigenomic analysis finds two groups of conventional high-grade MPNST sharing a frequent loss of neurofibromin. The larger of the two groups shows an additional loss of trimethylation of histone H3 at lysine 27 (H3K27me3). The smaller one retains H3K27me3 and is found in spinal locations. Sporadic MPNST with retained neurofibromin expression did not form an epigenetic group and most cases could be reclassified as cellular schwannomas or soft tissue sarcomas. Widespread immunohistochemical loss of H3K27me3 was exclusively seen in MPNST of the main methylation cluster, which defines it as an additional useful marker for the differentiation of cellular schwannoma and MPNST.

Type 2 diabetes and leucocyte DNA methylation: an epigenome-wide association study in over 1,500 older adults

AIMS/HYPOTHESIS

Development of type 2 diabetes depends on environmental and genetic factors. We investigated the epigenome-wide association of prevalent diabetes with DNA methylation (DNAm) in peripheral blood.

METHODS

DNAm was measured in whole blood with the Illumina Infinium HumanMethylation450 BeadChip in two subsamples of participants from the ESTHER cohort study. Cohort 1 included 988 participants, who were consecutively recruited between July and October 2000 and cohort 2 included 527 randomly selected participants. The association of DNAm with prevalent type 2 diabetes at recruitment was estimated using median regression analysis adjusting for sex, age, BMI, smoking behaviour, cell composition and batch at 361,922 CpG sites.

RESULTS

Type 2 diabetes was prevalent in 16% of the participants, and diabetes was poorly controlled in 45% of the diabetic patients. In cohort 1 (discovery) DNAm at 39 CpGs was significantly associated with prevalent diabetes after correction for multiple testing. In cohort 2 (replication) at one of these CpGs, DNAm was still significantly associated. Decreasing methylation levels at cg19693031 with increasing fasting glucose and HbA1c concentrations were observed using restricted cubic spline analysis. In diabetic patients with poorly controlled diabetes, the decrease in estimated DNAm levels was approximately 5% in comparison with participants free of diagnosed diabetes.

CONCLUSIONS/INTERPRETATION

Cg19693031, which is located within the 3'-untranslated region of TXNIP, might play a role in the pathophysiology of type 2 diabetes. This result appears biologically plausible given that thioredoxin-interacting protein is overexpressed in diabetic animals and humans and 3'-untranslated regions are known to play a regulatory role in gene expression.

MEST mediates the impact of prenatal bisphenol A exposure on long-term body weight development

Background

Exposure to endocrine-disrupting chemicals can alter normal physiology and increase susceptibility to non-communicable diseases like obesity. Especially the prenatal and early postnatal period is highly vulnerable to adverse effects by environmental exposure, promoting developmental reprogramming by epigenetic alterations. To obtain a deeper insight into the role of prenatal bisphenol A (BPA) exposure in children's overweight development, we combine epidemiological data with experimental models and BPA-dependent DNA methylation changes.

Methods

BPA concentrations were measured in maternal urine samples of the LINA mother-child-study obtained during pregnancy (n = 552), and BPA-associated changes in cord blood DNA methylation were analyzed by Illumina Infinium HumanMethylation450 BeadChip arrays (n = 472). Methylation changes were verified by targeted MassARRAY analyses, assessed for their functional translation by qPCR and correlated with children's body mass index (BMI) z scores at the age of 1 and 6 years. Further, female BALB/c mice were exposed to BPA from 1 week before mating until delivery, and weight development of their pups was monitored (n ≥ 8/group). Additionally, human adipose-derived mesenchymal stem cells were treated with BPA during the adipocyte differentiation period and assessed for exposure-related epigenetic, transcriptional and morphological changes (n = 4).

Results

In prenatally BPA-exposed children two CpG sites with deviating cord blood DNA-methylation profiles were identified, among them a hypo-methylated CpG in the promoter of the obesity-associated mesoderm-specific transcript (MEST). A mediator analysis suggested that prenatal BPA exposure was connected to cord blood MEST promoter methylation and MEST expression as well as BMI z scores in early infancy. This effect could be confirmed in mice in which prenatal BPA exposure altered Mest promoter methylation and transcription with a concomitant increase in the body weight of the juvenile offspring. An experimental model of in vitro differentiated human mesenchymal stem cells also revealed an epigenetically induced MEST expression and enhanced adipogenesis following BPA exposure.

Conclusions

Our study provides evidence that MEST mediates the impact of prenatal BPA exposure on long-term body weight development in offspring by triggering adipocyte differentiation.

Integrated molecular characterization of IDH-mutant glioblastomas

AIMS

Mutations of isocitrate dehydrogenase (IDH)1/2 affect almost all astrocytomas of WHO grade II and III. A subset of IDH-mutant astrocytic tumours progresses to IDH-mutant glioblastoma or presents with the histology of a glioblastoma at first presentation. We set out here to assess the molecular spectrum of IDH-mutant glioblastomas.

METHODS

We performed an integrated molecular analysis of a mono-centric cohort (n = 97); assessed through genome-wide DNA methylation analysis, copy-number profiling and targeted next generation sequencing using a neurooncology-tailored gene panel.

RESULTS

Of these 97 IDH-mutant glioblastomas, 68 had a glioblastoma at first presentation ('de novo' IDH-mutant glioblastoma) and 29 emerged from a prior low-grade lesion ('evolved' IDH-mutant glioblastoma). Unsupervised hierarchical clustering of DNA methylation data disclosed that IDH-mutant glioblastoma ('de novo' and 'evolved') formed a distinct group separate from other diffuse glioma subtypes. Homozygous deletions of CDKN2A/B were found to be associated with shorter survival.

CONCLUSIONS

This study demonstrates DNA methylation patterns in IDH-mutant glioblastoma to be distinct from lower-grade astrocytic counterparts but homogeneous within de novo and evolved IDH-mutant glioblastomas, and identifies CDKN2A as a marker for possible genetic sub-stratification.

Methylation profiling of choroid plexus tumors reveals 3 clinically distinct subgroups

BACKGROUND

Choroid plexus tumors are intraventricular neoplasms derived from the choroid plexus epithelium. A better knowledge of molecular factors involved in choroid plexus tumor biology may aid in identifying patients at risk for recurrence.

METHODS

Methylation profiles were examined in 29 choroid plexus papillomas (CPPs, WHO grade I), 32 atypical choroid plexus papillomas (aCPPs, WHO grade II), and 31 choroid plexus carcinomas (CPCs, WHO grade III) by Illumina Infinium HumanMethylation450 Bead Chip Array.

RESULTS

Unsupervised hierarchical clustering identified 3 subgroups: methylation cluster 1 (pediatric CPP and aCPP of mainly supratentorial location), methylation cluster 2 (adult CPP and aCPP of mainly infratentorial location), and methylation cluster 3 (pediatric CPP, aCPP, and CPC of supratentorial location). In methylation cluster 3, progression-free survival (PFS) accounted for a mean of 72 months (CI, 55-89 mo), whereas only 1 of 42 tumors of methylation clusters 1 and 2 progressed (P< .001). On stratification of outcome data according to WHO grade, all CPCs clustered within cluster 3 and were associated with shorter overall survival (mean, 105 mo [CI, 81-128 mo]) and PFS (mean, 55 mo [CI, 36-73 mo]). The aCPP of methylation cluster 3 also progressed frequently (mean, 69 mo [CI, 44-93 mo]), whereas no tumor progression was observed in aCPP of methylation clusters 1 and 2 (P< .05). Only 1 of 29 CPPs recurred.

CONCLUSIONS

Methylation profiling of choroid plexus tumors reveals 3 distinct subgroups (ie, pediatric low-risk choroid plexus tumors [cluster 1], adult low-risk choroid plexus tumors [cluster 2], and pediatric high-risk choroid plexus tumors [cluster 3]) and may provide useful prognostic information in addition to histopathology.

Mutational patterns and regulatory networks in epigenetic subgroups of meningioma

DNA methylation patterns delineate clinically relevant subgroups of meningioma. We previously established the six meningioma methylation classes (MC) benign 1-3, intermediate A and B, and malignant. Here, we set out to identify subgroup-specific mutational patterns and gene regulation. Whole genome sequencing was performed on 62 samples across all MCs and WHO grades from 62 patients with matched blood control, including 40 sporadic meningiomas and 22 meningiomas arising after radiation (Mrad). RNA sequencing was added for 18 of these cases and chromatin-immunoprecipitation for histone H3 lysine 27 acetylation (H3K27ac) followed by sequencing (ChIP-seq) for 16 samples. Besides the known mutations in meningioma, structural variants were found as the mechanism of NF2 inactivation in a small subset (5%) of sporadic meningiomas, similar to previous reports for Mrad. Aberrations of DMD were found to be enriched in MCs with NF2 mutations, and DMD was among the most differentially upregulated genes in NF2 mutant compared to NF2 wild-type cases. The mutational signature AC3, which has been associated with defects in homologous recombination repair (HRR), was detected in both sporadic meningioma and Mrad, but widely distributed across the genome in sporadic cases and enriched near genomic breakpoints in Mrad. Compared to the other MCs, the number of single nucleotide variants matching the AC3 pattern was significantly higher in the malignant MC, which also exhibited higher genomic instability, determined by the numbers of both large segments affected by copy number alterations and breakpoints between large segments. ChIP-seq analysis for H3K27ac revealed a specific activation of genes regulated by the transcription factor FOXM1 in the malignant MC. This analysis also revealed a super enhancer near the HOXD gene cluster in this MC, which, together with general upregulation of HOX genes in the malignant MC, indicates a role of HOX genes in meningioma aggressiveness. This data elucidates the biological mechanisms rendering different epigenetic subgroups of meningiomas, and suggests leveraging HRR as a novel therapeutic target.

Histone 3.3 hotspot mutations in conventional osteosarcomas: a comprehensive clinical and molecular characterization of six H3F3A mutated cases

BACKGROUND

Histone 3.3 (H3.3) hotspot mutations in bone tumors occur in the vast majority of giant cell tumors of bone (GCTBs; 96%), chondroblastomas (95%) and in a few cases of osteosarcomas. However, clinical presentation, histopathological features, and additional molecular characteristics of H3.3 mutant osteosarcomas are largely unknown.

METHODS

In this multicentre, retrospective study, a total of 106 conventional high-grade osteosarcomas, across all age groups were re-examined for hotspot mutations in the H3.3 coding genes H3F3A and H3F3B. H3.3 mutant osteosarcomas were re-evaluated in a multidisciplinary manner and analyzed for genome-wide DNA-methylation patterns and DNA copy number aberrations alongside H3.3 wild-type osteosarcomas and H3F3A G34W/L mutant GCTBs.

RESULTS

Six osteosarcomas (6/106) carried H3F3A hotspot mutations. No mutations were found in H3F3B. All patients with H3F3A mutant osteosarcoma were older than 30 years with a median age of 65 years. Copy number aberrations that are commonly encountered in high-grade osteosarcomas also occurred in H3F3A mutant osteosarcomas. Unlike a single osteosarcoma with a H3F3A K27M mutation, the DNA methylation profiles of H3F3A G34W/R mutant osteosarcomas were clearly different from H3.3 wild-type osteosarcomas, but more closely related to GCTBs. The most differentially methylated promoters between H3F3A G34W/R mutant and H3.3 wild-type osteosarcomas were in KLLN/PTEN (p < 0.00005) and HIST1H2BB (p < 0.0005).

CONCLUSIONS

H3.3 mutations in osteosarcomas may occur in H3F3A at mutational hotspots. They are overall rare, but become more frequent in osteosarcoma patients older than 30 years. Osteosarcomas carrying H3F3A G34W/R mutations are associated with epigenetic dysregulation of KLLN/PTEN and HIST1H2BB.

Genome-wide methylation profiling and copy number analysis in atypical fibroxanthomas and pleomorphic dermal sarcomas indicate a similar molecular phenotype

Background

Atypical fibroxanthomas (AFX) and pleomorphic dermal sarcomas (PDS) are lesions of the skin with overlapping histologic features and unspecific molecular traits. PDS behaves aggressive compared to AFX. Thus, a precise delineation, although challenging in some instances, is relevant.

Methods

We examined the value of DNA-methylation profiling and copy number analysis for separating these tumors. DNA-methylation data were generated from 17 AFX and 15 PDS using the Illumina EPIC array. These were compared with DNA-methylation data generated from 196 tumors encompassing potential histologic mimics like cutaneous squamous carcinomas (cSCC; n = 19), basal cell carcinomas (n = 10), melanoma metastases originating from the skin (n = 11), leiomyosarcomas (n = 11), angiosarcomas of the skin and soft tissue (n = 11), malignant peripheral nerve sheath tumors (n = 19), dermatofibrosarcomas protuberans (n = 13), extraskeletal myxoid chondrosarcomas (n = 9), myxoid liposarcomas (n = 14), schwannomas (n = 10), neurofibromas (n = 21), alveolar (n = 19) and embryonal (n = 17) rhabdomyosarcomas as well as undifferentiated pleomorphic sarcomas (n = 12).

Results

DNA-methylation profiling did not separate AFX from PDS. The DNA-methylation profiles of the other cases, however, were distinct from AFX/PDS. They reliably assigned to subtype-specific DNA-methylation clusters, although overlap occurred between some AFX/PDS and cSCC. Copy number profiling revealed alterations in a similar frequency and distribution between AFX and PDS. They involved losses of 9p (22/32) and 13q (25/32). Gains frequently involved 8q (8/32). Notably, a homozygous deletion of CDKN2A was more frequent in PDS (6/15) than in AFX (2/17), whereas amplifications were non-recurrent and overall rare (5/32).

Conclusions

Our findings support the concept that AFX and PDS belong to a common tumor spectrum. We could demonstrate the diagnostic value of DNA-methylation profiling to delineating AFX/PDS from potential mimics. However, the assessment of certain histologic features remains crucial for separating PDS from AFX.

Molecular differences in IDH wildtype glioblastoma according to MGMT promoter methylation

Background

O6-methylguanine-DNA-methyltransferase (MGMT) promoter methylation status is a predictive biomarker in glioblastoma. We investigated whether this marker furthermore defines a molecularly distinct tumor subtype with clinically different outcome.

Methods

We analyzed copy number variation (CNV) and methylation profiles of 1095 primary and 92 progressive isocitrate dehydrogenase wildtype glioblastomas, including paired samples from 49 patients. DNA mutation data from 182 glioblastoma samples of The Cancer Genome Atlas (TCGA) and RNA expression from 107 TCGA and 55 Chinese Glioma Genome Atlas samples were analyzed.

Results

Among untreated glioblastomas, MGMT promoter methylated (mMGMT) and unmethylated (uMGMT) tumors did not show different CNV or specific gene mutations, but a higher mutation count in mMGMT tumors. We identified 3 methylation clusters. Cluster 1 showed the highest average methylation and was enriched for mMGMT tumors. Seventeen genes including gastrulation brain homeobox 2 (GBX2) were found to be hypermethylated and downregulated on the mRNA level in mMGMT tumors. In progressive glioblastomas, platelet derived growth factor receptor alpha (PDGFRA) and GLI2 amplifications were enriched in mMGMT tumors. Methylated MGMT tumors gain PDGFRA amplification of PDGFRA, whereas uMGMT tumors with amplified PDGFRA frequently lose this amplification upon progression. Glioblastoma patients surviving <6 months and with mMGMT harbored less frequent epidermal growth factor receptor (EGFR) amplifications, more frequent TP53 mutations, and a higher tumor necrosis factor-nuclear factor-kappaB (TNF-NFκB) pathway activation compared with patients surviving >12 months.

Conclusions

MGMT promoter methylation status does not define a molecularly distinct glioblastoma subpopulation among untreated tumors. Progressive mMGMT glioblastomas and mMGMT tumors of patients with short survival tend to have more unfavorable molecular profiles.

Profiling of gallbladder carcinoma reveals distinct miRNA profiles and activation of STAT1 by the tumor suppressive miRNA-145-5p

Gallbladder carcinoma (GBC) is a biliary tract cancer with few treatment options and poor prognosis. Radical surgery is the only potentially curative treatment option but most patients diagnosed with GBC are unresectable. Thus, there is a great need for the development of new treatment options including targeted therapy. Here, we aimed at identifying deregulated miRNAs and affected pathways involved in GBC development and progression. We performed global miRNA profiling of 40 GBC and 8 normal gallbladder tissues and identified large differences with 30% of miRNAs being differentially expressed (false discovery rate: FDR < 0.001). We found 24 miRNAs to be differentially regulated in GBC with poor outcome (p < 0.05) of which miR-145-5p was the most downregulated miRNA. Overexpression of miR-145-5p significantly reduced cell proliferation and colony formation. Gene expression analysis of cells expressing miR-145-5p mimics revealed activation of the Signal transducer and activator of transcription 1 (STAT1) signaling pathway which is mainly tumor suppressive. Furthermore, the activation of STAT1 by miR-145-5p was specifically observed in gallbladder carcinoma and cholangiocarcinoma but not in hepatocellular carcinoma cells. The Protein Tyrosine Phosphatase Receptor Type F (PTPRF) is downregulated upon miR-145 expression and may be involved in STAT1 regulation. In addition, we found that the STAT1-regulated protein IRF7 is downregulated in GBC compared to normal gallbladder tissue and low IRF7 expression is associated with significantly lower overall survival of GBC patients. Thus, this study identified GBC patient subgroups and provides new mechanistic insights in the tumor suppressive function of miR-145-5p leading to activation of STAT1 signaling.

Systematic screening of isogenic cancer cells identifies DUSP6 as context-specific synthetic lethal target in melanoma

Next-generation sequencing has dramatically increased genome-wide profiling options and conceptually initiates the possibility for personalized cancer therapy. State-of-the-art sequencing studies yield large candidate gene sets comprising dozens or hundreds of mutated genes. However, few technologies are available for the systematic downstream evaluation of these results to identify novel starting points of future cancer therapies.

We improved and extended a site-specific recombination-based system for systematic analysis of the individual functions of a large number of candidate genes. This was facilitated by a novel system for the construction of isogenic constitutive and inducible gain- and loss-of-function cell lines. Additionally, we demonstrate the construction of isogenic cell lines with combinations of the traits for advanced functional in vitro analyses. In a proof-of-concept experiment, a library of 108 isogenic melanoma cell lines was constructed and 8 genes were identified that significantly reduced viability in a discovery screen and in an independent validation screen. Here, we demonstrate the broad applicability of this recombination-based method and we proved its potential to identify new drug targets via the identification of the tumor suppressor DUSP6 as potential synthetic lethal target in melanoma cell lines with BRAF V600E mutations and high DUSP6 expression.

DNA-methylome-assisted classification of patients with poor prognostic subventricular zone associated IDH-wildtype glioblastoma

Glioblastoma (GBM) derived from the “stem cell” rich subventricular zone (SVZ) may constitute a therapy-refractory subgroup of tumors associated with poor prognosis. Risk stratification for these cases is necessary but is curtailed by error prone imaging-based evaluation. Therefore, we aimed to establish a robust DNA methylome-based classification of SVZ GBM and subsequently decipher underlying molecular characteristics. MRI assessment of SVZ association was performed in a retrospective training set of IDH-wildtype GBM patients (n = 54) uniformly treated with postoperative chemoradiotherapy. DNA isolated from FFPE samples was subject to methylome and copy number variation (CNV) analysis using Illumina Platform and cnAnalysis450k package. Deep next-generation sequencing (NGS) of a panel of 130 GBM-related genes was conducted (Agilent SureSelect/Illumina). Methylome, transcriptome, CNV, MRI, and mutational profiles of SVZ GBM were further evaluated in a confirmatory cohort of 132 patients (TCGA/TCIA). A 15 CpG SVZ methylation signature (SVZM) was discovered based on clustering and random forest analysis. One third of CpG in the SVZM were associated with MAB21L2/LRBA. There was a 14.8% (n = 8) discordance between SVZM vs. MRI classification. Re-analysis of these patients favored SVZM classification with a hazard ratio (HR) for OS of 2.48 [95% CI 1.35–4.58], p = 0.004 vs. 1.83 [1.0–3.35], p = 0.049 for MRI classification. In the validation cohort, consensus MRI based assignment was achieved in 62% of patients with an intraclass correlation (ICC) of 0.51 and non-significant HR for OS (2.03 [0.81–5.09], p = 0.133). In contrast, SVZM identified two prognostically distinct subgroups (HR 3.08 [1.24–7.66], p = 0.016). CNV alterations revealed loss of chromosome 10 in SVZM– and gains on chromosome 19 in SVZM– tumors. SVZM– tumors were also enriched for differentially mutated genes (p < 0.001). In summary, SVZM classification provides a novel means for stratifying GBM patients with poor prognosis and deciphering molecular mechanisms governing aggressive tumor phenotypes.

A complementary study approach unravels novel players in the pathoetiology of Hirschsprung disease

Hirschsprung disease (HSCR, OMIM 142623) involves congenital intestinal obstruction caused by dysfunction of neural crest cells and their progeny during enteric nervous system (ENS) development. HSCR is a multifactorial disorder; pathogenetic variants accounting for disease phenotype are identified only in a minority of cases, and the identification of novel disease-relevant genes remains challenging. In order to identify and to validate a potential disease-causing relevance of novel HSCR candidate genes, we established a complementary study approach, combining whole exome sequencing (WES) with transcriptome analysis of murine embryonic ENS-related tissues, literature and database searches, in silico network analyses, and functional readouts using candidate gene-specific genome-edited cell clones. WES datasets of two patients with HSCR and their non-affected parents were analysed, and four novel HSCR candidate genes could be identified: ATP7A, SREBF1, ABCD1 and PIAS2. Further rare variants in these genes were identified in additional HSCR patients, suggesting disease relevance. Transcriptomics revealed that these genes are expressed in embryonic and fetal gastrointestinal tissues. Knockout of these genes in neuronal cells demonstrated impaired cell differentiation, proliferation and/or survival. Our approach identified and validated candidate HSCR genes and provided further insight into the underlying pathomechanisms of HSCR.