Segmental folding of chromosomes: a basis for structural and regulatory chromosomal neighborhoods?

Tracing the Shared Foundations of Gene Expression and Chromatin Structure

The three-dimensional organization of chromatin into topologically associating domains (TADs) may impact gene regulation by bringing distant genes into contact. However, many questions about TADs' function and their influence on transcription remain unresolved due to technical limitations in defining TAD boundaries and measuring the direct effect that TADs have on gene expression. Here, we develop consensus TAD maps for human and mouse with a novel "bag-of-genes" approach for defining the gene composition within TADs. This approach enables new functional interpretations of TADs by providing a way to capture species-level differences in chromatin organization. We also leverage a generative AI foundation model computed from 33 million transcriptomes to define contextual similarity, an embedding-based metric that is more powerful than co-expression at representing functional gene relationships. Our analytical framework directly leads to testable hypotheses about chromatin organization across cellular states. We find that TADs play an active role in facilitating gene co-regulation, possibly through a mechanism involving transcriptional condensates. We also discover that the TAD-linked enhancement of transcriptional context is strongest in early developmental stages and systematically declines with aging. Investigation of cancer cells show distinct patterns of TAD usage that shift with chemotherapy treatment, suggesting specific roles for TAD-mediated regulation in cellular development and plasticity. Finally, we develop "TAD signatures" to improve statistical analysis of single-cell transcriptomic data sets in predicting cancer cell-line drug response. These findings reshape our understanding of cellular plasticity in development and disease, indicating that chromatin organization acts through probabilistic mechanisms rather than deterministic rules.

Software availability

https://singhlab.net/tadmap.

Insulation between adjacent TADs is controlled by the width of their boundaries through distinct mechanisms

Topologically associating domains (TADs) are sub-Megabase regions in vertebrate genomes with enriched intradomain interactions that restrict enhancer-promoter contacts across their boundaries. However, the mechanisms that separate TADs remain incompletely understood. Most boundaries between TADs contain CTCF binding sites (CBSs), which individually contribute to the blocking of Cohesin-mediated loop extrusion. Using genome-wide classification, here we show that the width of TAD boundaries forms a continuum from narrow to highly extended and correlates with CBSs distribution, chromatin features, and gene regulatory elements. To investigate how these boundary widths emerge, we modified the random crosslinker polymer model to incorporate specific boundary configurations, enabling us to evaluate the differential impact of boundary composition on TAD insulation. Our analysis, using three generic boundary categories, identifies differential influence on TAD insulation, with varying local and distal effects on neighboring domains. Notably, we find that increasing boundary width reduces long-range inter-TAD contacts, as confirmed by Hi-C data. While blocking loop extrusion at boundaries indirectly promotes spurious intermingling of neighboring TADs, extended boundaries counteract this effect, emphasizing their role in establishing genome organization. In conclusion, TAD boundary width not only enhances the efficiency of loop extrusion blocking but may also modulate enhancer-promoter contacts over long distances across TAD boundaries, providing a further mechanism for transcriptional regulation.

The impact of the embryonic DNA methylation program on CTCF-mediated genome regulation

During mammalian embryogenesis, both the 5-cytosine DNA methylation (5meC) landscape and three dimensional (3D) chromatin architecture are profoundly remodeled during a process known as 'epigenetic reprogramming.' An understudied aspect of epigenetic reprogramming is how the 5meC flux, per se, affects the 3D genome. This is pertinent given the 5meC-sensitivity of DNA binding for a key regulator of chromosome folding: CTCF. We profiled the CTCF binding landscape using a mouse embryonic stem cell (ESC) differentiation protocol that models embryonic 5meC dynamics. Mouse ESCs lacking DNA methylation machinery are able to exit naive pluripotency, thus allowing for dissection of subtle effects of CTCF on gene expression. We performed CTCF HiChIP in both wild-type and mutant conditions to assess gained CTCF-CTCF contacts in the absence of 5meC. We performed H3K27ac HiChIP to determine the impact that ectopic CTCF binding has on cis-regulatory contacts. Using 5meC epigenome editing, we demonstrated that the methyl-mark is able to impair CTCF binding at select loci. Finally, a detailed dissection of the imprinted Zdbf2 locus showed how 5meC-antagonism of CTCF allows for proper gene regulation during differentiation. This work provides a comprehensive overview of how 5meC impacts the 3D genome in a relevant model for early embryonic events.

Permeable TAD boundaries and their impact on genome-associated functions

TAD boundaries are genomic elements that separate biological processes in neighboring domains by blocking DNA loops that are formed through Cohesin-mediated loop extrusion. Most TAD boundaries consist of arrays of binding sites for the CTCF protein, whose interaction with the Cohesin complex blocks loop extrusion. TAD boundaries are not fully impermeable though and allow a limited amount of inter-TAD loop formation. Based on the reanalysis of Nano-C data, a multicontact Chromosome Conformation Capture assay, we propose a model whereby clustered CTCF binding sites promote the successive stalling of Cohesin and subsequent dissociation from the chromatin. A fraction of Cohesin nonetheless achieves boundary read-through. Due to a constant rate of Cohesin dissociation elsewhere in the genome, the maximum length of inter-TAD loops is restricted though. We speculate that the DNA-encoded organization of stalling sites regulates TAD boundary permeability and discuss implications for enhancer-promoter loop formation and other genomic processes.

A genomic hotspot of diversifying selection and structural change in the hoary bat (Lasiurus cinereus)

Background

Previous work found that numerous genes positively selected within the hoary bat (Lasiurus cinereus) lineage are physically clustered in regions of conserved synteny. Here I further validate and expand on those finding utilizing an updated L. cinereus genome assembly and additional bat species as well as other tetrapod outgroups.

Methods

A chromosome-level assembly was generated by chromatin-contact mapping and made available by DNAZoo (www.dnazoo.org). The genomic organization of orthologous genes was extracted from annotation data for multiple additional bat species as well as other tetrapod clades for which chromosome-level assemblies were available from the National Center for Biotechnology Information (NCBI). Tests of branch-specific positive selection were performed for L. cinereus using PAML as well as with the HyPhy package for comparison.

Results

Twelve genes exhibiting significant diversifying selection in the L. cinereus lineage were clustered within a 12-Mb genomic window; one of these (Trpc4) also exhibited diversifying selection in bats generally. Ten of the 12 genes are landmarks of two distinct blocks of ancient synteny that are not linked in other tetrapod clades. Bats are further distinguished by frequent structural rearrangements within these synteny blocks, which are rarely observed in other Tetrapoda. Patterns of gene order and orientation among bat taxa are incompatible with phylogeny as presently understood, implying parallel evolution or subsequent reversals. Inferences of positive selection were found to be robust to alternative phylogenetic topologies as well as a strong shift in background nucleotide composition in some taxa.

Discussion

