Functional analysis of CTCF during mammalian limb development

Conserved Cis-Acting Range Extender Element Mediates Extreme Long-Range Enhancer Activity in Mammals

While most mammalian enhancers regulate their cognate promoters over moderate distances of tens of kilobases (kb), some enhancers act over distances in the megabase range. The sequence features enabling such extreme-distance enhancer-promoter interactions remain elusive. Here, we used in vivo enhancer replacement experiments in mice to show that short- and medium-range enhancers cannot initiate gene expression at extreme-distance range. We uncover a novel conserved cis-acting element, Range EXtender (REX), that confers extreme-distance regulatory activity and is located next to a long-range enhancer of Sall1. The REX element itself has no endogenous enhancer activity. However, addition of the REX to other short- and mid-range enhancers substantially increases their genomic interaction range. In the most extreme example observed, addition of the REX increased the range of an enhancer by an order of magnitude, from its native 71kb to 840kb. The REX element contains highly conserved [C/T]AATTA homeodomain motifs. These motifs are enriched around long-range limb enhancers genome-wide, including the ZRS, a benchmark long-range limb enhancer of Shh. Mutating the [C/T]AATTA motifs within the ZRS does not affect its limb-specific enhancer activity at short range, but selectively abolishes its long-range activity, resulting in severe limb reduction in knock-in mice. In summary, we identify a sequence signature globally associated with long-range enhancer-promoter interactions and describe a prototypical REX element that is necessary and sufficient to confer extreme-distance gene activation by remote enhancers.

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

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

Enhancer-promoter interactions can form independently of genomic distance and be functional across TAD boundaries

Topologically Associating Domains (TADs) have been suggested to facilitate and constrain enhancer-promoter interactions. However, the role of TAD boundaries in effectively restricting these interactions remains unclear. Here, we show that a significant proportion of enhancer-promoter interactions are established across TAD boundaries in Drosophila embryos, but that developmental genes are strikingly enriched in intra- but not inter-TAD interactions. We pursued this observation using the twist locus, a master regulator of mesoderm development, and systematically relocated one of its enhancers to various genomic locations. While this developmental gene can establish inter-TAD interactions with its enhancer, the functionality of these interactions remains limited, highlighting the existence of topological constraints. Furthermore, contrary to intra-TAD interactions, the formation of inter-TAD enhancer-promoter interactions is not solely driven by genomic distance, with distal interactions sometimes favored over proximal ones. These observations suggest that other general mechanisms must exist to establish and maintain specific enhancer-promoter interactions across large distances.

Establishing and maintaining Hox profiles during spinal cord development

The chromosomally-arrayed Hox gene family plays central roles in embryonic patterning and the specification of cell identities throughout the animal kingdom. In vertebrates, the relatively large number of Hox genes and pervasive expression throughout the body has hindered understanding of their biological roles during differentiation. Studies on the subtype diversification of spinal motor neurons (MNs) have provided a tractable system to explore the function of Hox genes during differentiation, and have provided an entry point to explore how neuronal fate determinants contribute to motor circuit assembly. Recent work, using both in vitro and in vivo models of MN subtype differentiation, have revealed how patterning morphogens and regulation of chromatin structure determine cell-type specific programs of gene expression. These studies have not only shed light on basic mechanisms of rostrocaudal patterning in vertebrates, but also have illuminated mechanistic principles of gene regulation that likely operate in the development and maintenance of terminal fates in other systems.

CTCF deletion alters the pluripotency and DNA methylation profile of human iPSCs

Pluripotent stem cells are characterized by their differentiation potential toward endoderm, mesoderm, and ectoderm. However, it is still largely unclear how these cell-fate decisions are mediated by epigenetic mechanisms. In this study, we explored the relevance of CCCTC-binding factor (CTCF), a zinc finger-containing DNA-binding protein, which mediates long-range chromatin organization, for directed cell-fate determination. We generated human induced pluripotent stem cell (iPSC) lines with deletions in the protein-coding region in exon 3 of CTCF, resulting in shorter transcripts and overall reduced protein expression. Chromatin immunoprecipitation showed a considerable loss of CTCF binding to target sites. The CTCF deletions resulted in slower growth and modest global changes in gene expression, with downregulation of a subset of pluripotency-associated genes and neuroectodermal genes. CTCF deletion also evoked DNA methylation changes, which were moderately associated with differential gene expression. Notably, CTCF-deletions lead to upregulation of endo-mesodermal associated marker genes and epigenetic signatures, whereas ectodermal differentiation was defective. These results indicate that CTCF plays an important role in the maintenance of pluripotency and differentiation, especially towards ectodermal lineages.

CTCF coordinates cell fate specification via orchestrating regulatory hubs with pioneer transcription factors

CCCTC-binding factor (CTCF), a ubiquitously expressed architectural protein, has emerged as a key regulator of cell identity gene transcription. However, the precise molecular mechanism underlying specialized functions of CTCF remains elusive. Here, we investigate the mechanism through integrative analyses of primary hepatocytes, myocytes, and B cells from mouse and human. We demonstrate that CTCF cooperates with lineage-specific pioneer transcription factors (TFs), including MyoD, FOXA, and PU.1, to control cell identity at 1D and 3D levels. At the 1D level, pioneer TFs facilitate lineage-specific CTCF occupancy via opening chromatin. At the 3D level, CTCF and pioneer TFs form regulatory hubs to govern the expression of cell identity genes. This mechanism is validated using MyoD-null mice, CTCF knockout mice, and CRISPR editing during myogenic differentiation. Collectively, these findings uncover a general mechanism whereby CTCF acts as a cell identity cofactor to control cell identity genes via orchestrating regulatory hubs with pioneer TFs.

Hnrnpk is essential for embryonic limb bud development as a transcription activator and a collaborator of insulator protein Ctcf

Proper development of the limb bud relies on the concordance of various signals, but its molecular mechanisms have not yet been fully illustrated. Here we report that heterogeneous nuclear ribonucleoprotein K (hnRNPK) is essential for limb bud development. Its ablation in the limb bud results in limbless forelimbs and severe deformities of the hindlimbs. In terms of mechanism, hnRNPK functions as a transcription activator for the vital genes involved in the three regulatory axes of limb bud development. Simultaneously, for the first time we elucidate that hnRNPK binds to and coordinates with the insulator protein CCCTC binding factor (CTCF) to maintain a three-dimensional chromatin architecture. Ablation of hnRNPK weakens the binding strength of CTCF to topologically associating domain (TAD) boundaries, then leading to the loose TADs, and decreased interactions between promoters and enhancers, and further decreased transcription of developmental genes. Our study establishes a fundamental and novel role of hnRNPK in regulating limb bud development.

Three-dimensional genome structure and function

