The genome-wide multi-layered architecture of chromosome pairing in early Drosophila embryos

The zinc-finger protein Z4 cooperates with condensin II to regulate somatic chromosome pairing and 3D chromatin organization

Chromosome pairing constitutes an important level of genome organization, yet the mechanisms that regulate pairing in somatic cells and the impact on 3D chromatin organization are still poorly understood. Here, we address these questions in Drosophila, an organism with robust somatic pairing. In Drosophila, pairing preferentially occurs at loci consisting of numerous architectural protein binding sites (APBSs), suggesting a role of architectural proteins (APs) in pairing regulation. Amongst these, the anti-pairing function of the condensin II subunit CAP-H2 is well established. However, the factors that regulate CAP-H2 localization and action at APBSs remain largely unknown. Here, we identify two factors that control CAP-H2 occupancy at APBSs and, therefore, regulate pairing. We show that Z4, interacts with CAP-H2 and is required for its localization at APBSs. We also show that hyperosmotic cellular stress induces fast and reversible unpairing in a Z4/CAP-H2 dependent manner. Moreover, by combining the opposite effects of Z4 depletion and osmostress, we show that pairing correlates with the strength of intrachromosomal 3D interactions, such as active (A) compartment interactions, intragenic gene-loops, and polycomb (Pc)-mediated chromatin loops. Altogether, our results reveal new players in CAP-H2-mediated pairing regulation and the intimate interplay between inter-chromosomal and intra-chromosomal 3D interactions.

Building Haplotype-Resolved 3D Genome Maps of Chicken Skeletal Muscle

Haplotype-resolved 3D chromatin architecture related to allelic differences in avian skeletal muscle development has not been addressed so far, although chicken husbandry for meat consumption has been prevalent feature of cultures on every continent for more than thousands of years. Here, high-resolution Hi-C diploid maps (1.2-kb maximum resolution) are generated for skeletal muscle tissues in chicken across three developmental stages (embryonic day 15 to day 30 post-hatching). The sequence features governing spatial arrangement of chromosomes and characterize homolog pairing in the nucleus, are identified. Multi-scale characterization of chromatin reorganization between stages from myogenesis in the fetus to myofiber hypertrophy after hatching show concordant changes in transcriptional regulation by relevant signaling pathways. Further interrogation of parent-of-origin-specific chromatin conformation supported that genomic imprinting is absent in birds. This study also reveals promoter-enhancer interaction (PEI) differences between broiler and layer haplotypes in skeletal muscle development-related genes are related to genetic variation between breeds, however, only a minority of breed-specific variations likely contribute to phenotypic divergence in skeletal muscle potentially via allelic PEI rewiring. Beyond defining the haplotype-specific 3D chromatin architecture in chicken, this study provides a rich resource for investigating allelic regulatory divergence among chicken breeds.

Modeling homologous chromosome recognition via nonspecific interactions

In many organisms, most notably Drosophila, homologous chromosomes associate in somatic cells, a phenomenon known as somatic pairing, which takes place without double strand breaks or strand invasion, thus requiring some other mechanism for homologs to recognize each other. Several studies have suggested a "specific button" model, in which a series of distinct regions in the genome, known as buttons, can associate with each other, mediated by different proteins that bind to these different regions. Here, we use computational modeling to evaluate an alternative "button barcode" model, in which there is only one type of recognition site or adhesion button, present in many copies in the genome, each of which can associate with any of the others with equal affinity. In this model, buttons are nonuniformly distributed, such that alignment of a chromosome with its correct homolog, compared with a nonhomolog, is energetically favored; since to achieve nonhomologous alignment, chromosomes would be required to mechanically deform in order to bring their buttons into mutual register. By simulating randomly generated nonuniform button distributions, many highly effective button barcodes can be easily found, some of which achieve virtually perfect pairing fidelity. This model is consistent with existing literature on the effect of translocations of different sizes on homolog pairing. We conclude that a button barcode model can attain highly specific homolog recognition, comparable to that seen in actual cells undergoing somatic homolog pairing, without the need for specific interactions. This model may have implications for how meiotic pairing is achieved.

A cytological F1 RNAi screen for defects in Drosophila melanogaster female meiosis

Genetic screens for recessive alleles induce mutations, make the mutated chromosomes homozygous, and then assay those homozygotes for the phenotype of interest. When screening for genes required for female meiosis, the phenotype of interest has typically been nondisjunction from chromosome segregation errors. As this requires that mutant females be viable and fertile, any mutants that are lethal or sterile when homozygous cannot be recovered by this approach. To overcome these limitations, we have screened the VALIUM22 collection of RNAi constructs that target germline-expressing genes in a vector optimized for germline expression by driving RNAi with GAL4 under control of a germline-specific promoter (nanos or mat-alpha4). This allowed us to test genes that would be lethal if knocked down in all cells, and by examining unfertilized metaphase-arrested mature oocytes, we could identify defects in sterile females. After screening >1,450 lines of the collection for two different defects (chromosome congression and the hypoxic sequestration of Mps1-GFP to ooplasmic filaments), we obtained multiple hits for both phenotypes, identified novel meiotic phenotypes for genes that had been previously characterized in other processes, and identified the first phenotypes to be associated with several previously uncharacterized genes.

