A genome-wide map of CTCF multivalency redefines the CTCF code

Nickel-induced alterations to chromatin structure and function

Nickel (Ni), a heavy metal is prevalent in the atmosphere due to both natural and anthropogenic activities. Ni is a carcinogen implicated in the development of lung and nasal cancers in humans. Furthermore, Ni exposure is associated with a number of chronic lung diseases in humans including asthma, chronic bronchitis, emphysema, pulmonary fibrosis, pulmonary edema and chronic obstructive pulmonary disease (COPD). While Ni compounds are weak mutagens, a number of studies have demonstrated the potential of Ni to alter the epigenome, suggesting epigenomic dysregulation as an important underlying cause for its pathogenicity. In the eukaryotic nucleus, the DNA is organized in a three-dimensional (3D) space through assembly of higher order chromatin structures. Such an organization is critically important for transcription and other biological activities. Accumulating evidence suggests that by negatively affecting various cellular regulatory processes, Ni could potentially affect chromatin organization. In this review, we discuss the role of Ni in altering the chromatin architecture, which potentially plays a major role in Ni pathogenicity.

"Stripe" transcription factors provide accessibility to co-binding partners in mammalian genomes

Regulatory elements activate promoters by recruiting transcription factors (TFs) to specific motifs. Notably, TF-DNA interactions often depend on cooperativity with colocalized partners, suggesting an underlying cis-regulatory syntax. To explore TF cooperativity in mammals, we analyze ∼500 mouse and human primary cells by combining an atlas of TF motifs, footprints, ChIP-seq, transcriptomes, and accessibility. We uncover two TF groups that colocalize with most expressed factors, forming stripes in hierarchical clustering maps. The first group includes lineage-determining factors that occupy DNA elements broadly, consistent with their key role in tissue-specific transcription. The second one, dubbed universal stripe factors (USFs), comprises ∼30 SP, KLF, EGR, and ZBTB family members that recognize overlapping GC-rich sequences in all tissues analyzed. Knockouts and single-molecule tracking reveal that USFs impart accessibility to colocalized partners and increase their residence time. Mammalian cells have thus evolved a TF superfamily with overlapping DNA binding that facilitate chromatin accessibility.

CTCF shapes chromatin structure and gene expression in health and disease

CCCTC-binding factor (CTCF) is an eleven zinc finger (ZF), multivalent transcriptional regulator, that recognizes numerous motifs thanks to the deployment of distinct combinations of its ZFs. The great majority of the ~50,000 genomic locations bound by the CTCF protein in a given cell type is intergenic, and a fraction of these sites overlaps with transcriptional enhancers. Furthermore, a proportion of the regions bound by CTCF intersect genes and promoters. This suggests multiple ways in which CTCF may impact gene expression. At promoters, CTCF can directly affect transcription. At more distal sites, CTCF may orchestrate interactions between regulatory elements and help separate eu- and heterochromatic areas in the genome, exerting a chromatin barrier function. In this review, we outline how CTCF contributes to the regulation of the three-dimensional structure of chromatin and the formation of chromatin domains. We discuss how CTCF binding and architectural functions are regulated. We examine the literature implicating CTCF in controlling gene expression in development and disease both by acting as an insulator and a factor facilitating regulatory elements to efficiently interact with each other in the nuclear space.

The role of chromatin loop extrusion in antibody diversification

Cohesin mediates chromatin loop formation across the genome by extruding chromatin between convergently oriented CTCF-binding elements. Recent studies indicate that cohesin-mediated loop extrusion in developing B cells presents immunoglobulin heavy chain (Igh) variable (V), diversity (D) and joining (J) gene segments to RAG endonuclease through a process referred to as RAG chromatin scanning. RAG initiates V(D)J recombinational joining of these gene segments to generate the large number of different Igh variable region exons that are required for immune responses to diverse pathogens. Antigen-activated mature B cells also use chromatin loop extrusion to mediate the synapsis, breakage and end joining of switch regions flanking Igh constant region exons during class-switch recombination, which allows for the expression of different antibody constant region isotypes that optimize the functions of antigen-specific antibodies to eliminate pathogens. Here, we review recent advances in our understanding of chromatin loop extrusion during V(D)J recombination and class-switch recombination at the Igh locus.

Single base-pair resolution analysis of DNA binding motif with MoMotif reveals an oncogenic function of CTCF zinc-finger 1 mutation

Defining the impact of missense mutations on the recognition of DNA motifs is highly dependent on bioinformatic tools that define DNA binding elements. However, classical motif analysis tools remain limited in their capacity to identify subtle changes in complex binding motifs between distinct conditions. To overcome this limitation, we developed a new tool, MoMotif, that facilitates a sensitive identification, at the single base-pair resolution, of complex, or subtle, alterations to core binding motifs, discerned from ChIP-seq data. We employed MoMotif to define the previously uncharacterized recognition motif of CTCF zinc-finger 1 (ZF1), and to further define the impact of CTCF ZF1 mutation on its association with chromatin. Mutations of CTCF ZF1 are exclusive to breast cancer and are associated with metastasis and therapeutic resistance, but the underlying mechanisms are unclear. Using MoMotif, we identified an extension of the CTCF core binding motif, necessitating a functional ZF1 to bind appropriately. Using a combination of ChIP-Seq and RNA-Seq, we discover that the inability to bind this extended motif drives an altered transcriptional program associated with the oncogenic phenotypes observed clinically. Our study demonstrates that MoMotif is a powerful new tool for comparative ChIP-seq analysis and characterising DNA-protein contacts.

CTCF DNA binding domain undergoes dynamic and selective protein–protein interactions

CTCF is a predominant insulator protein required for three-dimensional chromatin organization. However, the roles of its insulation of enhancers in a 3D nuclear organization have not been fully explained. Here, we found that the CTCF DNA-binding domain (DBD) forms dynamic self-interacting clusters. Strikingly, CTCF DBD clusters were found to incorporate other insulator proteins but are not coenriched with transcriptional activators in the nucleus. This property is not observed in other domains of CTCF or the DBDs of other transcription factors. Moreover, endogenous CTCF shows a phenotype consistent with the DBD by forming small protein clusters and interacting with CTCF motif arrays that have fewer transcriptional activators bound. Our results reveal an interesting phenomenon in which CTCF DBD interacts with insulator proteins and selectively localizes to nuclear positions with lower concentrations of transcriptional activators, providing insights into the insulation function of CTCF.

•

The CTCF DNA-binding domain forms protein clusters in vivo and in vitro

•

CTCF DBD clusters colocalize with insulator proteins but not with activators

•

Arginine residues of CTCF DBD are frequently mutated in cancers

•

Multiple transcription factor DBDs form protein clusters

NuRD complex recruitment to Thpok mediates CD4+ T cell lineage differentiation

