Quantitative differences in TAD border strength underly the TAD hierarchy in Drosophila chromosomes

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.

3D chromatin interactions involving Drosophila insulators are infrequent but preferential and arise before TADs and transcription

In mammals, insulators contribute to the regulation of loop extrusion to organize chromatin into topologically associating domains. In Drosophila the role of insulators in 3D genome organization is, however, under current debate. Here, we addressed this question by combining bioinformatics analysis and multiplexed chromatin imaging. We describe a class of Drosophila insulators enriched at regions forming preferential chromatin interactions genome-wide. Notably, most of these 3D interactions do not involve TAD borders. Multiplexed imaging shows that these interactions occur infrequently, and only rarely involve multiple genomic regions coalescing together in space in single cells. Finally, we show that non-border preferential 3D interactions enriched in this class of insulators are present before TADs and transcription during Drosophila development. Our results are inconsistent with insulators forming stable hubs in single cells, and instead suggest that they fine-tune existing 3D chromatin interactions, providing an additional regulatory layer for transcriptional regulation.

TADfit is a multivariate linear regression model for profiling hierarchical chromatin domains on replicate Hi-C data

Topologically associating domains (TADs) are fundamental building blocks of three dimensional genome, and organized into complex hierarchies. Identifying hierarchical TADs on Hi-C data helps to understand the relationship between genome architectures and gene regulation. Herein we propose TADfit, a multivariate linear regression model for profiling hierarchical chromatin domains, which tries to fit the interaction frequencies in Hi-C contact matrix with and without replicates using all-possible hierarchical TADs, and the significant ones can be determined by the regression coefficients obtained with the help of an online learning solver called Follow-The-Regularized-Leader (FTRL). Beyond the existing methods, TADfit has an ability to handle multiple contact matrix replicates and find partially overlapping TADs on them, which helps to find the comprehensive underlying TADs across replicates from different experiments. The comparative results tell that TADfit has better accuracy and reproducibility, and the hierarchical TADs called by it exhibit a reasonable biological relevance.

The 3D Genome: From Structure to Function

The genome is the most functional part of a cell, and genomic contents are organized in a compact three-dimensional (3D) structure. The genome contains millions of nucleotide bases organized in its proper frame. Rapid development in genome sequencing and advanced microscopy techniques have enabled us to understand the 3D spatial organization of the genome. Chromosome capture methods using a ligation approach and the visualization tool of a 3D genome browser have facilitated detailed exploration of the genome. Topologically associated domains (TADs), lamin-associated domains, CCCTC-binding factor domains, cohesin, and chromatin structures are the prominent identified components that encode the 3D structure of the genome. Although TADs are the major contributors to 3D genome organization, they are absent in Arabidopsis. However, a few research groups have reported the presence of TAD-like structures in the plant kingdom.

3D genome organisation in Drosophila

Ever since Thomas Hunt Morgan's discovery of the chromosomal basis of inheritance by using Drosophila melanogaster as a model organism, the fruit fly has remained an essential model system in studies of genome biology, including chromatin organisation. Very much as in vertebrates, in Drosophila, the genome is organised in territories, compartments and topologically associating domains (TADs). However, these domains might be formed through a slightly different mechanism than in vertebrates due to the presence of a large and potentially redundant set of insulator proteins and the minor role of dCTCF in TAD boundary formation. Here, we review the different levels of chromatin organisation in Drosophila and discuss mechanisms and factors that might be involved in TAD formation. The dynamics of TADs and enhancer-promoter interactions in the context of transcription are covered in the light of currently conflicting results. Finally, we illustrate the value of polymer modelling approaches to infer the principles governing the three-dimensional organisation of the Drosophila genome.

Multi-Scale Organization of the Drosophila melanogaster Genome

