Search

Genome organization involves cis and trans chromosomal interactions, both implicated in gene regulation, development, and disease. Here, we focus on trans interactions in Drosophila, where homologous chromosomes are paired in somatic cells from embryogenesis through adulthood. We first address long-standing questions regarding the structure of embryonic homolog pairing and, to this end, develop a haplotype-resolved Hi-C approach to minimize homolog misassignment and thus robustly distinguish trans-homolog from cis contacts. This computational approach, which we call Ohm, reveals pairing to be surprisingly structured genome-wide, with trans-homolog domains, compartments, and interaction peaks, many coinciding with analogous cis features. We also find a significant genome-wide correlation between pairing, transcription during zygotic genome activation, and binding of the pioneer factor Zelda. Our findings reveal a complex, highly structured organization underlying homolog pairing, first discovered a century ago in Drosophila. Finally, we demonstrate the versatility of our haplotype-resolved approach by applying it to mammalian embryos.

Compared to noncoding RNAs (ncRNAs), such as rRNAs and ribozymes, for which high-resolution structures abound, little is known about the tertiary structures of mRNAs. In eukaryotic cells, newly made mRNAs are packaged with proteins in highly compacted mRNA particles (mRNPs), but the manner of this mRNA compaction is unknown. Here, we developed and implemented RIPPLiT (RNA immunoprecipitation and proximity ligation in tandem), a transcriptome-wide method for probing the 3D conformations of RNAs stably associated with defined proteins, in this case, exon junction complex (EJC) core factors. EJCs multimerize with other mRNP components to form megadalton-sized complexes that protect large swaths of newly synthesized mRNAs from endonuclease digestion. Unlike ncRNPs, wherein strong locus-specific structures predominate, mRNPs behave more like flexible polymers. Polymer analysis of proximity ligation data for hundreds of mRNA species demonstrates that nascent and pre-translational mammalian mRNAs are compacted by their associated proteins into linear rod-like structures. Metkar et al. developed a new toolkit, a biochemical approach, RIPPLiT and a bioinformatics suite to capture and analyze higher-order organization of RNPs. They were able to enrich for interactions within mRNAs and further, using polymer analysis on 100 s of mRNAs, identified a unifying principle for mRNP packaging.

Fertilization triggers assembly of higher-order chromatin structure from a condensed maternal and a naïve paternal genome to generate a totipotent embryo. Chromatin loops and domains have been detected in mouse zygotes by single-nucleus Hi-C (snHi-C), but not bulk Hi-C. It is therefore unclear when and how embryonic chromatin conformations are assembled. Here, we investigated whether a mechanism of cohesin-dependent loop extrusion generates higher-order chromatin structures within the one-cell embryo. Using snHi-C of mouse knockout embryos, we demonstrate that the zygotic genome folds into loops and domains that critically depend on Scc1-cohesin and that are regulated in size and linear density by Wapl. Remarkably, we discovered distinct effects on maternal and paternal chromatin loop sizes, likely reflecting differences in loop extrusion dynamics and epigenetic reprogramming. Dynamic polymer models of chromosomes reproduce changes in snHi-C, suggesting a mechanism where cohesin locally compacts chromatin by active loop extrusion, whose processivity is controlled by Wapl. Our simulations and experimental data provide evidence that cohesin-dependent loop extrusion organizes mammalian genomes over multiple scales from the one-cell embryo onward., National Institutes of Health (U.S.) (Grant R01GM114190), National Institutes of Health (U.S.) (Grant U54DK107980), National Science Foundation (U.S.) (Grant 1504942)