This study confirms and further localizes a genomic hotspot of protein-coding divergence in the hoary bat, one that also exhibits an increased tempo of structural change in bats compared with other mammals. Most genes in the two synteny blocks have elevated expression in brain tissue in humans and model organisms, and genetic studies implicate the selected genes in cranial and neurological development, among other functions.

Hi-C sequencing unravels dynamic three-dimensional chromatin interactions in muntjac lineage: insights from chromosome fusions in Fea's muntjac genome

Eukaryotes have varying numbers and structures of characteristic chromosomes across lineages or species. The evolutionary trajectory of species may have been affected by spontaneous genome rearrangements. Chromosome fusion drastically alters karyotypes. However, the mechanisms and consequences of chromosome fusions, particularly in muntjac species, are poorly understood. Recent research-based advancements in three-dimensional (3D) genomics, particularly high-throughput chromatin conformation capture (Hi-C) sequencing, have allowed for the identification of chromosome fusions and provided mechanistic insights into three muntjac species: Muntiacus muntjak, M. reevesi, and M. crinifrons. This study aimed to uncover potential genome rearrangement patterns in the threatened species Fea's muntjac (Muntiacus feae), which have not been previously examined for such characteristics. Deep Hi-C sequencing (31.42 × coverage) was performed to reveal the 3D chromatin architecture of the Fea's muntjac genome. Patterns of repeated chromosome fusions that were potentially mediated by high-abundance transposable elements were identified. Comparative Hi-C maps demonstrated linkage homology between the sex chromosomes in Fea's muntjac and autosomes in M. reevesi, indicating that fusions may have played a crucial role in the evolution of the sex chromosomes of the lineage. The species-level dynamics of topologically associated domains (TADs) suggest that TAD organization could be altered by differential chromosome interactions owing to repeated chromosome fusions. However, research on the effect of TADs on muntjac genome evolution is insufficient. This study generated Hi-C data for the Fea's muntjac, providing a genomic resource for future investigations of the evolutionary patterns of chromatin conformation at the chromosomal level.

Meiosis: Dances Between Homologs

The raison d'être of meiosis is shuffling of genetic information via Mendelian segregation and, within individual chromosomes, by DNA crossing-over. These outcomes are enabled by a complex cellular program in which interactions between homologous chromosomes play a central role. We first provide a background regarding the basic principles of this program. We then summarize the current understanding of the DNA events of recombination and of three processes that involve whole chromosomes: homolog pairing, crossover interference, and chiasma maturation. All of these processes are implemented by direct physical interaction of recombination complexes with underlying chromosome structures. Finally, we present convergent lines of evidence that the meiotic program may have evolved by coupling of this interaction to late-stage mitotic chromosome morphogenesis.

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.

A model for cis-regulation of transcriptional condensates and gene expression by proximal lncRNAs

Long noncoding RNAs (lncRNAs) perform several important functions in cells including cis-regulation of transcription. Barring a few specific cases, the mechanisms underlying transcriptional regulation by lncRNAs remain poorly understood. Transcriptional proteins can form condensates via phase separation at protein-binding loci (BL) on the genome (e.g., enhancers and promoters). lncRNA-coding genes are present at loci in close genomic proximity of these BL and these RNAs can interact with transcriptional proteins via attractive heterotypic interactions mediated by their net charge. Motivated by these observations, we propose that lncRNAs can dynamically regulate transcription in cis via charge-based heterotypic interactions with transcriptional proteins in condensates. To study the consequences of this mechanism, we developed and studied a dynamical phase-field model. We find that proximal lncRNAs can promote condensate formation at the BL. Vicinally localized lncRNA can migrate to the BL to attract more protein because of favorable interaction free energies. However, increasing the distance beyond a threshold leads to a sharp decrease in protein recruitment to the BL. This finding could potentially explain why genomic distances between lncRNA-coding genes and protein-coding genes are conserved across metazoans. Finally, our model predicts that lncRNA transcription can fine-tune transcription from neighboring condensate-controlled genes, repressing transcription from highly expressed genes and enhancing transcription of genes expressed at a low level. This nonequilibrium effect can reconcile conflicting reports that lncRNAs can enhance or repress transcription from proximal genes.

POSTRE: a tool to predict the pathological effects of human structural variants

Understanding the pathological impact of non-coding genetic variation is a major challenge in medical genetics. Accumulating evidences indicate that a significant fraction of genetic alterations, including structural variants (SVs), can cause human disease by altering the function of non-coding regulatory elements, such as enhancers. In the case of SVs, described pathomechanisms include changes in enhancer dosage and long-range enhancer-gene communication. However, there is still a clear gap between the need to predict and interpret the medical impact of non-coding variants, and the existence of tools to properly perform these tasks. To reduce this gap, we have developed POSTRE (Prediction Of STRuctural variant Effects), a computational tool to predict the pathogenicity of SVs implicated in a broad range of human congenital disorders. By considering disease-relevant cellular contexts, POSTRE identifies SVs with either coding or long-range pathological consequences with high specificity and sensitivity. Furthermore, POSTRE not only identifies pathogenic SVs, but also predicts the disease-causative genes and the underlying pathological mechanism (e.g, gene deletion, enhancer disconnection, enhancer adoption, etc.). POSTRE is available at https://github.com/vicsanga/Postre.

Deciphering the chromatin spatial organization landscapes during BMMSC differentiation

The differentiation imbalance in bone marrow mesenchymal stem cells (BMMSCs) is critical for the development of bone density diseases as the population ages. BMMSCs are precursor cells for osteoblasts and adipocytes; however, the chromatin organization landscapes during BMMSC differentiation remain elusive. In this study, we systematically delineate the four-dimensional (4D) genome and dynamic epigenetic atlas of BMMSCs by RNA sequencing (RNA-seq), assay for transposase-accessible chromatin sequencing (ATAC-seq), and high-throughput chromosome conformation capture (Hi-C). The structure analyses reveal 17.5% common and 28.5%-30% specific loops among BMMSCs, osteoblasts, and adipocytes. The subsequent correlation of genome-wide association studies (GWAS) and expression quantitative trait locus (eQTL) data with multi-omics analysis reveal 274 genes and 3634 single nucleotide polymorphisms (SNPs) associated with bone degeneration and osteoporosis (OP). We hypothesize that SNP mutations affect transcription factor (TF) binding sites, thereby affecting changes in gene expression. Furthermore, 26 motifs, 260 TFs, and 291 SNPs are identified to affect the eQTL. Among these genes, DAAM2, TIMP2, and TMEM241 were found to be essential for diseases such as bone degeneration and OP and may serve as potential drug targets.

CLNN-loop: a deep learning model to predict CTCF-mediated chromatin loops in the different cell lines and CTCF-binding sites (CBS) pair types

MOTIVATION

Three-dimensional (3D) genome organization is of vital importance in gene regulation and disease mechanisms. Previous studies have shown that CTCF-mediated chromatin loops are crucial to studying the 3D structure of cells. Although various experimental techniques have been developed to detect chromatin loops, they have been found to be time-consuming and costly. Nowadays, various sequence-based computational methods can capture significant features of 3D genome organization and help predict chromatin loops. However, these methods have low performance and poor generalization ability in predicting chromatin loops.