Cooltools: Enabling high-resolution Hi-C analysis in Python

Chromosome conformation capture (3C) technologies reveal the incredible complexity of genome organization. Maps of increasing size, depth, and resolution are now used to probe genome architecture across cell states, types, and organisms. Larger datasets add challenges at each step of computational analysis, from storage and memory constraints to researchers' time; however, analysis tools that meet these increased resource demands have not kept pace. Furthermore, existing tools offer limited support for customizing analysis for specific use cases or new biology. Here we introduce cooltools (https://github.com/open2c/cooltools), a suite of computational tools that enables flexible, scalable, and reproducible analysis of high-resolution contact frequency data. Cooltools leverages the widely-adopted cooler format which handles storage and access for high-resolution datasets. Cooltools provides a paired command line interface (CLI) and Python application programming interface (API), which respectively facilitate workflows on high-performance computing clusters and in interactive analysis environments. In short, cooltools enables the effective use of the latest and largest genome folding datasets.

Pairtools: From sequencing data to chromosome contacts

The field of 3D genome organization produces large amounts of sequencing data from Hi-C and a rapidly-expanding set of other chromosome conformation protocols (3C+). Massive and heterogeneous 3C+ data require high-performance and flexible processing of sequenced reads into contact pairs. To meet these challenges, we present pairtools-a flexible suite of tools for contact extraction from sequencing data. Pairtools provides modular command-line interface (CLI) tools that can be flexibly chained into data processing pipelines. The core operations provided by pairtools are parsing of.sam alignments into Hi-C pairs, sorting and removal of PCR duplicates. In addition, pairtools provides auxiliary tools for building feature-rich 3C+ pipelines, including contact pair manipulation, filtration, and quality control. Benchmarking pairtools against popular 3C+ data pipelines shows advantages of pairtools for high-performance and flexible 3C+ analysis. Finally, pairtools provides protocol-specific tools for restriction-based protocols, haplotype-resolved contacts, and single-cell Hi-C. The combination of CLI tools and tight integration with Python data analysis libraries makes pairtools a versatile foundation for a broad range of 3C+ pipelines.

Haplotype-resolved 3D chromatin architecture of the hybrid pig

In diploid mammals, allele-specific three-dimensional (3D) genome architecture may lead to imbalanced gene expression. Through ultradeep in situ Hi-C sequencing of three representative somatic tissues (liver, skeletal muscle, and brain) from hybrid pigs generated by reciprocal crosses of phenotypically and physiologically divergent Berkshire and Tibetan pigs, we uncover extensive chromatin reorganization between homologous chromosomes across multiple scales. Haplotype-based interrogation of multi-omic data revealed the tissue dependence of 3D chromatin conformation, suggesting that parent-of-origin-specific conformation may drive gene imprinting. We quantify the effects of genetic variations and histone modifications on allelic differences of long-range promoter-enhancer contacts, which likely contribute to the phenotypic differences between the parental pig breeds. We also observe the fine structure of somatically paired homologous chromosomes in the pig genome, which has a functional implication genome-wide. This work illustrates how allele-specific chromatin architecture facilitates concomitant shifts in allele-biased gene expression, as well as the possible consequential phenotypic changes in mammals.

Uncharted territories: Solving the mysteries of male meiosis in flies

The segregation of homologous chromosomes during meiosis typically requires tight end-to-end chromosome pairing. However, in Drosophila spermatogenesis, male flies segregate their chromosomes without classic synaptonemal complex formation and without recombination, instead compartmentalizing homologs into subnuclear domains known as chromosome territories (CTs). How homologs find each other in the nucleus and are separated into CTs has been one of the biggest riddles in chromosome biology. Here, we discuss our current understanding of pairing and CT formation in flies and review recent data on how homologs are linked and partitioned during meiosis in male flies.

Transvection between nonallelic genomic positions in Drosophila

In Drosophila, pairing of maternal and paternal homologous chromosomes can permit trans-interactions between enhancers on one homolog and promoters on another, an example of transvection. Although trans-interactions have been observed at many loci in the Drosophila genome and in other organisms, the parameters that govern enhancer action in trans remain poorly understood. Using a transgenic reporter system, we asked whether enhancers and promoters at nonallelic, but nearby, genomic positions can communication in trans. Using one transgenic insertion carrying the synthetic enhancer GMR and another nearby insertion carrying the hsp70 promoter driving a fluorescent reporter, we show that transgenes separated by 2.6 kb of linear distance can support enhancer action in trans at the 53F8 locus. Furthermore, transvection between the nonallelic insertions can be augmented by a small deletion flanking one insert, likely via changes to the paired configuration of the homologs. Subsequent analyses of other insertions in 53F8 that carry different transgenic sequences demonstrate that the capacity to support transvection between nonallelic sites varies greatly, suggesting that factors beyond the linear distance between insertion sites play an important role. Finally, analysis of transvection between nearby nonallelic sites at other genomic locations shows evidence of position effects, where one locus supported GMR action in trans over a linear distance of over 10 kb, whereas another locus showed no evidence of transvection over a span <200 bp. Overall, our data demonstrate that transvection between nonallelic sites represents a complex interplay between genomic context, interallelic distance, and promoter identity.

