Structural Basis for the Versatile and Methylation-Dependent Binding of CTCF to DNA

One protein to rule them all: The role of CCCTC-binding factor in shaping human genome in health and disease

The eukaryotic genome, constituting several billion base pairs, must be contracted to fit within the volume of a nucleus where the diameter is on the scale of μm. The 3D structure and packing of such a long sequence cannot be left to pure chance, as DNA must be efficiently used for its primary roles as a matrix for transcription and replication. In recent years, methods like chromatin conformation capture (including 3C, 4C, Hi-C, ChIA-PET and Multi-ChIA) and optical microscopy have advanced substantially and have shed new light on how eukaryotic genomes are hierarchically organized; first into 10-nm fiber, next into DNA loops, topologically associated domains and finally into interphase or mitotic chromosomes. This knowledge has allowed us to revise our understanding regarding the mechanisms governing the process of DNA organization. Mounting experimental evidence suggests that the key element in the formation of loops is the binding of the CCCTC-binding factor (CTCF) to DNA; a protein that can be referred to as the chief organizer of the genome. However, CTCF does not work alone but in cooperation with other proteins, such as cohesin or Yin Yang 1 (YY1). In this short review, we briefly describe our current understanding of the structure of eukaryotic genomes, how they are established and how the formation of DNA loops can influence gene expression. We discuss the recent discoveries describing the 3D structure of the CTCF-DNA complex and the role of CTCF in establishing genome structure. Finally, we briefly explain how various genetic disorders might arise as a consequence of mutations in the CTCF target sequence or alteration of genomic imprinting.

A tour of 3D genome with a focus on CTCF

The complex three-dimensional (3D) structure of the genome plays critical roles in the maintenance of genome stability, organization, and dynamics and in regulation of gene expression for understanding molecular mechanisms and diseases. Chromatin maintains biological functions and transcriptional activities through long distance interaction and interactions between loops and enhancers-promoters. We firstly overview the architecture and biology of chromatin and loops, topologically associated domains (TADs) and interactions, and compartments and functions. We specifically focus on CCCTC-binding factor (CTCF) in 3D genome organization and function to furthermore understand the significance of CTCF biology, transcriptional regulations, interactions with cohesin, roles in DNA binding, influences of CTCF degradation, and communication with wings-apart like (Wapl) protein. We also summarize the advanced single cell approaches to further monitor dynamics of CTCF functions and structures in the maintenance of 3D genome organization and function at single cell level.

Identification of Cancer Drivers at CTCF Insulators in 1,962 Whole Genomes

Recent studies have shown that mutations at non-coding elements, such as promoters and enhancers, can act as cancer drivers. However, an important class of non-coding elements, namely CTCF insulators, has been overlooked in the previous driver analyses. We used insulator annotations from CTCF and cohesin ChIA-PET and analyzed somatic mutations in 1,962 whole genomes from 21 cancer types. Using the heterogeneous patterns of transcription-factor-motif disruption, functional impact, and recurrence of mutations, we developed a computational method that revealed 21 insulators showing signals of positive selection. In particular, mutations in an insulator in multiple cancer types, including 16% of melanoma samples, are associated with TGFB1 up-regulation. Using CRISPR-Cas9, we find that alterations at two of the most frequently mutated regions in this insulator increase cell growth by 40%-50%, supporting the role of this boundary element as a cancer driver. Thus, our study reveals several CTCF insulators as putative cancer drivers.

DNA methylation repels targeting of Arabidopsis REF6

RELATIVE OF EARLY FLOWERING 6 (REF6/JMJ12), a Jumonji C (JmjC)-domain-containing H3K27me3 histone demethylase, finds its target loci in Arabidopsis genome by directly recognizing the CTCTGYTY motif via its zinc-finger (ZnF) domains. REF6 tends to bind motifs located in active chromatin states that are depleted for heterochromatic modifications. However, the underlying mechanism remains unknown. Here, we show that REF6 preferentially bind to hypo-methylated CTCTGYTY motifs in vivo, and that CHG methylation decreases REF6 DNA binding affinity in vitro. In addition, crystal structures of ZnF-clusters in complex with DNA oligonucleotides reveal that 5-methylcytosine is unfavorable for REF6 binding. In drm1 drm2 cmt2 cmt3 (ddcc) quadruple mutants, in which non-CG methylation is significantly reduced, REF6 can ectopically bind a small number of new target loci, most of which are located in or neighbored with short TEs in euchromatic regions. Collectively, our findings reveal that DNA methylation, likely acting in combination with other epigenetic modifications, may partially explain why REF6 binding is depleted in heterochromatic loci.

