Multiple myeloma is an incurable malignancy of plasma cells, and its pathogenesis is poorly understood. Here we report the massively parallel sequencing of 38 tumour genomes and their comparison to matched normal DNAs. Several new and unexpected oncogenic mechanisms were suggested by the pattern of somatic mutation across the data set. These include the mutation of genes involved in protein translation (seen in nearly half of the patients), genes involved in histone methylation, and genes involved in blood coagulation. In addition, a broader than anticipated role of NF-κB signalling was indicated by mutations in 11 members of the NF-κB pathway. Of potential immediate clinical relevance, activating mutations of the kinase BRAF were observed in 4% of patients, suggesting the evaluation of BRAF inhibitors in multiple myeloma clinical trials. These results indicate that cancer genome sequencing of large collections of samples will yield new insights into cancer not anticipated by existing knowledge.

The transcription factor GATA-1 is required for terminal erythroid maturation and functions as an activator or repressor depending on gene context. Yet its in vivo site selectivity and ability to distinguish between activated versus repressed genes remain incompletely understood. In this study, we performed GATA-1 ChIP-seq in erythroid cells and compared it to GATA-1-induced gene expression changes. Bound and differentially expressed genes contain a greater number of GATA-binding motifs, a higher frequency of palindromic GATA sites, and closer occupancy to the transcriptional start site versus nondifferentially expressed genes. Moreover, we show that the transcription factor Zbtb7a occupies GATA-1-bound regions of some direct GATA-1 target genes, that the presence of SCL/TAL1 helps distinguish transcriptional activation versus repression, and that polycomb repressive complex 2 (PRC2) is involved in epigenetic silencing of a subset of GATA-1-repressed genes. These data provide insights into GATA-1-mediated gene regulation in vivo.

Activation of the RAS-RAF-MEK-ERK signaling pathway is thought to be the key driver of pediatric low-grade astrocytoma (PLGA) growth. Sorafenib is a multikinase inhibitor targeting BRAF, VEGFR, PDGFR, and c-kit. This multicenter phase II study was conducted to determine the response rate to sorafenib in patients with recurrent or progressive PLGA.

The SMARCA4 (also known as BRG1 in humans) chromatin remodeling factor is critical for establishing lineage-specific chromatin states during early mammalian development. However, the role of SMARCA4 in tissue-specific gene regulation during embryogenesis remains poorly defined. To investigate the genome-wide binding landscape of SMARCA4 in differentiating tissues, we engineered a Smarca4(FLAG) knock-in mouse line. Using ChIP-seq, we identified ∼51,000 SMARCA4-associated regions across six embryonic mouse tissues (forebrain, hindbrain, neural tube, heart, limb, and face) at mid-gestation (E11.5). The majority of these regions was distal from promoters and showed dynamic occupancy, with most distal SMARCA4 sites (73%) confined to a single or limited subset of tissues. To further characterize these regions, we profiled active and repressive histone marks in the same tissues and examined the intersection of informative chromatin states and SMARCA4 binding. This revealed distinct classes of distal SMARCA4-associated elements characterized by activating and repressive chromatin signatures that were associated with tissue-specific up- or down-regulation of gene expression and relevant active/repressed biological pathways. We further demonstrate the predicted active regulatory properties of SMARCA4-associated elements by retrospective analysis of tissue-specific enhancers and direct testing of SMARCA4-bound regions in transgenic mouse assays. Our results indicate a dual active/repressive function of SMARCA4 at distal regulatory sequences in vivo and support its role in tissue-specific gene regulation during embryonic development.

Translocation events are frequent in cancer and may create chimeric fusions or 'regulatory rearrangements' that drive oncogene overexpression. Here we identify super-enhancer translocations that drive overexpression of the oncogenic transcription factor MYB as a recurrent theme in adenoid cystic carcinoma (ACC). Whole-genome sequencing data and chromatin maps highlight distinct chromosomal rearrangements that juxtapose super-enhancers to the MYB locus. Chromosome conformation capture confirms that the translocated enhancers interact with the MYB promoter. Remarkably, MYB protein binds to the translocated enhancers, creating a positive feedback loop that sustains its expression. MYB also binds enhancers that drive different regulatory programs in alternate cell lineages in ACC, cooperating with TP63 in myoepithelial cells and a Notch program in luminal epithelial cells. Bromodomain inhibitors slow tumor growth in ACC primagraft models in vivo. Thus, our study identifies super-enhancer translocations that drive MYB expression and provides insight into downstream MYB functions in alternate ACC lineages.

