Co-binding by YY1 identifies the transcriptionally active, highly conserved set of CTCF-bound regions in primate genomes

Epigenetic Regulation of Mammalian Cardiomyocyte Development

The heart is the first organ formed during mammalian development and functions to distribute nutrients and oxygen to other parts of the developing embryo. Cardiomyocytes are the major cell types of the heart and provide both structural support and contractile function to the heart. The successful differentiation of cardiomyocytes during early development is under tight regulation by physical and molecular factors. We have reviewed current studies on epigenetic factors critical for cardiomyocyte differentiation, including DNA methylation, histone modifications, chromatin remodelers, and noncoding RNAs. This review also provides comprehensive details on structural and morphological changes associated with the differentiation of fetal and postnatal cardiomyocytes and highlights their differences. A holistic understanding of all aspects of cardiomyocyte development is critical for the successful in vitro differentiation of cardiomyocytes for therapeutic purposes.

BERT-TFBS: a novel BERT-based model for predicting transcription factor binding sites by transfer learning

Transcription factors (TFs) are proteins essential for regulating genetic transcriptions by binding to transcription factor binding sites (TFBSs) in DNA sequences. Accurate predictions of TFBSs can contribute to the design and construction of metabolic regulatory systems based on TFs. Although various deep-learning algorithms have been developed for predicting TFBSs, the prediction performance needs to be improved. This paper proposes a bidirectional encoder representations from transformers (BERT)-based model, called BERT-TFBS, to predict TFBSs solely based on DNA sequences. The model consists of a pre-trained BERT module (DNABERT-2), a convolutional neural network (CNN) module, a convolutional block attention module (CBAM) and an output module. The BERT-TFBS model utilizes the pre-trained DNABERT-2 module to acquire the complex long-term dependencies in DNA sequences through a transfer learning approach, and applies the CNN module and the CBAM to extract high-order local features. The proposed model is trained and tested based on 165 ENCODE ChIP-seq datasets. We conducted experiments with model variants, cross-cell-line validations and comparisons with other models. The experimental results demonstrate the effectiveness and generalization capability of BERT-TFBS in predicting TFBSs, and they show that the proposed model outperforms other deep-learning models. The source code for BERT-TFBS is available at https://github.com/ZX1998-12/BERT-TFBS.

The N-terminus of Stag1 is required to repress the 2C program by maintaining rRNA expression and nucleolar integrity

Our understanding of how STAG proteins contribute to cell identity and disease have largely been studied from the perspective of chromosome topology and protein-coding gene expression. Here, we show that STAG1 is the dominant paralog in mouse embryonic stem cells (mESCs) and is required for pluripotency. mESCs express a wide diversity of naturally occurring Stag1 isoforms, resulting in complex regulation of both the levels of STAG paralogs and the proportion of their unique terminal ends. Skewing the balance of these isoforms impacts cell identity. We define a novel role for STAG1, in particular its N-terminus, in regulating repeat expression, nucleolar integrity, and repression of the two-cell (2C) state to maintain mESC identity. Our results move beyond protein-coding gene regulation via chromatin loops to new roles for STAG1 in nucleolar structure and function, and offer fresh perspectives on how STAG proteins, known to be cancer targets, contribute to cell identity and disease.

Sequence-based machine learning reveals 3D genome differences between bonobos and chimpanzees

Phenotypic divergence between closely related species, including bonobos and chimpanzees (genus Pan), is largely driven by variation in gene regulation. The 3D structure of the genome mediates gene expression; however, genome folding differences in Pan are not well understood. Here, we apply machine learning to predict genome-wide 3D genome contact maps from DNA sequence for 56 bonobos and chimpanzees, encompassing all five extant lineages. We use a pairwise approach to estimate 3D divergence between individuals from the resulting contact maps in 4,420 1 Mb genomic windows. While most pairs were similar, ∼17% were predicted to be substantially divergent in genome folding. The most dissimilar maps were largely driven by single individuals with rare variants that produce unique 3D genome folding in a region. We also identified 89 genomic windows where bonobo and chimpanzee contact maps substantially diverged, including several windows harboring genes associated with traits implicated in Pan phenotypic divergence. We used in silico mutagenesis to identify 51 3D-modifying variants in these bonobo-chimpanzee divergent windows, finding that 34 or 66.67% induce genome folding changes via CTCF binding motif disruption. Our results reveal 3D genome variation at the population-level and identify genomic regions where changes in 3D folding may contribute to phenotypic differences in our closest living relatives.

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.

Heat shock induces alternative polyadenylation through dynamic DNA methylation-regulated chromatin looping

Alternative cleavage and polyadenylation (APA) is a gene regulatory mechanism used by cells under stress to upregulate proteostasis-promoting transcripts, but how cells achieve this remains poorly understood. Previously, we elucidated a DNA methylation-regulated APA mechanism, in which gene body DNA methylation enhances distal poly(A) isoform expression by blocking CTCF binding and chromatin loop formation at APA control regions. We hypothesized that DNA methylation-regulated APA is one mechanism cells employ to induce proteostasis-promoting poly(A) isoforms. At the DNAJB6 co-chaperone gene locus, acute heat shock resulted in binding of stress response transcription factors HSF1, ATF6, and YY1 at the APA control region and an increase in the expression of the proximal poly(A) isoform known to prevent protein aggregation. Furthermore, TET1 was recruited to rapidly demethylate DNA, facilitating CTCF binding and chromatin loop formation, thereby reinforcing preferential proximal poly(A) isoform expression. As cells recovered, the transcription factors vacated the APA control region, and DNMT1 was recruited to remethylate the region. This process resolved chromatin looping and reset the poly(A) isoform expression pattern. Our findings unveil an epigenetic mechanism enabling cells to dynamically modulate poly(A) isoforms in response to stress while shedding light on the interplay between DNA methylation, transcription factors, and chromatin looping.

Metabolic regulation of CTCF expression and chromatin association dictates starvation response in mice and flies

