Clustered CTCF binding is an evolutionary mechanism to maintain topologically associating domains

Evolutionary patterns and functional effects of 3D chromatin structures in butterflies with extensive genome rearrangements

Chromosome rearrangements may distort 3D chromatin architectures and thus change gene regulation, yet how 3D chromatin structures evolve in insects is largely unknown. Here, we obtain chromosome-level genomes for four butterfly species, Graphium cloanthus, Graphium sarpedon, Graphium eurypylus with 2n = 30, 40, and 60, respectively, and Papilio bianor with 2n = 60. Together with large-scale Hi-C data, we find that inter-chromosome rearrangements very rarely disrupted the pre-existing 3D chromatin structure of ancestral chromosomes. However, some intra-chromosome rearrangements changed 3D chromatin structures compared to the ancestral configuration. We find that new TADs and subTADs have emerged across the rearrangement sites where their adjacent compartments exhibit uniform types. Two intra-chromosome rearrangements altered Rel and lft regulation, potentially contributing to wing patterning differentiation and host plant choice. Notably, butterflies exhibited chromatin loops between Hox gene cluster ANT-C and BX-C, unlike Drosophila. Our CRISPR-Cas9 experiments in butterflies confirm that knocking out the CTCF binding site of the loops in BX-C affected the phenotypes regulated by Antp in ANT-C, resulting in legless larva. Our results reveal evolutionary patterns of insect 3D chromatin structures and provide evidence that 3D chromatin structure changes can play important roles in the evolution of traits.

Epigenetic alterations affecting hematopoietic regulatory networks as drivers of mixed myeloid/lymphoid leukemia

Leukemias with ambiguous lineage comprise several loosely defined entities, often without a clear mechanistic basis. Here, we extensively profile the epigenome and transcriptome of a subgroup of such leukemias with CpG Island Methylator Phenotype. These leukemias exhibit comparable hybrid myeloid/lymphoid epigenetic landscapes, yet heterogeneous genetic alterations, suggesting they are defined by their shared epigenetic profile rather than common genetic lesions. Gene expression enrichment reveals similarity with early T-cell precursor acute lymphoblastic leukemia and a lymphoid progenitor cell of origin. In line with this, integration of differential DNA methylation and gene expression shows widespread silencing of myeloid transcription factors. Moreover, binding sites for hematopoietic transcription factors, including CEBPA, SPI1 and LEF1, are uniquely inaccessible in these leukemias. Hypermethylation also results in loss of CTCF binding, accompanied by changes in chromatin interactions involving key transcription factors. In conclusion, epigenetic dysregulation, and not genetic lesions, explains the mixed phenotype of this group of leukemias with ambiguous lineage. The data collected here constitute a useful and comprehensive epigenomic reference for subsequent studies of acute myeloid leukemias, T-cell acute lymphoblastic leukemias and mixed-phenotype leukemias.

Single-molecule imaging reveals a direct role of CTCF's zinc fingers in SA interaction and cluster-dependent RNA recruitment

CTCF is a zinc finger protein associated with transcription regulation that also acts as a barrier factor for topologically associated domains (TADs) generated by cohesin via loop extrusion. These processes require different properties of CTCF-DNA interaction, and it is still unclear how CTCF's structural features may modulate its diverse roles. Here, we employ single-molecule imaging to study both full-length CTCF and truncation mutants. We show that CTCF enriches at CTCF binding sites (CBSs), displaying a longer lifetime than observed previously. We demonstrate that the zinc finger domains mediate CTCF clustering and that clustering enables RNA recruitment, possibly creating a scaffold for interaction with RNA-binding proteins like cohesin's subunit SA. We further reveal a direct recruitment and an increase of SA residence time by CTCF bound at CBSs, suggesting that CTCF-SA interactions are crucial for cohesin stability on chromatin at TAD borders. Furthermore, we establish a single-molecule T7 transcription assay and show that although a transcribing polymerase can remove CTCF from CBSs, transcription is impaired. Our study shows that context-dependent nucleic acid binding determines the multifaceted CTCF roles in genome organization and transcription regulation.

DNA methylation patterns of transcription factor binding regions characterize their functional and evolutionary contexts

BACKGROUND

DNA methylation is an important epigenetic modification which has numerous roles in modulating genome function. Its levels are spatially correlated across the genome, typically high in repressed regions but low in transcription factor (TF) binding sites and active regulatory regions. However, the mechanisms establishing genome-wide and TF binding site methylation patterns are still unclear.

RESULTS

Here we use a comparative approach to investigate the association of DNA methylation to TF binding evolution in mammals. Specifically, we experimentally profile DNA methylation and combine this with published occupancy profiles of five distinct TFs (CTCF, CEBPA, HNF4A, ONECUT1, FOXA1) in the liver of five mammalian species (human, macaque, mouse, rat, dog). TF binding sites are lowly methylated, but they often also have intermediate methylation levels. Furthermore, biding sites are influenced by the methylation status of CpGs in their wider binding regions even when CpGs are absent from the core binding motif. Employing a classification and clustering approach, we extract distinct and species-conserved patterns of DNA methylation levels at TF binding regions. CEBPA, HNF4A, ONECUT1, and FOXA1 share the same methylation patterns, while CTCF's differ. These patterns characterize alternative functions and chromatin landscapes of TF-bound regions. Leveraging our phylogenetic framework, we find DNA methylation gain upon evolutionary loss of TF occupancy, indicating coordinated evolution. Furthermore, each methylation pattern has its own evolutionary trajectory reflecting its genomic contexts.

CONCLUSIONS

Our epigenomic analyses indicate a role for DNA methylation in TF binding changes across species including that specific DNA methylation profiles characterize TF binding and are associated with their regulatory activity, chromatin contexts, and evolutionary trajectories.

Deciphering the role of a SINE-VNTR-Alu retrotransposon polymorphism as a biomarker of Parkinson's disease progression

SINE-VNTR-Alu (SVA) retrotransposons are transposable elements which represent a source of genetic variation. We previously demonstrated that the presence/absence of a human-specific SVA, termed SVA_67, correlated with the progression of Parkinson's disease (PD). In the present study, we demonstrate that SVA_67 acts as expression quantitative trait loci, thereby exhibiting a strong regulatory effect across the genome using whole genome and transcriptomic data from the Parkinson's progression markers initiative cohort. We further show that SVA_67 is polymorphic for its variable number tandem repeat domain which correlates with both regulatory properties in a luciferase reporter gene assay in vitro and differential expression of multiple genes in vivo. Additionally, this variation's utility as a biomarker is reflected in a correlation with a number of PD progression markers. These experiments highlight the plethora of transcriptomic and phenotypic changes associated with SVA_67 polymorphism which should be considered when investigating the missing heritability of neurodegenerative diseases.

