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

Varying levels of complexity in transcription factor binding motifs

Binding of transcription factors to DNA is one of the keystones of gene regulation. The existence of statistical dependencies between binding site positions is widely accepted, while their relevance for computational predictions has been debated. Building probabilistic models of binding sites that may capture dependencies is still challenging, since the most successful motif discovery approaches require numerical optimization techniques, which are not suited for selecting dependency structures. To overcome this issue, we propose sparse local inhomogeneous mixture (Slim) models that combine putative dependency structures in a weighted manner allowing for numerical optimization of dependency structure and model parameters simultaneously. We find that Slim models yield a substantially better prediction performance than previous models on genomic context protein binding microarray data sets and on ChIP-seq data sets. To elucidate the reasons for the improved performance, we develop dependency logos, which allow for visual inspection of dependency structures within binding sites. We find that the dependency structures discovered by Slim models are highly diverse and highly transcription factor-specific, which emphasizes the need for flexible dependency models. The observed dependency structures range from broad heterogeneities to sparse dependencies between neighboring and non-neighboring binding site positions.

Spotting the enemy within: Targeted silencing of foreign DNA in mammalian genomes by the Krüppel-associated box zinc finger protein family

Tandem C2H2-type zinc finger proteins (ZFPs) constitute the largest transcription factor family in animals. Tandem-ZFPs bind DNA in a sequence-specific manner through arrays of multiple zinc finger domains that allow high flexibility and specificity in target recognition. In tetrapods, a large proportion of tandem-ZFPs contain Krüppel-associated-box (KRAB) repression domains, which are able to induce epigenetic silencing through the KAP1 corepressor. The KRAB-ZFP family continuously amplified in tetrapods through segmental gene duplications, often accompanied by deletions, duplications, and mutations of the zinc finger domains. As a result, tetrapod genomes contain unique sets of KRAB-ZFP genes, consisting of ancient and recently evolved family members. Although several hundred human and mouse KRAB-ZFPs have been identified or predicted, the biological functions of most KRAB-ZFP family members have gone unexplored. Furthermore, the evolutionary forces driving the extraordinary KRAB-ZFP expansion and diversification have remained mysterious for decades. In this review, we highlight recent studies that associate KRAB-ZFPs with the repression of parasitic DNA elements in the mammalian germ line and discuss the hypothesis that the KRAB-ZFP family primarily evolved as an adaptive genomic surveillance system against foreign DNA. Finally, we comment on the computational, genetic, and biochemical challenges of studying KRAB-ZFPs and attempt to predict how these challenges may be soon overcome.

CTCF Recruits Centromeric Protein CENP-E to the Pericentromeric/Centromeric Regions of Chromosomes through Unusual CTCF-Binding Sites

The role of CTCF in stabilizing long-range interactions between chromatin sites essential for maintaining nuclear architecture is well established. Most of these interactions involve recruitment of the cohesin complex to chromatin via CTCF. We find that CTCF also interacts with the centromeric protein CENP-E both in vitro and in vivo. We identified CTCF sites in pericentric/centromeric DNA and found that, early in mitosis, CTCF binds and recruits CENP-E to these sites. Unlike most known CTCF genomic sites, the CTCF-binding sites in the pericentric/centromeric regions interact strongly with the C-terminal fingers of CTCF. Overexpression of a small CENP-E fragment, targeted to these CTCF sites, results in a delay in alignment of some chromosomes during mitosis, suggesting that the recruitment of CENP-E by CTCF is physiologically important. We conclude that CTCF helps recruit CENP-E to the centromere during mitosis and that it may do so through a structure stabilized by the CTCF/CENP-E complex.

Comparative analyses of CTCF and BORIS occupancies uncover two distinct classes of CTCF binding genomic regions

Background

CTCF and BORIS (CTCFL), two paralogous mammalian proteins sharing nearly identical DNA binding domains, are thought to function in a mutually exclusive manner in DNA binding and transcriptional regulation.

Results

Here we show that these two proteins co-occupy a specific subset of regulatory elements consisting of clustered CTCF binding motifs (termed 2xCTSes). BORIS occupancy at 2xCTSes is largely invariant in BORIS-positive cancer cells, with the genomic pattern recapitulating the germline-specific BORIS binding to chromatin. In contrast to the single-motif CTCF target sites (1xCTSes), the 2xCTS elements are preferentially found at active promoters and enhancers, both in cancer and germ cells. 2xCTSes are also enriched in genomic regions that escape histone to protamine replacement in human and mouse sperm. Depletion of the BORIS gene leads to altered transcription of a large number of genes and the differentiation of K562 cells, while the ectopic expression of this CTCF paralog leads to specific changes in transcription in MCF7 cells.

Conclusions

We discover two functionally and structurally different classes of CTCF binding regions, 2xCTSes and 1xCTSes, revealed by their predisposition to bind BORIS. We propose that 2xCTSes play key roles in the transcriptional program of cancer and germ cells.

Electronic supplementary material

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

CRISPR Inversion of CTCF Sites Alters Genome Topology and Enhancer/Promoter Function

CTCF and the associated cohesin complex play a central role in insulator function and higher-order chromatin organization of mammalian genomes. Recent studies identified a correlation between the orientation of CTCF-binding sites (CBSs) and chromatin loops. To test the functional significance of this observation, we combined CRISPR/Cas9-based genomic-DNA-fragment editing with chromosome-conformation-capture experiments to show that the location and relative orientations of CBSs determine the specificity of long-range chromatin looping in mammalian genomes, using protocadherin (Pcdh) and β-globin as model genes. Inversion of CBS elements within the Pcdh enhancer reconfigures the topology of chromatin loops between the distal enhancer and target promoters and alters gene-expression patterns. Thus, although enhancers can function in an orientation-independent manner in reporter assays, in the native chromosome context, the orientation of at least some enhancers carrying CBSs can determine both the architecture of topological chromatin domains and enhancer/promoter specificity. These findings reveal how 3D chromosome architecture can be encoded by linear genome sequences.

