RepEnTools: an automated repeat enrichment analysis package for ChIP-seq data reveals hUHRF1 Tandem-Tudor domain enrichment in young repeats

BACKGROUND

Repeat elements (REs) play important roles for cell function in health and disease. However, RE enrichment analysis in short-read high-throughput sequencing (HTS) data, such as ChIP-seq, is a challenging task.

RESULTS

Here, we present RepEnTools, a software package for genome-wide RE enrichment analysis of ChIP-seq and similar chromatin pulldown experiments. Our analysis package bundles together various software with carefully chosen and validated settings to provide a complete solution for RE analysis, starting from raw input files to tabular and graphical outputs. RepEnTools implementations are easily accessible even with minimal IT skills (Galaxy/UNIX). To demonstrate the performance of RepEnTools, we analysed chromatin pulldown data by the human UHRF1 TTD protein domain and discovered enrichment of TTD binding on young primate and hominid specific polymorphic repeats (SVA, L1PA1/L1HS) overlapping known enhancers and decorated with H3K4me1-K9me2/3 modifications. We corroborated these new bioinformatic findings with experimental data by qPCR assays using newly developed primate and hominid specific qPCR assays which complement similar research tools. Finally, we analysed mouse UHRF1 ChIP-seq data with RepEnTools and showed that the endogenous mUHRF1 protein colocalizes with H3K4me1-H3K9me3 on promoters of REs which were silenced by UHRF1. These new data suggest a functional role for UHRF1 in silencing of REs that is mediated by TTD binding to the H3K4me1-K9me3 double mark and conserved in two mammalian species.

CONCLUSIONS

RepEnTools improves the previously available programmes for RE enrichment analysis in chromatin pulldown studies by leveraging new tools, enhancing accessibility and adding some key functions. RepEnTools can analyse RE enrichment rapidly, efficiently, and accurately, providing the community with an up-to-date, reliable and accessible tool for this important type of analysis.

Refined read-out: The hUHRF1 Tandem-Tudor domain prefers binding to histone H3 tails containing K4me1 in the context of H3K9me2/3

UHRF1 is an essential chromatin protein required for DNA methylation maintenance, mammalian development, and gene regulation. We investigated the Tandem-Tudor domain (TTD) of human UHRF1 that is known to bind H3K9me2/3 histones and is a major driver of UHRF1 localization in cells. We verified binding to H3K9me2/3 but unexpectedly discovered stronger binding to H3 peptides and mononucleosomes containing K9me2/3 with additional K4me1. We investigated the combined binding of TTD to H3K4me1-K9me2/3 versus H3K9me2/3 alone, engineered mutants with specific and differential changes of binding, and discovered a novel read-out mechanism for H3K4me1 in an H3K9me2/3 context that is based on the interaction of R207 with the H3K4me1 methyl group and on counting the H-bond capacity of H3K4. Individual TTD mutants showed up to a 10,000-fold preference for the double-modified peptides, suggesting that after a conformational change, WT TTD could exhibit similar effects. The frequent appearance of H3K4me1-K9me2 regions in human chromatin demonstrated in our TTD chromatin pull-down and ChIP-western blot data suggests that it has specific biological roles. Chromatin pull-down of TTD from HepG2 cells and full-length murine UHRF1 ChIP-seq data correlate with H3K4me1 profiles indicating that the H3K4me1-K9me2/3 interaction of TTD influences chromatin binding of full-length UHRF1. We demonstrate the H3K4me1-K9me2/3 specific binding of UHRF1-TTD to enhancers and promoters of cell-type-specific genes at the flanks of cell-type-specific transcription factor binding sites, and provided evidence supporting an H3K4me1-K9me2/3 dependent and TTD mediated downregulation of these genes by UHRF1. All these findings illustrate the important physiological function of UHRF1-TTD binding to H3K4me1-K9me2/3 double marks in a cellular context.

The MECP2-TRD domain interacts with the DNMT3A-ADD domain at the H3-tail binding site

The DNMT3A DNA methyltransferase and MECP2 methylation reader are highly expressed in neurons. Both proteins interact via their DNMT3A-ADD and MECP2-TRD domains, and the MECP2 interaction regulates the activity and subnuclear localization of DNMT3A. Here, we mapped the interface of both domains using peptide SPOT array binding, protein pull-down, equilibrium peptide binding assays, and structural analyses. The region D529-D531 on the surface of the ADD domain was identified as interaction point with the TRD domain. This includes important residues of the histone H3 N-terminal tail binding site to the ADD domain, explaining why TRD and H3 binding to the ADD domain is competitive. On the TRD domain, residues 214-228 containing K219 and K223 were found to be essential for the ADD interaction. This part represents a folded patch within the otherwise largely disordered TRD domain. A crystal structure analysis of ADD revealed that the identified H3/TDR lysine binding pocket is occupied by an arginine residue from a crystallographic neighbor in the ADD apoprotein structure. Finally, we show that mutations in the interface of ADD and TRD domains disrupt the cellular interaction of both proteins in NIH3T3 cells. In summary, our data show that the H3 peptide binding cleft of the ADD domain also mediates the interaction with the MECP2-TRD domain suggesting that this binding site may have a broader role also in the interaction of DNMT3A with other proteins leading to complex regulation options by competitive and PTM specific binding.

