High-throughput robust single-cell DNA methylation profiling with sciMETv2.

DNA methylation is a key epigenetic property that drives gene regulatory programs in development and disease. Current single-cell methods that produce high quality methylomes are expensive and low throughput without the aid of extensive automation. We previously described a proof-of-principle technique that enabled high cell throughput; however, it produced only low-coverage profiles and was a difficult protocol that required custom sequencing primers and recipes and frequently produced libraries with excessive adapter contamination. Here, we describe a greatly improved version that generates high-coverage profiles (~15-fold increase) using a robust protocol that does not require custom sequencing capabilities, includes multiple stopping points, and exhibits minimal adapter contamination. We demonstrate two versions of sciMETv2 on primary human cortex, a high coverage and rapid version, identifying distinct cell types using CH methylation patterns. These datasets are able to be directly integrated with one another as well as with existing snmC-seq2 datasets with little discernible bias. Finally, we demonstrate the ability to determine cell types using CG methylation alone, which is the dominant context for DNA methylation in most cell types other than neurons and the most applicable analysis outside of brain tissue.

NELF Regulates a Promoter-Proximal Step Distinct from RNA Pol II Pause-Release.

RNA polymerase II (RNA Pol II) is generally paused at promoter-proximal regions in most metazoans, and based on in vitro studies, this function has been attributed to the negative elongation factor (NELF). Here, we show that upon rapid depletion of NELF, RNA Pol II fails to be released into gene bodies, stopping instead around the +1 nucleosomal dyad-associated region. The transition to the 2nd pause region is independent of positive transcription elongation factor P-TEFb. During the heat shock response, RNA Pol II is rapidly released from pausing at heat shock-induced genes, while most genes are paused and transcriptionally downregulated. Both of these aspects of the heat shock response remain intact upon NELF loss. We find that NELF depletion results in global loss of cap-binding complex from chromatin without global reduction of nascent transcript 5' cap stability. Thus, our studies implicate NELF functioning in early elongation complexes distinct from RNA Pol II pause-release.

Targeting Processive Transcription Elongation via SEC Disruption for MYC-Induced Cancer Therapy.

The super elongation complex (SEC) is required for robust and productive transcription through release of RNA polymerase II (Pol II) with its P-TEFb module and promoting transcriptional processivity with its ELL2 subunit. Malfunction of SEC contributes to multiple human diseases including cancer. Here, we identify peptidomimetic lead compounds, KL-1 and its structural homolog KL-2, which disrupt the interaction between the SEC scaffolding protein AFF4 and P-TEFb, resulting in impaired release of Pol II from promoter-proximal pause sites and a reduced average rate of processive transcription elongation. SEC is required for induction of heat-shock genes and treating cells with KL-1 and KL-2 attenuates the heat-shock response from Drosophila to human. SEC inhibition downregulates MYC and MYC-dependent transcriptional programs in mammalian cells and delays tumor progression in a mouse xenograft model of MYC-driven cancer, indicating that small-molecule disruptors of SEC could be used for targeted therapy of MYC-induced cancer.

Structural Analysis of the Ash2L/Dpy-30 Complex Reveals a Heterogeneity in H3K4 Methylation.

Dpy-30 is a regulatory subunit controlling the histone methyltransferase activity of the KMT2 enzymes in vivo. Paradoxically, in vitro methyltransferase assays revealed that Dpy-30 only modestly participates in the positive heterotypic allosteric regulation of these methyltransferases. Detailed genome-wide, molecular and structural studies reveal that an extensive network of interactions taking place at the interface between Dpy-30 and Ash2L are critical for the correct placement, genome-wide, of H3K4me2 and H3K4me3 but marginally contribute to the methyltransferase activity of KMT2 enzymes in vitro. Moreover, we show that H3K4me2 peaks persisting following the loss of Dpy-30 are found in regions of highly transcribed genes, highlighting an interplay between Complex of Proteins Associated with SET1 (COMPASS) kinetics and the cycling of RNA polymerase to control H3K4 methylation. Overall, our data suggest that Dpy-30 couples its modest positive heterotypic allosteric regulation of KMT2 methyltransferase activity with its ability to help the positioning of SET1/COMPASS to control epigenetic signaling.

A high-throughput method to identify trans-activation domains within transcription factor sequences.

Even though transcription factors (TFs) are central players of gene regulation and have been extensively studied, their regulatory trans-activation domains (tADs) often remain unknown and a systematic functional characterization of tADs is lacking. Here, we present a novel high-throughput approach tAD-seq to functionally test thousands of candidate tADs from different TFs in parallel. The tADs we identify by pooled screening validate in individual luciferase assays, whereas neutral regions do not. Interestingly, the tADs are found at arbitrary positions within the TF sequences and can contain amino acid (e.g., glutamine) repeat regions or overlap structured domains, including helix-loop-helix domains that are typically annotated as DNA-binding. We also identified tADs in the non-native reading frames, confirming that random sequences can function as tADs, albeit weakly. The identification of tADs as short protein sequences sufficient for transcription activation will enable the systematic study of TF function, which-particularly for TFs of different transcription activating functionalities-is still poorly understood.