A discrete chromatin loop in the mouse Tcra-Tcrd locus shapes the TCRδ and TCRα repertoires

The Tcra-Tcrd locus undergoes V(D)J recombination in CD4⁻CD8⁻double-negative thymocytes and CD4⁺CD8⁺ double-positive thymocytes to generate diverse TCRδ and TCRα repertoires, respectively. Here we reveal a Tcra-Tcrd locus chromatin interaction network in double-negative thymocytes that was formed by interactions between CTCF-binding elements. Disruption of a discrete chromatin loop encompassing Tcrd diversity, joining and constant gene segments allowed a single variable gene segment to frequently contact and rearrange to diversity and joining gene segments and dominate the adult TCRδ repertoire. Disruption of this loop also narrowed the TCRα repertoire, which, we believe, follows as a consequence of the restricted TCRδ repertoire. Hence, a single CTCF-mediated chromatin loop directly regulates TCRδ diversity and indirectly regulates TCRα diversity.

Role of DNA Methylation in Modulating Transcription Factor Occupancy

Although DNA methylation is commonly invoked as a mechanism for transcriptional repression, the extent to which it actively silences transcription factor (TF) occupancy sites in vivo is unknown. To study the role of DNA methylation in the active modulation of TF binding, we quantified the effect of DNA methylation depletion on the genomic occupancy patterns of CTCF, an abundant TF with known methylation sensitivity that is capable of autonomous binding to its target sites in chromatin. Here, we show that the vast majority (>98.5%) of the tens of thousands of unoccupied, methylated CTCF recognition sequences remain unbound upon abrogation of DNA methylation. The small fraction of sites that show methylation-dependent binding in vivo are in turn characterized by highly variable CTCF occupancy across cell types. Our results suggest that DNA methylation is not a primary groundskeeper of genomic TF landscapes, but rather a specialized mechanism for stabilizing intrinsically labile sites.

Targeting CTCF to Control Virus Gene Expression: A Common Theme amongst Diverse DNA Viruses

All viruses target host cell factors for successful life cycle completion. Transcriptional control of DNA viruses by host cell factors is important in the temporal and spatial regulation of virus gene expression. Many of these factors are recruited to enhance virus gene expression and thereby increase virus production, but host cell factors can also restrict virus gene expression and productivity of infection. CCCTC binding factor (CTCF) is a host cell DNA binding protein important for the regulation of genomic chromatin boundaries, transcriptional control and enhancer element usage. CTCF also functions in RNA polymerase II regulation and in doing so can influence co-transcriptional splicing events. Several DNA viruses, including Kaposi's sarcoma-associated herpesvirus (KSHV), Epstein-Barr virus (EBV) and human papillomavirus (HPV) utilize CTCF to control virus gene expression and many studies have highlighted a role for CTCF in the persistence of these diverse oncogenic viruses. CTCF can both enhance and repress virus gene expression and in some cases CTCF increases the complexity of alternatively spliced transcripts. This review article will discuss the function of CTCF in the life cycle of DNA viruses in the context of known host cell CTCF functions.

CCL8 and the Immune Control of Cytomegalovirus in Organ Transplant Recipients

Monitoring of cytomegalovirus cell-mediated immunity is a promising tool for the refinement of preventative and therapeutic strategies posttransplantation. Typically, the interferon-γ response to T cell stimulation is measured. We evaluated a broad range of cytokine and chemokines to better characterize the ex vivo host-response to CMV peptide stimulation. In a cohort of CMV viremic organ transplant recipients, chemokine expression-specifically CCL8 (AUC 0.849 95% CI 0.721-0.978; p = 0.003) and CXCL10 (AUC 0.841, 95% CI 0.707-0.974; p = 0.004)-was associated with control of viral replication. In a second cohort of transplant recipients at high-risk for CMV, the presence of a polymorphism in the CCL8 promoter conferred an increased risk of viral replication after discontinuation of antiviral prophylaxis (logrank hazard ratio 3.6; 95% CI 2.077-51.88). Using cell-sorting experiments, we determined that the primary cell type producing CCL8 in response to CMV peptide stimulation was the monocyte fraction. Finally, in vitro experiments using standard immunosuppressive agents demonstrated a dose-dependent reduction in CCL8 production. Chemokines appear to be important elements of the cell-mediated response to CMV infection posttransplant, as here suggested for CCL8, and translation of this knowledge may allow for the tailoring and improvement of preventative strategies.

Regulatory Elements in Vectors for Efficient Generation of Cell Lines Producing Target Proteins

To date, there has been an increasing number of drugs produced in mammalian cell cultures. In order to enhance the expression level and stability of target recombinant proteins in cell cultures, various regulatory elements with poorly studied mechanisms of action are used. In this review, we summarize and discuss the potential mechanisms of action of such regulatory elements.

A possible role of Drosophila CTCF in mitotic bookmarking and maintaining chromatin domains during the cell cycle

Background

The CCCTC-binding factor (CTCF) is a highly conserved insulator protein that plays various roles in many cellular processes. CTCF is one of the main architecture proteins in higher eukaryotes, and in combination with other architecture proteins and regulators, also shapes the three-dimensional organization of a genome. Experiments show CTCF partially remains associated with chromatin during mitosis. However, the role of CTCF in the maintenance and propagation of genome architectures throughout the cell cycle remains elusive.

Results