Linear DNA undergoes a series of compression and folding events, forming various three-dimensional (3D) structural units in mammalian cells, including chromosomal territory, compartment, topologically associating domain, and chromatin loop. These structures play crucial roles in regulating gene expression, cell differentiation, and disease progression. Deciphering the principles underlying 3D genome folding and the molecular mechanisms governing cell fate determination remains a challenge. With advancements in high-throughput sequencing and imaging techniques, the hierarchical organization and functional roles of higher-order chromatin structures have been gradually illuminated. This review systematically discussed the structural hierarchy of the 3D genome, the effects and mechanisms of cis-regulatory elements interaction in the 3D genome for regulating spatiotemporally specific gene expression, the roles and mechanisms of dynamic changes in 3D chromatin conformation during embryonic development, and the pathological mechanisms of diseases such as congenital developmental abnormalities and cancer, which are attributed to alterations in 3D genome organization and aberrations in key structural proteins. Finally, prospects were made for the research about 3D genome structure, function, and genetic intervention, and the roles in disease development, prevention, and treatment, which may offer some clues for precise diagnosis and treatment of related diseases.

Low-affinity CTCF binding drives transcriptional regulation whereas high-affinity binding encompasses architectural functions

CTCF is a DNA-binding protein which plays critical roles in chromatin structure organization and transcriptional regulation; however, little is known about the functional determinants of different CTCF-binding sites (CBS). Using a conditional mouse model, we have identified one set of CBSs that are lost upon CTCF depletion (lost CBSs) and another set that persists (retained CBSs). Retained CBSs are more similar to the consensus CTCF-binding sequence and usually span tandem CTCF peaks. Lost CBSs are enriched at enhancers and promoters and associate with active chromatin marks and higher transcriptional activity. In contrast, retained CBSs are enriched at TAD and loop boundaries. Integration of ChIP-seq and RNA-seq data has revealed that retained CBSs are located at the boundaries between distinct chromatin states, acting as chromatin barriers. Our results provide evidence that transient, lost CBSs are involved in transcriptional regulation, whereas retained CBSs are critical for establishing higher-order chromatin architecture.

Drosophila insulator proteins exhibit in vivo liquid-liquid phase separation properties

Drosophila insulator proteins and the cohesin subunit Rad21 coalesce in vivo to form liquid-droplet condensates, suggesting that liquid–liquid phase separation mediates their function in 3D genome organization.

Mounting evidence implicates liquid–liquid phase separation (LLPS), the condensation of biomolecules into liquid-like droplets in the formation and dissolution of membraneless intracellular organelles (MLOs). Cells use MLOs or condensates for various biological processes, including emergency signaling and spatiotemporal control over steady-state biochemical reactions and heterochromatin formation. Insulator proteins are architectural elements involved in establishing independent domains of transcriptional activity within eukaryotic genomes. In Drosophila, insulator proteins form nuclear foci known as insulator bodies in response to osmotic stress. However, the mechanism through which insulator proteins assemble into bodies is yet to be investigated. Here, we identify signatures of LLPS by insulator bodies, including high disorder tendency in insulator proteins, scaffold–client–dependent assembly, extensive fusion behavior, sphericity, and sensitivity to 1,6-hexanediol. We also show that the cohesin subunit Rad21 is a component of insulator bodies, adding to the known insulator protein constituents and γH2Av. Our data suggest a concerted role of cohesin and insulator proteins in insulator body formation and under physiological conditions. We propose a mechanism whereby these architectural proteins modulate 3D genome organization through LLPS.

A Tremendous Reorganization Journey for the 3D Chromatin Structure from Gametes to Embryos

The 3D chromatin structure within the nucleus is important for gene expression regulation and correct developmental programs. Recently, the rapid development of low-input chromatin conformation capture technologies has made it possible to study 3D chromatin structures in gametes, zygotes and early embryos in a variety of species, including flies, vertebrates and mammals. There are distinct 3D chromatin structures within the male and female gametes. Following the fertilization of male and female gametes, fertilized eggs undergo drastic epigenetic reprogramming at multi levels, including the 3D chromatin structure, to convert the terminally differentiated gamete state into the totipotent state, which can give rise to an individual. However, to what extent the 3D chromatin structure reorganization is evolutionarily conserved and what the underlying mechanisms are for the tremendous reorganization in early embryos remain elusive. Here, we review the latest findings on the 3D chromatin structure reorganization during embryogenesis, and discuss the convergent and divergent reprogramming patterns and key molecular mechanisms for the 3D chromatin structure reorganization from gametes to embryos in different species. These findings shed light on how the 3D chromatin structure reorganization contribute to embryo development in different species. The findings also indicate the role of the 3D chromatin structure on the acquisition of totipotent developmental potential.

CTCF acetylation at lysine 20 is required for the early cardiac mesoderm differentiation of embryonic stem cells

The CCCTC-binding factor (CTCF) protein and its modified forms regulate gene expression and genome organization. However, information on CTCF acetylation and its biological function is still lacking. Here, we show that CTCF can be acetylated at lysine 20 (CTCF-K20) by CREB-binding protein (CBP) and deacetylated by histone deacetylase 6 (HDAC6). CTCF-K20 is required for the CTCF interaction with CBP. A CTCF point mutation at lysine 20 had no effect on self-renewal but blocked the mesoderm differentiation of mouse embryonic stem cells (mESCs). The CTCF-K20 mutation reduced CTCF binding to the promoters and enhancers of genes associated with early cardiac mesoderm differentiation, resulting in diminished chromatin accessibility and decreased enhancer-promoter interactions, impairing gene expression. In summary, this study reveals the important roles of CTCF-K20 in regulating CTCF genomic functions and mESC differentiation into mesoderm.

Cohesin is required for long-range enhancer action at the Shh locus

The regulatory landscapes of developmental genes in mammals can be complex, with enhancers spread over many hundreds of kilobases. It has been suggested that three-dimensional genome organization, particularly topologically associating domains formed by cohesin-mediated loop extrusion, is important for enhancers to act over such large genomic distances. By coupling acute protein degradation with synthetic activation by targeted transcription factor recruitment, here we show that cohesin, but not CTCF, is required for activation of the target gene Shh by distant enhancers in mouse embryonic stem cells. Cohesin is not required for activation directly at the promoter or by an enhancer located closer to the Shh gene. Our findings support the hypothesis that chromatin compaction via cohesin-mediated loop extrusion allows for genes to be activated by enhancers that are located many hundreds of kilobases away in the linear genome and suggests that cohesin is dispensable for enhancers located more proximally.

3D or Not 3D: Shaping the Genome during Development

One of the most fundamental questions in developmental biology is how one fertilized cell can give rise to a fully mature organism and how gene regulation governs this process. Precise spatiotemporal gene expression is required for development and is believed to be achieved through a complex interplay of sequence-specific information, epigenetic modifications, trans-acting factors, and chromatin folding. Here we review the role of chromatin folding during development, the mechanisms governing 3D genome organization, and how it is established in the embryo. Furthermore, we discuss recent advances and debated questions regarding the contribution of the 3D genome to gene regulation during organogenesis. Finally, we describe the mechanisms that can reshape the 3D genome, including disease-causing structural variations and the emerging view that transposable elements contribute to chromatin organization.

