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

Crossed wires: 3D genome misfolding in human disease

Norton and Phillips-Cremins review the 3D architecture of the genome and discuss links between chromatin misfolding and human diseases.

Mammalian genomes are folded into unique topological structures that undergo precise spatiotemporal restructuring during healthy development. Here, we highlight recent advances in our understanding of how the genome folds inside the 3D nucleus and how these folding patterns are miswired during the onset and progression of mammalian disease states. We discuss potential mechanisms underlying the link among genome misfolding, genome dysregulation, and aberrant cellular phenotypes. We also discuss cases in which the endogenous 3D genome configurations in healthy cells might be particularly susceptible to mutation or translocation. Together, these data support an emerging model in which genome folding and misfolding is critically linked to the onset and progression of a broad range of human diseases.

Myc Regulates Chromatin Decompaction and Nuclear Architecture during B Cell Activation

50 years ago, Vincent Allfrey and colleagues discovered that lymphocyte activation triggers massive acetylation of chromatin. However, the molecular mechanisms driving epigenetic accessibility are still unknown. We here show that stimulated lymphocytes decondense chromatin by three differentially regulated steps. First, chromatin is repositioned away from the nuclear periphery in response to global acetylation. Second, histone nanodomain clusters decompact into mononucleosome fibers through a mechanism that requires Myc and continual energy input. Single-molecule imaging shows that this step lowers transcription factor residence time and non-specific collisions during sampling for DNA targets. Third, chromatin interactions shift from long range to predominantly short range, and CTCF-mediated loops and contact domains double in numbers. This architectural change facilitates cognate promoter-enhancer contacts and also requires Myc and continual ATP production. Our results thus define the nature and transcriptional impact of chromatin decondensation and reveal an unexpected role for Myc in the establishment of nuclear topology in mammalian cells.

MeCP2 and CTCF: enhancing the cross-talk of silencers

This paper provides a brief introductory review of the most recent advances in our knowledge about the structural and functional aspects of two transcriptional regulators: MeCP2, a protein whose mutated forms are involved in Rett syndrome; and CTCF, a constitutive transcriptional insulator. This is followed by a description of the PTMs affecting these two proteins and an analysis of their known interacting partners. A special emphasis is placed on the recent studies connecting these two proteins, focusing on the still poorly understood potential structural and functional interactions between the two of them on the chromatin substrate. An overview is provided for some of the currently known genes that are dually regulated by these two proteins. Finally, a model is put forward to account for their possible involvement in their regulation of gene expression.

The β-NAD+ salvage pathway and PKC-mediated signaling influence localized PARP-1 activity and CTCF Poly(ADP)ribosylation

Poly(ADP)ribosylation (PARylation) of the chromatin architectural protein CTCF is critical for CTCF-dependent regulation of chromatin boundary and insulator elements. Loss of CTCF PARylation results in epigenetic silencing of certain tumor suppressor genes through destabilization of nearby chromatin boundaries. We investigated the metabolic and mechanistic processes that regulate PARP-1-mediated CTCF PARylation in human cancer cell lines and discovered a key role for the expression and activity of β-NAD+ salvage enzymes, NAMPT and NMNAT-1. These enzymes are downregulated in cells that exhibit reduced CTCF PARylation, resulting in a decreased concentration of nuclear β-NAD+. In these cells, decreased NMNAT-1 expression is enforced by a proteasome-mediated feedback loop resulting in degradation of NMNAT-1, transcriptional repression of NAMPT, and suppression of PARP-1 activity. Interestingly, dePARylated CTCF is associated in a stable protein complex with PARP-1 and NMNAT-1 in cancer cells harboring silenced tumor suppressor genes. Although the metabolic context in these cells favors suppression of PARP-1 activity, CTCF PARylation can be restored by Protein Kinase C (PKC) signaling. PKC induces dissociation of the catalytically inactive PARP-1/NMNAT-1/CTCF protein complex and phosphorylation of NMNAT-1, which stimulates its proteasome-mediated degradation. Our findings suggest that CTCF PARylation is underpinned by a cellular metabolic context engendered by regulation of the β-NAD+ salvage pathway in which NMNAT-1 acts as a rheostat to control localized β-NAD+ synthesis at CTCF/PARP-1 complexes.

Tissue-specific CTCF-cohesin-mediated chromatin architecture delimits enhancer interactions and function in vivo

The genome is organized via CTCF-cohesin-binding sites, which partition chromosomes into 1-5 megabase (Mb) topologically associated domains (TADs), and further into smaller sub-domains (sub-TADs). Here we examined in vivo an ∼80 kb sub-TAD, containing the mouse α-globin gene cluster, lying within a ∼1 Mb TAD. We find that the sub-TAD is flanked by predominantly convergent CTCF-cohesin sites that are ubiquitously bound by CTCF but only interact during erythropoiesis, defining a self-interacting erythroid compartment. Whereas the α-globin regulatory elements normally act solely on promoters downstream of the enhancers, removal of a conserved upstream CTCF-cohesin boundary extends the sub-TAD to adjacent upstream CTCF-cohesin-binding sites. The α-globin enhancers now interact with the flanking chromatin, upregulating expression of genes within this extended sub-TAD. Rather than acting solely as a barrier to chromatin modification, CTCF-cohesin boundaries in this sub-TAD delimit the region of chromatin to which enhancers have access and within which they interact with receptive promoters.

Genome Organization Drives Chromosome Fragility

In this study, we show that evolutionarily conserved chromosome loop anchors bound by CCCTC-binding factor (CTCF) and cohesin are vulnerable to DNA double strand breaks (DSBs) mediated by topoisomerase 2B (TOP2B). Polymorphisms in the genome that redistribute CTCF/cohesin occupancy rewire DNA cleavage sites to novel loop anchors. While transcription- and replication-coupled genomic rearrangements have been well documented, we demonstrate that DSBs formed at loop anchors are largely transcription-, replication-, and cell-type-independent. DSBs are continuously formed throughout interphase, are enriched on both sides of strong topological domain borders, and frequently occur at breakpoint clusters commonly translocated in cancer. Thus, loop anchors serve as fragile sites that generate DSBs and chromosomal rearrangements. VIDEO ABSTRACT.

CTCF binding landscape in jawless fish with reference to Hox cluster evolution

