NucleoMap: A computational tool for identifying nucleosomes in ultra-high resolution contact maps

Although poorly positioned nucleosomes are ubiquitous in the eukaryotic genome, they are difficult to identify with existing nucleosome identification methods. Recently available enhanced high-throughput chromatin conformation capture techniques such as Micro-C, DNase Hi-C, and Hi-CO characterize nucleosome-level chromatin proximity, probing the positions of mono-nucleosomes and the spacing between nucleosome pairs at the same time, enabling nucleosome profiling in poorly positioned regions. Here we develop a novel computational approach, NucleoMap, to identify nucleosome positioning from ultra-high resolution chromatin contact maps. By integrating nucleosome read density, contact distances, and binding preferences, NucleoMap precisely locates nucleosomes in both prokaryotic and eukaryotic genomes and outperforms existing nucleosome identification methods in both precision and recall. We rigorously characterize genome-wide association in eukaryotes between the spatial organization of mono-nucleosomes and their corresponding histone modifications, protein binding activities, and higher-order chromatin functions. We also find evidence of two tetra-nucleosome folding structures in human embryonic stem cells and analyze their association with multiple structural and functional regions. Based on the identified nucleosomes, nucleosome contact maps are constructed, reflecting the inter-nucleosome distances and preserving the contact distance profiles in original contact maps.     

Divergence of Mammalian Higher Order Chromatin Structure Is Associated with Developmental Loci

Several recent studies have examined different aspects of mammalian higher order chromatin structure – replication timing, lamina association and Hi-C inter-locus interactions — and have suggested that most of these features of genome organisation are conserved over evolution. However, the extent of evolutionary divergence in higher order structure has not been rigorously measured across the mammalian genome, and until now little has been known about the characteristics of any divergent loci present. Here, we generate a dataset combining multiple measurements of chromatin structure and organisation over many embryonic cell types for both human and mouse that, for the first time, allows a comprehensive assessment of the extent of structural divergence between mammalian genomes. Comparison of orthologous regions confirms that all measurable facets of higher order structure are conserved between human and mouse, across the vast majority of the detectably orthologous genome. This broad similarity is observed in spite of many loci possessing cell type specific structures. However, we also identify hundreds of regions (from 100 Kb to 2.7 Mb in size) showing consistent evidence of divergence between these species, constituting at least 10% of the orthologous mammalian genome and encompassing many hundreds of human and mouse genes. These regions show unusual shifts in human GC content, are unevenly distributed across both genomes, and are enriched in human subtelomeric regions. Divergent regions are also relatively enriched for genes showing divergent expression patterns between human and mouse ES cells, implying these regions cause divergent regulation. Particular divergent loci are strikingly enriched in genes implicated in vertebrate development, suggesting important roles for structural divergence in the evolution of mammalian developmental programmes. These data suggest that, though relatively rare in the mammalian genome, divergence in higher order chromatin structure has played important roles during evolution.     

Dynamic distribution of linker histone h1.5 in cellular differentiation

Linker histones are essential components of chromatin, but the distributions and functions of many during cellular differentiation are not well understood. Here, we show that H1.5 binds to genic and intergenic regions, forming blocks of enrichment, in differentiated human cells from all three embryonic germ layers but not in embryonic stem cells. In differentiated cells, H1.5, but not H1.3, binds preferentially to genes that encode membrane and membrane-related proteins. Strikingly, 37% of H1.5 target genes belong to gene family clusters, groups of homologous genes that are located in proximity to each other on chromosomes. H1.5 binding is associated with gene repression and is required for SIRT1 binding, H3K9me2 enrichment, and chromatin compaction. Depletion of H1.5 results in loss of SIRT1 and H3K9me2, increased chromatin accessibility, deregulation of gene expression, and decreased cell growth. Our data reveal for the first time a specific and novel function for linker histone subtype H1.5 in maintenance of condensed chromatin at defined gene families in differentiated human cells.     

Comparison of the three-dimensional organization of sperm and fibroblast genomes using the Hi-C approach

BackgroundThe three-dimensional organization of the genome is tightly connected to its biological function. The Hi-C approach was recently introduced as a method that can be used to identify higher-order chromatin interactions genome-wide. The aim of this study was to determine genome-wide chromatin interaction frequencies using the Hi-C approach in mouse sperm cells and embryonic fibroblasts.ResultsThe obtained data demonstrate that the three-dimensional genome organizations of sperm and fibroblast cells show a high degree of similarity both with each other and with the previously described mouse embryonic stem cells. Both A- and B-compartments and topologically associated domains are present in spermatozoa and fibroblasts. Nevertheless, sperm cells and fibroblasts exhibit statistically significant differences between each other in the contact probabilities of defined loci. Tight packaging of the sperm genome results in an enrichment of long-range contacts compared with the fibroblasts. However, only 30% of the differences in the number of contacts are based on differences in the densities of their genome packages; the main source of the differences is the gain or loss of contacts that are specific for defined genome regions. We find that the dependence of the contact probability on genomic distance for sperm is close to the dependence predicted for the fractal globular folding of chromatin.ConclusionsOverall, we can conclude that the three-dimensional structure of the genome is passed through generations without being dramatically changed in sperm cells.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0642-0) contains supplementary material, which is available to authorized users.     