Coordinated temporal control of gene expression is essential for physiological homeostasis, especially during metabolic transitions. However, the interplay between chromatin architectural proteins and metabolism in regulating transcription is less understood. Here, we demonstrate a conserved bidirectional interplay between CTCF (CCCTC-binding factor) expression/function and metabolic inputs during feed-fast cycles. Our results indicate that its loci-specific functional diversity is associated with physiological plasticity in mouse hepatocytes. CTCF differential expression and long non-coding RNA-Jpx mediated changes in chromatin occupancy, unraveled its paradoxical yet tuneable functions, which are governed by metabolic inputs. We illustrate the key role of CTCF in controlling temporal cascade of transcriptional response, with effects on hepatic mitochondrial energetics and lipidome. Underscoring the evolutionary conservation of CTCF-dependent metabolic homeostasis, CTCF knockdown in flies abrogated starvation resistance. In summary, we demonstrate the interplay between CTCF and metabolic inputs that highlights the coupled plasticity of physiological responses and chromatin function.

Targeting Transcription Factor YY1 for Cancer Treatment: Current Strategies and Future Directions

Cancer represents a significant and persistent global health burden, with its impact underscored by its prevalence and devastating consequences. Whereas numerous oncogenes could contribute to cancer development, a group of transcription factors (TFs) are overactive in the majority of tumors. Targeting these TFs may also combat the downstream oncogenes activated by the TFs, making them attractive potential targets for effective antitumor therapeutic strategy. One such TF is yin yang 1 (YY1), which plays crucial roles in the development and progression of various tumors. In preclinical studies, YY1 inhibition has shown efficacy in inhibiting tumor growth, promoting apoptosis, and sensitizing tumor cells to chemotherapy. Recent studies have also revealed the potential of combining YY1 inhibition with immunotherapy for enhanced antitumor effects. However, clinical translation of YY1-targeted therapy still faces challenges in drug specificity and delivery. This review provides an overview of YY1 biology, its role in tumor development and progression, as well as the strategies explored for YY1-targeted therapy, with a focus on their clinical implications, including those using small molecule inhibitors, RNA interference, and gene editing techniques. Finally, we discuss the challenges and current limitations of targeting YY1 and the need for further research in this area.

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.

Characterizing collaborative transcription regulation with a graph-based deep learning approach

Human epigenome and transcription activities have been characterized by a number of sequence-based deep learning approaches which only utilize the DNA sequences. However, transcription factors interact with each other, and their collaborative regulatory activities go beyond the linear DNA sequence. Therefore leveraging the informative 3D chromatin organization to investigate the collaborations among transcription factors is critical. We developed ECHO, a graph-based neural network, to predict chromatin features and characterize the collaboration among them by incorporating 3D chromatin organization from 200-bp high-resolution Micro-C contact maps. ECHO predicted 2,583 chromatin features with significantly higher average AUROC and AUPR than the best sequence-based model. We observed that chromatin contacts of different distances affected different types of chromatin features’ prediction in diverse ways, suggesting complex and divergent collaborative regulatory mechanisms. Moreover, ECHO was interpretable via gradient-based attribution methods. The attributions on chromatin contacts identify important contacts relevant to chromatin features. The attributions on DNA sequences identify TF binding motifs and TF collaborative binding. Furthermore, combining the attributions on contacts and sequences reveals important sequence patterns in the neighborhood which are relevant to a target sequence’s chromatin feature prediction.

Inferring mammalian tissue-specific regulatory conservation by predicting tissue-specific differences in open chromatin

BACKGROUND

Evolutionary conservation is an invaluable tool for inferring functional significance in the genome, including regions that are crucial across many species and those that have undergone convergent evolution. Computational methods to test for sequence conservation are dominated by algorithms that examine the ability of one or more nucleotides to align across large evolutionary distances. While these nucleotide alignment-based approaches have proven powerful for protein-coding genes and some non-coding elements, they fail to capture conservation of many enhancers, distal regulatory elements that control spatial and temporal patterns of gene expression. The function of enhancers is governed by a complex, often tissue- and cell type-specific code that links combinations of transcription factor binding sites and other regulation-related sequence patterns to regulatory activity. Thus, function of orthologous enhancer regions can be conserved across large evolutionary distances, even when nucleotide turnover is high.

RESULTS

We present a new machine learning-based approach for evaluating enhancer conservation that leverages the combinatorial sequence code of enhancer activity rather than relying on the alignment of individual nucleotides. We first train a convolutional neural network model that can predict tissue-specific open chromatin, a proxy for enhancer activity, across mammals. Next, we apply that model to distinguish instances where the genome sequence would predict conserved function versus a loss of regulatory activity in that tissue. We present criteria for systematically evaluating model performance for this task and use them to demonstrate that our models accurately predict tissue-specific conservation and divergence in open chromatin between primate and rodent species, vastly out-performing leading nucleotide alignment-based approaches. We then apply our models to predict open chromatin at orthologs of brain and liver open chromatin regions across hundreds of mammals and find that brain enhancers associated with neuron activity have a stronger tendency than the general population to have predicted lineage-specific open chromatin.

CONCLUSION

The framework presented here provides a mechanism to annotate tissue-specific regulatory function across hundreds of genomes and to study enhancer evolution using predicted regulatory differences rather than nucleotide-level conservation measurements.

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.

OUP accepted manuscript

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 D₃ (1,25(OH)₂D₃). We performed CTCF HiChIP assays to produce the first genome-wide dataset of CTCF long-range interactions in 1,25(OH)₂D₃-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.

Variants on the UBE2L3/YDJC Autoimmune Disease Risk Haplotype Increase UBE2L3 Expression by Modulating CCCTC-Binding Factor and YY1 Binding

OBJECTIVE

Genetic variants spanning UBE2L3 are associated with increased expression of the UBE2L3-encoded E2 ubiquitin-conjugating enzyme H7 (UbcH7), which facilitates activation of proinflammatory NF-κB signaling and susceptibility to autoimmune diseases. We undertook this study to delineate how genetic variants carried on the UBE2L3/YDJC autoimmune risk haplotype function to drive hypermorphic UBE2L3 expression.

METHODS

We used bioinformatic analyses, electrophoretic mobility shift assays, and luciferase reporter assays to identify and functionally characterize allele-specific effects of risk variants positioned in chromatin accessible regions of immune cells. Chromatin conformation capture with quantitative polymerase chain reaction (3C-qPCR), chromatin immunoprecipitation (ChIP)-qPCR, and small interfering RNA (siRNA) knockdown assays were performed on patient-derived Epstein-Barr virus-transformed B cells homozygous for the UBE2L3/YDJC nonrisk or risk haplotype to determine if the risk haplotype increases UBE2L3 expression by altering the regulatory chromatin architecture in the region.

RESULTS