The nuclear protein CCCTC-binding factor (CTCF) contributes as an insulator to chromatin organization in animal genomes. Currently, our knowledge of its binding property is confined mainly to mammals. In this study, we identified CTCF homologs in extant jawless fishes and performed ChIP-seq for the CTCF protein in the Arctic lamprey. Our phylogenetic analysis suggests that the lamprey lineage experienced gene duplication that gave rise to its unique paralog, designated CTCF2, which is independent from the previously recognized duplication between CTCF and CTCFL. The ChIP-seq analysis detected comparable numbers of CTCF binding sites between lamprey, chicken, and human, and revealed that the lamprey CTCF protein binds to the two-part motif, consisting of core and upstream motifs previously reported for mammals. These findings suggest that this mode of CTCF binding was established in the last common ancestor of extant vertebrates (more than 500 million years ago). We analyzed CTCF binding inside Hox clusters, which revealed a reinforcement of CTCF binding in the region spanning Hox1-4 genes that is unique to lamprey. Our study provides not only biological insights into the antiquity of CTCF-based epigenomic regulation known in mammals but also a technical basis for comparative epigenomic studies encompassing the whole taxon Vertebrata.

CTCF orchestrates the germinal centre transcriptional program and prevents premature plasma cell differentiation

In germinal centres (GC) mature B cells undergo intense proliferation and immunoglobulin gene modification before they differentiate into memory B cells or long-lived plasma cells (PC). GC B-cell-to-PC transition involves a major transcriptional switch that promotes a halt in cell proliferation and the production of secreted immunoglobulins. Here we show that the CCCTC-binding factor (CTCF) is required for the GC reaction in vivo, whereas in vitro the requirement for CTCF is not universal and instead depends on the pathways used for B-cell activation. CTCF maintains the GC transcriptional programme, allows a high proliferation rate, and represses the expression of Blimp-1, the master regulator of PC differentiation. Restoration of Blimp-1 levels partially rescues the proliferation defect of CTCF-deficient B cells. Thus, our data reveal an essential function of CTCF in maintaining the GC transcriptional programme and preventing premature PC differentiation.

Activated B cells differentiate into antibody-producing plasma cells in the germinal centre in secondary lymphoid organs. Here the authors show that this differentiation process and related transcription programs are modulated by the transcription factor CTCF, partly by suppressing the premature expression of Blimp-1.

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

The multidomain CCCTC-binding factor (CTCF), containing a tandem array of 11 zinc fingers (ZFs), modulates the three-dimensional organization of chromatin. We crystallized the human CTCF DNA-binding domain in complex with a known CTCF-binding site. While ZF2 does not make sequence-specific contacts, each finger of ZF3-7 contacts three bases of the 15-bp consensus sequence. Each conserved nucleotide makes base-specific hydrogen bonds with a particular residue. Most of the variable base pairs within the core sequence also engage in interactions with the protein. These interactions compensate for deviations from the consensus sequence, allowing CTCF to adapt to sequence variations. CTCF is sensitive to cytosine methylation at position 2, but insensitive at position 12 of the 15-bp core sequence. These differences can be rationalized structurally. Although included in crystallizations, ZF10 and ZF11 are not visible, while ZF8 and ZF9 span the backbone of the DNA duplex, conferring no sequence specificity but adding to overall binding stability.

The C2H2-ZF transcription factor Zfp335 recognizes two consensus motifs using separate zinc finger arrays

Here, Han et al. show that transcription factor Zfp335 binds DNA and drives transcription via recognition of two distinct consensus motifs by separate ZF clusters and identify the specific motif interaction disrupted by the mutation R1092W. This study presents Zfp335 as a model for understanding how C2H2-ZF TFs may use multiple recognition motifs to control gene expression.

The complexities of DNA recognition by transcription factors (TFs) with multiple Cys2–His2 zinc fingers (C2H2-ZFs) remain poorly studied. We previously reported a mutation (R1092W) in the C2H2-ZF TF Zfp335 that led to selective loss of binding at a subset of targets, although the basis for this effect was unclear. We show that Zfp335 binds DNA and drives transcription via recognition of two distinct consensus motifs by separate ZF clusters and identify the specific motif interaction disrupted by R1092W. Our work presents Zfp335 as a model for understanding how C2H2-ZF TFs may use multiple recognition motifs to control gene expression.

Direct Observation of Cell-Cycle-Dependent Interactions between CTCF and Chromatin

The three-dimensional arrangement of chromatin encodes regulatory traits important for nuclear processes such as transcription and replication. Chromatin topology is in part mediated by the architectural protein CCCTC-binding factor (CTCF) that binds to the boundaries of topologically associating domains. Whereas sites of CTCF interactions are well characterized, little is known on how long CTCF binds to chromatin and how binding evolves during the cell cycle. We monitored CTCF-chromatin interactions by live cell single molecule tracking in different phases of the cell cycle. In G1-, S-, and G2-phases, a majority of CTCF molecules was bound transiently (∼0.2 s) to chromatin, whereas minor fractions were bound dynamically (∼4 s) or stably (>15 min). During mitosis, CTCF was mostly excluded from chromatin. Our data suggest that CTCF scans DNA in search for two different subsets of specific target sites and provide information on the timescales over which topologically associating domains might be restructured. During S-phase, dynamic and stable interactions decreased considerably compared to G1-phase, but were resumed in G2-phase, indicating that specific interactions need to be dissolved for replication to proceed.

CTCF and cohesin regulate chromatin loop stability with distinct dynamics

Folding of mammalian genomes into spatial domains is critical for gene regulation. The insulator protein CTCF and cohesin control domain location by folding domains into loop structures, which are widely thought to be stable. Combining genomic and biochemical approaches we show that CTCF and cohesin co-occupy the same sites and physically interact as a biochemically stable complex. However, using single-molecule imaging we find that CTCF binds chromatin much more dynamically than cohesin (~1-2 min vs. ~22 min residence time). Moreover, after unbinding, CTCF quickly rebinds another cognate site unlike cohesin for which the search process is long (~1 min vs. ~33 min). Thus, CTCF and cohesin form a rapidly exchanging 'dynamic complex' rather than a typical stable complex. Since CTCF and cohesin are required for loop domain formation, our results suggest that chromatin loops are dynamic and frequently break and reform throughout the cell cycle.

New factors in mammalian DNA repair-the chromatin connection

