Numerous large-scale genomic studies of matched tumor-normal samples have established the somatic landscapes of most cancer types. However, the downstream analysis of data from somatic mutations entails a number of computational and statistical approaches, requiring usage of independent software and numerous tools. Here, we describe an R Bioconductor package, Maftools, which offers a multitude of analysis and visualization modules that are commonly used in cancer genomic studies, including driver gene identification, pathway, signature, enrichment, and association analyses. Maftools only requires somatic variants in Mutation Annotation Format (MAF) and is independent of larger alignment files. With the implementation of well-established statistical and computational methods, Maftools facilitates data-driven research and comparative analysis to discover novel results from publicly available data sets. In the present study, using three of the well-annotated cohorts from The Cancer Genome Atlas (TCGA), we describe the application of Maftools to reproduce known results. More importantly, we show that Maftools can also be used to uncover novel findings through integrative analysis.

Glioblastoma multiforme (GBM) is a lethal brain tumour in adults and children. However, DNA copy number and gene expression signatures indicate differences between adult and paediatric cases. To explore the genetic events underlying this distinction, we sequenced the exomes of 48 paediatric GBM samples. Somatic mutations in the H3.3-ATRX-DAXX chromatin remodelling pathway were identified in 44% of tumours (21/48). Recurrent mutations in H3F3A, which encodes the replication-independent histone 3 variant H3.3, were observed in 31% of tumours, and led to amino acid substitutions at two critical positions within the histone tail (K27M, G34R/G34V) involved in key regulatory post-translational modifications. Mutations in ATRX (α-thalassaemia/mental retardation syndrome X-linked) and DAXX (death-domain associated protein), encoding two subunits of a chromatin remodelling complex required for H3.3 incorporation at pericentric heterochromatin and telomeres, were identified in 31% of samples overall, and in 100% of tumours harbouring a G34R or G34V H3.3 mutation. Somatic TP53 mutations were identified in 54% of all cases, and in 86% of samples with H3F3A and/or ATRX mutations. Screening of a large cohort of gliomas of various grades and histologies (n = 784) showed H3F3A mutations to be specific to GBM and highly prevalent in children and young adults. Furthermore, the presence of H3F3A/ATRX-DAXX/TP53 mutations was strongly associated with alternative lengthening of telomeres and specific gene expression profiles. This is, to our knowledge, the first report to highlight recurrent mutations in a regulatory histone in humans, and our data suggest that defects of the chromatin architecture underlie paediatric and young adult GBM pathogenesis.

Monozygous twins share a common genotype. However, most monozygotic twin pairs are not identical; several types of phenotypic discordance may be observed, such as differences in susceptibilities to disease and a wide range of anthropomorphic features. There are several possible explanations for these observations, but one is the existence of epigenetic differences. To address this issue, we examined the global and locus-specific differences in DNA methylation and histone acetylation of a large cohort of monozygotic twins. We found that, although twins are epigenetically indistinguishable during the early years of life, older monozygous twins exhibited remarkable differences in their overall content and genomic distribution of 5-methylcytosine DNA and histone acetylation, affecting their gene-expression portrait. These findings indicate how an appreciation of epigenetics is missing from our understanding of how different phenotypes can be originated from the same genotype.DNA methylation ͉ epigenetics ͉ histones H uman monozygotic (MZ) twins account for 1 in 250 live births (1). The origin of MZ twins is attributed to two or more daughter cells of a single zygote undergoing independent mitotic divisions, leading to independent development and births. They are considered genetically identical, but significant phenotypic discordance between them may exist. This quality is particularly noticeable for psychiatric diseases, such as schizophrenia and bipolar disorder (2). MZ twins have been used to demonstrate the role of environmental factors in determining complex diseases and phenotypes, but the true nature of the phenotypic discordance nevertheless remains extremely poorly understood. In this context, differences in the placenta, amniotic sac, and vascularization of the separate cell masses or even mosaicism in genetic and cytogenetic markers in MZ may exist (3), although the published studies are very few in number. Thus, the real causes for MZ twin discordance for common diseases and traits remain to be established. Epigenetic differences may be an important part of the solution to this puzzle. Indeed, epigenetic profiles may represent the link between an environmental factor and phenotypic differences in MZ twins. Cloned animals provide another example of how epigenetics may explain phenotypic differences in beings that have identical genetic sequences. In this case, inefficient epigenetic reprogramming of the transplanted nucleus is associated with aberrations in imprinting, aberrant growth, and lethality beyond a threshold of faulty epigenetic control (4). MZ twins are another phenomenon in which epigenetics can ''make the difference.'' To address this possibility, we have profiled the epigenetic patterns related to global and locus-specific DNA methylation and histone H3 and H4 acetylation in the largest series of MZ twins for which molecular studies have been reported. Materials and MethodsSubjects. Eighty volunteer Caucasian twins from Spain were recruited in the study, including 30 male and 50 female subje...

