Topological domains in mammalian genomes identified by analysis of chromatin interactions - PubMed (original) (raw)

Topological domains in mammalian genomes identified by analysis of chromatin interactions

Jesse R Dixon et al. Nature. 2012.

Abstract

The spatial organization of the genome is intimately linked to its biological function, yet our understanding of higher order genomic structure is coarse, fragmented and incomplete. In the nucleus of eukaryotic cells, interphase chromosomes occupy distinct chromosome territories, and numerous models have been proposed for how chromosomes fold within chromosome territories. These models, however, provide only few mechanistic details about the relationship between higher order chromatin structure and genome function. Recent advances in genomic technologies have led to rapid advances in the study of three-dimensional genome organization. In particular, Hi-C has been introduced as a method for identifying higher order chromatin interactions genome wide. Here we investigate the three-dimensional organization of the human and mouse genomes in embryonic stem cells and terminally differentiated cell types at unprecedented resolution. We identify large, megabase-sized local chromatin interaction domains, which we term 'topological domains', as a pervasive structural feature of the genome organization. These domains correlate with regions of the genome that constrain the spread of heterochromatin. The domains are stable across different cell types and highly conserved across species, indicating that topological domains are an inherent property of mammalian genomes. Finally, we find that the boundaries of topological domains are enriched for the insulator binding protein CTCF, housekeeping genes, transfer RNAs and short interspersed element (SINE) retrotransposons, indicating that these factors may have a role in establishing the topological domain structure of the genome.

PubMed Disclaimer

Conflict of interest statement

The authors declare no competing financial interests.

Figures

Figure 1

Figure 1. Topological Domains in the Mouse ES cell Genome

a, Normalized Hi-C interaction frequencies displayed as a 2D heatmap overlayed on ChIP-Seq data (from ref. 3), DI, HMM Bias State Calls, and domains. For both DI and HMM State calls, downstream bias (red) and upstream bias (green) are indicated. b, Schematic illustrating topological domains and resulting directional bias. c, Distribution of the DI (absolute value, in blue) compared to random (red). d, Mean interaction frequencies at all genomic distances between 40kb to 2Mb. Above 40kb, the intra-versus inter-domain interaction frequencies are significantly different (p < 0.005, Wilcoxan test). e, Box plot of all interaction frequencies at 80kb distance. Intra-domain interactions are enriched for high-frequency interactions. f–i, Diagram of “Intra-domain” (f) and “Inter-domain” FISH probes (g) and the genomic distance between pairs (h). i, Bar chart of the squared interprobe distance (from ref. 7) FISH probe pairs. Error bars indicate standard error (n = 100 for each probe pair).

Figure 2

Figure 2. Topological Boundaries Demonstrate Classical Insulator or Barrier Elements Features

a, 2D heatmap surrounding the HoxA locus and CS5 insulator in IMR90 cells. b, Enrichment of CTCF at boundary regions. c, The portion of CTCF binding sites that are considered “associated” with a boundary (within +/− 20kb window is used as the expected uncertainty due to 40kb binning). d, Heat maps of H3K9me3 at boundary sites in human and mouse. e, UCSC Genome Browser shot showing heterochromatin spreading in the human ES cells and IMR90 cells. The 2D heat map shows the interaction frequency in hES cells. f, Heat map of LADs (from ref. 15) surrounding the boundary regions.

Figure 3

Figure 3. Boundaries are shared across cell types and conserved in evolution

a, Overlap of boundaries between cell types. b, Genome browser shot of a cortex enriched dynamic interacting region that overlaps with the Foxg1 gene. c, Foxg1 expression in mouse ES cells and cortex as measured by RNA-seq. d, Heat map of the gene expression ratio between mouse ES cell and cortex of genes at dynamic interactions. e, Pie chart of inter- and intra-domain dynamic interactions. f, Overlap of boundaries between syntenic mouse and human sequences (p-value < 2.2*10−16 compared to random, Fisher’s exact test). g and h, Genome browser shots showing domain structure over a syntenic region in the mouse (g) and human ES cells (h). Note: the region in humans has been inverted from its normal UCSC coordinates for proper display purposes.

Figure 4

Figure 4. Boundary regions are enriched for housekeeping genes

a, Chromatin modifications, TSS, GRO-Seq, and SINE elements surrounding boundary regions in mESCs or IMR90. b, Boundaries associated with a CTCF binding site, housekeeping gene, or tRNA gene (purple) compared to expected at random (grey). c, Gene Ontology p-value chart. d, Enrichment of housekeeping genes (gold) and tissue specific genes (blue) as defined by on Shannon entropy scores near boundaries normalized for the number of genes in each class (TSS/10kb/total TSS). e, Percentage of boundaries with a given mark within 20kb of the boundaries.

Similar articles

Cited by

References

    1. Cremer T, Cremer M. Chromosome territories. Cold Spring Harb Perspect Biol. 2:a003889. - PMC - PubMed
    1. Lieberman-Aiden E, et al. Comprehensive mapping of long-range interactions reveals folding principles of the human genome. Science. 2009;326:289–93. - PMC - PubMed
    1. Shen Y, et al. A Map of cis-Regulatory Sequences in the Mouse Genome. 2012 in submission. - PMC - PubMed
    1. Yaffe E, Tanay A. Probabilistic modeling of Hi-C contact maps eliminates systematic biases to characterize global chromosomal architecture. Nat Genet. 43:1059–65. - PubMed
    1. Wang KC, et al. A long noncoding RNA maintains active chromatin to coordinate homeotic gene expression. Nature. 472:120–4. - PMC - PubMed

Publication types

MeSH terms

Substances

Grants and funding

LinkOut - more resources