The N-terminal dimerization domains of human and Drosophila CTCF have similar functionality

BACKGROUND

CTCF is highly likely to be the ancestor of proteins that contain large clusters of C2H2 zinc finger domains, and its conservation is observed across most bilaterian organisms. In mammals, CTCF is the primary architectural protein involved in organizing chromosome topology and mediating enhancer-promoter interactions over long distances. In Drosophila, CTCF (dCTCF) cooperates with other architectural proteins to establish long-range interactions and chromatin boundaries. CTCFs of various organisms contain an unstructured N-terminal dimerization domain (DD) and clusters comprising eleven zinc-finger domains of the C2H2 type. The Drosophila (dCTCF) and human (hCTCF) CTCFs share sequence homology in only five C2H2 domains that specifically bind to a conserved 15 bp motif.

RESULTS

Previously, we demonstrated that CTCFs from different organisms carry unstructured N-terminal dimerization domains (DDs) that lack sequence homology. Here we used the CTCFattP(mCh) platform to introduce desired changes in the Drosophila CTCF gene and generated a series of transgenic lines expressing dCTCF with different variants of the N-terminal domain. Our findings revealed that the functionality of dCTCF is significantly affected by the deletion of the N-terminal DD. Additionally, we observed a strong impact on the binding of the dCTCF mutant to chromatin upon deletion of the DD. However, chromatin binding was restored in transgenic flies expressing a chimeric CTCF protein with the DD of hCTCF. Although the chimeric protein exhibited lower expression levels than those of the dCTCF variants, it efficiently bound to chromatin similarly to the wild type (wt) protein.

CONCLUSIONS

Our findings suggest that one of the evolutionarily conserved functions of the unstructured N-terminal dimerization domain is to recruit dCTCF to its genomic sites in vivo.

Implications of the three-dimensional chromatin organization for genome evolution in a fungal plant pathogen

The spatial organization of eukaryotic genomes is linked to their biological functions, although it is not clear how this impacts the overall evolution of a genome. Here, we uncover the three-dimensional (3D) genome organization of the phytopathogen Verticillium dahliae, known to possess distinct genomic regions, designated adaptive genomic regions (AGRs), enriched in transposable elements and genes that mediate host infection. Short-range DNA interactions form clear topologically associating domains (TADs) with gene-rich boundaries that show reduced levels of gene expression and reduced genomic variation. Intriguingly, TADs are less clearly insulated in AGRs than in the core genome. At a global scale, the genome contains bipartite long-range interactions, particularly enriched for AGRs and more generally containing segmental duplications. Notably, the patterns observed for V. dahliae are also present in other Verticillium species. Thus, our analysis links 3D genome organization to evolutionary features conserved throughout the Verticillium genus.

Epigenetic regulation in cancer

Epigenetic modifications are defined as heritable changes in gene activity that do not involve changes in the underlying DNA sequence. The oncogenic process is driven by the accumulation of alterations that impact genome's structure and function. Genetic mutations, which directly disrupt the DNA sequence, are complemented by epigenetic modifications that modulate gene expression, thereby facilitating the acquisition of malignant characteristics. Principals among these epigenetic changes are shifts in DNA methylation and histone mark patterns, which promote tumor development and metastasis. Notably, the reversible nature of epigenetic alterations, as opposed to the permanence of genetic changes, positions the epigenetic machinery as a prime target in the discovery of novel therapeutics. Our review delves into the complexities of epigenetic regulation, exploring its profound effects on tumor initiation, metastatic behavior, metabolic pathways, and the tumor microenvironment. We place a particular emphasis on the dysregulation at each level of epigenetic modulation, including but not limited to, the aberrations in enzymes responsible for DNA methylation and histone modification, subunit loss or fusions in chromatin remodeling complexes, and the disturbances in higher-order chromatin structure. Finally, we also evaluate therapeutic approaches that leverage the growing understanding of chromatin dysregulation, offering new avenues for cancer treatment.

DARDN: A Deep-Learning Approach for CTCF Binding Sequence Classification and Oncogenic Regulatory Feature Discovery

Characterization of gene regulatory mechanisms in cancer is a key task in cancer genomics. CCCTC-binding factor (CTCF), a DNA binding protein, exhibits specific binding patterns in the genome of cancer cells and has a non-canonical function to facilitate oncogenic transcription programs by cooperating with transcription factors bound at flanking distal regions. Identification of DNA sequence features from a broad genomic region that distinguish cancer-specific CTCF binding sites from regular CTCF binding sites can help find oncogenic transcription factors in a cancer type. However, the presence of long DNA sequences without localization information makes it difficult to perform conventional motif analysis. Here, we present DNAResDualNet (DARDN), a computational method that utilizes convolutional neural networks (CNNs) for predicting cancer-specific CTCF binding sites from long DNA sequences and employs DeepLIFT, a method for interpretability of deep learning models that explains the model's output in terms of the contributions of its input features. The method is used for identifying DNA sequence features associated with cancer-specific CTCF binding. Evaluation on DNA sequences associated with CTCF binding sites in T-cell acute lymphoblastic leukemia (T-ALL) and other cancer types demonstrates DARDN's ability in classifying DNA sequences surrounding cancer-specific CTCF binding from control constitutive CTCF binding and identifying sequence motifs for transcription factors potentially active in each specific cancer type. We identify potential oncogenic transcription factors in T-ALL, acute myeloid leukemia (AML), breast cancer (BRCA), colorectal cancer (CRC), lung adenocarcinoma (LUAD), and prostate cancer (PRAD). Our work demonstrates the power of advanced machine learning and feature discovery approach in finding biologically meaningful information from complex high-throughput sequencing data.

TAD evolutionary and functional characterization reveals diversity in mammalian TAD boundary properties and function

Topological associating domains (TADs) are self-interacting genomic units crucial for shaping gene regulation patterns. Despite their importance, the extent of their evolutionary conservation and its functional implications remain largely unknown. In this study, we generate Hi-C and ChIP-seq data and compare TAD organization across four primate and four rodent species and characterize the genetic and epigenetic properties of TAD boundaries in correspondence to their evolutionary conservation. We find 14% of all human TAD boundaries to be shared among all eight species (ultraconserved), while 15% are human-specific. Ultraconserved TAD boundaries have stronger insulation strength, CTCF binding, and enrichment of older retrotransposons compared to species-specific boundaries. CRISPR-Cas9 knockouts of an ultraconserved boundary in a mouse model lead to tissue-specific gene expression changes and morphological phenotypes. Deletion of a human-specific boundary near the autism-related AUTS2 gene results in the upregulation of this gene in neurons. Overall, our study provides pertinent TAD boundary evolutionary conservation annotations and showcases the functional importance of TAD evolution.