In response to DNA damage mammalian cells activate a complex network of stress response pathways collectively termed DNA damage response (DDR). DDR involves a temporary arrest of the cell cycle to allow for the repair of the damage. DDR also attenuates gene expression by silencing global transcription and translation. Main function of DDR is, however, to prevent the fixation of debilitating changes to DNA by activation of various DNA repair pathways. Proper execution of DDR requires careful coordination between these interdependent cellular responses. Deregulation of some aspects of DDR orchestration is potentially pathological and could lead to various undesired outcomes such as DNA translocations, cellular transformation or acute cell death. It is thus critical to understand the regulation of DDR in cells especially in the light of a strong linkage between the DDR impairment and the occurrence of common human diseases such as cancer. In this review we focus on recent advances in understanding of mammalian DNA repair regulation and a on the function of PAXX/c9orf142 and ZNF281 proteins that recently had been discovered to play a role in that process. We focus on regulation of double-strand DNA break (DSB) repair via the non-homologous end joining pathway, as unrepaired DSBs are the primary cause of pathological cellular states after DNA damage. Interestingly these new factors operate at the level of chromatin, which reinforces a notion of a central role of chromatin structure in the regulation of cellular DDR regulation.

C2H2 Zinc Finger Proteins: The Largest but Poorly Explored Family of Higher Eukaryotic Transcription Factors

The emergence of whole-genome assays has initiated numerous genome-wide studies of transcription factor localizations at genomic regulatory elements (enhancers, promoters, silencers, and insulators), as well as facilitated the uncovering of some of the key principles of chromosomal organization. However, the proteins involved in the formation and maintenance of the chromosomal architecture and the organization of regulatory domains remain insufficiently studied. This review attempts to collate the available data on the abundant but still poorly understood family of proteins with clusters of the C2H2 zinc finger domains. One of the best known proteins of this family is a well conserved protein known as CTCF, which plays a key role in the establishment of the chromosomal architecture in vertebrates. The distinctive features of C2H2 zinc finger proteins include strong and specific binding to a long and unique DNA recognition target sequence and rapid expansion within various animal taxa during evolution. The reviewed data support a proposed model according to which many of the C2H2 proteins have functions that are similar to those of the CTCF in the organization of the chromatin architecture.

In vivo binding of PRDM9 reveals interactions with noncanonical genomic sites

In mouse and human meiosis, DNA double-strand breaks (DSBs) initiate homologous recombination and occur at specific sites called hotspots. The localization of these sites is determined by the sequence-specific DNA binding domain of the PRDM9 histone methyl transferase. Here, we performed an extensive analysis of PRDM9 binding in mouse spermatocytes. Unexpectedly, we identified a noncanonical recruitment of PRDM9 to sites that lack recombination activity and the PRDM9 binding consensus motif. These sites include gene promoters, where PRDM9 is recruited in a DSB-dependent manner. Another subset reveals DSB-independent interactions between PRDM9 and genomic sites, such as the binding sites for the insulator protein CTCF. We propose that these DSB-independent sites result from interactions between hotspot-bound PRDM9 and genomic sequences located on the chromosome axis.

CTCF genetic alterations in endometrial carcinoma are pro-tumorigenic

CTCF is a haploinsufficient tumour suppressor gene with diverse normal functions in genome structure and gene regulation. However the mechanism by which CTCF haploinsufficiency contributes to cancer development is not well understood. CTCF is frequently mutated in endometrial cancer. Here we show that most CTCF mutations effectively result in CTCF haploinsufficiency through nonsense-mediated decay of mutant transcripts, or loss-of-function missense mutation. Conversely, we identified a recurrent CTCF mutation K365T, which alters a DNA binding residue, and acts as a gain-of-function mutation enhancing cell survival. CTCF genetic deletion occurs predominantly in poor prognosis serous subtype tumours, and this genetic deletion is associated with poor overall survival. In addition, we have shown that CTCF haploinsufficiency also occurs in poor prognosis endometrial clear cell carcinomas and has some association with endometrial cancer relapse and metastasis. Using shRNA targeting CTCF to recapitulate CTCF haploinsufficiency, we have identified a novel role for CTCF in the regulation of cellular polarity of endometrial glandular epithelium. Overall, we have identified two novel pro-tumorigenic roles (promoting cell survival and altering cell polarity) for genetic alterations of CTCF in endometrial cancer.

CTCF participates in DNA damage response via poly(ADP-ribosyl)ation

CCCTC-binding factor (CTCF) plays an essential role in regulating the structure of chromatin by binding DNA strands for defining the boundary between active and heterochromatic DNA. However, the role of CTCF in DNA damage response remains elusive. Here, we show that CTCF is quickly recruited to the sites of DNA damage. The fast recruitment is mediated by the zinc finger domain and poly (ADP-ribose) (PAR). Further analyses show that only three zinc finger motifs are essential for PAR recognition. Moreover, CTCF-deficient cells are hypersensitive to genotoxic stress such as ionizing radiation (IR). Taken together, these results suggest that CTCF participate in DNA damage response via poly(ADP-ribosylation).

Cohesin biology meets the loop extrusion model

Extensive research has revealed that cohesin acts as a topological device, trapping chromosomal DNA within a large tripartite ring. In so doing, cohesin contributes to the formation of compact and organized genomes. How exactly the cohesin subunits interact, how it opens, closes, and translocates on chromatin, and how it actually tethers DNA strands together are still being elucidated. A comprehensive understanding of these questions will shed light on how cohesin performs its many functions, including its recently proposed role as a chromatid loop extruder. Here, we discuss this possibility in light of our understanding of the molecular properties of cohesin complexes.

Mitotic phosphorylation of CCCTC-binding factor (CTCF) reduces its DNA binding activity

During mitosis, higher order chromatin structures are disrupted and chromosomes are condensed to achieve accurate chromosome segregation. CCCTC‐binding factor (CTCF) is a highly conserved and ubiquitously expressed C2H2‐type zinc finger protein which is considered to be involved in epigenetic memory through regulation of higher order chromatin architecture. However, the regulatory mechanism of CTCF in mitosis is still unclear. Here we found that the DNA‐binding activity of CTCF is regulated in a phosphorylation‐dependent manner during mitosis. The linker domains of the CTCF zinc finger domain were found to be phosphorylated during mitosis. The phosphorylation of linker domains impaired the DNA‐binding activity in vitro. Mutation analyses showed that amino acid residues (Thr289, Thr317, Thr346, Thr374, Ser402, Ser461, and Thr518) located in the linker domains were phosphorylated during mitosis. Based on these results, we propose that the mitotic phosphorylation of the linker domains of CTCF is important for the dissociation of CTCF from mitotic chromatin.

A dynamic mode of mitotic bookmarking by transcription factors

