Does CTCF mediate between nuclear organization and gene expression?

C2H2 proteins: Evolutionary aspects of domain architecture and diversification

The largest group of transcription factors in higher eukaryotes are C2H2 proteins, which contain C2H2-type zinc finger domains that specifically bind to DNA. Few well-studied C2H2 proteins, however, demonstrate their key role in the control of gene expression and chromosome architecture. Here we review the features of the domain architecture of C2H2 proteins and the likely origin of C2H2 zinc fingers. A comprehensive investigation of proteomes for the presence of proteins with multiple clustered C2H2 domains has revealed a key difference between groups of organisms. Unlike plants, transcription factors in metazoans contain clusters of C2H2 domains typically separated by a linker with the TGEKP consensus sequence. The average size of C2H2 clusters varies substantially, even between genomes of higher metazoans, and with a tendency to increase in combination with SCAN, and especially KRAB domains, reflecting the increasing complexity of gene regulatory networks.

Atypical Modes of CTCF Binding Facilitate Tissue-Specific and Neuronal Activity-Dependent Gene Expression States

Multivalent binding of CTCF to variable DNA sequences is thought to underlie its ability to mediate diverse cellular functions. CTCF typically binds a 20 base-pair consensus DNA sequence, but the full diversity of CTCF binding sites (CBS) within the genome has not been interrogated. We assessed CTCF occupancy in cultured cortical neurons and observed surprisingly that ~ 22% of CBS lack the consensus CTCF motif. We report here that sequence diversity at most of these atypical CBS involves degeneracy at specific nucleotide positions within the consensus CTCF motif, which likely affect the binding of CTCF zinc fingers 6 and 7. This mode of atypical CTCF binding defines most CBS at gene promoters, as well as CBS that are dynamically altered during neural differentiation and following neuronal stimulation, revealing how atypical CTCF binding could influence gene activity. Dynamic CBS are distributed both within and outside loop anchors and TAD boundaries, suggesting both looping-dependent and independent roles for CTCF. Finally, we describe a second mode of atypical CTCF binding to DNA sequences that are completely unrelated to the consensus CTCF motif, which are enriched within the bodies of tissue-specific genes. These tissue-specific atypical CBS are also enriched in H3K27ac, which marks cis-regulatory elements within chromatin, including enhancers. Overall, these results indicate how atypical CBS could dynamically regulate gene activity patterns during differentiation, development, and in response to environmental cues.

Involvement of CCCTC-binding factor in epigenetic regulation of cancer

A major global health burden continues to be borne by the complex and multifaceted disease of cancer. Epigenetic changes, which are essential for the emergence and spread of cancer, have drawn a huge amount of attention recently. The CCCTC-binding factor (CTCF), which takes part in a wide range of cellular processes including genomic imprinting, X chromosome inactivation, 3D chromatin architecture, local modifications of histone, and RNA polymerase II-mediated gene transcription, stands out among the diverse array of epigenetic regulators. CTCF not only functions as an architectural protein but also modulates DNA methylation and histone modifications. Epigenetic regulation of cancer has already been the focus of plenty of studies. Understanding the role of CTCF in the cancer epigenetic landscape may lead to the development of novel targeted therapeutic strategies for cancer. CTCF has already earned its status as a tumor suppressor gene by acting like a homeostatic regulator of genome integrity and function. Moreover, CTCF has a direct effect on many important transcriptional regulators that control the cell cycle, apoptosis, senescence, and differentiation. As we learn more about CTCF-mediated epigenetic modifications and transcriptional regulations, the possibility of utilizing CTCF as a diagnostic marker and therapeutic target for cancer will also increase. Thus, the current review intends to promote personalized and precision-based therapeutics for cancer patients by shedding light on the complex interplay between CTCF and epigenetic processes.

Chromatin architectural factor CTCF is essential for progesterone-dependent uterine maturation

Receptors for estrogen and progesterone frequently interact, via Cohesin/CTCF loop extrusion, at enhancers distal from regulated genes. Loss-of-function CTCF mutation in >20% of human endometrial tumors indicates its importance in uterine homeostasis. To better understand how CTCF-mediated enhancer-gene interactions impact endometrial development and function, the Ctcf gene was selectively deleted in female reproductive tissues of mice. Prepubertal Ctcfd/d uterine tissue exhibited a marked reduction in the number of uterine glands compared to those without Ctcf deletion (Ctcff/f mice). Post-pubertal Ctcfd/d uteri were hypoplastic with significant reduction in both the amount of the endometrial stroma and number of glands. Transcriptional profiling revealed increased expression of stem cell molecules Lif, EOMES, and Lgr5, and enhanced inflammation pathways following Ctcf deletion. Analysis of the response of the uterus to steroid hormone stimulation showed that CTCF deletion affects a subset of progesterone-responsive genes. This finding indicates (1) Progesterone-mediated signaling remains functional following Ctcf deletion and (2) certain progesterone-regulated genes are sensitive to Ctcf deletion, suggesting they depend on gene-enhancer interactions that require CTCF. The progesterone-responsive genes altered by CTCF ablation included Ihh, Fst, and Errfi1. CTCF-dependent progesterone-responsive uterine genes enhance critical processes including anti-tumorigenesis, which is relevant to the known effectiveness of progesterone in inhibiting progression of early-stage endometrial tumors. Overall, our findings reveal that uterine Ctcf plays a key role in progesterone-dependent expression of uterine genes underlying optimal post-pubertal uterine development.

Epigenomic analysis of Alzheimer's disease brains reveals diminished CTCF binding on genes involved in synaptic organization

Epigenetic aberrations are suggested to play an important role in transcriptional alterations in Alzheimer's disease (AD). One of the key mechanisms of epigenetic regulation of gene expression is through the dynamic organization of chromatin structure via the master genome architecture protein, CCCTC-binding factor (CTCF). By forming chromatin loops, CTCF can influence gene transcription in a complex manner. To find out whether genome-wide DNA binding sites for CTCF are altered in AD, we compared CTCF chromatin immunoprecipitation sequencing (ChIP-Seq) data from frontal cortex of human AD patients and normal controls (n = 9 pairs, all females). We have revealed that CTCF-binding affinity on many genes is significantly reduced in AD patients, and these genes are enriched in synaptic organization, cell adhesion, and actin cytoskeleton, including synaptic scaffolding molecules and receptors, such as SHANK2, HOMER1, NRXN1, CNTNAP2 and GRIN2A, and protocadherin (PCDH) and cadherin (CDH) family members. By comparing transcriptomic data from AD patients, we have discovered that many of the synaptic and adhesion genes with reduced CTCF binding in AD are significantly reduced in their mRNA expression. Moreover, a significant overlap of genes with the diminished CTCF binding and the reduced H3K27ac is identified in AD, with the common genes enriched in synaptic organization. These data suggest that the CTCF-controlled 3D chromatin organization is perturbed in AD, which may be linked to the diminished expression of target genes, probably through changes in histone modification.