Increased long-distance and homo-trans interactions related to H3K27me3 in Arabidopsis hybrids

In plants, the genome structure of hybrids changes compared with their parents, but the effects of these changes in hybrids remain elusive. Comparing reciprocal crosses between Col × C24 and C24 × Col in Arabidopsis using high-throughput chromosome conformation capture assay (Hi-C) analysis, we found that hybrid three-dimensional (3D) chromatin organization had more long-distance interactions relative to parents, and this was mainly located in promoter regions and enriched in genes with heterosis-related pathways. The interactions between euchromatin and heterochromatin were increased, and the compartment strength decreased in hybrids. In compartment domain (CD) boundaries, the distal interactions were more in hybrids than their parents. In the hybrids of CURLY LEAF (clf) mutants clfCol  × clfC24 and clfC24  × clfCol , the heterosis phenotype was damaged, and the long-distance interactions in hybrids were fewer than in their parents with lower H3K27me3. ChIP-seq data revealed higher levels of H3K27me3 in the region adjacent to the CD boundary and the same interactional homo-trans sites in the wild-type (WT) hybrids, which may have led to more long-distance interactions. In addition, the differentially expressed genes (DEGs) located in the boundaries of CDs and loop regions changed obviously in WT, and the functional enrichment for DEGs was different between WT and clf in the long-distance interactions and loop regions. Our findings may therefore propose a new epigenetic explanation of heterosis in the Arabidopsis hybrids and provide new insights into crop breeding and yield increase.

A Cytological F1 RNAi Screen for Defects in Drosophila melanogaster Female Meiosis

Genetic screens for recessive alleles induce mutations, make the mutated chromosomes homozygous, and then assay those homozygotes for the phenotype of interest. When screening for genes required for female meiosis, the phenotype of interest has typically been nondisjunction from chromosome segregation errors. As this requires that mutant females be viable and fertile, any mutants that are lethal or sterile when homozygous cannot be recovered by this approach. To overcome these limitations, our lab has screened the VALIUM22 collection produced by the Harvard TRiP Project, which contains RNAi constructs targeting genes known to be expressed in the germline in a vector optimized for germline expression. By driving RNAi with GAL4 under control of a germline-specific promoter (nanos or mat-alpha4), we can test genes that would be lethal if knocked down in all cells, and by examining unfertilized metaphase-arrested mature oocytes, we can identify defects associated with genes whose knockdown results in sterility or causes other errors besides nondisjunction. We screened this collection to identify genes that disrupt either of two phenotypes when knocked down: the ability of meiotic chromosomes to congress to a single mass at the end of prometaphase, and the sequestration of Mps1-GFP to ooplasmic filaments in response to hypoxia. After screening >1450 lines of the collection, we obtained multiple hits for both phenotypes, identified novel meiotic phenotypes for genes that had been previously characterized in other processes, and identified the first phenotypes to be associated with several previously uncharacterized genes.

A spatial genome aligner for resolving chromatin architectures from multiplexed DNA FISH

Multiplexed fluorescence in situ hybridization (FISH) is a widely used approach for analyzing three-dimensional genome organization, but it is challenging to derive chromosomal conformations from noisy fluorescence signals, and tracing chromatin is not straightforward. Here we report a spatial genome aligner that parses true chromatin signal from noise by aligning signals to a DNA polymer model. Using genomic distances separating imaged loci, our aligner estimates spatial distances expected to separate loci on a polymer in three-dimensional space. Our aligner then evaluates the physical probability observed signals belonging to these loci are connected, thereby tracing chromatin structures. We demonstrate that this spatial genome aligner can efficiently model chromosome architectures from DNA FISH data across multiple scales and be used to predict chromosome ploidies de novo in interphase cells. Reprocessing of previous whole-genome chromosome tracing data with this method indicates the spatial aggregation of sister chromatids in S/G2 phase cells in asynchronous mouse embryonic stem cells and provides evidence for extranumerary chromosomes that remain tightly paired in postmitotic neurons of the adult mouse cortex.

Homologous chromosome recognition via nonspecific interactions