RESULTS

Here, we propose a novel deep learning model, called CLNN-loop, to predict chromatin loops in different cell lines and CTCF-binding sites (CBS) pair types by fusing multiple sequence-based features. The analysis of a series of examinations based on the datasets in the previous study shows that CLNN-loop has satisfactory performance and is superior to the existing methods in terms of predicting chromatin loops. In addition, we apply the SHAP framework to interpret the predictions of different models, and find that CTCF motif and sequence conservation are important signs of chromatin loops in different cell lines and CBS pair types.

AVAILABILITY AND IMPLEMENTATION

The source code of CLNN-loop is freely available at https://github.com/HaoWuLab-Bioinformatics/CLNN-loop and the webserver of CLNN-loop is freely available at http://hwclnn.sdu.edu.cn.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Complex chromosomal rearrangements induced by transposons in maize

Eukaryotic genomes are large and complex, and gene expression can be affected by multiple regulatory elements and their positions within the dynamic chromatin architecture. Transposable elements are known to play important roles in genome evolution, yet questions remain as to how transposable elements alter genome structure and affect gene expression. Previous studies have shown that genome rearrangements can be induced by Reversed Ends Transposition involving termini of Activator and related transposable elements in maize and other plants. Here, we show that complex alleles can be formed by the rapid and progressive accumulation of Activator-induced duplications and rearrangements. The p1 gene enhancer in maize can induce ectopic expression of the nearby p2 gene in pericarp tissue when placed near it via different structural rearrangements. By screening for p2 expression, we identified and studied 5 cases in which multiple sequential transposition events occurred and increased the p1 enhancer copy number. We see active p2 expression due to multiple copies of the p1 enhancer present near p2 in all 5 cases. The p1 enhancer effects are confirmed by the observation that loss of p2 expression is correlated with transposition-induced excision of the p1 enhancers. We also performed a targeted Chromosome Conformation Capture experiment to test the physical interaction between the p1 enhancer and p2 promoter region. Together, our results show that transposon-induced rearrangements can accumulate rapidly and progressively increase genetic variation important for genomic evolution.

Reconstituted TAD-size chromatin fibers feature heterogeneous nucleosome clusters

Large topologically associated domains (TADs) contain irregularly spaced nucleosome clutches, and interactions between such clutches are thought to aid the compaction of these domains. Here, we reconstituted TAD-sized chromatin fibers containing hundreds of nucleosomes on native source human and lambda-phage DNA and compared their mechanical properties at the single-molecule level with shorter ‘601’ arrays with various nucleosome repeat lengths. Fluorescent imaging showed increased compaction upon saturation of the DNA with histones and increasing magnesium concentration. Nucleosome clusters and their structural fluctuations were visualized in confined nanochannels. Force spectroscopy revealed not only similar mechanical properties of the TAD-sized fibers as shorter fibers but also large rupture events, consistent with breaking the interactions between distant clutches of nucleosomes. Though the arrays of native human DNA, lambda-phage and ‘601’ DNA featured minor differences in reconstitution yield and nucleosome stability, the fibers’ global structural and mechanical properties were similar, including the interactions between nucleosome clutches. These single-molecule experiments quantify the mechanical forces that stabilize large TAD-sized chromatin domains consisting of disordered, dynamically interacting nucleosome clutches and their effect on the condensation of large chromatin domains.

The Molecular and Nuclear Dynamics of X-Chromosome Inactivation

In female eutherian mammals, dosage compensation of X-linked gene expression is achieved during development through transcriptional silencing of one of the two X chromosomes. Following X chromosome inactivation (XCI), the inactive X chromosome remains faithfully silenced throughout somatic cell divisions. XCI is dependent on Xist, a long noncoding RNA that coats and silences the X chromosome from which it is transcribed. Xist coating triggers a cascade of chromosome-wide changes occurring at the levels of transcription, chromatin composition, chromosome structure, and spatial organization within the nucleus. XCI has emerged as a paradigm for the study of such crucial nuclear processes and the dissection of their functional interplay. In the past decade, the advent of tools to characterize and perturb these processes have provided an unprecedented understanding into their roles during XCI. The mechanisms orchestrating the initiation of XCI as well as its maintenance are thus being unraveled, although many questions still remain. Here, we introduce key aspects of the XCI process and review the recent discoveries about its molecular basis.

CTCF: a misguided jack-of-all-trades in cancer cells

The emergence and progression of cancers is accompanied by a dysregulation of transcriptional programs. The three-dimensional (3D) organization of the human genome has emerged as an important multi-level mediator of gene transcription and regulation. In cancer cells, this organization can be restructured, providing a framework for the deregulation of gene activity. The CTCF protein, initially identified as the product from a tumor suppressor gene, is a jack-of-all-trades for the formation of 3D genome organization in normal cells. Here, we summarize how CTCF is involved in the multi-level organization of the human genome and we discuss emerging insights into how perturbed CTCF function and DNA binding causes the activation of oncogenes in cancer cells, mostly through a process of enhancer hijacking. Moreover, we highlight non-canonical functions of CTCF that can be relevant for the emergence of cancers as well. Finally, we provide guidelines for the computational identification of perturbed CTCF binding and reorganized 3D genome structure in cancer cells.

ProbC: joint modeling of epigenome and transcriptome effects in 3D genome

BACKGROUND

Hi-C and its high nucleosome resolution variant Micro-C provide a window into the spatial packing of a genome in 3D within the cell. Even though both techniques do not directly depend on the binding of specific antibodies, previous work has revealed enriched interactions and domain structures around multiple chromatin marks; epigenetic modifications and transcription factor binding sites. However, the joint impact of chromatin marks in Hi-C and Micro-C interactions have not been globally characterized, which limits our understanding of 3D genome characteristics. An emerging question is whether it is possible to deduce 3D genome characteristics and interactions by integrative analysis of multiple chromatin marks and associate interactions to functionality of the interacting loci.

RESULT

We come up with a probabilistic method PROBC to decompose Hi-C and Micro-C interactions by known chromatin marks. PROBC is based on convex likelihood optimization, which can directly take into account both interaction existence and nonexistence. Through PROBC, we discover histone modifications (H3K27ac, H3K9me3, H3K4me3, H3K4me1) and CTCF as particularly predictive of Hi-C and Micro-C contacts across cell types and species. Moreover, histone modifications are more effective than transcription factor binding sites in explaining the genome's 3D shape through these interactions. PROBC can successfully predict Hi-C and Micro-C interactions in given species, while it is trained on different cell types or species. For instance, it can predict missing nucleosome resolution Micro-C interactions in human ES cells trained on mouse ES cells only from these 5 chromatin marks with above 0.75 AUC. Additionally, PROBC outperforms the existing methods in predicting interactions across almost all chromosomes.