REF6 is a H3K27me3 histone demethylase that targets regions of active chromatin in the Arabidopsis genome. Here Qiu et al. show that REF6 preferentially binds unmethylated DNA via its zinc-finger domains suggesting that DNA methylation might deter REF6 binding in heterochromatin.

DNA (de)methylation in embryonic stem cells controls CTCF-dependent chromatin boundaries

Coordinated changes of DNA (de)methylation, nucleosome positioning, and chromatin binding of the architectural protein CTCF play an important role for establishing cell-type-specific chromatin states during differentiation. To elucidate molecular mechanisms that link these processes, we studied the perturbed DNA modification landscape in mouse embryonic stem cells (ESCs) carrying a double knockout (DKO) of the Tet1 and Tet2 dioxygenases. These enzymes are responsible for the conversion of 5-methylcytosine (5mC) into its hydroxymethylated (5hmC), formylated (5fC), or carboxylated (5caC) forms. We determined changes in nucleosome positioning, CTCF binding, DNA methylation, and gene expression in DKO ESCs and developed biophysical models to predict differential CTCF binding. Methylation-sensitive nucleosome repositioning accounted for a significant portion of CTCF binding loss in DKO ESCs, whereas unmethylated and nucleosome-depleted CpG islands were enriched for CTCF sites that remained occupied. A number of CTCF sites also displayed direct correlations with the CpG modification state: CTCF was preferentially lost from sites that were marked with 5hmC in wild-type (WT) cells but not from 5fC-enriched sites. In addition, we found that some CTCF sites can act as bifurcation points defining the differential methylation landscape. CTCF loss from such sites, for example, at promoters, boundaries of chromatin loops, and topologically associated domains (TADs), was correlated with DNA methylation/demethylation spreading and can be linked to down-regulation of neighboring genes. Our results reveal a hierarchical interplay between cytosine modifications, nucleosome positions, and DNA sequence that determines differential CTCF binding and regulates gene expression.

An alternative CTCF isoform antagonizes canonical CTCF occupancy and changes chromatin architecture to promote apoptosis

CTCF plays key roles in gene regulation, chromatin insulation, imprinting, X chromosome inactivation and organizing the higher-order chromatin architecture of mammalian genomes. Previous studies have mainly focused on the roles of the canonical CTCF isoform. Here, we explore the functions of an alternatively spliced human CTCF isoform in which exons 3 and 4 are skipped, producing a shorter isoform (CTCF-s). Functionally, we find that CTCF-s competes with the genome binding of canonical CTCF and binds a similar DNA sequence. CTCF-s binding disrupts CTCF/cohesin binding, alters CTCF-mediated chromatin looping and promotes the activation of IFI6 that leads to apoptosis. This effect is caused by an abnormal long-range interaction at the IFI6 enhancer and promoter. Taken together, this study reveals a non-canonical function for CTCF-s that antagonizes the genomic binding of canonical CTCF and cohesin, and that modulates chromatin looping and causes apoptosis by stimulating IFI6 expression.

CTCF plays key roles in gene regulation, chromatin insulation and organizing the higher-order chromatin architecture of mammalian genomes. Here the authors investigate the function an alternatively spliced shorter CTCF isoform, finding that this isoform antagonizes canonical CTCF occupancy and changes chromatin architecture to promote apoptosis.

The anti-cancer drugs curaxins target spatial genome organization

Recently we characterized a class of anti-cancer agents (curaxins) that disturbs DNA/histone interactions within nucleosomes. Here, using a combination of genomic and in vitro approaches, we demonstrate that curaxins strongly affect spatial genome organization and compromise enhancer-promoter communication, which is necessary for the expression of several oncogenes, including MYC. We further show that curaxins selectively inhibit enhancer-regulated transcription of chromatinized templates in cell-free conditions. Genomic studies also suggest that curaxins induce partial depletion of CTCF from its binding sites, which contributes to the observed changes in genome topology. Thus, curaxins can be classified as epigenetic drugs that target the 3D genome organization.

Curaxins are a recently discovered class of anti-cancer agents that disturbs DNA/histone interactions within. Here the authors provide evidence that curaxins affect the spatial genome organization and compromise enhancer-promoter communication necessary for expression of several oncogenes, including MYC.

TET Enzymes and 5hmC in Adaptive and Innate Immune Systems