Outward-oriented sites within clustered CTCF boundaries are key for intra-TAD chromatin interactions and gene regulation

CTCF plays an important role in 3D genome organization by adjusting the strength of chromatin insulation at TAD boundaries, where clustered CBS (CTCF-binding site) elements are often arranged in a tandem array with a complex divergent or convergent orientation. Here, using Pcdh and HOXD loci as a paradigm, we look into the clustered CTCF TAD boundaries and find that, counterintuitively, outward-oriented CBS elements are crucial for inward enhancer-promoter interactions as well as for gene regulation. Specifically, by combinatorial deletions of a series of putative enhancer elements in mice in vivo or CBS elements in cultured cells in vitro, in conjunction with chromosome conformation capture and RNA-seq analyses, we show that deletions of outward-oriented CBS elements weaken the strength of long-distance intra-TAD promoter-enhancer interactions and enhancer activation of target genes. Our data highlight the crucial role of outward-oriented CBS elements within the clustered CTCF TAD boundaries in developmental gene regulation and have interesting implications on the organization principles of clustered CTCF sites within TAD boundaries.

Three-dimensional chromatin architecture in plants - General features and novelties

Research on the three-dimensional (3D) structure of the genome and its distribution within the nuclear space has made a big leap in the last two decades. Work in the animal field has led to significant advances in our general understanding on eukaryotic genome organization. This did not only bring along insights into how the 3D genome interacts with the epigenetic landscape and the transcriptional machinery but also how 3D genome architecture is relevant for fundamental developmental processes, such as cell differentiation. In parallel, the 3D organization of plant genomes have been extensively studied, which resulted in both congruent and novel findings, contributing to a more complete view on how eukaryotic genomes are organized in multiple dimensions. Plant genomes are remarkably diverse in size, composition, and ploidy. Furthermore, as intrinsically sessile organisms without the possibility to relocate to more favorable environments, plants have evolved an elaborate epigenetic repertoire to rapidly respond to environmental challenges. The diversity in genome organization and the complex epigenetic programs make plants ideal study subjects to acquire a better understanding on universal features and inherent constraints of genome organization. Furthermore, considering a wide range of species allows us to study the evolutionary crosstalk between the various levels of genome architecture. In this article, we aim at summarizing important findings on 3D genome architecture obtained in various plant species. These findings cover many aspects of 3D genome organization on a wide range of levels, from gene loops to topologically associated domains and to global 3D chromosome configurations. We present an overview on plant 3D genome organizational features that resemble those in animals and highlight facets that have only been observed in plants to date.

OSR1 disruption contributes to uterine factor infertility via impaired Müllerian duct development and endometrial receptivity

Three sisters, born from consanguineous parents, manifested a unique Müllerian anomaly characterized by uterine hypoplasia with thin estrogen-unresponsive endometrium and primary amenorrhea, but with spontaneous tubal pregnancies. Through whole-exome sequencing followed by comprehensive genetic analysis, a missense variant was identified in the OSR1 gene. We therefore investigated OSR1/OSR1 expression in postpubertal human uteri, and the prenatal and postnatal expression pattern of Osr1/Osr1 in murine developing Müllerian ducts (MDs) and endometrium, respectively. We then investigated whether Osr1 deletion would affect MD development, using WT and genetically engineered mice. Human uterine OSR1/OSR1 expression was found primarily in the endometrium. Mouse Osr1 was expressed prenatally in MDs and Wolffian ducts (WDs), from rostral to caudal segments, in E13.5 embryos. MDs and WDs were absent on the left side and MDs were rostrally truncated on the right side of E13.5 Osr1-/- embryos. Postnatally, Osr1 was expressed in mouse uteri throughout their lifespan, peaking at postnatal days 14 and 28. Osr1 protein was present primarily in uterine luminal and glandular epithelial cells and in the epithelial cells of mouse oviducts. Through this translational approach, we demonstrated that OSR1 in humans and mice is important for MD development and endometrial receptivity and may be implicated in uterine factor infertility.

In silico discovery of repetitive elements as key sequence determinants of 3D genome folding

Natural and experimental genetic variants can modify DNA loops and insulating boundaries to tune transcription, but it is unknown how sequence perturbations affect chromatin organization genome wide. We developed a deep-learning strategy to quantify the effect of any insertion, deletion, or substitution on chromatin contacts and systematically scored millions of synthetic variants. While most genetic manipulations have little impact, regions with CTCF motifs and active transcription are highly sensitive, as expected. Our unbiased screen and subsequent targeted experiments also point to noncoding RNA genes and several families of repetitive elements as CTCF-motif-free DNA sequences with particularly large effects on nearby chromatin interactions, sometimes exceeding the effects of CTCF sites and explaining interactions that lack CTCF. We anticipate that our disruption tracks may be of broad interest and utility as a measure of 3D genome sensitivity, and our computational strategies may serve as a template for biological inquiry with deep learning.

Transposable elements in mammalian chromatin organization

Transposable elements (TEs) are mobile DNA elements that comprise almost 50% of mammalian genomic sequence. TEs are capable of making additional copies of themselves that integrate into new positions in host genomes. This unique property has had an important impact on mammalian genome evolution and on the regulation of gene expression because TE-derived sequences can function as cis-regulatory elements such as enhancers, promoters and silencers. Now, advances in our ability to identify and characterize TEs have revealed that TE-derived sequences also regulate gene expression by both maintaining and shaping 3D genome architecture. Studies are revealing how TEs contribute raw sequence that can give rise to the structures that shape chromatin organization, and thus gene expression, allowing for species-specific genome innovation and evolutionary novelty.

Regulation of the three-dimensional chromatin organization by transposable elements in pig spleen

Like other mammalian species, the pig genome is abundant with transposable elements (TEs). The importance of TEs for three-dimensional (3D) chromatin organization has been observed in species like human and mouse, yet current understanding about pig TEs is absent. Here, we investigated the contribution of TEs for the 3D chromatin organization in three pig tissues, focusing on spleen which is crucial for both adaptive and innate immunity. We identified dozens of TE families overrepresented with CTCF binding sites, including LTR22_SS, LTR15_SS and LTR16_SSc which are pig-specific families of endogenous retroviruses (ERVs). Interestingly, LTR22_SS elements harbor a CTCF motif and create hundreds of CTCF binding sites that are associated with adaptive immunity. We further applied Hi-C to profile the 3D chromatin structure in spleen and found that TE-derived CTCF binding sites correlate with chromatin insulation and frequently overlap TAD borders and loop anchors. Notably, one LTR22_SS-derived CTCF binding site demarcate a TAD boundary upstream of XCL1, which is a spleen-enriched chemokine gene important for lymphocyte trafficking and inflammation. Overall, this study represents a first step toward understanding the function of TEs on 3D chromatin organization regulation in pigs and expands our understanding about the functional importance of TEs in mammals.

