Integrated genomic analyses of de novo pathways

Full text

Turn on search term navigation

ARTICLE

Received 2 Aug 2016 | Accepted 28 Dec 2016 | Published 14 Feb 2017

Akdes Serin Harmanc1,2, Mark W. Youngblood1,2,3, Victoria E. Clark1,2,3, Sleyman Coskun1,2, Octavian Henegariu1,2,3,4,5, Daniel Duran1,2, E. Zeynep Erson-Omay1,2, Leon D. Kaulen1,2, Tong Ihn Lee6, Brian J. Abraham6,

Matthias Simon7, Boris Krischek8, Marco Timmer8, Roland Goldbrunner8, S. Blent Omay1,2, Jacob Baranoski1,2,3, Burin Baran1,2, Geneive Carrin-Grant1,2, Hanwen Bai1,3, Ketu Mishra-Gorur1,2,3,4,5, Johannes Schramm7, Jennifer Moliterno1,2, Alexander O. Vortmeyer9, Kaya Bilgvar1,3,10, Katsuhito Yasuno1,2,Richard A. Young6,11 & Murat Gnel1,2,3,4,5,12

Meningiomas are mostly benign brain tumours, with a potential for becoming atypical or malignant. On the basis of comprehensive genomic, transcriptomic and epigenomic analyses, we compared benign meningiomas to atypical ones. Here, we show that the majority of primary (de novo) atypical meningiomas display loss of NF2, which co-occurs either with genomic instability or recurrent SMARCB1 mutations. These tumours harbour increased H3K27me3 signal and a hypermethylated phenotype, mainly occupying the polycomb repressive complex 2 (PRC2) binding sites in human embryonic stem cells, thereby phenocopying a more primitive cellular state. Consistent with this observation, atypical meningiomas exhibit upregulation of EZH2, the catalytic subunit of the PRC2 complex, as well as the E2F2 and FOXM1 transcriptional networks. Importantly, these primary atypical meningiomas do not harbour TERT promoter mutations, which have been reported in atypical tumours that progressed from benign ones. Our results establish the genomic landscape of primary atypical meningiomas and potential therapeutic targets.

1 Yale Program in Brain Tumor Research, Yale School of Medicine, New Haven, Connecticut 06510, USA. 2 Department of Neurosurgery, Yale School of Medicine, New Haven, Connecticut 06510, USA. 3 Department of Genetics, Yale School of Medicine, New Haven, Connecticut 06510, USA. 4 Department of Neurobiology, Yale School of Medicine, New Haven, Connecticut 06510, USA. 5 Yale Program on Neurogenetics, Yale School of Medicine, New Haven, Connecticut 06510, USA. 6 Whitehead Institute for Biomedical Research, 9 Cambridge Center, Cambridge, Massachusetts 02142, USA. 7 Department of Neurosurgery, University of Bonn Medical School, Bonn 53105, Germany. 8 Department of General Neurosurgery, University Hospital of Cologne, Cologne 50937, Germany. 9 Department of Pathology, Yale School of Medicine, New Haven, Connecticut 06510, USA. 10 Yale Center for Genome Analysis, Yale School of Medicine, Orange, Connecticut 06477, USA. 11 Department of Biology, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA. 12 Yale Comprehensive Cancer Center, Yale School of Medicine, New Haven, Connecticut 06510, USA. Correspondence and requests for materials should be addressed to M.G. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 8:14433 | DOI: 10.1038/ncomms14433 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

DOI: 10.1038/ncomms14433 OPEN

Integrated genomic analyses of de novo pathways underlying atypical meningiomas

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14433

Over 35% of all primary tumours that affect the central nervous system are meningiomas, which originate from the three-layer meningeal membrane ensheathing the

brain and spinal cord1. According to the World Health Organization (WHO) criteria, meningiomas are classied into three pathological grades, based mainly on morphologic ndings2,3. These include histological criteria such as mitotic activity, cellularity, cellular morphology and growth pattern, necrosis, and brain invasion. Approximately 7080% of meningiomas are grade I and benign, while grade II and III meningiomas are higher grade and classied as atypical (520%) or malignant (13%), respectively2,3. The grading of a tumour carries prognostic value, as higher grade lesions are more likely to recur and decrease the chances of long-term survival4.

Work by our lab and others has identied mutually exclusive molecular subgroups in benign meningioma, including loss of NF2 (occasionally with recurrent mutations in SMARCB1), mutations in the WD40 repeat region of TRAF7 (co-occurring with either PI3K activating mutations or recurrent KLFp.Lys409Gln mutation), activation of Hedgehog signalling (via SMO, SUFU or PRKAR1A), and recurrent p.Gln403Lys and p.Leu438_His439del mutations in the dock domain of POLR2A57. Interestingly, these subgroups were associated with distinct pathological and clinical ndings. For example, secretory meningiomas were driven exclusively by TRAF7/KLF4 co-mutations, while brous meningiomas were primarily associated with NF2 loss. The intracranial origin of a meningioma was also predicted by the underlying meningioma mutations, with non-NF2 mutant tumours being enriched in the neural crest derived anterior skull base region, while samples harbouring NF2 loss arose from the mesoderm-derived posterior regions6.

While these studies led to comprehensive genomic characterization of benign meningiomas, the genomic pathways that lead to formation of atypical cases are not well established. Primary atypical meningiomas form de novo, whereas secondary ones form due to recurrence and malignant progression of benign meningiomas8. Overall, atypical tumours comprise B520% of all meningiomas and are associated with poor prognosis and a 10-year survival less than 80%, mostly due to local recurrence9. As with other high-grade forms of neoplasia, these tumours show increased genomic instability, including loss of chromosomes 22q, 1p and 14q6,10. Other than identifying these large chromosomal events, studies to date have failed to identify specic molecular changes that distinguish atypical from benign meningiomas.

The genomic alterations that differentiate low-grade and high-grade tumours have been extensively studied in other forms of neoplasia11. While somatic coding mutations have been identied in many cases, epigenetic modications have also emerged as a potent mechanism to induce formation of malignant tumours. These include alterations in histone modications or DNA methylation, both of which remodel chromatin to affect changes in gene expression and alter the transcriptional prole of the cancer cells11,12.

Here, we used integrated genomic techniques to dissect the molecular landscape of primary atypical tumours compared with their benign counterparts. We show that the mutational background of primary atypical meningiomas to be comprised mainly of NF2 mutants, which frequently co-occur with either chromosomal instability or recurrent p.(Arg383Gln) or p.Arg386His mutations in SMARCB1 (co-occurrence P 1.2 10 7, Fishers

exact test). Meningiomas with these alterations carry a signicantly higher risk of being atypical as compared with non-NF2 meningiomas, including TRAF7 (with PI3K or KLF4 alterations), Hedgehog or POLR2A mutant tumours. Differences in the

number of coding mutations do not contribute to the risk of being atypical; while large-scale copy number variant (CNV) events, transcriptional and epigenetic changes as well as alterations in miRNAs show substantial association. Indeed, genomically unstable, hypermethylated NF2 mutant meningiomas, which display activation of the cell cycle as well as the PRC2 pathways, account for the majority of the genomic landscape of primary atypical meningiomas. These ndings dene novel therapeutic targets in atypical meningiomas, which continue to represent signicant treatment challenges due to a lack of effective chemotherapeutics.

ResultsNF2 and SMARCB1 mutations in atypical meningiomas. We hypothesized that similar to gliomas, in which malignancy (glioblastoma multiforme) can occur either de novo or through progression of a lower grade tumour11, the molecular pathways that underlie formation of primary atypical meningiomas would be different than those associated with the progression of benign tumours8. On the basis of this assumption, we choose to specically focus our analysis on primary atypical samples, as the molecular features of this tumour have not been described. Of the 775 meningiomas that we studied using next-generation exome (n 107) or targeted sequencing, we concentrated our initial

efforts on comparing histologically benign samples (n 468)

with de novo atypical tumours (n 88) (Supplementary Fig. 1;

Supplementary Data 1ac). This dataset did not include tumours that have undergone previous chemotherapy or radiation, as these treatments may induce exogenous genomic changes.

We rst searched for known meningioma driver mutations in each sample, dening the distribution of benign versus atypical meningiomas in each of the molecular subgroups including NF2, TRAF7 (co-mutated with PI3K pathway or KLF4), Hedgehog and POLR2A mutant tumours. We found signicant differences in the percentage of atypical versus benign tumours within these molecular subgroups. Of the 88 histologically atypical meningiomas, 75% contained NF2 mutations (Supplementary Data 1c), while the remaining 9% were TRAF7/PI3K mutant and 16% did not harbour a mutation in the previously established meningioma genes. In our sample set, we did not observe any atypical tumours that harboured mutations in TRAF7/KLF4, POLR2A or the Hedgehog pathway. Because of the high prevalence of NF2 alterations in the atypical cohort, a tumour harbouring NF2 loss has a 3.78 times greater risk to be atypical compared with a non-NF2 meningioma (P 2.2 10 10, Fishers exact test).