Of the 7 prioritized variants, 5 demonstrated allele-specific increases in nuclear protein binding affinity and regulatory activity. High-throughput sequencing of chromosome conformation capture coupled with ChIP (HiChIP) and 3C-qPCR uncovered a long-range interaction between the UBE2L3 promoter (rs140490, rs140491, rs11089620) and the downstream YDJC promoter (rs3747093) that was strengthened in the presence of the UBE2L3/YDJC risk haplotype, and correlated with the loss of CCCTC-binding factor (CTCF) and gain of YY1 binding at the risk alleles. Depleting YY1 by siRNA disrupted the long-range interaction between the 2 promoters and reduced UBE2L3 expression.

CONCLUSION

The UBE2L3/YDJC autoimmune risk haplotype increases UBE2L3 expression through strengthening a YY1-mediated interaction between the UBE2L3 and YDJC promoters.

Emerging roles for endogenous retroviruses in immune epigenetic regulation

In recent years, there has been significant progress toward understanding the transcriptional networks underlying mammalian immune responses, fueled by advances in regulatory genomic technologies. Epigenomic studies profiling immune cells have generated detailed genome-wide maps of regulatory elements that will be key to deciphering the regulatory networks underlying cellular immune responses and autoimmune disorders. Unbiased analyses of these genomic maps have uncovered endogenous retroviruses as an unexpected ally in the regulation of human immune systems. Despite their parasitic origins, studies are finding an increasing number of examples of retroviral sequences having been co-opted for beneficial immune function and regulation by the host cell. Here, we review how endogenous retroviruses have given rise to numerous regulatory elements that shape the epigenetic landscape of host immune responses. We will discuss the implications of these elements on the function, dysfunction, and evolution of innate immunity.

The genetic architecture of DNA replication timing in human pluripotent stem cells

DNA replication follows a strict spatiotemporal program that intersects with chromatin structure but has a poorly understood genetic basis. To systematically identify genetic regulators of replication timing, we exploited inter-individual variation in human pluripotent stem cells from 349 individuals. We show that the human genome's replication program is broadly encoded in DNA and identify 1,617 cis-acting replication timing quantitative trait loci (rtQTLs) - sequence determinants of replication initiation. rtQTLs function individually, or in combinations of proximal and distal regulators, and are enriched at sites of histone H3 trimethylation of lysines 4, 9, and 36 together with histone hyperacetylation. H3 trimethylation marks are individually repressive yet synergistically associate with early replication. We identify pluripotency-related transcription factors and boundary elements as positive and negative regulators of replication timing, respectively. Taken together, human replication timing is controlled by a multi-layered mechanism with dozens of effectors working combinatorially and following principles analogous to transcription regulation.

The emerging role of chromatin remodelers in neurodevelopmental disorders: a developmental perspective

Neurodevelopmental disorders (NDDs), including intellectual disability (ID) and autism spectrum disorders (ASD), are a large group of disorders in which early insults during brain development result in a wide and heterogeneous spectrum of clinical diagnoses. Mutations in genes coding for chromatin remodelers are overrepresented in NDD cohorts, pointing towards epigenetics as a convergent pathogenic pathway between these disorders. In this review we detail the role of NDD-associated chromatin remodelers during the developmental continuum of progenitor expansion, differentiation, cell-type specification, migration and maturation. We discuss how defects in chromatin remodelling during these early developmental time points compound over time and result in impaired brain circuit establishment. In particular, we focus on their role in the three largest cell populations: glutamatergic neurons, GABAergic neurons, and glia cells. An in-depth understanding of the spatiotemporal role of chromatin remodelers during neurodevelopment can contribute to the identification of molecular targets for treatment strategies.

CTCF As an Example of DNA-Binding Transcription Factors Containing Clusters of C2H2-Type Zinc Fingers

In mammals, most of the boundaries of topologically associating domains and all well-studied insulators are rich in binding sites for the CTCF protein. According to existing experimental data, CTCF is a key factor in the organization of the architecture of mammalian chromosomes. A characteristic feature of the CTCF is that the central part of the protein contains a cluster consisting of eleven domains of C2H2-type zinc fingers, five of which specifically bind to a long DNA sequence conserved in most animals. The class of transcription factors that carry a cluster of C2H2-type zinc fingers consisting of five or more domains (C2H2 proteins) is widely represented in all groups of animals. The functions of most C2H2 proteins still remain unknown. This review presents data on the structure and possible functions of these proteins, using the example of the vertebrate CTCF protein and several well- characterized C2H2 proteins in Drosophila and mammals.

The Why of YY1: Mechanisms of Transcriptional Regulation by Yin Yang 1

First described in 1991, Yin Yang 1 (YY1) is a transcription factor that is ubiquitously expressed throughout mammalian cells. It regulates both transcriptional activation and repression, in a seemingly context-dependent manner. YY1 has a well-established role in the development of the central nervous system, where it is involved in neurogenesis and maintenance of homeostasis in the developing brain. In neurodevelopmental and neurodegenerative disease, the crucial role of YY1 in cellular processes in the central nervous system is further underscored. In this mini-review, we discuss the various mechanisms leading to the transcriptional activating and repressing roles of YY1, including its role as a traditional transcription factor, its interactions with cofactors and chromatin modifiers, the role of YY1 in the non-coding genome and 3D chromatin organization and the possible implications of the phase-separation mechanism on YY1 function. We provide examples on how these processes can be involved in normal development and how alterations can lead to various diseases.

Yin Yang 1 is a potent activator of human T lymphotropic virus type 1 LTR-driven gene expression via RNA binding

Yin Yang 1 (YY1) is a DNA-binding transcription factor that either activates or represses gene expression. YY1 has previously been implicated in the transcriptional silencing of many retroviruses by binding to DNA sequences in the U3 region of the viral long terminal repeat (LTR). We here show that YY1 overexpression leads to profound activation, rather than repression, of human T lymphotropic virus type 1 (HTLV-1) expression, while YY1 down-regulation reduces HTLV-1 expression. The YY1 responsive element mapped not to YY1 DNA-binding sites in the HTLV-1 LTR but to the R region. The HTLV-1 R sequence alone is sufficient to provide YY1 responsiveness to a nonresponsive promoter, but only in the sense orientation and only when included as part of the mRNA. YY1 binds to the R region of HTLV-1 RNA in vitro and in vivo, leading to increased transcription initiation and elongation. The findings indicate that YY1 is a potent transactivator of HTLV-1 gene expression acting via binding viral RNA, rather than DNA.