Benefiting from the intrinsic role of epigenetics to predict patterns of CTCF binding

Motivation

One of the most relevant mechanisms involved in the determination of chromatin structure is the formation of structural loops that are also related with the conservation of chromatin states. Many of these loops are stabilized by CCCTC-binding factor (CTCF) proteins at their base. Despite the relevance of chromatin structure and the key role of CTCF, the role of the epigenetic factors that are involved in the regulation of CTCF binding, and thus, in the formation of structural loops in the chromatin, is not thoroughly understood.

Results

Here we describe a CTCF binding predictor based on Random Forest that employs different epigenetic data and genomic features. Importantly, given the ability of Random Forests to determine the relevance of features for the prediction, our approach also shows how the different types of descriptors impact the binding of CTCF, confirming previous knowledge on the relevance of chromatin accessibility and DNA methylation, but demonstrating the effect of epigenetic modifications on the activity of CTCF. We compared our approach against other predictors and found improved performance in terms of areas under PR and ROC curves (PRAUC-ROCAUC), outperforming current state-of-the-art methods.

Auxin-inducible degron 2 system deciphers functions of CTCF domains in transcriptional regulation

BACKGROUND

CTCF is a well-established chromatin architectural protein that also plays various roles in transcriptional regulation. While CTCF biology has been extensively studied, how the domains of CTCF function to regulate transcription remains unknown. Additionally, the original auxin-inducible degron 1 (AID1) system has limitations in investigating the function of CTCF.

RESULTS

We employ an improved auxin-inducible degron technology, AID2, to facilitate the study of acute depletion of CTCF while overcoming the limitations of the previous AID system. As previously observed through the AID1 system and steady-state RNA analysis, the new AID2 system combined with SLAM-seq confirms that CTCF depletion leads to modest nascent and steady-state transcript changes. A CTCF domain sgRNA library screening identifies the zinc finger (ZF) domain as the region within CTCF with the most functional relevance, including ZFs 1 and 10. Removal of ZFs 1 and 10 reveals genomic regions that independently require these ZFs for DNA binding and transcriptional regulation. Notably, loci regulated by either ZF1 or ZF10 exhibit unique CTCF binding motifs specific to each ZF.

CONCLUSIONS

By extensively comparing the AID1 and AID2 systems for CTCF degradation in SEM cells, we confirm that AID2 degradation is superior for achieving miniAID-tagged protein degradation without the limitations of the AID1 system. The model we create that combines AID2 depletion of CTCF with exogenous overexpression of CTCF mutants allows us to demonstrate how peripheral ZFs intricately orchestrate transcriptional regulation in a cellular context for the first time.

CTCF: an R/bioconductor data package of human and mouse CTCF binding sites

Summary

CTCF (CCCTC-binding factor) is an 11-zinc-finger DNA binding protein which regulates much of the eukaryotic genome's 3D structure and function. The diversity of CTCF binding motifs has led to a fragmented landscape of CTCF binding data. We collected position weight matrices of CTCF binding motifs and defined strand-oriented CTCF binding sites in the human and mouse genomes, including the recent Telomere to Telomere and mm39 assemblies. We included selected experimentally determined and predicted CTCF binding sites, such as CTCF-bound cis-regulatory elements from SCREEN ENCODE. We recommend filtering strategies for CTCF binding motifs and demonstrate that liftOver is a viable alternative to convert CTCF coordinates between assemblies. Our comprehensive data resource and usage recommendations can serve to harmonize and strengthen the reproducibility of genomic studies utilizing CTCF binding data.

Availability and implementation

https://bioconductor.org/packages/CTCF. Companion website: https://dozmorovlab.github.io/CTCF/; Code to reproduce the analyses: https://github.com/dozmorovlab/CTCF.dev.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

CTCF shapes chromatin structure and gene expression in health and disease

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

3D chromatin architecture and transcription regulation in cancer

Chromatin has distinct three-dimensional (3D) architectures important in key biological processes, such as cell cycle, replication, differentiation, and transcription regulation. In turn, aberrant 3D structures play a vital role in developing abnormalities and diseases such as cancer. This review discusses key 3D chromatin structures (topologically associating domain, lamina-associated domain, and enhancer-promoter interactions) and corresponding structural protein elements mediating 3D chromatin interactions [CCCTC-binding factor, polycomb group protein, cohesin, and Brother of the Regulator of Imprinted Sites (BORIS) protein] with a highlight of their associations with cancer. We also summarise the recent development of technologies and bioinformatics approaches to study the 3D chromatin interactions in gene expression regulation, including crosslinking and proximity ligation methods in the bulk cell population (ChIA-PET and HiChIP) or single-molecule resolution (ChIA-drop), and methods other than proximity ligation, such as GAM, SPRITE, and super-resolution microscopy techniques.

Canonical WNT signaling-dependent gating of MYC requires a noncanonical CTCF function at a distal binding site

Abnormal WNT signaling increases MYC expression in colon cancer cells in part via oncogenic super-enhancer-(OSE)-mediated gating of the active MYC to the nuclear pore in a poorly understood process. We show here that the principal tenet of the WNT-regulated MYC gating, facilitating nuclear export of the MYC mRNA, is regulated by a CTCF binding site (CTCFBS) within the OSE to confer growth advantage in HCT-116 cells. To achieve this, the CTCFBS directs the WNT-dependent trafficking of the OSE to the nuclear pore from intra-nucleoplasmic positions in a stepwise manner. Once the OSE reaches a peripheral position, which is triggered by a CTCFBS-mediated CCAT1 eRNA activation, its final stretch (≤0.7 μm) to the nuclear pore requires the recruitment of AHCTF1, a key nucleoporin, to the CTCFBS. Thus, a WNT/ß-catenin-AHCTF1-CTCF-eRNA circuit enables the OSE to promote pathological cell growth by coordinating the trafficking of the active MYC gene within the 3D nuclear architecture.

Gene-gating of a MYC oncogenic super-enhancer (OSE) increases its expression in colon cancer cells in a poorly understood process. Here the authors show that MYC gating requires a CTCF binding site (CTCFBS) within the OSE that directs the stepwise trafficking of the OSE to the nuclear pore to facilitate increased nuclear export of MYC mRNA, which results in a growth advantage.

Key roles of CCCTC-binding factor in cancer evolution and development