Multi-feature clustering of CTCF binding creates robustness for loop extrusion blocking and Topologically Associating Domain boundaries

Topologically Associating Domains (TADs) separate vertebrate genomes into insulated regulatory neighborhoods that focus genome-associated processes. TADs are formed by Cohesin-mediated loop extrusion, with many TAD boundaries consisting of clustered binding sites of the CTCF insulator protein. Here we determine how this clustering of CTCF binding contributes to the blocking of loop extrusion and the insulation between TADs. We identify enrichment of three features of CTCF binding at strong TAD boundaries, consisting of strongly bound and closely spaced CTCF binding peaks, with a further enrichment of DNA-binding motifs within these peaks. Using multi-contact Nano-C analysis in cells with normal and perturbed CTCF binding, we establish that individual CTCF binding sites contribute to the blocking of loop extrusion, but in an incomplete manner. When clustered, individual CTCF binding sites thus create a stepwise insulation between neighboring TADs. Based on these results, we propose a model whereby multiple instances of temporal loop extrusion blocking create strong insulation between TADs.

Inferring CTCF-binding patterns and anchored loops across human tissues and cell types

CCCTC-binding factor (CTCF) is a transcription regulator with a complex role in gene regulation. The recognition and effects of CTCF on DNA sequences, chromosome barriers, and enhancer blocking are not well understood. Existing computational tools struggle to assess the regulatory potential of CTCF-binding sites and their impact on chromatin loop formation. Here we have developed a deep-learning model, DeepAnchor, to accurately characterize CTCF binding using high-resolution genomic/epigenomic features. This has revealed distinct chromatin and sequence patterns for CTCF-mediated insulation and looping. An optimized implementation of a previous loop model based on DeepAnchor score excels in predicting CTCF-anchored loops. We have established a compendium of CTCF-anchored loops across 52 human tissue/cell types, and this suggests that genomic disruption of these loops could be a general mechanism of disease pathogenesis. These computational models and resources can help investigate how CTCF-mediated cis-regulatory elements shape context-specific gene regulation in cell development and disease progression.

TADs: Dynamic structures to create stable regulatory functions

Mammalian chromosomes are organized at different length scales within the cell nucleus. Topologically Associating Domains (TADs) are structural units of 3D genome organization with functions in gene regulation, DNA replication, recombination and repair. Whereas TADs were initially interpreted as insulated domains, recent studies are revealing that these domains should be interpreted as dynamic collections of actively extruding loops. This process of loop extrusion is subsequently blocked at dedicated TAD boundaries, thereby promoting intra-domain interactions over their surroundings. In this review, we discuss how mammalian TAD structure can emerge from this dynamic process and we discuss recent evidence that TAD boundaries can have regulatory functions.

CTCF and Its Multi-Partner Network for Chromatin Regulation

Architectural proteins are essential epigenetic regulators that play a critical role in organizing chromatin and controlling gene expression. CTCF (CCCTC-binding factor) is a key architectural protein responsible for maintaining the intricate 3D structure of chromatin. Because of its multivalent properties and plasticity to bind various sequences, CTCF is similar to a Swiss knife for genome organization. Despite the importance of this protein, its mechanisms of action are not fully elucidated. It has been hypothesized that its versatility is achieved through interaction with multiple partners, forming a complex network that regulates chromatin folding within the nucleus. In this review, we delve into CTCF's interactions with other molecules involved in epigenetic processes, particularly histone and DNA demethylases, as well as several long non-coding RNAs (lncRNAs) that are able to recruit CTCF. Our review highlights the importance of CTCF partners to shed light on chromatin regulation and pave the way for future exploration of the mechanisms that enable the finely-tuned role of CTCF as a master regulator of chromatin.

DNA architectural protein CTCF facilitates subset-specific chromatin interactions to limit the formation of memory CD8+ T cells

Although the importance of genome organization for transcriptional regulation of cell-fate decisions and function is clear, the changes in chromatin architecture and how these impact effector and memory CD8+ T cell differentiation remain unknown. Using Hi-C, we studied how genome configuration is integrated with CD8+ T cell differentiation during infection and investigated the role of CTCF, a key chromatin remodeler, in modulating CD8+ T cell fates through CTCF knockdown approaches and perturbation of specific CTCF-binding sites. We observed subset-specific changes in chromatin organization and CTCF binding and revealed that weak-affinity CTCF binding promotes terminal differentiation of CD8+ T cells through the regulation of transcriptional programs. Further, patients with de novo CTCF mutations had reduced expression of the terminal-effector genes in peripheral blood lymphocytes. Therefore, in addition to establishing genome architecture, CTCF regulates effector CD8+ T cell heterogeneity through altering interactions that regulate the transcription factor landscape and transcriptome.

Elevated enhancer-oncogene contacts and higher oncogene expression levels by recurrent CTCF inactivating mutations in acute T cell leukemia

Monoallelic inactivation of CCCTC-binding factor (CTCF) in human cancer drives altered methylated genomic states, altered CTCF occupancy at promoter and enhancer regions, and deregulated global gene expression. In patients with T cell acute lymphoblastic leukemia (T-ALL), we find that acquired monoallelic CTCF-inactivating events drive subtle and local genomic effects in nearly half of t(5; 14) (q35; q32.2) rearranged patients, especially when CTCF-binding sites are preserved in between the BCL11B enhancer and the TLX3 oncogene. These solitary intervening sites insulate TLX3 from the enhancer by inducing competitive looping to multiple binding sites near the TLX3 promoter. Reduced CTCF levels or deletion of the intervening CTCF site abrogates enhancer insulation by weakening competitive looping while favoring TLX3 promoter to BCL11B enhancer looping, which elevates oncogene expression levels and leukemia burden.

TAD Evolutionary and functional characterization reveals diversity in mammalian TAD boundary properties and function

