Wilms tumour is an embryonal tumour of childhood that closely resembles the developing kidney. Genomic changes responsible for the development of the majority of Wilms tumours remain largely unknown. Here we identify recurrent mutations within Wilms tumours that involve the highly conserved YEATS domain of MLLT1 (ENL), a gene known to be involved in transcriptional elongation during early development. The mutant MLLT1 protein shows altered binding to acetylated histone tails. Moreover, MLLT1-mutant tumours show an increase in MYC gene expression and HOX dysregulation. Patients with MLLT1-mutant tumours present at a younger age and have a high prevalence of precursor intralobar nephrogenic rests. These data support a model whereby activating MLLT1 mutations early in renal development result in the development of Wilms tumour.

Glioblastoma (GBM) is the most lethal and aggressive adult brain tumor, requiring the development of efficacious therapeutics. Towards this goal, we screened five genetically distinct patient-derived brain-tumor initiating cell lines (BTIC) with a unique collection of small molecule epigenetic modulators from the Structural Genomics Consortium (SGC). We identified multiple hits that inhibited the growth of BTICs in vitro, and further evaluated the therapeutic potential of EZH2 and HDAC inhibitors due to the high relevance of these targets for GBM. We found that the novel SAM-competitive EZH2 inhibitor UNC1999 exhibited low micromolar cytotoxicity in vitro on a diverse collection of BTIC lines, synergized with dexamethasone (DEX) and suppressed tumor growth in vivo in combination with DEX. In addition, a unique brain-penetrant class I HDAC inhibitor exhibited cytotoxicity in vitro on a panel of BTIC lines and extended survival in combination with TMZ in an orthotopic BTIC model in vivo. Finally, a combination of EZH2 and HDAC inhibitors demonstrated synergy in vitro by augmenting apoptosis and increasing DNA damage. Our findings identify key epigenetic modulators in GBM that regulate BTIC growth and survival and highlight promising combination therapies.

Posterior fossa ependymoma comprises two distinct molecular variants termed EPN_PFA and EPN_PFB that have a distinct biology and natural history. The therapeutic value of cytoreductive surgery and radiation therapy for posterior fossa ependymoma after accounting for molecular subgroup is not known.

DDX3X encodes a DEAD-box family RNA helicase (DDX3) commonly mutated in medulloblastoma, a highly aggressive cerebellar tumor affecting both children and adults. Despite being implicated in several facets of RNA metabolism, the nature and scope of DDX3's interactions with RNA remain unclear. Here, we show DDX3 collaborates extensively with the translation initiation machinery through direct binding to 5'UTRs of nearly all coding RNAs, specific sites on the 18S rRNA, and multiple components of the translation initiation complex. Impairment of translation initiation is also evident in primary medulloblastomas harboring mutations in DDX3X, further highlighting DDX3's role in this process. Arsenite-induced stress shifts DDX3 binding from the 5'UTR into the coding region of mRNAs concomitant with a general reduction of translation, and both the shift of DDX3 on mRNA and decreased translation are blunted by expression of a catalytically-impaired, medulloblastoma-associated DDX3R534H variant. Furthermore, despite the global repression of translation induced by arsenite, translation is preserved on select genes involved in chromatin organization in DDX3R534H-expressing cells. Thus, DDX3 interacts extensively with RNA and ribosomal machinery to help remodel the translation landscape in response to stress, while cancer-related DDX3 variants adapt this response to selectively preserve translation.

Developmental history shapes the epigenome and biological function of differentiated cells. Epigenomic patterns have been broadly attributed to the three embryonic germ layers. Here we investigate how developmental origin influences epigenomes. We compare key epigenomes of cell types derived from surface ectoderm (SE), including keratinocytes and breast luminal and myoepithelial cells, against neural crest-derived melanocytes and mesoderm-derived dermal fibroblasts, to identify SE differentially methylated regions (SE-DMRs). DNA methylomes of neonatal keratinocytes share many more DMRs with adult breast luminal and myoepithelial cells than with melanocytes and fibroblasts from the same neonatal skin. This suggests that SE origin contributes to DNA methylation patterning, while shared skin tissue environment has limited effect on epidermal keratinocytes. Hypomethylated SE-DMRs are in proximity to genes with SE relevant functions. They are also enriched for enhancer- and promoter-associated histone modifications in SE-derived cells, and for binding motifs of transcription factors important in keratinocyte and mammary gland biology. Thus, epigenomic analysis of cell types with common developmental origin reveals an epigenetic signature that underlies a shared gene regulatory network.