The processes of cancer and embryonic development have a partially overlapping effect. Several transcription factor families, which are highly conserved in the evolutionary history of biology, play a key role in the development of cancer and are often responsible for the pivotal developmental processes such as cell survival, expansion, senescence, and differentiation. As an evolutionary conserved and ubiquitously expression protein, CCCTC-binding factor (CTCF) has diverse regulatory functions, including gene regulation, imprinting, insulation, X chromosome inactivation, and the establishment of three-dimensional (3D) chromatin structure during human embryogenesis. In various cancers, CTCF is considered as a tumor suppressor gene and plays homeostatic roles in maintaining genome function and integrity. However, the mechanisms of CTCF in tumor development have not been fully elucidated. Here, this review will focus on the key roles of CTCF in cancer evolution and development (Cancer Evo-Dev) and embryogenesis.

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

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

CTCF as a regulator of alternative splicing: new tricks for an old player

Three decades of research have established the CCCTC-binding factor (CTCF) as a ubiquitously expressed chromatin organizing factor and master regulator of gene expression. A new role for CTCF as a regulator of alternative splicing (AS) has now emerged. CTCF has been directly and indirectly linked to the modulation of AS at the individual transcript and at the transcriptome-wide level. The emerging role of CTCF-mediated regulation of AS involves diverse mechanisms; including transcriptional elongation, DNA methylation, chromatin architecture, histone modifications, and regulation of splicing factor expression and assembly. CTCF thereby appears to not only co-ordinate gene expression regulation but contributes to the modulation of transcriptomic complexity. In this review, we highlight previous discoveries regarding the role of CTCF in AS. In addition, we summarize detailed mechanisms by which CTCF mediates AS regulation. We propose opportunities for further research designed to examine the possible fate of CTCF-mediated alternatively spliced genes and associated biological consequences. CTCF has been widely acknowledged as the 'master weaver of the genome'. Given its multiple connections, further characterization of CTCF's emerging role in splicing regulation might extend its functional repertoire towards a 'conductor of the splicing orchestra'.

In silico Characterization of the Heme Oxygenase 1 From Bottlenose Dolphin (Tursiops truncatus): Evidence of Changes in the Active Site and Purifying Selection

Cetacea is a clade well-adapted to the aquatic lifestyle, with diverse adaptations and physiological responses, as well as a robust antioxidant defense system. Serious injuries caused by boats and fishing nets are common in bottlenose dolphins (Tursiops truncatus); however, these animals do not show signs of serious infections. Evidence suggests an adaptive response to tissue damage and associated infections in cetaceans. Heme oxygenase (HO) is a cytoprotective protein that participates in the anti-inflammatory response. HO catalyzes the first step in the oxidative degradation of the heme group. Various stimuli, including inflammatory mediators, regulate the inducible HO-1 isoform. This study aims to characterize HO-1 of the bottlenose dolphin in silico and compare its structure to the terrestrial mammal protein. Upstream HO-1 sequence of the bottlenose dolphin was obtained from NCBI and Ensemble databases, and the gene structure was determined using bioinformatics tools. Five exons and four introns were identified, and proximal regulatory elements were detected in the upstream region. The presence of 10 α-helices, three 310 helices, the heme group lodged between the proximal and distal helices, and a histidine-25 in the proximal helix serving as a ligand to the heme group were inferred for T. truncatus. Amino acid sequence alignment suggests HO-1 is a conserved protein. The HO-1 "fingerprint" and histidine-25 appear to be fully conserved among all species analyzed. Evidence of positive selection within an α-helix configuration without changes in protein configuration and evidence of purifying selection were found, indicating evolutionary conservation of the coding sequence structure.

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

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

ChIAMM: A Mixture Model for Statistical Analysis of Long-Range Chromatin Interactions From ChIA-PET Experiments

Chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) is an important experimental method for detecting specific protein-mediated chromatin loops genome-wide at high resolution. Here, we proposed a new statistical approach with a mixture model, chromatin interaction analysis using mixture model (ChIAMM), to detect significant chromatin interactions from ChIA-PET data. The statistical model is cast into a Bayesian framework to consider more systematic biases: the genomic distance, local enrichment, mappability, and GC content. Using different ChIA-PET datasets, we evaluated the performance of ChIAMM and compared it with the existing methods, including ChIA-PET Tool, ChiaSig, Mango, ChIA-PET2, and ChIAPoP. The result showed that the new approach performed better than most top existing methods in detecting significant chromatin interactions in ChIA-PET experiments.

Oplr16 serves as a novel chromatin factor to control stem cell fate by modulating pluripotency-specific chromosomal looping and TET2-mediated DNA demethylation

Formation of a pluripotency-specific chromatin network is a critical event in reprogramming somatic cells into pluripotent status. To characterize the regulatory components in this process, we used ‘chromatin RNA in situ reverse transcription sequencing’ (CRIST-seq) to profile RNA components that interact with the pluripotency master gene Oct4. Using this approach, we identified a novel nuclear lncRNA Oplr16 that was closely involved in the initiation of reprogramming. Oplr16 not only interacted with the Oct4 promoter and regulated its activity, but it was also specifically activated during reprogramming to pluripotency. Active expression of Oplr16 was required for optimal maintenance of pluripotency in embryonic stem cells. Oplr16 was also able to enhance reprogramming of fibroblasts into pluripotent cells. RNA reverse transcription-associated trap sequencing (RAT-seq) indicated that Oplr16 interacted with multiple target genes related to stem cell self-renewal. Of note, Oplr16 utilized its 3′-fragment to recruit the chromatin factor SMC1 to orchestrate pluripotency-specific intrachromosomal looping. After binding to the Oct4 promoter, Oplr16 recruited TET2 to induce DNA demethylation and activate Oct4 in fibroblasts, leading to enhanced reprogramming. These data suggest that Oplr16 may act as a pivotal chromatin factor to control stem cell fate by modulating chromatin architecture and DNA demethylation.

Signal-induced PARP1-Erk synergism mediates IEG expression

A recently disclosed Erk-induced PARP1 activation mechanism mediates the expression of immediate early genes (IEGs) in response to a variety of extra- and intracellular signals implicated in memory acquisition, development and proliferation. Here, we review this mechanism, which is initiated by stimulation-induced binding of PARP1 to phosphorylated Erk translocated into the nucleus. This binding maintains long-lasting synergistic activity of these proteins, which offers a new pattern for targeted therapy.

Advances in culture, expansion and mechanistic studies of corneal endothelial cells: a systematic review

Human corneal endothelial cells are notorious for their restricted proliferative ability in vivo and in vitro. Hence, injury or dysfunction of these cells may easily result in blindness. Currently, the only treatment is to transplant a donor cornea that contains a healthy corneal endothelium. However there is a severe global shortage of donor corneas and there remains an unmet clinical need to engineer human corneal grafts with healthy corneal endothelium. In this review, we present current advances in the culture, expansion, and molecular understandings of corneal endothelial cells in vitro in order to help establish methods of engineering human corneal endothelial grafts.