Impact of 3D genome organization, guided by cohesin and CTCF looping, on sex-biased chromatin interactions and gene expression in mouse liver

Several thousand sex-differential distal enhancers have been identified in mouse liver; however, their links to sex-biased genes and the impact of any sex-differences in nuclear organization and chromatin interactions are unknown. To address these issues, we first characterized 1847 mouse liver genomic regions showing significant sex differential occupancy by cohesin and CTCF, two key 3D nuclear organizing factors. These sex-differential binding sites were primarily distal to sex-biased genes but rarely generated sex-differential TAD (topologically associating domain) or intra-TAD loop anchors, and were sometimes found in TADs without sex-biased genes. A substantial subset of sex-biased cohesin-non-CTCF binding sites, but not sex-biased cohesin-and-CTCF binding sites, overlapped sex-biased enhancers. Cohesin depletion reduced the expression of male-biased genes with distal, but not proximal, sex-biased enhancers by >10-fold, implicating cohesin in long-range enhancer interactions regulating sex-biased genes. Using circularized chromosome conformation capture-based sequencing (4C-seq), we showed that sex differences in distal sex-biased enhancer-promoter interactions are common. Intra-TAD loops with sex-independent cohesin-and-CTCF anchors conferred sex specificity to chromatin interactions indirectly, by insulating sex-biased enhancer-promoter contacts and by bringing sex-biased genes into closer proximity to sex-biased enhancers. Furthermore, sex-differential chromatin interactions involving sex-biased gene promoters, enhancers, and lncRNAs were associated with sex-biased binding of cohesin and/or CTCF. These studies elucidate how 3D genome organization impacts sex-biased gene expression in a non-reproductive tissue through both direct and indirect effects of cohesin and CTCF looping on distal enhancer interactions with sex-differentially expressed genes.

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

BACKGROUND

CTCF binding contributes to the establishment of a higher-order genome structure by demarcating the boundaries of large-scale topologically associating domains (TADs). However, despite the importance and conservation of TADs, the role of CTCF binding in their evolution and stability remains elusive.

RESULTS

We carry out an experimental and computational study that exploits the natural genetic variation across five closely related species to assess how CTCF binding patterns stably fixed by evolution in each species contribute to the establishment and evolutionary dynamics of TAD boundaries. We perform CTCF ChIP-seq in multiple mouse species to create genome-wide binding profiles and associate them with TAD boundaries. Our analyses reveal that CTCF binding is maintained at TAD boundaries by a balance of selective constraints and dynamic evolutionary processes. Regardless of their conservation across species, CTCF binding sites at TAD boundaries are subject to stronger sequence and functional constraints compared to other CTCF sites. TAD boundaries frequently harbor dynamically evolving clusters containing both evolutionarily old and young CTCF sites as a result of the repeated acquisition of new species-specific sites close to conserved ones. The overwhelming majority of clustered CTCF sites colocalize with cohesin and are significantly closer to gene transcription start sites than nonclustered CTCF sites, suggesting that CTCF clusters particularly contribute to cohesin stabilization and transcriptional regulation.

CONCLUSIONS

Dynamic conservation of CTCF site clusters is an apparently important feature of CTCF binding evolution that is critical to the functional stability of a higher-order chromatin structure.

AIControl: replacing matched control experiments with machine learning improves ChIP-seq peak identification

ChIP-seq is a technique to determine binding locations of transcription factors, which remains a central challenge in molecular biology. Current practice is to use a 'control' dataset to remove background signals from a immunoprecipitation (IP) 'target' dataset. We introduce the AIControl framework, which eliminates the need to obtain a control dataset and instead identifies binding peaks by estimating the distributions of background signals from many publicly available control ChIP-seq datasets. We thereby avoid the cost of running control experiments while simultaneously increasing the accuracy of binding location identification. Specifically, AIControl can (i) estimate background signals at fine resolution, (ii) systematically weigh the most appropriate control datasets in a data-driven way, (iii) capture sources of potential biases that may be missed by one control dataset and (iv) remove the need for costly and time-consuming control experiments. We applied AIControl to 410 IP datasets in the ENCODE ChIP-seq database, using 440 control datasets from 107 cell types to impute background signal. Without using matched control datasets, AIControl identified peaks that were more enriched for putative binding sites than those identified by other popular peak callers that used a matched control dataset. We also demonstrated that our framework identifies binding sites that recover documented protein interactions more accurately.

p63 cooperates with CTCF to modulate chromatin architecture in skin keratinocytes

The transcription factor p63 regulates epidermal genes and the enhancer landscape in skin keratinocytes. Its molecular function in controlling the chromatin structure is, however, not yet completely understood. Here, we integrated multi-omics profiles, including the transcriptome, transcription factor DNA-binding and chromatin accessibility, in skin keratinocytes isolated from EEC syndrome patients carrying p63 mutations, to examine the role of p63 in shaping the chromatin architecture. We found decreased chromatin accessibility in p63- and CTCF-bound open chromatin regions that potentially contributed to gene deregulation in mutant keratinocytes. Cooperation of p63 and CTCF seemed to assist chromatin interactions between p63-bound enhancers and gene promoters in skin keratinocytes. Our study suggests an intriguing model where cell type-specific transcription factors such as p63 cooperate with the genome organizer CTCF in the three-dimensional chromatin space to regulate the transcription program important for the proper cell identity.

Identification of Cancer Drivers at CTCF Insulators in 1,962 Whole Genomes

Recent studies have shown that mutations at non-coding elements, such as promoters and enhancers, can act as cancer drivers. However, an important class of non-coding elements, namely CTCF insulators, has been overlooked in the previous driver analyses. We used insulator annotations from CTCF and cohesin ChIA-PET and analyzed somatic mutations in 1,962 whole genomes from 21 cancer types. Using the heterogeneous patterns of transcription-factor-motif disruption, functional impact, and recurrence of mutations, we developed a computational method that revealed 21 insulators showing signals of positive selection. In particular, mutations in an insulator in multiple cancer types, including 16% of melanoma samples, are associated with TGFB1 up-regulation. Using CRISPR-Cas9, we find that alterations at two of the most frequently mutated regions in this insulator increase cell growth by 40%-50%, supporting the role of this boundary element as a cancer driver. Thus, our study reveals several CTCF insulators as putative cancer drivers.

YY1 negatively regulates the XAF1 gene expression in prostate cancer