Nuclear receptor activation shapes spatial genome organization essential for gene expression control: lessons learned from the vitamin D receptor

Spatial genome organization is tightly controlled by several regulatory mechanisms and is essential for gene expression control. Nuclear receptors are ligand-activated transcription factors that modulate physiological and pathophysiological processes and are primary pharmacological targets. DNA binding of the important loop-forming insulator protein CCCTC-binding factor (CTCF) was modulated by 1α,25-dihydroxyvitamin D3 (1,25(OH)2D3). We performed CTCF HiChIP assays to produce the first genome-wide dataset of CTCF long-range interactions in 1,25(OH)2D3-treated cells, and to determine whether dynamic changes of spatial chromatin interactions are essential for fine-tuning of nuclear receptor signaling. We detected changes in 3D chromatin organization upon vitamin D receptor (VDR) activation at 3.1% of all observed CTCF interactions. VDR binding was enriched at both differential loop anchors and within differential loops. Differential loops were observed in several putative functional roles including TAD border formation, promoter-enhancer looping, and establishment of VDR-responsive insulated neighborhoods. Vitamin D target genes were enriched in differential loops and at their anchors. Secondary vitamin D effects related to dynamic chromatin domain changes were linked to location of downstream transcription factors in differential loops. CRISPR interference and loop anchor deletion experiments confirmed the functional relevance of nuclear receptor ligand-induced adjustments of the chromatin 3D structure for gene expression regulation.

Implications of Dosage Deficiencies in CTCF and Cohesin on Genome Organization, Gene Expression, and Human Neurodevelopment

Properly organizing DNA within the nucleus is critical to ensure normal downstream nuclear functions. CTCF and cohesin act as major architectural proteins, working in concert to generate thousands of high-intensity chromatin loops. Due to their central role in loop formation, a massive research effort has been dedicated to investigating the mechanism by which CTCF and cohesin create these loops. Recent results lead to questioning the direct impact of CTCF loops on gene expression. Additionally, results of controlled depletion experiments in cell lines has indicated that genome architecture may be somewhat resistant to incomplete deficiencies in CTCF or cohesin. However, heterozygous human genetic deficiencies in CTCF and cohesin have illustrated the importance of their dosage in genome architecture, cellular processes, animal behavior, and disease phenotypes. Thus, the importance of considering CTCF or cohesin levels is especially made clear by these heterozygous germline variants that characterize genetic syndromes, which are increasingly recognized in clinical practice. Defined primarily by developmental delay and intellectual disability, the phenotypes of CTCF and cohesin deficiency illustrate the importance of architectural proteins particularly in neurodevelopment. We discuss the distinct roles of CTCF and cohesin in forming chromatin loops, highlight the major role that dosage of each protein plays in the amplitude of observed effects on gene expression, and contrast these results to heterozygous mutation phenotypes in murine models and clinical patients. Insights highlighted by this comparison have implications for future research into these newly emerging genetic syndromes.

CTCF supports preferentially short lamina-associated domains

More than one third of the mammalian genome is in a close association with the nuclear lamina, thus these genomic regions were termed lamina-associated domains (LADs). This association is fundamental for many aspects of chromatin biology including transcription, replication, and DNA damage repair. LADs association with the nuclear envelope is thought to be dependent on two major mechanisms: The first mechanism is the interaction between nuclear membrane proteins such as LBR with heterochromatin modifications that are enriched in LADs chromatin. The second mechanism is based on proteins that bind the borders of the LADs and support the association of the LADs with the nuclear envelope. Two factors were suggested to support the second mechanism: CCCTC-binding factor (CTCF) and YY1 based on their enriched binding to LADs borders. However, this mechanism has not been proven yet at a whole genome level. Here, to test if CTCF supports the LADs landscape, we generated melanoma cells with a partial loss of function (pLoF) of CTCF by the CRISPR-Cas9 system and determined the LADs landscape by lamin B ChIP-seq analysis. We found that under regular growth conditions, CTCF pLoF led to modest changes in the LADs landscape that included an increase in the signal of 2% of the LADs and a decrease in the signal of 8% of the LADs. However, CTCF importance for the LADs landscape was much higher upon induction of a chromatin stress. We induced chromatin stress by inhibiting RNA polymerase II, an intervention that is known to alter chromatin compaction and supercoiling. Notably, only in CTCF pLoF cells, the chromatin stress led to the dissociation of 7% of the LADs from the lamina. The CTCF-dependent LADs had almost three times shorter median length than the non-affected LADs, were enriched in CTCF binding at their borders, and were higher in their facultative-status (cell-type specific). Thus, it appears that CTCF is a key factor in facilitating the association of short facultative LADs with the nuclear lamina upon chromatin stress.

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

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

CTCF knockout in zebrafish induces alterations in regulatory landscapes and developmental gene expression

Coordinated chromatin interactions between enhancers and promoters are critical for gene regulation. The architectural protein CTCF mediates chromatin looping and is enriched at the boundaries of topologically associating domains (TADs), which are sub-megabase chromatin structures. In vitro CTCF depletion leads to a loss of TADs but has only limited effects over gene expression, challenging the concept that CTCF-mediated chromatin structures are a fundamental requirement for gene regulation. However, how CTCF and a perturbed chromatin structure impacts gene expression during development remains poorly understood. Here we link the loss of CTCF and gene regulation during patterning and organogenesis in a ctcf knockout zebrafish model. CTCF absence leads to loss of chromatin structure and affects the expression of thousands of genes, including many developmental regulators. Our results demonstrate the essential role of CTCF in providing the structural context for enhancer-promoter interactions, thus regulating developmental genes.

Liver-Specific Deletion of Mouse CTCF Leads to Hepatic Steatosis via Augmented PPARγ Signaling

BACKGROUND & AIMS

The liver is the major organ for metabolizing lipids, and malfunction of the liver leads to various diseases. Nonalcoholic fatty liver disease is rapidly becoming a major health concern worldwide and is characterized by abnormal retention of excess lipids in the liver. CCCTC-binding factor (CTCF) is a highly conserved zinc finger protein that regulates higher-order chromatin organization and is involved in various gene regulation processes. Here, we sought to determine the physiological role of CTCF in hepatic lipid metabolism.

METHODS

We generated liver-specific, CTCF-ablated and/or CD36 whole-body knockout mice. Overexpression or knockdown of peroxisome proliferator-activated receptor (PPAR)γ in the liver was achieved using adenovirus. Mice were examined for development of hepatic steatosis and inflammation. RNA sequencing was performed to identify genes affected by CTCF depletion. Genome-wide occupancy of H3K27 acetylation, PPARγ, and CTCF were analyzed by chromatin immunoprecipitation sequencing. Genome-wide chromatin interactions were analyzed by in situ Hi-C.

RESULTS