BORIS: a key regulator of cancer stemness

BORIS (CTCFL) is a DNA binding protein which is involved in tumorigenesis. Although, there are different opinions on the level of gene expression and function of BORIS in normal and cancer tissues, the results of many studies have classified BORIS as a protein belonging to cancer/testis (CT) genes, which are identified as a group of genes that are expressed normally in testis, and abnormally in various types of cancers. In testis, BORIS induces the expression of some male germ cell/testis specific genes, and plays crucial roles during spermatogenesis and production of sperm. In tumorigenesis, the role of BORIS in the expression induction of some CT genes and oncogenes, as well as increasing proliferation/viability of cancer cells has been demonstrated in many researches. In addition to cancer cells, some believe that BORIS is also expressed in normal conditions and plays a universal function in cell division and regulation of genes. The following is a comprehensive review on contradictory views on the expression pattern and biological function of BORIS in normal, as well as cancer cells/tissues, and presents some evidence that support the expression of BORIS in cancer stem cells (CSCs) and advanced stage/poorer differentiation grade of cancers. Boris is involved in the regulation of CSC cellular and molecular features such as self-renewal, chemo-resistance, tumorigenicity, sphere-forming ability, and migration capacity. Finally, the role of BORIS in regulating two important signaling pathways including Wnt/β-catenin and Notch in CSCs, and its ability in recruiting transcription factors or chromatin-remodeling proteins to induce tumorigenesis is discussed.

Poly(ADP-ribosyl)ation associated changes in CTCF-chromatin binding and gene expression in breast cells

CTCF is an evolutionarily conserved and ubiquitously expressed architectural protein regulating a plethora of cellular functions via different molecular mechanisms. CTCF can undergo a number of post-translational modifications which change its properties and functions. One such modifications linked to cancer is poly(ADP-ribosyl)ation (PARylation). The highly PARylated CTCF form has an apparent molecular mass of 180 kDa (referred to as CTCF180), which can be distinguished from hypo- and non-PARylated CTCF with the apparent molecular mass of 130 kDa (referred to as CTCF130). The existing data accumulated so far have been mainly related to CTCF130. However, the properties of CTCF180 are not well understood despite its abundance in a number of primary tissues. In this study we performed ChIP-seq and RNA-seq analyses in human breast cells 226LDM, which display predominantly CTCF130 when proliferating, but CTCF180 upon cell cycle arrest. We observed that in the arrested cells the majority of sites lost CTCF, whereas fewer sites gained CTCF or remain bound (i.e. common sites). The classical CTCF binding motif was found in the lost and common, but not in the gained sites. The changes in CTCF occupancies in the lost and common sites were associated with increased chromatin densities and altered expression from the neighboring genes. Based on these results we propose a model integrating the CTCF130/180 transition with CTCF-DNA binding and gene expression changes. This study also issues an important cautionary note concerning the design and interpretation of any experiments using cells and tissues where CTCF180 may be present.

Enhancer functions in three dimensions: beyond the flat world perspective

Transcriptional enhancers constitute a subclass of regulatory elements that facilitate transcription. Such regions are generally organized by short stretches of DNA enriched in transcription factor-binding sites but also can include very large regions containing clusters of enhancers, termed super-enhancers. These regions increase the probability or the rate (or both) of transcription generally in cis and sometimes over very long distances by altering chromatin states and the activity of Pol II machinery at promoters. Although enhancers were discovered almost four decades ago, their inner workings remain enigmatic. One important opening into the underlying principle has been provided by observations that enhancers make physical contacts with their target promoters to facilitate the loading of the RNA polymerase complex. However, very little is known about how such chromatin loops are regulated and how they govern transcription in the three-dimensional context of the nuclear architecture. Here, we present current themes of how enhancers may boost gene expression in three dimensions and we identify currently unresolved key questions.

Novel CTCF mutations in Chinese patients with ovarian endometriosis

Endometriosis is a common gynecological disease characterized by the outgrowth of the endometrium, however, the detailed molecular etiology remains largely uncharacterized. Recent studies have implicated that endometriosis is potentially a precancerous lesion, and that CCCTC‑binding factor (CTCF) mutations may be involved in the pathogenesis of this disorder. However, the detailed CTCF mutation spectrum in Chinese patients with ovarian endometriosis remains largely unknown. In the present study, a cohort of 92 patients with ovarian endometriosis were analyzed for the presence of CTCF mutations by sequencing the entire coding regions. In addition, 67 healthy eutopic endometrial tissues and 46 healthy ovarian tissues from control samples (without endometriosis) were also analyzed. In total, two CTCF missense mutations, p.K206E (c.616A>G) and p.H373L (c.1118A>T), were identified in 2/92 (2.2%) endometriotic lesions. The patient with the p.K206E mutation was 26 years old and diagnosed with primary infertility, whereas the patient with the p.H373L mutation was 37 years old and concurrently diagnosed with uterine leiomyoma. The p.H373L mutation was previously identified in endometrial cancer samples with low frequency, while the p.K206E mutation was novel. In addition, no CTCF mutations were detected in the 67 healthy eutopic endometrial and 46 healthy ovarian tissue samples. In silico prediction and evolutionary conservation analysis suggested that these CTCF mutations may be pathogenic. In summary, the present study identified 2 potential pathogenic CTCF mutations in endometriotic lesions from 2/92 patients with ovarian endometriosis. These results, together with a prior exome‑sequencing based study, suggest that CTCF mutations may be involved in the development of ovarian endometriosis.

Deciphering the DNA code for the function of the Drosophila polydactyl zinc finger protein Suppressor of Hairy-wing

Polydactyl zinc finger (ZF) proteins have prominent roles in gene regulation and often execute multiple regulatory functions. To understand how these proteins perform varied regulation, we studiedDrosophila Suppressor of Hairy-wing [Su(Hw)], an exemplar multifunctional polydactyl ZF protein. We identified separation-of-function (SOF) alleles that encode proteins disrupted in a single ZF that retain one of the Su(Hw) regulatory activities. Through extended in vitro analyses of the Su(Hw) ZF domain, we show that clusters of ZFs bind individual modules within a compound DNA consensus sequence. Through in vivo analysis of SOF mutants, we find that Su(Hw) genomic sites separate into sequence subclasses comprised of combinations of modules, with subclasses enriched for different chromatin features. These data suggest a Su(Hw) code, wherein DNA binding dictates its cofactor recruitment and regulatory output. We propose that similar DNA codes might be used to confer multiple regulatory functions of other polydactyl ZF proteins.

Senescence Mediated by p16INK4a Impedes Reprogramming of Human Corneal Endothelial Cells into Neural Crest Progenitors