We performed a comprehensive bioinformatics analysis on public datasets of Drosophila CTCF (dCTCF). We characterized dCTCF-binding sites according to their occupancy status during the cell cycle, and identified three classes: interphase-mitosis-common (IM), interphase-only (IO) and mitosis-only (MO) sites. Integrated function analysis showed dCTCF-binding sites of different classes might be involved in different biological processes, and IM sites were more conserved and more intensely bound. dCTCF-binding sites of the same class preferentially localized closer to each other, and were highly enriched at chromatin syntenic and topologically associating domains boundaries.

Conclusions

Our results revealed different functions of dCTCF during the cell cycle and suggested that dCTCF might contribute to the establishment of the three-dimensional architecture of the Drosophila genome by maintaining local chromatin compartments throughout the whole cell cycle.

Electronic supplementary material

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

The genomic landscape of high hyperdiploid childhood acute lymphoblastic leukemia

High hyperdiploid (51-67 chromosomes) acute lymphoblastic leukemia (ALL) is one of the most common childhood malignancies, comprising 30% of all pediatric B cell-precursor ALL. Its characteristic genetic feature is the nonrandom gain of chromosomes X, 4, 6, 10, 14, 17, 18 and 21, with individual trisomies or tetrasomies being seen in over 75% of cases, but the pathogenesis remains poorly understood. We performed whole-genome sequencing (WGS) (n = 16) and/or whole-exome sequencing (WES) (n = 39) of diagnostic and remission samples from 51 cases of high hyperdiploid ALL to further define the genomic landscape of this malignancy. The majority of cases showed involvement of the RTK-RAS pathway and of histone modifiers. No recurrent fusion gene-forming rearrangement was found, and an analysis of mutations on trisomic chromosomes indicated that the chromosomal gains were early events, strengthening the notion that the high hyperdiploid pattern is the main driver event in this common pediatric malignancy.

Role of CTCF poly(ADP-Ribosyl)ation in the regulation of apoptosis in breast cancer cells

INTRODUCTION

CTCF is a candidate tumor suppressor gene encoding a multifunctional transcription factor. CTCF function is controlled by posttranslational modification and interaction with other proteins. Research findings suggested that CTCF function can be regulated by poly(ADP-ribosyl)ation (PARlation) and has specific anti-apoptotic function in breast cancer cells. The aim of this study is to assess the effect of CTCF-wild type (WT) and CTCF complete mutant, which is deficient of PARlation in regulating apoptosis in breast cancer cells.

MATERIALS AND METHODS

The effect of CTCF-WT and CTCF complete mutant was expressed in breast cancer cell-lines by DNA-mediated transfection technique monitored by enhanced green fluorescent protein fluorescence. Evaluation of apoptotic cell death was carried out with immunohistochemical staining using 4'-6'-diamino-2 phenylindole and scoring by fluorescent microscopy.

RESULTS

CTCF-WT supports survival of breast cancer cells and was observed that CTCF complete mutant interferes with the functions of the CTCF-WT and there was a considerable apoptotic cell death in the breast cancer cell lines such as MDA-MB-435, CAMA-1 and MCF-7.

CONCLUSION

The study enlighten CTCF as a Biological Marker for breast cancer and the role of CTCF PARlation may be involved in breast carcinogenesis.

The lncRNA Firre anchors the inactive X chromosome to the nucleolus by binding CTCF and maintains H3K27me3 methylation

BACKGROUND

In mammals, X chromosome genes are present in one copy in males and two in females. To balance the dosage of X-linked gene expression between the sexes, one of the X chromosomes in females is silenced. X inactivation is initiated by upregulation of the lncRNA (long non-coding RNA) Xist and recruitment of specific chromatin modifiers. The inactivated X chromosome becomes heterochromatic and visits a specific nuclear compartment adjacent to the nucleolus.

RESULTS

Here, we show a novel role for the lncRNA Firre in anchoring the inactive mouse X chromosome and preserving one of its main epigenetic features, H3K27me3. Similar to Dxz4, Firre is X-linked and expressed from a macrosatellite repeat locus associated with a cluster of CTCF and cohesin binding sites, and is preferentially located adjacent to the nucleolus. CTCF binding present initially in both male and female mouse embryonic stem cells is lost from the active X during development. Knockdown of Firre disrupts perinucleolar targeting and H3K27me3 levels in mouse fibroblasts, demonstrating a role in maintenance of an important epigenetic feature of the inactive X chromosome. No X-linked gene reactivation is seen after Firre knockdown; however, a compensatory increase in the expression of chromatin modifier genes implicated in X silencing is observed. Further experiments in female embryonic stem cells suggest that Firre does not play a role in X inactivation onset.

CONCLUSIONS

The X-linked lncRNA Firre helps to position the inactive X chromosome near the nucleolus and to preserve one of its main epigenetic features.

The interaction of DNA with multi-Cys2His2 zinc finger proteins

The multi-Cys2His2 (mC2H2) zinc finger protein, like CTCF, plays a central role in the three-dimensional organization of chromatin and gene regulation. The interaction between DNA and mC2H2 zinc finger proteins becomes crucial to better understand how CTCF dynamically shapes the chromatin structure. Here, we study a coarse-grained model of the mC2H2 zinc finger proteins in complexes with DNA, and in particular, we study how a mC2H2 zinc finger protein binds to and searches for its target DNA loci. On the basis of coarse-grained molecular dynamics simulations, we present several interesting kinetic conformational properties of the proteins, such as the rotation-coupled sliding, the asymmetrical roles of different zinc fingers and the partial binding partial dangling mode. In addition, two kinds of studied mC2H2 zinc finger proteins, of CG-rich and AT-rich binding motif each, were able to recognize their target sites and slide away from their non-target sites, which shows a proper sequence specificity in our model and the derived force field for mC2H2-DNA interaction. A further application to CTCF shows that the protein binds to a specific DNA duplex only with its central zinc fingers. The zinc finger domains of CTCF asymmetrically bend the DNA, but do not form a DNA loop alone in our simulations.