Given these results, which suggest that primary atypical meningiomas were overwhelmingly associated with loss of NF2, we divided our cohort into two large subgroups: those with NF2 mutations, and those that were NF2 wild type (including TRAF7/PI3K//KLF4, Hedgehog and POLR2A mutant tumours as well as mutation unknown samples). We investigated the potential role of coding variation by comparing the number of somatic coding mutations in atypical and benign samples stratied by NF2 status, which allowed us to control for the underlying driver mutation. We did not nd a statistically signicant difference (Students t-test) (Fig. 1a,b; Supplementary Figs 2 and 3; Supplementary Data 2ae).

Among all samples that underwent whole-exome sequencing, the only somatically mutated gene that we found to be enriched in the atypical samples was NF2 (Supplementary Data 2be). We next studied an extended cohort that included not only the whole-exome sequenced dataset, but also an independent dataset of meningioma samples that underwent targeted sequencing (n 556). We calculated the signicance of association of the

2 NATURE COMMUNICATIONS | 8:14433 | DOI: 10.1038/ncomms14433 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14433 ARTICLE

a b

Deleterious predicted

somatic mutation count

Atypical/benign

Atypical

Benign

Atypical

Benign

chr22 loss

NF2

SMARCB1

AKT1

TRAF7/ PI3K-KLF4

Hedgehog

Nonsense

Frameshift

Splice site

Missense

Inframe InDel

PIK3CA

PIK3R1

KLF4

TRAF7

SMO

POLR2A

NF2 Non-NF2

Somatic mutation count

Non-synonymous

Synonymous

Utr

P= 0.02 P= 0.27

Atypical

Benign

%30

% Altered genome

%20

100%

Mutation

signature

75%

%10

50%

Transition Transversion

25%

NF2 Non-NF2

d e

14q

Gain

Loss

TNFRSF14

50%

10q

FDR= 0.04

FDR= 0.06

FDR= 0.07

FDR= 0.1

Atypical

Benign

% Gain

% Loss

25%

10p

25%

ARID1B

BCL11B

SMAD2

CAMTA1

SPEN SDHB

SFPQ MUTYH

CDKN2C

BCL10

TRIM33

50%

GATA3

PRF1

KAT6B FAS

DICER1

RPL22

FAM46C

NOTCH2

PRDM1

PTPRK

ESR1

CCDC6

BMPR1A

PTEN

MAX

KLF6

%0 %20 %40

% of deletion events in NF2 mutants

2 3 4 5 6 7 8 9 10

Chromosome

Figure 1 | Somatic mutations and copy number variations in benign and atypical meningiomas. (a) Oncoprint depicting the mutational prole of75 exome-sequenced meningioma samples is shown. Histological grade, chromosome 22 loss status and recurrently mutated genes, which are clustered based on mutually exclusive meningioma subgroups (left) are summarized at the top panel. The distributions of somatic mutations according to their functional consequences are shown in the middle panel, whereas the mutational signatures are shown at the bottom. The colour codes are explained on the right. (b) Atypical versus benign meningiomas harbour similar number of damaging somatic mutations (atypical n 10, benign n 65). The lines above the

bars indicate statistical analysis (Students t-test; ns: non-signicant). Lines depict the median values; boxes plot 25th to 75th percentiles, whereas separately plotted dots show the outliers. (c) The percentage of genome alteration is statistically signicantly different between atypical versus benign meningiomas (atypical n 55, benign n 153) (Students t-test). Lines depict the median values; boxes plot 25th to 75th percentiles, whereas separately

plotted dots show the outliers. (d) Large-scale genomic events associated with atypical NF2 samples are shown (atypical n 43, benign n 57). Losses of

chromosomes 14q, 10q, 10p, 1p and 6q signicantly associate with atypical tumours (FDR adjusted Fishers exact test). (e) Differences in genomic alterations between atypical NF2 and benign NF2 mutant samples are shown (atypical n 43, benign n 57). Along the horizontal axis, losses are depicted

in blue, whereas gains are shown in red. The vertical axis represents the genome. Signicantly altered tumour suppressor consensus cancer genes in atypical NF2 samples are noted.

driver mutations (NF2, SMARCB1, TRAF7/PI3K, TRAF7/KLF4, POLR2A, Hedgehog) with atypical meningiomas. We observed a signicant association of being atypical only with NF2 and SMARCB1 mutations (NF2: P 2.2 10 10, SMARCB1:

P 0.05, Fishers exact test). These ndings suggest that except

for mutations in NF2 and SMARCB1, coding variation does not signicantly contribute to the risk of being atypical.

We next calculated the clonality rate of each somatic mutation based on the variant allele frequency and ploidy at that site, considering the admixture rate of each tumour. In the majority of our samples (13/17), we calculated NF2 mutations to have a clonality rate near 100%, suggesting that NF2 mutations occurred early during tumour formation (Supplementary Data 2f; Supplementary Fig. 4).

Atypical NF2 mutants demonstrate chromosomal instability. Using whole-genome genotyping (WGG), we next compared

large-scale CNV events between benign and atypical tumours (n 153 versus 55, benign versus atypical, respectively)

(Supplementary Data 3a). We dened large-scale CNV events as affecting more than one-third of a chromosomal arm and calculated the percentage of genome alteration (PGA) as a general indicator of genomic aberrations. Considering a calculated mean PGA value of 4.8% across 208 meningiomas, we classied any meningioma with a higher or lower PGA value as CNV-high or CNV-low, respectively (Supplementary Fig. 5).

Overall, primary atypical NF2 tumours were signicantly more likely to be CNV-high, as compared with benign NF2 ones (P 0.02, Students t-test, Fig. 1c). On the basis of correction for

mutational background by considering only the NF2 mutant meningiomas, we calculated atypical tumours to have a 2.19 times greater risk to be CNV-high (P 0.001, Fishers exact test)

(Supplementary Data 3b). When we performed the same analysis considering non-NF2 mutant meningiomas (n 96 versus 12,

NATURE COMMUNICATIONS | 8:14433 | DOI: 10.1038/ncomms14433 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14433

benign versus atypical, respectively), we did not nd any statistically signicant difference in PGA, likely due to the small number of non-NF2 mutant meningiomas (P 0.27, Students

t-test, Fig. 1c). Not surprisingly, when we combined the benign and atypical samples together and simply compared non-NF2 mutants to NF2 mutant meningiomas (among which percentage of atypical tumours is much higher), we obtained a similar result, with NF2 mutants being genomically more unstable. These ndings demonstrate that atypical tumours with NF2 loss display increased chromosomal instability compared with their benign counterparts, while those without NF2 loss do not.

We next compared specic CNV events associated with atypical samples, initially considering the NF2 mutants only. As expected, all benign as well as atypical NF2 mutant samples demonstrated chromosome 22 loss. There were, however, signicant differences with respect to several other chromosomes. We found that chromosome 14q loss was the most signicant difference between atypical and benign NF2 meningiomas (18% versus 7%, respectively; P 0.0016 false discovery rate

(FDR) 0.04, Fishers exact test). This was followed by whole

arm losses of chromosomes 10q, 10p, 1p and 6q, (P 0.002

FDR 0.04, P 0.005 FDR 0.06, P 0.008 FDR 0.07,

P 0.013 FDR 0.1 respectively, Fishers exact test) (Fig. 1d).

No CNV events were signicantly associated with atypical non-NF2 meningiomas compared with benign non-NF2 samples.

We next identied tumour suppressor genes located within these genomic regions frequently deleted in atypical NF2 meningiomas. Deletions in phosphatase and tensin homologue (PTEN) (10q23) (13% versus 0% in atypical versus benign NF2 tumours, FDR 0.018), myc-associated factor X (MAX) (14q23)

(34% versus 10%, FDR 0.03), neurogenic locus notch homolog

protein 2 (NOTCH2) (1p13) (20% versus %3, FDR 0.043),

AT-rich interactive domain-containing protein 1B (ARID1B) (6q25) (30% versus 10%, FDR 0.058), a member of the

SWI/SNF-A chromatin-remodeling complex and cyclin-dependent kinase inhibitor 2C (CDKN2C) (1p32) (46% versus 24%, FDR 0.08) genes were signicantly more common in

atypical NF2 samples (Fig. 1e).

When we corrected both for the mutational background (NF2 versus non-NF2 mutant) and CNV status (CNV-high versus-low), we again did not identify any differences in the overall number of damaging somatic mutations, further suggesting that changes in coding variations did not signicantly contribute to the risk of being atypical.

mRNA expression signature separates atypicals from benigns. After investigating genomic alterations, we expanded our analysis to understand the transcriptional changes underlying atypical meningiomas. We studied the messenger RNA (mRNA) expression proles of 138 primary meningiomas, again focusing on the comparison between atypical versus benign samples (n 26 versus 112, respectively) (Supplementary Data 4a).

Principal component analysis, unsupervised hierarchical and consensus clustering clearly distinguished NF2 CNV-low, NF2 CNV-high and non-NF2 samples, however failed to completely separate atypical versus benign tumours (Fig. 2a). This result suggested that gene expression correlated more closely with the underlying driver mutation rather than the histological grade.

To characterize specic transcriptional changes associated with atypical tumours, we next corrected for the underlying meningioma driver mutation by dividing the cohort into NF2 and non-NF2 subgroups as before. We identied mRNA signature genes (n 483) that were differentially expressed between

atypical and benign samples when stratied by NF2 status (Supplementary Data 4b). Principal component (PC) analysis of

gene expression data using these signature genes correctly separated atypical samples from benign ones, even when NF2 status was not considered (Fig. 2b; Supplementary Figs 69). This analysis demonstrated that atypical tumours harbour a distinct expression prole of a set of transcripts, which were signicantly enriched for cell cycle processes, including upregulation of the E2F and FOXM1 transcription factor networks (hypergeometric test) (Fig. 2c,d; Supplementary Fig. 10; Supplementary Data 4c,d). Importantly, expression levels of the genes associated with atypical samples were not affected by the CNV events, suggesting that transcriptional regulation, and not genomic instability, was the primary driver of these processes.

We next used the top 25 most differentially expressed genes to build a random forest prediction model, aiming to associate gene expression with meningioma histological grade. Our prediction model had a 96% prediction accuracy (4% out of bag error rate) on this training set. When we used the dataset of an independent meningioma gene expression study13 as a validation set, we obtained a 91% accuracy rate, after considering only histologically atypical meningiomas with high or medium Ki-67 index. Interestingly, 15 of the top 25 most signicant genes were involved in cell cycle processes, including E2F, ASPM, aurora kinase, cyclin and centromic proteins (Fig. 2c; Supplementary Data 4b). Other top upregulated genes included POLQ, a DNA damage repair gene associated with high tumour grade and genomic instability in breast cancer14; RET, a receptor tyrosine kinase oncogene; and BCL2, which has been found to have increased expression in many types of cancer such as lymphoma, small cell lung and prostate cancer1517. Importantly, atypical meningiomas also showed increased expression of the EZH2 gene. This gene codes for the catalytic subunit of the PRC2 complex, known to play a key role in both tissue-specic stem cell maintenance and tumour development (Fig. 2d).

microRNA regulatory networks in atypical meningiomas. We next compared miRNA expression patterns in benign versus and atypical meningiomas (n 15 versus 17, respectively)

(Supplementary Data 5a). Unsupervised hierarchical clustering of miRNA expression proles correlated well with being benign versus atypical, CNV status, and also with the underlying meningioma driver mutation, clearly separating NF2 mutants from non-NF2 meningiomas (Fig. 3a).

On the basis of this observation, we searched for differentially expressed miRNAs between atypical and benign meningiomas (Fig. 3b, Supplementary Data 5b). Differential expression of miRNAs has been shown to play important roles in the control of cancer hallmark functions such as invasion, metastasis, proliferation and apoptosis18 (Fig. 3c). In atypical samples, we identied 67 differentially regulated miRNAs (54 downregulated, 13 upregulated). We next correlated the differentially expressed miRNAs with large-scale chromosomal events. We identied miRNAs clustering on a deleted region on chromosome 14q32 (n 28) to be downregulated, suggesting the expression of these

miRNAs were largely driven by CNVs.

To correlate downregulated miRNA clusters with gene expression patterns, we inferred miRNA:mRNA regulatory networks using samples that contained both datasets (n 22).

We calculated the observed number of negatively correlated targets for each miRNA, with the assumption that if the miRNA expression were playing a signicant biological role, it is expected that the expression levels of its target genes would be altered. Consistent with the expected biological action of elevated miRNA expression, we identied the number of negatively correlated targets to be signicantly higher than the positively correlated ones (Po2.2e 16, paired Wilcoxon test) (Supplementary

4 NATURE COMMUNICATIONS | 8:14433 | DOI: 10.1038/ncomms14433 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14433 ARTICLE

Atypical/benign

Atypical

Benign

Mutation typeNF2

Non-NF2

Atypical/benign

Mutation type

CNV status

CNV-high CNV-low Unknown

CNV status

Consensus

clustering

NF2 CNV-low

NF2 CNV-high

Non-NF2

b c

Mitotic prometaphase

Pathway GO

FOXM1 transcription factor network

Mitotic G2G2/M phases

EPHEphrin signalling

PC1 (33%)

E2F transcription factor network

Mitotic cell cycle process

Mitotic cell cycle

0 3 6 9log10 FDR

Atypical

Benign

PC2 (10%)

EZH2 (ILMN_1708105)

E2F2(ILMN_1777233)

NF2 non-NF2

6.0

Expression value

5.5

PC3 (8%)

Atypical

Benign

5.0

4.5

NF2

Non-NF2

NF2

Non-NF2

NF2

Non-NF2

NF2

Non-NF2

Figure 2 | mRNA expression prole in atypical meningiomas. (a) Unsupervised hierarchical clustering of 138 meningiomas by genome-wide expression proling is shown. Atypical versus benign histology, underlying meningioma driver mutations, copy number variations, which are colour coded, are shown on the left. Although the expression prole accurately clusters meningioma samples based on driver mutations, it does not fully differentiate atypical versus benign tumours. Consensus clustering also separates non-NF2 meningiomas (Hedgehog, POLR2A, TRAF7/PI3K, TRAF7/KLF4) into different subgroups. (b) PC analysis using mRNA signature genes separates atypical and benign samples. (c) GO Term and pathway enrichment for genes associated with the atypical phenotype compared with benign are shown (atypical n 26, benign n 112). Black line indicates an FDR of 0.05. (d) EZH2 and E2F2 gene

expressions, which are upregulated in atypical versus benign samples, are plotted (empirical Bayesian method). Lines depict the median values; boxes plot 25th to 75th percentiles, whereas separately plotted dots show the outliers.

Fig. 11). We found that loss of miRNAs on the frequently deleted chromosomal segment 14q32 was signicantly inversely correlated with a set of their known target genes (n 85)

(Supplementary Data 5c). These upregulated genes were enriched for receptor tyrosine kinase signalling pathways as well as the GO term regulation of programmed cell death (Fig. 3d; Supplementary Data 5d).

We also identied 5 members of let-7 family (let-7a, chr22, log FC 0.66; let-7b, chr22, log FC 0.8; let-7c, chr21, log

FC 0.6; let-7d, chr9, log FC 0.78; let-7e, chr19,

log FC 0.86) to have aberrant expression in atypical

meningiomas. This family of miRNAs has been shown to be deregulated in various cancers such as prostate cancer and neuroblastoma19,20, as well as malignant schwannoma and meningioma21, potentially through negative regulation of EZH2 (refs 19,20,22). Consistent with these observations, the down regulation of let-7 family in atypical meningiomas was correlated with upregulation of EZH2 mRNA expression in our dataset (let-7c: B 0.65, P 0.001, let-7d: B 0.62, P 0.002, let-7e:

B 0.49, P 0.02 correlation test) (Fig. 3e). Importantly,

expression levels of these miRNAs, other than let-7b (deleted in 75% versus 47% in atypical versus benign meningiomas, respectively) were not affected by the CNV events.

DNA methylation patterns in atypical meningioma subgroups. We next focused on epigenetic alterations, again comparing atypical meningiomas to benign ones. We initially analyzed the DNA methylation status of 60 samples, including meningiomas (n 46 versus 11, benign versus atypical, respectively) as well as

control tissues obtained from the normal meninges (n 3)

(Supplementary Data 6a,b). Unsupervised analysis of whole-genome methylation data using two different clustering methods, namely principal components and consensus clustering, revealed consistent results and identied distinct subgroups of meningiomas (Fig. 4a,b). Neither of these clusterings was affected by large-scale chromosomal events, since when we removed the sites in the regions that are affected by large-scale chromosomal events (chr1, chr14, chr22), clustering results did not change (Supplementary Fig. 12).

Integration of the methylation data with copy number and mutational analyses revealed the presence of distinct molecular signatures in atypical meningiomas. First, genome-wide DNA methylation patterns clearly separated the NF2 mutant meningiomas from relatively less methylated non-NF2 tumours, with atypical non-NF2 meningiomas clustering distinctly from benign non-NF2 ones (Fig. 4a,b). Second, for NF2 mutant meningiomas, which also formed 2 different sub-clusters, one of

NATURE COMMUNICATIONS | 8:14433 | DOI: 10.1038/ncomms14433 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14433

a b

CNV status

Mutation type

NF2

Non-NF2

Atypical/benign

Atypical

Benign

CNV-high

CNV-low

Cancer related pathways

let-7a miR-125a miR-125b miR-128 miR-204 miR-21 miR-25 miR-30d miR-483 miR-944

let-7b let-7c let-7d let-7e miR-126 miR-127 miR-130a miR-134 miR-136 miR-143 miR-145 miR-154 miR-23b miR-27b miR-320a miR-320c

miR-326 miR-329 miR-337 miR-33a miR-345 miR-376c miR-377 miR-379 miR-382 miR-411 miR-432 miR-495 miR-502 miR-574 miR-655 miR-656

Mutation type

NF2

Non-NF2

Atypical/benign

Atypical

Benign

Low

High

CNV status

CNV-high

CNV-low

r2 = 0.65

r2 = 0.62

let-7d expression

r2 = 0.49

5.8

Down regulated miRNA

EZH2(ILMN_1708105)

expression

Up regulated miRNA

5.4

5.0

Fibroblast growth factor receptor signalling pathway

Epidermal growth factor receptor signalling pathway

4.6

ERBB signalling pathway

Neurotrophin TRK receptor signalling pathway

Negative regulation of programmed cell death

0 1 2 7 8 9

log10 FDR

8.5 9 9.5 10 10.5 11

let-7c expression

6 7 8

let-7e expression

Figure 3 | miRNA expression prole in atypical meningiomas. (a) Unsupervised hierarchical clustering of genome-wide miRNA expression of 32 meningiomas. Integrated view of the miRNA expression clustering, combined with the underlying meningioma driver gene mutation, histological grade and CNV-high or -low status, which are all colour coded, are shown. (b) Heatmap visualization of differentially expressed miRNAs between atypical versus benign samples. The colour scale for log2 (fold change) is shown at the bottom. MiRNAs affected by CNVs are marked by an asterisk. (c) Differentially expressed miRNAs with either inhibitor or activator effects on cancer hallmark functions are shown. Red coloured miRNAs are upregulated whereas blue coloured miRNAs are downregulated in atypical samples. (d) The top gene ontology terms for the genes regulated by the 14q32 miRNA cluster. Black line indicates an FDR value of 0.05. (e) Expression of let-7 family, a known inhibitor of EZH2, negatively correlates with EZH2 mRNA expression (Spearman correlation test). Lines depict the median values; boxes plot 25th to 75th percentiles, whereas separately plotted dots show the outliers.

the subgroups revealed a distinct hypermethylated phenotype and was signicantly enriched for atypical CNV-high samples (P 0.03 (two NF2 cluster comparison), Fishers exact test)

(Fig. 4ac; Supplementary Data 6b). The majority of samples in this cluster revealed large-scale chromosomal aberrations including chromosome 1p loss, which was found in all samples (Fig. 4a; Supplementary Data 6b). Indeed, we observed a statistically signicant positive correlation between the degree of chromosomal alterations and the amount of genome-wide DNA hypermethylation (P 0.01 correlation test) (Fig. 4d).

Interestingly, we found that atypical samples harbouring both NF2 and SMARCB1 mutations, although genomically stable, still clustered closer to the CNV-high atypical NF2 samples in

principal component analysis and were also hypermethylated. These results suggest the presence of two distinct pathways underlying NF2 mutant, atypical meningiomas: one through acquiring genomic instability and the other through recurrent SMARCB1 mutations (Fig. 4b,c).

Hypermethylation of PRC2 binding sites in atypical tumours. Given that NF2 mutant atypical meningiomas displayed a hypermethylated phenotype, we next investigated which gene sets were more methylated in atypical versus benign meningiomas (GREAT tool), making use of the Molecular Signature Database (MSigDB)23. We found that Polycomb Repressive Complex 2

6 NATURE COMMUNICATIONS | 8:14433 | DOI: 10.1038/ncomms14433 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14433 ARTICLE

a b

Consensus cluster

CNV status

Cluster 1

Cluster 2

Cluster 3

Mutation type

Atypical/benign

PC1

Consensus cluster

Atypical/benign

Atypical

Benign

Mutation type

NF2

Non-NF2 10

CNV status

CNV-high

CNV-low

Cluster 4

10 5 0 5

PC2

Hyper Hypo

Atypical

Benign

% Hyper/hypo

methylated

sites compared to

benign non-NF2

%80

%60

CNV-high

CNV-low NF2/SMARCB1

NF2

Non-NF2

%40

%20

25%

% of genome alteration

No. differentially hypermethylated

sites / total no. of differentially

NF2 CNV-high

SMARCB1 Atypical

NF2

SMARCB1 CNV-low

Benign

NF2 mutant samples

Atypical Non-NF2

0.76

0.64

0.47

methylated sites (compared to

benign Non-NF2)

20%

15%

10%

Atypical vs benign

Benign non-NF2 vs control

H3K27me3 targets in ES cells

NF2/SMARCB1

NF2 CNV-highvs benign non-NF2

NF2 CNV-lowvs benign non-NF2

PRC2 targets in ES cells

NF2CNV-high

(Cluster 1)

NF2CNV-low

(Cluster 2)

Non-NF2 Atypical

(Cluster 3)

Atypical non-NF2 vs benign non-NF2

NF2/SMARCB1vs benign non-NF2

0 50 100 150 200

log10 (adjusted P-value)

Figure 4 | DNA-methylation patterns across benign versus atypical meningiomas. (a) Integrated view of DNA methylation clustering combined with the underlying meningioma driver gene mutations, histological grades and genome copy numbers are shown (atypical n 11, benign n 46). Consensus

clustering denes four main subgroups based on methylation status. The CNV-high or -low categories, as well as NF2 and SMARCB1 mutations are marked. (b) PC analysis of meningioma methylation data, which separate NF2 CNV-high, -low and non-NF2 samples into different groups, is shown (atypical n 11,

benign n 46). (c) Percentages of hyper- or hypo-methylated sites among differentially methylated sites in various subgroups, as compared with benign

non-NF2 samples are plotted. All atypical as well as benign CNV-high NF2 mutant meningiomas are hypermethylated when compared with benign non-NF2 samples (atypical n 11, benign n 46). (d) Percentage of altered genome signicantly correlates with the extent of genome-wide DNA hypermethylation

(atypical n 11, benign n 46). Along the horizontal axis, the meningioma samples are grouped according to the consensus cluster numbers as in a.

The number of differentially hypermethylated sites divided by the total number of sites is shown as the green dashed line and quantied along the vertical axis on the right. Lines depict the median values; boxes plot 25th to 75th percentiles, whereas separately plotted dots show the outliers. (e) Comparisons for atypical versus benign as well as various other molecular subgroups, which are colour coded, are shown for Molecular Signatures Database (MSigDB) gene set enrichment summary plots. The most signicant gene sets are plotted.

(PRC2) target sites in human embryonic stem cells (PRC2-hESCs)24, as well as Homeobox domain sites25, represented the most differentially methylated regions (FDR 1.15 10 71 and

6.6 10 55, respectively) (Fig. 4e; Supplementary Data 7a).

Interestingly, we also found enriched methylation of these sites across benign meningiomas compared with control meninges, albeit less than that observed in atypical samples (FDR 5.8 10 17). This nding suggests that the degree of

methylation at PRC2-hESCs and Homeobox domain sites may establish a spectrum of meningioma severity, with normal meninges, benign meningioma, and atypical meningioma having progressive amounts of increased DNA methylation at these sites.

To further explore this nding, we next considered the methylation status of PRC2-hESC targets among the various

subgroups identied in our genomic, gene expression, and methylation analyses. We observed the highest enrichment for PRC2-hESCs hypermethylation in the atypical and benign NF2 CNV-high samples (FDR 9.3 10 142 when compared

with benign non-NF2 mutant samples), followed, in order, by the NF2/SMARCB1, NF2 CNV-low and atypical non-NF2 meningiomas (FDR values 1.2 10 91, 4.1 10 77 and

2.8 10 75 respectively) (Fig. 4e; Supplementary Data 7bf).

It has been previously shown that EZH2, the catalytic subunit of PRC2, is a recruitment platform for DNA methyltransferases (DNMTs), acting as a direct controller of DNA methylation at PRC2 binding sites26. On the basis of this knowledge and our observations, transcriptional upregulation of EZH2 and down regulation of its regulator let-7 in atypical samples suggest that increased methylation in atypical tumours might be related to

NATURE COMMUNICATIONS | 8:14433 | DOI: 10.1038/ncomms14433 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14433

deregulated PRC2 activity. However, this association was not experimentally tested in our study. (Supplementary Figs 13 and 14; Supplementary Data 8a,b).

H3K27me3 ChIP-seq conrms silencing of PRC2 targets. In addition to DNA methylation, the PRC2 complex also plays an important role in regulation of histone state. In particular, EZH2 is involved in trimethylation of histone 3 lysine 27 (H3K27me3) to drive long-term silencing of gene expression during development and other processes27. To conrm increased activity of EZH2 in atypical tumours, we next performed H3K27me3 ChIP-seq using both atypical and benign meningioma samples (n 3 each) (Supplementary Data 9a). H3K27me3 proles clearly

separated atypical meningiomas from benign ones (Fig. 5a). Interestingly, differential binding analysis between atypical and benign samples showed an overall increase in H3K27me3 binding in atypical samples (Fig. 5b; Supplementary Data 9a).

GO Term enrichment analysis of H3K27me3 differentially bound and underexpressed genes in atypical as compared with benign samples revealed neuron projection morphogenesis, neuron differentiation, and neurogenesis GO terms, suggesting that atypical meningiomas, while turning on embryonically active pathways, repress pathways involved in differentiation (Fig. 5c, Supplementary Data 9b). On the basis of GREAT enrichment analysis, H3K27me3 differentially bound and underexpressed genes in atypicals as compared with benign ones were enriched for PRC2-hESC targets (FDR 0.03). Since GREAT analysis

assesses the enrichment using hypergeometric test, the enrichment was signicant even though there was an overall increase in H3K27me3 signal in atypicals.

H3K27ac ChIP-seq identies activated regions in atypicals. Given the role of epigenetic silencing in atypical samples, we investigated if activating epigenetic alterations may also be associated with these tumours. To investigate this question, we performed ChIP-seq targeting histone 3, lysine 27 acetylation (H3K27ac) using 3 atypical tumours and 15 benign tumours, as well as control meningeal tissues (2 dura samples) (Supplementary Data 10a). We identied transcriptionally active regions by focusing on broad genomic loci that harboured dense clustering of H3K27ac signals, previously described as super-enhancer regions. These features have been shown to be important in determining not only cell identity, but also cancer cell properties28. Indeed, super-enhancer driven expression of particular oncogenes has been shown to be fundamental in formation of specic tumours28. Although active genomic regions as dened by H3K27ac binding were highly correlated among different meningioma samples (with a minimum correlation coefcient of 0.7), we were able to clearly classify meningiomas into atypical versus benign samples as well as into various molecular subgroups, including NF2 CNV-high, NF2 CNV-low and non-NF2 samples (Fig. 5d).

We identied 19 super-enhancers with concordant changes in gene expression between atypical and benign samples (Fig. 5e, Supplementary Data 10b). Notably, in the atypical group we identied a differentially active super-enhancer near the ZIC Family Member 1 (ZIC1) transcription factor, a regulator of neural crest differentiation29, that is shown to play an essential role in the proliferation of meningeal cell progenitors30. This gene showed increased H3K27ac binding in atypical versus benign meningiomas (FDR 0.007, fold change log FC 4.95) with

increased expression in atypical NF2 samples (gene expression FDR 0.001) (Fig. 5f,g).

To characterize super-enhancers associated with atypical tumours independent of the mutational background, we next

compared atypical NF2 mutant samples to benign NF2 mutants and identied ve super-enhancers with concordant changes in gene expression. We found decreased H3K27ac signal near the GDNF family receptor alpha-1 (GFRA1) gene that is involved in neuron survival and differentiation31. Importantly, decreased H3K27ac binding was associated with decreased GFRA1 expression in atypical NF2 versus benign NF2 meningiomas (FDR 6.2 10 5 fold change log FC 9.54).

Benign meningiomas undergoing atypical progression. We next expanded our analysis of primary atypical meningioma to include recurrent (or progressed) atypical tumours. On the basis of our ndings as well as previously published studies11, we hypothesized that distinct molecular mechanisms may be involved in formation of de novo atypical meningioma versus progression of benign tumours. Our analysis focused on a small cohort of paired meningioma samples that included the original benign sample as well as the progressed atypical counterpart from the same patient (n 4). On the basis of exome sequencing, we

identied all samples to be NF2 mutant (Supplementary Data 11).

Consistent with the previously published literature32, we found Telomerase Reverse Transcriptase (TERT) gene promoter mutations (C228T/-124G-A chr5:1,295,228; C250T/-146G-A chr5:1,295,250) in 2 of the 4 progressed atypical samples. A larger screen of 27 recurrent atypical meningiomas (non-paired) identied 2 additional samples, suggesting approximately 13% (4/31) of progressed atypical samples harbour mutations in the TERT promoter (Supplementary Data 11). Importantly, when we screened our cohort of 110 primary, non-recurrent meningiomas (n 66 NF2 atypical, n 12 NF2 benign CNV-high, n 4

non-NF2 benign, CNV-high, n 12 NF2 benign CNV-low,

n 12 non-NF2 benign CNV-low and n 4 NF2/SMARCB1),

we did not identify any TERT promoter mutations (Supplementary Data 11). This nding further supports the presence of two distinct molecular pathways that underlie formation of atypical meningiomas: either through a de novo pathway as detailed in this study, or due to transformation of benign meningiomas, partially through acquisition of an activating TERT promoter mutation.

DiscussionAmong the three pathological grades of meningiomas, grade II atypical tumours show increased mitotic activity as well as high cellularity, small cells with a high nuclear to cytoplasmic ratio, prominent nucleoli, uninterrupted patternless or sheet-like growth and occasional foci of spontaneous necrosis2,3. The 2007 revision of the WHO system classied meningiomas with brain invasion also as grade II, even if they were histologically benign, as the presence of brain invasion worsened the clinical outcome and often necessitating the need for adjuvant radiotherapy. Thus, the WHO grading system correlates well with clinical parameters including the likelihood and time to recurrence after initial treatment (that is, surgery) as well as overall survival, which worsens with increasing grade. Indeed, atypical meningiomas are associated with up to a 40% recurrence rate at 5 years following total resection4. Besides dictating the clinical course, pathological grading of meningiomas further guides treatment such that surgical resection is the standard rst-line treatment for all accessible, symptomatic meningiomas, with radiation being reserved for the post-operative management of higher grade II and III tumours. There are no standard chemotherapy options.

Unlike grade I benign meningiomas, the genomic landscape of atypical meningiomas is not well understood. Here we studied the genomic landscape of histologically atypical, primary grade II

8 NATURE COMMUNICATIONS | 8:14433 | DOI: 10.1038/ncomms14433 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14433 ARTICLE

meningiomas only, excluding those tumours that progressed from a benign, recurrent meningioma or those that were classied as grade II based solely on brain invasion. Utilizing comprehensive next-generation genomic approaches, including exome sequencing, mRNA and miRNA expression proling, as well as H3K27 acetylation and trimethylation ChIP-seq and DNA methylation analyses, we identied the molecular features associated with

atypical samples independent of their mutational background and CNVs. We found that primary atypical meningiomas were enriched for NF2 loss. On the basis of exome sequencing, we did not observe any signicant difference in the number of somatic protein altering mutations in atypical versus benign samples (even when stratied by NF2 status), and also did not nd novel recurrent driver genes.

a b c

Increased H3K27me3

Organ development

0 2 4

Atypical vs benign (-log 10FDR)

System development

Neuron projection morphogenesis

Anatomical structure morphogenesis

Cell projection morphogenesis

Cell morphogenesis involved in neuron differentiation

Nervous system development

hESC PRC2 targets

Neurogenesis

log10 FDR

Atypical

Benign

5 0 5 10

0.85

0.9

0.95

Atypical vs benign H3K27me3 log2 fold change

d e

Atypical/benign

Atypical

Benign

Mutation type

NF2

Non-NF2

Control

Atypical vs benign gene expression

log 2fold change

EPB41L4A

IL1RAP

KCNJ8

PTPRC

PSD3

ZIC1

Control

PLXNA4

C2orf71

ADAMTS2

EPHA2

FGD5

ZBTB7C

RSG1

LPHN2

PLEKHH2

LDLRAD3

LMO4

PTHLH

CNV status

CNV-high

CNV-low

TIMP3

5 0 5

Atypical vs benign H3K27ac ChIP-Seq log2 fold change

4.2

0.6 0.8 1

f g

ZIC1(ILMN_3213167)

100 kb hg19

147,100,000 147,150,000 147,200,000 147,250,000 147,300,000 147,350,000

FDR=0.001

Atypical

Benign

Super enhancer

4.4

Expression value

rpm/bp

NF2 atypical

NF2 benign

4.3

Non-NF2 benign

Non-NF2

ZIC1

NF2

NATURE COMMUNICATIONS | 8:14433 | DOI: 10.1038/ncomms14433 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 9

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14433

However, we did nd signicant differences in the extent of chromosomal instability, with genomically unstable samples having the greatest chance of being NF2 mutant atypical. The CNV status associated signicantly with levels of DNA methylation, such that CNV-high NF2 mutant meningiomas displayed a hypermethylated phenotype, particularly affecting the PRC2 targets in embryonic stem cells. These sites also revealed increased H3K27me3 signals based on ChIP-seq. These observations raise the possibility that differentiation-related genes may be suppressed in atypical meningiomas, similar to PRC2-mediated inhibition of these genes during embryonic development33. Indeed, silencing of PRC2 targets has been shown to play a role in the formation and maintenance of other forms of cancer by locking the cells into stem-like cellular state11,34. We also identied the Homeobox genes, which similar to PRC2 signalling play a role in embryogenesis and differentiation35, to be hypermethylated in NF2 CNV-high meningiomas.

Consistent with these ndings, we showed transcriptional upregulation of the catalytic subgroup of the PRC2 complex, EZH2, in atypical meningiomas as compared with benign ones. Besides functioning as a histone methyltransferase, EZH2 has been shown to be a recruitment platform for DNA methyltransferases (DNMTs), acting as a direct controller of DNA methylation at PRC2 binding sites36. Previous reports have also identied this enzyme as a marker for aggressiveness, particularly in glioblastoma and renal cell carcinoma37,38. The association of EZH2 expression with atypical meningiomas suggests a role for EZH2 as a marker in higher grade meningiomas.

A potential mechanism for increased EZH2 expression is the loss of the miRNA let-7, a known negative regulator of EZH2, in atypical samples (let-7c: FDR 0.03, let-7d: FDR 0.01, let-7e:

FDR 0.01). Indeed, samples with let-7 loss revealed increased

EZH2 mRNA expression (let-7c: B 0.65, P 0.001, let-7d:

B 0.62, P 0.002, let-7e: B 0.49, P 0.02, correlation test)

in atypical compared with benign samples. This nding suggests that loss of let-7 contributes to increased EZH2 expression and deregulated PRC2 activity in atypical meningiomas. We also identied miRNA cluster at 14q32 to be downregulated in atypical samples, affecting genes in receptor tyrosine kinase signalling pathways. Deregulation of miRNA cluster at 14q32 has been reported in various cancer types39.

In a distinct subgroup of CNV-low but hypermethylated NF2 mutant atypical meningiomas, we observed recurrent SMARCB1 mutations. Interestingly, loss of the tumour suppressor SMARCB1 has previously been reported in various malignant rhabdoid tumours40, as well as in families with multiple meningiomas with schwannomatosis41. Our group has also recently described recurrent SMARCB1 mutations in a subset of benign

meningioma7, however enrichment of this alteration in a distinct molecular subgroup of atypical samples has not previously been reported.

Taken together, we report a unique set of genomic and epigenomic events that distinguish primary atypical tumours from other types of meningiomas. Similar to the molecular pathways underlying formation of malignant gliomas, we propose that distinct molecular mechanisms underlie the formation of de novo (primary) and progressed (recurrent) atypical meningiomas. While previous groups have reported TERT promoter mutations in atypical samples, we nd that these

Chromosomal instability14q, 10q, 10p, 1p and 6q deletions

SMARCB1p.(Arg383Gln)/p.Arg386His

Proliferation upregulation of E2F2 and FOXM1 cell cycle pathways

promoter mutations

De novo atypical pathway

Figure 6 | De novo pathway in the formation of atypical meningiomas.

Samples with NF2 loss were highly enriched among primary (de novo) atypical cases, which co-occur either with genomic instability or recurrent mutations in SMARCB1. These tumours have hypermethylation of the PRC2 binding sites in hESCs, thereby phenocopying a more primitive cellular state. Epigenetic changes acting on the PRC2 networks are associated with upregulation of the catalytic subunit EZH2, and down regulation of its regulator let-7. Primary atypical tumours harbour increased expression of cell-cycle related genes, including the E2F and FOXM1 transcription factor networks. TERTpromoter mutations are not present in de novo tumours, and limited to recurrent samples only.

Figure 5 | Histone modications in meningiomas. (a) Unsupervised hierarchical clustering of meningioma samples using H3K27me3 proles is shown (Pearson correlation) (atypical n 3, benign n 3). Grades are colour coded. (b) Starburst plot for atypical versus benign H3K27me3 ChIP-seq signal fold

change (horizontal axis) and FDR (vertical axis) is shown. Red and blue circles indicate regions that have increased binding in atypical versus benign meningiomas, respectively. PRC2-hESC targets are marked in red colour. These analyses reveal increased H3K27me3 binding in atypical samples, which shows particular enrichment of PRC2-hESC targets. (c) GO Term enrichment plot for genes with increased H3K27me3 binding and decreased expression is shown (atypical n 3, benign n 3). Black line indicates an FDR of 0.05 (d) Unsupervised hierarchical clustering of 18 meningiomas and 2 controls based

on H3K27ac prole in super-enhancer sites is shown (Pearson correlation). Super-enhancer proles are highly correlated among different meningioma subtypes. All samples are colour coded, which are shown on the right. (e) H3K27ac ChIP-seq results correlate with gene expression proles (n 18).

Atypical versus benign H3K27ac ChIP-seq signal fold change is plotted along the horizontal axis, whereas atypical versus benign gene expression fold change is shown along the vertical axis. Red dots indicate genes that have both activated super-enhancers and are overexpressed in atypical versus benign meningiomas, whereas blue points indicate genes that are associated with de-activated enhancers and are underexpressed. An FDR threshold of 0.05 is used both for gene expression and H3K27ac ChIP-seq data (empirical bayes and DiffBind methods). (f) H3K27ac ChIP-seq occupancy at a super-enhancer near ZIC1 (n 18). The horizontal axis shows genomic position whereas the vertical axis shows signal of ChIP-seq occupancy in units of reads per million

(r.p.m.). Super-enhancer region that is differentially bound in atypical samples is depicted as a black line over the gene track and highlighted. (g) ZIC1 gene expression across atypical and benign samples is plotted (n 138) (empirical bayes method). Lines depict the median values; boxes plot 25th to 75th

percentiles, whereas separately plotted dots show the outliers.

10 NATURE COMMUNICATIONS | 8:14433 | DOI: 10.1038/ncomms14433 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14433 ARTICLE

events were limited to recurrent samples only, and not present in de novo tumours. Our results show that primary atypical lesions arise in response to synergistic signals that stimulate proliferation (through upregulation of the E2F2 and FOXM1 regulated cell-cycle pathways), as well as epigenetic changes (through EZH2/PRC2 target activation), that have previously been shown to phenocopy more primitive cellular states34 (Fig. 6).

In conclusion, we report comprehensive genomic, transcriptomic and epigenomic analyses of primary atypical meningiomas, identifying distinct molecular pathways that dene the landscape of these tumours and separate them from recurrent atypical samples. Our results have biological as well as clinical implications through enhanced molecular classication and therapeutic intervention in atypical meningiomas.

Methods

Clinical materials. Institutional Review Board approvals for genetic studies, along with written consent from all study subjects, were obtained at the participating institutions. The specic approval committees included the Human Investigation Committee at Yale University, Ethikkommission der Medizinischen Fakultat der Universitat zu Kln and Ethikkommission an der Medizinischen Fakultat der Rheinischen Friedrich-Wilhelms-Universitat Bonn.

Selective tissue dissection. For each frozen specimen submitted for whole-exome sequencing, sections were re-reviewed to conrm the diagnosis and assess the adequacy of the frozen tissue for experimental analysis. On H&E-stained sections from frozen tissue blocks, areas of interest were identied and microscopically dissected to ensure that each sample consisted of 470% tumour cells;

unwanted regions such as inammatory and necrotic areas were excluded. Tumours in the replication cohorts did not undergo selective tissue dissection. DNA/RNA was prepared using the Allprep DNA/RNA Mini Kit (Qiagen) with the assistance of a QIAcube.

Exome capture and sequencing. Nimblegen/Roche human solution-capture exome array (Roche Nimblegen, Inc.) was used to capture the exomes of blood and tumour samples according to the manufacturers protocol. Sequencing of the library was performed on Illumina HiSeq instruments using paired-end 74 basepair reads by multiplexing two tumour samples or three blood samples per lane. Image analysis and base calling was performed by Illumina Pipeline with default parameters, installed on Yale Universitys High Performance Computing Cluster.

Whole-exome sequence analysis. We rst ltered the reads based on Illumina quality score. The low-quality 30-end of the reads using FASTX-Toolkit (http://hannonlab.cshl.edu/fastx_toolkit/index.html

Web End =http://hannonlab.cshl.edu/fastx_toolkit/index.html) and PCR primer-contaminated sequences, that are considered to lead to alignment artifacts, were trimmed using cutadapt. The sequences were kept only if both reads in a pair had more than35 bases remaining after the above trimming and ltering quality measures.

The reads were aligned to the human reference genome (version GRCh37) using Stampy (version 1.0.21) in a hybrid mode with BWA (version 0.5.9-r16)42,43. MarkDuplicates algorithm from Picard was used for agging PCR duplicates. Alignment quality metrics were calculated using CollectAlignmentSummaryMetrics and CalculateHsMetrics utilities of Picard (http://picard.sourceforge.net/

Web End =http://picard.sourceforge.net/). Mean target coverage was 246 and 160 for tumour and blood respectively. The average percentage of reads with at least 20 coverage was 93% and 90% for tumour and

blood respectively. We performed multi-sequence local realignment around putative and known insertion/deletion sites. This was followed by the base quality score recalibration using the Genome Analysis Toolkit (GATK, version 2.52)44. We detected variant sites (point mutations and small indels) for tumour and matched blood pairs using the HaplotypeCaller algorithm from GATK. The genotype likelihood-based somatic score, which was proposed by Li, was used45. We ltered out the variants according to the classes of genotype changes in tumour with respect to the blood46. We used various quality metrics to lter out variants:(1) somatic score less than 20, (2) overlapping a RepeatMasker or segmental duplication annotated region, (3) low quality (o30) and low quality-by-depth values (o1) (4) high mapping quality zero reads, (5) strand bias, (6) mutation cluster of size 42 (7) homopolymer runs of length 4 10 base pairs within / 5 base pairs around the mutation or from the right of the mutation or

(8) ClippingRankSum (calculated by GATK)o 3 or 43. We also excluded the

sites that have more than 1% frequency in the NHLBI Exome Variant Server Database (http://evs.gs.washington.edu/EVS/

Web End =http://evs.gs.washington.edu/EVS/) and 1000 Genome Database. In addition we used our internal database of 2216 exomes to compare the variant allele frequencies of each gene and excluded the variants in the genes that have greater than 150 variant alleles. Finally, we annotated variant alleles using Ensembl database (version 69) with the help of Variant Effect Predictor (v2.7) tool47. From these functional annotations, we selected the most-deleterious consequence

for each variant site and considered the variant allele to be deleterious if its consequence was annotated as transcript ablation, splice acceptor/donor,stop gained, frame shift, stop lost and splice region, initiator codon, nonsynonymous codon predicted to be deleterious/damaging, or in framecodon loss/gain. We also called somatic point mutations and indels with MuTect48/Indelocator (http://archive.broadinstitute.org/cancer/cga/indelocator

Web End =http://archive.broadinstitute.org/cancer/cga/indelocator) and Strelka49. Since MuTect can only detect somatic point mutations it is used together with short indel variant calling method Indelocator. The co-occurrence and mutual exclusivity of genes were assessed using one-sided Fishers exact test. The MutSig algorithm was used for determining the genes that were mutated more often than expected by chance with FDR o5% (ref. 50). Mutational signatures of 6 main categories, T4C/A4G, A4C/T4G, G4C/C4G, C4T/G4A, A4T/T4A, G4T/C4A) are calculated based on ltered somatic variants.

Clonality analysis. Clonality rate of each somatic mutation was calculated based on the variant allele frequency and ploidy at that site, taking into account the admixture rate of each tumour.

CNV identication from exome data. The log ratio of depth of coverage between tumour and blood was calculated using GATK-Depth Of Coverage tool. CNV segments were then called from the log ratio of depth of coverage using ExomeCNV R package51. False positive CNV events were corrected by calculating minor allele frequencies (BAF) in each CNV segment. In each CNV segment, B-allele frequencies (BAF) at heterozygous sites should deviate from 0.5 by at least0.05 units. We estimated the admixture rate based on CNV analysis of paired tumour and blood samples. Copy number loss regions were extracted and for those regions the BAF of each tumour SNP that was heterozygous in blood was calculated. Finally the admixture rate was estimated from the degree of deviation from homozygosity using the qpure R package52.

Custom molecular inversion probe sequencing and analysis. Targeted sequencing of exomic regions and exon-intron boundaries of NF2, SMARCB1, TRAF7, PIK3CA, PIK3R1, PRKAR1A, SMO and SUFU plus the recurrent variants AKT1 p.Glu17Lys and KLF4 p.Lys409Gln was performed using molecular inversion probes (MIPs). Recurrent mutations in POLR2A were assessed with Sanger screening.

Custom amplicon sequencing and analysis. Libraries consisting of the coding exons from TRAF7, NF2, SMO, and the recurrent mutations for AKT1 p.Glu17Lys and KLF4 p.Lys409Gln were created using the TargetRich custom amplicon kit (Kailos Genetics).

Sanger sequencing. Coding variants detected by whole-exome sequencing or targeted next-generation sequencing were conrmed by Sanger sequencing using standard protocols.

Whole-genome genotyping. The Illumina Platform was used for WGG and CNV analyses of the samples. Human OmniExpress-12v1.0 BeadChips that contain 733,202 markers were used according to the manufacturers protocol (Illumina, San Diego, CA, USA). CNVs were detected by comparing the normalized signal intensity between tumour and matched blood or tumour and the average of all blood samples. Segmentation was performed on log intensity (R) ratios using DNACopy algorithm53. Large-scale chromosomal deletion or amplication was dened as affecting more than one-third of the chromosomal arm, whereas focal event deletion or amplication was dened as affecting less than one-third and more than one-tenth of the chromosomal arm with accompanying log ratio of signal intensities o 0.1 or 40.1 and B-allele frequencies (BAF) at heterozygous

sites deviating from 0.5 by at least 0.05 units. Large-scale copy neutral LOH was dened similarly, with the exception of log ratio of signal intensities being between0.1 and 0.1. Percentage of genome alteration (PGA) was dened as the percentage of loss or gain base pairs relative to the entire genome. To identify signicantly altered genes we applied Wilcoxon test between atypical and benign tumours using intensity log ratios. Cancer genes were identied based on the COSMIC database54.

Gene expression data. We used Illumina HumanHT12.v4 chips on the gene expression data. Data was normalized using normal-exponential convolution model-based background correction and quantile normalization using the limma R package55. All batches were normalized at once after excluding probes with low quality. Samples estimated to have zero proportion of expressed probes, mean signal intensity being less than 5.54, or RIN value o3 were excluded. Hierarchical clustering of the gene expression data showed batch effect in the data. The batch effect was removed using ComBat in sva R package56. We performed unsupervised hierarchical clustering based on a Euclidean distance metric and average linking clustering on the probes that showed the top 1,000 most variable gene expression levels. GO term enrichment analysis was performed using the Cytoscape Reactome Plugin. Differentially expressed genes were identied using an empirical Bayesian method ebayes implemented in limma R package55. Genes were considered differentially expressed with adjusted P-valueo0.05.

NATURE COMMUNICATIONS | 8:14433 | DOI: 10.1038/ncomms14433 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 11

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14433

Random forest prediction model using gene expression data. We used the randomForest R package (https://cran.r-project.org/web/packages/randomForest/index.html

Web End =https://cran.r-project.org/web/packages/randomForest/ https://cran.r-project.org/web/packages/randomForest/index.html

Web End =index.html ) for building a model to predict the histological grade57. We rst performed feature selection by identifying the top 25 genes associated with histological grade from differential gene expression analysis. We trained the model on our top 25 differentially expressed genes using randomForest function with ntree 5,000 parameter. Out of bag error rate (OOB), which is the error rate of the

trained algorithm on a left out dataset that is not used in trained algorithm, was calculated. We next validated our prediction model on an independent dataset13. Raw Affymetrix data of the validation set was processed using the RMA method in affy R package (https://www.bioconductor.org/packages/devel/bioc/manuals/affy/man/affy.pdf

Web End =https://www.bioconductor.org/packages/devel/bioc/manuals/affy/ https://www.bioconductor.org/packages/devel/bioc/manuals/affy/man/affy.pdf

Web End =man/affy.pdf ).

miRNA sequencing analyses. Sequencing was performed on Illumina HiSeq instruments using 75 basepair, single-end reads and multiplexing 8 tumour samples per lane. Adaptor sequences were trimmed using cutadapt58. Reads shorter than 18 basepair and reads with quality less than 20 were ltered. Per base quality scores were assessed using FASTX tool. Adaptor trimmed reads were aligned to miRBase using mirDeep2 tool59.

The set of known human miRNA precursors were downloaded from miRBase version 20 (ref. 60). We reported total read counts for 5p and 3p strands. Read counts for each sample were normalized to reads per million reads (RPM) and log2-transformed. Batch effect was corrected using ComBat in sva R package56.

We then performed unsupervised hierarchical clustering based on a Euclidean distance metric and complete linking. Differentially expressed miRNAs were identied using limma R package. An FDR threshold of 0.05 was used.

Six algorithms were used for miRNA target prediction: Miranda, Mirbase, Mirtarget2, Pictar, Tarbase and TargetScan using the RmiR.Hs.miRNA R package61. Regulatory targets of individual miRNAs were dened as those genes having signicant negative correlation with the miRNA (spearman correlation test Po0.05) and prediction support in at least three databases. GO term enrichment analysis was performed using the GOStats R package62.

DNA methylation data. We performed DNA methylation proling on 57 tumour samples and 3 control samples using the Illumina Innium HumanMethylation450 BeadChip, which assesses the level of methylation at over 450,000 CpG sites across the entire genome (covering 99% of RefSeq genes and 96% of CpG islands).

We rst processed raw intensity les (*.idat) and obtained ratio between Illumina methylated probe intensity and total probe intensities (beta-values) using the champ package63. Then, we assessed the quality of the methylation samples and probes. Finally batch effect is corrected using combat algorithm in champ R package63.

Quality control and preprocessing for DNA methylation data. We removed sites containing missing values. The probes targeting a CpG with a SNP were also removed from analysis. Probes targeting the X and Y chromosomes were excluded. The sites having at least 50% samples with detection P-value40.05 are removed.

The detection P-value is calculated using the reported background signal levels of both the methylated and unmethylated channels. After preprocessing the raw data, we performed beta mixture quantile normalization (BMIQ), a normalization correction for the technical differences between the Type I and Type II array designs23. The DNA methylation score for each probe is described as the ratio of intensities between methylated and unmethylated alleles: M/(M U).

Unsupervised clustering for DNA methylation. We performed consensus clustering on methylation b-values with 80% subsampling over 1,000 iterations of hierarchical clustering based on a Pearson correlation distance metric and average linking64.

Identifying differentially methylated sites. Differentially methylated sites were calculated using the empirical bayes method called limma55. Sites were considered statistically signicantly different between the groups if they had an adjusted P-valueo0.05 and median b-value difference 40.1 or o( 0.1). We used the

GREAT tool, which internally maps genomic regions to genes and statistically controls for the fact that genes differ in size and their relative distance to each other23. Raw reads of EZH2 H1hesc ChIP-seq were downloaded from UCSC Genome Browser ENCODE database65. EZH2, ChIP-seq peaks were downloaded from Cistrome website66.

Immunouorescence staining. Meningioma frozen tissue sections were washed with phosphate buffer saline (PBS) for 5 min, then placed in 4% formaldehyde in PBS for 3 min (xation), rinsed in PBS and permeabilized with 0.3% Triton-X100 in PBS for 40 s. The sections were rinsed with PBS 3 times for 5 min each and were incubated in BSDSGS blocking solution (PBS with 1% bovine serum albmin, 5% donkey serum, 5% goat serum, 0.1% glycine, 0.1% lysine) with 0.1% Tween 20. Afterwards, the sections were stained with 1:100 rabbit polyclonal anti-EZH2 antibody (5246P, Cell Signaling, Danvers, MA, USA). After washing, the slides were

incubated 30 min with 1:200 donkey-anti rabbit Alexa Fluor 555 (A315772, Life Technologies, Grand Island, NY, USA). After washing, slides were mounted with Vectashield DAPI medium (H1200, Vector Labs, Burlingame, CA, USA). Images were analyzed under an inverted microscope (Axio Vert A1, Zeiss, Obekochen, Germany) with uorescent light source (X-Cite 120Q, Lumen Dynamics).

H3K27ac and H3K27me3 ChIP-seq. Briey, 10 to 15 frozen sections of each tumour block or dura sample were collected for ChIP-Seq experiments. Tissue was crosslinked with 1% formaldehyde, quenched with glycine, and washed with PBS. Nuclei were extracted by dounce homogenization and resuspended in nuclear lysis buffer containing 0.3% SDS. Chromatin was sheared by sonication with a Q800R2 Sonicator by QSonica (60 min total, amplitude 30, 10 s pulses, 10 s rest). Soluble chromatin was incubated with magnetic beads coated with either H3K27ac antibody (ab4729) or H3K27me3 antibody (ab6002) overnight at 4 C. Chromatin was precipitated using a magnet, washed extensively, and eluted with TE 1%

SDS. Crosslinks were reversed, puried, and subjected to standard Illumina paired-end multiplexed library construction. H3K27ac and H3K27me3 ChIP and input samples were sequenced for each tumour (1 75 bp, HiSeq 2000). H3K27ac and

H3K27me3 reads were aligned uniquely with bowtie (0.12.7)67 to the human genome (hg19). In H3K27ac, regions of enrichment were identied with MACS (v1.4)68 whereas in H3K27me3, regions of enrichment were identied with MUSIC69. We used DiffBind to estimate signicance of super-enhancer read change (adjusted P-value threshold 0.05) between meningioma subtypes70.

Identifying super-enhancers from ChIP-seq data. To identify super-enhancers, rst regions enriched in ChIP-seq reads for H3K27ac were identied using MACS with input control. Super-enhancers were separated from typical enhancers using ROSE pipeline (https://bitbucket.org/young_computation/rose

Web End =https://bitbucket.org/young_computation/rose) with parameters -s (stitching) 12,500, -t (promoter exclusion zone) 2000. We used DiffBindto estimate signicance of super-enhancer read change (adjusted P-value threshold 0.05) between meningioma subtypes70. Cancer genes are dened

based on COSMIC databases54.

Data availability. All somatic mutations identied through exome sequencingof meningiomas were submitted to the COSMIC database previously7. Gene expression data is deposited in GEO database (accession: GSE84263)7. Our novel datasets including DNA methylation, miRNA sequencing, H3K27ac and H3K27me3 ChIP-seq are deposited in GEO database (accession: GSE91376).

References

1. Wiemels, J., Wrensch, M. & Claus, E. B. Epidemiology and etiology of meningioma. J. Neurooncol. 99, 307314 (2010).

2. Louis, D. N. et al. The 2007 WHO classication of tumours of the central nervous system. Acta Neuropathol. 114, 97109 (2007).

3. Louis, D. N. et al. The 2016 World Health Organization Classication of tumors of the central nervous system: a summary. Acta Neuropathol. 131, 803820 (2016).

4. Riemenschneider, M. J., Perry, A. & Reifenberger, G. Histological classication and molecular genetics of meningiomas. Lancet Neurol. 5, 10451054 (2006).

5. Brastianos, P. K. et al. Genomic sequencing of meningiomas identies oncogenic SMO and AKT1 mutations. Nat. Genet. 45, 285289 (2013).

6. Clark, V. E. et al. Genomic analysis of non-NF2 meningiomas reveals mutations in TRAF7, KLF4, AKT1, and SMO. Science 339, 10771080 (2013).

7. Clark, V. E. et al. Recurrent somatic mutations in POLR2A dene a distinct subset of meningiomas. Nat. Genet. 48, 12531259 (2016).

8. Moliterno, J. et al. Survival in patients treated for anaplastic meningioma.J. Neurosurg. 123, 2330 (2015).9. Maier, H., Ofner, D., Hittmair, A., Kitz, K. & Budka, H. Classic, atypical, and anaplastic meningioma: three histopathological subtypes of clinical relevance.J. Neurosurg. 77, 616623 (1992).10. Weber, R. G. et al. Analysis of genomic alterations in benign, atypical, and anaplastic meningiomas: toward a genetic model of meningioma progression. Proc. Natl Acad. Sci. USA 94, 1471914724 (1997).

11. Bai, H. et al. Integrated genomic characterization of IDH1-mutant glioma malignant progression. Nat. Genet. 48, 5966 (2016).

12. Jones, P. A. & Baylin, S. B. The fundamental role of epigenetic events in cancer. Nat. Rev. Genet. 3, 415428 (2002).

13. Lee, Y. et al. Genomic landscape of meningiomas. Brain Pathol. 20, 751762 (2010).

14. Lemee, F. et al. DNA polymerase theta up-regulation is associated with poor survival in breast cancer, perturbs DNA replication, and promotes genetic instability. Proc. Natl Acad. Sci. USA 107, 1339013395 (2010).

15. Rantanen, S., Monni, O., Joensuu, H., Franssila, K. & Knuutila, S. Causes and consequences of BCL2 overexpression in diffuse large B-cell lymphoma. Leuk. Lymphoma 42, 10891098 (2001).

16. Krajewska, M. et al. Immunohistochemical analysis of bcl-2, bax, bcl-X, and mcl-1 expression in prostate cancers. Am. J. Pathol. 148, 15671576 (1996).

12 NATURE COMMUNICATIONS | 8:14433 | DOI: 10.1038/ncomms14433 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14433 ARTICLE

17. Ikegaki, N., Katsumata, M., Minna, J. & Tsujimoto, Y. Expression of bcl-2 in small cell lung carcinoma cells. Cancer Res. 54, 68 (1994).

18. Xie, B., Ding, Q., Han, H. & Wu, D. miRCancer: a microRNA-cancer association database constructed by text mining on literature. Bioinformatics 29, 638644 (2013).

19. Boyerinas, B., Park, S. M., Hau, A., Murmann, A. E. & Peter, M. E. The role of let-7 in cell differentiation and cancer. Endocr. Relat. Cancer 17, F19F36 (2010).

20. Kong, D. et al. Loss of let-7 up-regulates EZH2 in prostate cancer consistent with the acquisition of cancer stem cell signatures that are attenuated by BR-DIM. PLoS ONE 7, e33729 (2012).

21. Erkan, E. P., Breakeeld, X. O. & Saydam, O. miRNA signature of schwannomas: possible role(s) of tumor suppressor miRNAs in benign tumors. Oncotarget 2, 265270 (2011).

22. Powers, J. T. et al. Multiple mechanisms disrupt the let-7 microRNA family in neuroblastoma. Nature 535, 246251 (2016).

23. McLean, C. Y. et al. GREAT improves functional interpretation of cis-regulatory regions. Nat. Biotechnol. 28, 495501 (2010).

24. Ben-Porath, I. et al. An embryonic stem cell-like gene expression signature in poorly differentiated aggressive human tumors. Nat. Genet. 40, 499507 (2008).

25. Mitchell, A. et al. The InterPro protein families database: the classication resource after 15 years. Nucleic Acids Res. 43, D213D221 (2015).

26. Vire, E. et al. The polycomb group protein EZH2 directly controls DNA methylation. Nature 439, 871874 (2006).

27. Ding, X. et al. The polycomb protein Ezh2 impacts on induced pluripotent stem cell generation. Stem Cells Dev. 23, 931940 (2014).

28. Hnisz, D. et al. Super-enhancers in the control of cell identity and disease. Cell 155, 934947 (2013).

29. Aruga, J. The role of Zic genes in neural development. Mol. Cell Neurosci. 26, 205221 (2004).

30. Aruga, J., Nozaki, Y., Hatayama, M., Odaka, Y. S. & Yokota, N. Expression of ZIC family genes in meningiomas and other brain tumors. BMC Cancer 10, 79 (2010).

31. Airaksinen, M. S. & Saarma, M. The GDNF family: signalling, biological functions and therapeutic value. Nat. Rev. Neurosci. 3, 383394 (2002).

32. Goutagny, S. et al. High incidence of activating TERT promoter mutations in meningiomas undergoing malignant progression. Brain Pathol. 24, 184189 (2014).

33. Lee, T. I. et al. Control of developmental regulators by polycomb in human embryonic stem cells. Cell 125, 301313 (2006).

34. Widschwendter, M. et al. Epigenetic stem cell signature in cancer. Nat. Genet. 39, 157158 (2006).

35. Garcia-Fernandez, J. The genesis and evolution of homeobox gene clusters. Nat. Rev. Genet. 6, 881892 (2005).

36. Vir, E. et al. The polycomb group protein EZH2 directly controls DNA methylation. Nature 439, 871874 (2005).

37. Sato, Y. et al. Integrated molecular analysis of clear-cell renal cell carcinoma. Nat. Genet. 45, 860867 (2013).

38. Suv, M.-L. et al. EZH2 is essential for glioblastoma cancer stem cell maintenance. Cancer Res. 69, 92119218 (2009).

39. Lucon, D. R. et al. Downregulation of 14q32 microRNAs in primary human desmoplastic medulloblastoma. Front. Oncol. 3, 254 (2013).

40. Biegel, J. A. et al. Germ-line and acquired mutations of INI1 in atypical teratoid and rhabdoid tumors. Cancer Res. 59, 7479 (1999).

41. van den Munckhof, P., Christiaans, I., Kenter, S. B., Baas, F. & Hulsebos, T. J. Germline SMARCB1 mutation predisposes to multiple meningiomas and schwannomas with preferential location of cranial meningiomas at the falx cerebri. Neurogenetics 13, 17 (2012).

42. Li, H. & Durbin, R. Fast and accurate short read alignment with BurrowsWheeler transform. Bioinformatics 25, 17541760 (2009).

43. Lunter, G. & Goodson, M. Stampy: a statistical algorithm for sensitive and fast mapping of Illumina sequence reads. Genome Res. 21, 936939 (2011).44. DePristo, M. A. et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat. Genet. 43, 491498 (2011).

45. Li, H. A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics 27, 29872993 (2011).

46. Erson-Omay, E. Z. et al. Somatic POLE mutations cause an ultramutated giant cell high-grade glioma subtype with better prognosis. Neurooncology 17, 13561364 (2015).

47. McLaren, W. et al. Deriving the consequences of genomic variants with the Ensembl API and SNP Effect Predictor. Bioinformatics 26, 20692070 (2010).

48. Cibulskis, K. et al. Sensitive detection of somatic point mutations in impure and heterogeneous cancer samples. Nat. Biotechnol. 31, 213219 (2013).

49. Saunders, C. T. et al. Strelka: accurate somatic small-variant calling from sequenced tumor-normal sample pairs. Bioinformatics 28, 18111817 (2012).

50. Lawrence, M. S. et al. Mutational heterogeneity in cancer and the search for new cancer-associated genes. Nature 499, 214218 (2013).

51. Sathirapongsasuti, J. F. et al. Exome sequencing-based copy-number variation and loss of heterozygosity detection: ExomeCNV. Bioinformatics 27, 26482654 (2011).

52. Song, S. et al. qpure: A tool to estimate tumor cellularity from genome-wide single-nucleotide polymorphism proles. PLoS ONE 7, e45835 (2012).

53. Olshen, A. B., Venkatraman, E. S., Lucito, R. & Wigler, M. Circular binary segmentation for the analysis of array-based DNA copy number data. Biostatistics 5, 557572 (2004).

54. Forbes, S. A. et al. COSMIC: mining complete cancer genomes in the catalogue of somatic mutations in cancer. Nucleic Acids Res. 39, D945D950 (2011).55. Smyth, G. K. Limma: linear models for microarray data. 397420 (Springer, 2005).56. Leek, J. T., Johnson, W. E., Parker, H. S., Jaffe, A. E. & Storey, J. D. The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics 28, 882883 (2012).

57. Wiener, A. L. A. M. Classication and regression by randomForest. R News 2, 1822 (2002).

58. Martin, M. Cutadapt removes adapter sequences from high throughput sequencing reads. EMBnet J. 17, 10 (2011).

59. Friedlander, M. R. et al. Discovering microRNAs from deep sequencing data using miRDeep. Nat. Biotechnol. 26, 407415 (2008).

60. Grifths-Jones, S. miRBase: microRNA sequences, targets and gene nomenclature. Nucleic Acids Res. 34, D140D144 (2006).

61. Favero, F. RmiR: package to work with miRNAs and miRNA targets with R, R package version 1.10.0. Available at https://www.bioconductor.org/packages/release/bioc/html/RmiR.html

Web End =https://www.bioconductor.org/packages/ https://www.bioconductor.org/packages/release/bioc/html/RmiR.html

Web End =release/bioc/html/RmiR.html (2006).

62. Falcon, S. & Gentleman, R. Using GOstats to test gene lists for GO term association. Bioinformatics 23, 257258 (2007).

63. Morris, T. J. et al. ChAMP: 450k chip analysis methylation pipeline. Bioinformatics. 30, 428430 (2014).

64. Monti, S., Tamayo, P., Mesirov, J. & Golub, T. Consensus clustering: a resampling-based method for class discovery and visualization of gene expression microarray data. Mach. Learn. 52, 91118 (2003).

65. Consortium TEP. The ENCODE (ENCyclopedia Of DNA Elements) Project. Science 306, 636640 (2004).

66. Liu, T. et al. Cistrome: an integrative platform for transcriptional regulation studies. Genome Biol. 12, R83 (2011).

67. Langmead, B., Trapnell, C., Pop, M. & Salzberg, S. L. Ultrafast and memory-efcient alignment of short DNA sequences to the human genome. Genome Biol. 10, R25 (2009).

68. Zhang, Y. et al. Model-based analysis of ChIP-Seq (MACS). Genome Biol. 9, R137 (2008).

69. Harmanci, A., Rozowsky, J. & Gerstein, M. MUSIC: identication of enriched regions in ChIP-Seq experiments using a mappability-corrected multiscale signal processing framework. Genome Biol. 15, 474 (2014).

70. Ross-Innes, C. S. et al. Differential oestrogen receptor binding is associated with clinical outcome in breast cancer. Nature 481, 389393 (2012).

Acknowledgements

This study was supported by Gregory M. Kiez and Mehmet Kutman Foundation and Yale School of Medicine funds. Partial funding was provided through a research agreement between Gilead Sciences, Inc and Yale University. V.E.C. and M.W.Y. are supported by NIH Medical Scientist Training Program Grant T32GM007205. We gratefully acknowledge Arif O. Harmanci for constructive comments and discussions on the study.

Author contributions

A.S.H. performed exome CNA, SNP array, DNA methylation array, ChIP-seq and miRNA expression analysis. A.S.H., E.Z.E.-O. and K.Y. performed mRNA expression analysis. A.S.H., E.Z.E.-O., V.E.C. and K.Y. performed whole-exome sequencing analysis.V.E.C., D.D.M., L.K. and M.W.Y. performed targeted next-generation sequencing, Sanger conrmation and mutation prole/clinical characterization. E.Z.E.-O. performed clonality analysis. V.E.C., O.H., G.C.-G., M.W.Y. aided in the selection of and prepared samples. A.S.H. and B.J.A. performed super-enhancer analysis. M.S., J.S. and B.K. provided samples and clinical data. M.W.Y., S.C. and V.E.C. performed ChIP-Seq experiments. K.B. supervised genomic experiments. O.H. and B.B. conducted immunouorescence staining experiment. A.S.H. and M.G. wrote the manuscript. M.G. designed and oversaw the project.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing nancial interests: Partial funding for sequencing of the tumour samples and partial salary support for S.C. and K.M.-G. was provided through a research agreement between Gilead Sciences, Inc. and Yale University. All other authors declare no competing nancial interests.

NATURE COMMUNICATIONS | 8:14433 | DOI: 10.1038/ncomms14433 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 13

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14433

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Harmanc, A. S. et al. Integrated genomic analyses of de novo pathways underlying atypical meningiomas. Nat. Commun. 8, 14433 doi: 10.1038/ncomms14433 (2017).

Publishers note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional afliations.

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/

r The Author(s) 2017

14 NATURE COMMUNICATIONS | 8:14433 | DOI: 10.1038/ncomms14433 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

Word count: 12208

Show less

Abstract

Translate

Details

Title

Integrated genomic analyses of de novo pathways underlying atypical meningiomas

Author

Harmanci, Akdes Serin; Youngblood, Mark W; Clark, Victoria E; Coskun, Süleyman; Henegariu, Octavian; Duran, Daniel; Erson-omay, E Zeynep; Kaulen, Leon D; Lee, Tong Ihn; Abraham, Brian J; Simon, Matthias; Krischek, Boris; Timmer, Marco; Goldbrunner, Roland; Omay, S Bülent; Baranoski, Jacob; Baran, Burçin; Carrión-grant, Geneive; Bai, Hanwen; Mishra-gorur, Ketu; Schramm, Johannes; Moliterno, Jennifer; Vortmeyer, Alexander O; Bilgüvar, Kaya; Yasuno, Katsuhito; Young, Richard A; Günel, Murat

Pages

14433

Publication year

2017

Publication date

Feb 2017

Publisher

Nature Publishing Group

e-ISSN

20411723

Source type

Scholarly Journal

Language of publication

English

DOI

https://doi.org/10.1038/ncomms14433

ProQuest document ID

1867965900

Integrated genomic analyses of de novo pathways underlying atypical meningiomas

Jump to:

Full text

Abstract

Details

Suggested sources