Human corneal endothelial cells (HCECs) have limited proliferative capacity due to “contact-inhibition” at G1 phase. Such contact-inhibition can be delayed from Day 21 to Day 42 by switching EGF-containing SHEM to LIF/bFGF-containing MESCM through transient activation of LIF-JAK1-STAT3 signaling that delays eventual nuclear translocation of p16^INK4a. Using the latter system, we have reported a novel tissue engineering technique by implementing 5 weekly knockdowns with p120 catenin (p120) and Kaiso siRNAs since Day 7 to achieve effective expansion of HCEC monolayers to a transplantable size with a normal HCEC density, through reprogramming of HCECs into neural crest progenitors by activating p120-Kaiso-RhoA-ROCK-canonical BMP signaling. Herein, we noted that a single knockdown with p120-Kaiso siRNAs at Day 42 failed to achieve such reprogramming when contact inhibition transitioned to senescence with nuclear translocation of p16^INK4a. In contrast, 5 weekly knockdowns with p120-Kaiso siRNAs since Day 7 precluded senescence mediated by p16^INK4a by inducing nuclear translocation of Bmi1 because of sustained activation of JAK2-STAT3 signaling downstream of p120-Kaiso-RhoA-ROCK signaling. STAT3 or Bmi1 siRNA impeded nuclear exclusion of p16^INK4a and suppressed the reprogramming induced by p120-Kaiso siRNAs, suggesting that another important engineering strategy of HCEC lies in prevention of senescence mediated by nuclear translocation of p16^INK4a.

Motif oriented high-resolution analysis of ChIP-seq data reveals the topological order of CTCF and cohesin proteins on DNA

BACKGROUND

ChIP-seq provides a wealth of information on the approximate location of DNA-binding proteins genome-wide. It is known that the targeted motifs in most cases can be found at the peak centers. A high resolution mapping of ChIP-seq peaks could in principle allow the fine mapping of the protein constituents within protein complexes, but the current ChIP-seq analysis pipelines do not target the basepair resolution strand specific mapping of peak summits.

RESULTS

The approach proposed here is based on i) locating regions that are bound by a sufficient number of proteins constituting a complex; ii) determining the position of the underlying motif using either a direct or a de novo motif search approach; and iii) determining the exact location of the peak summits with respect to the binding motif in a strand specific manner. We applied this method for analyzing the CTCF/cohesin complex, which holds together DNA loops. The relative positions of the constituents of the complex were determined with one-basepair estimated accuracy. Mapping the positions on a 3D model of DNA made it possible to deduce the approximate local topology of the complex that allowed us to predict how the CTCF/cohesin complex locks the DNA loops. As the positioning of the proteins was not compatible with previous models of loop closure, we proposed a plausible "double embrace" model in which the DNA loop is held together by two adjacent cohesin rings in such a way that the ring anchored by CTCF to one DNA duplex encircles the other DNA double helix and vice versa.

CONCLUSIONS

A motif-centered, strand specific analysis of ChIP-seq data improves the accuracy of determining peak positions. If a genome contains a large number of binding sites for a given protein complex, such as transcription factor heterodimers or transcription factor/cofactor complexes, the relative position of the constituent proteins on the DNA can be established with an accuracy that allow one to deduce the local topology of the protein complex. The proposed high resolution mapping approach of ChIP-seq data is applicable for detecting the contact topology of DNA-binding protein complexes.

Functional Dissection of the Blocking and Bypass Activities of the Fab-8 Boundary in the Drosophila Bithorax Complex

Functionally autonomous regulatory domains direct the parasegment-specific expression of the Drosophila Bithorax complex (BX-C) homeotic genes. Autonomy is conferred by boundary/insulator elements that separate each regulatory domain from its neighbors. For six of the nine parasegment (PS) regulatory domains in the complex, at least one boundary is located between the domain and its target homeotic gene. Consequently, BX-C boundaries must not only block adventitious interactions between neighboring regulatory domains, but also be permissive (bypass) for regulatory interactions between the domains and their gene targets. To elucidate how the BX-C boundaries combine these two contradictory activities, we have used a boundary replacement strategy. We show that a 337 bp fragment spanning the Fab-8 boundary nuclease hypersensitive site and lacking all but 83 bp of the 625 bp Fab-8 PTS (promoter targeting sequence) fully rescues a Fab-7 deletion. It blocks crosstalk between the iab-6 and iab-7 regulatory domains, and has bypass activity that enables the two downstream domains, iab-5 and iab-6, to regulate Abdominal-B (Abd-B) transcription in spite of two intervening boundary elements. Fab-8 has two dCTCF sites and we show that they are necessary both for blocking and bypass activity. However, CTCF sites on their own are not sufficient for bypass. While multimerized dCTCF (or Su(Hw)) sites have blocking activity, they fail to support bypass. Moreover, this bypass defect is not rescued by the full length PTS. Finally, we show that orientation is critical for the proper functioning the Fab-8 replacement. Though the inverted Fab-8 boundary still blocks crosstalk, it disrupts the topology of the Abd-B regulatory domains and does not support bypass. Importantly, altering the orientation of the Fab-8 dCTCF sites is not sufficient to disrupt bypass, indicating that orientation dependence is conferred by other factors.

Boundary elements in the Bithorax complex have two seemingly contradictory activities. They must block crosstalk between neighboring regulatory domains, but at the same time be permissive (insulator bypass) for regulatory interactions between the domains and the BX-C homeotic genes. We have used a replacement strategy to investigate how they carry out these two functions. We show that a 337 bp fragment spanning the Fab-8 boundary nuclease hypersensitive site is sufficient to fully rescue a Fab-7 boundary deletion. It blocks crosstalk and supports bypass. As has been observed in transgene assays, blocking activity requires the Fab-8 dCTCF sites, while full bypass activity requires the dCTCF sites plus a small part of PTS. In transgene assays, bypass activity typically depends on the orientation of the two insulators relative to each other. A similar orientation dependence is observed for the Fab-8 replacement in BX-C. When the orientation of the Fab-8 boundary is reversed, bypass activity is lost, while blocking is unaffected. Interestingly, unlike what has been observed in mammals, reversing the orientation of only the Fab-8 dCTCF sites does not affect boundary function. This finding indicates that other Fab-8 factors must play a critical role in determining orientation. Taken together, our findings argue that carrying out the paradoxical functions of the BX-C boundaries does not require any unusual or special properties; rather BX-C boundaries utilize generic blocking and insulator bypass activities that are appropriately adapted to their regulatory context. Thus making them a good model for studying the functional properties of boundaries/insulators in their native setting.

A PARP1-ERK2 synergism is required for the induction of LTP