Topological associating domains (TADs) are self-interacting genomic units crucial for shaping gene regulation patterns. Despite their importance, the extent of their evolutionary conservation and its functional implications remain largely unknown. In this study, we generate Hi-C and ChIP-seq data and compare TAD organization across four primate and four rodent species, and characterize the genetic and epigenetic properties of TAD boundaries in correspondence to their evolutionary conservation. We find that only 14% of all human TAD boundaries are shared among all eight species (ultraconserved), while 15% are human-specific. Ultraconserved TAD boundaries have stronger insulation strength, CTCF binding, and enrichment of older retrotransposons, compared to species-specific boundaries. CRISPR-Cas9 knockouts of two ultraconserved boundaries in mouse models leads to tissue-specific gene expression changes and morphological phenotypes. Deletion of a human-specific boundary near the autism-related AUTS2 gene results in upregulation of this gene in neurons. Overall, our study provides pertinent TAD boundary evolutionary conservation annotations, and showcase the functional importance of TAD evolution.

The circadian demethylation of a unique intronic deoxymethylCpG-rich island boosts the transcription of its cognate circadian clock output gene

We demonstrate that there is a tight functional relationship between two highly evolutionary conserved cell processes, i.e., the circadian clock (CC) and the circadian DNA demethylation-methylation of cognate deoxyCpG-rich islands. We have discovered that every circadian clock-controlled output gene (CCG), but not the core clock nor its immediate-output genes, contains a single cognate intronic deoxyCpG-rich island, the demethylation-methylation of which is controlled by the CC. During the transcriptional activation period, these intronic islands are demethylated and, upon dimerization of two YY1 protein binding sites located upstream to the transcriptional enhancer and downstream from the deoxyCpG-rich island, store activating components initially assembled on a cognate active enhancer (a RORE, a D-box or an E-box), in keeping with the generation of a transcriptionally active condensate that boosts the initiation of transcription of their cognate pre-mRNAs. We report how these single intronic deoxyCpG-rich islands are instrumental in such a circadian activation/repression transcriptional process.

Widespread contribution of transposable elements to the rewiring of mammalian 3D genomes

Transposable elements (TEs) are major contributors of genetic material in mammalian genomes. These often include binding sites for architectural proteins, including the multifarious master protein, CTCF, which shapes the 3D genome by creating loops, domains, compartment borders, and RNA-DNA interactions. These play a role in the compact packaging of DNA and have the potential to facilitate regulatory function. In this study, we explore the widespread contribution of TEs to mammalian 3D genomes by quantifying the extent to which they give rise to loops and domain border differences across various cell types and species using several 3D genome mapping technologies. We show that specific families and subfamilies of TEs have contributed to lineage-specific 3D chromatin structures across mammalian species. In many cases, these loops may facilitate sustained interaction between distant cis-regulatory elements and target genes, and domains may segregate chromatin state to impact gene expression in a lineage-specific manner. An experimental validation of our analytical findings using CRISPR-Cas9 to delete a candidate TE resulted in disruption of species-specific 3D chromatin structure. Taken together, we comprehensively quantify and selectively validate our finding that TEs contribute to shaping 3D genome organization and may, in some cases, impact gene regulation during the course of mammalian evolution.

Evolution of 3D chromatin organization at different scales

Most animal genomes fold in 3D chromatin domains called topologically associated domains (TADs) that facilitate interactions between cis-regulatory elements (CREs) and promoters. Owing to their critical role in the control of developmental gene expression, we explore how TADs have shaped animal evolution. In the light of recent studies that profile TADs in disparate animal lineages, we discuss their phylogenetic distribution and the mechanisms that underlie their formation. We present evidence indicating that TADs are plastic entities composed of genomic strata of different ages: ancient cores are combined with newer regions and brought into extant TADs through genomic rearrangements. We highlight that newly incorporated TAD strata enable the establishment of new CRE-promoter interactions and in turn new expression patterns that can drive phenotypical innovation. We further highlight how subtle changes in chromatin folding may fine-tune the expression levels of developmental genes and hold a potential for evolutionary significance.

Profiling human pathogenic repeat expansion regions by synergistic and multi-level impacts on molecular connections

Whilst DNA repeat expansions cause numerous heritable human disorders, their origins and underlying pathological mechanisms are often unclear. We collated a dataset comprising 224 human repeat expansions encompassing 203 different genes, and performed a systematic analysis with respect to key topological features at the DNA, RNA and protein levels. Comparison with controls without known pathogenicity and genomic regions lacking repeats, allowed the construction of the first tool to discriminate repeat regions harboring pathogenic repeat expansions (DPREx). At the DNA level, pathogenic repeat expansions exhibited stronger signals for DNA regulatory factors (e.g. H3K4me3, transcription factor-binding sites) in exons, promoters, 5'UTRs and 5'genes but were not significantly different from controls in introns, 3'UTRs and 3'genes. Additionally, pathogenic repeat expansions were also found to be enriched in non-B DNA structures. At the RNA level, pathogenic repeat expansions were characterized by lower free energy for forming RNA secondary structure and were closer to splice sites in introns, exons, promoters and 5'genes than controls. At the protein level, pathogenic repeat expansions exhibited a preference to form coil rather than other types of secondary structure, and tended to encode surface-located protein domains. Guided by these features, DPREx ( http://biomed.nscc-gz.cn/zhaolab/geneprediction/# ) achieved an Area Under the Curve (AUC) value of 0.88 in a test on an independent dataset. Pathogenic repeat expansions are thus located such that they exert a synergistic influence on the gene expression pathway involving inter-molecular connections at the DNA, RNA and protein levels.

Active enhancers strengthen insulation by RNA-mediated CTCF binding at chromatin domain boundaries

Vertebrate genomes are partitioned into chromatin domains or topologically associating domains (TADs), which are typically bound by head-to-head pairs of CTCF binding sites. Transcription at domain boundaries correlates with better insulation; however, it is not known whether the boundary transcripts themselves contribute to boundary function. Here we characterize boundary-associated RNAs genome-wide, focusing on the disease-relevant INK4a/ARF and MYC TAD. Using CTCF site deletions and boundary-associated RNA knockdowns, we observe that boundary-associated RNAs facilitate recruitment and clustering of CTCF at TAD borders. The resulting CTCF enrichment enhances TAD insulation, enhancer-promoter interactions, and TAD gene expression. Importantly, knockdown of boundary-associated RNAs results in loss of boundary insulation function. Using enhancer deletions and CRISPRi of promoters, we show that active TAD enhancers, but not promoters, induce boundary-associated RNA transcription, thus defining a novel class of regulatory enhancer RNAs.

The sound of silence: Transgene silencing in mammalian cell engineering