CONCLUSION

Via our proposed method, we optimally decompose Hi-C interactions in terms of these chromatin marks at genome and chromosome levels. We find a subset of histone modifications and transcription factor binding sites to be predictive of both Hi-C and Micro-C interactions and TADs across human, mouse, and different cell types. Through learned models, we can predict interactions on species just from chromatin marks for which Hi-C data may be limited.

StackTADB: a stacking-based ensemble learning model for predicting the boundaries of topologically associating domains (TADs) accurately in fruit flies

Chromosome is composed of many distinct chromatin domains, referred to variably as topological domains or topologically associating domains (TADs). The domains are stable across different cell types and highly conserved across species, thus these chromatin domains have been considered as the basic units of chromosome folding and regarded as an important secondary structure in chromosome organization. However, the identification of TAD boundaries is still a great challenge due to the high cost and low resolution of Hi-C data or experiments. In this study, we propose a novel ensemble learning framework, termed as StackTADB, for predicting the boundaries of TADs. StackTADB integrates four base classifiers including Random Forest, Logistic Regression, K-NearestNeighbor and Support Vector Machine. From the analysis of a series of examinations on the data set in the previous study, it is concluded that StackTADB has optimal performance in six metrics, AUC, Accuracy, MCC, Precision, Recall and F1 score, and it is superior to the existing methods. In addition, the comparison of the performance of multiple features shows that Kmers-based features play an essential role in predicting TADs boundaries of fruit flies, and we also apply the SHapley Additive exPlanations (SHAP) framework to interpret the predictions of StackTADB to identify the reason why Kmers-based features are vital. The experimental results show that the subsequences matching the BEAF-32 motif play a crucial role in predicting the boundaries of TADs. The source code is freely available at https://github.com/HaoWuLab-Bioinformatics/StackTADB and the webserver of StackTADB is freely available at http://hwtad.sdu.edu.cn:8002/StackTADB.

Mechanisms of sex determination and X-chromosome dosage compensation

Abnormalities in chromosome number have the potential to disrupt the balance of gene expression and thereby decrease organismal fitness and viability. Such abnormalities occur in most solid tumors and also cause severe developmental defects and spontaneous abortions. In contrast to the imbalances in chromosome dose that cause pathologies, the difference in X-chromosome dose used to determine sexual fate across diverse species is well tolerated. Dosage compensation mechanisms have evolved in such species to balance X-chromosome gene expression between the sexes, allowing them to tolerate the difference in X-chromosome dose. This review analyzes the chromosome counting mechanism that tallies X-chromosome number to determine sex (XO male and XX hermaphrodite) in the nematode Caenorhabditis elegans and the associated dosage compensation mechanism that balances X-chromosome gene expression between the sexes. Dissecting the molecular mechanisms underlying X-chromosome counting has revealed how small quantitative differences in intracellular signals can be translated into dramatically different fates. Dissecting the process of X-chromosome dosage compensation has revealed the interplay between chromatin modification and chromosome structure in regulating gene expression over vast chromosomal territories.

Position effects at the FGF8 locus are associated with femoral hypoplasia

Copy-number variations (CNVs) are a common cause of congenital limb malformations and are interpreted primarily on the basis of their effect on gene dosage. However, recent studies show that CNVs also influence the 3D genome chromatin organization. The functional interpretation of whether a phenotype is the result of gene dosage or a regulatory position effect remains challenging. Here, we report on two unrelated families with individuals affected by bilateral hypoplasia of the femoral bones, both harboring de novo duplications on chromosome 10q24.32. The ∼0.5 Mb duplications include FGF8, a key regulator of limb development and several limb enhancer elements. To functionally characterize these variants, we analyzed the local chromatin architecture in the affected individuals' cells and re-engineered the duplications in mice by using CRISPR-Cas9 genome editing. We found that the duplications were associated with ectopic chromatin contacts and increased FGF8 expression. Transgenic mice carrying the heterozygous tandem duplication including Fgf8 exhibited proximal shortening of the limbs, resembling the human phenotype. To evaluate whether the phenotype was a result of gene dosage, we generated another transgenic mice line, carrying the duplication on one allele and a concurrent Fgf8 deletion on the other allele, as a control. Surprisingly, the same malformations were observed. Capture Hi-C experiments revealed ectopic interaction with the duplicated region and Fgf8, indicating a position effect. In summary, we show that duplications at the FGF8 locus are associated with femoral hypoplasia and that the phenotype is most likely the result of position effects altering FGF8 expression rather than gene dosage effects.

Seeing the forest through the trees: prioritising potentially functional interactions from Hi-C

Eukaryotic genomes are highly organised within the nucleus of a cell, allowing widely dispersed regulatory elements such as enhancers to interact with gene promoters through physical contacts in three-dimensional space. Recent chromosome conformation capture methodologies such as Hi-C have enabled the analysis of interacting regions of the genome providing a valuable insight into the three-dimensional organisation of the chromatin in the nucleus, including chromosome compartmentalisation and gene expression. Complicating the analysis of Hi-C data, however, is the massive amount of identified interactions, many of which do not directly drive gene function, thus hindering the identification of potentially biologically functional 3D interactions. In this review, we collate and examine the downstream analysis of Hi-C data with particular focus on methods that prioritise potentially functional interactions. We classify three groups of approaches: structural-based discovery methods, e.g. A/B compartments and topologically associated domains, detection of statistically significant chromatin interactions, and the use of epigenomic data integration to narrow down useful interaction information. Careful use of these three approaches is crucial to successfully identifying potentially functional interactions within the genome.

Mapping the micro-proteome of the nuclear lamina and lamina-associated domains

The nuclear lamina is a proteinaceous network of filaments that provide both structural and gene regulatory functions by tethering proteins and large domains of DNA, the so-called lamina-associated domains (LADs), to the periphery of the nucleus. LADs are a large fraction of the mammalian genome that are repressed, in part, by their association to the nuclear periphery. The genesis and maintenance of LADs is poorly understood as are the proteins that participate in these functions. In an effort to identify proteins that reside at the nuclear periphery and potentially interact with LADs, we have taken a two-pronged approach. First, we have undertaken an interactome analysis of the inner nuclear membrane bound LAP2β to further characterize the nuclear lamina proteome. To accomplish this, we have leveraged the BioID system, which previously has been successfully used to characterize the nuclear lamina proteome. Second, we have established a system to identify proteins that bind to LADs by developing a chromatin-directed BioID system. We combined the BioID system with the m6A-tracer system which binds to LADs in live cells to identify both LAD proximal and nuclear lamina proteins. In combining these datasets, we have further characterized the protein network at the nuclear lamina, identified putative LAD proximal proteins and found several proteins that appear to interface with both micro-proteomes. Importantly, several proteins essential for LAD function, including heterochromatin regulating proteins related to H3K9 methylation, were identified in this study.

Dysregulation of chromatin organization in pediatric and adult brain tumors: oncoepigenomic contributions to tumorigenesis and cancer stem cell properties

