Superenhancer Analysis Defines Novel Epigenomic Subtypes of Non-APL AML, Including an RARα Dependency Targetable by SY-1425, a Potent and Selective RARα Agonist.

We characterized the enhancer landscape of 66 patients with acute myeloid leukemia (AML), identifying 6 novel subgroups and their associated regulatory loci. These subgroups are defined by their superenhancer (SE) maps, orthogonal to somatic mutations, and are associated with distinct leukemic cell states. Examination of transcriptional drivers for these epigenomic subtypes uncovers a subset of patients with a particularly strong SE at the retinoic acid receptor alpha (RARA) gene locus. The presence of a RARA SE and concomitant high levels of RARA mRNA predisposes cell lines and ex vivo models to exquisite sensitivity to a selective agonist of RARα, SY-1425 (tamibarotene). Furthermore, only AML patient-derived xenograft (PDX) models with high RARA mRNA were found to respond to SY-1425. Mechanistically, we show that the response to SY-1425 in RARA-high AML cells is similar to that of acute promyelocytic leukemia treated with retinoids, characterized by the induction of known retinoic acid response genes, increased differentiation, and loss of proliferation.Significance: We use the SE landscape of primary human AML to elucidate transcriptional circuitry and identify novel cancer vulnerabilities. A subset of patients were found to have an SE at RARA, which is predictive for response to SY-1425, a potent and selective RARα agonist, in preclinical models, forming the rationale for its clinical investigation in biomarker-selected patients. Cancer Discov; 7(10); 1136-53. ©2017 AACR.See related commentary by Wang and Aifantis, p. 1065.This article is highlighted in the In This Issue feature, p. 1047.

Noncoding Transcription Is a Driving Force for Nucleosome Instability in spt16 Mutant Cells.

FACT (facilitates chromatin transcription) consists of two essential subunits, Spt16 and Pob3, and functions as a histone chaperone. Mutation of spt16 results in a global loss of nucleosomes as well as aberrant transcription. Here, we show that the majority of nucleosome changes upon Spt16 depletion are alterations in nucleosome fuzziness and position shift. Most nucleosomal changes are suppressed by the inhibition of RNA polymerase II (Pol II) activity. Surprisingly, a small subgroup of nucleosome changes is resistant to transcriptional inhibition. Notably, Spt16 and distinct histone modifications are enriched at this subgroup of nucleosomes. We also report 1,037 Spt16-suppressed noncoding transcripts (SNTs) and found that the SNT start sites are enriched with the subgroup of nucleosomes resistant to Pol II inhibition. Finally, the nucleosomes at genes overlapping SNTs are more susceptible to changes upon Spt16 depletion than those without SNTs. Taken together, our results support a model in which Spt16 has a role in maintaining local nucleosome stability to inhibit initiation of SNT transcription, which once initiated drives additional nucleosome loss upon Spt16 depletion.

Alzheimer's loci: epigenetic associations and interaction with genetic factors.

OBJECTIVE: We explore the role of DNA methylation in Alzheimer's disease (AD). To elucidate where DNA methylation falls along the causal pathway linking risk factors to disease, we examine causal models to assess its role in the pathology of AD. METHODS: DNA methylation profiles were generated in 740 brain samples using the Illumina HumanMet450K beadset. We focused our analysis on CpG sites from 11 AD susceptibility gene regions. The primary outcome was a quantitative measure of neuritic amyloid plaque (NP), a key early element of AD pathology. We tested four causal models: (1) independent associations, (2) CpG mediating the association of a variant, (3) reverse causality, and (4) genetic variant by CpG interaction. RESULTS: Six genes regions (17 CpGs) showed evidence of CpG associations with NP, independent of genetic variation - BIN1 (5), CLU (5), MS4A6A (3), ABCA7 (2), CD2AP (1), and APOE (1). Together they explained 16.8% of the variability in NP. An interaction effect was seen in the CR1 region for two CpGs, cg10021878 (P = 0.01) and cg05922028 (P = 0.001), in relation to NP. In both cases, subjects with the risk allele rs6656401(AT) (/) (AA) display more methylation being associated with more NP burden, whereas subjects with the rs6656401(TT) protective genotype have an inverse association with more methylation being associated with less NP. INTERPRETATION: These observations suggest that, within known AD susceptibility loci, methylation is related to pathologic processes of AD and may play a largely independent role by influencing gene expression in AD susceptibility loci.