Interphase chromatin, despite its appearance, is a highly organized framework of loops and bends. Chromosomes are folded into topologically associating domains, or TADs, and each chromosome and its homolog occupy a distinct territory within the nucleus. In Drosophila, genome organization is exceptional because homologous chromosome pairing is in both germline and somatic tissues, which promote interhomolog interactions such as transvection that can affect gene expression in trans. In this review, we focus on what is known about genome organization in Drosophila and discuss it from TADs to territory. We start by examining intrachromosomal organization at the sub-chromosome level into TADs, followed by a comprehensive analysis of the known proteins that play a key role in TAD formation and boundary establishment. We then zoom out to examine interhomolog interactions such as pairing and transvection that are abundant in Drosophila but rare in other model systems. Finally, we discuss chromosome territories that form within the nucleus, resulting in a complete picture of the multi-scale organization of the Drosophila genome.

Order and stochasticity in the folding of individual Drosophila genomes

Mammalian and Drosophila genomes are partitioned into topologically associating domains (TADs). Although this partitioning has been reported to be functionally relevant, it is unclear whether TADs represent true physical units located at the same genomic positions in each cell nucleus or emerge as an average of numerous alternative chromatin folding patterns in a cell population. Here, we use a single-nucleus Hi-C technique to construct high-resolution Hi-C maps in individual Drosophila genomes. These maps demonstrate chromatin compartmentalization at the megabase scale and partitioning of the genome into non-hierarchical TADs at the scale of 100 kb, which closely resembles the TAD profile in the bulk in situ Hi-C data. Over 40% of TAD boundaries are conserved between individual nuclei and possess a high level of active epigenetic marks. Polymer simulations demonstrate that chromatin folding is best described by the random walk model within TADs and is most suitably approximated by a crumpled globule build of Gaussian blobs at longer distances. We observe prominent cell-to-cell variability in the long-range contacts between either active genome loci or between Polycomb-bound regions, suggesting an important contribution of stochastic processes to the formation of the Drosophila 3D genome.

Genomes are partitioned into topologically associating domains (TADs). Here the authors present single-nucleus Hi-C maps in Drosophila at 10 kb resolution, demonstrating the presence of chromatin compartments in individual nuclei, and partitioning of the genome into non-hierarchical TADs at the scale of 100 kb, which resembles population TAD profiles.

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.

SpectralTAD: an R package for defining a hierarchy of topologically associated domains using spectral clustering

BACKGROUND

The three-dimensional (3D) structure of the genome plays a crucial role in gene expression regulation. Chromatin conformation capture technologies (Hi-C) have revealed that the genome is organized in a hierarchy of topologically associated domains (TADs), sub-TADs, and chromatin loops. Identifying such hierarchical structures is a critical step in understanding genome regulation. Existing tools for TAD calling are frequently sensitive to biases in Hi-C data, depend on tunable parameters, and are computationally inefficient.

METHODS

To address these challenges, we developed a novel sliding window-based spectral clustering framework that uses gaps between consecutive eigenvectors for TAD boundary identification.

RESULTS

Our method, implemented in an R package, SpectralTAD, detects hierarchical, biologically relevant TADs, has automatic parameter selection, is robust to sequencing depth, resolution, and sparsity of Hi-C data. SpectralTAD outperforms four state-of-the-art TAD callers in simulated and experimental settings. We demonstrate that TAD boundaries shared among multiple levels of the TAD hierarchy were more enriched in classical boundary marks and more conserved across cell lines and tissues. In contrast, boundaries of TADs that cannot be split into sub-TADs showed less enrichment and conservation, suggesting their more dynamic role in genome regulation.

CONCLUSION

SpectralTAD is available on Bioconductor, http://bioconductor.org/packages/SpectralTAD/ .

C-InterSecture-a computational tool for interspecies comparison of genome architecture

MOTIVATION

Recent development of Hi-C technique, a biochemical method to study 3D genome architecture, provided large amount of information describing spatial organization of chromosomes in different cell types and species. While multiple tools are available for analysis and comparison of Hi-C data of different cell types, there are almost no resources for systematic interspecies comparison.

RESULTS