The three-dimensional (3D) organization of the genome is a crucial enabler of cell fate, identity, and function. In this review, we will focus on the emerging role of altered 3D genome organization in the etiology of disease, with a special emphasis on brain cancers. We discuss how different genetic alterations can converge to disrupt the epigenome in childhood and adult brain tumors, by causing aberrant DNA methylation and by affecting the amounts and genomic distribution of histone post-translational modifications. We also highlight examples that illustrate how epigenomic alterations have the potential to affect 3D genome architecture in brain tumors. Finally, we will propose the concept of "epigenomic erosion" to explain the transition from stem-like cells to differentiated cells in hierarchically organized brain cancers.

Rare or Overlooked? Structural Disruption of Regulatory Domains in Human Neurocristopathies

In the last few years, the role of non-coding regulatory elements and their involvement in human disease have received great attention. Among the non-coding regulatory sequences, enhancers are particularly important for the proper establishment of cell type-specific gene-expression programs. Furthermore, the disruption of enhancers can lead to human disease through two main mechanisms: (i) Mutations or copy number variants can directly alter the enhancer sequences and thereby affect expression of their target genes; (ii) structural variants can provoke changes in 3-D chromatin organization that alter neither the enhancers nor their target genes, but rather the physical communication between them. In this review, these pathomechanisms are mostly discussed in the context of neurocristopathies, congenital disorders caused by defects that occur during neural crest development. We highlight why, due to its contribution to multiple tissues and organs, the neural crest represents an important, yet understudied, cell type involved in multiple congenital disorders. Moreover, we discuss currently available resources and experimental models for the study of human neurocristopathies. Last, we provide some practical guidelines that can be followed when investigating human neurocristopathies caused by structural variants. Importantly, these guidelines can be useful not only to uncover the etiology of human neurocristopathies, but also of other human congenital disorders in which enhancer disruption is involved.

Computing 3D Chromatin Configurations from Contact Probability Maps by Inverse Brownian Dynamics

The three-dimensional (3D) organization of chromatin, on the length scale of a few genes, is crucial in determining the functional state-accessibility and amount of gene expression-of the chromatin. Recent advances in chromosome conformation capture experiments provide partial information on the chromatin organization in a cell population, namely the contact count between any segment pairs, but not on the interaction strength that leads to these contact counts. However, given the contact matrix, determining the complete 3D organization of the whole chromatin polymer is an inverse problem. In this work, a novel inverse Brownian dynamics method based on a coarse-grained bead-spring chain model has been proposed to compute the optimal interaction strengths between different segments of chromatin such that the experimentally measured contact count probability constraints are satisfied. Applying this method to the α-globin gene locus in two different cell types, we predict the 3D organizations corresponding to active and repressed states of chromatin at the locus. We show that the average distance between any two segments of the region has a broad distribution and cannot be computed as a simple inverse relation based on the contact probability alone. The results presented for multiple normalization methods suggest that all measurable quantities may crucially depend on the nature of normalization. We argue that by experimentally measuring predicted quantities, one may infer the appropriate form of normalization.

Programming PAM antennae for efficient CRISPR-Cas9 DNA editing

Bacterial CRISPR-Cas9 nucleases have been repurposed as powerful genome editing tools. Whereas engineering guide RNAs or Cas nucleases have proven to improve the efficiency of CRISPR editing, modulation of protospacer-adjacent motif (PAM), indispensable for CRISPR, has been less explored. Here, we develop a DNA origami-based platform to program a PAM antenna microenvironment and address its performance at the single-molecule level with submolecular resolution. To mimic spatially controlled in vivo PAM distribution as may occur in chromatin, we investigate the effect of PAM antennae surrounding target DNA. We find that PAM antennae effectively sensitize the DNA cleavage by recruiting Cas9 molecules. Super-resolution tracking of single single-guide RNA/Cas9s reveals localized translocation of Cas9 among spatially proximal PAMs. We find that the introduction of the PAM antennae effectively modulates the microenvironment for enhanced target cleavage (up to ~50%). These results provide insight into factors that promote more efficient genome editing.

White spot syndrome viral protein VP9 alters the cellular higher-order chromatin structure

Viral protein 9 (VP9) is a non-structural protein of white spot syndrome virus (WSSV) highly expressed during the early stage of infection. The crystal structure of VP9 suggests that the polymers of VP9 dimers resemble a DNA mimic, but its function remains elusive. In this study, we demonstrated that VP9 impedes histones binding to DNA via single-molecule manipulation. We established VP9 expression in HeLa cells due to the lack of a WSSV-susceptible cell line, and observed abundant VP9 in the nucleus, which mirrors its distribution in the hemocytes of WSSV-infected shrimp. VP9 expression increased the dynamics and rotational mobility of histones in stable H3-GFP HeLa cells as revealed by fluorescent recovery after photobleaching and fluorescence anisotropy imaging, which suggested a loosened compaction of chromatin structure. Successive salt fractionation showed that a prominent population of histones was solubilized in high salt concentrations, which implies alterations of bulk chromatin structure. Southern blotting identified that VP9 alters juxtacentromeric chromatin structures to be more accessible to micrococcal nuclease digestion. RNA microarray revealed that VP9 expression also leads to significant changes of cellular gene expression. Our findings provide evidence that VP9 alters the cellular higher-order chromatin structure, uncovering a potential strategy adopted by WSSV to facilitate its replication.

Forged by DXZ4, FIRRE, and ICCE: How Tandem Repeats Shape the Active and Inactive X Chromosome

Recent efforts in mapping spatial genome organization have revealed three evocative and conserved structural features of the inactive X in female mammals. First, the chromosomal conformation of the inactive X reveals a loss of topologically associated domains (TADs) present on the active X. Second, the macrosatellite DXZ4 emerges as a singular boundary that suppresses physical interactions between two large TAD-depleted "megadomains." Third, DXZ4 reaches across several megabases to form "superloops" with two other X-linked tandem repeats, FIRRE and ICCE, which also loop to each other. Although all three structural features are conserved across rodents and primates, deletion of mouse and human orthologs of DXZ4 and FIRRE from the inactive X have revealed limited impact on X chromosome inactivation (XCI) and escape in vitro. In contrast, loss of Xist or SMCHD1 have been shown to impair TAD erasure and gene silencing on the inactive X. In this perspective, we summarize these results in the context of new research describing disruption of X-linked tandem repeats in vivo, and discuss their possible molecular roles through the lens of evolutionary conservation and clinical genetics. As a null hypothesis, we consider whether the conservation of some structural features on the inactive X may reflect selection for X-linked tandem repeats on account of necessary cis- and trans-regulatory roles they may play on the active X, rather than the inactive X. Additional hypotheses invoking a role for X-linked tandem repeats on X reactivation, for example in the germline or totipotency, remain to be assessed in multiple developmental models spanning mammalian evolution.