Although BTB-zinc finger (BTB-ZF) transcription factors control the differentiation of multiple hematopoietic and immune lineages, how they function is poorly understood. The BTB-ZF factor Thpok controls intrathymic CD4⁺ T cell development and the expression of most CD4⁺ and CD8⁺ lineage genes. Here, we identify the nucleosome remodeling and deacetylase (NuRD) complex as a critical Thpok cofactor. Using mass spectrometry and coimmunoprecipitation in primary T cells, we show that Thpok binds NuRD components independently of DNA association. We locate three amino acid residues within the Thpok BTB domain that are required for both NuRD binding and Thpok functions. Conversely, a chimeric protein merging the NuRD component Mta2 to a BTB-less version of Thpok supports CD4⁺ T cell development, indicating that NuRD recruitment recapitulates the functions of the Thpok BTB domain. We found that NuRD mediates Thpok repression of CD8⁺ lineage genes, including the transcription factor Runx3, but is dispensable for Cd4 expression. We show that these functions cannot be performed by the BTB domain of the Thpok-related factor Bcl6, which fails to bind NuRD. Thus, cofactor binding critically contributes to the functional specificity of BTB-ZF factors, which control the differentiation of most hematopoietic subsets.

Mitochondrial base editor induces substantial nuclear off-target mutations

DddA-derived cytosine base editors (DdCBEs), which are fusions of the split-DddA halves and transcription activator-like effector (TALE) array proteins, enable targeted C·G-to- T·A conversions in mitochondrial DNA¹. However, its genome-wide specificity is poorly understood. Here we show that the mitochondrial base editor induces extensive off-target editing in the nuclear genome. Genome-wide, unbiased analysis of its editome reveals hundreds of off-target sites that are TALE array sequence (TAS)-dependent or -independent. TAS-dependent off-target sites in the nuclear DNA (nDNA) are often specified by only one of the two TALE repeats, challenging the principle that DdCBEs are guided by a paired TALE proteins positioned in close proximity. TAS-independent nDNA off-target sites are frequently shared among DdCBEs with distinct TALE arrays. Notably, they co-localize strongly with CTCF-binding sites and are enriched in TAD boundaries. We also engineered DdCBE to alleviate such off-target effect. Collectively, our results have implications for the use of DdCBEs in basic research and therapeutic applications, and suggest the need to thoroughly define and evaluate the off-target effects of base editing tools.

CTCF: a misguided jack-of-all-trades in cancer cells

The emergence and progression of cancers is accompanied by a dysregulation of transcriptional programs. The three-dimensional (3D) organization of the human genome has emerged as an important multi-level mediator of gene transcription and regulation. In cancer cells, this organization can be restructured, providing a framework for the deregulation of gene activity. The CTCF protein, initially identified as the product from a tumor suppressor gene, is a jack-of-all-trades for the formation of 3D genome organization in normal cells. Here, we summarize how CTCF is involved in the multi-level organization of the human genome and we discuss emerging insights into how perturbed CTCF function and DNA binding causes the activation of oncogenes in cancer cells, mostly through a process of enhancer hijacking. Moreover, we highlight non-canonical functions of CTCF that can be relevant for the emergence of cancers as well. Finally, we provide guidelines for the computational identification of perturbed CTCF binding and reorganized 3D genome structure in cancer cells.

Can aggressive cancers be identified by the "aggressiveness" of their chromatin?

Phenotypic plasticity is a crucial feature of aggressive cancer, providing the means for cancer progression. Stochastic changes in tumor cell transcriptional programs increase the chances of survival under any condition. I hypothesize that unstable chromatin permits stochastic transitions between transcriptional programs in aggressive cancers and supports non-genetic heterogeneity of tumor cells as a basis for their adaptability. I present a mechanistic model for unstable chromatin which includes destabilized nucleosomes, mobile chromatin fibers and random enhancer-promoter contacts, resulting in stochastic transcription. I suggest potential markers for "unsettled" chromatin in tumors associated with poor prognosis. Although many of the characteristics of unstable chromatin have been described, they were mostly used to explain changes in the transcription of individual genes. I discuss approaches to evaluate the role of unstable chromatin in non-genetic tumor cell heterogeneity and suggest using the degree of chromatin instability and transcriptional noise in tumor cells to predict cancer aggressiveness.

Neural network modeling of differential binding between wild-type and mutant CTCF reveals putative binding preferences for zinc fingers 1-2

BACKGROUND

Many transcription factors (TFs), such as multi zinc-finger (ZF) TFs, have multiple DNA binding domains (DBDs), and deciphering the DNA binding motifs of individual DBDs is a major challenge. One example of such a TF is CCCTC-binding factor (CTCF), a TF with eleven ZFs that plays a variety of roles in transcriptional regulation, most notably anchoring DNA loops. Previous studies found that CTCF ZFs 3-7 bind CTCF's core motif and ZFs 9-11 bind a specific upstream motif, but the motifs of ZFs 1-2 have yet to be identified.

RESULTS

We developed a new approach to identifying the binding motifs of individual DBDs of a TF through analyzing chromatin immunoprecipitation sequencing (ChIP-seq) experiments in which a single DBD is mutated: we train a deep convolutional neural network to predict whether wild-type TF binding sites are preserved in the mutant TF dataset and interpret the model. We applied this approach to mouse CTCF ChIP-seq data and identified the known binding preferences of CTCF ZFs 3-11 as well as a putative GAG binding motif for ZF 1. We analyzed other CTCF datasets to provide additional evidence that ZF 1 is associated with binding at the motif we identified, and we found that the presence of the motif for ZF 1 is associated with CTCF ChIP-seq peak strength.

CONCLUSIONS

Our approach can be applied to any TF for which in vivo binding data from both the wild-type and mutated versions of the TF are available, and our findings provide new potential insights binding preferences of CTCF's DBDs.

Rice RS2-9, which is bound by transcription factor OSH1, blocks enhancer-promoter interactions in plants

Insulators characterized in Drosophila and mammals have been shown to play a key role in the restriction of promiscuous enhancer-promoter interactions, as well as reshaping the topological landscape of chromosomes. Yet the role of insulators in plants remains poorly understood, in large part because of a lack of well-characterized insulators and binding factor(s). In this study, we isolated a 1.2-kb RS2-9 insulator from the Oryza sativa (rice) genome that can, when interposed between an enhancer and promoter, efficiently block the activation function of both constitutive and floral organ-specific enhancers in transgenic Arabidopsis and Nicotiana tabacum (tobacco). In the rice genome, the genes flanking RS2-9 exhibit an absence of mutual transcriptional interactions, as well as a lack of histone modification spread. We further determined that O. sativa Homeobox 1 (OSH1) bound two regions of RS2-9, as well as over 50 000 additional sites in the rice genome, the majority of which resided in intergenic regions. Mutation of one of the two OSH1-binding sites in RS2-9 impaired insulation activity by up to 60%, whereas the mutation of both binding sites virtually abolished insulator function. We also demonstrated that OSH1 binding sites were associated with 72% of the boundaries of topologically associated domains (TADs) identified in the rice genome, which is comparable to the 77% of TAD boundaries bound by the insulator CCCTC-binding factor (CTCF) in mammals. Taken together, our findings indicate that OSH1-RS2-9 acts as a true insulator in plants, and highlight a potential role for OSH1 in gene insulation and topological organization in plant genomes.