Liver-specific, CTCF-deficient mice developed hepatic steatosis and inflammation when fed a standard chow diet. Global analysis of the transcriptome and enhancer landscape revealed that CTCF-depleted liver showed enhanced accumulation of PPARγ in the nucleus, which leads to increased expression of its downstream target genes, including fat storage-related gene CD36, which is involved in the lipid metabolic process. Hepatic steatosis developed in liver-specific, CTCF-deficient mice was ameliorated by repression of PPARγ via pharmacologic blockade or adenovirus-mediated knockdown, but hardly rescued by additional knockout of CD36.

CONCLUSIONS

Our data indicate that liver-specific deletion of CTCF leads to hepatosteatosis through augmented PPARγ DNA-binding activity, which up-regulates its downstream target genes associated with the lipid metabolic process.

Drosophila architectural protein CTCF is not essential for fly survival and is able to function independently of CP190

CTCF is the most likely ancestor of proteins that contain large clusters of C2H2 zinc finger domains (C2H2) and is conserved among most bilateral organisms. In mammals, CTCF functions as the main architectural protein involved in the organization of topology-associated domains (TADs). In vertebrates and Drosophila, CTCF is involved in the regulation of homeotic genes. Previously, it was found that null mutations in the dCTCF gene died as pharate adults, which failed to eclose from their pupal case, or shortly after hatching of adults. Here, we obtained several new null dCTCF mutations and found that the complete inactivation of dCTCF appears is limited mainly to phenotypic manifestations of the Abd-B gene and fertility of adult flies. Many modifiers that are not associated with an independent phenotypic manifestation can significantly enhance the expressivity of the null dCTCF mutations, indicating that other architectural proteins are able to functionally compensate for dCTCF inactivation in Drosophila. We also mapped the 715-735 aa region of dCTCF as being essential for the interaction with the BTB (Broad-Complex, Tramtrack, and Bric a brac) and microtubule-targeting (M) domains of the CP190 protein, which binds to many architectural proteins. However, the mutational analysis showed that the interaction with CP190 was not important for the functional activity of dCTCF in vivo.

Induction of a chromatin boundary in vivo upon insertion of a TAD border

Mammalian genomes are partitioned into sub-megabase to megabase-sized units of preferential interactions called topologically associating domains or TADs, which are likely important for the proper implementation of gene regulatory processes. These domains provide structural scaffolds for distant cis regulatory elements to interact with their target genes within the three-dimensional nuclear space and architectural proteins such as CTCF as well as the cohesin complex participate in the formation of the boundaries between them. However, the importance of the genomic context in providing a given DNA sequence the capacity to act as a boundary element remains to be fully investigated. To address this question, we randomly relocated a topological boundary functionally associated with the mouse HoxD gene cluster and show that it can indeed act similarly outside its initial genomic context. In particular, the relocated DNA segment recruited the required architectural proteins and induced a significant depletion of contacts between genomic regions located across the integration site. The host chromatin landscape was re-organized, with the splitting of the TAD wherein the boundary had integrated. These results provide evidence that topological boundaries can function independently of their site of origin, under physiological conditions during mouse development.

During development, enhancer sequences tightly regulate the spatio-temporal expression of target genes often located hundreds of kilobases away. This complex process is made possible by the folding of chromatin into domains, which are separated from one another by specific genomic regions referred to as boundaries. In order to understand whether such boundary sequences require their particular genomic contexts to achieve their isolating effect, we analyzed the impact of introducing one such boundary, taken from the HoxD locus, into a distinct topological domain. We show that this ectopic boundary splits the host domain into two sub-domains and affects the expression levels of a neighboring gene. We conclude that this sequence can work independently from its genomic context and thus carries all the information necessary to act as a boundary element.

Genome folding through loop extrusion by SMC complexes

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

The combined action of CTCF and its testis-specific paralog BORIS is essential for spermatogenesis

CTCF is a key organizer of the 3D genome. Its specialized paralog, BORIS, heterodimerizes with CTCF but is expressed only in male germ cells and in cancer states. Unexpectedly, BORIS-null mice have only minimal germ cell defects. To understand the CTCF-BORIS relationship, mouse models with varied CTCF and BORIS levels were generated. Whereas Ctcf+/+Boris+/+, Ctcf+/-Boris+/+, and Ctcf+/+Boris-/- males are fertile, Ctcf+/-Boris-/- (Compound Mutant; CM) males are sterile. Testes with combined depletion of both CTCF and BORIS show reduced size, defective meiotic recombination, increased apoptosis, and malformed spermatozoa. Although CM germ cells exhibit only 25% of CTCF WT expression, chromatin binding of CTCF is preferentially lost from CTCF-BORIS heterodimeric sites. Furthermore, CM testes lose the expression of a large number of spermatogenesis genes and gain the expression of developmentally inappropriate genes that are "toxic" to fertility. Thus, a combined action of CTCF and BORIS is required to both repress pre-meiotic genes and activate post-meiotic genes for a complete spermatogenesis program.

MACMIC Reveals A Dual Role of CTCF in Epigenetic Regulation of Cell Identity Genes

Numerous studies of relationship between epigenomic features have focused on their strong correlation across the genome, likely because such relationship can be easily identified by many established methods for correlation analysis. However, two features with little correlation may still colocalize at many genomic sites to implement important functions. There is no bioinformatic tool for researchers to specifically identify such feature pairs. Here, we develop a method to identify feature pairs in which two features have maximal colocalization minimal correlation (MACMIC) across the genome. By MACMIC analysis of 3306 feature pairs in 16 human cell types, we reveal a dual role of CCCTC-binding factor (CTCF) in epigenetic regulation of cell identity genes. Although super-enhancers are associated with activation of target genes, only a subset of super-enhancers colocalized with CTCF regulate cell identity genes. At super-enhancers colocalized with CTCF, CTCF is required for the active marker H3K27ac in cell types requiring the activation, and also required for the repressive marker H3K27me3 in other cell types requiring repression. Our work demonstrates the biological utility of the MACMIC analysis and reveals a key role for CTCF in epigenetic regulation of cell identity. The code for MACMIC is available at https://github.com/bxia888/MACMIC.

CTCF loss has limited effects on global genome architecture in Drosophila despite critical regulatory functions

Vertebrate genomes are partitioned into contact domains defined by enhanced internal contact frequency and formed by two principal mechanisms: compartmentalization of transcriptionally active and inactive domains, and stalling of chromosomal loop-extruding cohesin by CTCF bound at domain boundaries. While Drosophila has widespread contact domains and CTCF, it is currently unclear whether CTCF-dependent domains exist in flies. We genetically ablate CTCF in Drosophila and examine impacts on genome folding and transcriptional regulation in the central nervous system. We find that CTCF is required to form a small fraction of all domain boundaries, while critically controlling expression patterns of certain genes and supporting nervous system function. We also find that CTCF recruits the pervasive boundary-associated factor Cp190 to CTCF-occupied boundaries and co-regulates a subset of genes near boundaries together with Cp190. These results highlight a profound difference in CTCF-requirement for genome folding in flies and vertebrates, in which a large fraction of boundaries are CTCF-dependent and suggest that CTCF has played mutable roles in genome architecture and direct gene expression control during metazoan evolution.