Glioblastoma (GBM) is a brain tumor that carries a dismal prognosis and displays considerable heterogeneity. We have recently identified recurrent H3F3A mutations affecting two critical amino acids (K27 and G34) of histone H3.3 in one-third of pediatric GBM. Here, we show that each H3F3A mutation defines an epigenetic subgroup of GBM with a distinct global methylation pattern, and that they are mutually exclusive with IDH1 mutations, which characterize a third mutation-defined subgroup. Three further epigenetic subgroups were enriched for hallmark genetic events of adult GBM and/or established transcriptomic signatures. We also demonstrate that the two H3F3A mutations give rise to GBMs in separate anatomic compartments, with differential regulation of transcription factors OLIG1, OLIG2, and FOXG1, possibly reflecting different cellular origins.

The pan-cancer analysis of whole genomes The expansion of whole-genome sequencing studies from individual ICGC and TCGA working groups presented the opportunity to undertake a meta-analysis of genomic features across tumour types. To achieve this, the PCAWG Consortium was established. A Technical Working Group implemented the informatics analyses by aggregating the raw sequencing data from different working groups that studied individual tumour types, aligning the sequences to the human genome and delivering a set of high-quality somatic mutation calls for downstream analysis (Extended Data Fig. 1). Given the recent meta-analysis

CpG islands frequently contain gene promoters or exons and are usually unmethylated in normal cells. Methylation of CpG islands is associated with delayed replication, condensed chromatin and inhibition of transcription initiation. The investigation of aberrant CpG-island methylation in human cancer has primarily taken a candidate gene approach, and has focused on less than 15 of the estimated 45,000 CpG islands in the genome. Here we report a global analysis of the methylation status of 1,184 unselected CpG islands in each of 98 primary human tumours using restriction landmark genomic scanning (RLGS). We estimate that an average of 600 CpG islands (range of 0 to 4,500) of the 45,000 in the genome were aberrantly methylated in the tumours, including early stage tumours. We identified patterns of CpG-island methylation that were shared within each tumour type, together with patterns and targets that displayed distinct tumour-type specificity. The expression of many of these genes was reactivated by experimental demethylation in cultured tumour cells. Thus, the methylation of particular subsets of CpG islands may have consequences for specific tumour types.

Two recurrent mutations, K27M and G34R/V, within histone variant H3.3 were recently identified in ∼50% of pHGGs. Both mutations define clinically and biologically distinct subgroups of pHGGs. Here, we provide further insight about the dominant-negative effect of K27M mutant H3.3, leading to a global reduction of the repressive histone mark H3K27me3. We demonstrate that this is caused by aberrant recruitment of the PRC2 complex to K27M mutant H3.3 and enzymatic inhibition of the H3K27me3-establishing methyltransferase EZH2. By performing chromatin immunoprecipitation followed by next-generation sequencing and whole-genome bisulfite sequencing in primary pHGGs, we show that reduced H3K27me3 levels and DNA hypomethylation act in concert to activate gene expression in K27M mutant pHGGs.

In this study, we present integrated quantitative proteome, transcriptome, and methylome analyses of hematopoietic stem cells (HSCs) and four multipotent progenitor (MPP) populations. From the characterization of more than 6,000 proteins, 27,000 transcripts, and 15,000 differentially methylated regions (DMRs), we identified coordinated changes associated with early differentiation steps. DMRs show continuous gain or loss of methylation during differentiation, and the overall change in DNA methylation correlates inversely with gene expression at key loci. Our data reveal the differential expression landscape of 493 transcription factors and 682 lncRNAs and highlight specific expression clusters operating in HSCs. We also found an unexpectedly dynamic pattern of transcript isoform regulation, suggesting a critical regulatory role during HSC differentiation, and a cell cycle/DNA repair signature associated with multipotency in MPP2 cells. This study provides a comprehensive genome-wide resource for the functional exploration of molecular, cellular, and epigenetic regulation at the top of the hematopoietic hierarchy.