Structures and biological functions of zinc finger proteins and their roles in hepatocellular carcinoma

Zinc finger proteins are transcription factors with the finger domain, which plays a significant role in gene regulation. As the largest family of transcription factors in the human genome, zinc finger (ZNF) proteins are characterized by their different DNA binding motifs, such as C2H2 and Gag knuckle. Different kinds of zinc finger motifs exhibit a wide variety of biological functions. Zinc finger proteins have been reported in various diseases, especially in several cancers. Hepatocellular carcinoma (HCC) is the third leading cause of cancer-associated death worldwide, especially in China. Most of HCC patients have suffered from hepatitis B virus (HBV) and hepatitis C virus (HCV) injection for a long time. Although the surgical operation of HCC has been extremely developed, the prognosis of HCC is still very poor, and the underlying mechanisms in HCC tumorigenesis are still not completely understood. Here, we summarize multiple functions and recent research of zinc finger proteins in HCC tumorigenesis and progression. We also discuss the significance of zinc finger proteins in HCC diagnosis and prognostic evaluation.

Epigenetic Regulation of Human T-Cell Leukemia Virus Gene Expression

Viral and cellular gene expression are regulated by epigenetic alterations, including DNA methylation, histone modifications, nucleosome positioning, and chromatin looping. Human T-cell leukemia virus type 1 (HTLV-1) is a pathogenic retrovirus associated with inflammatory disorders and T-cell lymphoproliferative malignancy. The transforming activity of HTLV-1 is driven by the viral oncoprotein Tax, which acts as a transcriptional activator of the cAMP response element-binding protein (CREB) and nuclear factor kappa B (NFκB) pathways. The epigenetic effects of Tax and the induction of lymphoproliferative malignancy include alterations in DNA methylation and histone modifications. In addition, alterations in nucleosome positioning and DNA looping also occur in HTLV-1-induced malignant cells. A mechanistic definition of these effects will pave the way to new therapies for HTLV-1-associated disorders.

Jpx RNA regulates CTCF anchor site selection and formation of chromosome loops

Chromosome loops shift dynamically during development, homeostasis, and disease. CCCTC-binding factor (CTCF) is known to anchor loops and construct 3D genomes, but how anchor sites are selected is not yet understood. Here, we unveil Jpx RNA as a determinant of anchor selectivity. Jpx RNA targets thousands of genomic sites, preferentially binding promoters of active genes. Depleting Jpx RNA causes ectopic CTCF binding, massive shifts in chromosome looping, and downregulation of >700 Jpx target genes. Without Jpx, thousands of lost loops are replaced by de novo loops anchored by ectopic CTCF sites. Although Jpx controls CTCF binding on a genome-wide basis, it acts selectively at the subset of developmentally sensitive CTCF sites. Specifically, Jpx targets low-affinity CTCF motifs and displaces CTCF protein through competitive inhibition. We conclude that Jpx acts as a CTCF release factor and shapes the 3D genome by regulating anchor site usage.

Structure-function relationships explain CTCF zinc finger mutation phenotypes in cancer

CCCTC-binding factor (CTCF) plays fundamental roles in transcriptional regulation and chromatin architecture maintenance. CTCF is also a tumour suppressor frequently mutated in cancer, however, the structural and functional impact of mutations have not been examined. We performed molecular and structural characterisation of five cancer-specific CTCF missense zinc finger (ZF) mutations occurring within key intra- and inter-ZF residues. Functional characterisation of CTCF ZF mutations revealed a complete (L309P, R339W, R377H) or intermediate (R339Q) abrogation as well as an enhancement (G420D) of the anti-proliferative effects of CTCF. DNA binding at select sites was disrupted and transcriptional regulatory activities abrogated. Molecular docking and molecular dynamics confirmed that mutations in residues specifically contacting DNA bases or backbone exhibited loss of DNA binding. However, R339Q and G420D were stabilised by the formation of new primary DNA bonds, contributing to gain-of-function. Our data confirm that a spectrum of loss-, change- and gain-of-function impacts on CTCF zinc fingers are observed in cell growth regulation and gene regulatory activities. Hence, diverse cellular phenotypes of mutant CTCF are clearly explained by examining structure-function relationships.

Sequential in cis mutagenesis in vivo reveals various functions for CTCF sites at the mouse HoxD cluster

Mammalian Hox gene clusters contain a range of CTCF binding sites. In addition to their importance in organizing a TAD border, which isolates the most posterior genes from the rest of the cluster, the positions and orientations of these sites suggest that CTCF may be instrumental in the selection of various subsets of contiguous genes, which are targets of distinct remote enhancers located in the flanking regulatory landscapes. We examined this possibility by producing an allelic series of cumulative in cis mutations in these sites, up to the abrogation of CTCF binding in the five sites located on one side of the TAD border. In the most impactful alleles, the global chromatin architecture of the locus was modified, yet not drastically, illustrating that CTCF sites located on one side of a strong TAD border are sufficient to organize at least part of this insulation. Spatial colinearity in the expression of these genes along the major body axis was nevertheless maintained, despite abnormal expression boundaries. In contrast, strong effects were scored in the selection of target genes responding to particular enhancers, leading to the misregulation of Hoxd genes in specific structures. Altogether, while most enhancer-promoter interactions can occur in the absence of this series of CTCF sites, the binding of CTCF in the Hox cluster is required to properly transform a rather unprecise process into a highly discriminative mechanism of interactions, which is translated into various patterns of transcription accompanied by the distinctive chromatin topology found at this locus. Our allelic series also allowed us to reveal the distinct functional contributions for CTCF sites within this Hox cluster, some acting as insulator elements, others being necessary to anchor or stabilize enhancer-promoter interactions, and some doing both, whereas they all together contribute to the formation of a TAD border. This variety of tasks may explain the amazing evolutionary conservation in the distribution of these sites among paralogous Hox clusters or between various vertebrates.

Decoding the Equine Genome: Lessons from ENCODE

The horse reference genome assemblies, EquCab2.0 and EquCab3.0, have enabled great advancements in the equine genomics field, from tools to novel discoveries. However, significant gaps of knowledge regarding genome function remain, hindering the study of complex traits in horses. In an effort to address these gaps and with inspiration from the Encyclopedia of DNA Elements (ENCODE) project, the equine Functional Annotation of Animal Genome (FAANG) initiative was proposed to bridge the gap between genome and gene expression, providing further insights into functional regulation within the horse genome. Three years after launching the initiative, the equine FAANG group has generated data from more than 400 experiments using over 50 tissues, targeting a variety of regulatory features of the equine genome. In this review, we examine how valuable lessons learned from the ENCODE project informed our decisions in the equine FAANG project. We report the current state of the equine FAANG project and discuss how FAANG can serve as a template for future expansion of functional annotation in the equine genome and be used as a reference for studies of complex traits in horse. A well-annotated reference functional atlas will also help advance equine genetics in the pan-genome and precision medicine era.

Long-Range Control of Class Switch Recombination by Transcriptional Regulatory Elements