Three-dimensional regulation of transcription

Cells can adapt to environment and development by reconstructing their transcriptional networks to regulate diverse cellular processes without altering the underlying DNA sequences. These alterations, namely epigenetic changes, occur during cell division, differentiation and cell death. Numerous evidences demonstrate that epigenetic changes are governed by various types of determinants, including DNA methylation patterns, histone posttranslational modification signatures, histone variants, chromatin remodeling, and recently discovered chromosome conformation characteristics and non-coding RNAs (ncRNAs). Here, we highlight recent efforts on how the two latter epigenetic factors participate in the sophisticated transcriptional process and describe emerging techniques which permit us to uncover and gain insights into the fascinating genomic regulation.

Genome-wide targeting of the epigenetic regulatory protein CTCF to gene promoters by the transcription factor TFII-I

CCCTC-binding factor (CTCF) is a key regulator of nuclear chromatin structure and gene regulation. The impact of CTCF on transcriptional output is highly varied, ranging from repression to transcriptional pausing and transactivation. The multifunctional nature of CTCF may be directed solely through remodeling chromatin architecture. However, another hypothesis is that the multifunctional nature of CTCF is mediated, in part, through differential association with protein partners having unique functions. Consistent with this hypothesis, our mass spectrometry analyses of CTCF interacting partners reveal a previously undefined association with the transcription factor general transcription factor II-I (TFII-I). Biochemical fractionation of CTCF indicates that a distinct CTCF complex incorporating TFII-I is assembled on DNA. Unexpectedly, we found that the interaction between CTCF and TFII-I is essential for directing CTCF to the promoter proximal regulatory regions of target genes across the genome, particularly at genes involved in metabolism. At genes coregulated by CTCF and TFII-I, we find knockdown of TFII-I results in diminished CTCF binding, lack of cyclin-dependent kinase 8 (CDK8) recruitment, and an attenuation of RNA polymerase II phosphorylation at serine 5. Phenotypically, knockdown of TFII-I alters the cellular response to metabolic stress. Our data indicate that TFII-I directs CTCF binding to target genes, and in turn the two proteins cooperate to recruit CDK8 and enhance transcription initiation.

Locus-specific targeting to the X chromosome revealed by the RNA interactome of CTCF

CTCF is a master regulator that plays important roles in genome architecture and gene expression. How CTCF is recruited in a locus-specific manner is not fully understood. Evidence from epigenetic processes, such as X chromosome inactivation (XCI), indicates that CTCF associates functionally with RNA. Using genome-wide approaches to investigate the relationship between its RNA interactome and epigenomic landscape, here we report that CTCF binds thousands of transcripts in mouse embryonic stem cells, many in close proximity to CTCF's genomic binding sites. CTCF is a specific and high-affinity RNA-binding protein (Kd < 1 nM). During XCI, CTCF differentially binds the active and inactive X chromosomes and interacts directly with Tsix, Xite, and Xist RNAs. Tsix and Xite RNAs target CTCF to the X inactivation center, thereby inducing homologous X chromosome pairing. Our work elucidates one mechanism by which CTCF is recruited in a locus-specific manner and implicates CTCF-RNA interactions in long-range chromosomal interactions.

Genomic discovery of potent chromatin insulators for human gene therapy

Insertional mutagenesis and genotoxicity, which usually manifest as hematopoietic malignancy, represent major barriers to realizing the promise of gene therapy. Although insulator sequences that block transcriptional enhancers could mitigate or eliminate these risks, so far no human insulators with high functional potency have been identified. Here we describe a genomic approach for the identification of compact sequence elements that function as insulators. These elements are highly occupied by the insulator protein CTCF, are DNase I hypersensitive and represent only a small minority of the CTCF recognition sequences in the human genome. We show that the elements identified acted as potent enhancer blockers and substantially decreased the risk of tumor formation in a cancer-prone animal model. The elements are small, can be efficiently accommodated by viral vectors and have no detrimental effects on viral titers. The insulators we describe here are expected to increase the safety of gene therapy for genetic diseases.

DNA recombination. Recombination initiation maps of individual human genomes

DNA double-strand breaks (DSBs) are introduced in meiosis to initiate recombination and generate crossovers, the reciprocal exchanges of genetic material between parental chromosomes. Here, we present high-resolution maps of meiotic DSBs in individual human genomes. Comparing DSB maps between individuals shows that along with DNA binding by PRDM9, additional factors may dictate the efficiency of DSB formation. We find evidence for both GC-biased gene conversion and mutagenesis around meiotic DSB hotspots, while frequent colocalization of DSB hotspots with chromosome rearrangement breakpoints implicates the aberrant repair of meiotic DSBs in genomic disorders. Furthermore, our data indicate that DSB frequency is a major determinant of crossover rate. These maps provide new insights into the regulation of meiotic recombination and the impact of meiotic recombination on genome function.

MACE: model based analysis of ChIP-exo

Understanding the role of a given transcription factor (TF) in regulating gene expression requires precise mapping of its binding sites in the genome. Chromatin immunoprecipitation-exo, an emerging technique using λ exonuclease to digest TF unbound DNA after ChIP, is designed to reveal transcription factor binding site (TFBS) boundaries with near-single nucleotide resolution. Although ChIP-exo promises deeper insights into transcription regulation, no dedicated bioinformatics tool exists to leverage its advantages. Most ChIP-seq and ChIP-chip analytic methods are not tailored for ChIP-exo, and thus cannot take full advantage of high-resolution ChIP-exo data. Here we describe a novel analysis framework, termed MACE (model-based analysis of ChIP-exo) dedicated to ChIP-exo data analysis. The MACE workflow consists of four steps: (i) sequencing data normalization and bias correction; (ii) signal consolidation and noise reduction; (iii) single-nucleotide resolution border peak detection using the Chebyshev Inequality and (iv) border matching using the Gale-Shapley stable matching algorithm. When applied to published human CTCF, yeast Reb1 and our own mouse ONECUT1/HNF6 ChIP-exo data, MACE is able to define TFBSs with high sensitivity, specificity and spatial resolution, as evidenced by multiple criteria including motif enrichment, sequence conservation, direct sequence pileup, nucleosome positioning and open chromatin states. In addition, we show that the fundamental advance of MACE is the identification of two boundaries of a TFBS with high resolution, whereas other methods only report a single location of the same event. The two boundaries help elucidate the in vivo binding structure of a given TF, e.g. whether the TF may bind as dimers or in a complex with other co-factors.