During mitosis, transcription is shut off, chromatin condenses, and most transcription factors (TFs) are reported to be excluded from chromosomes. How do daughter cells re-establish the original transcription program? Recent discoveries that a select set of TFs remain bound on mitotic chromosomes suggest a potential mechanism for maintaining transcriptional programs through the cell cycle termed mitotic bookmarking. Here we report instead that many TFs remain associated with chromosomes in mouse embryonic stem cells, and that the exclusion previously described is largely a fixation artifact. In particular, most TFs we tested are significantly enriched on mitotic chromosomes. Studies with Sox2 reveal that this mitotic interaction is more dynamic than in interphase and is facilitated by both DNA binding and nuclear import. Furthermore, this dynamic mode results from lack of transcriptional activation rather than decreased accessibility of underlying DNA sequences in mitosis. The nature of the cross-linking artifact prompts careful re-examination of the role of TFs in mitotic bookmarking.

DOI:http://dx.doi.org/10.7554/eLife.22280.001

A kidney cell functions differently from a skin cell despite the fact that all the cells in one organism share the same DNA. This is because not all of the genes encoded within the DNA are active in the cells. Instead, cells can turn on just those genes that are specific to how that cell type works. One way that cells can regulate their genes is by using proteins called transcription factors that can bind to DNA to turn nearby genes on and off.

When cells divide to form new cells, the DNA is condensed and gene activity is turned off. However, each dividing cell also has to ‘remember’ the program of genes that specifies its identity. After division, how do the cells know which genes to turn on and which ones to keep off?

It was thought that the transcription factors attached to the DNA were all detached from it during cell division. Through studies in mouse embryonic stem cells, Teves et al. now show that this finding is largely an artifact of the methods used to study the process. In fact, many transcription factors still bind to and interact with DNA during cell division. This provides an efficient way for the newly formed cells to quickly reset to the pattern of gene activity appropriate for their cell type.

Having found that many key transcription factors are still bound to DNA during cell division, the next challenge is to find out what role this binding plays in allowing cells to ‘remember’ their identity.

DOI:http://dx.doi.org/10.7554/eLife.22280.002

CTCF and Cohesin in Genome Folding and Transcriptional Gene Regulation

Genome function, replication, integrity, and propagation rely on the dynamic structural organization of chromosomes during the cell cycle. Genome folding in interphase provides regulatory segmentation for appropriate transcriptional control, facilitates ordered genome replication, and contributes to genome integrity by limiting illegitimate recombination. Here, we review recent high-resolution chromosome conformation capture and functional studies that have informed models of the spatial and regulatory compartmentalization of mammalian genomes, and discuss mechanistic models for how CTCF and cohesin control the functional architecture of mammalian chromosomes.

Topoisomerase II beta interacts with cohesin and CTCF at topological domain borders

BACKGROUND

Type II DNA topoisomerases (TOP2) regulate DNA topology by generating transient double stranded breaks during replication and transcription. Topoisomerase II beta (TOP2B) facilitates rapid gene expression and functions at the later stages of development and differentiation. To gain new insight into the genome biology of TOP2B, we used proteomics (BioID), chromatin immunoprecipitation, and high-throughput chromosome conformation capture (Hi-C) to identify novel proximal TOP2B protein interactions and characterize the genomic landscape of TOP2B binding at base pair resolution.

RESULTS

Our human TOP2B proximal protein interaction network included members of the cohesin complex and nucleolar proteins associated with rDNA biology. TOP2B associates with DNase I hypersensitivity sites, allele-specific transcription factor (TF) binding, and evolutionarily conserved TF binding sites on the mouse genome. Approximately half of all CTCF/cohesion-bound regions coincided with TOP2B binding. Base pair resolution ChIP-exo mapping of TOP2B, CTCF, and cohesin sites revealed a striking structural ordering of these proteins along the genome relative to the CTCF motif. These ordered TOP2B-CTCF-cohesin sites flank the boundaries of topologically associating domains (TADs) with TOP2B positioned externally and cohesin internally to the domain loop.

CONCLUSIONS

TOP2B is positioned to solve topological problems at diverse cis-regulatory elements and its occupancy is a highly ordered and prevalent feature of CTCF/cohesin binding sites that flank TADs.

CTCF: making the right connections

Ghirlando and Felsenfeld review recent major advances in understanding the multiple roles of CTCF in gene regulation and genome organization and especially in how such domains are generated and organized.

The role of the zinc finger protein CTCF in organizing the genome within the nucleus is now well established. Widely separated sites on DNA, occupied by both CTCF and the cohesin complex, make physical contacts that create large loop domains. Additional contacts between loci within those domains, often also mediated by CTCF, tend to be favored over contacts between loci in different domains. A large number of studies during the past 2 years have addressed the questions: How are these loops generated? What are the effects of disrupting them? Are there rules governing large-scale genome organization? It now appears that the strongest and evolutionarily most conserved of these CTCF interactions have specific rules for the orientation of the paired CTCF sites, implying the existence of a nonequilibrium mechanism of generation. Recent experiments that invert, delete, or inactivate one of a mating CTCF pair result in major changes in patterns of organization and gene expression in the surrounding regions. What remain to be determined are the detailed molecular mechanisms for re-establishing loop domains and maintaining them after replication and mitosis. As recently published data show, some mechanisms may involve interactions with noncoding RNAs as well as protein cofactors. Many CTCF sites are also involved in other functions such as modulation of RNA splicing and specific regulation of gene expression, and the relationship between these activities and loop formation is another unanswered question that should keep investigators occupied for some time.

PAtCh-Cap: input strategy for improving analysis of ChIP-exo data sets and beyond

Recently, a number of advances have been implemented into the core ChIP-seq (chromatin immunoprecipitation coupled with next-generation sequencing) methodology to streamline the process, reduce costs or improve data resolution. Several of these emerging ChIP-based methods perform additional chemical steps on bead-bound immunoprecipitated chromatin, posing a challenge for generating similarly treated input controls required for artifact removal during bioinformatics analyses. Here we present a versatile method for producing technique-specific input controls for ChIP-based methods that utilize additional bead-bound processing steps. This reported method, termed protein attached chromatin capture (PAtCh-Cap), relies on the non-specific capture of chromatin-bound proteins via their carboxylate groups, leaving the DNA accessible for subsequent chemical treatments in parallel with chromatin separately immunoprecipitated for the target protein. Application of this input strategy not only significantly enhanced artifact removal from ChIP-exo data, increasing confidence in peak identification and allowing for de novo motif searching, but also afforded discovery of a novel CTCF binding motif.

Evaluating the impact of single nucleotide variants on transcription factor binding