Immunoglobulin class switch recombination (CSR) plays a crucial role in adaptive immune responses through a change of the effector functions of antibodies and is triggered by T-cell-dependent as well as T-cell-independent antigens. Signals generated following encounter with each type of antigen direct CSR to different isotypes. At the genomic level, CSR occurs between highly repetitive switch sequences located upstream of the constant gene exons of the immunoglobulin heavy chain locus. Transcription of switch sequences is mandatory for CSR and is induced in a stimulation-dependent manner. Switch transcription takes place within dynamic chromatin domains and is regulated by long-range regulatory elements which promote alignment of partner switch regions in CSR centers. Here, we review recent work and models that account for the function of long-range transcriptional regulatory elements and the chromatin-based mechanisms involved in the control of CSR.

Is developmental synchrony enabled by CTCF residence time?

Depletion of CTCF in cultured cells has minor effects on transcription whereas its mutation leads to embryonic lethality and developmental defects. In a recent issue of Nature Cell Biology, Soochit et al. (2021) show that the residence time of CTCF on DNA may explain its critical role in cell differentiation.

CTCF chromatin residence time controls three-dimensional genome organization, gene expression and DNA methylation in pluripotent cells

The 11 zinc finger (ZF) protein CTCF regulates topologically associating domain formation and transcription through selective binding to thousands of genomic sites. Here, we replaced endogenous CTCF in mouse embryonic stem cells with green-fluorescent-protein-tagged wild-type or mutant proteins lacking individual ZFs to identify additional determinants of CTCF positioning and function. While ZF1 and ZF8-ZF11 are not essential for cell survival, ZF8 deletion strikingly increases the DNA binding off-rate of mutant CTCF, resulting in reduced CTCF chromatin residence time. Loss of ZF8 results in widespread weakening of topologically associating domains, aberrant gene expression and increased genome-wide DNA methylation. Thus, important chromatin-templated processes rely on accurate CTCF chromatin residence time, which we propose depends on local sequence and chromatin context as well as global CTCF protein concentration.

Power-law behavior of transcription factor dynamics at the single-molecule level implies a continuum affinity model

Single-molecule tracking (SMT) allows the study of transcription factor (TF) dynamics in the nucleus, giving important information regarding the diffusion and binding behavior of these proteins in the nuclear environment. Dwell time distributions obtained by SMT for most TFs appear to follow bi-exponential behavior. This has been ascribed to two discrete populations of TFs-one non-specifically bound to chromatin and another specifically bound to target sites, as implied by decades of biochemical studies. However, emerging studies suggest alternate models for dwell-time distributions, indicating the existence of more than two populations of TFs (multi-exponential distribution), or even the absence of discrete states altogether (power-law distribution). Here, we present an analytical pipeline to evaluate which model best explains SMT data. We find that a broad spectrum of TFs (including glucocorticoid receptor, oestrogen receptor, FOXA1, CTCF) follow a power-law distribution of dwell-times, blurring the temporal line between non-specific and specific binding, suggesting that productive binding may involve longer binding events than previously believed. From these observations, we propose a continuum of affinities model to explain TF dynamics, that is consistent with complex interactions of TFs with multiple nuclear domains as well as binding and searching on the chromatin template.

Dynamic CTCF binding directly mediates interactions among cis-regulatory elements essential for hematopoiesis

While constitutive CCCTC-binding factor (CTCF)-binding sites are needed to maintain relatively invariant chromatin structures, such as topologically associating domains, the precise roles of CTCF to control cell-type-specific transcriptional regulation remain poorly explored. We examined CTCF occupancy in different types of primary blood cells derived from the same donor to elucidate a new role for CTCF in gene regulation during blood cell development. We identified dynamic, cell-type-specific binding sites for CTCF that colocalize with lineage-specific transcription factors. These dynamic sites are enriched for single-nucleotide polymorphisms that are associated with blood cell traits in different linages, and they coincide with the key regulatory elements governing hematopoiesis. CRISPR-Cas9-based perturbation experiments demonstrated that these dynamic CTCF-binding sites play a critical role in red blood cell development. Furthermore, precise deletion of CTCF-binding motifs in dynamic sites abolished interactions of erythroid genes, such as RBM38, with their associated enhancers and led to abnormal erythropoiesis. These results suggest a novel, cell-type-specific function for CTCF in which it may serve to facilitate interaction of distal regulatory emblements with target promoters. Our study of the dynamic, cell-type-specific binding and function of CTCF provides new insights into transcriptional regulation during hematopoiesis.

Resolving diverse protein-DNA footprints from exonuclease-based ChIP experiments

MOTIVATION

High-throughput chromatin immunoprecipitation (ChIP) sequencing-based assays capture genomic regions associated with the profiled transcription factor (TF). ChIP-exo is a modified protocol, which uses lambda exonuclease to digest DNA close to the TF-DNA complex, in order to improve on the positional resolution of the TF-DNA contact. Because the digestion occurs in the 5'-3' orientation, the protocol produces directional footprints close to the complex, on both sides of the double stranded DNA. Like all ChIP-based methods, ChIP-exo reports a mixture of different regions associated with the TF: those bound directly to the TF as well as via intermediaries. However, the distribution of footprints are likely to be indicative of the complex forming at the DNA.

RESULTS

We present ExoDiversity, which uses a model-based framework to learn a joint distribution over footprints and motifs, thus resolving the mixture of ChIP-exo footprints into diverse binding modes. It uses no prior motif or TF information and automatically learns the number of different modes from the data. We show its application on a wide range of TFs and organisms/cell-types. Because its goal is to explain the complete set of reported regions, it is able to identify co-factor TF motifs that appear in a small fraction of the dataset. Further, ExoDiversity discovers small nucleotide variations within and outside canonical motifs, which co-occur with variations in footprints, suggesting that the TF-DNA structural configuration at those regions is likely to be different. Finally, we show that detected modes have specific DNA shape features and conservation signals, giving insights into the structure and function of the putative TF-DNA complexes.

AVAILABILITY AND IMPLEMENTATION

The code for ExoDiversity is available on https://github.com/NarlikarLab/exoDIVERSITY.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genome folding through loop extrusion by SMC complexes

Genomic DNA is folded into loops and topologically associating domains (TADs), which serve important structural and regulatory roles. It has been proposed that these genomic structures are formed by a loop extrusion process, which is mediated by structural maintenance of chromosomes (SMC) protein complexes. Recent single-molecule studies have shown that the SMC complexes condensin and cohesin are indeed able to extrude DNA into loops. In this Review, we discuss how the loop extrusion hypothesis can explain key features of genome architecture; cellular functions of loop extrusion, such as separation of replicated DNA molecules, facilitation of enhancer-promoter interactions and immunoglobulin gene recombination; and what is known about the mechanism of loop extrusion and its regulation, for example, by chromatin boundaries that depend on the DNA binding protein CTCF. We also discuss how the loop extrusion hypothesis has led to a paradigm shift in our understanding of both genome architecture and the functions of SMC complexes.

CTCF mediates dosage- and sequence-context-dependent transcriptional insulation by forming local chromatin domains