Recurrent mutations affecting the histone H3.3 residues Lys27 or indirectly Lys36 are frequent drivers of pediatric high-grade gliomas (over 30% of HGGs). To identify additional driver mutations in HGGs, we investigated a cohort of 60 pediatric HGGs using whole-exome sequencing (WES) and compared them to 543 exomes from non-cancer control samples. We identified mutations in SETD2, a H3K36 trimethyltransferase, in 15% of pediatric HGGs, a result that was genome-wide significant (FDR = 0.029). Most SETD2 alterations were truncating mutations. Sequencing the gene in this cohort and another validation cohort (123 gliomas from all ages and grades) showed SETD2 mutations to be specific to high-grade tumors affecting 15% of pediatric HGGs (11/73) and 8% of adult HGGs (5/65) while no SETD2 mutations were identified in low-grade diffuse gliomas (0/45). Furthermore, SETD2 mutations were mutually exclusive with H3F3A mutations in HGGs (P = 0.0492) while they partly overlapped with IDH1 mutations (4/14), and SETD2-mutant tumors were found exclusively in the cerebral hemispheres (P = 0.0055). SETD2 is the only H3K36 trimethyltransferase in humans, and SETD2-mutant tumors showed a substantial decrease in H3K36me3 levels (P < 0.001), indicating that the mutations are loss-of-function. These data suggest that loss-of-function SETD2 mutations occur in older children and young adults and are specific to HGG of the cerebral cortex, similar to the H3.3 G34R/V and IDH mutations. Taken together, our results suggest that mutations disrupting the histone code at H3K36, including H3.3 G34R/V, IDH1 and/or SETD2 mutations, are central to the genesis of hemispheric HGGs in older children and young adults.

Understanding the principles governing mammalian gene regulation has been hampered by the difficulty in measuring in vivo binding dynamics of large numbers of transcription factors (TF) to DNA. Here, we develop a high-throughput Chromatin ImmunoPrecipitation (HT-ChIP) method to systematically map protein-DNA interactions. HT-ChIP was applied to define the dynamics of DNA binding by 25 TFs and 4 chromatin marks at 4 time-points following pathogen stimulus of dendritic cells. Analyzing over 180,000 TF-DNA interactions we find that TFs vary substantially in their temporal binding landscapes. This data suggests a model for transcription regulation whereby TF networks are hierarchically organized into cell differentiation factors, factors that bind targets prior to stimulus to prime them for induction, and factors that regulate specific gene programs. Overlaying HT-ChIP data on gene-expression dynamics shows that many TF-DNA interactions are established prior to the stimuli, predominantly at immediate-early genes, and identified specific TF ensembles that coordinately regulate gene-induction.

Reprogramming of cellular metabolism is a key event during tumorigenesis. Despite being known for decades (Warburg effect), the molecular mechanisms regulating this switch remained unexplored. Here, we identify SIRT6 as a tumor suppressor that regulates aerobic glycolysis in cancer cells. Importantly, loss of SIRT6 leads to tumor formation without activation of known oncogenes, whereas transformed SIRT6-deficient cells display increased glycolysis and tumor growth, suggesting that SIRT6 plays a role in both establishment and maintenance of cancer. By using a conditional SIRT6 allele, we show that SIRT6 deletion in vivo increases the number, size, and aggressiveness of tumors. SIRT6 also functions as a regulator of ribosome metabolism by corepressing MYC transcriptional activity. Lastly, Sirt6 is selectively downregulated in several human cancers, and expression levels of SIRT6 predict prognosis and tumor-free survival rates, highlighting SIRT6 as a critical modulator of cancer metabolism. Our studies reveal SIRT6 to be a potent tumor suppressor acting to suppress cancer metabolism.