Chromatin topology and the timing of enhancer function at the HoxD locus

The HoxD gene cluster is critical for proper limb formation in tetrapods. In the emerging limb buds, different subgroups of Hoxd genes respond first to a proximal regulatory signal, then to a distal signal that organizes digits. These two regulations are exclusive from one another and emanate from two distinct topologically associating domains (TADs) flanking HoxD, both containing a range of appropriate enhancer sequences. The telomeric TAD (T-DOM) contains several enhancers active in presumptive forearm cells and is divided into two sub-TADs separated by a CTCF-rich boundary, which defines two regulatory submodules. To understand the importance of this particular regulatory topology to control Hoxd gene transcription in time and space, we either deleted or inverted this sub-TAD boundary, eliminated the CTCF binding sites, or inverted the entire T-DOM to exchange the respective positions of the two sub-TADs. The effects of such perturbations on the transcriptional regulation of Hoxd genes illustrate the requirement of this regulatory topology for the precise timing of gene activation. However, the spatial distribution of transcripts was eventually resumed, showing that the presence of enhancer sequences, rather than either their exact topology or a particular chromatin architecture, is the key factor. We also show that the affinity of enhancers to find their natural target genes can overcome the presence of both a strong TAD border and an unfavorable orientation of CTCF sites.

Robust and efficient gene regulation through localized nuclear microenvironments

Developmental enhancers drive gene expression in specific cell types during animal development. They integrate signals from many different sources mediated through the binding of transcription factors, producing specific responses in gene expression. Transcription factors often bind low-affinity sequences for only short durations. How brief, low-affinity interactions drive efficient transcription and robust gene expression is a central question in developmental biology. Localized high concentrations of transcription factors have been suggested as a possible mechanism by which to use these enhancer sites effectively. Here, we discuss the evidence for such transcriptional microenvironments, mechanisms for their formation and the biological consequences of such sub-nuclear compartmentalization for developmental decisions and evolution.

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

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

Functional signatures of evolutionarily young CTCF binding sites

Background

The introduction of novel CTCF binding sites in gene regulatory regions in the rodent lineage is partly the effect of transposable element expansion, particularly in the murine lineage. The exact mechanism and functional impact of evolutionarily novel CTCF binding sites are not yet fully understood. We investigated the impact of novel subspecies-specific CTCF binding sites in two Mus genus subspecies, Mus musculus domesticus and Mus musculus castaneus, that diverged 0.5 million years ago.

Results

CTCF binding site evolution is influenced by the action of the B2-B4 family of transposable elements independently in both lineages, leading to the proliferation of novel CTCF binding sites. A subset of evolutionarily young sites may harbour transcriptional functionality as evidenced by the stability of their binding across multiple tissues in M. musculus domesticus (BL6), while overall the distance of subspecies-specific CTCF binding to the nearest transcription start sites and/or topologically associated domains (TADs) is largely similar to musculus-common CTCF sites. Remarkably, we discovered a recurrent regulatory architecture consisting of a CTCF binding site and an interferon gene that appears to have been tandemly duplicated to create a 15-gene cluster on chromosome 4, thus forming a novel BL6 specific immune locus in which CTCF may play a regulatory role.

Conclusions

Our results demonstrate that thousands of CTCF binding sites show multiple functional signatures rapidly after incorporation into the genome.

Loss of U11 small nuclear RNA in the developing mouse limb results in micromelia

Disruption of the minor spliceosome due to mutations in RNU4ATAC is linked to primordial dwarfism in microcephalic osteodysplastic primordial dwarfism type 1, Roifman syndrome, and Lowry-Wood syndrome. Similarly, primordial dwarfism in domesticated animals is linked to positive selection in minor spliceosome components. Despite being vital for limb development and size regulation, its role remains unexplored. Here, we disrupt minor spliceosome function in the developing mouse limb by ablating one of its essential components, U11 small nuclear RNA, which resulted in micromelia. Notably, earlier loss of U11 corresponded to increased severity. We find that limb size is reduced owing to elevated minor intron retention in minor intron-containing genes that regulate cell cycle. As a result, limb progenitor cells experience delayed prometaphase-to-metaphase transition and prolonged S-phase. Moreover, we observed death of rapidly dividing, distally located progenitors. Despite cell cycle defects and cell death, the spatial expression of key limb patterning genes was maintained. Overall, we show that the minor spliceosome is required for limb development via size control potentially shared in disease and domestication.

The emergence of genome architecture and zygotic genome activation

The fusion of two transcriptionally silent gametes, egg and sperm, generates a totipotent zygote that activates zygotic transcription to support further development. Although the molecular details of zygotic genome activation (ZGA) are not well understood in most species, an emerging concept is that one or more pioneer transcription factors trigger zygotic transcription. Concomitantly, extensive changes in 3D chromatin organization occur during development. In this review, we discuss recent advances in understanding when and how genome architecture emerges in early metazoan embryos, how the zygotic genome is activated, and how these events might be coordinated. We also highlight some of the unknowns that may be critical to address in the future.

A complex regulatory landscape involved in the development of mammalian external genitals

Developmental genes are often controlled by large regulatory landscapes matching topologically associating domains (TADs). In various contexts, the associated chromatin backbone is modified by specific enhancer-enhancer and enhancer-promoter interactions. We used a TAD flanking the mouse HoxD cluster to study how these regulatory architectures are formed and deconstructed once their function achieved. We describe this TAD as a functional unit, with several regulatory sequences acting together to elicit a transcriptional response. With one exception, deletion of these sequences didn't modify the transcriptional outcome, a result at odds with a conventional view of enhancer function. The deletion and inversion of a CTCF site located near these regulatory sequences did not affect transcription of the target gene. Slight modifications were nevertheless observed, in agreement with the loop extrusion model. We discuss these unexpected results considering both conventional and alternative explanations relying on the accumulation of poorly specific factors within the TAD backbone.

Interplay between genetics and epigenetics in osteoarthritis

Research into the molecular genetics of osteoarthritis (OA) has been substantially bolstered in the past few years by the implementation of powerful genome-wide scans that have revealed a large number of novel risk loci associated with the disease. This refreshing wave of discovery has occurred concurrently with epigenetic studies of joint tissues that have examined DNA methylation, histone modifications and regulatory RNAs. These epigenetic analyses have involved investigations of joint development, homeostasis and disease and have used both human samples and animal models. What has become apparent from a comparison of these two complementary approaches is that many OA genetic risk signals interact with, map to or correlate with epigenetic mediators. This discovery implies that epigenetic mechanisms, and their effect on gene expression, are a major conduit through which OA genetic risk polymorphisms exert their functional effects. This observation is particularly exciting as it provides mechanistic insight into OA susceptibility. Furthermore, this knowledge reveals avenues for attenuating the negative effect of risk-conferring alleles by exposing the epigenome as an exploitable target for therapeutic intervention in OA.