Insulators play a critical role in spatiotemporal gene regulation in animals. The evolutionarily conserved CCCTC-binding factor (CTCF) is required for insulator function in mammals, but not all of its binding sites act as insulators. Here we explore the sequence requirements of CTCF-mediated transcriptional insulation using a sensitive insulator reporter in mouse embryonic stem cells. We find that insulation potency depends on the number of CTCF-binding sites in tandem. Furthermore, CTCF-mediated insulation is dependent on upstream flanking sequences at its binding sites. CTCF-binding sites at topologically associating domain boundaries are more likely to function as insulators than those outside topologically associating domain boundaries, independently of binding strength. We demonstrate that insulators form local chromatin domain boundaries and weaken enhancer-promoter contacts. Taken together, our results provide genetic, molecular and structural evidence connecting chromatin topology to the action of insulators in the mammalian genome.

On the choreography of genome folding: A grand pas de deux of cohesin and CTCF

Cohesin and CTCF are key to the 3D folding of interphase chromosomes. Cohesin forms chromatin loops via loop extrusion, a process that involves the formation and enlargement of DNA loops. The architectural protein CTCF controls this process by acting as an anchor for chromatin looping. How CTCF controls cohesin has long been a mystery. Recent work shows that CTCF dictates chromatin looping via a direct interaction of its N-terminus with cohesin. CTCF's ability to regulate chromatin looping turns out to also be partially dependent on several RNA-binding domains. In this review, we discuss recent insights and consider how cohesin and CTCF together may orchestrate the folding of the genome into chromosomal loops.

How DNA loop extrusion mediated by cohesin enables V(D)J recombination

'Structural maintenance of chromosomes' (SMC) complexes are required for the folding of genomic DNA into loops. Theoretical considerations and single-molecule experiments performed with the SMC complexes cohesin and condensin indicate that DNA folding occurs via loop extrusion. Recent work indicates that this process is essential for the assembly of antigen receptor genes by V(D)J recombination in developing B and T cells of the vertebrate immune system. Here, I review how recent studies of the mouse immunoglobulin heavy chain locus Igh have provided evidence for this hypothesis and how the formation of chromatin loops by cohesin and regulation of this process by CTCF and Wapl might ensure that all variable gene segments in this locus (VH segments) participate in recombination with a re-arranged DJH segment, to ensure generation of a maximally diverse repertoire of B-cell receptors and antibodies.

Dual DNA and protein tagging of open chromatin unveils dynamics of epigenomic landscapes in leukemia

The architecture of chromatin regulates eukaryotic cell states by controlling transcription factor access to sites of gene regulation. Here we describe a dual transposase-peroxidase approach, integrative DNA and protein tagging (iDAPT), which detects both DNA (iDAPT-seq) and protein (iDAPT-MS) associated with accessible regions of chromatin. In addition to direct identification of bound transcription factors, iDAPT enables the inference of their gene regulatory networks, protein interactors and regulation of chromatin accessibility. We applied iDAPT to profile the epigenomic consequences of granulocytic differentiation of acute promyelocytic leukemia, yielding previously undescribed mechanistic insights. Our findings demonstrate the power of iDAPT as a platform for studying the dynamic epigenomic landscapes and their transcription factor components associated with biological phenomena and disease.

A Newly Assigned Role of CTCF in Cellular Response to Broken DNAs

Best known as a transcriptional factor, CCCTC-binding factor (CTCF) is a highly conserved multifunctional DNA-binding protein with 11 zinc fingers. It functions in diverse genomic processes, including transcriptional activation/repression, insulation, genome imprinting and three-dimensional genome organization. A big surprise has recently emerged with the identification of CTCF engaging in the repair of DNA double-strand breaks (DSBs) and in the maintenance of genome fidelity. This discovery now adds a new dimension to the multifaceted attributes of this protein. CTCF facilitates the most accurate DSB repair via homologous recombination (HR) that occurs through an elaborate pathway, which entails a chain of timely assembly/disassembly of various HR-repair complexes and chromatin modifications and coordinates multistep HR processes to faithfully restore the original DNA sequences of broken DNA sites. Understanding the functional crosstalks between CTCF and other HR factors will illuminate the molecular basis of various human diseases that range from developmental disorders to cancer and arise from impaired repair. Such knowledge will also help understand the molecular mechanisms underlying the diverse functions of CTCF in genome biology. In this review, we discuss the recent advances regarding this newly assigned versatile role of CTCF and the mechanism whereby CTCF functions in DSB repair.

Distinct properties and functions of CTCF revealed by a rapidly inducible degron system

CCCTC-binding factor (CTCF) is a conserved zinc finger transcription factor implicated in a wide range of functions, including genome organization, transcription activation, and elongation. To explore the basis for CTCF functional diversity, we coupled an auxin-induced degron system with precision nuclear run-on. Unexpectedly, oriented CTCF motifs in gene bodies are associated with transcriptional stalling in a manner independent of bound CTCF. Moreover, CTCF at different binding sites (CBSs) displays highly variable resistance to degradation. Motif sequence does not significantly predict degradation behavior, but location at chromatin boundaries and chromatin loop anchors, as well as co-occupancy with cohesin, are associated with delayed degradation. Single-molecule tracking experiments link chromatin residence time to CTCF degradation kinetics, which has ramifications regarding architectural CTCF functions. Our study highlights the heterogeneity of CBSs, uncovers properties specific to architecturally important CBSs, and provides insights into the basic processes of genome organization and transcription regulation.

On the prediction of DNA-binding preferences of C2H2-ZF domains using structural models: application on human CTCF

Cis2-His2 zinc finger (C2H2-ZF) proteins are the largest family of transcription factors in human and higher metazoans. To date, the DNA-binding preferences of many members of this family remain unknown. We have developed a computational method to predict their DNA-binding preferences. We have computed theoretical position weight matrices (PWMs) of proteins composed by C2H2-ZF domains, with the only requirement of an input structure. We have predicted more than two-third of a single zinc-finger domain binding site for about 70% variants of Zif268, a classical member of this family. We have successfully matched between 60 and 90% of the binding-site motif of examples of proteins composed by three C2H2-ZF domains in JASPAR, a standard database of PWMs. The tests are used as a proof of the capacity to scan a DNA fragment and find the potential binding sites of transcription-factors formed by C2H2-ZF domains. As an example, we have tested the approach to predict the DNA-binding preferences of the human chromatin binding factor CTCF. We offer a server to model the structure of a zinc-finger protein and predict its PWM.

Genome-Wide Histone Modifications and CTCF Enrichment Predict Gene Expression in Sheep Macrophages