Molecular mechanism for the multiple sclerosis risk variant rs17594362

Multiple sclerosis (MS) is known as an autoimmune demyelinating disease of the central nervous system. However, its cause remains elusive. Given previous studies suggesting that dysfunctional oligodendrocytes (OLs) may trigger MS, we tested whether single nucleotide polymorphisms (SNPs) associated with MS affect OL enhancers, potentially increasing MS risk by dysregulating gene expression of OL lineage cells. We found that two closely spaced OL enhancers, which are 3 Kb apart on chromosome 13, overlap two MS SNPs in linkage disequilibrium-rs17594362 and rs12429256. Our data revealed that the two MS SNPs significantly up-regulate the associated OL enhancers, which we have named as Rgcc-E1 and Rgcc-E2. Analysis of Hi-C data and epigenome editing experiments shows that Rgcc is the primary target of Rgcc-E1 and Rgcc-E2. Collectively, these data indicate that the molecular mechanism of rs17594362 and rs12429256 is to induce Rgcc overexpression by potentiating the enhancer activity of Rgcc-E1 and Rgcc-E2. Importantly, the dosage of the rs17594362/rs12429256 risk allele is positively correlated with the expression level of Rgcc in the human population, confirming our molecular mechanism. Our study also suggests that Rgcc overexpression in OL lineage cells may be a key cellular mechanism of rs17594362 and rs12429256 for MS.

A principled strategy for mapping enhancers to genes

Mapping enhancers to genes is a fundamental goal of modern biology. We have developed an innovative strategy that maps enhancers to genes in a principled manner. We illustrate its power by applying it to Myrf. Despite being a master regulator of oligodendrocytes, oligodendrocyte enhancers governing Myrf expression remain elusive. Since chromatin conformation capture studies have shown that a gene and its enhancer tend to be found in the same topologically associating domain (TAD), we started with the delineation of the Myrf TAD. A genome-wide map of putative oligodendrocyte enhancers uncovered 6 putative oligodendrocyte enhancers in the Myrf TAD, narrowing down the search space for Myrf enhancers from the entire genome to 6 loci in a principled manner. Epigenome editing experiments revealed that two of them govern Myrf expression for oligodendrocyte development. Our new method is simple, principled, and powerful, providing a systematic way to find enhancers that regulate the expression of a gene of interest. Since it can be applied to most cell types, it would greatly facilitate our effort to unravel transcriptional regulatory networks of diverse cell types.

Structural variation in the 3D genome

Structural and quantitative chromosomal rearrangements, collectively referred to as structural variation (SV), contribute to a large extent to the genetic diversity of the human genome and thus are of high relevance for cancer genetics, rare diseases and evolutionary genetics. Recent studies have shown that SVs can not only affect gene dosage but also modulate basic mechanisms of gene regulation. SVs can alter the copy number of regulatory elements or modify the 3D genome by disrupting higher-order chromatin organization such as topologically associating domains. As a result of these position effects, SVs can influence the expression of genes distant from the SV breakpoints, thereby causing disease. The impact of SVs on the 3D genome and on gene expression regulation has to be considered when interpreting the pathogenic potential of these variant types.

Genome sequencing-the dawn of a game-changing era

The development of genome sequencing technologies has revolutionized the biological sciences in ways which could not have been imagined at the time. This article sets out to document the dawning of the age of genomics and to consider the impact of this revolution on biological investigation, our understanding of life, and the relationship between science and society.

Highly interacting regions of the human genome are enriched with enhancers and bound by DNA repair proteins

In specific cases, chromatin clearly forms long-range loops that place distant regulatory elements in close proximity to transcription start sites, but we have limited understanding of many loops identified by Chromosome Conformation Capture (such as Hi-C) analyses. In efforts to elucidate their characteristics and functions, we have identified highly interacting regions (HIRs) using intra-chromosomal Hi-C datasets with a new computational method based on looking at the eigenvector that corresponds to the smallest eigenvalue (here unity). Analysis of these regions using ENCODE data shows that they are in general enriched in bound factors involved in DNA damage repair and have actively transcribed genes. However, both highly transcribed regions as well as transcriptionally inactive regions can form HIRs. The results also indicate that enhancers and super-enhancers in particular form long-range interactions within the same chromosome. The accumulation of DNA repair factors in most identified HIRs suggests that protection from DNA damage in these regions is essential for avoidance of detrimental rearrangements.

Semi-nonparametric modeling of topological domain formation from epigenetic data

Background

Hi-C experiments capturing the 3D genome architecture have led to the discovery of topologically-associated domains (TADs) that form an important part of the 3D genome organization and appear to play a role in gene regulation and other functions. Several histone modifications have been independently associated with TAD formation, but their combinatorial effects on domain formation remain poorly understood at a global scale.

Results

We propose a convex semi-nonparametric approach called nTDP based on Bernstein polynomials to explore the joint effects of histone markers on TAD formation as well as predict TADs solely from the histone data. We find a small subset of modifications to be predictive of TADs across species. By inferring TADs using our trained model, we are able to predict TADs across different species and cell types, without the use of Hi-C data, suggesting their effect is conserved. This work provides the first comprehensive joint model of the effect of histone markers on domain formation.

Conclusions

Our approach, nTDP, can form the basis of a unified, explanatory model of the relationship between epigenetic marks and topological domain structures. It can be used to predict domain boundaries for cell types, species, and conditions for which no Hi-C data is available. The model may also be of use for improving Hi-C-based domain finders.

Analysis of 3D genomic interactions identifies candidate host genes that transposable elements potentially regulate

BACKGROUND

The organization of chromatin in the nucleus plays an essential role in gene regulation. About half of the mammalian genome comprises transposable elements. Given their repetitive nature, reads associated with these elements are generally discarded or randomly distributed among elements of the same type in genome-wide analyses. Thus, it is challenging to identify the activities and properties of individual transposons. As a result, we only have a partial understanding of how transposons contribute to chromatin folding and how they impact gene regulation.

RESULTS

Using PCR and Capture-based chromosome conformation capture (3C) approaches, collectively called 4Tran, we take advantage of the repetitive nature of transposons to capture interactions from multiple copies of endogenous retrovirus (ERVs) in the human and mouse genomes. With 4Tran-PCR, reads are selectively mapped to unique regions in the genome. This enables the identification of transposable element interaction profiles for individual ERV families and integration events specific to particular genomes. With this approach, we demonstrate that transposons engage in long-range intra-chromosomal interactions guided by the separation of chromosomes into A and B compartments as well as topologically associated domains (TADs). In contrast to 4Tran-PCR, Capture-4Tran can uniquely identify both ends of an interaction that involve retroviral repeat sequences, providing a powerful tool for uncovering the individual transposable element insertions that interact with and potentially regulate target genes.

CONCLUSIONS

4Tran provides new insight into the manner in which transposons contribute to chromosome architecture and identifies target genes that transposable elements can potentially control.

Evolutionary stability of topologically associating domains is associated with conserved gene regulation

