Multi-feature clustering of CTCF binding creates robustness for loop extrusion blocking and Topologically Associating Domain boundaries

Systematic screening of enhancer-blocking insulators in Drosophila identifies their DNA sequence determinants

Long-range transcriptional activation of gene promoters by abundant enhancers in animal genomes calls for mechanisms to limit inappropriate regulation. DNA elements called insulators serve this purpose by shielding promoters from an enhancer when interposed. Unlike promoters and enhancers, insulators have not been systematically characterized due to lacking high-throughput screening assays, and questions regarding how insulators are distributed and encoded in the genome remain. Here, we establish "insulator-seq" as a plasmid-based massively parallel reporter assay in Drosophila cultured cells to perform a systematic insulator screen of selected genomic loci. Screening developmental gene loci showed that not all insulator protein binding sites effectively block enhancer-promoter communication. Deep insulator mutagenesis identified sequences flexibly positioned around the CTCF insulator protein binding motif that are critical for functionality. The ability to screen millions of DNA sequences without positional effect has enabled functional mapping of insulators and provided further insights into the determinants of insulators.

The impact of the embryonic DNA methylation program on CTCF-mediated genome regulation

During mammalian embryogenesis, both the 5-cytosine DNA methylation (5meC) landscape and three dimensional (3D) chromatin architecture are profoundly remodeled during a process known as 'epigenetic reprogramming.' An understudied aspect of epigenetic reprogramming is how the 5meC flux, per se, affects the 3D genome. This is pertinent given the 5meC-sensitivity of DNA binding for a key regulator of chromosome folding: CTCF. We profiled the CTCF binding landscape using a mouse embryonic stem cell (ESC) differentiation protocol that models embryonic 5meC dynamics. Mouse ESCs lacking DNA methylation machinery are able to exit naive pluripotency, thus allowing for dissection of subtle effects of CTCF on gene expression. We performed CTCF HiChIP in both wild-type and mutant conditions to assess gained CTCF-CTCF contacts in the absence of 5meC. We performed H3K27ac HiChIP to determine the impact that ectopic CTCF binding has on cis-regulatory contacts. Using 5meC epigenome editing, we demonstrated that the methyl-mark is able to impair CTCF binding at select loci. Finally, a detailed dissection of the imprinted Zdbf2 locus showed how 5meC-antagonism of CTCF allows for proper gene regulation during differentiation. This work provides a comprehensive overview of how 5meC impacts the 3D genome in a relevant model for early embryonic events.

Permeable TAD boundaries and their impact on genome-associated functions

TAD boundaries are genomic elements that separate biological processes in neighboring domains by blocking DNA loops that are formed through Cohesin-mediated loop extrusion. Most TAD boundaries consist of arrays of binding sites for the CTCF protein, whose interaction with the Cohesin complex blocks loop extrusion. TAD boundaries are not fully impermeable though and allow a limited amount of inter-TAD loop formation. Based on the reanalysis of Nano-C data, a multicontact Chromosome Conformation Capture assay, we propose a model whereby clustered CTCF binding sites promote the successive stalling of Cohesin and subsequent dissociation from the chromatin. A fraction of Cohesin nonetheless achieves boundary read-through. Due to a constant rate of Cohesin dissociation elsewhere in the genome, the maximum length of inter-TAD loops is restricted though. We speculate that the DNA-encoded organization of stalling sites regulates TAD boundary permeability and discuss implications for enhancer-promoter loop formation and other genomic processes.

Prediction and functional interpretation of inter-chromosomal genome architecture from DNA sequence with TwinC

Three-dimensional nuclear DNA architecture comprises well-studied intra-chromosomal (cis) folding and less characterized inter-chromosomal (trans) interfaces. Current predictive models of 3D genome folding can effectively infer pairwise cis-chromatin interactions from the primary DNA sequence but generally ignore trans contacts. There is an unmet need for robust models of trans-genome organization that provide insights into their underlying principles and functional relevance. We present TwinC, an interpretable convolutional neural network model that reliably predicts trans contacts measurable through genome-wide chromatin conformation capture (Hi-C). TwinC uses a paired sequence design from replicate Hi-C experiments to learn single base pair relevance in trans interactions across two stretches of DNA. The method achieves high predictive accuracy (AUROC=0.80) on a cross-chromosomal test set from Hi-C experiments in heart tissue. Mechanistically, the neural network learns the importance of compartments, chromatin accessibility, clustered transcription factor binding and G-quadruplexes in forming trans contacts. In summary, TwinC models and interprets trans genome architecture, shedding light on this poorly understood aspect of gene regulation.

The characteristics of CTCF binding sequences contribute to enhancer blocking activity