Metabolic adaptation is essential for cell survival during nutrient deprivation. We report that eukaryotic elongation factor 2 kinase (eEF2K), which is activated by AMP-kinase (AMPK), confers cell survival under acute nutrient depletion by blocking translation elongation. Tumor cells exploit this pathway to adapt to nutrient deprivation by reactivating the AMPK-eEF2K axis. Adaptation of transformed cells to nutrient withdrawal is severely compromised in cells lacking eEF2K. Moreover, eEF2K knockdown restored sensitivity to acute nutrient deprivation in highly resistant human tumor cell lines. In vivo, overexpression of eEF2K rendered murine tumors remarkably resistant to caloric restriction. Expression of eEF2K strongly correlated with overall survival in human medulloblastoma and glioblastoma multiforme. Finally, C. elegans strains deficient in efk-1, the eEF2K ortholog, were severely compromised in their response to nutrient depletion. Our data highlight a conserved role for eEF2K in protecting cells from nutrient deprivation and in conferring tumor cell adaptation to metabolic stress. PAPERCLIP:

Glioblastoma (GBM) is thought to be driven by a subpopulation of cancer stem cells (CSCs) that self-renew and recapitulate tumor heterogeneity yet remain poorly understood. Here, we present a comparative analysis of chromatin state in GBM CSCs that reveals widespread activation of genes normally held in check by Polycomb repressors. These activated targets include a large set of developmental transcription factors (TFs) whose coordinated activation is unique to the CSCs. We demonstrate that a critical factor in the set, ASCL1, activates Wnt signaling by repressing the negative regulator DKK1. We show that ASCL1 is essential for the maintenance and in vivo tumorigenicity of GBM CSCs. Genome-wide binding profiles for ASCL1 and the Wnt effector LEF-1 provide mechanistic insight and suggest widespread interactions between the TF module and the signaling pathway. Our findings demonstrate regulatory connections among ASCL1, Wnt signaling, and collaborating TFs that are essential for the maintenance and tumorigenicity of GBM CSCs.

A report on the genomics workshop 'Identification of Functional Elements in Mammalian Genomes', Cold Spring Harbor, New York, 11-13 November 2004.

SOX9 is a master transcription factor that regulates development and stem cell programs. However, its potential oncogenic activity and regulatory mechanisms that control SOX9 protein stability are poorly understood. Here, we show that SOX9 is a substrate of FBW7, a tumor suppressor, and a SCF (SKP1/CUL1/F-box)-type ubiquitin ligase. FBW7 recognizes a conserved degron surrounding threonine 236 (T236) in SOX9 that is phosphorylated by GSK3 kinase and consequently degraded by SCFFBW7α Failure to degrade SOX9 promotes migration, metastasis, and treatment resistance in medulloblastoma, one of the most common childhood brain tumors. FBW7 is either mutated or downregulated in medulloblastoma, and in cases where FBW7 mRNA levels are low, SOX9 protein is significantly elevated and this phenotype is associated with metastasis at diagnosis and poor patient outcome. Transcriptional profiling of medulloblastoma cells expressing a degradation-resistant SOX9 mutant reveals activation of pro-metastatic genes and genes linked to cisplatin resistance. Finally, we show that pharmacological inhibition of PI3K/AKT/mTOR pathway activity destabilizes SOX9 in a GSK3/FBW7-dependent manner, rendering medulloblastoma cells sensitive to cytostatic treatment.