BRCA1 recruitment to transcriptional pause sites is required for R-loop-driven DNA damage repair.

The mechanisms contributing to transcription-associated genomic instability are both complex and incompletely understood. Although R-loops are normal transcriptional intermediates, they are also associated with genomic instability. Here, we show that BRCA1 is recruited to R-loops that form normally over a subset of transcription termination regions. There it mediates the recruitment of a specific, physiological binding partner, senataxin (SETX). Disruption of this complex led to R-loop-driven DNA damage at those loci as reflected by adjacent γ-H2AX accumulation and ssDNA breaks within the untranscribed strand of relevant R-loop structures. Genome-wide analysis revealed widespread BRCA1 binding enrichment at R-loop-rich termination regions (TRs) of actively transcribed genes. Strikingly, within some of these genes in BRCA1 null breast tumors, there are specific insertion/deletion mutations located close to R-loop-mediated BRCA1 binding sites within TRs. Thus, BRCA1/SETX complexes support a DNA repair mechanism that addresses R-loop-based DNA damage at transcriptional pause sites.

Integrative analysis of 111 reference human epigenomes.

The reference human genome sequence set the stage for studies of genetic variation and its association with human disease, but epigenomic studies lack a similar reference. To address this need, the NIH Roadmap Epigenomics Consortium generated the largest collection so far of human epigenomes for primary cells and tissues. Here we describe the integrative analysis of 111 reference human epigenomes generated as part of the programme, profiled for histone modification patterns, DNA accessibility, DNA methylation and RNA expression. We establish global maps of regulatory elements, define regulatory modules of coordinated activity, and their likely activators and repressors. We show that disease- and trait-associated genetic variants are enriched in tissue-specific epigenomic marks, revealing biologically relevant cell types for diverse human traits, and providing a resource for interpreting the molecular basis of human disease. Our results demonstrate the central role of epigenomic information for understanding gene regulation, cellular differentiation and human disease.

Alzheimer's disease: early alterations in brain DNA methylation at ANK1, BIN1, RHBDF2 and other loci.

We used a collection of 708 prospectively collected autopsied brains to assess the methylation state of the brain's DNA in relation to Alzheimer's disease (AD). We found that the level of methylation at 71 of the 415,848 interrogated CpGs was significantly associated with the burden of AD pathology, including CpGs in the ABCA7 and BIN1 regions, which harbor known AD susceptibility variants. We validated 11 of the differentially methylated regions in an independent set of 117 subjects. Furthermore, we functionally validated these CpG associations and identified the nearby genes whose RNA expression was altered in AD: ANK1, CDH23, DIP2A, RHBDF2, RPL13, SERPINF1 and SERPINF2. Our analyses suggest that these DNA methylation changes may have a role in the onset of AD given that we observed them in presymptomatic subjects and that six of the validated genes connect to a known AD susceptibility gene network.

Developmental control of gene copy number by repression of replication initiation and fork progression.

Precise DNA replication is crucial for genome maintenance, yet this process has been inherently difficult to study on a genome-wide level in untransformed differentiated metazoan cells. To determine how metazoan DNA replication can be repressed, we examined regions selectively under-replicated in Drosophila polytene salivary glands, and found they are transcriptionally silent and enriched for the repressive H3K27me3 mark. In the first genome-wide analysis of binding of the origin recognition complex (ORC) in a differentiated metazoan tissue, we find that ORC binding is dramatically reduced within these large domains, suggesting reduced initiation as one mechanism leading to under-replication. Inhibition of replication fork progression by the chromatin protein SUUR is an additional repression mechanism to reduce copy number. Although repressive histone marks are removed when SUUR is mutated and copy number restored, neither transcription nor ORC binding is reinstated. Tethering of the SUUR protein to a specific site is insufficient to block replication, however. These results establish that developmental control of DNA replication, at both the initiation and elongation stages, is a mechanism to change gene copy number during differentiation.