DNA methylation is an abundant and stable epigenetic modification that allows inheritance of information from parental to daughter cells. At active genomic regions, DNA methylation can be reversed by TET (Ten-eleven translocation) enzymes, which are responsible for fine-tuning methylation patterns. TET enzymes oxidize the methyl group of 5-methylcytosine (5mC) to yield 5-hydroxymethylcytosine (5hmC) and other oxidized methylcytosines, facilitating both passive and active demethylation. Increasing evidence has demonstrated the essential functions of TET enzymes in regulating gene expression, promoting cell differentiation, and suppressing tumor formation. In this review, we will focus on recent discoveries of the functions of TET enzymes in the development and function of lymphoid and myeloid cells. How TET activity can be modulated by metabolites, including vitamin C and 2-hydroxyglutarate, and its potential application in shaping the course of immune response will be discussed.

DNA methylation in mice is influenced by genetics as well as sex and life experience

DNA methylation is an essential epigenetic process in mammals, intimately involved in gene regulation. Here we address the extent to which genetics, sex, and pregnancy influence genomic DNA methylation by intercrossing 2 inbred mouse strains, C57BL/6N and C3H/HeN, and analyzing DNA methylation in parents and offspring using whole-genome bisulfite sequencing. Differential methylation across genotype is detected at thousands of loci and is preserved on parental alleles in offspring. In comparison of autosomal DNA methylation patterns across sex, hundreds of differentially methylated regions are detected. Comparison of animals with different histories of pregnancy within our study reveals a CpG methylation pattern that is restricted to female animals that had borne offspring. Collectively, our results demonstrate the stability of CpG methylation across generations, clarify the interplay of epigenetics with genetics and sex, and suggest that CpG methylation may serve as an epigenetic record of life events in somatic tissues at loci whose expression is linked to the relevant biology.

DNA methylation is an epigenetic mark involved in gene regulation. Here the authors investigate the extent to which genetics, sex and pregnancy influence genomic DNA methylation in mice, providing evidence of the stability of CpG methylation across generation and suggest that CpG methylation may serve as an epigenetic record of life events in somatic tissues at loci whose expression is linked to the relevant biology.

DNA Methylation and Transcription Factors Competitively Regulate SIRT4 Promoter Activity in Bovine Adipocytes: Roles of NRF1 and CMYB

Sirtuin 4 (SIRT4) belongs to the mitochondrial sirtuin protein family, a class of NAD⁺-dependent protein deacylases that remove post-translational acyl modifications from cellular substrates during the regulation of various biological pathways. SIRT4 has been shown to regulate lipid homeostasis. However, the mechanism by which the bovine SIRT4 gene is transcriptionally regulated remains unknown. To explore the molecular mechanism of SIRT4 expression, we obtained a 400-kb fragment of the 5'-regulatory region of bovine SIRT4 by molecular cloning, which contained a CpG island. Electrophoretic mobility shift assays and luciferase reporter gene assays identified the nuclear respiratory factor 1 (NRF1) and myb proto-oncogene protein (CMYB) binding sites as transcriptional repression and activation sites in the SIRT4 promoter region, respectively. We further verified that NRF1 and CMYB bind to the SIRT4 promoter using chromatin immunoprecipitation assays. In addition, from DNA methylation and reporter gene assays, results revealed that SIRT4 promoter activity was enhanced by demethylation. Further, NRF1-mediated transcriptional inhibition and CMYB-mediated transcriptional activation of SIRT4 expression were strengthened by demethylation during bovine adipocyte differentiation. Taken together, our results shed light on the mechanism underlying the promoter methylation and transcriptional regulation of SIRT4 expression in bovine adipocytes.

Molecular Processes Connecting DNA Methylation Patterns with DNA Methyltransferases and Histone Modifications in Mammalian Genomes

DNA methylation is an essential part of the epigenome chromatin modification network, which also comprises several covalent histone protein post-translational modifications. All these modifications are highly interconnected, because the writers and erasers of one mark, DNA methyltransferases (DNMTs) and ten eleven translocation enzymes (TETs) in the case of DNA methylation, are directly or indirectly targeted and regulated by other marks. Here, we have collected information about the genomic distribution and variability of DNA methylation in human and mouse DNA in different genomic elements. After summarizing the impact of DNA methylation on genome evolution including CpG depletion, we describe the connection of DNA methylation with several important histone post-translational modifications, including methylation of H3K4, H3K9, H3K27, and H3K36, but also with nucleosome remodeling. Moreover, we present the mechanistic features of mammalian DNA methyltransferases and their associated factors that mediate the crosstalk between DNA methylation and chromatin modifications. Finally, we describe recent advances regarding the methylation of non-CpG sites, methylation of adenine residues in human cells and methylation of mitochondrial DNA. At several places, we highlight controversial findings or open questions demanding future experimental work.

Structural insights into methylated DNA recognition by the C-terminal zinc fingers of the DNA reader protein ZBTB38