Mammalian chromatin is spatially organized at many scales showing two prominent features in interphase: (i) alternating regions (1–10 Mb) of active and inactive chromatin that spatially segregate into different compartments, and (II) domains (<1 Mb), that is, regions that preferentially interact internally [topologically associating domains (TADs)] and are central to gene regulation. There is growing evidence that TADs are formed by active extrusion of chromatin loops by cohesin, whereas compartmentalization is established according to local chromatin states. Here, we use polymer simulations to examine how loop extrusion and compartmental segregation work collectively and potentially interfere in shaping global chromosome organization. A model with differential attraction between euchromatin and heterochromatin leads to phase separation and reproduces compartmentalization as observed in Hi-C. Loop extrusion, essential for TAD formation, in turn, interferes with compartmentalization. Our integrated model faithfully reproduces Hi-C data from puzzling experimental observations where altering loop extrusion also led to changes in compartmentalization. Specifically, depletion of chromatin-associated cohesin reduced TADs and revealed finer compartments, while increased processivity of cohesin strengthened large TADs and reduced compartmentalization; and depletion of the TAD boundary protein CTCF weakened TADs while leaving compartments unaffected. We reveal that these experimental perturbations are special cases of a general polymer phenomenon of active mixing by loop extrusion. Our results suggest that chromatin organization on the megabase scale emerges from competition of nonequilibrium active loop extrusion and epigenetically defined compartment structure. Keywords: chromatin; genome architecture; Hi-C; polymer physics; active matter, National Science Foundation (U.S.) (Grant 1504942)

Metazoan genomes are spatially organized at multiple scales, from packaging of DNA around individual nucleosomes to segregation of whole chromosomes into distinct territories. At the intermediate scale of kilobases to megabases, which encompasses the sizes of genes, gene clusters and regulatory domains, the three-dimensional (3D) organization of DNA is implicated in multiple gene regulatory mechanisms, but understanding this organization remains a challenge. At this scale, the genome is partitioned into domains of different epigenetic states that are essential for regulating gene expression. Here we investigate the 3D organization of chromatin in different epigenetic states using super-resolution imaging. We classified genomic domains in Drosophila cells into transcriptionally active, inactive or Polycomb-repressed states, and observed distinct chromatin organizations for each state. All three types of chromatin domains exhibit power-law scaling between their physical sizes in 3D and their domain lengths, but each type has a distinct scaling exponent. Polycomb-repressed domains show the densest packing and most intriguing chromatin folding behaviour, in which chromatin packing density increases with domain length. Distinct from the self-similar organization displayed by transcriptionally active and inactive chromatin, the Polycomb-repressed domains are characterized by a high degree of chromatin intermixing within the domain. Moreover, compared to inactive domains, Polycomb-repressed domains spatially exclude neighbouring active chromatin to a much stronger degree. Computational modelling and knockdown experiments suggest that reversible chromatin interactions mediated by Polycomb-group proteins play an important role in these unique packaging properties of the repressed chromatin. Taken together, our super-resolution images reveal distinct chromatin packaging for different epigenetic states at the kilobase-to-megabase scale, a length scale that is directly relevant

Chromatin is reprogrammed after fertilization to produce a totipotent zygote with the potential to generate a new organism. The maternal genome inherited from the oocyte and the paternal genome provided by sperm coexist as separate haploid nuclei in the zygote. How these two epigenetically distinct genomes are spatially organized is poorly understood. Existing chromosome conformation capture-based methods are not applicable to oocytes and zygotes owing to a paucity of material. To study three-dimensional chromatin organization in rare cell types, we developed a single-nucleus Hi-C (high-resolution chromosome conformation capture) protocol that provides greater than tenfold more contacts per cell than the previous method. Here we show that chromatin architecture is uniquely reorganized during the oocyte-to-zygote transition in mice and is distinct in paternal and maternal nuclei within single-cell zygotes. Features of genomic organization including compartments, topologically associating domains (TADs) and loops are present in individual oocytes when averaged over the genome, but the presence of each feature at a locus varies between cells. At the sub-megabase level, we observed stochastic clusters of contacts that can occur across TAD boundaries but average into TADs. Notably, we found that TADs and loops, but not compartments, are present in zygotic maternal chromatin, suggesting that these are generated by different mechanisms. Our results demonstrate that the global chromatin organization of zygote nuclei is fundamentally different from that of other interphase cells. An understanding of this zygotic chromatin 'ground state' could potentially provide insights into reprogramming cells to a state of totipotency.