High-resolution digital profiling of the epigenome

The widespread adoption of short-read DNA sequencing as a digital epigenomic readout platform has motivated the development of genome-wide tools that achieve base-pair resolution. New methods for footprinting and affinity purification of nucleosomes, RNA polymerases, chromatin remodellers and transcription factors have increased the resolution of epigenomic profiling by two orders of magnitude, leading to new insights into how the chromatin landscape affects gene regulation. These digital epigenomic tools have also been applied to directly profile both turnover kinetics and transcription in situ. In this Review, we describe how these new genome-wide tools allow interrogation of diverse aspects of the epigenome.

The role of CCCTC-binding factor (CTCF) in genomic imprinting, development, and reproduction

CCCTC-binding factor (CTCF) is the major protein involved in insulator activity in vertebrates, with widespread DNA binding sites in the genome. CTCF participates in many processes related to global chromatin organization and remodeling, contributing to the repression or activation of gene transcription. It is also involved in epigenetic reprogramming and is essential during gametogenesis and embryo development. Abnormal DNA methylation patterns at CTCF motifs may impair CTCF binding to DNA, and are related to fertility disorders in mammals. Therefore, CTCF and its binding sites are important candidate regions to be investigated as molecular markers for gamete and embryo quality. This article reviews the role of CTCF in genomic imprinting, gametogenesis, and early embryo development and, moreover, highlights potential opportunities for environmental influences associated with assisted reproductive techniques (ARTs) to affect CTCF-mediated processes. We discuss the potential use of CTCF as a molecular marker for assessing gamete and embryo quality in the context of improving the efficiency and safety of ARTs.

Epigenetic dysregulation by nickel through repressive chromatin domain disruption

Investigations into the genomic landscape of histone modifications in heterochromatic regions have revealed histone H3 lysine 9 dimethylation (H3K9me2) to be important for differentiation and maintaining cell identity. H3K9me2 is associated with gene silencing and is organized into large repressive domains that exist in close proximity to active genes, indicating the importance of maintenance of proper domain structure. Here we show that nickel, a nonmutagenic environmental carcinogen, disrupted H3K9me2 domains, resulting in the spreading of H3K9me2 into active regions, which was associated with gene silencing. We found weak CCCTC-binding factor (CTCF)-binding sites and reduced CTCF binding at the Ni-disrupted H3K9me2 domain boundaries, suggesting a loss of CTCF-mediated insulation function as a potential reason for domain disruption and spreading. We furthermore show that euchromatin islands, local regions of active chromatin within large H3K9me2 domains, can protect genes from H3K9me2-spreading-associated gene silencing. These results have major implications in understanding H3K9me2 dynamics and the consequences of chromatin domain disruption during pathogenesis.

DNase footprint signatures are dictated by factor dynamics and DNA sequence

Genomic footprinting has emerged as an unbiased discovery method for transcription factor (TF) occupancy at cognate DNA in vivo. A basic premise of footprinting is that sequence-specific TF-DNA interactions are associated with localized resistance to nucleases, leaving observable signatures of cleavage within accessible chromatin. This phenomenon is interpreted to imply protection of the critical nucleotides by the stably bound protein factor. However, this model conflicts with previous reports of many TFs exchanging with specific binding sites in living cells on a timescale of seconds. We show that TFs with short DNA residence times have no footprints at bound motif elements. Moreover, the nuclease cleavage profile within a footprint originates from the DNA sequence in the factor-binding site, rather than from the protein occupying specific nucleotides. These findings suggest a revised understanding of TF footprinting and reveal limitations in comprehensive reconstruction of the TF regulatory network using this approach.

Genetic variants within the second intron of the KCNQ1 gene affect CTCF binding and confer a risk of Beckwith-Wiedemann syndrome upon maternal transmission

BACKGROUND

Disruption of 11p15 imprinting results in two fetal growth disorders with opposite phenotypes: the Beckwith-Wiedemann (BWS; MIM 130650) and the Silver-Russell (SRS; MIM 180860) syndromes. DNA methylation defects account for 60% of BWS and SRS cases and, in most cases, occur without any identified mutation in a cis-acting regulatory sequence or a trans-acting factor.

METHODS

We investigated whether 11p15 cis-acting sequence variants account for primary DNA methylation defects in patients with SRS and BWS with loss of DNA methylation at ICR1 and ICR2, respectively.

RESULTS

We identified a 4.5 kb haplotype that, upon maternal transmission, is associated with a risk of ICR2 loss of DNA methylation in patients with BWS. This novel region is located within the second intron of the KCNQ1 gene, 170 kb upstream of the ICR2 imprinting centre and encompasses two CTCF binding sites. We showed that, within the 4.5 kb region, two SNPs (rs11823023 and rs179436) affect CTCF occupancy at DNA motifs flanking the CTCF 20 bp core motif.

CONCLUSIONS

This study shows that genetic variants confer a risk of DNA methylation defect with a parent-of-origin effect and highlights the crucial role of CTCF for the regulation of genomic imprinting of the CDKN1C/KCNQ1 domain.

CTCF haploinsufficiency destabilizes DNA methylation and predisposes to cancer