Synovial sarcoma (SS) is a rare sarcoma with distinct morphologic and genetic features, which, despite its name, does not arise from synovium. While most SSs (>80%) arise in the deep soft tissue of the extremities, up to 5% of these tumors are encountered in the body axis including the spine, mediastinum, retroperitoneum, and head/neck regions. Reports of SS located within the spinal axis have been rare to date.

Next generation sequencing has now enabled a cost-effective enumeration of the full mutational complement of a tumor genome-in particular single nucleotide variants (SNVs). Most current computational and statistical models for analyzing next generation sequencing data, however, do not account for cancer-specific biological properties, including somatic segmental copy number alterations (CNAs)-which require special treatment of the data. Here we present CoNAn-SNV (Copy Number Annotated SNV): a novel algorithm for the inference of single nucleotide variants (SNVs) that overlap copy number alterations. The method is based on modelling the notion that genomic regions of segmental duplication and amplification induce an extended genotype space where a subset of genotypes will exhibit heavily skewed allelic distributions in SNVs (and therefore render them undetectable by methods that assume diploidy). We introduce the concept of modelling allelic counts from sequencing data using a panel of Binomial mixture models where the number of mixtures for a given locus in the genome is informed by a discrete copy number state given as input. We applied CoNAn-SNV to a previously published whole genome shotgun data set obtained from a lobular breast cancer and show that it is able to discover 21 experimentally revalidated somatic non-synonymous mutations in a lobular breast cancer genome that were not detected using copy number insensitive SNV detection algorithms. Importantly, ROC analysis shows that the increased sensitivity of CoNAn-SNV does not result in disproportionate loss of specificity. This was also supported by analysis of a recently published lymphoma genome with a relatively quiescent karyotype, where CoNAn-SNV showed similar results to other callers except in regions of copy number gain where increased sensitivity was conferred. Our results indicate that in genomically unstable tumors, copy number annotation for SNV detection will be critical to fully characterize the mutational landscape of cancer genomes.

Effective treatment of diffuse large B-cell lymphoma (DLBCL) is plagued by heterogeneous responses to standard therapy, and molecular mechanisms underlying unfavorable outcomes in lymphoma patients remain elusive. Here, we profiled 148 genomes with 91 matching transcriptomes in a DLBCL cohort treated with rituximab plus cyclophosphamide, doxorubicin, vincristine, and prednisolone (R-CHOP) to uncover molecular subgroups linked to treatment failure. Systematic integration of high-resolution genotyping arrays and RNA sequencing data revealed novel deletions in RCOR1 to be associated with unfavorable progression-free survival (P = .001). Integration of expression data from the clinical samples with data from RCOR1 knockdowns in the lymphoma cell lines KM-H2 and Raji yielded an RCOR1 loss-associated gene signature comprising 233 genes. This signature identified a subgroup of patients with unfavorable overall survival (P = .023). The prognostic significance of the 233-gene signature for overall survival was reproduced in an independent cohort comprising 195 R-CHOP-treated patients (P = .039). Additionally, we discovered that within the International Prognostic Index low-risk group, the gene signature provides additional prognostic value that was independent of the cell-of-origin phenotype. We present a novel and reproducible molecular subgroup of DLBCL that impacts risk-stratification of R-CHOP-treated DLBCL patients and reveals a possible new avenue for therapeutic intervention strategies.

The role of intermediate methylation states in DNA is unclear. Here, to comprehensively identify regions of intermediate methylation and their quantitative relationship with gene activity, we apply integrative and comparative epigenomics to 25 human primary cell and tissue samples. We report 18,452 intermediate methylation regions located near 36% of genes and enriched at enhancers, exons and DNase I hypersensitivity sites. Intermediate methylation regions average 57% methylation, are predominantly allele-independent and are conserved across individuals and between mouse and human, suggesting a conserved function. These regions have an intermediate level of active chromatin marks and their associated genes have intermediate transcriptional activity. Exonic intermediate methylation correlates with exon inclusion at a level between that of fully methylated and unmethylated exons, highlighting gene context-dependent functions. We conclude that intermediate DNA methylation is a conserved signature of gene regulation and exon usage.