In many organisms, most notably Drosophila, homologous chromosomes in somatic cells associate with each other, a phenomenon known as somatic homolog pairing. Unlike in meiosis, where homology is read out at the level of DNA sequence complementarity, somatic homolog pairing takes place without double strand breaks or strand invasion, thus requiring some other mechanism for homologs to recognize each other. Several studies have suggested a "specific button" model, in which a series of distinct regions in the genome, known as buttons, can associate with each other, presumably mediated by different proteins that bind to these different regions. Here we consider an alternative model, which we term the "button barcode" model, in which there is only one type of recognition site or adhesion button, present in many copies in the genome, each of which can associate with any of the others with equal affinity. An important component of this model is that the buttons are non-uniformly distributed, such that alignment of a chromosome with its correct homolog, compared with a non-homolog, is energetically favored; since to achieve nonhomologous alignment, chromosomes would be required to mechanically deform in order to bring their buttons into mutual register. We investigated several types of barcodes and examined their effect on pairing fidelity. We found that high fidelity homolog recognition can be achieved by arranging chromosome pairing buttons according to an actual industrial barcode used for warehouse sorting. By simulating randomly generated non-uniform button distributions, many highly effective button barcodes can be easily found, some of which achieve virtually perfect pairing fidelity. This model is consistent with existing literature on the effect of translocations of different sizes on homolog pairing. We conclude that a button barcode model can attain highly specific homolog recognition, comparable to that seen in actual cells undergoing somatic homolog pairing, without the need for specific interactions. This model may have implications for how meiotic pairing is achieved.

Somatic chromosome pairing has a determinant impact on 3D chromatin organization

In the nucleus, chromatin is intricately structured into multiple layers of 3D organization important for genome activity. How distinct layers influence each other is not well understood. In particular, the contribution of chromosome pairing to 3D chromatin organization has been largely neglected. Here, we address this question in Drosophila, an organism that shows robust chromosome pairing in interphasic somatic cells. The extent of chromosome pairing depends on the balance between pairing and anti-pairing factors, with the anti-pairing activity of the CAP-H2 condensin II subunit being the best documented. Here, we identify the zinc-finger protein Z4 as a strong anti-pairer that interacts with and mediates the chromatin binding of CAP-H2. We also report that hyperosmotic cellular stress induces fast and reversible chromosome unpairing that depends on Z4/CAP-H2. And, most important, by combining Z4 depletion and osmostress, we show that chromosome pairing reinforces intrachromosomal 3D interactions. On the one hand, pairing facilitates RNAPII occupancy that correlates with enhanced intragenic gene-loop interactions. In addition, acting at a distance, pairing reinforces chromatin-loop interactions mediated by Polycomb (Pc). In contrast, chromosome pairing does not affect which genomic intervals segregate to active (A) and inactive (B) compartments, with only minimal effects on the strength of A-A compartmental interactions. Altogether, our results unveil the intimate interplay between inter-chromosomal and intra-chromosomal 3D interactions, unraveling the interwoven relationship between different layers of chromatin organization and the essential contribution of chromosome pairing.

C-DNA may facilitate homologous DNA pairing

Recombination-independent homologous pairing represents a prominent yet largely enigmatic feature of chromosome biology. As suggested by studies in the fungus Neurospora crassa, this process may be based on the direct pairing of homologous DNA molecules. Theoretical search for the DNA structures consistent with those genetic results has led to an all-atom model in which the B-DNA conformation of the paired double helices is strongly shifted toward C-DNA. Coincidentally, C-DNA also features a very shallow major groove that could permit initial homologous contacts without atom-atom clashes. The hereby conjectured role of C-DNA in homologous pairing should encourage the efforts to discover its biological functions and may also clarify the mechanism of recombination-independent recognition of DNA homology.

Pairtools: from sequencing data to chromosome contacts

The field of 3D genome organization produces large amounts of sequencing data from Hi-C and a rapidly-expanding set of other chromosome conformation protocols (3C+). Massive and heterogeneous 3C+ data require high-performance and flexible processing of sequenced reads into contact pairs. To meet these challenges, we present pairtools - a flexible suite of tools for contact extraction from sequencing data. Pairtools provides modular command-line interface (CLI) tools that can be flexibly chained into data processing pipelines. Pairtools provides both crucial core tools as well as auxiliary tools for building feature-rich 3C+ pipelines, including contact pair manipulation, filtration, and quality control. Benchmarking pairtools against popular 3C+ data pipelines shows advantages of pairtools for high-performance and flexible 3C+ analysis. Finally, pairtools provides protocol-specific tools for multi-way contacts, haplotype-resolved contacts, and single-cell Hi-C. The combination of CLI tools and tight integration with Python data analysis libraries makes pairtools a versatile foundation for a broad range of 3C+ pipelines.

A global high-density chromatin interaction network reveals functional long-range and trans-chromosomal relationships

BACKGROUND

Chromatin contacts are essential for gene-expression regulation; however, obtaining a high-resolution genome-wide chromatin contact map is still prohibitively expensive owing to large genome sizes and the quadratic scale of pairwise data. Chromosome conformation capture (3C)-based methods such as Hi-C have been extensively used to obtain chromatin contacts. However, since the sparsity of these maps increases with an increase in genomic distance between contacts, long-range or trans-chromatin contacts are especially challenging to sample.

RESULTS

Here, we create a high-density reference genome-wide chromatin contact map using a meta-analytic approach. We integrate 3600 human, 6700 mouse, and 500 fly Hi-C experiments to create species-specific meta-Hi-C chromatin contact maps with 304 billion, 193 billion, and 19 billion contacts in respective species. We validate that meta-Hi-C contact maps are uniquely powered to capture functional chromatin contacts in both cis and trans. We find that while individual dataset Hi-C networks are largely unable to predict any long-range coexpression (median 0.54 AUC), meta-Hi-C networks perform comparably in both cis and trans (0.65 AUC vs 0.64 AUC). Similarly, for long-range expression quantitative trait loci (eQTL), meta-Hi-C contacts outperform all individual Hi-C experiments, providing an improvement over the conventionally used linear genomic distance-based association. Assessing between species, we find patterns of chromatin contact conservation in both cis and trans and strong associations with coexpression even in species for which Hi-C data is lacking.

