Neva C. Durand

Chromatin extrusion explains key features of loop and domain formation in wild-type and engineered genomes

Significance When the human genome folds up inside the cell nucleus, it is spatially partitioned into numerous loops and contact domains. How these structures form is unknown. Here, we show that data from high-resolution spatial proximity maps are consistent with a model in which a complex, including the proteins CCCTC-binding factor (CTCF) and cohesin, mediates the formation of loops by a process of extrusion. Contact domains form as a byproduct of this process. The model accurately predicts how the genome will fold, using only information about the locations at which CTCF is bound. We demonstrate the ability to reengineer loops and domains in a predictable manner by creating highly targeted mutations, some as small as a single base pair, at CTCF sites. We recently used in situ Hi-C to create kilobase-resolution 3D maps of mammalian genomes. Here, we combine these maps with new Hi-C, microscopy, and genome-editing experiments to study the physical structure of chromatin fibers, domains, and loops. We find that the observed contact domains are inconsistent with the equilibrium state for an ordinary condensed polymer. Combining Hi-C data and novel mathematical theorems, we show that contact domains are also not consistent with a fractal globule. Instead, we use physical simulations to study two models of genome folding. In one, intermonomer attraction during polymer condensation leads to formation of an anisotropic “tension globule.” In the other, CCCTC-binding factor (CTCF) and cohesin act together to extrude unknotted loops during interphase. Both models are consistent with the observed contact domains and with the observation that contact domains tend to form inside loops. However, the extrusion model explains a far wider array of observations, such as why loops tend not to overlap and why the CTCF-binding motifs at pairs of loop anchors lie in the convergent orientation. Finally, we perform 13 genome-editing experiments examining the effect of altering CTCF-binding sites on chromatin folding. The convergent rule correctly predicts the affected loops in every case. Moreover, the extrusion model accurately predicts in silico the 3D maps resulting from each experiment using only the location of CTCF-binding sites in the WT. Thus, we show that it is possible to disrupt, restore, and move loops and domains using targeted mutations as small as a single base pair.

Juicer Provides a One-Click System for Analyzing Loop-Resolution Hi-C Experiments

Hi-C experiments explore the 3D structure of the genome, generating terabases of data to create high-resolution contact maps. Here, we introduce Juicer, an open-source tool for analyzing terabase-scale Hi-C datasets. Juicer allows users without a computational background to transform raw sequence data into normalized contact maps with one click. Juicer produces a hic file containing compressed contact matrices at many resolutions, facilitating visualization and analysis at multiple scales. Structural features, such as loops and domains, are automatically annotated. Juicer is available as open source software at http://aidenlab.org/juicer/.

Deletion of DXZ4 on the human inactive X chromosome alters higher-order genome architecture.

Significance In human females, one of the two X chromosomes is inactive (Xi) and adopts an unusual 3D conformation. The Xi chromosome contains superloops, large chromatin loops that are often anchored at the macrosatellite repeat DXZ4, and is partitioned into two large intervals, called superdomains, whose boundary lies at DXZ4. Here, we use spatial proximity mapping, microscopy, and genome editing to study the Xi. We find that superloops and superdomains are conserved across humans, macaque, and mouse. By mapping proximity between three or more loci, we show that superloops tend to occur simultaneously. Deletion of DXZ4 from the human Xi disrupts superloops, eliminates superdomains, and alters chromatin modifications. Finally, we show that a model in which CCCTC-binding factor (CTCF) and cohesin extrude chromatin can explain the formation of superloops and superdomains. During interphase, the inactive X chromosome (Xi) is largely transcriptionally silent and adopts an unusual 3D configuration known as the “Barr body.” Despite the importance of X chromosome inactivation, little is known about this 3D conformation. We recently showed that in humans the Xi chromosome exhibits three structural features, two of which are not shared by other chromosomes. First, like the chromosomes of many species, Xi forms compartments. Second, Xi is partitioned into two huge intervals, called “superdomains,” such that pairs of loci in the same superdomain tend to colocalize. The boundary between the superdomains lies near DXZ4, a macrosatellite repeat whose Xi allele extensively binds the protein CCCTC-binding factor. Third, Xi exhibits extremely large loops, up to 77 megabases long, called “superloops.” DXZ4 lies at the anchor of several superloops. Here, we combine 3D mapping, microscopy, and genome editing to study the structure of Xi, focusing on the role of DXZ4. We show that superloops and superdomains are conserved across eutherian mammals. By analyzing ligation events involving three or more loci, we demonstrate that DXZ4 and other superloop anchors tend to colocate simultaneously. Finally, we show that deleting DXZ4 on Xi leads to the disappearance of superdomains and superloops, changes in compartmentalization patterns, and changes in the distribution of chromatin marks. Thus, DXZ4 is essential for proper Xi packaging.