Manipulation of human minichromosomes to carry greater than megabase-sized chromosome inserts

For introducing regions of human chromosomes greater than a megabase into cells or animals, we have developed a chromosome-cloning system in which defined regions of human chromosomes can be cloned into a stable human minichromosome vector in homologous recombination-proficient chicken DT40 cells. The stable minichromosome vector allowed a 10 Mb-sized region of the mitotically unstable human chromosome 22 to be stably maintained in mouse embryonic stem (ES) cells, and in mice. Furthermore, we demonstrated functional expression of human genes from the HAC in mice. This study describes a stable cloning and expression system for greater than megabase-sized regions of human chromosomes.     

Evidence for Local Regulatory Control of Escape from Imprinted X Chromosome Inactivation

X chromosome inactivation (XCI) is an epigenetic process that almost completely inactivates one of two X chromosomes in somatic cells of mammalian females. A few genes are known to escape XCI and the mechanism for this escape remains unclear. Here, using mouse trophoblast stem (TS) cells, we address whether particular chromosomal interactions facilitate escape from imprinted XCI. We demonstrate that promoters of genes escaping XCI do not congregate to any particular region of the genome in TS cells. Further, the escape status of a gene was uncorrelated with the types of genomic features and gene activity located in contacted regions. Our results suggest that genes escaping imprinted XCI do so by using the same regulatory sequences as their expressed alleles on the active X chromosome. We suggest a model where regulatory control of escape from imprinted XCI is mediated by genomic elements located in close linear proximity to escaping genes.     

Comparative 3D Genome Structure Analysis of the Fission and the Budding Yeast

We studied the 3D structural organization of the fission yeast genome, which emerges from the tethering of heterochromatic regions in otherwise randomly configured chromosomes represented as flexible polymer chains in an nuclear environment. This model is sufficient to explain in a statistical manner many experimentally determined distinctive features of the fission yeast genome, including chromatin interaction patterns from Hi-C experiments and the co-locations of functionally related and co-expressed genes, such as genes expressed by Pol-III. Our findings demonstrate that some previously described structure-function correlations can be explained as a consequence of random chromatin collisions driven by a few geometric constraints (mainly due to centromere-SPB and telomere-NE tethering) combined with the specific gene locations in the chromosome sequence. We also performed a comparative analysis between the fission and budding yeast genome structures, for which we previously detected a similar organizing principle. However, due to the different chromosome sizes and numbers, substantial differences are observed in the 3D structural genome organization between the two species, most notably in the nuclear locations of orthologous genes, and the extent of nuclear territories for genes and chromosomes. However, despite those differences, remarkably, functional similarities are maintained, which is evident when comparing spatial clustering of functionally related genes in both yeasts. Functionally related genes show a similar spatial clustering behavior in both yeasts, even though their nuclear locations are largely different between the yeast species.     

Large-scale reconstruction of 3D structures of human chromosomes from chromosomal contact data

Chromosomes are not positioned randomly within a nucleus, but instead, they adopt preferred spatial conformations to facilitate necessary long-range gene–gene interactions and regulations. Thus, obtaining the 3D shape of chromosomes of a genome is critical for understanding how the genome folds, functions and how its genes interact and are regulated. Here, we describe a method to reconstruct preferred 3D structures of individual chromosomes of the human genome from chromosomal contact data generated by the Hi-C chromosome conformation capturing technique. A novel parameterized objective function was designed for modeling chromosome structures, which was optimized by a gradient descent method to generate chromosomal structural models that could satisfy as many intra-chromosomal contacts as possible. We applied the objective function and the corresponding optimization method to two Hi-C chromosomal data sets of both a healthy and a cancerous human B-cell to construct 3D models of individual chromosomes at resolutions of 1 MB and 200 KB, respectively. The parameters used with the method were calibrated according to an independent fluorescence in situ hybridization experimental data. The structural models generated by our method could satisfy a high percentage of contacts (pairs of loci in interaction) and non-contacts (pairs of loci not in interaction) and were compatible with the known two-compartment organization of human chromatin structures. Furthermore, structural models generated at different resolutions and from randomly permuted data sets were consistent.     

Ten families of variant genes encoded in subtelomeric regions of multiple chromosomes of Plasmodium chabaudi,a malaria species that undergoes antigenic variation in the laboratory mouse