Unexpectedly, a post-translational modification of DNA-binding proteins, initiating the cell response to single-strand DNA damage, was also required for long-term memory acquisition in a variety of learning paradigms. Our findings disclose a molecular mechanism based on PARP1-Erk synergism, which may underlie this phenomenon. A stimulation induced PARP1 binding to phosphorylated Erk2 in the chromatin of cerebral neurons caused Erk-induced PARP1 activation, rendering transcription factors and promoters of immediate early genes (IEG) accessible to PARP1-bound phosphorylated Erk2. Thus, Erk-induced PARP1 activation mediated IEG expression implicated in long-term memory. PARP1 inhibition, silencing, or genetic deletion abrogated stimulation-induced Erk-recruitment to IEG promoters, gene expression and LTP generation in hippocampal CA3-CA1-connections. Moreover, a predominant binding of PARP1 to single-strand DNA breaks, occluding its Erk binding sites, suppressed IEG expression and prevented the generation of LTP. These findings outline a PARP1-dependent mechanism required for LTP generation, which may be implicated in long-term memory acquisition and in its deterioration in senescence.

A PARP1-Erk2 synergism is required for stimulation-induced expression of immediate early genes

A PARP1-Erk2 synergism was required to generate synaptic long-term potentiation in the CA3-CA1 hippocampal connections. This molecular mechanism was associated with the recently identified pivotal role of polyADP-ribosylation in learning. High frequency electrical stimulation of cortical and hippocampal neurons induced binding of phosphorylated Erk2 (transported into the nucleus) to the nuclear protein PARP1. PARP1-Erk2 binding induced PARP1 activation and polyADP-ribosylation of its prominent substrate, linker histone H1. A facilitated access of PARP1-bound phosphorylated Erk2 to its substrates, transcription factors Elk1 and CREB was attributed to the release of polyADP-ribosylated H1 from the DNA, causing local DNA relaxation. Erk-induced phosphorylation of transcription factors activating the HAT activity of CBP (CREB binding protein), recruited acetylated histone H4 to the promoters of immediate early genes (IEG) cfos, zif268 and arc, which are implicated in synaptic plasticity. In accordance, their induced expression was suppressed after PARP1 genetic deletion in PARP1-KO mice, or after PARP1 inhibition or silencing. Moreover, under these conditions, long-term synaptic potentiation (LTP) (indicating synaptic plasticity) was not generation in the hippocampal CA3-CA1 connections, and learning abilities were impaired. Furthermore, both IEG expression and LTP generation failed when cerebral neurons accumulated single strand DNA breaks, due to a predominant binding of PARP1 to nicked DNA, occluding its Erk binding sites. Thus, a declined synaptic plasticity is anticipated when aged cerebral neurons accumulate DNA single-strand breaks during life span.

Evidence for a common evolutionary rate in metazoan transcriptional networks

Genome sequences diverge more rapidly in mammals than in other animal lineages, such as birds or insects. However, the effect of this rapid divergence on transcriptional evolution remains unclear. Recent reports have indicated a faster divergence of transcription factor binding in mammals than in insects, but others found the reverse for mRNA expression. Here, we show that these conflicting interpretations resulted from differing methodologies. We performed an integrated analysis of transcriptional network evolution by examining mRNA expression, transcription factor binding and cis-regulatory motifs across >25 animal species, including mammals, birds and insects. Strikingly, we found that transcriptional networks evolve at a common rate across the three animal lineages. Furthermore, differences in rates of genome divergence were greatly reduced when restricting comparisons to chromatin-accessible sequences. The evolution of transcription is thus decoupled from the global rate of genome sequence evolution, suggesting that a small fraction of the genome regulates transcription.

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

The genetic information that makes each individual unique is encoded in DNA molecules. Cells read this molecular instruction manual by a process called transcription, in which proteins called transcription factors bind to DNA in specific places and regulate which sections of the DNA will be expressed. These 'transcripts' are active molecules that determine the cell’s – and ultimately the individual’s – characteristics. However, it is not well understood how alterations in the DNA of different individuals or species can lead to changes in where the transcription factors bind, and in which transcripts are expressed.

Carvunis, Wang, Skola et al. set out to determine if there is a relationship between how often DNA changes and how often transcription changes during the evolution of animals. The experiments examined the abundance of transcripts in the cells of a variety of animal species with close or distant evolutionary relationships. For example, the house mouse was compared to a close relative called the Algerian mouse, to another species of rodent (rat) and to humans.

The experiments show that the changes in transcript abundances are happening at similar rates in mammals, birds and insects, even though DNA changes at very different rates in these groups of animals. This similarity was also observed for other aspects of transcription, such as in changes to where transcription factors bind to DNA. The next challenges are to find out what makes transcription evolve at such similar rates in these groups of animals, and whether these findings extend to other species and to other processes in cells.

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

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.

Copy number variations alter methylation and parallel IGF2 overexpression in adrenal tumors

Overexpression of insulin growth factor 2 (IGF2) is a hallmark of adrenocortical carcinomas and pheochromocytomas. Previous studies investigating the IGF2/H19 locus have mainly focused on a single molecular level such as genomic alterations or altered DNA methylation levels and the causal changes underlying IGF2 overexpression are still not fully established. In the current study, we analyzed 62 tumors of the adrenal gland from patients with Conn's adenoma (CA, n=12), pheochromocytomas (PCC, n=10), adrenocortical benign tumors (ACBT, n=20), and adrenocortical carcinomas (ACC, n=20). Gene expression, somatic copy number variation of chr11p15.5, and DNA methylation status of three differential methylated regions of the IGF2/H19 locus including the H19 imprinting control region were integratively analyzed. IGF2 overexpression was found in 85% of the ACCs and 100% of the PCCs compared to 23% observed in CAs and ACBTs. Copy number aberrations of chr11p15.5 were abundant in both PCCs and ACCs but while PCCs retained a diploid state, ACCs were frequently tetraploid (7/19). Loss of either a single allele or loss of two alleles of the same parental origin in tetraploid samples resulted in a uniparental disomy-like genotype. These copy number changes correlated with hypermethylation of the H19 ICR suggesting that the lost alleles were the unmethylated maternal alleles. Our data provide conclusive evidence that loss of the maternal allele correlates with IGF2 overexpression in adrenal tumors and that hypermethylation of the H19 ICR is a consequence thereof.

CTCF regulates NELF, DSIF and P-TEFb recruitment during transcription