Alveolar macrophages function in innate and adaptive immunity, wound healing, and homeostasis in the lungs dependent on tissue-specific gene expression under epigenetic regulation. The functional diversity of tissue resident macrophages, despite their common myeloid lineage, highlights the need to study tissue-specific regulatory elements that control gene expression. Increasing evidence supports the hypothesis that subtle genetic changes alter sheep macrophage response to important production pathogens and zoonoses, for example, viruses like small ruminant lentiviruses and bacteria like Coxiella burnetii. Annotation of transcriptional regulatory elements will aid researchers in identifying genetic mutations of immunological consequence. Here we report the first genome-wide survey of regulatory elements in any sheep immune cell, utilizing alveolar macrophages. We assayed histone modifications and CTCF enrichment by chromatin immunoprecipitation with deep sequencing (ChIP-seq) in two sheep to determine cis-regulatory DNA elements and chromatin domain boundaries that control immunity-related gene expression. Histone modifications included H3K4me3 (denoting active promoters), H3K27ac (active enhancers), H3K4me1 (primed and distal enhancers), and H3K27me3 (broad silencers). In total, we identified 248,674 reproducible regulatory elements, which allowed assignment of putative biological function in macrophages to 12% of the sheep genome. Data exceeded the FAANG and ENCODE standards of 20 million and 45 million useable fragments for narrow and broad marks, respectively. Active elements showed consensus with RNA-seq data and were predictive of gene expression in alveolar macrophages from the publicly available Sheep Gene Expression Atlas. Silencer elements were not enriched for expressed genes, but rather for repressed developmental genes. CTCF enrichment enabled identification of 11,000 chromatin domains with mean size of 258 kb. To our knowledge, this is the first report to use immunoprecipitated CTCF to determine putative topological domains in sheep immune cells. Furthermore, these data will empower phenotype-associated mutation discovery since most causal variants are within regulatory elements.

CTCF As an Example of DNA-Binding Transcription Factors Containing Clusters of C2H2-Type Zinc Fingers

In mammals, most of the boundaries of topologically associating domains and all well-studied insulators are rich in binding sites for the CTCF protein. According to existing experimental data, CTCF is a key factor in the organization of the architecture of mammalian chromosomes. A characteristic feature of the CTCF is that the central part of the protein contains a cluster consisting of eleven domains of C2H2-type zinc fingers, five of which specifically bind to a long DNA sequence conserved in most animals. The class of transcription factors that carry a cluster of C2H2-type zinc fingers consisting of five or more domains (C2H2 proteins) is widely represented in all groups of animals. The functions of most C2H2 proteins still remain unknown. This review presents data on the structure and possible functions of these proteins, using the example of the vertebrate CTCF protein and several well- characterized C2H2 proteins in Drosophila and mammals.

Dynamic regulation of CTCF stability and sub-nuclear localization in response to stress

The nuclear protein CCCTC-binding factor (CTCF) has diverse roles in chromatin architecture and gene regulation. Functionally, CTCF associates with thousands of genomic sites and interacts with proteins, such as cohesin, or non-coding RNAs to facilitate specific transcriptional programming. In this study, we examined CTCF during the cellular stress response in human primary cells using immune-blotting, quantitative real time-PCR, chromatin immunoprecipitation-sequence (ChIP-seq) analysis, mass spectrometry, RNA immunoprecipitation-sequence analysis (RIP-seq), and Airyscan confocal microscopy. Unexpectedly, we found that CTCF is exquisitely sensitive to diverse forms of stress in normal patient-derived human mammary epithelial cells (HMECs). In HMECs, a subset of CTCF protein forms complexes that localize to Serine/arginine-rich splicing factor (SC-35)-containing nuclear speckles. Upon stress, this species of CTCF protein is rapidly downregulated by changes in protein stability, resulting in loss of CTCF from SC-35 nuclear speckles and changes in CTCF-RNA interactions. Our ChIP-seq analysis indicated that CTCF binding to genomic DNA is largely unchanged. Restoration of the stress-sensitive pool of CTCF protein abundance and re-localization to nuclear speckles can be achieved by inhibition of proteasome-mediated degradation. Surprisingly, we observed the same characteristics of the stress response during neuronal differentiation of human pluripotent stem cells (hPSCs). CTCF forms stress-sensitive complexes that localize to SC-35 nuclear speckles during a specific stage of neuronal commitment/development but not in differentiated neurons. We speculate that these particular CTCF complexes serve a role in RNA processing that may be intimately linked with specific genes in the vicinity of nuclear speckles, potentially to maintain cells in a certain differentiation state, that is dynamically regulated by environmental signals. The stress-regulated activity of CTCF is uncoupled in persistently stressed, epigenetically re-programmed "variant" HMECs and certain cancer cell lines. These results reveal new insights into CTCF function in cell differentiation and the stress-response with implications for oxidative damage-induced cancer initiation and neuro-degenerative diseases.

CTCF as a boundary factor for cohesin-mediated loop extrusion: evidence for a multi-step mechanism

Mammalian genome structure is closely linked to function. At the scale of kilobases to megabases, CTCF and cohesin organize the genome into chromatin loops. Mechanistically, cohesin is proposed to extrude chromatin loops bidirectionally until it encounters occupied CTCF DNA-binding sites. Curiously, loops form predominantly between CTCF binding sites in a convergent orientation. How CTCF interacts with and blocks cohesin extrusion in an orientation-specific manner has remained a mechanistic mystery. Here, we review recent papers that have shed light on these processes and suggest a multi-step interaction between CTCF and cohesin. This interaction may first involve a pausing step, where CTCF halts cohesin extrusion, followed by a stabilization step of the CTCF-cohesin complex, resulting in a chromatin loop. Finally, we discuss our own recent studies on an internal RNA-Binding Region (RBRi) in CTCF to elucidate its role in regulating CTCF clustering, target search mechanisms and chromatin loop formation and future challenges.

Epigenomic Reprogramming as a Driver of Malignant Glioma

Malignant gliomas are central nervous system tumors and remain among the most treatment-resistant cancers. Exome sequencing has revealed significant heterogeneity and important insights into the molecular pathogenesis of gliomas. Mutations in chromatin modifiers-proteins that shape the epigenomic landscape through remodeling and regulation of post-translational modifications on chromatin-are very frequent and often define specific glioma subtypes. This suggests that epigenomic reprogramming may be a fundamental driver of glioma. Here, we describe the key chromatin regulatory pathways disrupted in gliomas, delineating their physiological function and our current understanding of how their dysregulation may contribute to gliomagenesis.

Predicting 3D genome folding from DNA sequence with Akita

In interphase, the human genome sequence folds in three dimensions into a rich variety of locus-specific contact patterns. Cohesin and CTCF (CCCTC-binding factor) are key regulators; perturbing the levels of either greatly disrupts genome-wide folding as assayed by chromosome conformation capture methods. Still, how a given DNA sequence encodes a particular locus-specific folding pattern remains unknown. Here we present a convolutional neural network, Akita, that accurately predicts genome folding from DNA sequence alone. Representations learned by Akita underscore the importance of an orientation-specific grammar for CTCF binding sites. Akita learns predictive nucleotide-level features of genome folding, revealing effects of nucleotides beyond the core CTCF motif. Once trained, Akita enables rapid in silico predictions. Accounting for this, we demonstrate how Akita can be used to perform in silico saturation mutagenesis, interpret eQTLs, make predictions for structural variants and probe species-specific genome folding. Collectively, these results enable decoding genome function from sequence through structure.