Resolving systematic errors in widely used enhancer activity assays in human cells.

The identification of transcriptional enhancers in the human genome is a prime goal in biology. Enhancers are typically predicted via chromatin marks, yet their function is primarily assessed with plasmid-based reporter assays. Here, we show that such assays are rendered unreliable by two previously reported phenomena relating to plasmid transfection into human cells: (i) the bacterial plasmid origin of replication (ORI) functions as a conflicting core promoter and (ii) a type I interferon (IFN-I) response is activated. These cause confounding false positives and negatives in luciferase assays and STARR-seq screens. We overcome both problems by employing the ORI as core promoter and by inhibiting two IFN-I-inducing kinases, enabling genome-wide STARR-seq screens in human cells. In HeLa-S3 cells, we uncover strong enhancers, IFN-I-induced enhancers, and enhancers endogenously silenced at the chromatin level. Our findings apply to all episomal enhancer activity assays in mammalian cells and are key to the characterization of human enhancers.

Therapeutic targeting of polycomb and BET bromodomain proteins in diffuse intrinsic pontine gliomas.

Diffuse intrinsic pontine glioma (DIPG) is a highly aggressive pediatric brainstem tumor characterized by rapid and uniform patient demise. A heterozygous point mutation of histone H3 occurs in more than 80% of these tumors and results in a lysine-to-methionine substitution (H3K27M). Expression of this histone mutant is accompanied by a reduction in the levels of polycomb repressive complex 2 (PRC2)-mediated H3K27 trimethylation (H3K27me3), and this is hypothesized to be a driving event of DIPG oncogenesis. Despite a major loss of H3K27me3, PRC2 activity is still detected in DIPG cells positive for H3K27M. To investigate the functional roles of H3K27M and PRC2 in DIPG pathogenesis, we profiled the epigenome of H3K27M-mutant DIPG cells and found that H3K27M associates with increased H3K27 acetylation (H3K27ac). In accordance with previous biochemical data, the majority of the heterotypic H3K27M-K27ac nucleosomes colocalize with bromodomain proteins at the loci of actively transcribed genes, whereas PRC2 is excluded from these regions; this suggests that H3K27M does not sequester PRC2 on chromatin. Residual PRC2 activity is required to maintain DIPG proliferative potential, by repressing neuronal differentiation and function. Finally, to examine the therapeutic potential of blocking the recruitment of bromodomain proteins by heterotypic H3K27M-K27ac nucleosomes in DIPG cells, we performed treatments in vivo with BET bromodomain inhibitors and demonstrate that they efficiently inhibit tumor progression, thus identifying this class of compounds as potential therapeutics in DIPG.

Therapeutic Targeting of MLL Degradation Pathways in MLL-Rearranged Leukemia.

Chromosomal translocations of the mixed-lineage leukemia (MLL) gene with various partner genes result in aggressive leukemia with dismal outcomes. Despite similar expression at the mRNA level from the wild-type and chimeric MLL alleles, the chimeric protein is more stable. We report that UBE2O functions in regulating the stability of wild-type MLL in response to interleukin-1 signaling. Targeting wild-type MLL degradation impedes MLL leukemia cell proliferation, and it downregulates a specific group of target genes of the MLL chimeras and their oncogenic cofactor, the super elongation complex. Pharmacologically inhibiting this pathway substantially delays progression, and it improves survival of murine leukemia through stabilizing wild-type MLL protein, which displaces the MLL chimera from some of its target genes and, therefore, relieves the cellular oncogenic addiction to MLL chimeras. Stabilization of MLL provides us with a paradigm in the development of therapies for aggressive MLL leukemia and perhaps for other cancers caused by translocations.

An Evolutionary Conserved Epigenetic Mark of Polycomb Response Elements Implemented by Trx/MLL/COMPASS.

Polycomb response elements (PREs) are specific DNA sequences that stably maintain the developmental pattern of gene expression. Drosophila PREs are well characterized, whereas the existence of PREs in mammals remains debated. Accumulating evidence supports a model in which CpG islands recruit Polycomb group (PcG) complexes; however, which subset of CGIs is selected to serve as PREs is unclear. Trithorax (Trx) positively regulates gene expression in Drosophila and co-occupies PREs to antagonize Polycomb-dependent silencing. Here we demonstrate that Trx-dependent H3K4 dimethylation (H3K4me2) marks Drosophila PREs and maintains the developmental expression pattern of nearby genes. Similarly, the mammalian Trx homolog, MLL1, deposits H3K4me2 at CpG-dense regions that could serve as PREs. In the absence of MLL1 and H3K4me2, H3K27me3 levels, a mark of Polycomb repressive complex 2 (PRC2), increase at these loci. By inhibiting PRC2-dependent H3K27me3 in the absence of MLL1, we can rescue expression of these loci, demonstrating a functional balance between MLL1 and PRC2 activities at these sites. Thus, our study provides rules for identifying cell-type-specific functional mammalian PREs within the human genome.