Methyl-CpG-binding proteins (MBPs) are selective readers of DNA methylation that play an essential role in mediating cellular transcription processes in both normal and diseased cells. This physiological function of MBPs has generated significant interest in understanding the mechanisms by which these proteins read and interpret DNA methylation signals. Zinc finger and BTB domain-containing 38 (ZBTB38) represents one member of the zinc finger (ZF) family of MBPs. We recently demonstrated that the C-terminal ZFs of ZBTB38 exhibit methyl-selective DNA binding within the ((A/G)TmCG(G/A)(mC/T)(G/A)) context both in vitro and within cells. Here we report the crystal structure of the first four C-terminal ZBTB38 ZFs (ZFs 6-9) in complex with the previously identified methylated consensus sequence at 1.75 Å resolution. From the structure, methyl-selective binding is preferentially localized at the 5' mCpG site of the bound DNA, which is facilitated through a series of base-specific interactions from residues within the α-helices of ZF7 and ZF8. ZF6 and ZF9 primarily stabilize ZF7 and ZF8 to facilitate the core base-specific interactions. Further structural and biochemical analyses, including solution NMR spectroscopy and electrophoretic mobility gel shift assays, revealed that the C-terminal ZFs of ZBTB38 utilize an alternative mode of mCpG recognition from the ZF MBPs structurally evaluated to date. Combined, these findings provide insight into the mechanism by which this ZF domain of ZBTB38 selectively recognizes methylated CpG sites and expands our understanding of how ZF-containing proteins can interpret this essential epigenetic mark.

Competitive regulation by transcription factors and DNA methylation in the bovine SIRT5 promoter: Roles of E2F4 and KLF6

Sirtuin 5 (SIRT5) belongs to the mitochondrial sirtuin family, which constitutes a highly conserved family of nicotinamide adenine dinucleotide NAD⁺-dependent deacetylases and ADP-ribosyltransferases that play important regulatory roles in stress resistance and metabolic homeostasis. SIRT5 was shown to have deacetylase, desuccinylase, and demalonylase activities. However, the mechanisms regulating SIRT5 transcription remain unclear. To explore the molecular regulation of bovine SIRT5 expression, we obtained a 500-base pair fragment of the 5'-regulatory region of bovine SIRT5 by molecular cloning, which contained a region with 3 CpG islands. Electrophoretic mobility shift assays and luciferase reporter assays revealed the E2F transcription factor 4 (E2F4) and Kruppel-like factor 6 (KLF6) binding sites as transcriptional activators or repressors in the promoter region of SIRT5. We further verified that E2F4 and KLF6 bind to the SIRT5 promoter by chromatin immunoprecipitation assays. Additionally, methylation and luciferase report assays showed that SIRT5 promoter activity was enhanced by demethylation, and transcriptional activation by E2F4 and transcriptional inhibition by KLF6 of SIRT5 expression was strengthened by demethylation during adipocytes differentiation. This study focused on the mechanism underlying the methylation and transcriptional regulation of SIRT5 expression in bovine adipocytes.

Base-pair resolution detection of transcription factor binding site by deep deconvolutional network

Motivation

Transcription factor (TF) binds to the promoter region of a gene to control gene expression. Identifying precise TF binding sites (TFBSs) is essential for understanding the detailed mechanisms of TF-mediated gene regulation. However, there is a shortage of computational approach that can deliver single base pair resolution prediction of TFBS.

Results

In this paper, we propose DeepSNR, a Deep Learning algorithm for predicting TF binding location at Single Nucleotide Resolution de novo from DNA sequence. DeepSNR adopts a novel deconvolutional network (deconvNet) model and is inspired by the similarity to image segmentation by deconvNet. The proposed deconvNet architecture is constructed on top of 'DeepBind' and we trained the entire model using TF-specific data from ChIP-exonuclease (ChIP-exo) experiments. DeepSNR has been shown to outperform motif search-based methods for several evaluation metrics. We have also demonstrated the usefulness of DeepSNR in the regulatory analysis of TFBS as well as in improving the TFBS prediction specificity using ChIP-seq data.

Availability and implementation