Diseases and phenotypes caused by disrupted transcription factor (TF) binding are being identified, but progress is hampered by our limited capacity to predict such functional alterations. Improving predictions may be dependent on expanding the set of bona fide TF binding alterations. Allele-specific binding (ASB) events, where TFs preferentially bind to one of the two alleles at heterozygous sites, reveal the impact of sequence variations in altered TF binding. Here, we present the largest ASB compilation to our knowledge, 10 765 ASB events retrieved from 45 ENCODE ChIP-Seq data sets. Our analysis showed that ASB events were frequently associated with motif alterations of the ChIP'ed TF and potential partner TFs, allelic difference of DNase I hypersensitivity and allelic difference of histone modifications. For TF dimers bound symmetrically to DNA, ASB data revealed that central positions of the TF binding motifs were disproportionately important for binding. Lastly, the impact of variation on TF binding was predicted by a classification model incorporating all the investigated features of ASB events. Classification models using only DNase I hypersensitivity and sequence data exhibited predictive accuracy approaching the models with substantially more features. Taken together, the combination of ASB data and the classification model represents an important step toward elucidating regulatory variants across the human genome.

Genome-wide footprinting: ready for prime time?

High-throughput sequencing technologies have allowed many gene locus-level molecular biology assays to become genome-wide profiling methods. DNA-cleaving enzymes such as DNase I have been used to probe accessible chromatin. The accessible regions contain functional regulatory sites, including promoters, insulators and enhancers. Deep sequencing of DNase-seq libraries and computational analysis of the cut profiles have been used to infer protein occupancy in the genome at the nucleotide level, a method introduced as 'digital genomic footprinting'. The approach has been proposed as an attractive alternative to the analysis of transcription factors (TFs) by chromatin immunoprecipitation followed by sequencing (ChIP-seq), and in theory it should overcome antibody issues, poor resolution and batch effects. Recent reports point to limitations of the DNase-based genomic footprinting approach and call into question the scope of detectable protein occupancy, especially for TFs with short-lived chromatin binding. The genomics community is grappling with issues concerning the utility of genomic footprinting and is reassessing the proposed approaches in terms of robust deliverables. Here we summarize the consensus as well as different views emerging from recent reports, and we describe the remaining issues and hurdles for genomic footprinting.

Bayesian Markov models consistently outperform PWMs at predicting motifs in nucleotide sequences

Position weight matrices (PWMs) are the standard model for DNA and RNA regulatory motifs. In PWMs nucleotide probabilities are independent of nucleotides at other positions. Models that account for dependencies need many parameters and are prone to overfitting. We have developed a Bayesian approach for motif discovery using Markov models in which conditional probabilities of order k - 1 act as priors for those of order k This Bayesian Markov model (BaMM) training automatically adapts model complexity to the amount of available data. We also derive an EM algorithm for de-novo discovery of enriched motifs. For transcription factor binding, BaMMs achieve significantly (P    =  1/16) higher cross-validated partial AUC than PWMs in 97% of 446 ChIP-seq ENCODE datasets and improve performance by 36% on average. BaMMs also learn complex multipartite motifs, improving predictions of transcription start sites, polyadenylation sites, bacterial pause sites, and RNA binding sites by 26-101%. BaMMs never performed worse than PWMs. These robust improvements argue in favour of generally replacing PWMs by BaMMs.

CRISPR/Cas9 genome editing throws descriptive 3-D genome folding studies for a loop

CRISPR/Cas9 genome editing studies have recently shed new light into the causal link between the linear DNA sequence and 3-D chromatin architecture. Here we describe current models for the folding of genomes into a nested hierarchy of higher-order structures and discuss new insights into the organizing principles governing genome folding at each length scale. WIREs Syst Biol Med 2016, 8:286-299. doi: 10.1002/wsbm.1338 For further resources related to this article, please visit the WIREs website.

Chromatin Scanning by Dynamic Binding of Pioneer Factors

Pioneer factors such as FoxA target nucleosomal DNA and initiate cooperative interactions at silent genes during development, cellular reprogramming, and steroid hormone induction. Biophysical studies previously showed that the nuclear mobility of FoxA1 is slower than for many other transcription factors, whereas a new single molecule study (Swinstead et al., 2016, Cell) shows comparable chromatin residence times for FoxA1 and steroid receptors. Despite that steroid receptors engage nucleosome-remodeling complexes, the vast majority of co-bound sites with FoxA are dependent upon FoxA, not vice versa. Taken together, the distinguishing feature of pioneer factors remains nucleosomal access rather than an exceptional residence time in chromatin.

Formation of Chromosomal Domains by Loop Extrusion

Topologically associating domains (TADs) are fundamental structural and functional building blocks of human interphase chromosomes, yet the mechanisms of TAD formation remain unclear. Here, we propose that loop extrusion underlies TAD formation. In this process, cis-acting loop-extruding factors, likely cohesins, form progressively larger loops but stall at TAD boundaries due to interactions with boundary proteins, including CTCF. Using polymer simulations, we show that this model produces TADs and finer-scale features of Hi-C data. Each TAD emerges from multiple loops dynamically formed through extrusion, contrary to typical illustrations of single static loops. Loop extrusion both explains diverse experimental observations-including the preferential orientation of CTCF motifs, enrichments of architectural proteins at TAD boundaries, and boundary deletion experiments-and makes specific predictions for the depletion of CTCF versus cohesin. Finally, loop extrusion has potentially far-ranging consequences for processes such as enhancer-promoter interactions, orientation-specific chromosomal looping, and compaction of mitotic chromosomes.

Defining the multivalent functions of CTCF from chromatin state and three-dimensional chromatin interactions

CCCTC-binding factor (CTCF) is a multi-functional protein that is assigned various, even contradictory roles in the genome. High-throughput sequencing-based technologies such as ChIP-seq and Hi-C provided us the opportunity to assess the multivalent functions of CTCF in the human genome. The location of CTCF-binding sites with respect to genomic features provides insights into the possible roles of this protein. Here we present the first genome-wide survey and characterization of three important functions of CTCF: enhancer insulator, chromatin barrier and enhancer linker. We developed a novel computational framework to discover the multivalent functions of CTCF based on chromatin state and three-dimensional chromatin architecture. We applied our method to five human cell lines and identified ∼46 000 non-redundant CTCF sites related to the three functions. Disparate effects of these functions on gene expression were found and distinct genomic features of these CTCF sites were characterized in GM12878 cells. Finally, we investigated the cell-type specificities of CTCF sites related to these functions across five cell types. Our study provides new insights into the multivalent functions of CTCF in the human genome.

Steroid Receptors Reprogram FoxA1 Occupancy through Dynamic Chromatin Transitions