We report a case of early-onset pancreatic ductal adenocarcinoma in a patient harboring biallelic MUTYH germline mutations, whose tumor featured somatic mutational signatures consistent with defective MUTYH-mediated base excision repair and the associated driver KRAS transversion mutation p.Gly12Cys. Analysis of an additional 730 advanced cancer cases (N = 731) was undertaken to determine whether the mutational signatures were also present in tumors from germline MUTYH heterozygote carriers or if instead the signatures were only seen in those with biallelic loss of function. We identified two patients with breast cancer each carrying a pathogenic germline MUTYH variant with a somatic MUTYH copy loss leading to the germline variant being homozygous in the tumor and demonstrating the same somatic signatures. Our results suggest that monoallelic inactivation of MUTYH is not sufficient for C:G>A:T transversion signatures previously linked to MUTYH deficiency to arise (N = 9), but that biallelic complete loss of MUTYH function can cause such signatures to arise even in tumors not classically seen in MUTYH-associated polyposis (N = 3). Although defective MUTYH is not the only determinant of these signatures, MUTYH germline variants may be present in a subset of patients with tumors demonstrating elevated somatic signatures possibly suggestive of MUTYH deficiency (e.g., COSMIC Signature 18, SigProfiler SBS18/SBS36, SignatureAnalyzer SBS18/SBS36).

Temozolomide (TMZ) is active against glioblastomas (GBM) in which the O6-methylguanine-DNA methyltransferase (MGMT) gene is silenced. However, even in responsive cases, its beneficial effect is undermined by the emergence of drug resistance. Here, we tested whether inhibition of poly (ADP-ribose) polymerase-1 and -2 (PARP) enhanced the effectiveness of TMZ.

Telomerase reverse transcriptase (TERT) promoter mutations were recently shown to drive telomerase activity in various cancer types, including medulloblastoma. However, the clinical and biological implications of TERT mutations in medulloblastoma have not been described. Hence, we sought to describe these mutations and their impact in a subgroup-specific manner. We analyzed the TERT promoter by direct sequencing and genotyping in 466 medulloblastomas. The mutational distributions were determined according to subgroup affiliation, demographics, and clinical, prognostic, and molecular features. Integrated genomics approaches were used to identify specific somatic copy number alterations in TERT promoter-mutated and wild-type tumors. Overall, TERT promoter mutations were identified in 21 % of medulloblastomas. Strikingly, the highest frequencies of TERT mutations were observed in SHH (83 %; 55/66) and WNT (31 %; 4/13) medulloblastomas derived from adult patients. Group 3 and Group 4 harbored this alteration in <5 % of cases and showed no association with increased patient age. The prognostic implications of these mutations were highly subgroup-specific. TERT mutations identified a subset with good and poor prognosis in SHH and Group 4 tumors, respectively. Monosomy 6 was mostly restricted to WNT tumors without TERT mutations. Hallmark SHH focal copy number aberrations and chromosome 10q deletion were mutually exclusive with TERT mutations within SHH tumors. TERT promoter mutations are the most common recurrent somatic point mutation in medulloblastoma, and are very highly enriched in adult SHH and WNT tumors. TERT mutations define a subset of SHH medulloblastoma with distinct demographics, cytogenetics, and outcomes.