DeepSNR is available open source in the GitHub repository (https://github.com/sirajulsalekin/DeepSNR).

Supplementary information

Supplementary data are available at Bioinformatics online.

Zinc Finger Readers of Methylated DNA

DNA methylation is a prevalent epigenetic modification involved in regulating a number of essential cellular processes, including genomic accessibility and transcriptional outcomes. As such, aberrant alterations in global DNA methylation patterns have been associated with a growing number of disease conditions. Nevertheless, the full mechanisms by which DNA methylation information is interpreted and translated into genomic responses is not yet fully understood. Methyl-CpG binding proteins (MBPs) function as important mediators of this essential process by selectively reading DNA methylation signals and translating this information into down-stream cellular outcomes. The Cys₂His₂ zinc finger scaffold is one of the most abundant DNA binding motifs found within human transcription factors, yet only a few zinc finger containing proteins capable of conferring selectivity for mCpG over CpG sites have been characterized. This review summarizes our current structural understanding for the mechanisms by which the zinc finger MBPs evaluated to date read this essential epigenetic mark. Further, some of the biological implications for mCpG readout elicited by this family of MBPs are discussed.

DNA recognition patterns of the multi-zinc-finger protein CTCF: a mutagenesis study

CCCTC-binding factor (CTCF) is a zinc-finger protein, serving an important part in the genome architecture as well as some biochemical processes. Over 70,000 CTCF binding DNA sites have been detected genome-wide, and most anchors of chromatin loops are demarcated with the CTCF binding. Various protein or RNA molecules interact with DNA-bound CTCF to conduct different biological functions, and potentially the interfaces between CTCF and its cofactors can be targets for drug development. Here we identify the effective region of CTCF in DNA recognition, which defines the exposed CTCF surface feature for the interaction of cofactors. While the zinc-finger region contributes the most in DNA association, its binding affinity varies based on different DNA sequences. To investigate the effectiveness of individual zinc-fingers, the key residues are mutated to inactivate the DNA binding ability, while the finger configuration and the spacing between fingers are preserved. The strategy is proved to be successful, while clear differences are observed in the DNA binding affinities among the 11 finger mutants and the result is consistent to previous studies in general. With the help of inactivated finger mutants, we identify the ineffective fingers and the dominant effective fingers, which form distinctive patterns on different DNA targets.

Dynamic DNA methylation: In the right place at the right time

The classical model of cytosine DNA methylation (the presence of 5-methylcytosine, 5mC) regulation depicts this covalent modification as a stable repressive regulator of promoter activity. However, whole-genome analysis of 5mC reveals widespread tissue- and cell type-specific patterns and pervasive dynamics during mammalian development. Here we review recent findings that delineate 5mC functions in developmental stages and diverse genomic compartments as well as discuss the molecular mechanisms that connect transcriptional regulation and 5mC. Beyond the newly appreciated dynamics, regulatory roles for 5mC have been suggested in new biological contexts, such as learning and memory or aging. The use of new single-cell measurement techniques and precise editing tools will enable functional analyses of 5mC in gene expression, clarifying its role in various biological processes.

Dynamic Nature of CTCF Tandem 11 Zinc Fingers in Multivalent Recognition of DNA As Revealed by NMR Spectroscopy

The 11 zinc fingers (ZFs) of the transcription factor CTCF play a versatile role in the regulation of gene expression. CTCF binds to numerous genomic sites to form chromatin loops and topologically associated domains and thus mediates the 3D architecture of chromatin. Although CTCF inter-ZF plasticity is essential for the recognition of multiple genomic sites, the dynamic nature of its 11 ZFs remains unknown. We assigned the chemical shifts of the CTCF ZFs 1-11 and solved the solution structures of each ZF. NMR backbone dynamics, residual dipolar couplings, and small-angle X-ray scattering experiments suggest a high inter-ZF plasticity of the free-form ZFs 1-11. As exemplified by two different protocadherin DNA sequences, the titration of DNAs to ¹⁵N-labeled CTCF ZFs 1-11 enabled systematic mapping of binding of CTCF ZFs to various chromatin sites. Our work paves the way for illustrating the molecular basis of the versatile DNA recognized by CTCF and has interesting implications for its conformational transition during DNA binding.

Detecting and interpreting DNA methylation marks

The generation, alteration, recognition, and erasure of epigenetic modifications of DNA are fundamental to controlling gene expression in mammals. These covalent DNA modifications include cytosine methylation by AdoMet-dependent methyltransferases and 5-methylcytosine oxidation by Fe(II)-dependent and α-ketoglutarate-dependent dioxygenases. Sequence-specific transcription factors are responsible for interpreting the modification status of specific regions of chromatin. This review focuses on recent developments in characterizing the functional and structural links between the modification status of two DNA bases: 5-methylcytosine and 5-methyluracil (thymine).

An integrated -omics analysis of the epigenetic landscape of gene expression in human blood cells

Background

Gene expression can be influenced by DNA methylation 1) distally, at regulatory elements such as enhancers, as well as 2) proximally, at promoters. Our current understanding of the influence of distal DNA methylation changes on gene expression patterns is incomplete. Here, we characterize genome-wide methylation and expression patterns for ~ 13 k genes to explore how DNA methylation interacts with gene expression, throughout the genome.

Results