XAF1 is a tumor suppressor gene with low or absent expression in cancer. Since transcriptional reactivation or ectopic-mediated expression of XAF1 inhibits tumor growth, it is of great interest to elucidate the molecular mechanisms leading to XAF1 silencing. YY1 is a transcription factor that acts as a repressor or an activator to modulate several cancer-associated cellular processes. Both YY1 and XAF1 have key roles in prostate cancer (PCa) progression and are associated with worse clinical outcomes. To assess whether YY1 regulates the transcriptional activation of the XAF1 gene, we performed gene-reporter assays coupled with site-directed mutagenesis, which showed that YY1 is able to mediate XAF1 silencing. Concordantly, ChIP-qPCR assays showed that YY1 interacts with the XAF1 promoter in PC3 cells that lacks XAF1 expression. This association was lost after exposure to epigenetic modulators that induce XAF1 expression. Further supporting the YY1's repressive role, we found transcriptional reactivation of the XAF1 gene by YY1 downregulation. As expected by previous reports showing that HDAC1 is needed for YY1-mediated repressive actions, we observed XAF1 re-expression after either inhibition or downregulation of the HDAC1 gene. Finally, expression data retrieved from the TCGA consortium showed that PCa samples presented lower XAF1 and higher HDAC expression levels than normal tissues. Thus, our results support a model in which YY1 is able to silence tumor suppressor genes such as XAF1 through HDAC1 in PCa.

Disruption of CTCF-YY1-dependent looping of the human papillomavirus genome activates differentiation-induced viral oncogene transcription

The complex life cycle of oncogenic human papillomavirus (HPV) initiates in undifferentiated basal epithelial keratinocytes where expression of the E6 and E7 oncogenes is restricted. Upon epithelial differentiation, E6/E7 transcription is increased through unknown mechanisms to drive cellular proliferation required to support virus replication. We report that the chromatin-organising CCCTC-binding factor (CTCF) promotes the formation of a chromatin loop in the HPV genome that epigenetically represses viral enhancer activity controlling E6/E7 expression. CTCF-dependent looping is dependent on the expression of the CTCF-associated Yin Yang 1 (YY1) transcription factor and polycomb repressor complex (PRC) recruitment, resulting in trimethylation of histone H3 at lysine 27. We show that viral oncogene up-regulation during cellular differentiation results from YY1 down-regulation, disruption of viral genome looping, and a loss of epigenetic repression of viral enhancer activity. Our data therefore reveal a key role for CTCF-YY1-dependent looping in the HPV life cycle and identify a regulatory mechanism that could be disrupted in HPV carcinogenesis.

Characterization of the ICCE Repeat in Mammals Reveals an Evolutionary Relationship with the DXZ4 Macrosatellite through Conserved CTCF Binding Motifs

Appreciation is growing for how chromosomes are organized in three-dimensional space at interphase. Microscopic and high throughput sequence-based studies have established that the mammalian inactive X chromosome (Xi) adopts an alternate conformation relative to the active X chromosome. The Xi is organized into several multi-megabase chromatin loops called superloops. At the base of these loops are superloop anchors, and in humans three of these anchors are composed of large tandem repeat DNA that include DXZ4, Functional Intergenic Repeating RNA Element, and Inactive-X CTCF-binding Contact Element (ICCE). Each repeat contains a high density of binding sites for the architectural organization protein CCCTC-binding factor (CTCF) which exclusively associates with the Xi allele in normal cells. Removal of DXZ4 from the Xi compromises proper folding of the chromosome. In this study, we report the characterization of the ICCE tandem repeat, for which very little is known. ICCE is embedded within an intron of the Nobody (NBDY) gene locus at Xp11.21. We find that primary DNA sequence conservation of ICCE is only retained in higher primates, but that ICCE orthologs exist beyond the primate lineage. Like DXZ4, what is conserved is organization of the underlying DNA into a large tandem repeat, physical location within the NBDY locus and conservation of short DNA sequences corresponding to specific CTCF and Yin Yang 1 binding motifs that correlate with female-specific DNA hypomethylation. Unlike DXZ4, ICCE is not common to all eutherian mammals. Analysis of certain ICCE CTCF motifs reveal striking similarity with the DXZ4 motif and support an evolutionary relationship between DXZ4 and ICCE.

Species-Specific Changes in a Primate Transcription Factor Network Provide Insights into the Molecular Evolution of the Primate Prefrontal Cortex

The human prefrontal cortex (PFC) differs from that of other primates with respect to size, histology, and functional abilities. Here, we analyzed genome-wide expression data of humans, chimpanzees, and rhesus macaques to discover evolutionary changes in transcription factor (TF) networks that may underlie these phenotypic differences. We determined the co-expression networks of all TFs with species-specific expression including their potential target genes and interaction partners in the PFC of all three species. Integrating these networks allowed us inferring an ancestral network for all three species. This ancestral network as well as the networks for each species is enriched for genes involved in forebrain development, axonogenesis, and synaptic transmission. Our analysis allows us to directly compare the networks of each species to determine which links have been gained or lost during evolution. Interestingly, we detected that most links were gained on the human lineage, indicating increase TF cooperativity in humans. By comparing network changes between different tissues, we discovered that in brain tissues, but not in the other tissues, the human networks always had the highest connectivity. To pinpoint molecular changes underlying species-specific phenotypes, we analyzed the sub-networks of TFs derived only from genes with species-specific expression changes in the PFC. These sub-networks differed significantly in structure and function between the human and chimpanzee. For example, the human-specific sub-network is enriched for TFs implicated in cognitive disorders and for genes involved in synaptic plasticity and cognitive functions. Our results suggest evolutionary changes in TF networks that might have shaped morphological and functional differences between primate brains, in particular in the human PFC.

The chromatin basis of neurodevelopmental disorders: Rethinking dysfunction along the molecular and temporal axes