To fill this gap, we developed C-InterSecture, a computational pipeline allowing systematic comparison of genome architecture between species. C-InterSecture allows statistical comparison of contact frequencies of individual pairs of loci, as well as interspecies comparison of contacts pattern within defined genomic regions, i.e. topologically associated domains. We employed C-InterSecture to compare mammalian and avian genome organization and showed how evolutionary changes of genomic distance affect 3D architecture of vertebrate's genome.

AVAILABILITY AND IMPLEMENTATION

C-InterSecture is implemented as a collection of python scripts freely available on GitHub repository at https://github.com/NuriddinovMA/C-InterSecture. Jucebox-compatible .hic files produced by C-InterSecture are available at http://genedev.bionet.nsc.ru/site/CIntersecture.html.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Tight DNA-protein complexes isolated from barley seedlings are rich in potential guanine quadruplex sequences

Background

The concept of chromatin domains attached to the nuclear matrix is being revisited, with nucleus described as a set of topologically associating domains. The significance of the tightly bound to DNA proteins (TBP), a protein group that remains attached to DNA after its deproteinization should be also revisited, as the existence of these interactions is in good agreement with the concept of the topologically associating domain. The work aimed to characterize the DNA component of TBP isolated from barley seedlings.

Methods

The tight DNA-protein complexes from the first leaves, coleoptiles, and roots of barley seedlings were isolated by purification with chromatography on nitrocellulose or exhaustive digestion of DNA with DNase I. Cloning and transformation were performed using pMOSBBlue Blunt Ended Cloning Kit. Inserts were amplified by PCR, and sequencing was performed on the MegaBace 1000 Sequencing System. The BLAST search was performed using sequence databases at NCBI, CR-EST, and TREP and Ensembl Plants databases. Comparison to MAR/SAR sequences was performed using http://smartdb.bioinf.med.uni-goettingen.de/cgi-bin/SMARtDB/smar.cgi database. The prediction of G quadruplexes (GQ) was performed with the aid of R-studio library pqsfinder. CD spectra were recorded on a Chirascan CS/3D spectrometer.

Results

Although the barley genome is AT-rich (43% of GC pairs), most DNA fragments associated with TBP were GC-rich (up to 70% in some fractions). Both fractionation procedures yielded a high proportion of CT-motif sequences presented predominantly by the 16-bp CC(TCTCCC)₂ TC fragment present in clones derived from the TBP-bound DNA and absent in free DNA. BLAST analysis revealed alignment with different barley repeats. Some clones, however, aligned with both nuclear and chloroplast structural genes. Alignments with MAR/SAR motifs were very few. The analysis produced by the pqsfinder program revealed numerous potential quadruplex-forming sites in the TBP-bound sequences. A set of oligonucleotides containing sites of possible GQs were designed and ordered. Three of them represented the minus strand of the CT-repeat. Two were derived from sequences of two clones of nitrocellulose retained fraction from leaves and contained GC-rich motifs different from the CT motif. Circular dichroism spectroscopy revealed profound changes in spectra when oligonucleotides were incubated with 100 mM KCl. There was either an increase of positive band in the area of 260 nm or the formation of a positive band at 290 nm. In the former case, changes are typical for parallel G-quadruplexes and, in the latter, 3 + 1 structures.

Discussion

The G-quadruplexes anchor proteins are probably involved in the maintenance of the topologically associated domain structure.

The Many Roles of Cohesin in Drosophila Gene Transcription

The cohesin protein complex mediates sister chromatid cohesion to ensure accurate chromosome segregation, and also influences gene transcription in higher eukaryotes. Modest deficits in cohesin function that do not alter chromosome segregation cause significant birth defects. The mechanisms by which cohesin participates in gene regulation have been studied in Drosophila, revealing that it is involved in gene activation by transcriptional enhancers and epigenetic gene silencing mediated by Polycomb group proteins. Recent studies reveal that early DNA replication origins are important for determining which genes associate with cohesin, and suggest that cohesin at replication origins is important for establishing both sister chromatid cohesion and enhancer-promoter communication.

Principles of genome folding into topologically associating domains

This review discusses the features of TADs across species, and their role in chromosome organization, genome function, and evolution.