We used a linear mixed model framework to assess the correlation of DNA methylation at ~ 400 k CpGs with gene expression changes at ~ 13 k transcripts in two independent datasets from human blood cells. Among CpGs at which methylation significantly associates with transcription (eCpGs), > 50% are distal (> 50 kb) or trans (different chromosome) to the correlated gene. Many eCpG-transcript pairs are consistent between studies and ~ 90% of neighboring eCpGs associate with the same gene, within studies. We find that enhancers (P < 5e-18) and microRNA genes (P = 9e-3) are overrepresented among trans eCpGs, and insulators and long intergenic non-coding RNAs are enriched among cis and distal eCpGs. Intragenic-eCpG-transcript correlations are negative in 60–70% of occurrences and are enriched for annotated gene promoters and enhancers (P < 0.002), highlighting the importance of intragenic regulation. Gene Ontology analysis indicates that trans eCpGs are enriched for transcription factor genes and chromatin modifiers, suggesting that some trans eCpGs represent the influence of gene networks and higher-order transcriptional control.

Conclusions

This work sheds new light on the interplay between epigenetic changes and gene expression, and provides useful data for mining biologically-relevant results from epigenome-wide association studies.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4842-3) contains supplementary material, which is available to authorized users.

A modular dCas9-SunTag DNMT3A epigenome editing system overcomes pervasive off-target activity of direct fusion dCas9-DNMT3A constructs

Detection of DNA methylation in the genome has been possible for decades; however, the ability to deliberately and specifically manipulate local DNA methylation states in the genome has been extremely limited. Consequently, this has impeded our understanding of the direct effect of DNA methylation on transcriptional regulation and transcription factor binding in the native chromatin context. Thus, highly specific targeted epigenome editing tools are needed to address this. Recent adaptations of genome editing technologies, including fusion of the DNMT3A DNA methyltransferase catalytic domain to catalytically inactive Cas9 (dC9-D3A), have aimed to alter DNA methylation at desired loci. Here, we show that these tools exhibit consistent off-target DNA methylation deposition in the genome, limiting their capabilities to unambiguously assess the functional consequences of DNA methylation. To address this, we developed a modular dCas9-SunTag (dC9Sun-D3A) system that can recruit multiple DNMT3A catalytic domains to a target site for editing DNA methylation. dC9Sun-D3A is tunable, specific, and exhibits much higher induction of DNA methylation at target sites than the dC9-D3A direct fusion protein. Importantly, genome-wide characterization of dC9Sun-D3A binding sites and DNA methylation revealed minimal off-target protein binding and induction of DNA methylation with dC9Sun-D3A, compared to pervasive off-target methylation by dC9-D3A. Furthermore, we used dC9Sun-D3A to demonstrate the binding sensitivity to DNA methylation for CTCF and NRF1 in situ. Overall, this modular dC9Sun-D3A system enables precise DNA methylation deposition with the lowest off-target DNA methylation levels reported to date, allowing accurate functional determination of the role of DNA methylation at single loci.

Putative bovine topological association domains and CTCF binding motifs can reduce the search space for causative regulatory variants of complex traits

BACKGROUND

Topological association domains (TADs) are chromosomal domains characterised by frequent internal DNA-DNA interactions. The transcription factor CTCF binds to conserved DNA sequence patterns called CTCF binding motifs to either prohibit or facilitate chromosomal interactions. TADs and CTCF binding motifs control gene expression, but they are not yet well defined in the bovine genome. In this paper, we sought to improve the annotation of bovine TADs and CTCF binding motifs, and assess whether the new annotation can reduce the search space for cis-regulatory variants.

RESULTS

We used genomic synteny to map TADs and CTCF binding motifs from humans, mice, dogs and macaques to the bovine genome. We found that our mapped TADs exhibited the same hallmark properties of those sourced from experimental data, such as housekeeping genes, transfer RNA genes, CTCF binding motifs, short interspersed elements, H3K4me3 and H3K27ac. We showed that runs of genes with the same pattern of allele-specific expression (ASE) (either favouring paternal or maternal allele) were often located in the same TAD or between the same conserved CTCF binding motifs. Analyses of variance showed that when averaged across all bovine tissues tested, TADs explained 14% of ASE variation (standard deviation, SD: 0.056), while CTCF explained 27% (SD: 0.078). Furthermore, we showed that the quantitative trait loci (QTLs) associated with gene expression variation (eQTLs) or ASE variation (aseQTLs), which were identified from mRNA transcripts from 141 lactating cows' white blood and milk cells, were highly enriched at putative bovine CTCF binding motifs. The linearly-furthermost, and most-significant aseQTL and eQTL for each genic target were located within the same TAD as the gene more often than expected (Chi-Squared test P-value < 0.001).