The complexity of the human brain emerges from a long and finely tuned developmental process orchestrated by the crosstalk between genome and environment. Vis à vis other species, the human brain displays unique functional and morphological features that result from this extensive developmental process that is, unsurprisingly, highly vulnerable to both genetically and environmentally induced alterations. One of the most striking outcomes of the recent surge of sequencing-based studies on neurodevelopmental disorders (NDDs) is the emergence of chromatin regulation as one of the two domains most affected by causative mutations or Copy Number Variations besides synaptic function, whose involvement had been largely predicted for obvious reasons. These observations place chromatin dysfunction at the top of the molecular pathways hierarchy that ushers in a sizeable proportion of NDDs and that manifest themselves through synaptic dysfunction and recurrent systemic clinical manifestation. Here we undertake a conceptual investigation of chromatin dysfunction in NDDs with the aim of systematizing the available evidence in a new framework: first, we tease out the developmental vulnerabilities in human corticogenesis as a structuring entry point into the causation of NDDs; second, we provide a much needed clarification of the multiple meanings and explanatory frameworks revolving around "epigenetics", highlighting those that are most relevant for the analysis of these disorders; finally we go in-depth into paradigmatic examples of NDD-causing chromatin dysregulation, with a special focus on human experimental models and datasets.

Silencing of transposable elements may not be a major driver of regulatory evolution in primate iPSCs

Transposable elements (TEs) comprise almost half of primate genomes and their aberrant regulation can result in deleterious effects. In pluripotent stem cells, rapidly evolving KRAB-ZNF genes target TEs for silencing by H3K9me3. To investigate the evolution of TE silencing, we performed H3K9me3 ChIP-seq experiments in induced pluripotent stem cells from 10 human and 7 chimpanzee individuals. We identified four million orthologous TEs and found the SVA and ERV families to be marked most frequently by H3K9me3. We found little evidence of inter-species differences in TE silencing, with as many as 82% of putatively silenced TEs marked at similar levels in humans and chimpanzees. TEs that are preferentially silenced in one species are a similar age to those silenced in both species and are not more likely to be associated with expression divergence of nearby orthologous genes. Our data suggest limited species-specificity of TE silencing across 6 million years of primate evolution.

Repeat associated mechanisms of genome evolution and function revealed by the Mus caroli and Mus pahari genomes

Understanding the mechanisms driving lineage-specific evolution in both primates and rodents has been hindered by the lack of sister clades with a similar phylogenetic structure having high-quality genome assemblies. Here, we have created chromosome-level assemblies of the Mus caroli and Mus pahari genomes. Together with the Mus musculus and Rattus norvegicus genomes, this set of rodent genomes is similar in divergence times to the Hominidae (human-chimpanzee-gorilla-orangutan). By comparing the evolutionary dynamics between the Muridae and Hominidae, we identified punctate events of chromosome reshuffling that shaped the ancestral karyotype of Mus musculus and Mus caroli between 3 and 6 million yr ago, but that are absent in the Hominidae. Hominidae show between four- and sevenfold lower rates of nucleotide change and feature turnover in both neutral and functional sequences, suggesting an underlying coherence to the Muridae acceleration. Our system of matched, high-quality genome assemblies revealed how specific classes of repeats can play lineage-specific roles in related species. Recent LINE activity has remodeled protein-coding loci to a greater extent across the Muridae than the Hominidae, with functional consequences at the species level such as reproductive isolation. Furthermore, we charted a Muridae-specific retrotransposon expansion at unprecedented resolution, revealing how a single nucleotide mutation transformed a specific SINE element into an active CTCF binding site carrier specifically in Mus caroli, which resulted in thousands of novel, species-specific CTCF binding sites. Our results show that the comparison of matched phylogenetic sets of genomes will be an increasingly powerful strategy for understanding mammalian biology.

RTFAdb: A database of computationally predicted associations between retrotransposons and transcription factors in the human and mouse genomes

In recent years, retrotransposons have gained increasing attention as a source of binding motifs for transcription factors (TFs). Despite the substantial roles of these mobile genetic elements in the regulation of gene expression, a comprehensive resource enabling the investigation of retrotransposon species that are bound by TFs is still lacking. Herein, I introduce for the first time a novel database called RTFAdb, which allows exploring computationally predicted associations between retrotransposons and TFs in diverse cell lines and tissues of human and mouse. My database, using over 3.000 TF ChIP-seq binding profiles collected from human and mouse samples, makes possible searching more than 1.500 retrotransposon species in the binding sites of a total of 596 TFs. RTFAdb is freely available at http://tools.ibg.deu.edu.tr/rtfa/ and has the potential to offer novel insights into mammalian transcriptional networks by providing an additional layer of information regarding the regulatory roles of retrotransposons.

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

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

YY1 and CTCF orchestrate a 3D chromatin looping switch during early neural lineage commitment

CTCF is an architectural protein with a critical role in connecting higher-order chromatin folding in pluripotent stem cells. Recent reports have suggested that CTCF binding is more dynamic during development than previously appreciated. Here, we set out to understand the extent to which shifts in genome-wide CTCF occupancy contribute to the 3D reconfiguration of fine-scale chromatin folding during early neural lineage commitment. Unexpectedly, we observe a sharp decrease in CTCF occupancy during the transition from naïve/primed pluripotency to multipotent primary neural progenitor cells (NPCs). Many pluripotency gene-enhancer interactions are anchored by CTCF, and its occupancy is lost in parallel with loop decommissioning during differentiation. Conversely, CTCF binding sites in NPCs are largely preexisting in pluripotent stem cells. Only a small number of CTCF sites arise de novo in NPCs. We identify another zinc finger protein, Yin Yang 1 (YY1), at the base of looping interactions between NPC-specific genes and enhancers. Putative NPC-specific enhancers exhibit strong YY1 signal when engaged in 3D contacts and negligible YY1 signal when not in loops. Moreover, siRNA knockdown of Yy1 specifically disrupts interactions between key NPC enhancers and their target genes. YY1-mediated interactions between NPC regulatory elements are often nested within constitutive loops anchored by CTCF. Together, our results support a model in which YY1 acts as an architectural protein to connect developmentally regulated looping interactions; the location of YY1-mediated interactions may be demarcated in development by a preexisting topological framework created by constitutive CTCF-mediated interactions.

Conserved expression of transposon-derived non-coding transcripts in primate stem cells

BACKGROUND

A significant portion of expressed non-coding RNAs in human cells is derived from transposable elements (TEs). Moreover, it has been shown that various long non-coding RNAs (lncRNAs), which come from the human endogenous retrovirus subfamily H (HERVH), are not only expressed but required for pluripotency in human embryonic stem cells (hESCs).

RESULTS