To elucidate principles operating in native biological systems and to develop novel biotechnologies, synthetic biology aims to build and integrate synthetic gene circuits within native transcriptional networks. The utility of synthetic gene circuits for cell engineering relies on the ability to control the expression of all constituent transgene components. Transgene silencing, defined as the loss of expression over time, persists as an obstacle for engineering primary cells and stem cells with transgenic cargos. In this review, we highlight the challenge that transgene silencing poses to the robust engineering of mammalian cells, outline potential molecular mechanisms of silencing, and present approaches for preventing transgene silencing. We conclude with a perspective identifying future research directions for improving the performance of synthetic gene circuits.

Pattern recognition of topologically associating domains using deep learning

BACKGROUND

Recent increasing evidence indicates that three-dimensional chromosome structure plays an important role in genomic function. Topologically associating domains (TADs) are self-interacting regions that have been shown to be a chromosomal structural unit. During evolution, these are conserved based on checking synteny block cross species. Are there common TAD patterns across species or cell lines?

RESULTS

To address the above question, we propose a novel task-TAD recognition-as opposed to traditional TAD identification. Specifically, we treat Hi-C maps as images, thus re-casting TAD recognition as image pattern recognition, for which we use a convolutional neural network and a residual neural network. In addition, we propose an elegant way to generate non-TAD data for binary classification. We demonstrate deep learning performance which is quite promising, AUC > 0.80, through cross-species and cell-type validation.

CONCLUSIONS

TADs have been shown to be conserved during evolution. Interestingly, our results confirm that the TAD recognition model is practical across species, which indicates that TADs between human and mouse show common patterns from an image classification point of view. Our approach could be a new way to identify TAD variations or patterns among Hi-C maps. For example, TADs of two Hi-C maps are conserved if the two classification models are exchangeable.

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.

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

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

In vivo dissection of a clustered-CTCF domain boundary reveals developmental principles of regulatory insulation

Vertebrate genomes organize into topologically associating domains, delimited by boundaries that insulate regulatory elements from nontarget genes. However, how boundary function is established is not well understood. Here, we combine genome-wide analyses and transgenic mouse assays to dissect the regulatory logic of clustered-CCCTC-binding factor (CTCF) boundaries in vivo, interrogating their function at multiple levels: chromatin interactions, transcription and phenotypes. Individual CTCF binding site (CBS) deletions revealed that the characteristics of specific sites can outweigh other factors such as CBS number and orientation. Combined deletions demonstrated that CBSs cooperate redundantly and provide boundary robustness. We show that divergent CBS signatures are not strictly required for effective insulation and that chromatin loops formed by nonconvergently oriented sites could be mediated by a loop interference mechanism. Further, we observe that insulation strength constitutes a quantitative modulator of gene expression and phenotypes. Our results highlight the modular nature of boundaries and their control over developmental processes.

Association between Triplex-Forming Sites of Cardiac Long Noncoding RNA GATA6-AS1 and Chromatin Organization

This study explored the relationship between 3D genome organization and RNA–DNA triplex-forming sites of long noncoding RNAs (lncRNAs), a group of RNAs that do not code for proteins but are important factors regulating different aspects of genome activity. The triplex-forming sites of anti-sense cardiac lncRNA GATA6-AS1 derived from DBD-Capture-Seq were examined and compared to modular features of 3D genome organization called topologically associated domains (TADs) obtained from Hi-C data. It was found that GATA6-AS1 triplex-forming sites are positioned non-randomly in TADs and their boundaries. The triplex sites showed a preference for TAD boundaries over internal regions of TADs. Computational prediction analysis indicated that CTCF, the key protein involved in TAD specification, may interact with GATA6-AS1, and their binding sites correlate with each other. Examining locations of repeat elements in the genome suggests that the ability of lncRNA GATA6-AS1 to form triplex sites with many genomic locations may be achieved by the rapid expansion of different repeat elements. Some of the triplex-forming sites were found to be positioned in regions that undergo dynamic chromatin organization events such as loss/gain of TAD boundaries during cardiac differentiation. These observed associations suggest that lncRNA–DNA triplex formation may contribute to the specification of TADs in 3D genome organization.

Genome-wide analysis of cis-regulatory changes underlying metabolic adaptation of cavefish

Cis-regulatory changes are key drivers of adaptative evolution. However, their contribution to the metabolic adaptation of organisms is not well understood. Here, we used a unique vertebrate model, Astyanax mexicanus-different morphotypes of which survive in nutrient-rich surface and nutrient-deprived cave waters-to uncover gene regulatory networks underlying metabolic adaptation. We performed genome-wide epigenetic profiling in the liver tissues of Astyanax and found that many of the identified cis-regulatory elements (CREs) have genetically diverged and have differential chromatin features between surface and cave morphotypes, while retaining remarkably similar regulatory signatures between independently derived cave populations. One such CRE in the hpdb gene harbors a genomic deletion in cavefish that abolishes IRF2 repressor binding and derepresses enhancer activity in reporter assays. Selection of this mutation in multiple independent cave populations supports its importance in cave adaptation, and provides novel molecular insights into the evolutionary trade-off between loss of pigmentation and adaptation to food-deprived caves.

Methylation pattern of polymorphically imprinted nc886 is not conserved across mammalia

Background

In humans, the nc886 locus is a polymorphically imprinted metastable epiallele. Periconceptional conditions have an effect on the methylation status of nc886, and further, this methylation status is associated with health outcomes in later life, in line with the Developmental Origins of Health and Disease (DOHaD) hypothesis. Animal models would offer opportunities to study the associations between periconceptional conditions, nc886 methylation status and metabolic phenotypes further. Thus, we set out to investigate the methylation pattern of the nc886 locus in non-human mammals.

Data

We obtained DNA methylation data from the data repository GEO for mammals, whose nc886 gene included all three major parts of nc886 and had sequency similarity of over 80% with the human nc886. Our final sample set consisted of DNA methylation data from humans, chimpanzees, bonobos, gorillas, orangutangs, baboons, macaques, vervets, marmosets and guinea pigs.

Results

In human data sets the methylation pattern of nc886 locus followed the expected bimodal distribution, indicative of polymorphic imprinting. In great apes, we identified a unimodal DNA methylation pattern with 50% methylation level in all individuals and in all subspecies. In Old World monkeys, the between individual variation was greater and methylation on average was close to 60%. In guinea pigs the region around the nc886 homologue was non-methylated. Results obtained from the sequence comparison of the CTCF binding sites flanking the nc886 gene support the results on the DNA methylation data.

Conclusions

Our results indicate that unlike in humans, nc886 is not a polymorphically imprinted metastable epiallele in non-human primates or in guinea pigs, thus implying that animal models are not applicable for nc886 research. The obtained data suggests that the nc886 region may be classically imprinted in great apes, and potentially also in Old World monkeys, but not in guinea pigs.