CTCF is a versatile transcription factor with well-established roles in chromatin organization and insulator function. Recent findings also implicate CTCF in the control of elongation by RNA polymerase (RNAP) II. Here we show that CTCF knockdown abrogates RNAP II pausing at the early elongation checkpoint of c-myc by affecting recruitment of DRB-sensitivity-inducing factor (DSIF). CTCF knockdown also causes a termination defect on the U2 snRNA genes (U2), by affecting recruitment of negative elongation factor (NELF). In addition, CTCF is required for recruitment of positive elongation factor b (P-TEFb), which phosphorylates NELF, DSIF, and Ser2 of the RNAP II CTD to activate elongation of transcription of c-myc and recognition of the snRNA gene-specific 3' box RNA processing signal. These findings implicate CTCF in a complex network of protein:protein/protein:DNA interactions and assign a key role to CTCF in controlling RNAP II transcription through the elongation checkpoint of the protein-coding c-myc and the termination site of the non-coding U2, by regulating the recruitment and/or activity of key players in these processes.

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.

Functional role of dimerization and CP190 interacting domains of CTCF protein in Drosophila melanogaster

Background

Insulators play a central role in gene regulation, chromosomal architecture and genome function in higher eukaryotes. To learn more about how insulators carry out their diverse functions, we have begun an analysis of the Drosophila CTCF (dCTCF). CTCF is one of the few insulator proteins known to be conserved from flies to man.

Results

In the studies reported here we have focused on the identification and characterization of two dCTCF protein interaction modules. The first mediates dCTCF multimerization, while the second mediates dCTCF–CP190 interactions. The multimerization domain maps in the N-terminus of the dCTCF protein and likely mediates the formation of tetrameric complexes. The CP190 interaction module encompasses a sequence ~200 amino acids long that spans the C-terminal and mediates interactions with the N-terminal BTB domain of the CP190 protein. Transgene rescue experiments showed that a dCTCF protein lacking sequences critical for CP190 interactions was almost as effective as wild type in rescuing the phenotypic effects of a dCTCF null allele. The mutation did, however, affect CP190 recruitment to specific Drosophila insulator elements and had a modest effect on dCTCF chromatin association. A protein lacking the N-terminal dCTCF multimerization domain incompletely rescued the zygotic and maternal effect lethality of the null and did not rescue the defects in Abd-B regulation evident in surviving adult dCTCF mutant flies. Finally, we show that elimination of maternally contributed dCTCF at the onset of embryogenesis has quite different effects on development and Abd-B regulation than is observed when the homozygous mutant animals develop in the presence of maternally derived dCTCF activity.

Conclusions

Our results indicate that dCTCF–CP190 interactions are less critical for the in vivo functions of the dCTCF protein than the N-terminal dCTCF–dCTCF interaction domain. We also show that the phenotypic consequences of dCTCF mutations differ depending upon when and how dCTCF activity is lost.

Electronic supplementary material

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

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.

Structural organization of human replication timing domains

Recent analysis of genome-wide epigenetic modification data, mean replication timing (MRT) profiles and chromosome conformation data in mammals have provided increasing evidence that flexibility in replication origin usage is regulated locally by the epigenetic landscape and over larger genomic distances by the 3D chromatin architecture. Here, we review the recent results establishing some link between replication domains and chromatin structural domains in pluripotent and various differentiated cell types in human. We reconcile the originally proposed dichotomic picture of early and late constant timing regions that replicate by multiple rather synchronous origins in separated nuclear compartments of open and closed chromatins, with the U-shaped MRT domains bordered by "master" replication origins specified by a localized (∼200-300 kb) zone of open and transcriptionally active chromatin from which a replication wave likely initiates and propagates toward the domain center via a cascade of origin firing. We discuss the relationships between these MRT domains, topologically associated domains and lamina-associated domains. This review sheds a new light on the epigenetically regulated global chromatin reorganization that underlies the loss of pluripotency and the determination of differentiation properties.

Embryonic stem cell specific "master" replication origins at the heart of the loss of pluripotency

Epigenetic regulation of the replication program during mammalian cell differentiation remains poorly understood. We performed an integrative analysis of eleven genome-wide epigenetic profiles at 100 kb resolution of Mean Replication Timing (MRT) data in six human cell lines. Compared to the organization in four chromatin states shared by the five somatic cell lines, embryonic stem cell (ESC) line H1 displays (i) a gene-poor but highly dynamic chromatin state (EC4) associated to histone variant H2AZ rather than a HP1-associated heterochromatin state (C4) and (ii) a mid-S accessible chromatin state with bivalent gene marks instead of a polycomb-repressed heterochromatin state. Plastic MRT regions (≲ 20% of the genome) are predominantly localized at the borders of U-shaped timing domains. Whereas somatic-specific U-domain borders are gene-dense GC-rich regions, 31.6% of H1-specific U-domain borders are early EC4 regions enriched in pluripotency transcription factors NANOG and OCT4 despite being GC poor and gene deserts. Silencing of these ESC-specific “master” replication initiation zones during differentiation corresponds to a loss of H2AZ and an enrichment in H3K9me3 mark characteristic of late replicating C4 heterochromatin. These results shed a new light on the epigenetically regulated global chromatin reorganization that underlies the loss of pluripotency and lineage commitment.

During development, embryonic stem cell (ESC) enter a program of cell differentiation eventually leading to all the necessary differentiated cell types. Understanding the mechanisms responsible for the underlying modifications of the gene expression program is of fundamental importance, as it will likely have strong impact on the development of regenerative medicine. We show that besides some epigenetic regulation, ubiquitous master replication origins at replication timing U-domain borders shared by 6 human cell types are transcriptionally active open chromatin regions specified by a local enrichment in nucleosome free regions encoded in the DNA sequence suggesting that they have been selected during evolution. In contrast, ESC specific master replication origins bear a unique epigenetic signature (enrichment in CTCF, H2AZ, NANOG, OCT4, …) likely contributing to maintain ESC chromatin in a highly dynamic and accessible state that is refractory to polycomb and HP1 heterochromatin spreading. These ESC specific master origins thus appear as key genomic regions where epigenetic control of chromatin organization is at play to maintain pluripotency of stem cell lineages and to guide lineage commitment to somatic cell types.

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.

Experimental strategies to manipulate the cellular levels of the multifunctional factor CTCF

Cellular homeostasis is the result of an intricate and coordinated combinatorial of biochemical and molecular processes. Among them is the control of gene expression in the context of the chromatin structure which is central for cell survival. Interdependent action of transcription factors, cofactors, chromatin remodeling activities, and three-dimensional organization of the genome are responsible to reach exquisite levels of gene expression. Among such transcription factors there is a subset of highly specialized nuclear factors with features resembling master regulators with a large variety of functions. This is turning to be the case of the multifunctional nuclear factor CCCTC-binding protein (CTCF) which is involved in gene regulation, chromatin organization, and three-dimensional conformation of the genome inside the cell nucleus. Technically its study has turned to be challenging, in particular its posttranscriptional interference by small interference RNAs. Here we describe three main strategies to downregulate the overall abundance of CTCF in culture cell lines.

