GAMIBHEAR: whole-genome haplotype reconstruction from genome architecture mapping data

Motivation: Understanding haplotype-specific regulatory mechanisms becomes increasingly important in genomics and medical research. Investigating differences in allele-specific gene expression, epigenetic changes and their causal variants greatly benefits from haplotype reconstruction or phasing of genetic variants, but direct evidence for the haplotype structure is difficult to obtain from standard short-read sequencing data. Chromatin conformation data obtained from 3C experiments allows inference of haplotypes because inter-chromosomal contacts are more frequent than homologous intra-chromosomal contacts, but these data suffer from technical biases owing to the digestion and ligation process of the 3C technique. Genome Architecture Mapping (GAM) is a novel digestion- and ligation-free method for the inference of chromatin conformation from nuclear cryosections. Due to its high resolution and independence of enzymatic digestion it is well-suited for haplotype reconstruction and for detecting haplotype-specific chromatin contacts. Results: Here, we present GAMIBHEAR, a tool for accurate haplotype reconstruction from GAM data. GAMIBHEAR aggregates allelic co-observation frequencies across multiple nuclear slices and employs a GAM-specific probabilistic model of haplotype capture to optimise phasing accuracy. Using a hybrid mouse embryonic stem cell line with known haplotype structure as a benchmark dataset, we assess correctness and completeness of the reconstructed haplotypes, and demonstrate the power of GAM data and the accuracy of GAMIBHEAR to infer genome-wide haplotypes. Availability: GAMIBHEAR is available as an R package under the open source GPL-2 license at https://bitbucket.org/schwarzlab/gamibhear Maintainer: julia.markowski{at}mdc-berlin.d

Comparison of the Hi-C, GAM and SPRITE methods by use of polymer models of chromatin

Powerful technologies have been developed to probe chromatin 3D physical interactions genome-wide, such as Hi-C, GAM and SPRITE. Due to their intrinsic differences and without a benchmarking reference, it is currently difficult to assess how well each method represents the genome 3D structure and their relative performance. Here, we develop a computational approach to implement Hi-C, GAM and SPRITE in-silico to compare the three methods in a simplified, yet controlled framework against known polymer 3D structures. We test our approach on models of three 6-Mb genomic regions, around the Sox9 and the HoxD genes in mouse ES cells, and around the Epha4 gene in mouse CHLX-12 cells. The model-derived contact matrices consistently match Hi-C, GAM and SPRITE experiments. We show that in-silico Hi-C, GAM and SPRITE average data are overall faithful to the 3D structures of the polymer models. We find that the inherent variability of model single-molecule 3D conformations and experimental efficiency differently affect the contact data of the different methods. Similarly, the noise-to-signal levels vary with genomic distance differently in in-silico Hi-C, SPRITE and GAM. We benchmark the performance of each technology in bulk and in single-cell experiments, and identify the minimal number of cells required for replicates to return statistically consistent chromatin contact measures. Under the same experimental conditions, SPRITE requires the lowest number of cells, Hi-C is close to SPRITE, while GAM is the most reproducible method to capture interactions at large genomic distances

Single-cell-resolved dynamics of chromatin architecture delineate cell and regulatory states in wildtype and cloche/npas4l mutant zebrafish embryos

DNA accessibility of cis regulatory elements (CREs) dictates transcriptional activity and drives cell differentiation during development. While many of the genes that regulate embryonic development have been described, the underlying CRE dynamics controlling their expression remain largely unknown. To address this, we applied single-cell combinatorial indexing ATAC-seq (sci-ATAC-seq) to whole 24 hours post fertilization (hpf) stage zebrafish embryos and developed a new computational tool, ScregSeg, that selects informative genome segments and classifies complex accessibility dynamics. We integrated the ScregSeg output with bulk measurements for histone post-translational modifications and 3D genome organization, expanding knowledge of regulatory principles between chromatin modalities. Sci-ATAC-seq profiling of npas4l/cloche mutant embryos revealed novel cellular roles for this hemato-vascular transcriptional master regulator and suggests an intricate mechanism regulating its expression. Our work constitutes a valuable resource for future studies in developmental, molecular, and computational biology

Investigating field-induced magnetic order in Han purple by neutron scattering up to 25.9 T

BaCuSi2O6 is a quasi-two-dimensional (2D) quantum antiferromagnet containing three different types of stacked, square-lattice bilayer hosting spin-1/2 dimers. Although this compound has been studied extensively over the last two decades, the critical applied magnetic field required to close the dimer spin gap and induce magnetic order, which exceeds 23 T, has to date precluded any kind of neutron scattering investigation. However, the HFM/EXED instrument at the Helmholtz-Zentrum Berlin made this possible at magnetic fields up to 25.9 T. Thus we have used HFM/EXED to investigate the field-induced ordered phase, in particular to look for quasi-2D physics arising from the layered structure and from the different bilayer types. From neutron diffraction data, we determined the global dependence of the magnetic order parameter on both magnetic field and temperature, finding a form consistent with 3D quantum critical scaling; from this we deduce that the quasi-2D interactions and nonuniform layering of BaCuSi2O6 are not anisotropic enough to induce hallmarks of 2D physics. From neutron spectroscopy data, we measured the dispersion of the strongly Zeeman-split magnetic excitations, finding good agreement with the zero-field interaction parameters of BaCuSi2O6. We conclude that HFM/EXED allowed a significant extension in the application of neutron scattering techniques to the field range above 20 T and in particular opened previously unavailable possibilities in the study of field-induced magnetic quantum phase transitions