CONCLUSIONS

Our results suggest that genomic synteny can be used to functionally annotate conserved transcriptional components, and provides a tool to reduce the search space for causative regulatory variants in the bovine genome.

Emerging Evidence of Chromosome Folding by Loop Extrusion

Chromosome organization poses a remarkable physical problem with many biological consequences: How can molecular interactions between proteins at the nanometer scale organize micron-long chromatinized DNA molecules, insulating or facilitating interactions between specific genomic elements? The mechanism of active loop extrusion holds great promise for explaining interphase and mitotic chromosome folding, yet remains difficult to assay directly. We discuss predictions from our polymer models of loop extrusion with barrier elements and review recent experimental studies that provide strong support for loop extrusion, focusing on perturbations to CTCF and cohesin assayed via Hi-C in interphase. Finally, we discuss a likely molecular mechanism of loop extrusion by structural maintenance of chromosomes complexes.

CH···O Hydrogen Bonds Mediate Highly Specific Recognition of Methylated CpG Sites by the Zinc Finger Protein Kaiso

Many eukaryotic transcription factors recognize the epigenetic marker 5-methylcytosine (mC) at CpG sites in DNA. Despite their structural diversity, methyl-CpG-binding proteins (MBPs) share a common mode of recognition of mC methyl groups that involves hydrophobic pockets and weak hydrogen bonds of the CH···O type. The zinc finger protein Kaiso possesses a remarkably high specificity for methylated over unmethylated CpG sites. A key contribution to this specificity is provided by glutamate 535 (E535), which is optimally positioned to form multiple interactions with mCpG, including direct CH···O hydrogen bonds. To examine the role of E535 and CH···O hydrogen bonding in the preferential recognition of mCpG sites, we determined the structures of wild type Kaiso (WT) and E535 mutants and characterized their interactions with methylated DNA by nuclear magnetic resonance spectroscopy (NMR), X-ray crystallography, and in vitro protein-DNA binding assays. Our data show that Kaiso favors an mCpG over a CpG site by 2 orders of magnitude in affinity and that an important component of this effect is the presence of hydrophobic and CH···O contacts involving E535. Moreover, we present the first direct evidence for formation of a CH···O hydrogen bond between an MBP and 5-methylcytosine by using experimental (NMR) and quantum mechanical chemical shift analysis of the mC methyl protons. Together, our findings uncover a critical function of methyl-specific interactions, including CH···O hydrogen bonds, that optimize the specificity and affinity of MBPs for methylated DNA and contribute to the precise control of gene expression.

Developing in 3D: the role of CTCF in cell differentiation

CTCF is a highly conserved zinc-finger DNA-binding protein that mediates interactions between distant sequences in the genome. As a consequence, CTCF regulates enhancer-promoter interactions and contributes to the three-dimensional organization of the genome. Recent studies indicate that CTCF is developmentally regulated, suggesting that it plays a role in cell type-specific genome organization. Here, we review these studies and discuss how CTCF functions during the development of various cell and tissue types, ranging from embryonic stem cells and gametes, to neural, muscle and cardiac cells. We propose that the lineage-specific control of CTCF levels, and its partnership with lineage-specific transcription factors, allows for the control of cell type-specific gene expression via chromatin looping.

DNA Conformation Induces Adaptable Binding by Tandem Zinc Finger Proteins

Tandem zinc finger (ZF) proteins are the largest and most rapidly diverging family of DNA-binding transcription regulators in mammals. ZFP568 represses a transcript of placental-specific insulin like growth factor 2 (Igf2-P0) in mice. ZFP568 binds a 24-base pair sequence-specific element upstream of Igf2-P0 via the eleven-ZF array. Both DNA and protein conformations deviate from the conventional one finger-three bases recognition, with individual ZFs contacting 2, 3, or 4 bases and recognizing thymine on the opposite strand. These interactions arise from a shortened minor groove caused by an AT-rich stretch, suggesting adaptability of ZF arrays to sequence variations. Despite conservation in mammals, mutations at Igf2 and ZFP568 reduce their binding affinity in chimpanzee and humans. Our studies provide important insights into the evolutionary and structural dynamics of ZF-DNA interactions that play a key role in mammalian development and evolution.

CTCF-Induced Circular DNA Complexes Observed by Atomic Force Microscopy