Epigenetic alterations, particularly in DNA methylation, are ubiquitous in cancer, yet the molecular origins and the consequences of these alterations are poorly understood. CTCF, a DNA-binding protein that regulates higher-order chromatin organization, is frequently altered by hemizygous deletion or mutation in human cancer. To date, a causal role for CTCF in cancer has not been established. Here, we show that Ctcf hemizygous knockout mice are markedly susceptible to spontaneous, radiation-, and chemically induced cancer in a broad range of tissues. Ctcf(+/-) tumors are characterized by increased aggressiveness, including invasion, metastatic dissemination, and mixed epithelial/mesenchymal differentiation. Molecular analysis of Ctcf(+/-) tumors indicates that Ctcf is haploinsufficient for tumor suppression. Tissues with hemizygous loss of CTCF exhibit increased variability in CpG methylation genome wide. These findings establish CTCF as a prominent tumor-suppressor gene and point to CTCF-mediated epigenetic stability as a major barrier to neoplastic progression.

CTCF regulates the human p53 gene through direct interaction with its natural antisense transcript, Wrap53

The multifunctional CCCTC-binding factor (CTCF) protein exhibits a broad range of functions, including that of insulator and higher-order chromatin organizer. We found that CTCF comprises a previously unrecognized region that is necessary and sufficient to bind RNA (RNA-binding region [RBR]) and is distinct from its DNA-binding domain. Depletion of cellular CTCF led to a decrease in not only levels of p53 mRNA, as expected, but also those of Wrap53 RNA, an antisense transcript originated from the p53 locus. PAR-CLIP-seq (photoactivatable ribonucleoside-enhanced cross-linking and immunoprecipitation [PAR-CLIP] combined with deep sequencing) analyses indicate that CTCF binds a multitude of transcripts genome-wide as well as to Wrap53 RNA. Apart from its established role at the p53 promoter, CTCF regulates p53 expression through its physical interaction with Wrap53 RNA. Cells harboring a CTCF mutant in its RBR exhibit a defective p53 response to DNA damage. Moreover, the RBR facilitates CTCF multimerization in an RNA-dependent manner, which may bear directly on its role in establishing higher-order chromatin structures in vivo.

The clustering of CpG islands may constitute an important determinant of the 3D organization of interphase chromosomes

We used the 4C-Seq technique to characterize the genome-wide patterns of spatial contacts of several CpG islands located on chromosome 14 in cultured chicken lymphoid and erythroid cells. We observed a clear tendency for the spatial clustering of CpG islands present on the same and different chromosomes, regardless of the presence or absence of promoters within these CpG islands. Accordingly, we observed preferential spatial contacts between Sp1 binding motifs and other GC-rich genomic elements, including the DNA sequence motifs capable of forming G-quadruplexes. However, an anchor placed in a gene/CpG island-poor area formed spatial contacts with other gene/CpG island-poor areas on chromosome 14 and other chromosomes. These results corroborate the two-compartment model of the spatial organization of interphase chromosomes and suggest that the clustering of CpG islands constitutes an important determinant of the 3D organization of the eukaryotic genome in the cell nucleus. Using the ChIP-Seq technique, we mapped the genome-wide CTCF deposition sites in the chicken lymphoid and erythroid cells that were used for the 4C analysis. We observed a good correlation between the density of CTCF deposition sites and the level of 4C signals for the anchors located in CpG islands but not for an anchor located in a gene desert. It is thus possible that CTCF contributes to the clustering of CpG islands observed in our experiments.

Dual effect of CTCF loss on neuroprogenitor differentiation and survival

An increasing number of proteins involved in genome organization have been implicated in neurodevelopmental disorders, highlighting the importance of chromatin architecture in the developing CNS. The CCCTC-binding factor (CTCF) is a zinc finger DNA binding protein involved in higher-order chromatin organization, and mutations in the human CTCF gene cause an intellectual disability syndrome associated with microcephaly. However, information on CTCF function in vivo in the developing brain is lacking. To address this gap, we conditionally inactivated the Ctcf gene at early stages of mouse brain development. Cre-mediated Ctcf deletion in the telencephalon and anterior retina at embryonic day 8.5 triggered upregulation of the p53 effector PUMA (p53 upregulated modulator of apoptosis), resulting in massive apoptosis and profound ablation of telencephalic structures. Inactivation of Ctcf several days later at E11 also resulted in PUMA upregulation and increased apoptotic cell death, and the Ctcf-null forebrain was hypocellular and disorganized at birth. Although deletion of both Ctcf and Puma in the embryonic brain efficiently rescued Ctcf-null progenitor cell apoptosis, it failed to improve neonatal hypocellularity due to decreased proliferative capacity of rescued apical and outer radial glia progenitor cells. This was exacerbated by an independent effect of CTCF loss that resulted in depletion of the progenitor pool due to premature neurogenesis earlier in development. Our findings demonstrate that CTCF activities are required for two distinct events in early cortex formation: first, to correctly regulate the balance between neuroprogenitor cell proliferation and differentiation, and second, for the survival of neuroprogenitor cells, providing new clues regarding the contributions of CTCF in microcephaly/intellectual disability syndrome pathologies.

CTCF: an architectural protein bridging genome topology and function

The eukaryotic genome is organized in the three-dimensional nuclear space in a specific manner that is both a cause and a consequence of its function. This organization is partly established by a special class of architectural proteins, of which CCCTC-binding factor (CTCF) is the best characterized. Although CTCF has been assigned various roles that are often contradictory, new results now help to draw a unifying model to explain the many functions of this protein. CTCF creates boundaries between topologically associating domains in chromosomes and, within these domains, facilitates interactions between transcription regulatory sequences. Thus, CTCF links the architecture of the genome to its function.

Mechanisms and proteins involved in long-distance interactions