To identify additional TE-derived functional non-coding transcripts, we generated RNA-seq data from induced pluripotent stem cells (iPSCs) of four primate species (human, chimpanzee, gorilla, and rhesus) and searched for transcripts whose expression was conserved. We observed that about 30% of TE instances expressed in human iPSCs had orthologous TE instances that were also expressed in chimpanzee and gorilla. Notably, our analysis revealed a number of repeat families with highly conserved expression profiles including HERVH but also MER53, which is known to be the source of a placental-specific family of microRNAs (miRNAs). We also identified a number of repeat families from all classes of TEs, including MLT1-type and Tigger families, that contributed a significant amount of sequence to primate lncRNAs whose expression was conserved.

CONCLUSIONS

Together, these results describe TE families and TE-derived lncRNAs whose conserved expression patterns can be used to identify what are likely functional TE-derived non-coding transcripts in primate iPSCs.

The cancer-associated CTCFL/BORIS protein targets multiple classes of genomic repeats, with a distinct binding and functional preference for humanoid-specific SVA transposable elements

Background

A common aberration in cancer is the activation of germline-specific proteins. The DNA-binding proteins among them could generate novel chromatin states, not found in normal cells. The germline-specific transcription factor BORIS/CTCFL, a paralog of chromatin architecture protein CTCF, is often erroneously activated in cancers and rewires the epigenome for the germline-like transcription program. Another common feature of malignancies is the changed expression and epigenetic states of genomic repeats, which could alter the transcription of neighboring genes and cause somatic mutations upon transposition. The role of BORIS in transposable elements and other repeats has never been assessed.

Results

The investigation of BORIS and CTCF binding to DNA repeats in the K562 cancer cells dependent on BORIS for self-renewal by ChIP-chip and ChIP-seq revealed three classes of occupancy by these proteins: elements cohabited by BORIS and CTCF, CTCF-only bound, or BORIS-only bound. The CTCF-only enrichment is characteristic for evolutionary old and inactive repeat classes, while BORIS and CTCF co-binding predominately occurs at uncharacterized tandem repeats. These repeats form staggered cluster binding sites, which are a prerequisite for CTCF and BORIS co-binding. At the same time, BORIS preferentially occupies a specific subset of the evolutionary young, transcribed, and mobile genomic repeat family, SVA. Unlike CTCF, BORIS prominently binds to the VNTR region of the SVA repeats in vivo. This suggests a role of BORIS in SVA expression regulation. RNA-seq analysis indicates that BORIS largely serves as a repressor of SVA expression, alongside DNA and histone methylation, with the exception of promoter capture by SVA.

Conclusions

Thus, BORIS directly binds to, and regulates SVA repeats, which are essentially movable CpG islands, via clusters of BORIS binding sites. This finding uncovers a new function of the global germline-specific transcriptional regulator BORIS in regulating and repressing the newest class of transposable elements that are actively transposed in human genome when activated. This function of BORIS in cancer cells is likely a reflection of its roles in the germline.

Electronic supplementary material

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

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Mitochondrial heteroplasmy in vertebrates using ChIP-sequencing data

BACKGROUND

Mitochondrial heteroplasmy, the presence of more than one mitochondrial DNA (mtDNA) variant in a cell or individual, is not as uncommon as previously thought. It is mostly due to the high mutation rate of the mtDNA and limited repair mechanisms present in the mitochondrion. Motivated by mitochondrial diseases, much focus has been placed into studying this phenomenon in human samples and in medical contexts. To place these results in an evolutionary context and to explore general principles of heteroplasmy, we describe an integrated cross-species evaluation of heteroplasmy in mammals that exploits previously reported NGS data. Focusing on ChIP-seq experiments, we developed a novel approach to detect heteroplasmy from the concomitant mitochondrial DNA fraction sequenced in these experiments.

RESULTS

We first demonstrate that the sequencing coverage of mtDNA in ChIP-seq experiments is sufficient for heteroplasmy detection. We then describe a novel detection method for accurate detection of heteroplasmies, which also accounts for the error rate of NGS technology. Applying this method to 79 individuals from 16 species resulted in 107 heteroplasmic positions present in a total of 45 individuals. Further analysis revealed that the majority of detected heteroplasmies occur in intergenic regions.

CONCLUSION

In addition to documenting the prevalence of mtDNA in ChIP-seq data, the results of our mitochondrial heteroplasmy detection method suggest that mitochondrial heteroplasmies identified across vertebrates share similar characteristics as found for human heteroplasmies. Although largely consistent with previous studies in individual vertebrates, our integrated cross-species analysis provides valuable insights into the evolutionary dynamics of mitochondrial heteroplasmy.

3D Chromosome Regulatory Landscape of Human Pluripotent Cells

In this study, we describe the 3D chromosome regulatory landscape of human naive and primed embryonic stem cells. To devise this map, we identified transcriptional enhancers and insulators in these cells and placed them within the context of cohesin-associated CTCF-CTCF loops using cohesin ChIA-PET data. The CTCF-CTCF loops we identified form a chromosomal framework of insulated neighborhoods, which in turn form topologically associating domains (TADs) that are largely preserved during the transition between the naive and primed states. Regulatory changes in enhancer-promoter interactions occur within insulated neighborhoods during cell state transition. The CTCF anchor regions we identified are conserved across species, influence gene expression, and are a frequent site of mutations in cancer cells, underscoring their functional importance in cellular regulation. These 3D regulatory maps of human pluripotent cells therefore provide a foundation for future interrogation of the relationships between chromosome structure and gene control in development and disease.

Towards a predictive model of chromatin 3D organization

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

Modeling co-occupancy of transcription factors using chromatin features

Regulation of gene expression requires both transcription factor (TFs) and epigenetic modifications, and interplays between the two types of factors have been discovered. However study of relationships between chromatin features and TF–TF co-occupancy remains limited. Here, we revealed the relationship by first illustrating distinct profile patterns of chromatin features related to different binding events, including single TF binding and TF–TF co-occupancy of 71 TFs from five human cell lines. We further implemented statistical analyses to demonstrate the relationship by accurately predicting co-occupancy genome-widely using chromatin features including DNase I hypersensitivity, 11 histone modifications (HMs) and GC content. Remarkably, our results showed that the combination of chromatin features enables accurate predictions across the five cells. For individual chromatin features, DNase I enables high and consistent predictions. H3K27ac, H3K4me 2, H3K4me3 and H3K9ac are more reliable predictors than other HMs. Although the combination of 11 HMs achieves accurate predictions, their predictive ability varies considerably when a model obtained from one cell is applied to others, indicating relationship between HMs and TF–TF co-occupancy is cell type dependent. GC content is not a reliable predictor, but the addition of GC content to any other features enhances their predictive ability. Together, our results elucidate a strong relationship between TF–TF co-occupancy and chromatin features.