Functional genomic approaches to elucidate the role of enhancers during development

Successful development depends on the precise tissue-specific regulation of genes by enhancers, genetic elements that act as switches to control when and where genes are expressed. Because enhancers are critical for development, and the majority of disease-associated mutations reside within enhancers, it is essential to understand which sequences within enhancers are important for function. Advances in sequencing technology have enabled the rapid generation of genomic data that predict putative active enhancers, but functionally validating these sequences at scale remains a fundamental challenge. Herein, we discuss the power of genome-wide strategies used to identify candidate enhancers, and also highlight limitations and misconceptions that have arisen from these data. We discuss the use of massively parallel reporter assays to test enhancers for function at scale. We also review recent advances in our ability to study gene regulation during development, including CRISPR-based tools to manipulate genomes and single-cell transcriptomics to finely map gene expression. Finally, we look ahead to a synthesis of complementary genomic approaches that will advance our understanding of enhancer function during development. This article is categorized under: Physiology > Mammalian Physiology in Health and Disease Developmental Biology > Developmental Processes in Health and Disease Laboratory Methods and Technologies > Genetic/Genomic Methods.

Why a Constant Number of Vertebrae? Digital Control of Segmental Identity during Vertebrate Development: The Somite Cycle Controls a Digital, Chromatin-Based Counter That Defines Segmental Identity and Body Plans in Vertebrate Animals

It is not understood how the numbers and identities of vertebrae are controlled during mammalian development. The remarkable robustness and conservation of segmental numbers may suggest the digital nature of the underlying process. The study proposes a mechanism that allows cells to obtain and store the segmental information in digital form, and to produce a pattern of chromatin accessibility that in turn regulates Hox gene expression specific to the metameric segment. The model requires that a regulatory element be present such that the number of occurrences of the motif between two consecutive Hox genes equals the number of segments under the control of the anterior gene. This is true for the recently discovered hydroxyl radical cleavage 3bp-periodic (HRC3) motif, associated with histone modifications and developmental genes. The finding not only allows the correct prediction of the numbers of segments using only sequence information, but also resolves the 40-year-old enigma of the function of temporal and spatial collinearity of Hox genes. The logic of the mechanism is illustrated in the attached animated video. How different aspects of the proposed mechanism can be tested experimentally is also discussed.

Developmentally regulated Shh expression is robust to TAD perturbations

Topologically associating domains (TADs) have been proposed to both guide and constrain enhancer activity. Shh is located within a TAD known to contain all its enhancers. To investigate the importance of chromatin conformation and TAD integrity on developmental gene regulation, we have manipulated the Shh TAD - creating internal deletions, deleting CTCF sites, and deleting and inverting sequences at TAD boundaries. Chromosome conformation capture and fluorescence in situ hybridisation assays were used to investigate the changes in chromatin conformation that result from these manipulations. Our data suggest that these substantial alterations in TAD structure have no readily detectable effect on Shh expression patterns or levels of Shh expression during development - except where enhancers are deleted - and result in no detectable phenotypes. Only in the case of a larger deletion at one TAD boundary could ectopic influence of the Shh limb enhancer be detected on a gene (Mnx1) in the neighbouring TAD. Our data suggests that, contrary to expectations, the developmental regulation of Shh expression is remarkably robust to TAD perturbations.

The Role of Insulation in Patterning Gene Expression

Development is orchestrated by regulatory elements that turn genes ON or OFF in precise spatial and temporal patterns. Many safety mechanisms prevent inappropriate action of a regulatory element on the wrong gene promoter. In flies and mammals, dedicated DNA elements (insulators) recruit protein factors (insulator binding proteins, or IBPs) to shield promoters from regulatory elements. In mammals, a single IBP called CCCTC-binding factor (CTCF) is known, whereas genetic and biochemical analyses in Drosophila have identified a larger repertoire of IBPs. How insulators function at the molecular level is not fully understood, but it is currently thought that they fold chromosomes into conformations that affect regulatory element-promoter communication. Here, we review the discovery of insulators and describe their properties. We discuss recent genetic studies in flies and mice to address the question: Is gene insulation important for animal development? Comparing and contrasting observations in these two species reveal that they have different requirements for insulation, but that insulation is a conserved and critical gene regulation strategy.

Impact of genome architecture on the functional activation and repression of Hox regulatory landscapes

BACKGROUND

The spatial organization of the mammalian genome relies upon the formation of chromatin domains of various scales. At the level of gene regulation in cis, collections of enhancer sequences define large regulatory landscapes that usually match with the presence of topologically associating domains (TADs). These domains often contain ranges of enhancers displaying similar or related tissue specificity, suggesting that in some cases, such domains may act as coherent regulatory units, with a global on or off state. By using the HoxD gene cluster, which specifies the topology of the developing limbs via highly orchestrated regulation of gene expression, as a paradigm, we investigated how the arrangement of regulatory domains determines their activity and function.

RESULTS

Proximal and distal cells in the developing limb express different levels of Hoxd genes, regulated by flanking 3' and 5' TADs, respectively. We characterized the effect of large genomic rearrangements affecting these two TADs, including their fusion into a single chromatin domain. We show that, within a single hybrid TAD, the activation of both proximal and distal limb enhancers globally occurred as when both TADs are intact. However, the activity of the 3' TAD in distal cells is generally increased in the fused TAD, when compared to wild type where it is silenced. Also, target gene activity in distal cells depends on whether or not these genes had previously responded to proximal enhancers, which determines the presence or absence of H3K27me3 marks. We also show that the polycomb repressive complex 2 is mainly recruited at the Hox gene cluster and can extend its coverage to far-cis regulatory sequences as long as confined to the neighboring TAD structure.

CONCLUSIONS

We conclude that antagonistic limb proximal and distal enhancers can exert their specific effects when positioned into the same TAD and in the absence of their genuine target genes. We also conclude that removing these target genes reduced the coverage of a regulatory landscape by chromatin marks associated with silencing, which correlates with its prolonged activity in time.

The constrained architecture of mammalian Hox gene clusters

Vertebrate Hox genes are clustered. This organization has a functional relevance, as the transcription of each gene in time and space depends upon its relative position within the gene cluster. Hox clusters display a high organization, and all genes are transcribed from the same DNA strand. Here, we investigate the importance of this uniform transcriptional polarity by engineering alleles where one or several transcription units are inverted, with or without a CTCF site. We observe that inversions are likely detrimental to the proper implementation of this genetic system. We propose that the enhanced organization of Hox clusters in vertebrates evolved in conjunction with the emergence of global gene regulation to optimize a coordinated response of selected subsets of target genes.