Due to advances in genome-wide technologies, consistent distant interactions within chromosomes of higher eukaryotes have been revealed. In particular, it has been shown that enhancers can specifically and directly interact with promoters by looping out intervening sequences, which can be up to several hundred kilobases long. This review is focused on transcription factors that are supposed to be involved in long-range interactions. Available data are in agreement with the model that several known transcription factors and insulator proteins belong to an abundant but poorly studied class of proteins that are responsible for chromosomal architecture.

On the value of intra-motif dependencies of human insulator protein CTCF

The binding affinity of DNA-binding proteins such as transcription factors is mainly determined by the base composition of the corresponding binding site on the DNA strand. Most proteins do not bind only a single sequence, but rather a set of sequences, which may be modeled by a sequence motif. Algorithms for de novo motif discovery differ in their promoter models, learning approaches, and other aspects, but typically use the statistically simple position weight matrix model for the motif, which assumes statistical independence among all nucleotides. However, there is no clear justification for that assumption, leading to an ongoing debate about the importance of modeling dependencies between nucleotides within binding sites. In the past, modeling statistical dependencies within binding sites has been hampered by the problem of limited data. With the rise of high-throughput technologies such as ChIP-seq, this situation has now changed, making it possible to make use of statistical dependencies effectively. In this work, we investigate the presence of statistical dependencies in binding sites of the human enhancer-blocking insulator protein CTCF by using the recently developed model class of inhomogeneous parsimonious Markov models, which is capable of modeling complex dependencies while avoiding overfitting. These findings lead to a more detailed characterization of the CTCF binding motif, which is only poorly represented by independent nucleotide frequencies at several positions, predominantly at the 3' end.

T-KDE: a method for genome-wide identification of constitutive protein binding sites from multiple ChIP-seq data sets

BACKGROUND

A protein may bind to its target DNA sites constitutively, i.e., regardless of cell type. Intuitively, constitutive binding sites should be biologically functional. A prerequisite for understanding their functional relevance is knowing all their locations for a protein of interest. Genome-wide discovery of constitutive binding sites requires robust and efficient computational methods to integrate results from numerous binding experiments. Such methods are lacking, however.

RESULTS

To locate constitutive binding sites for a protein using ChIP-seq data for that protein from multiple cell lines, we developed a method, T-KDE, which combines a binary range tree with a kernel density estimator. Using 132 CTCF (CCCTC-binding factor) ChIP-seq datasets, we showed that the number of constitutive sites identified by T-KDE is robust to the choice of tuning parameter and that T-KDE identifies binding site locations more accurately than a binning approach. Furthermore, T-KDE can identify constitutive sites that are missed by a motif-based approach either because a bound site failed to reach the motif significance cutoff or because the peak sequence scanned was too short. By studying sites declared constitutive by T-KDE but not by the motif-based approach, we discovered two new CTCF motif variants. Using ENCODE data on 22 transcription factors (TF) in 132 cell lines, we identified constitutive binding sites for each TF and provide evidence that, for some TFs, they may be biologically meaningful.

CONCLUSIONS

T-KDE is an efficient and effective method to predict constitutive protein binding sites using ChIP-seq peaks from multiple cell lines. Besides constitutive binding sites for a given protein, T-KDE can identify genomic "hot spots" where several different proteins bind and, conversely, cell-type-specific sites bound by a given protein.

Sequence specificity incompletely defines the genome-wide occupancy of Myc

BACKGROUND

The Myc-Max heterodimer is a transcription factor that regulates expression of a large number of genes. Genome occupancy of Myc-Max is thought to be driven by Enhancer box (E-box) DNA elements, CACGTG or variants, to which the heterodimer binds in vitro.

RESULTS

By analyzing ChIP-Seq datasets, we demonstrate that the positions occupied by Myc-Max across the human genome correlate with the RNA polymerase II, Pol II, transcription machinery significantly better than with E-boxes. Metagene analyses show that in promoter regions, Myc is uniformly positioned about 100 bp upstream of essentially all promoter proximal paused polymerases with Max about 15 bp upstream of Myc. We re-evaluate the DNA binding properties of full length Myc-Max proteins. Electrophoretic mobility shift assay results demonstrate Myc-Max heterodimers display significant sequence preference, but have high affinity for any DNA. Quantification of the relative affinities of Myc-Max for all possible 8-mers using universal protein-binding microarray assays shows that sequences surrounding core 6-mers significantly affect binding. Compared to the in vitro sequence preferences,Myc-Max genomic occupancy measured by ChIP-Seq is largely, although not completely, independent of sequence specificity.

CONCLUSIONS

We quantified the affinity of Myc-Max to all possible 8-mers and compared this with the sites of Myc binding across the human genome. Our results indicate that the genomic occupancy of Myc cannot be explained by its intrinsic DNA specificity and suggest that the transcription machinery and associated promoter accessibility play a predominant role in Myc recruitment.

CTCF and BORIS in genome regulation and cancer

CTCF plays a vital role in chromatin structure and function. CTCF is ubiquitously expressed and plays diverse roles in gene regulation, imprinting, insulation, intra/interchromosomal interactions, nuclear compartmentalisation, and alternative splicing. CTCF has a single paralogue, the testes-specific CTCF-like gene (CTCFL)/BORIS. CTCF and BORIS can be deregulated in cancer. The tumour suppressor gene CTCF can be mutated or deleted in cancer, or CTCF DNA binding can be altered by epigenetic changes. BORIS is aberrantly expressed frequently in cancer, leading some to propose a pro-tumourigenic role for BORIS. However, BORIS can inhibit cell proliferation, and is mutated in cancer similarly to CTCF suggesting BORIS activation in cancer may be due to global genetic or epigenetic changes typical of malignant transformation.