Background

The human genome is highly organized in the three-dimensional nucleus. Chromosomes fold locally into topologically associating domains (TADs) defined by increased intra-domain chromatin contacts. TADs contribute to gene regulation by restricting chromatin interactions of regulatory sequences, such as enhancers, with their target genes. Disruption of TADs can result in altered gene expression and is associated to genetic diseases and cancers. However, it is not clear to which extent TAD regions are conserved in evolution and whether disruption of TADs by evolutionary rearrangements can alter gene expression.

Results

Here, we hypothesize that TADs represent essential functional units of genomes, which are stable against rearrangements during evolution. We investigate this using whole-genome alignments to identify evolutionary rearrangement breakpoints of different vertebrate species. Rearrangement breakpoints are strongly enriched at TAD boundaries and depleted within TADs across species. Furthermore, using gene expression data across many tissues in mouse and human, we show that genes within TADs have more conserved expression patterns. Disruption of TADs by evolutionary rearrangements is associated with changes in gene expression profiles, consistent with a functional role of TADs in gene expression regulation.

Conclusions

Together, these results indicate that TADs are conserved building blocks of genomes with regulatory functions that are often reshuffled as a whole instead of being disrupted by rearrangements.

Electronic supplementary material

The online version of this article (10.1186/s12915-018-0556-x) contains supplementary material, which is available to authorized users.

Comparing and Contrasting the Effects of Drosophila Condensin II Subunit dCAP-D3 Overexpression and Depletion in Vivo

The Condensin II complex plays important, conserved roles in genome organization throughout the cell cycle and in the regulation of gene expression. Previous studies have linked decreased Condensin II subunit expression with a variety of diseases. Here, we show that elevated levels of Condensin II subunits are detected in somatic cancers. To evaluate potential biological effects of elevated Condensin II levels, we overexpressed the Condensin II subunit, dCAP-D3 in Drosophila melanogaster larval tissues and examined the effects on the mitotic- and interphase-specific functions of Condensin II. Interestingly, while ubiquitous overexpression resulted in pupal lethality, tissue specific overexpression of dCAP-D3 caused formation of nucleoplasmic protein aggregates which slowed mitotic prophase progression, mimicking results observed when dCAP-D3 levels are depleted. Surprisingly, dCAP-D3 aggregate formation resulted in faster transitions from metaphase to anaphase. Overexpressed dCAP-D3 protein failed to precipitate other Condensin II subunits in nondividing tissues, but did cause changes to gene expression which occurred in a manner opposite of what was observed when dCAP-D3 levels were depleted in both dividing and nondividing tissues. Our findings show that altering dCAP-D3 levels in either direction has detrimental effects on mitotic timing, the regulation of gene expression, and organism development. Taken together, these data suggest that the different roles for Condensin II throughout the cell cycle may be independent of each other and/or that dCAP-D3 may possess functions that are separate from those involving its association with the Condensin II complex. If conserved, these findings could have implications for tumors harboring elevated CAP-D3 levels.

Multiple Functional Variants at 13q14 Risk Locus for Osteoporosis Regulate RANKL Expression Through Long-Range Super-Enhancer

RANKL is a key regulator involved in bone metabolism, and a drug target for osteoporosis. The clinical diagnosis and assessment of osteoporosis are mainly based on bone mineral density (BMD). Previous powerful genomewide association studies (GWASs) have identified multiple intergenic single-nucleotide polymorphisms (SNPs) located over 100 kb upstream of RANKL and 65 kb downstream of AKAP11 at 13q14.11 for osteoporosis. Whether these SNPs exert their roles on osteoporosis through RANKL is unknown. In this study, we conducted integrative analyses combining expression quantitative trait locus (eQTL), genomic chromatin interaction (high-throughput chromosome conformation capture [Hi-C]), epigenetic annotation, and a series of functional assays. The eQTL analysis identified six potential functional SNPs (rs9533090, rs9594738, r8001611, rs9533094, rs9533095, and rs9594759) exclusively correlated with RANKL gene expression (p < 0.001) at 13q14.11. Co-localization analyses suggested that eQTL signal for RANKL and BMD-GWAS signal shared the same causal variants. Hi-C analysis and functional annotation further validated that the first five osteoporosis SNPs are located in a super-enhancer region to regulate the expression of RANKL via long-range chromosomal interaction. Particularly, dual-luciferase assay showed that the region harboring rs9533090 in the super-enhancer has the strongest enhancer activity, and rs9533090 is an allele-specific regulatory SNP. Furthermore, deletion of the region harboring rs9533090 using CRISPR/Cas9 genome editing significantly reduced RANKL expression in both mRNA level and protein level. Finally, we found that the rs9533090-C robustly recruits transcription factor NFIC, which efficiently elevates the enhancer activity and increases the RANKL expression. In summary, we provided a feasible method to identify regulatory noncoding SNPs to distally regulate their target gene underlying the pathogenesis of osteoporosis by using bioinformatics data analyses and experimental validation. Our findings would be a potential and promising therapeutic target for precision medicine in osteoporosis. © 2018 American Society for Bone and Mineral Research.

Computational prediction of CTCF/cohesin-based intra-TAD loops that insulate chromatin contacts and gene expression in mouse liver

CTCF and cohesin are key drivers of 3D-nuclear organization, anchoring the megabase-scale Topologically Associating Domains (TADs) that segment the genome. Here, we present and validate a computational method to predict cohesin-and-CTCF binding sites that form intra-TAD DNA loops. The intra-TAD loop anchors identified are structurally indistinguishable from TAD anchors regarding binding partners, sequence conservation, and resistance to cohesin knockdown; further, the intra-TAD loops retain key functional features of TADs, including chromatin contact insulation, blockage of repressive histone mark spread, and ubiquity across tissues. We propose that intra-TAD loops form by the same loop extrusion mechanism as the larger TAD loops, and that their shorter length enables finer regulatory control in restricting enhancer-promoter interactions, which enables selective, high-level expression of gene targets of super-enhancers and genes located within repressive nuclear compartments. These findings elucidate the role of intra-TAD cohesin-and-CTCF binding in nuclear organization associated with widespread insulation of distal enhancer activity.

Structure and function of the Nppa-Nppb cluster locus during heart development and disease

Atrial natriuretic factor and brain natriuretic peptide are two important biomarkers in clinical cardiology. These two natriuretic peptide hormones are encoded by the paralogous genes Nppa and Nppb, which are evolutionary conserved. Both genes are predominantly expressed by the heart muscle during the embryonic and fetal stages, and in particular Nppa expression is strongly reduced in the ventricles after birth. Upon cardiac stress, Nppa and Nppb are strongly upregulated in the ventricular myocardium. Much is known about the molecular and physiological ques inducing Nppa and Nppb expression; however, the transcriptional regulatory mechanisms of the Nppa-Nppb cluster in vivo has proven to be quite complex and is not well understood. In this review, we will provide recent insights into the dynamic and complex regulation of Nppa and Nppb during heart development and hypertrophy, and the association of this gene cluster with the cardiomyocyte-intrinsic program of heart regeneration.