Regulation of Sox8 through lncRNA Mrhl-Mediated Chromatin Looping in Mouse Spermatogonia

Sox8 is a developmentally important transcription factor that plays an important role in sex maintenance and fertility of adult mice. In the B-type spermatogonial cells, Sox8 is regulated by the long noncoding RNAs (lncRNA) Mrhl in a p68-dependant manner under the control of the Wnt signaling pathway. The downregulation of Mrhl leads to the meiotic commitment of the spermatogonial cells in a Sox8-dependant manner. While the molecular players involved in the regulation of transcription at the Sox8 promoter have been worked out, our current study points to the involvement of the architectural proteins CTCF and cohesin in mediating a chromatin loop that brings the Sox8 promoter in contact with a silencer element present within the gene body in the presence of lncRNA Mrhl concomitant with transcriptional repression. Further, lncRNA Mrhl interacts with the Sox8 locus through the formation of a DNA:DNA:RNA triplex, which is necessary for the recruitment of PRC2 to the locus. The downregulation of lncRNA Mrhl results in the promoter-silencer loop giving way to a promoter-enhancer loop. This active transcription-associated chromatin loop is mediated by YY1 and brings the promoter in contact with the enhancer present downstream of the gene.

Regulation of Il6 expression by single CpG methylation in downstream of Il6 transcription initiation site

The innate immune system is an immediate defense against infectious pathogens by the production of inflammatory cytokines and other mediators. Deficiencies of epigenetic regulatory enzymes, such as Tet1 and Dnmt1, cause dysregulation of cytokine expression. However, it is unclear if DNA methylation at a single CpG dinucleotide in a specific gene locus can regulate gene expression. In this study, we demonstrated that CpG+286 and CpG+348 in exon 2 of the Il6 gene are similar in various primary mouse cells. In lipopolysaccharide-stimulated condition, hypomethylated CpG+286 promoted Il6 expression whereas deletion of CpG+348 led to a reduction in Il6 expression associated with enhanced CTCF binding to the Il6 locus. Moreover, hypomethylation at CpG+286 in alveolar macrophages from aged mice led to higher Il6 expression in response to LPS compared with young mice. Thus, DNA methylation at specific CpG dinucleotides plays an important regulatory role in Il6 expression.

•

TET enzymes regulate DNA methylation at Il6 TIS region and control Il6 expression

•

DNA methylation level at CpG+286 in the Il6 locus directly controls Il6 expression

•

Methylated CpG+348 regulates Il6 expression and CTCF recruitment to the Il6 locus

•

Methylation of CpG+286 in alveolar macrophages controls age-related Il6 expression

Characterisation of the Function of a SINE-VNTR-Alu Retrotransposon to Modulate Isoform Expression at the MAPT Locus

SINE-VNTR-Alu retrotransposons represent one class of transposable elements which contribute to the regulation and evolution of the primate genome and have the potential to be involved in genetic instability and disease progression. However, these polymorphic elements have not been extensively analysed when addressing the missing heritability of neurodegenerative diseases, including Parkinson’s disease (PD) and amyotrophic lateral sclerosis (ALS). SVA_67, a retrotransposon insertion polymorphism, is located in a 1.8 Mb region of high linkage disequilibrium, called the MAPT locus, which is known to contribute to increased risk of developing PD, frontotemporal dementia and other tauopathies. To investigate the role of SVA_67 in directing differential gene expression at this locus, we characterised the impact of SVA_67 allele dosage on isoform expression of several genes in the MAPT locus using the datasets from both the Parkinson’s Progression Markers Initiative and New York Genome Center Consortium Target ALS cohort. The Parkinson’s data was from gene expression in the blood and the ALS data from a variety of CNS regions and allowed us to demonstrate that SVA_67 presence or absence correlated with both isoform- and tissue-specific expression of multiple genes at this locus. This study highlights the importance of addressing SVA polymorphism in disease genetics to gain insight into a better understanding of the role of these regulatory domains to a variety of neurodegenerative diseases.

Dynamic Runx1 chromatin boundaries affect gene expression in hematopoietic development

The transcription factor RUNX1 is a critical regulator of developmental hematopoiesis and is frequently disrupted in leukemia. Runx1 is a large, complex gene that is expressed from two alternative promoters under the spatiotemporal control of multiple hematopoietic enhancers. To dissect the dynamic regulation of Runx1 in hematopoietic development, we analyzed its three-dimensional chromatin conformation in mouse embryonic stem cell (ESC) differentiation cultures. Runx1 resides in a 1.1 Mb topologically associating domain (TAD) demarcated by convergent CTCF motifs. As ESCs differentiate to mesoderm, chromatin accessibility, Runx1 enhancer-promoter (E-P) interactions, and CTCF-CTCF interactions increase in the TAD, along with initiation of Runx1 expression from the P2 promoter. Differentiation to hematopoietic progenitor cells is associated with the formation of tissue-specific sub-TADs over Runx1, a shift in E-P interactions, P1 promoter demethylation, and robust expression from both Runx1 promoters. Deletion of promoter-proximal CTCF sites at the sub-TAD boundaries has no obvious effects on E-P interactions but leads to partial loss of domain structure, mildly affects gene expression, and delays hematopoietic development. Together, our analysis of gene regulation at a large multi-promoter developmental gene reveals that dynamic sub-TAD chromatin boundaries play a role in establishing TAD structure and coordinated gene expression.

Analysis of epigenetic features characteristic of L1 loci expressed in human cells

Only a select few L1 loci in the human genome are expressed in any given cell line or organ, likely to minimize damage done to the genome. The epigenetic features and requirements of expressed L1 loci are currently unknown. Using human cells and comprehensive epigenetic analysis of individual expressed and unexpressed L1 loci, we determined that endogenous L1 transcription depends on a combination of epigenetic factors, including open chromatin, activating histone modifications, and hypomethylation at the L1 promoter. We demonstrate that the L1 promoter seems to require interaction with enhancer elements for optimal function. We utilize epigenetic context to predict the expression status of L1Hs loci that are poorly mappable with RNA-Seq. Our analysis identified a population of ‘transitional’ L1 loci that likely have greater potential to be activated during the epigenetic dysregulation seen in tumors and during aging because they are the most responsive to targeted CRISPR-mediated delivery of trans-activating domains. We demonstrate that an engineered increase in endogenous L1 mRNA expression increases Alu mobilization. Overall, our findings present the first global and comprehensive analysis of epigenetic status of individual L1 loci based on their expression status and demonstrate the importance of epigenetic context for L1 expression heterogeneity.

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.