CONCLUSIONS

We have generated an integrated chromatin interaction network which complements a large number of methodological and analytic approaches focused on improved specificity or interpretation. This high-depth "super-experiment" is surprisingly powerful in capturing long-range functional relationships of chromatin interactions, which are now able to predict coexpression, eQTLs, and cross-species relationships. The meta-Hi-C networks are available at https://labshare.cshl.edu/shares/gillislab/resource/HiC/ .

The 3D genome landscape: Diverse chromosomal interactions and their functional implications

Genome organization includes contacts both within a single chromosome and between distinct chromosomes. Thus, regulatory organization in the nucleus may include interplay of these two types of chromosomal interactions with genome activity. Emerging advances in omics and single-cell imaging technologies have allowed new insights into chromosomal contacts, including those of homologs and sister chromatids, and their significance to genome function. In this review, we highlight recent studies in this field and discuss their impact on understanding the principles of chromosome organization and associated functional implications in diverse cellular processes. Specifically, we describe the contributions of intra-chromosomal, inter-homolog, and inter-sister chromatid contacts to genome organization and gene expression.

Application of the 3C Method to Study the Developmental Genes in Drosophila Larvae

A transition from one developmental stage to another is accompanied by activation of developmental programs and corresponding gene ensembles. Changes in the spatial conformation of the corresponding loci are associated with this activation and can be investigated with the help of the Chromosome Conformation Capture (3C) methodology. Application of 3C to specific developmental stages is a sophisticated task. Here, we describe the use of the 3C method to study the spatial organization of developmental loci in Drosophila larvae. We critically analyzed the existing protocols and offered our own solutions and the optimized protocol to overcome limitations. To demonstrate the efficiency of our procedure, we studied the spatial organization of the developmental locus Dad in 3rd instar Drosophila larvae. Differences in locus conformation were found between embryonic cells and living wild-type larvae. We also observed the establishment of novel regulatory interactions in the presence of an adjacent transgene upon activation of its expression in larvae. Our work fills the gap in the application of the 3C method to Drosophila larvae and provides a useful guide for establishing 3C on an animal model.

Cas9/Nickase-induced allelic conversion by homologous chromosome-templated repair in Drosophila somatic cells

Repair of double-strand breaks (DSBs) in somatic cells is primarily accomplished by error-prone nonhomologous end joining and less frequently by precise homology-directed repair preferentially using the sister chromatid as a template. Here, a Drosophila system performs efficient somatic repair of both DSBs and single-strand breaks (SSBs) using intact sequences from the homologous chromosome in a process we refer to as homologous chromosome-templated repair (HTR). Unexpectedly, HTR-mediated allelic conversion at the white locus was more efficient (40 to 65%) in response to SSBs induced by Cas9-derived nickases D10A or H840A than to DSBs induced by fully active Cas9 (20 to 30%). Repair phenotypes elicited by Nickase versus Cas9 differ in both developmental timing (late versus early stages, respectively) and the production of undesired mutagenic events (rare versus frequent). Nickase-mediated HTR represents an efficient and unanticipated mechanism for allelic correction, with far-reaching potential applications in the field of gene editing.

Sister chromatid-sensitive Hi-C to map the conformation of replicated genomes

Chromosome conformation capture (Hi-C) techniques map the 3D organization of entire genomes. How sister chromatids fold in replicated chromosomes, however, cannot be determined with conventional Hi-C because of the identical DNA sequences of sister chromatids. Here, we present a protocol for sister chromatid-sensitive Hi-C (scsHi-C) that enables the distinction of DNA contacts within individual sister chromatids (cis sister contacts) from those between sister chromatids (trans sister contacts), thereby allowing investigation of the organization of replicated genomes. scsHi-C is based on live-cell labeling of nascent DNA by the synthetic nucleoside 4-thio-thymidine (4sT), which incorporates into a distinct DNA strand on each sister chromatid because of semi-conservative DNA replication. After purification of genomic DNA and in situ Hi-C library preparation, 4sT is chemically converted into 5-methyl-cytosine in the presence of OsO₄/NH₄Cl to introduce T-to-C signature point mutations on 4sT-labeled DNA. The Hi-C library is then sequenced, and ligated fragments are assigned to sister chromatids on the basis of strand orientation and the presence of signature mutations. The ensemble of scsHi-C contacts thereby represents genome-wide contact probabilities within and across sister chromatids. scsHi-C can be completed in 2 weeks, has been successfully applied in HeLa cells and can potentially be established for any cell type that allows proper cell cycle synchronization and incorporation of sufficient amounts of 4sT. The genome-wide maps of replicated chromosomes detected by scsHi-C enable investigation of the molecular mechanisms shaping sister chromatid topologies and the relevance of sister chromatid conformation in crucial processes like DNA repair, mitotic chromosome formation and potentially other biological processes.

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.