De novo assembly of the Aedes aegypti genome using Hi-C yields chromosome-length scaffolds

Hi-C for mosquito genomes Most genomes sequenced today are determined through the generation of short sequenced bits of DNA that are computationally pieced together like a jigsaw puzzle. This has resulted in the need for funds and additional data to fill in gaps in order to fully assemble the many chromosomes that make up a eukaryotic genome. Dudchenko et al. used the Hi-C method, which measures the distance between contact points within and between chromosomes for scaffold validation, together with correction and ordering to more completely determine the arrangement of short sequencing reads for genome mapping. They validated their approach through the de novo generation of a complete human genome. A comparative analysis of mosquito genomes was made possible by improving the Culex quinquefasciatus genome assembly and generating the genome of Aedes aegypti, the vector of Zika virus. Science, this issue p. 92 The DNA proximity ligation method Hi-C was used to create a genome assembly for the mosquito Aedes aegypti. The Zika outbreak, spread by the Aedes aegypti mosquito, highlights the need to create high-quality assemblies of large genomes in a rapid and cost-effective way. Here we combine Hi-C data with existing draft assemblies to generate chromosome-length scaffolds. We validate this method by assembling a human genome, de novo, from short reads alone (67× coverage). We then combine our method with draft sequences to create genome assemblies of the mosquito disease vectors Ae. aegypti and Culex quinquefasciatus, each consisting of three scaffolds corresponding to the three chromosomes in each species. These assemblies indicate that almost all genomic rearrangements among these species occur within, rather than between, chromosome arms. The genome assembly procedure we describe is fast, inexpensive, and accurate, and can be applied to many species.We present an end-to-end genome assembly of a female Aedes aegypti mosquito, which spreads viral diseases such as yellow fever, dengue, chikungunya, and Zika to humans. The assembly is based on an earlier genome published in 2007 and improved in 2013. The new assembly has a scaffold N50 of 419Mb, with 96.9% of the ungapped sequence anchored to chromosomes. We used the new assembly to examine the conservation of A. aegypti chromosomes. Our results suggest that synteny is strongly conserved between Ae. aegypti and An. gambiae. Comparison to D. melanogaster highlights the extent to which the identity of entire chromosome arms is preserved across dipterans. Main text: Due to its role in the spread of the Zika virus in the Americas, A. aegypti – an important mosquito vector of many human diseases – is causing a new wave of widespread concern (Harmon 2016). The lack of an end-to-end genome assembly is limiting our understanding of the biology of this major arbovirus vector and hinders efforts at disease control. To aid in the response, we present an improved assembly. Our assembly is based on AaegL2 (Nene et al. 2007), which was generated using Sanger reads (8X coverage) assembled using ARACHNE (Jaffe et al. 2003). The AaegL2 assembly consists of 4756 scaffolds spanning 1.3Gb of sequence, with a contig N50 of 82 Kb and scaffold N50 of 1.5Mb. Our effort to improve AaegL2 resulted in the first end-to-end genome assembly of A. aegypti. We refer to our new assembly as AaegL4. (AaegL3 is frequently used to refer to a variant of AaegL2 which includes the mitochondrial genome.) In addition to providing an end-to-end assembly, AaegL4 also improves the quality of anchoring and scaffolding. The AaegL4 assembly is shared at https://tinyurl.com/AaegL4. Given the public health relevance of this work, we choose to share this assembly with the scientific community without delay. A description of the methodology used to generate AaegL4 is currently in preparation.

Juicebox.js Provides a Cloud-Based Visualization System for Hi-C Data

SUMMARY Contact mapping experiments such as Hi-C explore how genomes fold in 3D. Here, we introduce Juicebox.js, a cloud-based web application for exploring the resulting datasets. Like the original Juicebox application, Juicebox.js allows users to zoom in and out of such datasets using an interface similar to Google Earth. Juicebox.js also has many features designed to facilitate data reproducibility and sharing. Furthermore, Juicebox.js encodes the exact state of the browser in a shareable URL. Creating a public browser for a new Hi-C dataset does not require coding and can be accomplished in under a minute. The web app also makes it possible to create interactive figures online that can complement or replace ordinary journal figures. When combined with Juicer, this makes the entire process of data analysis transparent, insofar as every step from raw reads to published figure is publicly available as open source code.

Genetic determinants of chromatin accessibility and gene regulation in T cell activation across human individuals