The estrogen receptor (ER), glucocorticoid receptor (GR), and forkhead box protein 1 (FoxA1) are significant factors in breast cancer progression. FoxA1 has been implicated in establishing ER-binding patterns though its unique ability to serve as a pioneer factor. However, the molecular interplay between ER, GR, and FoxA1 requires further investigation. Here we show that ER and GR both have the ability to alter the genomic distribution of the FoxA1 pioneer factor. Single-molecule tracking experiments in live cells reveal a highly dynamic interaction of FoxA1 with chromatin in vivo. Furthermore, the FoxA1 factor is not associated with detectable footprints at its binding sites throughout the genome. These findings support a model wherein interactions between transcription factors and pioneer factors are highly dynamic. Moreover, at a subset of genomic sites, the role of pioneer can be reversed, with the steroid receptors serving to enhance binding of FoxA1.

Analysis of computational footprinting methods for DNase sequencing experiments

DNase-seq allows nucleotide-level identification of transcription factor binding sites on the basis of a computational search of footprint-like DNase I cleavage patterns on the DNA. Frequently in high-throughput methods, experimental artifacts such as DNase I cleavage bias affect the computational analysis of DNase-seq experiments. Here we performed a comprehensive and systematic study on the performance of computational footprinting methods. We evaluated ten footprinting methods in a panel of DNase-seq experiments for their ability to recover cell-specific transcription factor binding sites. We show that three methods--HINT, DNase2TF and PIQ--consistently outperformed the other evaluated methods and that correcting the DNase-seq signal for experimental artifacts significantly improved the accuracy of computational footprints. We also propose a score that can be used to detect footprints arising from transcription factors with potentially short residence times.

CRISPR Double Cutting through the Labyrinthine Architecture of 3D Genomes

The genomes are organized into ordered and hierarchical topological structures in interphase nuclei. Within discrete territories of each chromosome, topologically associated domains (TADs) play important roles in various nuclear processes such as gene regulation. Inside TADs separated by relatively constitutive boundaries, distal elements regulate their gene targets through specific chromatin-looping contacts such as long-distance enhancer-promoter interactions. High-throughput sequencing studies have revealed millions of potential regulatory DNA elements, which are much more abundant than the mere ∼20,000 genes they control. The recently emerged CRISPR-Cas9 genome editing technologies have enabled efficient and precise genetic and epigenetic manipulations of genomes. The multiplexed and high-throughput CRISPR capabilities facilitate the discovery and dissection of gene regulatory elements. Here, we describe the applications of CRISPR for genome, epigenome, and 3D genome editing, focusing on CRISPR DNA-fragment editing with Cas9 and a pair of sgRNAs to investigate topological folding of chromatin TADs and developmental gene regulation.

CTCF-mediated chromatin loops enclose inducible gene regulatory domains

Background

The CCTC-binding factor (CTCF) protein is involved in genome organization, including mediating three-dimensional chromatin interactions. Human patient lymphocytes with mutations in a single copy of the CTCF gene have reduced expression of enhancer-associated genes involved in response to stimuli. We hypothesize that CTCF interactions stabilize enhancer-promoter chromatin interaction domains, facilitating increased expression of genes in response to stimuli. Here we systematically investigate this model using computational analyses.

Results

We use CTCF ChIA-PET data from the ENCODE project to show that CTCF-associated chromatin loops have a tendency to enclose regions of enhancer-regulated stimulus responsive genes, insulating them from neighboring regions of constitutively expressed housekeeping genes. To facilitate cell type-specific CTCF loop identification, we develop an algorithm to predict CTCF loops from ChIP-seq data alone by exploiting the CTCF motif directionality in loop anchors. We apply this algorithm to a hundred ENCODE cell line datasets, confirming the universality of our observations as well as identifying a general distinction between primary and immortal cells in loop-enclosed gene content. Finally, we combine the existing evidence to propose a model for the formation of CTCF loops in which partner sites are brought together by chromatin template reeling through stationary RNA polymerases, consistent with the transcription factory hypothesis.

Conclusions

We provide computational evidence that CTCF-mediated chromatin interactions enclose domains of stimulus responsive enhancer-regulated genes, insulating them from nearby housekeeping genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2516-6) contains supplementary material, which is available to authorized users.

Epigenetic silencing of miR-181c by DNA methylation in glioblastoma cell lines

Background

Post-transcriptional regulation by microRNAs is recognized as one of the major pathways for the control of cellular homeostasis. Less well understood is the transcriptional and epigenetic regulation of genes encoding microRNAs. In the present study we addressed the epigenetic regulation of the miR-181c in normal and malignant brain cells.

Methods

To explore the epigenetic regulation of the miR-181c we evaluated its expression using RT-qPCR and the in vivo binding of the CCCTC-binding factor (CTCF) to its regulatory region in different glioblastoma cell lines. DNA methylation survey, chromatin immunoprecipitation and RNA interference assays were used to assess the role of CTCF in the miR-181c epigenetic silencing.

Results

We found that miR-181c is downregulated in glioblastoma cell lines, as compared to normal brain tissues. Loss of expression correlated with a notorious gain of DNA methylation at the miR-181c promoter region and the dissociation of the multifunctional nuclear factor CTCF. Taking advantage of the genomic distribution of CTCF in different cell types we propose that CTCF has a local and cell type specific regulatory role over the miR-181c and not an architectural one through chromatin loop formation. This is supported by the depletion of CTCF in glioblastoma cells affecting the expression levels of NOTCH2 as a target of miR-181c.

Conclusion

Together, our results point to the epigenetic role of CTCF in the regulation of microRNAs implicated in tumorigenesis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12885-016-2273-6) contains supplementary material, which is available to authorized users.

Insulated Foamy Viral Vectors

Retroviral vector-mediated gene therapy is promising, but genotoxicity has limited its use in the clinic. Genotoxicity is highly dependent on the retroviral vector used, and foamy viral (FV) vectors appear relatively safe. However, internal promoters may still potentially activate nearby genes. We developed insulated FV vectors, using four previously described insulators: a version of the well-studied chicken hypersensitivity site 4 insulator (650cHS4), two synthetic CCCTC-binding factor (CTCF)-based insulators, and an insulator based on the CCAAT box-binding transcription factor/nuclear factor I (7xCTF/NF1). We directly compared these insulators for enhancer-blocking activity, effect on FV vector titer, and fidelity of transfer to both proviral long terminal repeats. The synthetic CTCF-based insulators had the strongest insulating activity, but reduced titers significantly. The 7xCTF/NF1 insulator did not reduce titers but had weak insulating activity. The 650cHS4-insulated FV vector was identified as the overall most promising vector. Uninsulated and 650cHS4-insulated FV vectors were both significantly less genotoxic than gammaretroviral vectors. Integration sites were evaluated in cord blood CD34(+) cells and the 650cHS4-insulated FV vector had fewer hotspots compared with an uninsulated FV vector. These data suggest that insulated FV vectors are promising for hematopoietic stem cell gene therapy.