A cis-regulatory map of the Drosophila genome.

Systematic annotation of gene regulatory elements is a major challenge in genome science. Direct mapping of chromatin modification marks and transcriptional factor binding sites genome-wide has successfully identified specific subtypes of regulatory elements. In Drosophila several pioneering studies have provided genome-wide identification of Polycomb response elements, chromatin states, transcription factor binding sites, RNA polymerase II regulation and insulator elements; however, comprehensive annotation of the regulatory genome remains a significant challenge. Here we describe results from the modENCODE cis-regulatory annotation project. We produced a map of the Drosophila melanogaster regulatory genome on the basis of more than 300 chromatin immunoprecipitation data sets for eight chromatin features, five histone deacetylases and thirty-eight site-specific transcription factors at different stages of development. Using these data we inferred more than 20,000 candidate regulatory elements and validated a subset of predictions for promoters, enhancers and insulators in vivo. We identified also nearly 2,000 genomic regions of dense transcription factor binding associated with chromatin activity and accessibility. We discovered hundreds of new transcription factor co-binding relationships and defined a transcription factor network with over 800 potential regulatory relationships.

Chromatin signatures of the Drosophila replication program.

DNA replication initiates from thousands of start sites throughout the Drosophila genome and must be coordinated with other ongoing nuclear processes such as transcription to ensure genetic and epigenetic inheritance. Considerable progress has been made toward understanding how chromatin modifications regulate the transcription program; in contrast, we know relatively little about the role of the chromatin landscape in defining how start sites of DNA replication are selected and regulated. Here, we describe the Drosophila replication program in the context of the chromatin and transcription landscape for multiple cell lines using data generated by the modENCODE consortium. We find that while the cell lines exhibit similar replication programs, there are numerous cell line-specific differences that correlate with changes in the chromatin architecture. We identify chromatin features that are associated with replication timing, early origin usage, and ORC binding. Primary sequence, activating chromatin marks, and DNA-binding proteins (including chromatin remodelers) contribute in an additive manner to specify ORC-binding sites. We also generate accurate and predictive models from the chromatin data to describe origin usage and strength between cell lines. Multiple activating chromatin modifications contribute to the function and relative strength of replication origins, suggesting that the chromatin environment does not regulate origins of replication as a simple binary switch, but rather acts as a tunable rheostat to regulate replication initiation events.

Identification of functional elements and regulatory circuits by Drosophila modENCODE.

To gain insight into how genomic information is translated into cellular and developmental programs, the Drosophila model organism Encyclopedia of DNA Elements (modENCODE) project is comprehensively mapping transcripts, histone modifications, chromosomal proteins, transcription factors, replication proteins and intermediates, and nucleosome properties across a developmental time course and in multiple cell lines. We have generated more than 700 data sets and discovered protein-coding, noncoding, RNA regulatory, replication, and chromatin elements, more than tripling the annotated portion of the Drosophila genome. Correlated activity patterns of these elements reveal a functional regulatory network, which predicts putative new functions for genes, reveals stage- and tissue-specific regulators, and enables gene-expression prediction. Our results provide a foundation for directed experimental and computational studies in Drosophila and related species and also a model for systematic data integration toward comprehensive genomic and functional annotation.

WNT11 expression is induced by estrogen-related receptor alpha and beta-catenin and acts in an autocrine manner to increase cancer cell migration.

Elevated expression of the orphan nuclear receptor estrogen-related receptor α (ERRα) has been associated with a negative outcome in several cancers, although the mechanism(s) by which this receptor influences the pathophysiology of this disease and how its activity is regulated remain unknown. Using a chemical biology approach, it was determined that compounds, previously shown to inhibit canonical Wnt signaling, also inhibited the transcriptional activity of ERRα. The significance of this association was revealed in a series of biochemical and genetic experiments that show that (a) ERRα, ß-catenin (ß-cat), and lymphoid enhancer-binding factor-1 form macromolecular complexes in cells, (b) ERRα transcriptional activity is enhanced by ß-cat expression and vice versa, and (c) there is a high level of overlap among genes previously shown to be regulated by ERRα or ß-cat. Furthermore, silencing of ERRα and ß-cat expression individually or together dramatically reduced the migratory capacity of breast, prostate, and colon cancer cells in vitro. This increased migration could be attributed to the ERRα/ß-cat-dependent induction of WNT11. Specifically, using (a) conditioned medium from cells overexpressing recombinant WNT11 or (b) WNT11 neutralizing antibodies, we were able to show that this protein was the key mediator of the promigratory activities of ERRα/ß-cat. Together, these data provide evidence for an autocrine regulatory loop involving transcriptional upregulation of WNT11 by ERRα and ß-cat that influences the migratory capacity of cancer cells.