The cohesin complex regulates immunoglobulin class switch recombination

Immunoglobulin (Ig) class switch recombination (CSR) is initiated by the transcription-coupled recruitment of activation-induced cytidine deaminase (AID) to switch regions and by the subsequent generation of double-stranded DNA breaks (DSBs). These DNA breaks are ultimately resolved through the nonhomologous end joining (NHEJ) pathway. We show that during CSR, AID associates with subunits of cohesin, a complex previously implicated in sister chromatid cohesion, DNA repair, and the formation of DNA loops between enhancers and promoters. Furthermore, we implicate the cohesin complex in the mechanism of CSR by showing that cohesin is dynamically recruited to the Sμ-Cμ region of the IgH locus during CSR and that knockdown of cohesin or its regulatory subunits results in impaired CSR and increased usage of microhomology-based end joining.

CTCF and CTCFL mRNA expression in 17β-estradiol-treated MCF7 cells

Estrogens play a key role in breast cancer, with 60-70% of the cases expressing estrogen receptors (ERs), which are encoded by the ESR1 gene. CTCFL, a paralogue of the chromatin organizer CTCF, is a potential biomarker of breast cancer, but its expression in this disease is currently controversial. A positive correlation has been reported between CTCFL and ERs in breast tumors and there also exists a coordinated interaction between CTCF and ERs in breast cancer cells. Therefore, there appears to be an association between CTCF, CTCFL and estrogens in breast cancer; however, there has been no report on the effects of estrogens on CTCF and CTCFL expression. The aim of this study was to determine the effect of 17β-estradiol (E2) on the CTCF and CTCFL mRNA expression in the MCF7 breast cancer cell line. The promoter methylation status of CTCFL and data mining for estrogen response elements in promoters of the CTCF and CTCFL genes were also determined. The transcription of CTCF and CTCFL was performed by quantitative polymerase chain reaction (qPCR) and the promoter methylation status of CTCFL was determined by methylation-specific PCR. The MCF7 cells exhibited basal transcription of CTCF, which was significantly downregulated to 0.68 by 1 μM E2; basal or E2-regulated transcription of CTCFL was not detected. Under basal conditions, the CTCFL promoter was methylated. Through data mining, an estrogen response element was identified in the CTCF promoter, but no such element was found in CTCFL. These results suggested that estrogens may modulate CTCF expression, although there was no apparent association between ERs and CTCFL.

CTCF binding site sequence differences are associated with unique regulatory and functional trends during embryonic stem cell differentiation

CTCF (CCCTC-binding factor) is a highly conserved multifunctional DNA-binding protein with thousands of binding sites genome-wide. Our previous work suggested that differences in CTCF’s binding site sequence may affect the regulation of CTCF recruitment and its function. To investigate this possibility, we characterized changes in genome-wide CTCF binding and gene expression during differentiation of mouse embryonic stem cells. After separating CTCF sites into three classes (LowOc, MedOc and HighOc) based on similarity to the consensus motif, we found that developmentally regulated CTCF binding occurs preferentially at LowOc sites, which have lower similarity to the consensus. By measuring the affinity of CTCF for selected sites, we show that sites lost during differentiation are enriched in motifs associated with weaker CTCF binding in vitro. Specifically, enrichment for T at the 18^th position of the CTCF binding site is associated with regulated binding in the LowOc class and can predictably reduce CTCF affinity for binding sites. Finally, by comparing changes in CTCF binding with changes in gene expression during differentiation, we show that LowOc and HighOc sites are associated with distinct regulatory functions. Our results suggest that the regulatory control of CTCF is dependent in part on specific motifs within its binding site.

Selective regulation of lymphopoiesis and leukemogenesis by individual zinc fingers of Ikaros

C2H2 zinc fingers are found in several transcriptional regulators in the immune system. However, these proteins usually contain more fingers than are needed for stable DNA binding, suggesting that different fingers regulate different genes and functions. Mice lacking finger 1 or finger 4 of Ikaros exhibited distinct subsets of the phenotypes of Ikaros-null mice. Most notably, the two fingers controlled different stages of lymphopoiesis and finger 4 was selectively required for tumor suppression. The distinct phenotypes suggest that only a small number of Ikaros target genes are critical for each of its biological functions. Subdivision of phenotypes and targets by mutagenesis of individual fingers will facilitate efforts to understand how members of this prevalent family regulate development, immunity and disease.

Characterization of constitutive CTCF/cohesin loci: a possible role in establishing topological domains in mammalian genomes

BACKGROUND

Recent studies suggested that human/mammalian genomes are divided into large, discrete domains that are units of chromosome organization. CTCF, a CCCTC binding factor, has a diverse role in genome regulation including transcriptional regulation, chromosome-boundary insulation, DNA replication, and chromatin packaging. It remains unclear whether a subset of CTCF binding sites plays a functional role in establishing/maintaining chromatin topological domains.

RESULTS

We systematically analysed the genomic, transcriptomic and epigenetic profiles of the CTCF binding sites in 56 human cell lines from ENCODE. We identified ~24,000 CTCF sites (referred to as constitutive sites) that were bound in more than 90% of the cell lines. Our analysis revealed: 1) constitutive CTCF loci were located in constitutive open chromatin and often co-localized with constitutive cohesin loci; 2) most constitutive CTCF loci were distant from transcription start sites and lacked CpG islands but were enriched with the full-spectrum CTCF motifs: a recently reported 33/34-mer and two other potentially novel (22/26-mer); 3) more importantly, most constitutive CTCF loci were present in CTCF-mediated chromatin interactions detected by ChIA-PET and these pair-wise interactions occurred predominantly within, but not between, topological domains identified by Hi-C.

CONCLUSIONS

Our results suggest that the constitutive CTCF sites may play a role in organizing/maintaining the recently identified topological domains that are common across most human cells.