While the elements encoding enhancers and promoters have been relatively well studied, the full spectrum of insulator elements which bind the CCCTC binding factor (CTCF), is relatively poorly characterized. This is partly due to the genomic context of CTCF sites greatly influencing their roles and activity. Here we have developed an experimental system to determine the ability of minimal, consistently sized, individual CTCF elements to interpose between enhancers and promoters and thereby reduce gene expression during differentiation. Importantly, each element is tested in the identical location thereby minimising the effect of genomic context. We found no correlation between the ability of CTCF elements to block enhancer-promoter activity with the degree of evolutionary conservation; their resemblance to the consensus core sequences; or the number of CTCF core motifs harboured in the element. Nevertheless, we have shown that the strongest enhancer-promoter blockers include a previously described bound element lying upstream of the CTCF core motif. In addition, we found other uncharacterised DNaseI footprints located close to the core motif that may affect function. We have developed an assay of CTCF sequences which will enable researchers to sub-classify individual CTCF elements in a uniform and unbiased way.

Cooperative insulation of regulatory domains by CTCF-dependent physical insulation and promoter competition

The specificity of gene expression during development requires the insulation of regulatory domains to avoid inappropriate enhancer-gene interactions. In vertebrates, this insulator function is mostly attributed to clusters of CTCF sites located at topologically associating domain (TAD) boundaries. However, TAD boundaries allow some physical crosstalk across regulatory domains, which is at odds with the specific and precise expression of developmental genes. Here we show that developmental genes and nearby clusters of CTCF sites cooperatively foster the robust insulation of regulatory domains. By genetically dissecting a couple of representative loci in mouse embryonic stem cells, we show that CTCF sites prevent undesirable enhancer-gene contacts (i.e. physical insulation), while developmental genes preferentially contribute to regulatory insulation through non-structural mechanisms involving promoter competition rather than enhancer blocking. Overall, our work provides important insights into the insulation of regulatory domains, which in turn might help interpreting the pathological consequences of certain structural variants.

Structural perturbation of chromatin domains with multiple developmental regulators can severely impact gene regulation and development

Chromatin domain boundaries delimited by CTCF motifs can restrict the range of enhancer action. However, disruption of domain structure often results in mild gene dysregulation and thus predicting the impact of boundary rearrangements on animal development remains challenging. Here, we tested whether structural perturbation of a chromatin domain with multiple developmental regulators can result in more acute gene dysregulation and severe developmental phenotypes. We targeted clusters of CTCF motifs in a domain of the mouse genome containing three FGF ligand genes-Fgf3, Fgf4, and Fgf15-that regulate several developmental processes. Deletion of the 23.9kb cluster that defines the centromeric boundary of this domain resulted in ectopic interactions of the FGF genes with enhancers located across the deleted boundary that are active in the developing brain. This caused strong induction of FGF expression and perinatal lethality with encephalocele and orofacial cleft phenotypes. Heterozygous boundary deletion was sufficient to cause these fully penetrant phenotypes, and strikingly, loss of a single CTCF motif within the cluster also recapitulated ectopic FGF expression and caused encephalocele. However, such phenotypic sensitivity to perturbation of domain structure did not extend to all CTCF clusters of this domain, nor to all developmental processes controlled by these three FGF genes-for example, the ability to undergo lineage specification in the blastocyst and pre-implantation development were not affected. By tracing the impact of different chromosomal rearrangements throughout mouse development, we start to uncover the determinants of phenotypic robustness and sensitivity to perturbation of chromatin boundaries. Our data show how small sequence variants at certain domain boundaries can have a surprisingly outsized effect and must be considered as potential sources of gene dysregulation during development and disease.

Single-molecule imaging reveals a direct role of CTCF's zinc fingers in SA interaction and cluster-dependent RNA recruitment

CTCF is a zinc finger protein associated with transcription regulation that also acts as a barrier factor for topologically associated domains (TADs) generated by cohesin via loop extrusion. These processes require different properties of CTCF-DNA interaction, and it is still unclear how CTCF's structural features may modulate its diverse roles. Here, we employ single-molecule imaging to study both full-length CTCF and truncation mutants. We show that CTCF enriches at CTCF binding sites (CBSs), displaying a longer lifetime than observed previously. We demonstrate that the zinc finger domains mediate CTCF clustering and that clustering enables RNA recruitment, possibly creating a scaffold for interaction with RNA-binding proteins like cohesin's subunit SA. We further reveal a direct recruitment and an increase of SA residence time by CTCF bound at CBSs, suggesting that CTCF-SA interactions are crucial for cohesin stability on chromatin at TAD borders. Furthermore, we establish a single-molecule T7 transcription assay and show that although a transcribing polymerase can remove CTCF from CBSs, transcription is impaired. Our study shows that context-dependent nucleic acid binding determines the multifaceted CTCF roles in genome organization and transcription regulation.

Pairtools: From sequencing data to chromosome contacts

The field of 3D genome organization produces large amounts of sequencing data from Hi-C and a rapidly-expanding set of other chromosome conformation protocols (3C+). Massive and heterogeneous 3C+ data require high-performance and flexible processing of sequenced reads into contact pairs. To meet these challenges, we present pairtools-a flexible suite of tools for contact extraction from sequencing data. Pairtools provides modular command-line interface (CLI) tools that can be flexibly chained into data processing pipelines. The core operations provided by pairtools are parsing of.sam alignments into Hi-C pairs, sorting and removal of PCR duplicates. In addition, pairtools provides auxiliary tools for building feature-rich 3C+ pipelines, including contact pair manipulation, filtration, and quality control. Benchmarking pairtools against popular 3C+ data pipelines shows advantages of pairtools for high-performance and flexible 3C+ analysis. Finally, pairtools provides protocol-specific tools for restriction-based protocols, haplotype-resolved contacts, and single-cell Hi-C. The combination of CLI tools and tight integration with Python data analysis libraries makes pairtools a versatile foundation for a broad range of 3C+ pipelines.