The conserved bromo-adjacent homology domain of yeast Orc1 functions in the selection of DNA replication origins within chromatin.

The origin recognition complex (ORC) binds to the specific positions on chromosomes that serve as DNA replication origins. Although ORC is conserved from yeast to humans, the DNA sequence elements that specify ORC binding are not. In particular, metazoan ORC shows no obvious DNA sequence specificity, whereas yeast ORC binds to a specific DNA sequence within all yeast origins. Thus, whereas chromatin must play an important role in metazoan ORC's ability to recognize origins, it is unclear whether chromatin plays a role in yeast ORC's recognition of origins. This study focused on the role of the conserved N-terminal bromo-adjacent homology domain of yeast Orc1 (Orc1BAH). Recent studies indicate that BAH domains are chromatin-binding modules. We show that the Orc1BAH domain was necessary for ORC's stable association with yeast chromosomes, and was physiologically relevant to DNA replication in vivo. This replication role was separable from the Orc1BAH domain's previously defined role in transcriptional silencing. Genome-wide analyses of ORC binding in ORC1 and orc1bahDelta cells revealed that the Orc1BAH domain contributed to ORC's association with most yeast origins, including a class of origins highly dependent on the Orc1BAH domain for ORC association (orc1bahDelta-sensitive origins). Orc1bahDelta-sensitive origins required the Orc1BAH domain for normal activity on chromosomes and plasmids, and were associated with a distinct local nucleosome structure. These data provide molecular insights into how the Orc1BAH domain contributes to ORC's selection of replication origins, as well as new tools for examining conserved mechanisms governing ORC's selection of origins within eukaryotic chromosomes.

Conserved nucleosome positioning defines replication origins.

The origin recognition complex (ORC) specifies replication origin location. The Saccharomyces cerevisiae ORC recognizes the ARS (autonomously replicating sequence) consensus sequence (ACS), but only a subset of potential genomic sites are bound, suggesting other chromosomal features influence ORC binding. Using high-throughput sequencing to map ORC binding and nucleosome positioning, we show that yeast origins are characterized by an asymmetric pattern of positioned nucleosomes flanking the ACS. The origin sequences are sufficient to maintain a nucleosome-free origin; however, ORC is required for the precise positioning of nucleosomes flanking the origin. These findings identify local nucleosomes as an important determinant for origin selection and function.

Multimodal regulation of E2F1 gene expression by progestins.

An analysis of mRNA expression in T47D breast cancer cells treated with the synthetic progestin R5020 revealed a subset of progesterone receptor (PR) target genes that are enriched for E2F binding sites. Following up on this observation, we determined that PR-B acts in both direct and indirect manners to positively upregulate E2F1 expression in T47D cells. The direct effects of PR on E2F1 expression were confirmed by chromatin immunoprecipitation (ChIP) analysis, which indicated that the agonist-bound receptor was recruited to several enhancer elements proximal to the E2F1 transcript. However, we also noted that cycloheximide partially inhibits R5020 induction of E2F1 expression, indicating that the ligand-dependent actions of PR on this gene may involve additional indirect regulatory pathways. In support of this hypothesis, we demonstrated that treatment with R5020 significantly increases both hyperphosphorylation of Rb and recruitment of E2F1 to its own promoter, thus activating a positive feedback loop that further amplifies its transcription. Furthermore, we established that PR-mediated induction of Krüppel-like factor 15 (KLF15), which can bind to GC-rich DNA within the E2F1 promoter, is required for maximal induction of E2F1 expression by progestins. Taken together, these results suggest a new paradigm for multimodal regulation of target gene expression by PR.