Mutational bias in spermatogonia impacts the anatomy of regulatory sites in the human genome

Mutation in the germline is the ultimate source of genetic variation, but little is known about the influence of germline chromatin structure on mutational processes. Using ATAC-seq, we profile the open chromatin landscape of human spermatogonia, the most proliferative cell type of the germline, identifying transcription factor binding sites (TFBSs) and PRDM9 binding sites, a subset of which will initiate meiotic recombination. We observe an increase in rare structural variant (SV) breakpoints at PRDM9-bound sites, implicating meiotic recombination in the generation of structural variation. Many germline TFBSs, such as NRF1, are also associated with increased rates of SV breakpoints, apparently independent of recombination. Singleton short insertions (≥5 bp) are highly enriched at TFBSs, particularly at sites bound by testis active TFs, and their rates correlate with those of structural variant breakpoints. Short insertions often duplicate the TFBS motif, leading to clustering of motif sites near regulatory regions in this male-driven evolutionary process. Increased mutation loads at germline TFBSs disproportionately affect neural enhancers with activity in spermatogonia, potentially altering neurodevelopmental regulatory architecture. Local chromatin structure in spermatogonia is thus pervasive in shaping both evolution and disease.

Epigenetic DNA Modifications Upregulate SPRY2 in Human Colorectal Cancers

Conventional wisdom is that Sprouty2 (SPRY2), a suppressor of Receptor Tyrosine Kinase (RTK) signaling, functions as a tumor suppressor and is downregulated in many solid tumors. We reported, for the first time, that increased expression of SPRY2 augments cancer phenotype and Epithelial-Mesenchymal-Transition (EMT) in colorectal cancer (CRC). In this report, we assessed epigenetic DNA modifications that regulate SPRY2 expression in CRC. A total of 4 loci within SPRY2 were evaluated for 5mC using Combined Bisulfite Restriction Analysis (COBRA). Previously sequenced 5hmC nano-hmC seal data within SPRY2 promoter and gene body were evaluated in CRC. Combined bioinformatics analyses of SPRY2 CRC transcripts by RNA-seq/microarray and 450K methyl-array data archived in The Cancer Genome Atlas (TCGA) and GEO database were performed. SPRY2 protein in CRC tumors and cells was measured by Western blotting. Increased SPRY2 mRNA was observed across several CRC datasets and increased protein expression was observed among CRC patient samples. For the first time, SPRY2 hypomethylation was identified in adenocarcinomas in the promoter and gene body. We also revealed, for the first time, increases of 5hmC deposition in the promoter region of SPRY2 in CRC. SPRY2 promoter hypomethylation and increased 5hmC may play an influential role in upregulating SPRY2 in CRC.

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

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

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.

Topologically Associating Domains and Regulatory Landscapes in Development, Evolution and Disease

Animal genomes are folded in topologically associating domains (TADs) that have been linked to the regulation of the genes they contain by constraining regulatory interactions between cis-regulatory elements and promoters. Therefore, TADs are proposed as structural scaffolds for the establishment of regulatory landscapes (RLs). In this review, we discuss recent advances in the connection between TADs and gene regulation, their relationship with gene RLs and their dynamics during development and differentiation. Moreover, we describe how restructuring TADs may lead to pathological conditions, which explains their high evolutionary conservation, but at the same time it provides a substrate for the emergence of evolutionary innovations that lay at the origin of vertebrates and other phylogenetic clades.

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

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

Purifying selection on noncoding deletions of human regulatory loci detected using their cellular pleiotropy

Genomic deletions provide a powerful loss-of-function model in noncoding regions to assess the role of purifying selection on genetic variation. Regulatory element function is characterized by nonuniform tissue and cell type activity, necessarily linking the study of fitness consequences from regulatory variants to their corresponding cellular activity. We generated a callset of deletions from genomes in the Alzheimer's Disease Neuroimaging Initiative (ADNI) and used deletions from The 1000 Genomes Project Consortium (1000GP) in order to examine whether purifying selection preserves noncoding sites of chromatin accessibility marked by DNase I hypersensitivity (DHS), histone modification (enhancer, transcribed, Polycomb-repressed, heterochromatin), and chromatin loop anchors. To examine this in a cellular activity-aware manner, we developed a statistical method, pleiotropy ratio score (PlyRS), which calculates a correlation-adjusted count of "cellular pleiotropy" for each noncoding base pair by analyzing shared regulatory annotations across tissues and cell types. By comparing real deletion PlyRS values to simulations in a length-matched framework and by using genomic covariates in analyses, we found that purifying selection acts to preserve both DHS and enhancer noncoding sites. However, we did not find evidence of purifying selection for noncoding transcribed, Polycomb-repressed, or heterochromatin sites beyond that of the noncoding background. Additionally, we found evidence that purifying selection is acting on chromatin loop integrity by preserving colocalized CTCF binding sites. At regions of DHS, enhancer, and CTCF within chromatin loop anchors, we found evidence that both sites of activity specific to a particular tissue or cell type and sites of cellularly pleiotropic activity are preserved by selection.

LncRNA:DNA triplex-forming sites are positioned at specific areas of genome organization and are predictors for Topologically Associated Domains

BACKGROUND

Chromosomes are organized into units called topologically associated domains (TADs). TADs dictate regulatory landscapes and other DNA-dependent processes. Even though various factors that contribute to the specification of TADs have been proposed, the mechanism is not fully understood. Understanding the process for specification and maintenance of these units is essential in dissecting cellular processes and disease mechanisms.

RESULTS

In this study, we report a genome-wide study that considers the idea of long noncoding RNAs (lncRNAs) mediating chromatin organization using lncRNA:DNA triplex-forming sites (TFSs). By analyzing the TFSs of expressed lncRNAs in multiple cell lines, we find that they are enriched in TADs, their boundaries, and loop anchors. However, they are evenly distributed across different regions of a TAD showing no preference for any specific portions within TADs. No relationship is observed between the locations of these TFSs and CTCF binding sites. However, TFSs are located not just in promoter regions but also in intronic, intergenic, and 3'UTR regions. We also show these triplex-forming sites can be used as predictors in machine learning models to discriminate TADs from other genomic regions. Finally, we compile a list of important "TAD-lncRNAs" which are top predictors for TADs identification.

CONCLUSIONS

Our observations advocate the idea that lncRNA:DNA TFSs are positioned at specific areas of the genome organization and are important predictors for TADs. LncRNA:DNA triplex formation most likely is a general mechanism of action exhibited by some lncRNAs, not just for direct gene regulation but also to mediate 3D chromatin organization.