Integrative modeling reveals the principles of multi-scale chromatin boundary formation in human nuclear organization

Background

Interphase chromosomes adopt a hierarchical structure, and recent data have characterized their chromatin organization at very different scales, from sub-genic regions associated with DNA-binding proteins at the order of tens or hundreds of bases, through larger regions with active or repressed chromatin states, up to multi-megabase-scale domains associated with nuclear positioning, replication timing and other qualities. However, we have lacked detailed, quantitative models to understand the interactions between these different strata.

Results

Here we collate large collections of matched locus-level chromatin features and Hi-C interaction data, representing higher-order organization, across three human cell types. We use quantitative modeling approaches to assess whether locus-level features are sufficient to explain higher-order structure, and identify the most influential underlying features. We identify structurally variable domains between cell types and examine the underlying features to discover a general association with cell-type-specific enhancer activity. We also identify the most prominent features marking the boundaries of two types of higher-order domains at different scales: topologically associating domains and nuclear compartments. We find parallel enrichments of particular chromatin features for both types, including features associated with active promoters and the architectural proteins CTCF and YY1.

Conclusions

We show that integrative modeling of large chromatin dataset collections using random forests can generate useful insights into chromosome structure. The models produced recapitulate known biological features of the cell types involved, allow exploration of the antecedents of higher-order structures and generate testable hypotheses for further experimental studies.

Electronic supplementary material

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

Explicit DNase sequence bias modeling enables high-resolution transcription factor footprint detection

DNaseI footprinting is an established assay for identifying transcription factor (TF)-DNA interactions with single base pair resolution. High-throughput DNase-seq assays have recently been used to detect in vivo DNase footprints across the genome. Multiple computational approaches have been developed to identify DNase-seq footprints as predictors of TF binding. However, recent studies have pointed to a substantial cleavage bias of DNase and its negative impact on predictive performance of footprinting. To assess the potential for using DNase-seq to identify individual binding sites, we performed DNase-seq on deproteinized genomic DNA and determined sequence cleavage bias. This allowed us to build bias corrected and TF-specific footprint models. The predictive performance of these models demonstrated that predicted footprints corresponded to high-confidence TF-DNA interactions. DNase-seq footprints were absent under a fraction of ChIP-seq peaks, which we show to be indicative of weaker binding, indirect TF-DNA interactions or possible ChIP artifacts. The modeling approach was also able to detect variation in the consensus motifs that TFs bind to. Finally, cell type specific footprints were detected within DNase hypersensitive sites that are present in multiple cell types, further supporting that footprints can identify changes in TF binding that are not detectable using other strategies.

CHD8 regulates neurodevelopmental pathways associated with autism spectrum disorder in neural progenitors

Truncating mutations of chromodomain helicase DNA-binding protein 8 (CHD8), and of many other genes with diverse functions, are strong-effect risk factors for autism spectrum disorder (ASD), suggesting multiple mechanisms of pathogenesis. We explored the transcriptional networks that CHD8 regulates in neural progenitor cells (NPCs) by reducing its expression and then integrating transcriptome sequencing (RNA sequencing) with genome-wide CHD8 binding (ChIP sequencing). Suppressing CHD8 to levels comparable with the loss of a single allele caused altered expression of 1,756 genes, 64.9% of which were up-regulated. CHD8 showed widespread binding to chromatin, with 7,324 replicated sites that marked 5,658 genes. Integration of these data suggests that a limited array of direct regulatory effects of CHD8 produced a much larger network of secondary expression changes. Genes indirectly down-regulated (i.e., without CHD8-binding sites) reflect pathways involved in brain development, including synapse formation, neuron differentiation, cell adhesion, and axon guidance, whereas CHD8-bound genes are strongly associated with chromatin modification and transcriptional regulation. Genes associated with ASD were strongly enriched among indirectly down-regulated loci (P < 10(-8)) and CHD8-bound genes (P = 0.0043), which align with previously identified coexpression modules during fetal development. We also find an intriguing enrichment of cancer-related gene sets among CHD8-bound genes (P < 10(-10)). In vivo suppression of chd8 in zebrafish produced macrocephaly comparable to that of humans with inactivating mutations. These data indicate that heterozygous disruption of CHD8 precipitates a network of gene-expression changes involved in neurodevelopmental pathways in which many ASD-associated genes may converge on shared mechanisms of pathogenesis.

Gibbon genome and the fast karyotype evolution of small apes

Gibbons are small arboreal apes that display an accelerated rate of evolutionary chromosomal rearrangement and occupy a key node in the primate phylogeny between Old World monkeys and great apes. Here we present the assembly and analysis of a northern white-cheeked gibbon (Nomascus leucogenys) genome. We describe the propensity for a gibbon-specific retrotransposon (LAVA) to insert into chromosome segregation genes and alter transcription by providing a premature termination site, suggesting a possible molecular mechanism for the genome plasticity of the gibbon lineage. We further show that the gibbon genera (Nomascus, Hylobates, Hoolock and Symphalangus) experienced a near-instantaneous radiation ∼5 million years ago, coincident with major geographical changes in southeast Asia that caused cycles of habitat compression and expansion. Finally, we identify signatures of positive selection in genes important for forelimb development (TBX5) and connective tissues (COL1A1) that may have been involved in the adaptation of gibbons to their arboreal habitat.

The genome of the gibbon, a tree-dwelling ape from Asia positioned between Old World monkeys and the great apes, is presented, providing insights into the evolutionary history of gibbon species and their accelerated karyotypes, as well as evidence for selection of genes such as those for forelimb development and connective tissue that may be important for locomotion through trees.

The many species of gibbons are small, tree-living apes from Southeast Asia, most of them listed as 'endangered' or 'critically endangered' on the IUCN list. In their presentation of the genome of the northern white-cheeked gibbon (Nomascus leucogenys) , Lucia Carbone and colleagues provide intriguing insights into the biology and evolutionary history of a group that straddles the divide between Old World monkeys and the great apes. The authors investigate how a novel gibbon-specific retrotransposon might be the source of gibbons' genome plasticity. Rapid karyotype evolution combined with multiple episodes of climate and environmental change might explain the almost instantaneous divergence of the four gibbon genera. Positive selection on genes involved in forelimb development and connective tissue might have been related to gibbons' unique mode of locomotion in the tropical canopy.