Amplification of the C19MC oncogenic miRNA cluster and high LIN28 expression has been linked to a distinctly aggressive group of cerebral CNS-PNETs (group 1 CNS-PNETs) arising in young children. In this study, we sought to evaluate the diagnostic specificity of C19MC and LIN28, and the clinical and biological spectra of C19MC amplified and/or LIN28+ CNS-PNETs. We interrogated 450 pediatric brain tumors using FISH and IHC analyses and demonstrate that C19MC alteration is restricted to a sub-group of CNS-PNETs with high LIN28 expression; however, LIN28 immunopositivity was not exclusive to CNS-PNETs but was also detected in a proportion of other malignant pediatric brain tumors including rhabdoid brain tumors and malignant gliomas. C19MC amplified/LIN28+ group 1 CNS-PNETs arose predominantly in children <4 years old; a majority arose in the cerebrum but 24 % (13/54) of tumors had extra-cerebral origins. Notably, group 1 CNS-PNETs encompassed several histologic classes including embryonal tumor with abundant neuropil and true rosettes (ETANTR), medulloepithelioma, ependymoblastoma and CNS-PNETs with variable differentiation. Strikingly, gene expression and methylation profiling analyses revealed a common molecular signature enriched for primitive neural features, high LIN28/LIN28B and DNMT3B expression for all group 1 CNS-PNETs regardless of location or tumor histology. Our collective findings suggest that current known histologic categories of CNS-PNETs which include ETANTRs, medulloepitheliomas, ependymoblastomas in various CNS locations, comprise a common molecular and diagnostic entity and identify inhibitors of the LIN28/let7/PI3K/mTOR axis and DNMT3B as promising therapeutics for this distinct histogenetic entity.

While molecular subgrouping has revolutionized medulloblastoma classification, the extent of heterogeneity within subgroups is unknown. Similarity network fusion (SNF) applied to genome-wide DNA methylation and gene expression data across 763 primary samples identifies very homogeneous clusters of patients, supporting the presence of medulloblastoma subtypes. After integration of somatic copy-number alterations, and clinical features specific to each cluster, we identify 12 different subtypes of medulloblastoma. Integrative analysis using SNF further delineates group 3 from group 4 medulloblastoma, which is not as readily apparent through analyses of individual data types. Two clear subtypes of infants with Sonic Hedgehog medulloblastoma with disparate outcomes and biology are identified. Medulloblastoma subtypes identified through integrative clustering have important implications for stratification of future clinical trials.

Loss of heterozygosity (LOH) and copy number alteration (CNA) feature prominently in the somatic genomic landscape of tumors. As such, karyotypic aberrations in cancer genomes have been studied extensively to discover novel oncogenes and tumor-suppressor genes. Advances in sequencing technology have enabled the cost-effective detection of tumor genome and transcriptome mutation events at single-base-pair resolution; however, computational methods for predicting segmental regions of LOH in this context are not yet fully explored. Consequently, whole transcriptome, nucleotide-level resolution analysis of monoallelic expression patterns associated with LOH has not yet been undertaken in cancer. We developed a novel approach for inference of LOH from paired tumor/normal sequence data and applied it to a cohort of 23 triple-negative breast cancer (TNBC) genomes. Following extensive benchmarking experiments, we describe the nucleotide-resolution landscape of LOH in TNBC and assess the consequent effect of LOH on the transcriptomes of these tumors using RNA-seq-derived measurements of allele-specific expression. We show that the majority of monoallelic expression in the transcriptomes of triple-negative breast cancer can be explained by genomic regions of LOH and establish an upper bound for monoallelic expression that may be explained by other tumor-specific modifications such as epigenetics or mutations. Monoallelically expressed genes associated with LOH reveal that cell cycle, homologous recombination and actin-cytoskeletal functions are putatively disrupted by LOH in TNBC. Finally, we show how inference of LOH can be used to interpret allele frequencies of somatic mutations and postulate on temporal ordering of mutations in the evolutionary history of these tumors.

DNA methylation plays key roles in diverse biological processes such as X chromosome inactivation, transposable element repression, genomic imprinting, and tissue-specific gene expression. Sequencing-based DNA methylation profiling provides an unprecedented opportunity to map and compare complete DNA methylomes. This includes one of the most widely applied technologies for measuring DNA methylation: methylated DNA immunoprecipitation followed by sequencing (MeDIP-seq), coupled with a complementary method, methylation-sensitive restriction enzyme sequencing (MRE-seq). A computational approach that integrates data from these two different but complementary assays and predicts methylation differences between samples has been unavailable. Here, we present a novel integrative statistical framework M&M (for integration of MeDIP-seq and MRE-seq) that dynamically scales, normalizes, and combines MeDIP-seq and MRE-seq data to detect differentially methylated regions. Using sample-matched whole-genome bisulfite sequencing (WGBS) as a gold standard, we demonstrate superior accuracy and reproducibility of M&M compared to existing analytical methods for MeDIP-seq data alone. M&M leverages the complementary nature of MeDIP-seq and MRE-seq data to allow rapid comparative analysis between whole methylomes at a fraction of the cost of WGBS. Comprehensive analysis of nineteen human DNA methylomes with M&M reveals distinct DNA methylation patterns among different tissue types, cell types, and individuals, potentially underscoring divergent epigenetic regulation at different scales of phenotypic diversity. We find that differential DNA methylation at enhancer elements, with concurrent changes in histone modifications and transcription factor binding, is common at the cell, tissue, and individual levels, whereas promoter methylation is more prominent in reinforcing fundamental tissue identities.