Recent advances enabled by the Hi-C technique have unraveled many principles of chromosomal folding that were subsequently linked to disease and gene regulation. In particular, Hi-C revealed that chromosomes of animals are organized into topologically associating domains (TADs), evolutionary conserved compact chromatin domains that influence gene expression. Mechanisms that underlie partitioning of the genome into TADs remain poorly understood. To explore principles of TAD folding in Drosophila melanogaster, we performed Hi-C and poly(A)+ RNA-seq in four cell lines of various origins (S2, Kc167, DmBG3-c2, and OSC). Contrary to previous studies, we find that regions between TADs (i.e., the inter-TADs and TAD boundaries) in Drosophila are only weakly enriched with the insulator protein dCTCF, while another insulator protein Su(Hw) is preferentially present within TADs. However, Drosophila inter-TADs harbor active chromatin and constitutively transcribed (housekeeping) genes. Accordingly, we find that binding of insulator proteins dCTCF and Su(Hw) predicts TAD boundaries much worse than active chromatin marks do. Interestingly, inter-TADs correspond to decompacted inter-bands of polytene chromosomes, whereas TADs mostly correspond to densely packed bands. Collectively, our results suggest that TADs are condensed chromatin domains depleted in active chromatin marks, separated by regions of active chromatin. We propose the mechanism of TAD self-assembly based on the ability of nucleosomes from inactive chromatin to aggregate, and lack of this ability in acetylated nucleosomal arrays. Finally, we test this hypothesis by polymer simulations and find that TAD partitioning may be explained by different modes of inter-nucleosomal interactions for active and inactive chromatin.

Recently, Chromosome Conformation Capture (3C) based experiments have highlighted the importance of computational models for the study of chromosome organization. In this review, we propose that current computational models can be grouped into roughly four classes, with two classes of data-driven models: consensus structures and data-driven ensembles, and two classes of de novo models: structural ensembles and mechanistic ensembles. Finally, we highlight specific questions mechanistic ensembles can address., National Institutes of Health (U.S.) (Grant R01HG003143), National Institutes of Health (U.S.) (Grant R01 GM114190)

Mitotic chromosomes are among the most recognizable structures in the cell, yet for over a century their internal organization remains largely unsolved. We applied chromosome conformation capture methods, 5C and Hi-C, across the cell cycle and revealed two alternative three-dimensional folding states of the human genome. We show that the highly compartmentalized and cell-type-specific organization described previously for non-synchronous cells is restricted to interphase. In metaphase, we identify a homogenous folding state, which is locus-independent, common to all chromosomes, and consistent among cell types, suggesting a general principle of metaphase chromosome organization. Using polymer simulations, we find that metaphase Hi-C data are inconsistent with classic hierarchical models, and is instead best described by a linearly-organized longitudinally compressed array of consecutive chromatin loops., National Cancer Institute (U.S.) (Grant U54CA143874)

Mammalian interphase chromosomes interact with the nuclear lamina (NL) through hundreds of large lamina-associated domains (LADs). We report a method to map NL contacts genome-wide in single human cells. Analysis of nearly 400 maps reveals a core architecture consisting of gene-poor LADs that contact the NL with high cell-to-cell consistency, interspersed by LADs with more variable NL interactions. The variable contacts tend to be cell-type specific and are more sensitive to changes in genome ploidy than the consistent contacts. Single-cell maps indicate that NL contacts involve multivalent interactions over hundreds of kilobases. Moreover, we observe extensive intra-chromosomal coordination of NL contacts, even over tens of megabases. Such coordinated loci exhibit preferential interactions as detected by Hi-C. Finally, the consistency of NL contacts is inversely linked to gene activity in single cells and correlates positively with the heterochromatic histone modification H3K9me3. These results highlight fundamental principles of single-cell chromatin organization., National Institutes of Health (U.S.) (Grant R01 GM114190), National Human Genome Research Institute (U.S.) (Grant R01 HG003143)

Topological constraints can affect both equilibrium and dynamic properties of polymer systems and can play a role in the organization of chromosomes. Despite many theoretical studies, the effects of topological constraints on the equilibrium state of a single compact polymer have not been systematically studied. Here we use simulations to address this longstanding problem. We find that sufficiently long unknotted polymers differ from knotted ones in the spatial and topological states of their subchains. The unknotted globule has subchains that are mostly unknotted and form asymptotically compact RG(s) ∼ s1/3 crumples. However, crumples display a high fractal dimension of the surface db = 2.8, forming excessive contacts and interpenetrating each other. We conclude that this topologically constrained equilibrium state resembles a conjectured crumpled globule [Grosberg et al., Journal de Physique, 1988, 49, 2095], but differs from its idealized hierarchy of self-similar, isolated and compact crumples., MIT-France Seed Fund, National Cancer Institute (U.S.). Physical Sciences-Oncology Center (MIT, (U54CA143874)), National Research University Higher School of Economics (Program for Basic Research)