The retrovirus HTLV-1 inserts an ectopic CTCF-binding site into the human genome

Human T-lymphotropic virus type 1 (HTLV-1) is a retrovirus that causes malignant and inflammatory diseases in ∼10% of infected people. A typical host has between 10(4) and 10(5) clones of HTLV-1-infected T lymphocytes, each clone distinguished by the genomic integration site of the single-copy HTLV-1 provirus. The HTLV-1 bZIP (HBZ) factor gene is constitutively expressed from the minus strand of the provirus, whereas plus-strand expression, required for viral propagation to uninfected cells, is suppressed or intermittent in vivo, allowing escape from host immune surveillance. It remains unknown what regulates this pattern of proviral transcription and latency. Here, we show that CTCF, a key regulator of chromatin structure and function, binds to the provirus at a sharp border in epigenetic modifications in the pX region of the HTLV-1 provirus in T cells naturally infected with HTLV-1. CTCF is a zinc-finger protein that binds to an insulator region in genomic DNA and plays a fundamental role in controlling higher order chromatin structure and gene expression in vertebrate cells. We show that CTCF bound to HTLV-1 acts as an enhancer blocker, regulates HTLV-1 mRNA splicing, and forms long-distance interactions with flanking host chromatin. CTCF-binding sites (CTCF-BSs) have been propagated throughout the genome by transposons in certain primate lineages, but CTCF binding has not previously been described in present-day exogenous retroviruses. The presence of an ectopic CTCF-BS introduced by the retrovirus in tens of thousands of genomic locations has the potential to cause widespread abnormalities in host cell chromatin structure and gene expression.

Insulators to Improve the Safety of Retroviral Vectors for HIV Gene Therapy

Retroviral vector gene therapy is a promising approach to treating HIV-1. However, integrated vectors are mutagens with the potential to dysregulate nearby genes and cause severe adverse side effects. Leukemia has already been a documented severe adverse event in gene therapy clinical trials for the treatment of primary immunodeficiencies. These side effects will need to be reduced or avoided if retroviral vectors are to be used clinically for HIV-1 treatment. The addition of chromatin insulators to retroviral vectors is a potential strategy for reducing adverse side effects. Insulators have already been effectively used in retroviral vectors to reduce genotoxicity in pre-clinical studies. Here, we will review how insulators function, genotoxicity in gene therapy clinical trials, the design of insulated retroviral vectors, promising results from insulated retroviral vector studies, and considerations for the development of insulated retroviral treatment vectors for HIV-1 gene therapy.

The evolutionary turnover of recombination hot spots contributes to speciation in mice

Meiotic recombination is required for the segregation of homologous chromosomes and is essential for fertility. In most mammals, the DNA double-strand breaks (DSBs) that initiate meiotic recombination are directed to a subset of genomic loci (hot spots) by sequence-specific binding of the PRDM9 protein. Rapid evolution of the DNA-binding specificity of PRDM9 and gradual erosion of PRDM9-binding sites by gene conversion will alter the recombination landscape over time. To better understand the evolutionary turnover of recombination hot spots and its consequences, we mapped DSB hot spots in four major subspecies of Mus musculus with different Prdm9 alleles and in their F1 hybrids. We found that hot spot erosion governs the preferential usage of some Prdm9 alleles over others in hybrid mice and increases sequence diversity specifically at hot spots that become active in the hybrids. As crossovers are disfavored at such hot spots, we propose that sequence divergence generated by hot spot turnover may create an impediment for recombination in hybrids, potentially leading to reduced fertility and, eventually, speciation.

TGF-β/β2-spectrin/CTCF-regulated tumor suppression in human stem cell disorder Beckwith-Wiedemann syndrome

Beckwith-Wiedemann syndrome (BWS) is a human stem cell disorder, and individuals with this disease have a substantially increased risk (~800-fold) of developing tumors. Epigenetic silencing of β2-spectrin (β2SP, encoded by SPTBN1), a SMAD adaptor for TGF-β signaling, is causally associated with BWS; however, a role of TGF-β deficiency in BWS-associated neoplastic transformation is unexplored. Here, we have reported that double-heterozygous Sptbn1+/- Smad3+/- mice, which have defective TGF-β signaling, develop multiple tumors that are phenotypically similar to those of BWS patients. Moreover, tumorigenesis-associated genes IGF2 and telomerase reverse transcriptase (TERT) were overexpressed in fibroblasts from BWS patients and TGF-β-defective mice. We further determined that chromatin insulator CCCTC-binding factor (CTCF) is TGF-β inducible and facilitates TGF-β-mediated repression of TERT transcription via interactions with β2SP and SMAD3. This regulation was abrogated in TGF-β-defective mice and BWS, resulting in TERT overexpression. Imprinting of the IGF2/H19 locus and the CDKN1C/KCNQ1 locus on chromosome 11p15.5 is mediated by CTCF, and this regulation is lost in BWS, leading to aberrant overexpression of growth-promoting genes. Therefore, we propose that loss of CTCF-dependent imprinting of tumor-promoting genes, such as IGF2 and TERT, results from a defective TGF-β pathway and is responsible at least in part for BWS-associated tumorigenesis as well as sporadic human cancers that are frequently associated with SPTBN1 and SMAD3 mutations.

Binding of Protein Factor CTCF within Chicken Genome Alpha-Globin Locus

A systematic search for DNA fragments containing potential CTCF transcription factor binding sites in the chicken alpha-globin domain and its flanking regions was performed by means of the two-dimension electrophoretic mobility shift assay. For the alpha-globin domain fragments selected, the occupancy by the CTCF in erythroid and lymphoid chicken cells was tested by chromatin immunoprecipitation. Only one of 13 DNA fragments capable of CTCF binding in vitro was efficiently bound to this protein in vivo in erythroid cells, and somewhat less efficiently - in lymphoid cells. So, binding of CTCF to the DNA fragment in vitro in most cases does not mean that this fragment will be occupied by CTCF in the cell nucleus. Yet, CTCF binding in vivo, as a rule, is accompanied by the binding of the protein to this DNA region in vitro. During the erythroid differentiation, no significant changes in CTCF binding to the DNA fragments studied were detected.

Towards a predictive model of chromatin 3D organization

Architectural proteins mediate interactions between distant regions in the genome to bring together different regulatory elements while establishing a specific three-dimensional organization of the genetic material. Depletion of specific architectural proteins leads to miss regulation of gene expression and alterations in nuclear organization. The specificity of interactions mediated by architectural proteins depends on the nature, number, and orientation of their binding site at individual genomic locations. Knowledge of the mechanisms and rules governing interactions among architectural proteins may provide a code to predict the 3D organization of the genome.