In cancer, recurrent somatic single-nucleotide variants-which are rare in most paediatric cancers-are confined largely to protein-coding genes1-3. Here we report highly recurrent hotspot mutations (r.3A>G) of U1 spliceosomal small nuclear RNAs (snRNAs) in about 50% of Sonic hedgehog (SHH) medulloblastomas. These mutations were not present across other subgroups of medulloblastoma, and we identified these hotspot mutations in U1 snRNA in only <0.1% of 2,442 cancers, across 36 other tumour types. The mutations occur in 97% of adults (subtype SHHδ) and 25% of adolescents (subtype SHHα) with SHH medulloblastoma, but are largely absent from SHH medulloblastoma in infants. The U1 snRNA mutations occur in the 5' splice-site binding region, and snRNA-mutant tumours have significantly disrupted RNA splicing and an excess of 5' cryptic splicing events. Alternative splicing mediated by mutant U1 snRNA inactivates tumour-suppressor genes (PTCH1) and activates oncogenes (GLI2 and CCND2), and represents a target for therapy. These U1 snRNA mutations provide an example of highly recurrent and tissue-specific mutations of a non-protein-coding gene in cancer.

H3K27M mutations resulting in epigenetic dysfunction are frequently observed in diffuse intrinsic pontine glioma (DIPGs), an incurable pediatric cancer. We conduct a CRISPR screen revealing that knockout of KDM1A encoding lysine-specific demethylase 1 (LSD1) sensitizes DIPG cells to histone deacetylase (HDAC) inhibitors. Consistently, Corin, a bifunctional inhibitor of HDACs and LSD1, potently inhibits DIPG growth in vitro and in xenografts. Mechanistically, Corin increases H3K27me3 levels suppressed by H3K27M histones, and simultaneously increases HDAC-targeted H3K27ac and LSD1-targeted H3K4me1 at differentiation-associated genes. Corin treatment induces cell death, cell-cycle arrest, and a cellular differentiation phenotype and drives transcriptional changes correlating with increased survival time in DIPG patients. These data suggest a strategy for treating DIPG by simultaneously inhibiting LSD1 and HDACs.

Development of effective targeted cancer therapies is fundamentally limited by our molecular understanding of disease pathogenesis. Diffuse intrinsic pontine glioma (DIPG) is a fatal malignancy of the childhood pons characterized by a unique substitution to methionine in histone H3 at lysine 27 (H3K27M) that results in globally altered epigenetic marks and oncogenic transcription. Through primary DIPG tumor characterization and isogenic oncohistone expression, we show that the same H3K27M mutation displays distinct modes of oncogenic reprogramming and establishes distinct enhancer architecture depending upon both the variant of histone H3 and the cell context in which the mutation occurs. Compared with non-malignant pediatric pontine tissue, we identify and functionally validate both shared and variant-specific pathophysiology. Altogether, we provide a powerful resource of epigenomic data in 25 primary DIPG samples and 5 rare normal pediatric pontine tissue samples, revealing clinically relevant functional distinctions previously unidentified in DIPG.

Most pancreatic neuroendocrine tumors (PNETs) do not produce excess hormones and are therefore considered 'non-functional'1-3. As clinical behaviors vary widely and distant metastases are eventually lethal2,4, biological classifications might guide treatment. Using enhancer maps to infer gene regulatory programs, we find that non-functional PNETs fall into two major subtypes, with epigenomes and transcriptomes that partially resemble islet α- and β-cells. Transcription factors ARX and PDX1 specify these normal cells, respectively5,6, and 84% of 142 non-functional PNETs expressed one or the other factor, occasionally both. Among 103 cases, distant relapses occurred almost exclusively in patients with ARX+PDX1- tumors and, within this subtype, in cases with alternative lengthening of telomeres. These markedly different outcomes belied similar clinical presentations and histology and, in one cohort, occurred irrespective of MEN1 mutation. This robust molecular stratification provides insight into cell lineage correlates of non-functional PNETs, accurately predicts disease course and can inform postoperative clinical decisions.

Acute myeloid leukemia (AML) is a heterogeneous disease that resides within a complex microenvironment, complicating efforts to understand how different cell types contribute to disease progression. We combined single-cell RNA sequencing and genotyping to profile 38,410 cells from 40 bone marrow aspirates, including 16 AML patients and five healthy donors. We then applied a machine learning classifier to distinguish a spectrum of malignant cell types whose abundances varied between patients and between subclones in the same tumor. Cell type compositions correlated with prototypic genetic lesions, including an association of FLT3-ITD with abundant progenitor-like cells. Primitive AML cells exhibited dysregulated transcriptional programs with co-expression of stemness and myeloid priming genes and had prognostic significance. Differentiated monocyte-like AML cells expressed diverse immunomodulatory genes and suppressed T cell activity in vitro. In conclusion, we provide single-cell technologies and an atlas of AML cell states, regulators, and markers with implications for precision medicine and immune therapies. VIDEO ABSTRACT.