Follicular lymphoma (FL) and diffuse large B-cell lymphoma (DLBCL) are the two most common non-Hodgkin lymphomas (NHLs). Here we sequenced tumour and matched normal DNA from 13 DLBCL cases and one FL case to identify genes with mutations in B-cell NHL. We analysed RNA-seq data from these and another 113 NHLs to identify genes with candidate mutations, and then re-sequenced tumour and matched normal DNA from these cases to confirm 109 genes with multiple somatic mutations. Genes with roles in histone modification were frequent targets of somatic mutation. For example, 32% of DLBCL and 89% of FL cases had somatic mutations in MLL2, which encodes a histone methyltransferase, and 11.4% and 13.4% of DLBCL and FL cases, respectively, had mutations in MEF2B, a calcium-regulated gene that cooperates with CREBBP and EP300 in acetylating histones. Our analysis suggests a previously unappreciated disruption of chromatin biology in lymphomagenesis.

To determine the maximum tolerated dose (MTD) and safety, and explore efficacy and biomarkers of vandetanib with cetuximab and irinotecan in second-line metastatic colorectal cancer.

A crucial step in the development of muscle cells in all metazoan animals is the assembly and anchorage of the sarcomere, the essential repeat unit responsible for muscle contraction. In Caenorhabditis elegans, many of the critical proteins involved in this process have been uncovered through mutational screens focusing on uncoordinated movement and embryonic arrest phenotypes. We propose that additional sarcomeric proteins exist for which there is a less severe, or entirely different, mutant phenotype produced in their absence. We have used Serial Analysis of Gene Expression (SAGE) to generate a comprehensive profile of late embryonic muscle gene expression. We generated two replicate long SAGE libraries for sorted embryonic muscle cells, identifying 7,974 protein-coding genes. A refined list of 3,577 genes expressed in muscle cells was compiled from the overlap between our SAGE data and available microarray data. Using the genes in our refined list, we have performed two separate RNA interference (RNAi) screens to identify novel genes that play a role in sarcomere assembly and/or maintenance in either embryonic or adult muscle. To identify muscle defects in embryos, we screened specifically for the Pat embryonic arrest phenotype. To visualize muscle defects in adult animals, we fed dsRNA to worms producing a GFP-tagged myosin protein, thus allowing us to analyze their myofilament organization under gene knockdown conditions using fluorescence microscopy. By eliminating or severely reducing the expression of 3,300 genes using RNAi, we identified 122 genes necessary for proper myofilament organization, 108 of which are genes without a previously characterized role in muscle. Many of the genes affecting sarcomere integrity have human homologs for which little or nothing is known.

The cause of mental retardation in one-third to one-half of all affected individuals is unknown. Microscopically detectable chromosomal abnormalities are the most frequently recognized cause, but gain or loss of chromosomal segments that are too small to be seen by conventional cytogenetic analysis has been found to be another important cause. Array-based methods offer a practical means of performing a high-resolution survey of the entire genome for submicroscopic copy-number variants. We studied 100 children with idiopathic mental retardation and normal results of standard chromosomal analysis, by use of whole-genome sampling analysis with Affymetrix GeneChip Human Mapping 100K arrays. We found de novo deletions as small as 178 kb in eight cases, de novo duplications as small as 1.1 Mb in two cases, and unsuspected mosaic trisomy 9 in another case. This technology can detect at least twice as many potentially pathogenic de novo copy-number variants as conventional cytogenetic analysis can in people with mental retardation.