Topologically associating domains (TADs) are fundamental structural and functional building blocks of human interphase chromosomes, yet the mechanisms of TAD formation remain unclear. Here, we propose that loop extrusion underlies TAD formation. In this process, cis-acting loop-extruding factors, likely cohesins, form progressively larger loops but stall at TAD boundaries due to interactions with boundary proteins, including CTCF. Using polymer simulations, we show that this model produces TADs and finer-scale features of Hi-C data. Each TAD emerges from multiple loops dynamically formed through extrusion, contrary to typical illustrations of single static loops. Loop extrusion both explains diverse experimental observations—including the preferential orientation of CTCF motifs, enrichments of architectural proteins at TAD boundaries, and boundary deletion experiments—and makes specific predictions for the depletion of CTCF versus cohesin. Finally, loop extrusion has potentially far-ranging consequences for processes such as enhancer-promoter interactions, orientation-specific chromosomal looping, and compaction of mitotic chromosomes., National Institutes of Health (U.S.) (grant R01 GM114190), National Institutes of Health (U.S.) (grant U54 DK107980), National Institutes of Health (U.S.) (grant R01 G003143), National Science Foundation (U.S.) (1504942)

The mechanism by which chromatids and chromosomes are segregated during mitosis and meiosis is a major puzzle of biology and biophysics. Using polymer simulations of chromosome dynamics, we show that a single mechanism of loop extrusion by condensins can robustly compact, segregate and disentangle chromosomes, arriving at individualized chromatids with morphology observed in vivo. Our model resolves the paradox of topological simplification concomitant with chromosome 'condensation', and explains how enzymes a few nanometers in size are able to control chromosome geometry and topology at micron length scales. We suggest that loop extrusion is a universal mechanism of genome folding that mediates functional interactions during interphase and compacts chromosomes during mitosis., National Science Foundation (U.S.) (Grant DMR-1206868), National Science Foundation (U.S.) (Grant MCB-1022117), National Institutes of Health (U.S.) (NIH Grant GM105847), National Institutes of Health (U.S.) (NIH Grant CA193419), National Institutes of Health (U.S.) (NIH grant GM114190), National Institutes of Health (U.S.) (NIH grant R01HG003143), National Institutes of Health (U.S.) (NIH grant DK107980)

Extracting biologically meaningful information from chromosomal interactions obtained with genome-wide chromosome conformation capture (3C) analyses requires the elimination of systematic biases. We present a computational pipeline that integrates a strategy to map sequencing reads with a data-driven method for iterative correction of biases, yielding genome-wide maps of relative contact probabilities. We validate this ICE (iterative correction and eigenvector decomposition) technique on published data obtained by the high-throughput 3C method Hi-C, and we demonstrate that eigenvector decomposition of the obtained maps provides insights into local chromatin states, global patterns of chromosomal interactions, and the conserved organization of human and mouse chromosomes., National Cancer Institute (U.S.). Physical Sciences-Oncology Center (U54CA143874)

Chromosomes must be highly compacted and organized within cells, but how this is achieved in vivo remains poorly understood. We report the use of chromosome conformation capture coupled with deep sequencing (Hi-C) to map the structure of bacterial chromosomes. Analysis of Hi-C data and polymer modeling indicates that the Caulobacter crescentus chromosome consists of multiple, largely independent spatial domains that are probably composed of supercoiled plectonemes arrayed into a bottle brush–like fiber. These domains are stable throughout the cell cycle and are reestablished concomitantly with DNA replication. We provide evidence that domain boundaries are established by highly expressed genes and the formation of plectoneme-free regions, whereas the histone-like protein HU and SMC (structural maintenance of chromosomes) promote short-range compaction and the colinearity of chromosomal arms, respectively. Collectively, our results reveal general principles for the organization and structure of chromosomes in vivo., National Institutes of Health (U.S.) (Grant R01GM082899), National Institutes of Health (U.S.) (Grant U54CA143874)