Understanding three-dimensional chromatin organization in diploid genomes

The three-dimensional (3D) organization of chromatin in the nucleus of diploid eukaryotic organisms has fascinated biologists for many years. Despite major progress in chromatin conformation studies, current knowledge regarding the spatial organization of diploid (maternal and paternal) genomes is still limited. Recent advances in Hi-C technology and data processing approaches have enabled construction of diploid Hi-C contact maps. These maps greatly accelerated the pace of novel discoveries in haplotype-resolved 3D genome studies, revealing the role of allele biased chromatin conformation in transcriptional regulation. Here, we review emerging concepts and haplotype phasing strategies of Hi-C data in 3D diploid genome studies. We discuss new insights on homologous chromosomal organization and the interplay between allelic biased chromatin architecture and several nuclear functions, explaining how haplotype-resolved Hi-C technologies have been used to resolve important biological questions.

Live imaging and biophysical modeling support a button-based mechanism of somatic homolog pairing in Drosophila

Three-dimensional eukaryotic genome organization provides the structural basis for gene regulation. In Drosophila melanogaster, genome folding is characterized by somatic homolog pairing, where homologous chromosomes are intimately paired from end to end; however, how homologs identify one another and pair has remained mysterious. Recently, this process has been proposed to be driven by specifically interacting 'buttons' encoded along chromosomes. Here, we turned this hypothesis into a quantitative biophysical model to demonstrate that a button-based mechanism can lead to chromosome-wide pairing. We tested our model using live-imaging measurements of chromosomal loci tagged with the MS2 and PP7 nascent RNA labeling systems. We show solid agreement between model predictions and experiments in the pairing dynamics of individual homologous loci. Our results strongly support a button-based mechanism of somatic homolog pairing in Drosophila and provide a theoretical framework for revealing the molecular identity and regulation of buttons.

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.

Conformation of sister chromatids in the replicated human genome

The three-dimensional organization of the genome supports regulated gene expression, recombination, DNA repair, and chromosome segregation during mitosis. Chromosome conformation capture (Hi-C)1,2 analysis has revealed a complex genomic landscape of internal chromosomal structures in vertebrate cells3-7, but the identical sequence of sister chromatids has made it difficult to determine how they topologically interact in replicated chromosomes. Here we describe sister-chromatid-sensitive Hi-C (scsHi-C), which is based on labelling of nascent DNA with 4-thio-thymidine and nucleoside conversion chemistry. Genome-wide conformation maps of human chromosomes reveal that sister-chromatid pairs interact most frequently at the boundaries of topologically associating domains (TADs). Continuous loading of a dynamic cohesin pool separates sister-chromatid pairs inside TADs and is required to focus sister-chromatid contacts at TAD boundaries. We identified a subset of TADs that are overall highly paired and are characterized by facultative heterochromatin and insulated topological domains that form separately within individual sister chromatids. The rich pattern of sister-chromatid topologies and our scsHi-C technology will make it possible to investigate how physical interactions between identical DNA molecules contribute to DNA repair, gene expression, chromosome segregation, and potentially other biological processes.

Phase separation drives pairing of homologous chromosomes

Pairing of homologous chromosomes is crucial for ensuring accurate segregation of chromosomes during meiosis. Molecular mechanisms of homologous chromosome pairing in meiosis have been extensively studied in the fission yeast Schizosaccharomyces pombe. In this organism, meiosis-specific noncoding RNA transcribed from specific genes accumulates at the respective gene loci, and chromosome-associated RNA-protein complexes mediate meiotic pairing of homologous loci through phase separation. Pairing of homologous chromosomes also occurs in somatic diploid cells in certain situations. For example, somatic pairing of homologous chromosomes occurs during the early embryogenesis in diptera, and relies on the transcription-associated chromatin architecture. Earlier models also suggest that transcription factories along the chromosome mediate pairing of homologous chromosomes in plants. These studies suggest that RNA bodies formed on chromosomes mediate the pairing of homologous chromosomes. This review summarizes lessons from S. pombe to provide general insights into mechanisms of homologous chromosome pairing mediated by phase separation of chromosome-associated RNA-protein complexes.

Pericentromeric heterochromatin is hierarchically organized and spatially contacts H3K9me2 islands in euchromatin

Membraneless pericentromeric heterochromatin (PCH) domains play vital roles in chromosome dynamics and genome stability. However, our current understanding of 3D genome organization does not include PCH domains because of technical challenges associated with repetitive sequences enriched in PCH genomic regions. We investigated the 3D architecture of Drosophila melanogaster PCH domains and their spatial associations with the euchromatic genome by developing a novel analysis method that incorporates genome-wide Hi-C reads originating from PCH DNA. Combined with cytogenetic analysis, we reveal a hierarchical organization of the PCH domains into distinct “territories.” Strikingly, H3K9me2-enriched regions embedded in the euchromatic genome show prevalent 3D interactions with the PCH domain. These spatial contacts require H3K9me2 enrichment, are likely mediated by liquid-liquid phase separation, and may influence organismal fitness. Our findings have important implications for how PCH architecture influences the function and evolution of both repetitive heterochromatin and the gene-rich euchromatin.