The chromosome ends of human malaria parasites harbour many genes encoding proteins that are exported to the surface of infected red cells, often being involved in host-parasite interactions and immune evasion. Unlike other murine malaria parasites Plasmodium chabaudi undergoes antigenic variation during passage in the laboratory mouse and hence is a model suitable for investigation of switching mechanisms. However, little is known about the subtelomeric regions of P. chabaudi chromosomes and its variable antigens. Here we report 80 kb of sequence from an end of one P. chabaudi chromosome. Hybridization of probes spanning this region to two dimensional pulsed field gels of the genome revealed 10 multicopy gene families located exclusively in subtelomeric regions of multiple P. chabaudi chromosomes, interspersed amongst multicopy intergenic regions. Hence all chromosomes share a common subtelomeric structure, presumably playing a similar role in spatial positioning as the P. falciparum Rep20 sequence. Expression in blood stages, domains characteristic of surface antigens and copy numbers between four and several hundred per genome, indicate a functional role in antigenic variation for some of these families. We identify members of the cir family, as well as novel genes, that although clearly homologous to cir have large low complexity regions in the predicted extracellular domains. Although all families have homologues in other rodent Plasmodium species, four were previously not known to be subtelomeric. Six have homologues in human and simian malarias.     

Common themes and cell type specific variations of higher order chromatin arrangements in the mouse

Background  

Similarities as well as differences in higher order chromatin arrangements of human cell types were previously reported. For an evolutionary comparison, we now studied the arrangements of chromosome territories and centromere regions in six mouse cell types (lymphocytes, embryonic stem cells, macrophages, fibroblasts, myoblasts and myotubes) with fluorescence in situ hybridization and confocal laser scanning microscopy. Both species evolved pronounced differences in karyotypes after their last common ancestors lived about 87 million years ago and thus seem particularly suited to elucidate common and cell type specific themes of higher order chromatin arrangements in mammals.     

Relaxed 3D genome conformation facilitates the pluripotent to totipotent-like state transition in embryonic stem cells

The 3D genome organization is crucial for gene regulation. Although recent studies have revealed a uniquely relaxed genome conformation in totipotent early blastomeres of both fertilized and cloned embryos, how weakened higher-order chromatin structure is functionally linked to totipotency acquisition remains elusive. Using low-input Hi-C, ATAC-seq and ChIP-seq, we systematically examined the dynamics of 3D genome and epigenome during pluripotent to totipotent-like state transition in mouse embryonic stem cells (ESCs). The spontaneously converted 2-cell-embryo-like cells (2CLCs) exhibited more relaxed chromatin architecture compared to ESCs, including global weakening of both enhancer-promoter interactions and TAD insulation. While the former correlated with inactivation of ESC enhancers and down-regulation of pluripotent genes, the latter might facilitate contacts between the putative new enhancers arising in 2CLCs and neighboring 2C genes. Importantly, disruption of chromatin organization by depleting CTCF or the cohesin complex promoted the ESC to 2CLC transition. Our results thus establish a critical role of 3D genome organization in totipotency acquisition.     

Innate structure of DNA foci restricts the mixing of DNA from different chromosome territories

The distribution of chromatin within the mammalian nucleus is constrained by its organization into chromosome territories (CTs). However, recent studies have suggested that promiscuous intra- and inter-chromosomal interactions play fundamental roles in regulating chromatin function and so might define the spatial integrity of CTs. In order to test the extent of DNA mixing between CTs, DNA foci of individual CTs were labeled in living cells following incorporation of Alexa-488 and Cy-3 conjugated replication precursor analogues during consecutive cell cycles. Uniquely labeled chromatin domains, resolved following random mitotic segregation, were visualized as discrete structures with defined borders. At the level of resolution analysed, evidence for mixing of chromatin from adjacent domains was only apparent within the surface volumes where neighboring CTs touched. However, while less than 1% of the nuclear volume represented domains of inter-chromosomal mixing, the dynamic plasticity of DNA foci within individual CTs allows continual transformation of CT structure so that different domains of chromatin mixing evolve over time. Notably, chromatin mixing at the boundaries of adjacent CTs had little impact on the innate structural properties of DNA foci. However, when TSA was used to alter the extent of histone acetylation changes in chromatin correlated with increased chromatin mixing. We propose that DNA foci maintain a structural integrity that restricts widespread mixing of DNA and discuss how the potential to dynamically remodel genome organization might alter during cell differentiation.     

Three-dimensional analysis of nuclear heterochromatin distribution during early development in the rabbit

Changes to the spatial organization of specific chromatin domains such as constitutive heterochromatin have been studied extensively in somatic cells. During early embryonic development, drastic epigenetic reprogramming of both the maternal and paternal genomes, followed by chromatin remodeling at the time of embryonic genome activation (EGA), have been observed in the mouse. Very few studies have been performed in other mammalian species (human, bovine, or rabbit) and the data are far from complete. During this work, we studied the three-dimensional organization of pericentromeric regions during the preimplantation period in the rabbit using specific techniques (3D-FISH) and tools (semi-automated image analysis). We observed that the pericentromeric regions (identified with specific probes for Rsat I and Rsat II genomic sequences) changed their shapes (from pearl necklaces to clusters), their nuclear localizations (from central to peripheral), as from the 4-cell stage. This reorganization goes along with histone modification changes and reduced amount of interactions with nucleolar precursor body surface. Altogether, our results suggest that the 4-cell stage may be a crucial window for events necessary before major EGA, which occurs during the 8-cell stage in the rabbit.