Mechanistic basis of the increased methylation activity of the SETD2 protein lysine methyltransferase towards a designed super-substrate peptide

Protein lysine methyltransferases have important regulatory functions in cells, but mechanisms determining their activity and specificity are incompletely understood. Naturally, SETD2 introduces H3K36me3, but previously an artificial super-substrate (ssK36) was identified, which is methylated >100-fold faster. The ssK36-SETD2 complex structure cannot fully explain this effect. We applied molecular dynamics (MD) simulations and biochemical experiments to unravel the mechanistic basis of the increased methylation of ssK36, considering peptide conformations in solution, association of peptide and enzyme, and formation of transition-state (TS) like conformations of the enzyme-peptide complex. We observed in MD and FRET experiments that ssK36 adopts a hairpin conformation in solution with V35 and K36 placed in the loop. The hairpin conformation has easier access into the active site of SETD2 and it unfolds during the association process. Peptide methylation experiments revealed that introducing a stable hairpin conformation in the H3K36 peptide increased its methylation by SETD2. In MD simulations of enzyme-peptide complexes, the ssK36 peptide approached TS-like structures more frequently than H3K36 and distinct, substrate-specific TS-like structures were observed. Hairpin association, hairpin unfolding during association, and substrate-specific catalytically competent conformations may also be relevant for other PKMTs and hairpins could represent a promising starting point for SETD2 inhibitor development.

A functional LSD1 coregulator screen reveals a novel transcriptional regulatory cascade connecting R-loop homeostasis with epigenetic regulation

The lysine specific demethylase 1 (LSD1) plays a pivotal role in cellular differentiation by regulating the expression of key developmental genes in concert with different coregulatory proteins. This process is impaired in different cancer types and incompletely understood. To comprehensively identify functional coregulators of LSD1, we established a novel tractable fluorescent reporter system to monitor LSD1 activity in living cells. Combining this reporter system with a state-of-the-art multiplexed RNAi screen, we identify the DEAD-box helicase 19A (DDX19A) as a novel coregulator and demonstrate that suppression of Ddx19a results in an increase of R-loops and reduced LSD1-mediated gene silencing. We further show that DDX19A binds to tri-methylated lysine 27 of histone 3 (H3K27me3) and it regulates gene expression through the removal of transcription promoting R-loops. Our results uncover a novel transcriptional regulatory cascade where the downregulation of genes is dependent on the LSD1 mediated demethylation of histone H3 lysine 4 (H3K4). This allows the polycomb repressive complex 2 (PRC2) to methylate H3K27, which serves as a binding site for DDX19A. Finally, the binding of DDX19A leads to the efficient removal of R-loops at active promoters, which further de-represses LSD1 and PRC2, establishing a positive feedback loop leading to a robust repression of the target gene.

H3K36me2/3 Binding and DNA Binding of the DNA Methyltransferase DNMT3A PWWP Domain Both Contribute to its Chromatin Interaction

The PWWP domain of DNMT3 DNA methyltransferases binds to histone H3 tails containing methylated K36, and this activity is important for heterochromatic targeting. Here, we show that the PWWP domain of mouse DNMT3A binds to H3K36me2 and H3K36me3 with a slight preference for H3K36me2. PWWP domains have also been reported to bind to DNA, and the close proximity of H3K36 and nucleosomal DNA suggests a combined binding to H3K36me2/3 and DNA. We show here that the DNMT3A PWWP domain binds to DNA with a weak preference for AT-rich sequences and that the designed charge reversal R362E mutation disrupts DNA binding. The K295E mutation, as well as K295I recently identified in paraganglioma, a rare neuroendocrine neoplasm, disrupts both DNA and H3K36me2/3 binding, which is in agreement with the proximity of K295 to residues involved in K36me2/3 methyllysine binding. Nucleosome pulldown experiments show that DNA binding and H3K36me2/3 binding are important for the interaction of the DNMT3A PWWP domain with nucleosomes. Localization studies of transiently transfected fluorescently-tagged wild-type and PWWP-mutated full-length DNMT3A indicate that both interactions contribute to the subnuclear localization of DNMT3A in mouse cells. In summary, our data demonstrate that the combined binding of the DNMT3A PWWP domain to the H3 tail containing K36me2/3 and to the nucleosomal or linker DNA is important for its chromatin interaction and subnuclear targeting of DNMT3A in living cells.