The three dimensional (3D) organization of genomes in cell nuclei can influence a wide variety of genome functions. However, most of our understanding of this critical architecture has been limited to the gene-rich euchromatin, and largely ignores the gene-poor and repeat-rich pericentromeric heterochromatin, or PCH. PCH comprises a large part of most eukaryotic genomes, forms 3D membraneless PCH domains in nuclei, and plays a vital role in chromosome dynamics and genome stability. In this study, we developed a new method that overcomes the technical challenges imposed by the highly repetitive PCH DNA, and generated a comprehensive picture of its 3D organization. Combined with image analyses, we reveal a hierarchical organization of the PCH domains. Surprisingly, we showed that distant euchromatic regions enriched for repressive epigenetic marks also dynamically interact with the main PCH domains. These 3D interactions are likely mediated by liquid-liquid phase separation (similar to how oil and vinegar separate in salad dressing) and the resulting liquid-like fusion events, and can influence the fitness of individuals. Our discoveries have strong implications for how seemingly “junk” DNA could impact functions in the gene-rich euchromatin.

Three-dimensional Organization of Polytene Chromosomes in Somatic and Germline Tissues of Malaria Mosquitoes

Spatial organization of chromosome territories and interactions between interphase chromosomes themselves, as well as with the nuclear periphery, play important roles in epigenetic regulation of the genome function. However, the interplay between inter-chromosomal contacts and chromosome-nuclear envelope attachments in an organism’s development is not well-understood. To address this question, we conducted microscopic analyses of the three-dimensional chromosome organization in malaria mosquitoes. We employed multi-colored oligonucleotide painting probes, spaced 1 Mb apart along the euchromatin, to quantitatively study chromosome territories in larval salivary gland cells and adult ovarian nurse cells of Anopheles gambiae, An. coluzzii, and An. merus. We found that the X chromosome territory has a significantly smaller volume and is more compact than the autosomal arm territories. The number of inter-chromosomal, and the percentage of the chromosome–nuclear envelope, contacts were conserved among the species within the same cell type. However, the percentage of chromosome regions located at the nuclear periphery was typically higher, while the number of inter-chromosomal contacts was lower, in salivary gland cells than in ovarian nurse cells. The inverse correlation was considerably stronger for the autosomes. Consistent with previous theoretical arguments, our data indicate that, at the genome-wide level, there is an inverse relationship between chromosome-nuclear envelope attachments and chromosome–chromosome interactions, which is a key feature of the cell type-specific nuclear architecture.

Fine-Mapping Complex Inversion Breakpoints and Investigating Somatic Pairing in the Anopheles gambiae Species Complex Using Proximity-Ligation Sequencing

Chromosomal inversions are fundamental drivers of genome evolution. In the main Afrotropical malaria vector species, belonging to the Anopheles gambiae species complex, inversions play an important role in local adaptation and have a rich history of cytological study. Despite the importance and ubiquity of some chromosomal inversions across the species complex, inversion breakpoints are often challenging to map molecularly due to the presence of large repetitive regions. Here, we develop an approach that uses Hi-C sequencing data to molecularly fine-map the breakpoints of inversions. We demonstrate that this approach is robust and likely to be widely applicable for both identification and fine-mapping inversion breakpoints in species whose inversions have heretofore been challenging to characterize. We apply our method to interrogate the previously unknown inversion breakpoints of 2Rbc and 2Rd in An. coluzzii We found that inversion breakpoints occur in large repetitive regions, and, strikingly, among three inversions analyzed, two breakpoints appear to be reused in two separate inversions. These breakpoint-adjacent regions are strongly enriched for the presence of a 30 bp satellite repeat sequence. Because low frequency inversion breakpoints are not correlated with genomic regions containing this satellite, we suggest that interrupting this particular repeat may result in arrangements with higher relative fitness. Additionally, we use heterozygous individuals to quantitatively investigate the impacts of somatic pairing in the regions immediately surrounding inversion breakpoints. Finally, we discuss important considerations for possible applications of this approach for inversion breakpoint identification in a range of organisms.

Global chromatin conformation differences in the Drosophila dosage compensated chromosome X

In Drosophila melanogaster the single male chromosome X undergoes an average twofold transcriptional upregulation for balancing the transcriptional output between sexes. Previous literature hypothesised that a global change in chromosome structure may accompany this process. However, recent studies based on Hi-C failed to detect these differences. Here we show that global conformational differences are specifically present in the male chromosome X and detectable using Hi-C data on sex-sorted embryos, as well as male and female cell lines, by leveraging custom data analysis solutions. We find the male chromosome X has more mid-/long-range interactions. We also identify differences at structural domain boundaries containing BEAF-32 in conjunction with CP190 or Chromator. Weakening of these domain boundaries in male chromosome X co-localizes with the binding of the dosage compensation complex and its co-factor CLAMP, reported to enhance chromatin accessibility. Together, our data strongly indicate that chromosome X dosage compensation affects global chromosome structure.