In many animal species with a bilateral symmetry, Hox genes are clustered either at one or at several genomic loci. This organization has a functional relevance, as the transcriptional control applied to each gene depends upon its relative position within the gene cluster. It was previously noted that vertebrate Hox clusters display a much higher level of genomic organization than their invertebrate counterparts. The former are always more compact than the latter, they are generally devoid of repeats and of interspersed genes, and all genes are transcribed by the same DNA strand, suggesting that particular factors constrained these clusters toward a tighter structure during the evolution of the vertebrate lineage. Here, we investigate the importance of uniform transcriptional orientation by engineering several alleles within the HoxD cluster, such as to invert one or several transcription units, with or without a neighboring CTCF site. We observe that the association between the tight structure of mammalian Hox clusters and their regulation makes inversions likely detrimental to the proper implementation of this complex genetic system. We propose that the consolidation of Hox clusters in vertebrates, including transcriptional polarity, evolved in conjunction with the emergence of global gene regulation via the flanking regulatory landscapes, to optimize a coordinated response of selected subsets of target genes in cis.

Genome-wide meta-analysis identifies genetic locus on chromosome 9 associated with Modic changes

BACKGROUND

Low back pain (LBP) is a common disabling condition. Lumbar disc degeneration (LDD) may be a contributing factor for LBP. Modic change (MC), a distinct phenotype of LDD, is presented as a pathological bone marrow signal change adjacent to vertebral endplate on MRI. It is strongly associated with LBP and has heritability around 30%. Our objective was to identify genetic loci associated with MC using a genome-wide meta-analysis.

METHODS

Presence of MC was evaluated in lumbar MRI in the Northern Finland Birth Cohort 1966 (n=1182) and TwinsUK (n=647). Genome-wide association analyses were carried out using linear regression model. Inverse-variance weighting approach was used in the meta-analysis.

RESULTS

A locus associated with MC (p<5e-8) was found on chromosome 9 with the lead SNP rs1934268 in an intron of the PTPRD gene. It is located in the binding region of BCL11A, SPI1 and PBX3 transcription factors. The SNP was nominally associated with LBP in TwinsUK (p=0.001) but not associated in the UK Biobank (p=0.914). Suggestive signals (p<1e-5) were identified near XKR4, SCIN, MGMT, DLG2, ZNF184 and OPRK1.

CONCLUSION

PTPRD is a novel candidate gene for MC that may act via the development of cartilage or nervous system; further work is needed to define the mechanisms underlying the pathways leading to development of MC. This is the first genome-wide meta-analysis of MC, and the results pave the way for further studies on the genetic factors underlying the various features of spine degeneration and LBP.

The interdependence of gene-regulatory elements and the 3D genome

Imaging studies, high-resolution chromatin conformation maps, and genome-wide occupancy data of architectural proteins have revealed that genome topology is tightly intertwined with gene expression. Cross-talk between gene-regulatory elements is often organized within insulated neighborhoods, and regulatory cues that induce transcriptional changes can reshape chromatin folding patterns and gene positioning within the nucleus. The cause-consequence relationship of genome architecture and gene expression is intricate, and its molecular mechanisms are under intense investigation. Here, we review the interdependency of transcription and genome organization with emphasis on enhancer-promoter contacts in gene regulation.

Targeting CTCFL/BORIS for the immunotherapy of cancer

Cancer vaccines have great potential in the fight against metastatic malignancies. Current anti-tumor immunotherapy is hindered by existing tolerance to tumor-associated antigens (TAA) and tumor escape using various mechanisms, highlighting the need for improved targets for immunotherapy. The cancer-testis antigen CTCFL/BORIS was discovered 16 years ago and possesses all features necessary for an ideal TAA. Recently CTCFL/BORIS has received additional attention as a target expressed in cancer stem cells (CSC). These cells drive tumor growth recurrence, metastasis, and treatment resistance. CTCFL/BORIS silencing leads to senescence and death of CSC. Therefore, an immunotherapeutic strategy that targets CTCFL/BORIS may lead to the selective destruction of CSC and potential eradication of metastatic disease. The high immunotherapeutic potential of CTCFL/BORIS antigen was shown in a stringent 4T1 mouse model of breast cancer. Using these highly metastatic, poorly immunogenic carcinoma cells inoculated into T-helper2 prone mice, we showed that DC fed with recombinant CTCFL/BORIS as an immunogen inhibited tumor growth and reduced the number of metastases in distant organs. About 20% of CTCFL/BORIS immunized animals were tumor free. 50% of animals remained metastasis free. Those having metastasis showed at least tenfold fewer metastases compared to controls. In a rat model of breast cancer, we showed that alphavirus-based CTCFL/BORIS immunotherapy was capable of cancer elimination as we were able to cure 50% of animals. Based on the above data, we believe that translation of CTCFL/BORIS-targeting immunotherapeutic strategies to the clinic will provide new avenues for improving survival of breast cancer patients with advanced metastatic disease.

CTCF Expression is Essential for Somatic Cell Viability and Protection Against Cancer

CCCTC-binding factor (CTCF) is a conserved transcription factor that performs diverse roles in transcriptional regulation and chromatin architecture. Cancer genome sequencing reveals diverse acquired mutations in CTCF, which we have shown functions as a tumour suppressor gene. While CTCF is essential for embryonic development, little is known of its absolute requirement in somatic cells and the consequences of CTCF haploinsufficiency. We examined the consequences of CTCF depletion in immortalised human and mouse cells using shRNA knockdown and CRISPR/Cas9 genome editing as well as examined the growth and development of heterozygous Ctcf (Ctcf^+/-) mice. We also analysed the impact of CTCF haploinsufficiency by examining gene expression changes in CTCF-altered endometrial carcinoma. Knockdown and CRISPR/Cas9-mediated editing of CTCF reduced the cellular growth and colony-forming ability of K562 cells. CTCF knockdown also induced cell cycle arrest and a pro-survival response to apoptotic insult. However, in p53 shRNA-immortalised Ctcf^+/- MEFs we observed the opposite: increased cellular proliferation, colony formation, cell cycle progression, and decreased survival after apoptotic insult compared to wild-type MEFs. CRISPR/Cas9-mediated targeting in Ctcf^+/- MEFs revealed a predominance of in-frame microdeletions in Ctcf in surviving clones, however protein expression could not be ablated. Examination of CTCF mutations in endometrial cancers showed locus-specific alterations in gene expression due to CTCF haploinsufficiency, in concert with downregulation of tumour suppressor genes and upregulation of estrogen-responsive genes. Depletion of CTCF expression imparts a dramatic negative effect on normal cell function. However, CTCF haploinsufficiency can have growth-promoting effects consistent with known cancer hallmarks in the presence of additional genetic hits. Our results confirm the absolute requirement for CTCF expression in somatic cells and provide definitive evidence of CTCF's role as a haploinsufficient tumour suppressor gene. CTCF genetic alterations in endometrial cancer indicate that gene dysregulation is a likely consequence of CTCF loss, contributing to, but not solely driving cancer growth.

Similarities and differences in the regulation of HoxD genes during chick and mouse limb development