A CTCF-dependent mechanism underlies the Hox timer: relation to a segmented body plan

During gastrulation, Hox genes are activated in a time-sequence that follows the order of the genes along their clusters. This property, which is observed in all animals that develop following a progressive rostral-to-caudal morphogenesis, is associated with changes in the chromatin structure and epigenetic profiles of Hox clusters, suggesting a process at least partly based on sequential gene accessibility. Here, we discuss recent work on this issue, as well as a possible mechanism based on the surprising conservation in both the distribution and orientation of CTCF sites inside vertebrate Hox clusters.

Identification of a distal enhancer regulating hedgehog interacting protein gene in human lung epithelial cells

BACKGROUND

An intergenic region at chromosome 4q31 is one of the most significant regions associated with COPD susceptibility and lung function in GWAS. In this region, the implicated causal gene HHIP has a unique epithelial expression pattern in adult human lungs, in contrast to dominant expression in fibroblasts in murine lungs. However, the mechanism underlying the species-dependent cell type-specific regulation of HHIP remains largely unknown.

METHODS

We employed snATAC-seq analysis to identify open chromatin regions within the COPD GWAS region in various human lung cell types. ChIP-quantitative PCR, reporter assays, chromatin conformation capture assays and Hi-C assays were conducted to characterize the regulatory element in this region. CRISPR/Cas9-editing was performed in BEAS-2B cells to generate single colonies with stable knockout of the regulatory element. RT-PCR and Western blot assays were used to evaluate expression of HHIP and epithelial-mesenchymal transition (EMT)-related marker genes.

FINDINGS

We identified a distal enhancer within the COPD 4q31 GWAS locus that regulates HHIP transcription at baseline and after TGFβ treatment in a SMAD3-dependent, but Hedgehog-independent manner in human bronchial epithelial cells. The distal enhancer also maintains chromatin topological domains near 4q31 locus and HHIP gene. Reduced HHIP expression led to increased EMT induced by TGFβ in human bronchial epithelial cells.

INTERPRETATION

A distal enhancer regulates HHIP expression both under homeostatic condition and upon TGFβ treatment in human bronchial epithelial cells. The interaction between HHIP and TGFβ signalling possibly contributes to COPD pathogenesis.

FUNDING

Supported by NIH grants R01HL127200, R01HL148667 and R01HL162783 (to X. Z).

The impact of DNA methylation on CTCF-mediated 3D genome organization

Cytosine DNA methylation is a highly conserved epigenetic mark in eukaryotes. Although the role of DNA methylation at gene promoters and repetitive elements has been extensively studied, the function of DNA methylation in other genomic contexts remains less clear. In the nucleus of mammalian cells, the genome is spatially organized at different levels, and strongly influences myriad genomic processes. There are a number of factors that regulate the three-dimensional (3D) organization of the genome, with the CTCF insulator protein being among the most well-characterized. Pertinently, CTCF binding has been reported as being DNA methylation-sensitive in certain contexts, perhaps most notably in the process of genomic imprinting. Therefore, it stands to reason that DNA methylation may play a broader role in the regulation of chromatin architecture. Here we summarize the current understanding that is relevant to both the mammalian DNA methylation and chromatin architecture fields and attempt to assess the extent to which DNA methylation impacts the folding of the genome. The focus is in early embryonic development and cellular transitions when the epigenome is in flux, but we also describe insights from pathological contexts, such as cancer, in which the epigenome and 3D genome organization are misregulated.

Outward-oriented sites within clustered CTCF boundaries are key for intra-TAD chromatin interactions and gene regulation

CTCF plays an important role in 3D genome organization by adjusting the strength of chromatin insulation at TAD boundaries, where clustered CBS (CTCF-binding site) elements are often arranged in a tandem array with a complex divergent or convergent orientation. Here, using Pcdh and HOXD loci as a paradigm, we look into the clustered CTCF TAD boundaries and find that, counterintuitively, outward-oriented CBS elements are crucial for inward enhancer-promoter interactions as well as for gene regulation. Specifically, by combinatorial deletions of a series of putative enhancer elements in mice in vivo or CBS elements in cultured cells in vitro, in conjunction with chromosome conformation capture and RNA-seq analyses, we show that deletions of outward-oriented CBS elements weaken the strength of long-distance intra-TAD promoter-enhancer interactions and enhancer activation of target genes. Our data highlight the crucial role of outward-oriented CBS elements within the clustered CTCF TAD boundaries in developmental gene regulation and have interesting implications on the organization principles of clustered CTCF sites within TAD boundaries.