In Drosophila, dosage compensation involves a twofold transcriptional upregulation of the single male chromosome X. Here the authors show that global conformational differences are specifically present in the male X chromosome and detectable using Hi-C data, indicating that dosage compensation affects global chromosome structure.

Direct Homologous dsDNA-dsDNA Pairing: How, Where, and Why?

The ability of homologous chromosomes (or selected chromosomal loci) to pair specifically in the apparent absence of DNA breakage and recombination represents a prominent feature of eukaryotic biology. The mechanism of homology recognition at the basis of such recombination-independent pairing has remained elusive. A number of studies have supported the idea that sequence homology can be sensed between intact DNA double helices in vivo. In particular, recent analyses of the two silencing phenomena in fungi, known as "repeat-induced point mutation" (RIP) and "meiotic silencing by unpaired DNA" (MSUD), have provided genetic evidence for the existence of the direct homologous dsDNA-dsDNA pairing. Both RIP and MSUD likely rely on the same search strategy, by which dsDNA segments are matched as arrays of interspersed base-pair triplets. This process is general and very efficient, yet it proceeds normally without the RecA/Rad51/Dmc1 proteins. Further studies of RIP and MSUD may yield surprising insights into the function of DNA in the cell.

Characterization of Button Loci that Promote Homologous Chromosome Pairing and Cell-Type-Specific Interchromosomal Gene Regulation

Homologous chromosomes colocalize to regulate gene expression in processes including genomic imprinting, X-inactivation, and transvection. In Drosophila, homologous chromosomes pair throughout development, promoting transvection. The "button" model of pairing proposes that specific regions along chromosomes pair with high affinity. Here, we identify buttons interspersed across the fly genome that pair with their homologous sequences, even when relocated to multiple positions in the genome. A majority of transgenes that span a full topologically associating domain (TAD) function as buttons, but not all buttons contain TADs. Additionally, buttons are enriched for insulator protein clusters. Fragments of buttons do not pair, suggesting that combinations of elements within a button are required for pairing. Pairing is necessary but not sufficient for transvection. Additionally, pairing and transvection are stronger in some cell types than in others, suggesting that pairing strength regulates transvection efficiency between cell types. Thus, buttons pair homologous chromosomes to facilitate cell-type-specific interchromosomal gene regulation.

Two major mechanisms of chromosome organization

The spatial organization of chromosomes has long been connected to their polymeric nature and is believed to be important for their biological functions, including the control of interactions between genomic elements, the maintenance of genetic information, and the compaction and safe transfer of chromosomes to cellular progeny. chromosome conformation capture techniques, particularly Hi-C, have provided a comprehensive picture of spatial chromosome organization and revealed new features and elements of chromosome folding. Furthermore, recent advances in microscopy have made it possible to obtain distance maps for extensive regions of chromosomes (Bintu et al., 2018; Nir et al., 2018 [2^••,3]), providing information complementary to, and in excellent agreement with, Hi-C maps. Not only has the resolution of both techniques advanced significantly, but new perturbation data generated in the last two years have led to the identification of molecular mechanisms behind large-scale genome organization. Two major mechanisms that have been proposed to govern chromosome organization are (i) the active (ATP-dependent) process of loop extrusion by Structural Maintenance of Chromosomes (SMC) complexes, and (ii) the spatial compartmentalization of the genome, which is likely mediated by affinity interactions between heterochromatic regions (Falk et al., 2019 [76^••]) rather than by ATP-dependent processes. Here, we review existing evidence that these two processes operate together to fold chromosomes in interphase and that loop extrusion alone drives mitotic compaction. We discuss possible implications of these mechanisms for chromosome function.

The Role of Insulators in Transgene Transvection in Drosophila

Transvection is the phenomenon where a transcriptional enhancer activates a promoter located on the homologous chromosome. It has been amply documented in Drosophila where homologs are closely paired in most, if not all, somatic nuclei, but it has been known to rarely occur in mammals as well. We have taken advantage of site-directed transgenesis to insert reporter constructs into the same genetic locus in Drosophila and have evaluated their ability to engage in transvection by testing many heterozygous combinations. We find that transvection requires the presence of an insulator element on both homologs. Homotypic trans-interactions between four different insulators can support transvection: the gypsy insulator (GI), Wari, Fab-8 and 1A2; GI and Fab-8 are more effective than Wari or 1A2 We show that, in the presence of insulators, transvection displays the characteristics that have been previously described: it requires homolog pairing, but can happen at any of several loci in the genome; a solitary enhancer confronted with an enhancerless reporter is sufficient to drive transcription; it is weaker than the action of the same enhancer-promoter pair in cis, and it is further suppressed by cis-promoter competition. Though necessary, the presence of homotypic insulators is not sufficient for transvection; their position, number and orientation matters. A single GI adjacent to both enhancer and promoter is the optimal configuration. The identity of enhancers and promoters in the vicinity of a trans-interacting insulator pair is also important, indicative of complex insulator-enhancer-promoter interactions.