High morbidity and mortality are common traits of malignant tumours and identification of the cells responsible is a focus of on-going research. Many studies are now reporting the use of antibodies specific to Clusters of Differentiation (CD) cell surface antigens to identify tumour-initiating cell (TIC) populations in neural tumours. Medulloblastoma is one of the most common malignant brain tumours in children and despite a considerable amount of research investigating this tumour, the identity of the TICs, and the means by which such cells can be targeted remain largely unknown. Current prognostication and stratification of medulloblastoma using clinical factors, histology and genetic profiling have classified this tumour into four main subgroups: WNT, Sonic hedgehog (SHH), Group 3 and Group 4. Of these subgroups, SHH remains one of the most studied tumour groups due to the ability to model medulloblastoma formation through targeted deletion of the Shh pathway inhibitor Patched1 (Ptch1). Here we sought to utilise CD antibody expression to identify and isolate TIC populations in Ptch1 deleted medulloblastoma, and determine if these antibodies can help classify the identity of human medulloblastoma subgroups. Using a fluorescence-activated cell sorted (FACS) CD antibody panel, we identified CD24 as a marker of TICs in Ptch1 deleted medulloblastoma. CD24 expression was not correlated with markers of astrocytes or oligodendrocytes, but co-labelled with markers of neural progenitor cells. In conjunction with CD15, proliferating CD24+/CD15+ granule cell precursors (GCPs) were identified as a TIC population in Ptch1 deleted medulloblastoma. On human medulloblastoma, CD24 was found to be highly expressed on Group 3, Group 4 and SHH subgroups compared with the WNT subgroup, which was predominantly positive for CD15, suggesting CD24 is an important marker of non-WNT medulloblastoma initiating cells and a potential therapeutic target in human medulloblastoma. This study reports the use of CD24 and CD15 to isolate a GCP-like TIC population in Ptch1 deleted medulloblastoma, and suggests CD24 expression as a marker to help stratify human WNT tumours from other medulloblastoma subgroups.

The HERBY trial was a phase II open-label, randomized, multicenter trial evaluating bevacizumab (BEV) in addition to temozolomide/radiotherapy in patients with newly diagnosed non-brainstem high-grade glioma (HGG) between the ages of 3 and 18 years. We carried out comprehensive molecular analysis integrated with pathology, radiology, and immune profiling. In post-hoc subgroup analysis, hypermutator tumors (mismatch repair deficiency and somatic POLE/POLD1 mutations) and those biologically resembling pleomorphic xanthoastrocytoma ([PXA]-like, driven by BRAF_V600E or NF1 mutation) had significantly more CD8+ tumor-infiltrating lymphocytes, and longer survival with the addition of BEV. Histone H3 subgroups (hemispheric G34R/V and midline K27M) had a worse outcome and were immune cold. Future clinical trials will need to take into account the diversity represented by the term "HGG" in the pediatric population.

Pilocytic astrocytoma is the most common childhood brain tumor, characterized by constitutive MAPK activation. MAPK signaling induces oncogene-induced senescence (OIS), which may cause unpredictable growth behavior of pilocytic astrocytomas. The senescence-associated secretory phenotype (SASP) has been shown to regulate OIS, but its role in pilocytic astrocytoma remains unknown.Experimental Design: The patient-derived pilocytic astrocytoma cell culture model, DKFZ-BT66, was used to demonstrate presence of the SASP and analyze its impact on OIS in pilocytic astrocytoma. The model allows for doxycycline-inducible switching between proliferation and OIS. Both states were studied using gene expression profiling (GEP), Western blot, ELISA, and cell viability testing. Primary pilocytic astrocytoma tumors were analyzed by GEP and multiplex assay.