The vast majority of genetic variants associated with complex human traits map to non-coding regions, but little is understood about how they modulate gene regulation in health and disease. Here, we analyzed Assay for Transposase-Accessible Chromatin (ATAC-seq) profiles from activated primary CD4+ T cells of 105 healthy donors to identify ATAC-QTLs: genetic variants that affect chromatin accessibility. We found that ATAC-QTLs are widespread, disrupt binding sites for transcription factors known to be important for CD4+ T cell differentiation and activation, overlap and mediate expression QTLs from the same cells and are enriched for SNPs associated with autoimmune diseases. We also identified numerous pairs of ATAC-peaks with highly correlated chromatin accessibility. When we characterize 3D chromosome organization in primary CD4+ T cells by in situ-Hi-C, we found that correlated peaks tend to reside in the same chromatin contact domains, span super-enhancers, and are more impacted by ATAC-QTLs than single peaks. Thus, variability in chromatin accessibility in primary CD4+ T cells is heritable, is determined by genetic variation in a manner affected by the 3D organization of the genome, and mediates genetic effects on gene expression. Our results provide insights into how genetic variants modulate chromatin state and gene expression in primary immune cells that play a key role in many human diseases.

Hybrid de novo genome assembly and centromere characterization of the gray mouse lemur (Microcebus murinus)

BackgroundThe de novo assembly of repeat-rich mammalian genomes using only high-throughput short read sequencing data typically results in highly fragmented genome assemblies that limit downstream applications. Here, we present an iterative approach to hybrid de novo genome assembly that incorporates datasets stemming from multiple genomic technologies and methods. We used this approach to improve the gray mouse lemur (Microcebus murinus) genome from early draft status to a near chromosome-scale assembly.MethodsWe used a combination of advanced genomic technologies to iteratively resolve conflicts and super-scaffold the M. murinus genome.ResultsWe improved the M. murinus genome assembly to a scaffold N50 of 93.32 Mb. Whole genome alignments between our primary super-scaffolds and 23 human chromosomes revealed patterns that are congruent with historical comparative cytogenetic data, thus demonstrating the accuracy of our de novo scaffolding approach and allowing assignment of scaffolds to M. murinus chromosomes. Moreover, we utilized our independent datasets to discover and characterize sequences associated with centromeres across the mouse lemur genome. Quality assessment of the final assembly found 96% of mouse lemur canonical transcripts nearly complete, comparable to other published high-quality reference genome assemblies.ConclusionsWe describe a new assembly of the gray mouse lemur (Microcebus murinus) genome with chromosome-scale scaffolds produced using a hybrid bioinformatic and sequencing approach. The approach is cost effective and produces superior results based on metrics of contiguity and completeness. Our results show that emerging genomic technologies can be used in combination to characterize centromeres of non-model species and to produce accurate de novo chromosome-scale genome assemblies of complex mammalian genomes.

The Juicebox Assembly Tools module facilitates de novo assembly of mammalian genomes with chromosome-length scaffolds for under

Hi-C contact maps are valuable for genome assembly (Lieberman-Aiden, van Berkum et al. 2009; Burton et al. 2013; Dudchenko et al. 2017). Recently, we developed Juicebox, a system for the visual exploration of Hi-C data (Durand, Robinson et al. 2016), and 3D-DNA, an automated pipeline for using Hi-C data to assemble genomes (Dudchenko et al. 2017). Here, we introduce “Assembly Tools,” a new module for Juicebox, which provides a point-and-click interface for using Hi-C heatmaps to identify and correct errors in a genome assembly. Together, 3D-DNA and the Juicebox Assembly Tools greatly reduce the cost of accurately assembling complex eukaryotic genomes. To illustrate, we generated de novo assemblies with chromosome-length scaffolds for three mammals: the wombat, Vombatus ursinus (3.3Gb), the Virginia opossum, Didelphis virginiana (3.3Gb), and the raccoon, Procyon lotor (2.5Gb). The only inputs for each assembly were Illumina reads from a short insert DNA-Seq library (300 million Illumina reads, maximum length 2x150 bases) and an in situ Hi-C library (100 million Illumina reads, maximum read length 2x150 bases), which cost <

1000

1000.

CTCF looping is established during gastrulation in medaka embryos

Genome architecture plays a critical role in gene regulation, but how the structures seen in mature cells emerge during embryonic development remains poorly understood. Here, we study early development in medaka (the Japanese killifish, Oryzias latipes) at 12 time points before, during, and after gastrulation which is the most dramatic event in early embryogenesis, and characterize transcription, protein binding, and genome architecture. We find that gastrulation is most associated with drastic changes in genome architecture, including the formation of the first loops between sites bound by the insulator protein CTCF and great increase in the size of contact domains. However, the position of CTCF is fixed throughout medaka embryogenesis. Interestingly, genome-wide transcription precedes the emergence of mature domains and CTCF-CTCF loops.