Chromatin associations in Arabidopsis interphase nuclei

The arrangement of chromatin within interphase nuclei seems to be caused by topological constraints and related to gene expression depending on tissue and developmental stage. In yeast and animals it was found that homologous and heterologous chromatin association are required to realize faithful expression and DNA repair. To test whether such associations are present in plants we analyzed Arabidopsis thaliana interphase nuclei by FISH using probes from different chromosomes. We found that chromatin fiber movement and variable associations, although in general relatively seldom, may occur between euchromatin segments along chromosomes, sometimes even over large distances. The combination of euchromatin segments bearing high or low co-expressing genes did not reveal different association frequencies probably due to adjacent genes of deviating expression patterns. Based on previous data and on FISH analyses presented here, we conclude that the global interphase chromatin organization in A. thaliana is relatively stable, due to the location of its 10 centromeres at the nuclear periphery and of the telomeres mainly at the centrally localized nucleolus. Nevertheless, chromatin movement enables a flexible spatial genome arrangement in plant nuclei.

Expression analysis of BORIS during pluripotent, differentiated, cancerous, and non-cancerous cell states

BORIS/CTCFL is an 11 zinc finger protein, which is the paralog of CTCF, a ubiquitously expressed protein with diverse roles in gene expression and chromatin organization. Several studies have shown that the expression of BORIS is restricted to normal adult testis, pluripotent cells, and diverse cancer cell lines. Thus, it is known as a cancer-testis (CT) gene that has been hypothesized to exhibit oncogenic properties and to be involved in cancer cell proliferation. On the contrary, other reports have shown that its expression is more widespread and can be detected in differentiated and normal somatic cells; hence, it might have roles in general cellular functions. The present study was aimed to analyze the expression of BORIS in different cell states of pluripotent, differentiated, cancerous and non-cancerous.We found that the two cell states of pluripotency and differentiation are not accompanied with significant variations of BORIS expression. Furthermore, Boris transcripts were detected at approximately the same level in cancer and non-cancer cell lines. These findings suggest that, in contrast to some previous reports, the expression of mouse BORIS is not limited to only cancerous cells or pluripotent cell states.

The epigenetic factor BORIS/CTCFL regulates the NOTCH3 gene expression in cancer cells

Aberrant upregulation of NOTCH3 gene plays a critical role in cancer pathogenesis. However, the underlying mechanisms are still unknown. We tested here the hypothesis that aberrant epigenetic modifications in the NOTCH3 promoter region might account for its upregulation in cancer cells. We compared DNA and histone methylation status of NOTCH3 promoter region in human normal blood cells and T cell acute lymphoblastic leukemia (T-ALL) cell lines, differentially expressing NOTCH3. We found that histone methylation, rather than DNA hypomethylation, contributes towards establishing an active chromatin status of NOTCH3 promoter in NOTCH3 overexpressing cancer cells. We discovered that the chromatin regulator protein BORIS/CTCFL plays an important role in regulating NOTCH3 gene expression. We observed that BORIS is present in T-ALL cell lines as well as in cell lines derived from several solid tumors overexpressing NOTCH3. Moreover, BORIS targets NOTCH3 promoter in cancer cells and it is able to induce and to maintain a permissive/active chromatin conformation. Importantly, the association between NOTCH3 overexpression and BORIS presence was confirmed in primary T-ALL samples from patients at the onset of the disease. Overall, our results provide novel insights into the determinants of NOTCH3 overexpression in cancer cells, by revealing a key role for BORIS as the main mediator of transcriptional deregulation of NOTCH3.

Subtelomeric CTCF and cohesin binding site organization using improved subtelomere assemblies and a novel annotation pipeline

Mapping genome-wide data to human subtelomeres has been problematic due to the incomplete assembly and challenges of low-copy repetitive DNA elements. Here, we provide updated human subtelomere sequence assemblies that were extended by filling telomere-adjacent gaps using clone-based resources. A bioinformatic pipeline incorporating multiread mapping for annotation of the updated assemblies using short-read data sets was developed and implemented. Annotation of subtelomeric sequence features as well as mapping of CTCF and cohesin binding sites using ChIP-seq data sets from multiple human cell types confirmed that CTCF and cohesin bind within 3 kb of the start of terminal repeat tracts at many, but not all, subtelomeres. CTCF and cohesin co-occupancy were also enriched near internal telomere-like sequence (ITS) islands and the nonterminal boundaries of subtelomere repeat elements (SREs) in transformed lymphoblastoid cell lines (LCLs) and human embryonic stem cell (ES) lines, but were not significantly enriched in the primary fibroblast IMR90 cell line. Subtelomeric CTCF and cohesin sites predicted by ChIP-seq using our bioinformatics pipeline (but not predicted when only uniquely mapping reads were considered) were consistently validated by ChIP-qPCR. The colocalized CTCF and cohesin sites in SRE regions are candidates for mediating long-range chromatin interactions in the transcript-rich SRE region. A public browser for the integrated display of short-read sequence-based annotations relative to key subtelomere features such as the start of each terminal repeat tract, SRE identity and organization, and subtelomeric gene models was established.

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.

A novel mechanism for CTCF in the epigenetic regulation of Bax in breast cancer cells

We previously reported the association of elevated levels of the multifunctional transcription factor, CCCTC binding factor (CTCF), in breast cancer cells with the specific anti-apoptotic function of CTCF. To understand the molecular mechanisms of this phenomenon, we investigated regulation of the human Bax gene by CTCF in breast and non-breast cells. Two CTCF binding sites (CTSs) within the Bax promoter were identified. In all cells, breast and non-breast, active histone modifications were present at these CTSs, DNA harboring this region was unmethylated, and levels of Bax mRNA and protein were similar. Nevertheless, up-regulation of Bax mRNA and protein and apoptotic cell death were observed only in breast cancer cells depleted of CTCF. We proposed that increased CTCF binding to the Bax promoter in breast cancer cells, by comparison with non-breast cells, may be mechanistically linked to the specific apoptotic phenotype in CTCF-depleted breast cancer cells. In this study, we show that CTCF binding was enriched at the Bax CTSs in breast cancer cells and tumors; in contrast, binding of other transcription factors (SP1, WT1, EGR1, and c-Myc) was generally increased in non-breast cells and normal breast tissues. Our findings suggest a novel mechanism for CTCF in the epigenetic regulation of Bax in breast cancer cells, whereby elevated levels of CTCF support preferential binding of CTCF to the Bax CTSs. In this context, CTCF functions as a transcriptional repressor counteracting influences of positive regulatory factors; depletion of breast cancer cells from CTCF therefore results in the activation of Bax and apoptosis.