The three-dimensional folding of chromosomes compartmentalizes the genome and and can bring distant functional elements, such as promoters and enhancers, into close spatial proximity. Deciphering the relationship between chromosome organization and genome activity will aid in understanding genomic processes, like transcription and replication. However, little is known about how chromosomes fold. Microscopy is unable to distinguish large numbers of loci simultaneously or at high resolution. To date, the detection of chromosomal interactions using chromosome conformation capture (3C) and its subsequent adaptations required the choice of a set of target loci, making genome-wide studies impossible. We developed Hi-C, an extension of 3C that is capable of identifying long range interactions in an unbiased, genome-wide fashion. In Hi-C, cells are fixed with formaldehyde, causing interacting loci to be bound to one another by means of covalent DNA-protein cross-links. When the DNA is subsequently fragmented with a restriction enzyme, these loci remain linked. A biotinylated residue is incorporated as the 5' overhangs are filled in. Next, blunt-end ligation is performed under dilute conditions that favor ligation events between cross-linked DNA fragments. This results in a genome-wide library of ligation products, corresponding to pairs of fragments that were originally in close proximity to each other in the nucleus. Each ligation product is marked with biotin at the site of the junction. The library is sheared, and the junctions are pulled-down with streptavidin beads. The purified junctions can subsequently be analyzed using a high-throughput sequencer, resulting in a catalog of interacting fragments. Direct analysis of the resulting contact matrix reveals numerous features of genomic organization, such as the presence of chromosome territories and the preferential association of small gene-rich chromosomes. Correlation analysis can be applied to the contact matrix, demo

Chromosomes must be highly compacted and organized within cells, but how this is achieved in vivo remains poorly understood. We report the use of chromosome conformation capture coupled with deep sequencing (Hi-C) to map the structure of bacterial chromosomes. Analysis of Hi-C data and polymer modeling indicates that the Caulobacter crescentus chromosome consists of multiple, largely independent spatial domains that are probably composed of supercoiled plectonemes arrayed into a bottle brush–like fiber. These domains are stable throughout the cell cycle and are reestablished concomitantly with DNA replication. We provide evidence that domain boundaries are established by highly expressed genes and the formation of plectoneme-free regions, whereas the histone-like protein HU and SMC (structural maintenance of chromosomes) promote short-range compaction and the colinearity of chromosomal arms, respectively. Collectively, our results reveal general principles for the organization and structure of chromosomes in vivo., National Institutes of Health (U.S.) (Grant R01GM082899), National Institutes of Health (U.S.) (Grant U54CA143874)

The classic model of eukaryotic gene expression requires direct spatial contact between a distal enhancer and a proximal promoter. Recent Chromosome Conformation Capture (3C) studies show that enhancers and promoters are embedded in a complex network of looping interactions. Here we use a polymer model of chromatin fiber to investigate whether, and to what extent, looping interactions between elements in the vicinity of an enhancer-promoter pair can influence their contact frequency. Our equilibrium polymer simulations show that a chromatin loop, formed by elements flanking either an enhancer or a promoter, suppresses enhancer-promoter interactions, working as an insulator. A loop formed by elements located in the region between an enhancer and a promoter, on the contrary, facilitates their interactions. We find that different mechanisms underlie insulation and facilitation; insulation occurs due to steric exclusion by the loop, and is a global effect, while facilitation occurs due to an effective shortening of the enhancer-promoter genomic distance, and is a local effect. Consistently, we find that these effects manifest quite differently for in silico 3C and microscopy. Our results show that looping interactions that do not directly involve an enhancer-promoter pair can nevertheless significantly modulate their interactions. This phenomenon is analogous to allosteric regulation in proteins, where a conformational change triggered by binding of a regulatory molecule to one site affects the state of another site., National Cancer Institute (U.S.). Physical Sciences-Oncology Center (U54-CA143874-04), Massachusetts Institute of Technology. Undergraduate Research Opportunities Program, National Science Foundation (U.S.). Massachusetts Institute of Technology. Program for Research in Mathematics, Engineering and Science for High School Students (PRIMES)