The CTCF protein has emerged as a key architectural protein involved in genome organization. Although hypothesized to initiate DNA looping, direct evidence of CTCF-induced DNA loop formation is still missing. Several studies have shown that the 11 zinc finger (11 ZF) domain of CTCF is actively involved in DNA binding. We here use atomic force microscopy to examine the effect of the 11 ZF domain comprising residues 266-579 (11 ZF CTCF) and the 3 ZF domain comprising residues 402-494 (6-8 ZF CTCF) of human CTCF on the DNA morphology. Our results show that both domains alter the DNA architecture from the relaxed morphology observed in control DNA samples to compact circular complexes, meshes, and networks, offering important insights into the multivalent character of the 11 ZF CTCF domain. Atomic force microscopy images reveal quasi-circular DNA/CTCF complexes, which are destabilized upon replacing the 11 ZF CTCF by the 6-8 ZF CTCF domain, highlighting the role of the 11 ZF motif in loop formation. Intriguingly, the formation of circular DNA/CTCF complexes is dominated by non-specific binding, whereby contour length and height profiles suggest a single DNA molecule twice wrapped around the protein.

Readers of DNA methylation, the MBD family as potential therapeutic targets

DNA methylation represents a fundamental epigenetic modification that regulates chromatin architecture and gene transcription. Many diseases, including cancer, show aberrant methylation patterns that contribute to the disease phenotype. DNA methylation inhibitors have been used to block methylation dependent gene silencing to treat hematopoietic neoplasms and to restore expression of developmentally silenced genes. However, these inhibitors disrupt methylation globally and show significant off-target toxicities. As an alternative approach, we have been studying readers of DNA methylation, the 5-methylcytosine binding domain family of proteins, as potential therapeutic targets to restore expression of aberrantly and developmentally methylated and silenced genes. In this review, we discuss the role of DNA methylation in gene regulation and cancer development, the structure and function of the 5-methylcytosine binding domain family of proteins, and the possibility of targeting the complexes these proteins form to treat human disease.

Molecular mechanism of directional CTCF recognition of a diverse range of genomic sites

CTCF, a conserved 3D genome architecture protein, determines proper genome-wide chromatin looping interactions through directional binding to specific sequence elements of four modules within numerous CTCF-binding sites (CBSs) by its 11 zinc fingers (ZFs). Here, we report four crystal structures of human CTCF in complex with CBSs of the protocadherin (Pcdh) clusters. We show that directional CTCF binding to cognate CBSs of the Pcdh enhancers and promoters is achieved through inserting its ZF3, ZFs 4-7, and ZFs 9-11 into the major groove along CBSs, resulting in a sequence-specific recognition of module 4, modules 3 and 2, and module 1, respectively; and ZF8 serves as a spacer element for variable distances between modules 1 and 2. In addition, the base contact with the asymmetric "A" in the central position of modules 2-3, is essential for directional recognition of the CBSs with symmetric core sequences but lacking module 1. Furthermore, CTCF tolerates base changes at specific positions within the degenerated CBS sequences, permitting genome-wide CTCF binding to a diverse range of CBSs. Together, these complex structures provide important insights into the molecular mechanisms for the directionality, diversity, flexibility, dynamics, and conservation of multivalent CTCF binding to its cognate sites across the entire human genome.

Structural basis of human PR/SET domain 9 (PRDM9) allele C-specific recognition of its cognate DNA sequence

PRDM9 is the only mammalian gene that has been associated with speciation. The PR/SET domain 9 (PRDM9) protein is a major determinant of meiotic recombination hot spots and acts through sequence-specific DNA binding via its C2H2 zinc finger (ZF) tandem array, which is highly polymorphic within and between species. The most common human variant, PRDM9 allele A (PRDM9a), contains 13 fingers (ZF1-13). Allele C (PRDM9c) is the second-most common among African populations and differs from PRDM9a by an arginine-to-serine change (R764S) in ZF9 and by replacement of ZF11 with two other fingers, yielding 14 fingers in PRDM9c. Here we co-crystallized the six-finger fragment ZF8-13 of PRDM9c, in complex with an oligonucleotide representing a known PRDM9c-specific hot spot sequence, and compared the structure with that of a characterized PRDM9a-specific complex. There are three major differences. First, Ser⁷⁶⁴ in ZF9 allows PRDM9c to accommodate a variable base, whereas PRDM9a Arg⁷⁶⁴ recognizes a conserved guanine. Second, the two-finger expansion of ZF11 allows PRDM9c to recognize three-base-pair-longer sequences. A tryptophan in the additional ZF interacts with a conserved thymine methyl group. Third, an Arg-Asp dipeptide immediately preceding the ZF helix, conserved in two PRDM9a fingers and three PRDM9c fingers, permits adaptability to variations from a C:G base pair (G-Arg interaction) to a G:C base pair (C-Asp interaction). This Arg-Asp conformational switch allows identical ZF modules to recognize different sequences. Our findings illuminate the molecular mechanisms for flexible and conserved binding of human PRDM9 alleles to their cognate DNA sequences.