Predicting 3D genome folding from DNA sequence with Akita

In interphase, the human genome sequence folds in three dimensions into a rich variety of locus-specific contact patterns. Cohesin and CTCF (CCCTC-binding factor) are key regulators; perturbing the levels of either greatly disrupts genome-wide folding as assayed by chromosome conformation capture methods. Still, how a given DNA sequence encodes a particular locus-specific folding pattern remains unknown. Here we present a convolutional neural network, Akita, that accurately predicts genome folding from DNA sequence alone. Representations learned by Akita underscore the importance of an orientation-specific grammar for CTCF binding sites. Akita learns predictive nucleotide-level features of genome folding, revealing effects of nucleotides beyond the core CTCF motif. Once trained, Akita enables rapid in silico predictions. Accounting for this, we demonstrate how Akita can be used to perform in silico saturation mutagenesis, interpret eQTLs, make predictions for structural variants and probe species-specific genome folding. Collectively, these results enable decoding genome function from sequence through structure.

The Functions and Mechanisms of Action of Insulators in the Genomes of Higher Eukaryotes

The mechanisms underlying long-range interactions between chromatin regions and the principles of chromosomal architecture formation are currently under extensive scrutiny. A special class of regulatory elements known as insulators is believed to be involved in the regulation of specific long-range interactions between enhancers and promoters. This review focuses on the insulators of Drosophila and mammals, and it also briefly characterizes the proteins responsible for their functional activity. It was initially believed that the main properties of insulators are blocking of enhancers and the formation of independent transcription domains. We present experimental data proving that the chromatin loops formed by insulators play only an auxiliary role in enhancer blocking. The review also discusses the mechanisms involved in the formation of topologically associating domains and their role in the formation of the chromosomal architecture and regulation of gene transcription.

CTCF orchestrates long-range cohesin-driven V(D)J recombinational scanning

The RAG endonuclease initiates Igh locus V(D)J recombination in progenitor (pro)-B cells1. Upon binding a recombination centre-based JH, RAG scans upstream chromatin via loop extrusion, potentially mediated by cohesin, to locate Ds and assemble a DJH-based recombination centre2. CTCF looping factor-bound elements (CBEs) within IGCR1 upstream of Ds impede RAG scanning3-5; however, their inactivation allows scanning to proximal VHs, where additional CBEs activate rearrangement and impede scanning any further upstream5. Distal VH utilization is thought to involve diffusional access to the recombination centre following large-scale Igh locus contraction6-8. Here we test the potential of linear RAG scanning to mediate distal VH usage in G1-arrested v-Abl pro-B cell lines9, which undergo robust D-to-JH but little VH-to-DJH rearrangements, presumably owing to lack of locus contraction2,5. Through an auxin-inducible approach10, we degraded the cohesin component RAD2110-12 or CTCF12,13 in these G1-arrested lines. Degradation of RAD21 eliminated all V(D)J recombination and interactions associated with RAG scanning, except for reecombination centre-located DQ52-to-JH joining, in which synapsis occurs by diffusion2. Remarkably, while degradation of CTCF suppressed most CBE-based chromatin interactions, it promoted robust recombination centre interactions with, and robust VH-to-DJH joining of, distal VHs, with patterns similar to those of 'locus-contracted' primary pro-B cells. Thus, downmodulation of CTCF-bound scanning-impediment activity promotes cohesin-driven RAG scanning across the 2.7-Mb Igh locus.

Cancer-specific CTCF binding facilitates oncogenic transcriptional dysregulation

BACKGROUND

The three-dimensional genome organization is critical for gene regulation and can malfunction in diseases like cancer. As a key regulator of genome organization, CCCTC-binding factor (CTCF) has been characterized as a DNA-binding protein with important functions in maintaining the topological structure of chromatin and inducing DNA looping. Among the prolific binding sites in the genome, several events with altered CTCF occupancy have been reported as associated with effects in physiology or disease. However, hitherto there is no comprehensive survey of genome-wide CTCF binding patterns across different human cancers.

RESULTS

To dissect functions of CTCF binding, we systematically analyze over 700 CTCF ChIP-seq profiles across human tissues and cancers and identify cancer-specific CTCF binding patterns in six cancer types. We show that cancer-specific lost and gained CTCF binding events are associated with altered chromatin interactions, partially with DNA methylation changes, and rarely with sequence mutations. While lost bindings primarily occur near gene promoters, most gained CTCF binding events exhibit enhancer activities and are induced by oncogenic transcription factors. We validate these findings in T cell acute lymphoblastic leukemia cell lines and patient samples and show that oncogenic NOTCH1 induces specific CTCF binding and they cooperatively activate expression of target genes, indicating transcriptional condensation phenomena.

CONCLUSIONS

Specific CTCF binding events occur in human cancers. Cancer-specific CTCF binding can be induced by other transcription factors to regulate oncogenic gene expression. Our results substantiate CTCF binding alteration as a functional epigenomic signature of cancer.

LncRNAs are regulated by chromatin states and affect the skeletal muscle cell differentiation

Objective

This study aims to clarify the mechanisms underlying transcriptional regulation and regulatory roles of lncRNAs in skeletal muscle cell differentiation.

Methods

We analysed the expression patterns of lncRNAs via time‐course RNA‐seq. Then, we further combined the ATAC‐seq and ChIP‐seq to investigate the governing mechanisms of transcriptional regulation of differentially expressed (DE) lncRNAs. Weighted correlation network analysis and GO analysis were conducted to identify the transcription factor (TF)‐lncRNA pairs related to skeletal muscle cell differentiation.

Results

We identified 385 DE lncRNAs during C2C12 differentiation, the transcription of which is determined by chromatin states around their transcriptional start sites. The TF‐lncRNA correlation network showed substantially concordant changes in DE lncRNAs between C2C12 differentiation and satellite cell rapid growth stages. Moreover, the up‐regulated lncRNAs showed a significant decrease following the differentiation capacity of satellite cells, which gradually declines during skeletal muscle development. Notably, inhibition of the lncRNA Atcayos and Trp53cor1 led to the delayed differentiation of satellite cells. Those lncRNAs were significantly up‐regulated during the rapid growth stage of satellite cells (4‐6 weeks) and down‐regulated with reduced differentiation capacity (8‐12 weeks). It confirms that these lncRNAs are positively associated with myogenic differentiation of satellite cells during skeletal muscle development.

Conclusions

This study extends the understanding of mechanisms governing transcriptional regulation of lncRNAs and provides a foundation for exploring their functions in skeletal muscle cell differentiation.

Many facades of CTCF unified by its coding for three-dimensional genome architecture