Chromatin 3D – will it make understanding of cancer transformation finally possible?

Despite enormous progress in molecular analysis of cancer cell genomes, the mechanism of tumorigenesis remains unclear. The information present in the genome is not limited to the DNA sequence itself. Indeed, a significant portion of this information is concealed in the spatial structure of chromatin. Ongoing scientific studies that focus on the three-dimensional structure of chromatin raise hopes of arriving at a general explanation of the cancer transformation phenomenon.

4Cin: A computational pipeline for 3D genome modeling and virtual Hi-C analyses from 4C data

The use of 3C-based methods has revealed the importance of the 3D organization of the chromatin for key aspects of genome biology. However, the different caveats of the variants of 3C techniques have limited their scope and the range of scientific fields that could benefit from these approaches. To address these limitations, we present 4Cin, a method to generate 3D models and derive virtual Hi-C (vHi-C) heat maps of genomic loci based on 4C-seq or any kind of 4C-seq-like data, such as those derived from NG Capture-C. 3D genome organization is determined by integrative consideration of the spatial distances derived from as few as four 4C-seq experiments. The 3D models obtained from 4C-seq data, together with their associated vHi-C maps, allow the inference of all chromosomal contacts within a given genomic region, facilitating the identification of Topological Associating Domains (TAD) boundaries. Thus, 4Cin offers a much cheaper, accessible and versatile alternative to other available techniques while providing a comprehensive 3D topological profiling. By studying TAD modifications in genomic structural variants associated to disease phenotypes and performing cross-species evolutionary comparisons of 3D chromatin structures in a quantitative manner, we demonstrate the broad potential and novel range of applications of our method.

Chromatin conformation capture (3C) methods have revealed the importance of the 3D organization of the chromatin, which is key to understand many aspects of genome biology. But each of these methods have their own limitations. Here we present 4Cin, a software that generates 3D models of the chromatin from a small number of 4C-seq experiments, a 3C-based method that provides the frequency of contacts between one fragments and the genome (one vs all). These 3D models are used to infer all chromosomal contacts within a given genomic region (many vs many). The contact maps facilitate the identification of Topological Associating Domains boundaries. Our software offers a much cheaper, accessible and versatile alternative to other available techniques while providing a comprehensive 3D topological profiling. We applied our software to two different loci to study modifications in genomic structural variants associated to disease phenotypes and to compare the chromatin organization in two different species in a quantitative manner.

Dynamic Control of Chromosome Topology and Gene Expression by a Chromatin Modification

The function of chromatin modification in establishing higher-order chromosome structure during gene regulation has been elusive. We dissected the machinery and mechanism underlying the enrichment of histone modification H4K20me1 on hermaphrodite X chromosomes during Caenorhabditis elegans dosage compensation and discovered a key role for H4K20me1 in regulating X-chromosome topology and chromosome-wide gene expression. Structural and functional analysis of the dosage compensation complex (DCC) subunit DPY-21 revealed a novel Jumonji C demethylase subfamily that converts H4K20me2 to H4K20me1 in worms and mammals. Inactivation of demethylase activity in vivo by genome editing eliminated H4K20me1 enrichment on X chromosomes of somatic cells, increased X-linked gene expression, reduced X-chromosome compaction, and disrupted X-chromosome conformation by diminishing the formation of topologically associated domains. H4K20me1 is also enriched on the inactive X of female mice, making our studies directly relevant to mammalian development. Unexpectedly, DPY-21 also associates specifically with autosomes of nematode germ cells in a DCC-independent manner to enrich H4K20me1 and trigger chromosome compaction. Thus, DPY-21 is an adaptable chromatin regulator. Its H4K20me2 demethylase activity can be harnessed during development for distinct biological functions by targeting it to diverse genomic locations through different mechanisms. In both somatic cells and germ cells, H4K20me1 enrichment modulates three-dimensional chromosome architecture, demonstrating the direct link between chromatin modification and higher-order chromosome structure.

Alterations of specific chromatin conformation affect ATRA-induced leukemia cell differentiation

Chromatin conformation plays a key role in regulating gene expression and controlling cell differentiation. However, the whole-genome chromatin conformation changes that occur during leukemia cell differentiation are poorly understood. Here, we characterized the changes in chromatin conformation, histone states, chromatin accessibility, and gene expression using an all-trans retinoic acid (ATRA)-induced HL-60 cell differentiation model. The results showed that the boundaries of topological associated domains (TADs) were stable during differentiation; however, the chromatin conformations within several specific TADs were obviously changed. By combining H3K4me3, H3K27ac, and Hi-C signals, we annotated the differential gene-regulatory chromatin interactions upon ATRA induction. The gains and losses of the gene-regulatory chromatin interactions are significantly correlated with gene expression and chromatin accessibility. Finally, we found that the loss of GATA2 expression and DNA binding are crucial for the differentiation process, and changes in the chromatin structure around the GATA2 regulate its expression upon ATRA induction. This study provided both statistical insights and experimental details regarding the relationship between chromatin conformation changes and transcription regulation during leukemia cell differentiation, and the results suggested that the chromatin conformation is a new type of potential drug target for cancer therapy.

3D Genome Organization Influences the Chromosome Translocation Pattern

Recent imaging, molecular, and computational modeling studies have greatly enhanced our knowledge of how eukaryotic chromosomes are folded in the nuclear space. This work has begun to reveal how 3D genome structure contributes to various DNA-mediated metabolic activities such as replication, transcription, recombination, and repair. Failure of proper DNA repair can lead to the chromosomal translocations observed in human cancers and other diseases. Questions about the role of 3D genome structure in translocation mechanisms have interested scientists for decades. Recent applications of imaging and Chromosome Conformation Capture approaches have clarified the influence of proximal positioning of chromosomal domains and gene loci on the formation of chromosomal translocations. These approaches have revealed the importance of 3D genome structure not only in translocation partner selection, but also in repair efficiency, likelihood of DNA damage, and the biological implications of translocations. This chapter focuses on our current understanding of the role of 3D genome structure in chromosome translocation formation and its potential implications in disease outcome.

Subnuclear distribution of proteins: Links with genome architecture

Metazoan genomes have a hierarchal 3-dimensional (3D) organization scaling from nucleosomes, loops, topologically associating domains (TADs), compartments, to chromosome territories. The 3D organization of genome has been linked with development, differentiation and disease. However, the principles governing the 3D chromatin architecture are just beginning to get unraveled. The nucleus has very high concentration of proteins and these proteins are either diffusely distributed throughout the nucleus, or aggregated in the form of foci/bodies/clusters/speckles or in combination of both. Several evidences suggest that the distribution of proteins within the nuclear space is linked to the organization and function of genome. Here, we describe advances made in understanding the relationship between subnuclear distribution of proteins and genome architecture.