Regulation of the imprinted Dlk1-Dio3 locus by allele-specific enhancer activity.

Genomic imprinting is a critical developmental process characteristic of parent of origin-specific gene expression. It is well accepted that differentially DNA-methylated regions (DMRs) and enhancers are two major classes of cis-elements determining parent of origin-specific gene expression, with each recruiting different sets of transcription factors. Previously, we identified the AF4/FMR2 (AFF) family protein AFF3 within the transcription elongation complex SEC-L3. Here, we report that AFF3 can specifically bind both gametic DMRs (gDMRs) and enhancers within imprinted loci in an allele-specific manner. We identify the molecular regulators involved in the recruitment of AFF3 to gDMRs and provide mechanistic insight into the requirement of AFF3 at an enhancer for the expression of an ∼200-kb polycistronic transcript within the imprinted Dlk1-Dio3 locus. Our data suggest that the heterochromatic environment at the gDMR reinforces silencing of its related enhancer by controlling the binding and activity of AFF3 in an allele-specific manner. In summary, this study provides molecular details about the regulation of dosage-critical imprinted gene expression through the regulated binding of the transcription elongation factor AFF3 between a DMR and an enhancer.

Mitotic Transcriptional Activation: Clearance of Actively Engaged Pol II via Transcriptional Elongation Control in Mitosis.

Although it is established that some general transcription factors are inactivated at mitosis, many details of mitotic transcription inhibition (MTI) and its underlying mechanisms are largely unknown. We have identified mitotic transcriptional activation (MTA) as a key regulatory step to control transcription in mitosis for genes with transcriptionally engaged RNA polymerase II (Pol II) to activate and transcribe until the end of the gene to clear Pol II from mitotic chromatin, followed by global impairment of transcription reinitiation through MTI. Global nascent RNA sequencing and RNA fluorescence in situ hybridization demonstrate the existence of transcriptionally engaged Pol II in early mitosis. Both genetic and chemical inhibition of P-TEFb in mitosis lead to delays in the progression of cell division. Together, our study reveals a mechanism for MTA and MTI whereby transcriptionally engaged Pol II can progress into productive elongation and finish transcription to allow proper cellular division.

PAF1, a Molecular Regulator of Promoter-Proximal Pausing by RNA Polymerase II.

The control of promoter-proximal pausing and the release of RNA polymerase II (Pol II) is a widely used mechanism for regulating gene expression in metazoans, especially for genes that respond to environmental and developmental cues. Here, we identify that Pol-II-associated factor 1 (PAF1) possesses an evolutionarily conserved function in metazoans in the regulation of promoter-proximal pausing. Reduction in PAF1 levels leads to an increased release of paused Pol II into gene bodies at thousands of genes. PAF1 depletion results in increased nascent and mature transcripts and increased levels of phosphorylation of Pol II's C-terminal domain on serine 2 (Ser2P). These changes can be explained by the recruitment of the Ser2P kinase super elongation complex (SEC) effecting increased release of paused Pol II into productive elongation, thus establishing PAF1 as a regulator of promoter-proximal pausing by Pol II.

Context dependency of Set1/COMPASS-mediated histone H3 Lys4 trimethylation.

The stimulation of trimethylation of histone H3 Lys4 (H3K4) by H2B monoubiquitination (H2Bub) has been widely studied, with multiple mechanisms having been proposed for this form of histone cross-talk. Cps35/Swd2 within COMPASS (complex of proteins associated with Set1) is considered to bridge these different processes. However, a truncated form of Set1 (762-Set1) is reported to function in H3K4 trimethylation (H3K4me3) without interacting with Cps35/Swd2, and such cross-talk is attributed to the n-SET domain of Set1 and its interaction with the Cps40/Spp1 subunit of COMPASS. Here, we used biochemical, structural, in vivo, and chromatin immunoprecipitation (ChIP) sequencing (ChIP-seq) approaches to demonstrate that Cps40/Spp1 and the n-SET domain of Set1 are required for the stability of Set1 and not the cross-talk. Furthermore, the apparent wild-type levels of H3K4me3 in the 762-Set1 strain are due to the rogue methylase activity of this mutant, resulting in the mislocalization of H3K4me3 from the promoter-proximal regions to the gene bodies and intergenic regions. We also performed detailed screens and identified yeast strains lacking H2Bub but containing intact H2Bub enzymes that have normal levels of H3K4me3, suggesting that monoubiquitination may not directly stimulate COMPASS but rather works in the context of the PAF and Rad6/Bre1 complexes. Our study demonstrates that the monoubiquitination machinery and Cps35/Swd2 function to focus COMPASS's H3K4me3 activity at promoter-proximal regions in a context-dependent manner.