CCCTC-binding factor (CTCF) is a multifunctional zinc finger protein that is conserved in metazoan species. CTCF is consistently found to play an important role in many diverse biological processes. CTCF/cohesin-mediated active chromatin 'loop extrusion' architects three-dimensional (3D) genome folding. The 3D architectural role of CTCF underlies its multifarious functions, including developmental regulation of gene expression, protocadherin (Pcdh) promoter choice in the nervous system, immunoglobulin (Ig) and T-cell receptor (Tcr) V(D)J recombination in the immune system, homeobox (Hox) gene control during limb development, as well as many other aspects of biology. Here, we review the pleiotropic functions of CTCF from the perspective of its essential role in 3D genome architecture and topological promoter/enhancer selection. We envision the 3D genome as an enormous complex architecture, with tens of thousands of CTCF sites as connecting nodes and CTCF proteins as mysterious bonds that glue together genomic building parts with distinct articulation joints. In particular, we focus on the internal mechanisms by which CTCF controls higher order chromatin structures that manifest its many façades of physiological and pathological functions. We also discuss the dichotomic role of CTCF sites as intriguing 3D genome nodes for seemingly contradictory 'looping bridges' and 'topological insulators' to frame a beautiful magnificent house for a cell's nuclear home.

MEDEA: analysis of transcription factor binding motifs in accessible chromatin

Deciphering the interplay between chromatin accessibility and transcription factor (TF) binding is fundamental to understanding transcriptional regulation, control of cellular states, and the establishment of new phenotypes. Recent genome-wide chromatin accessibility profiling studies have provided catalogs of putative open regions, where TFs can recognize their motifs and regulate gene expression programs. Here, we present motif enrichment in differential elements of accessibility (MEDEA), a computational tool that analyzes high-throughput chromatin accessibility genomic data to identify cell-type-specific accessible regions and lineage-specific motifs associated with TF binding therein. To benchmark MEDEA, we used a panel of reference cell lines profiled by ENCODE and curated by the ENCODE Project Consortium for the ENCODE-DREAM Challenge. By comparing results with RNA-seq data, ChIP-seq peaks, and DNase-seq footprints, we show that MEDEA improves the detection of motifs associated with known lineage specifiers. We then applied MEDEA to 610 ENCODE DNase-seq data sets, where it revealed significant motifs even when absolute enrichment was low and where it identified novel regulators, such as NRF1 in kidney development. Finally, we show that MEDEA performs well on both bulk and single-cell ATAC-seq data. MEDEA is publicly available as part of our Glossary-GENRE suite for motif enrichment analysis.

Defining the relative and combined contribution of CTCF and CTCFL to genomic regulation

BACKGROUND

Ubiquitously expressed CTCF is involved in numerous cellular functions, such as organizing chromatin into TAD structures. In contrast, its paralog, CTCFL, is normally only present in the testis. However, it is also aberrantly expressed in many cancers. While it is known that shared and unique zinc finger sequences in CTCF and CTCFL enable CTCFL to bind competitively to a subset of CTCF binding sites as well as its own unique locations, the impact of CTCFL on chromosome organization and gene expression has not been comprehensively analyzed in the context of CTCF function. Using an inducible complementation system, we analyze the impact of expressing CTCFL and CTCF-CTCFL chimeric proteins in the presence or absence of endogenous CTCF to clarify the relative and combined contribution of CTCF and CTCFL to chromosome organization and transcription.

RESULTS

We demonstrate that the N terminus of CTCF interacts with cohesin which explains the requirement for convergent CTCF binding sites in loop formation. By analyzing CTCF and CTCFL binding in tandem, we identify phenotypically distinct sites with respect to motifs, targeting to promoter/intronic intergenic regions and chromatin folding. Finally, we reveal that the N, C, and zinc finger terminal domains play unique roles in targeting each paralog to distinct binding sites to regulate transcription, chromatin looping, and insulation.

CONCLUSION

This study clarifies the unique and combined contribution of CTCF and CTCFL to chromosome organization and transcription, with direct implications for understanding how their co-expression deregulates transcription in cancer.

S3norm: simultaneous normalization of sequencing depth and signal-to-noise ratio in epigenomic data

Quantitative comparison of epigenomic data across multiple cell types or experimental conditions is a promising way to understand the biological functions of epigenetic modifications. However, differences in sequencing depth and signal-to-noise ratios in the data from different experiments can hinder our ability to identify real biological variation from raw epigenomic data. Proper normalization is required prior to data analysis to gain meaningful insights. Most existing methods for data normalization standardize signals by rescaling either background regions or peak regions, assuming that the same scale factor is applicable to both background and peak regions. While such methods adjust for differences in sequencing depths, they do not address differences in the signal-to-noise ratios across different experiments. We developed a new data normalization method, called S3norm, that normalizes the sequencing depths and signal-to-noise ratios across different data sets simultaneously by a monotonic nonlinear transformation. We show empirically that the epigenomic data normalized by our method, compared to existing methods, can better capture real biological variation, such as impact on gene expression regulation.

High Sensitivity Profiling of Chromatin Structure by MNase-SSP

A complete view of eukaryotic gene regulation requires that we accurately delineate how transcription factors (TFs) and nucleosomes are arranged along linear DNA in a sensitive, unbiased manner. Here we introduce MNase-SSP, a single-stranded sequencing library preparation method for nuclease-digested chromatin that enables simultaneous mapping of TF and nucleosome positions. As a proof of concept, we apply MNase-SSP toward the genome-wide, high-resolution mapping of nucleosome and TF occupancy in murine embryonic stem cells (mESCs). Compared with existing MNase-seq protocols, MNase-SSP markedly enriches for short DNA fragments, enabling detection of binding by subnucleosomal particles and TFs, in addition to nucleosomes. From these same data, we identify multiple, sequence-dependent binding modes of the architectural TF Ctcf and extend this analysis to the TF Nrsf/ Rest. Looking forward, we anticipate that single stranded protocol (SSP) adaptations of any protein-DNA interaction mapping technique (e.g., ChIP-exo and CUT&RUN) will enhance the information content of the resulting data.

Ramani et al. describe MNase-SSP, a single-stranded DNA sequencing library preparation method for profiling chromatin structure. MNase-SSP libraries harbor diminished sequence bias and capture shorter DNA fragments compared to traditional MNase-seq libraries. Applying MNase-SSP to murine embryonic stem cells enables simultaneous analysis of nucleosomal, subnucleosomal, and transcription factor-DNA interactions genome-wide.

A WIZ/Cohesin/CTCF Complex Anchors DNA Loops to Define Gene Expression and Cell Identity

Chromosome structure is a key regulator of gene expression. CTCF and cohesin play critical roles in structuring chromosomes by mediating physical interactions between distant genomic sites. The resulting DNA loops often contain genes and their cis-regulatory elements. Despite the importance of DNA loops in maintaining proper transcriptional regulation and cell identity, there is limited understanding of the molecular mechanisms that regulate their dynamics and function. We report a previously unrecognized role for WIZ (widely interspaced zinc finger-containing protein) in DNA loop architecture and regulation of gene expression. WIZ forms a complex with cohesin and CTCF that occupies enhancers, promoters, insulators, and anchors of DNA loops. Aberrant WIZ function alters cohesin occupancy and increases the number of DNA loop structures in the genome. WIZ is required for proper gene expression and transcriptional insulation. Our results uncover an unexpected role for WIZ in DNA loop architecture, transcriptional control, and maintenance of cell identity.