The ubiquitous transcription factor CTCF promotes lineage-specific epigenomic remodeling and establishment of transcriptional networks driving cell differentiation

Cell differentiation relies on tissue-specific transcription factors (TFs) that cooperate to establish unique transcriptomes and phenotypes. However, the role of ubiquitous TFs in these processes remains poorly defined. Recently, we have shown that the CCCTC-binding factor (CTCF) is required for adipocyte differentiation through epigenomic remodelling of adipose tissue-specific enhancers and transcriptional activation of Peroxisome proliferator-activated receptor gamma (PPARG), the main driver of the adipogenic program (PPARG), and its target genes. Here, we discuss how these findings, together with the recent literature, illuminate a functional role for ubiquitous TFs in lineage-determining transcriptional networks.

CTCF Binding Polarity Determines Chromatin Looping

CCCTC-binding factor (CTCF) is an architectural protein involved in the three-dimensional (3D) organization of chromatin. In this study, we assayed the 3D genomic contact profiles of a large number of CTCF binding sites with high-resolution 4C-seq. As recently reported, our data also suggest that chromatin loops preferentially form between CTCF binding sites oriented in a convergent manner. To directly test this, we used CRISPR/Cas9 genome editing to delete core CTCF binding sites in three loci, including the CTCF site in the Sox2 super-enhancer. In all instances, CTCF and cohesin recruitment were lost, and chromatin loops with distal, convergent CTCF sites were disrupted or destabilized. Re-insertion of oppositely oriented CTCF recognition sequences restored CTCF and cohesin recruitment, but did not re-establish chromatin loops. We conclude that CTCF binding polarity plays a functional role in the formation of higher-order chromatin structure.

DiffLogo: a comparative visualization of sequence motifs

BACKGROUND

For three decades, sequence logos are the de facto standard for the visualization of sequence motifs in biology and bioinformatics. Reasons for this success story are their simplicity and clarity. The number of inferred and published motifs grows with the number of data sets and motif extraction algorithms. Hence, it becomes more and more important to perceive differences between motifs. However, motif differences are hard to detect from individual sequence logos in case of multiple motifs for one transcription factor, highly similar binding motifs of different transcription factors, or multiple motifs for one protein domain.

RESULTS

Here, we present DiffLogo, a freely available, extensible, and user-friendly R package for visualizing motif differences. DiffLogo is capable of showing differences between DNA motifs as well as protein motifs in a pair-wise manner resulting in publication-ready figures. In case of more than two motifs, DiffLogo is capable of visualizing pair-wise differences in a tabular form. Here, the motifs are ordered by similarity, and the difference logos are colored for clarity. We demonstrate the benefit of DiffLogo on CTCF motifs from different human cell lines, on E-box motifs of three basic helix-loop-helix transcription factors as examples for comparison of DNA motifs, and on F-box domains from three different families as example for comparison of protein motifs.

CONCLUSIONS

DiffLogo provides an intuitive visualization of motif differences. It enables the illustration and investigation of differences between highly similar motifs such as binding patterns of transcription factors for different cell types, treatments, and algorithmic approaches.

Inferring intra-motif dependencies of DNA binding sites from ChIP-seq data

Background

Statistical modeling of transcription factor binding sites is one of the classical fields in bioinformatics. The position weight matrix (PWM) model, which assumes statistical independence among all nucleotides in a binding site, used to be the standard model for this task for more than three decades but its simple assumptions are increasingly put into question. Recent high-throughput sequencing methods have provided data sets of sufficient size and quality for studying the benefits of more complex models. However, learning more complex models typically entails the danger of overfitting, and while model classes that dynamically adapt the model complexity to data have been developed, effective model selection is to date only possible for fully observable data, but not, e.g., within de novo motif discovery.

Results

To address this issue, we propose a stochastic algorithm for performing robust model selection in a latent variable setting. This algorithm yields a solution without relying on hyperparameter-tuning via massive cross-validation or other computationally expensive resampling techniques. Using this algorithm for learning inhomogeneous parsimonious Markov models, we study the degree of putative higher-order intra-motif dependencies for transcription factor binding sites inferred via de novo motif discovery from ChIP-seq data. We find that intra-motif dependencies are prevalent and not limited to first-order dependencies among directly adjacent nucleotides, but that second-order models appear to be the significantly better choice.

Conclusions

The traditional PWM model appears to be indeed insufficient to infer realistic sequence motifs, as it is on average outperformed by more complex models that take into account intra-motif dependencies. Moreover, using such models together with an appropriate model selection procedure does not lead to a significant performance loss in comparison with the PWM model for any of the studied transcription factors. Hence, we find it worthwhile to recommend that any modern motif discovery algorithm should attempt to take into account intra-motif dependencies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0797-4) contains supplementary material, which is available to authorized users.

Genetic Tailors: CTCF and Cohesin Shape the Genome During Evolution

Research into chromosome structure and organization is an old field that has seen some fascinating progress in recent years. Modern molecular methods that can describe the shape of chromosomes have begun to revolutionize our understanding of genome organization and the mechanisms that regulate gene activity. A picture is beginning to emerge of chromatin loops representing a widespread organizing principle of the chromatin fiber and the proteins cohesin and CCCTC-binding factor (CTCF) as key players anchoring such chromatin loops. Here we review our current understanding of the features of CTCF- and cohesin-mediated genome organization and how their evolution may have helped to shape genome structure.

Chromosomal Loop Domains Direct the Recombination of Antigen Receptor Genes

RAG initiates antibody V(D)J recombination in developing lymphocytes by generating "on-target" DNA breaks at matched pairs of bona fide recombination signal sequences (RSSs). We employ bait RAG-generated breaks in endogenous or ectopically inserted RSS pairs to identify huge numbers of RAG "off-target" breaks. Such breaks occur at the simple CAC motif that defines the RSS cleavage site and are largely confined within convergent CTCF-binding element (CBE)-flanked loop domains containing bait RSS pairs. Marked orientation dependence of RAG off-target activity within loops spanning up to 2 megabases implies involvement of linear tracking. In this regard, major RAG off-targets in chromosomal translocations occur as convergent RSS pairs at enhancers within a loop. Finally, deletion of a CBE-based IgH locus element disrupts V(D)J recombination domains and, correspondingly, alters RAG on- and off-target distributions within IgH. Our findings reveal how RAG activity is developmentally focused and implicate mechanisms by which chromatin domains harness biological processes within them.