Comparison of the Hi-C, GAM and SPRITE methods using polymer models of chromatin

Hi-C, split-pool recognition of interactions by tag extension (SPRITE) and genome architecture mapping (GAM) are powerful technologies utilized to probe chromatin interactions genome wide, but how faithfully they capture three-dimensional (3D) contacts and how they perform relative to each other is unclear, as no benchmark exists. Here, we compare these methods in silico in a simplified, yet controlled, framework against known 3D structures of polymer models of murine and human loci, which can recapitulate Hi-C, GAM and SPRITE experiments and multiplexed fluorescence in situ hybridization (FISH) single-molecule conformations. We find that in silico Hi-C, GAM and SPRITE bulk data are faithful to the reference 3D structures whereas single-cell data reflect strong variability among single molecules. The minimal number of cells required in replicate experiments to return statistically similar contacts is different across the technologies, being lowest in SPRITE and highest in GAM under the same conditions. Noise-to-signal levels follow an inverse power law with detection efficiency and grow with genomic distance differently among the three methods, being lowest in GAM for genomic separations >1 Mb

Multiplex-GAM: genome-wide identification of chromatin contacts yields insights overlooked by Hi-C

Technology for measuring 3D genome topology is increasingly important for studying gene regulation, for genome assembly and for mapping of genome rearrangements. Hi-C and other ligation-based methods have become routine but have specific biases. Here, we develop multiplex-GAM, a faster and more affordable version of genome architecture mapping (GAM), a ligation-free technique that maps chromatin contacts genome-wide. We perform a detailed comparison of multiplex-GAM and Hi-C using mouse embryonic stem cells. When examining the strongest contacts detected by either method, we find that only one-third of these are shared. The strongest contacts specifically found in GAM often involve ‘active’ regions, including many transcribed genes and super-enhancers, whereas in Hi-C they more often contain ‘inactive’ regions. Our work shows that active genomic regions are involved in extensive complex contacts that are currently underestimated in ligation-based approaches, and highlights the need for orthogonal advances in genome-wide contact mapping technologies

Single-cell-resolved dynamics of chromatin architecture delineate cell and regulatory states in zebrafish embryos

DNA accessibility of cis-regulatory elements (CREs) dictates transcriptional activity and drives cell differentiation during development. While many genes regulating embryonic development have been identified, the underlying CRE dynamics controlling their expression remain largely uncharacterized. To address this, we produced a multimodal resource and genomic regulatory map for the zebrafish community, which integrates single-cell combinatorial indexing assay for transposase-accessible chromatin with high-throughput sequencing (sci-ATAC-seq) with bulk histone PTMs and Hi-C data to achieve a genome-wide classification of the regulatory architecture determining transcriptional activity in the 24-h post-fertilization (hpf) embryo. We characterized the genome-wide chromatin architecture at bulk and single-cell resolution, applying sci-ATAC-seq on whole 24-hpf stage zebrafish embryos, generating accessibility profiles for ∼23,000 single nuclei. We developed a genome segmentation method, ScregSeg (single-cell regulatory landscape segmentation), for defining regulatory programs, and candidate CREs, specific to one or more cell types. We integrated the ScregSeg output with bulk measurements for histone post-translational modifications and 3D genome organization and identified new regulatory principles between chromatin modalities prevalent during zebrafish development. Sci-ATAC-seq profiling of npas4l/cloche mutant embryos identified novel cellular roles for this hematovascular transcriptional master regulator and suggests an intricate mechanism regulating its expression. Our work defines regulatory architecture and principles in the zebrafish embryo and establishes a resource of cell-type-specific genome-wide regulatory annotations and candidate CREs, providing a valuable open resource for genomics, developmental, molecular, and computational biology

Methods for mapping 3D chromosome architecture

Determining how chromosomes are positioned and folded within the nucleus is critical to understanding the role of chromatin topology in gene regulation. Several methods are available for studying chromosome architecture, each with different strengths and limitations. Established imaging approaches and proximity ligation-based chromosome conformation capture (3C) techniques (such as DNA-FISH and Hi-C, respectively) have revealed the existence of chromosome territories, functional nuclear landmarks (such as splicing speckles and the nuclear lamina) and topologically associating domains. Improvements to these methods and the recent development of ligation-free approaches, including GAM, SPRITE and ChIA-Drop, are now helping to uncover new aspects of 3D genome topology that confirm the nucleus to be a complex, highly organized organelle