In all tetrapods examined thus far, the development and patterning of limbs require the activation of gene members of the HoxD cluster. In mammals, they are regulated by a complex bimodal process that controls first the proximal patterning and then the distal structure. During the shift from the former to the latter regulation, this bimodal regulatory mechanism allows the production of a domain with low Hoxd gene expression, at which both telomeric (T-DOM) and centromeric regulatory domains (C-DOM) are silent. These cells generate the future wrist and ankle articulations. We analyzed the implementation of this regulatory mechanism in chicken, i.e., in an animal for which large morphological differences exist between fore- and hindlimbs. We report that although this bimodal regulation is globally conserved between the mouse and the chick, some important modifications evolved at least between these two model systems, in particular regarding the activity of specific enhancers, the width of the TAD boundary separating the two regulations, and the comparison between the forelimb versus hindlimb regulatory controls. At least one aspect of these regulations seems to be more conserved between chick and bats than with mouse, which may relate to the extent to which forelimbs and hindlimbs of these various animals differ in their morphologies.

A comparison of Hox gene regulation during the development of limbs in birds and mammals reveals that whereas the characteristic bimodal regulatory system, based on large chromatin domains, is largely conserved between these morphologically distinct structures, some differences are revealed in the way this is implemented in various vertebrates.

The shapes of limbs vary greatly among tetrapod species, even between the forelimbs and hindlimbs of the same animal. Hox genes regulate the proper growth and patterning of tetrapod limbs. In order to evaluate whether variations in the complex regulation of a cluster of Hox genes—the Hoxd genes—during limb development contribute to the differences in limb shape, we compared their transcriptional control during limb bud development in the forelimbs and hindlimbs of mouse and chicken embryos. We found that the regulatory mechanism underlying Hoxd gene expression is highly conserved, but some clear differences exist. For instance, we observed a variation in the topologically associating domain (TAD; a self-interacting genomic region) boundary interval between the mouse and the chick, as well as differences in the activity of a conserved enhancer element situated within the telomeric regulatory domain. In contrast to the mouse, the chicken enhancer has a stronger activity in the forelimb buds than in the hindlimb buds, which is correlated with the striking differences in the mRNA levels of the genes. We conclude that differences in both the timing and duration of TAD activities and in the width of their boundary may parallel the important decrease in Hoxd gene transcription in chick hindlimb buds versus forelimb buds. These differences may also account for the slightly distinct regulatory strategies implemented by mammals and birds at this locus.

Heterogeneous combinatorial expression of Hoxd genes in single cells during limb development

BACKGROUND

Global analyses of gene expression during development reveal specific transcription patterns associated with the emergence of various cell types, tissues, and organs. These heterogeneous patterns are instrumental to ensure the proper formation of the different parts of our body, as shown by the phenotypic effects generated by functional genetic approaches. However, variations at the cellular level can be observed within each structure or organ. In the developing mammalian limbs, expression of Hox genes from the HoxD cluster is differentially controlled in space and time, in cells that will pattern the digits and the forearms. While the Hoxd genes broadly share a common regulatory landscape and large-scale analyses have suggested a homogenous Hox gene transcriptional program, it has not previously been clear whether Hoxd genes are expressed together at the same levels in the same cells.

RESULTS

We report a high degree of heterogeneity in the expression of the Hoxd11 and Hoxd13 genes. We analyzed single-limb bud cell transcriptomes and show that Hox genes are expressed in specific combinations that appear to match particular cell types. In cells giving rise to digits, we find that the expression of the five relevant Hoxd genes (Hoxd9 to Hoxd13) is unbalanced, despite their control by known global enhancers. We also report that specific combinatorial expression follows a pseudo-time sequence, which is established based on the transcriptional diversity of limb progenitors.

CONCLUSIONS

Our observations reveal the existence of distinct combinations of Hoxd genes at the single-cell level during limb development. In addition, we document that the increasing combinatorial expression of Hoxd genes in this developing structure is associated with specific transcriptional signatures and that these signatures illustrate a temporal progression in the differentiation of these cells.

The Insulator Protein CTCF Is Required for Correct Hox Gene Expression, but Not for Embryonic Development in Drosophila

Insulator binding proteins (IBPs) play an important role in regulating gene expression by binding to specific DNA sites to facilitate appropriate gene regulation. There are several IBPs in Drosophila, each defined by their ability to insulate target gene promoters in transgenic assays from the activating or silencing effects of neighboring regulatory elements. Of these, only CCCTC-binding factor (CTCF) has an obvious ortholog in mammals. CTCF is essential for mammalian cell viability and is an important regulator of genome architecture. In flies, CTCF is both maternally deposited and zygotically expressed. Flies lacking zygotic CTCF die as young adults with homeotic defects, suggesting that specific Hox genes are misexpressed in inappropriate body segments. The lack of any major embryonic defects was assumed to be due to the maternal supply of CTCF protein, as maternally contributed factors are often sufficient to progress through much of embryogenesis. Here, we definitively determined the requirement of CTCF for developmental progression in Drosophila We generated animals that completely lack both maternal and zygotic CTCF and found that, contrary to expectation, these mutants progress through embryogenesis and larval life. They develop to pharate adults, which fail to eclose from their pupal case. These mutants show exacerbated homeotic defects compared to zygotic mutants, misexpressing the Hox gene Abdominal-B outside of its normal expression domain early in development. Our results indicate that loss of Drosophila CTCF is not accompanied by widespread effects on gene expression, which may be due to redundant functions with other IBPs. Rather, CTCF is required for correct Hox gene expression patterns and for the viability of adult Drosophila.

CTCF is required for maintenance of auditory hair cells and hearing function in the mouse cochlea

Auditory hair cells play an essential role in hearing. These cells convert sound waves, mechanical stimuli, into electrical signals that are conveyed to the brain via spiral ganglion neurons. The hair cells are located in the organ of Corti within the cochlea. They assemble in a special arrangement with three rows of outer hair cells and one row of inner hair cells. The proper differentiation and preservation of auditory hair cells are essential for acquiring and maintaining hearing function, respectively. Many genetic regulatory mechanisms underlying hair-cell differentiation and maintenance have been elucidated to date. However, the role of epigenetic regulation in hair-cell differentiation and maintenance has not been definitively demonstrated. CTCF is an essential epigenetic component that plays a primary role in the organization of global chromatin architecture. To determine the role of CTCF in mammalian hair cells, we specifically deleted Ctcf in developing hair cells by crossing Ctcf^fl/fl mice with Gfi1^Cre/+ mice. Gfi1^Cre; Ctcf^fl/fl mice did not exhibit obvious developmental defects in hair cells until postnatal day 8. However, at 3 weeks, the Ctcf deficiency caused intermittent degeneration of the stereociliary bundles of outer hair cells, resulting in profound hearing impairment. At 5 weeks, most hair cells were degenerated in Gfi1^Cre; Ctcf^fl/fl mice, and defects in other structures of the organ of Corti, such as the tunnel of Corti and Nuel's space, became apparent. These results suggest that CTCF plays an essential role in maintaining hair cells and hearing function in mammalian cochlea.