The triploid endosperm genome of Arabidopsis adopts a peculiar, parental-dosage-dependent chromatin organization

The endosperm is a seed tissue unique to flowering plants. Due to its central role in nourishing and protecting the embryo, endosperm development is subject to parental conflicts and adaptive processes, which led to the evolution of parent-of-origin-dependent gene regulation. The role of higher-order chromatin organization in regulating the endosperm genome was long ignored due to technical hindrance. We developed a combination of approaches to analyze nuclear structure and chromatin organization in Arabidopsis thaliana endosperm. Endosperm nuclei showed a less condensed chromatin than other types of nuclei and a peculiar heterochromatin organization, with smaller chromocenters and additional heterochromatic foci interspersed in euchromatin. This is accompanied by a redistribution of the heterochromatin mark H3K9me1 from chromocenters toward euchromatin and interspersed heterochromatin. Thus, endosperm nuclei have a specific nuclear architecture and organization, which we interpret as a relaxed chromocenter-loop model. The analysis of endosperm with altered parental genome dosage indicated that the additional heterochromatin may be predominantly of maternal origin, suggesting differential regulation of maternal and paternal genomes, possibly linked to genome dosage regulation.     

The INTACT method for cell type-specific gene expression and chromatin profiling in Arabidopsis thaliana

Genomic studies of cell differentiation and function within a whole organism depend on the ability to isolate specific cell types from a tissue, but this is technically difficult. We developed a method called INTACT (isolation of nuclei tagged in specific cell types) that allows affinity-based isolation of nuclei from individual cell types of a tissue, thereby circumventing the problems associated with mechanical purification techniques. In this method nuclei are affinity-labeled through transgenic expression of a biotinylated nuclear envelope protein in the cell type of interest. Total nuclei are isolated from transgenic plants and biotin-labeled nuclei are then purified using streptavidin-coated magnetic beads, without the need for specialized equipment. INTACT gives high yield and purity of nuclei from the desired cell types, which can be used for genome-wide analysis of gene expression and chromatin features. The entire procedure, from nuclei purification through cDNA preparation or chromatin immunoprecipitation (ChIP), can be completed within 2 d. The protocol we present assumes that transgenic lines are already available, and includes procedural details for amplification of cDNA or ChIP DNA prior to microarray or deep sequencing analysis.     

Integrative modeling reveals the principles of multi-scale chromatin boundary formation in human nuclear organization

BackgroundInterphase chromosomes adopt a hierarchical structure, and recent data have characterized their chromatin organization at very different scales, from sub-genic regions associated with DNA-binding proteins at the order of tens or hundreds of bases, through larger regions with active or repressed chromatin states, up to multi-megabase-scale domains associated with nuclear positioning, replication timing and other qualities. However, we have lacked detailed, quantitative models to understand the interactions between these different strata.ResultsHere we collate large collections of matched locus-level chromatin features and Hi-C interaction data, representing higher-order organization, across three human cell types. We use quantitative modeling approaches to assess whether locus-level features are sufficient to explain higher-order structure, and identify the most influential underlying features. We identify structurally variable domains between cell types and examine the underlying features to discover a general association with cell-type-specific enhancer activity. We also identify the most prominent features marking the boundaries of two types of higher-order domains at different scales: topologically associating domains and nuclear compartments. We find parallel enrichments of particular chromatin features for both types, including features associated with active promoters and the architectural proteins CTCF and YY1.ConclusionsWe show that integrative modeling of large chromatin dataset collections using random forests can generate useful insights into chromosome structure. The models produced recapitulate known biological features of the cell types involved, allow exploration of the antecedents of higher-order structures and generate testable hypotheses for further experimental studies.

Electronic supplementary material

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

MaxHiC: A robust background correction model to identify biologically relevant chromatin interactions in Hi-C and capture Hi-C experiments

Hi-C is a genome-wide chromosome conformation capture technology that detects interactions between pairs of genomic regions and exploits higher order chromatin structures. Conceptually Hi-C data counts interaction frequencies between every position in the genome and every other position. Biologically functional interactions are expected to occur more frequently than transient background and artefactual interactions. To identify biologically relevant interactions, several background models that take biases such as distance, GC content and mappability into account have been proposed. Here we introduce MaxHiC, a background correction tool that deals with these complex biases and robustly identifies statistically significant interactions in both Hi-C and capture Hi-C experiments. MaxHiC uses a negative binomial distribution model and a maximum likelihood technique to correct biases in both Hi-C and capture Hi-C libraries. We systematically benchmark MaxHiC against major Hi-C background correction tools including Hi-C significant interaction callers (SIC) and Hi-C loop callers using published Hi-C, capture Hi-C, and Micro-C datasets. Our results demonstrate that 1) Interacting regions identified by MaxHiC have significantly greater levels of overlap with known regulatory features (e.g. active chromatin histone marks, CTCF binding sites, DNase sensitivity) and also disease-associated genome-wide association SNPs than those identified by currently existing models, 2) the pairs of interacting regions are more likely to be linked by eQTL pairs and 3) more likely to link known regulatory features including known functional enhancer-promoter pairs validated by CRISPRi than any of the existing methods. We also demonstrate that interactions between different genomic region types have distinct distance distributions only revealed by MaxHiC. MaxHiC is publicly available as a python package for the analysis of Hi-C, capture Hi-C and Micro-C data.     

THUNDER: A reference-free deconvolution method to infer cell type proportions from bulk Hi-C data

Hi-C data provide population averaged estimates of three-dimensional chromatin contacts across cell types and states in bulk samples. Effective analysis of Hi-C data entails controlling for the potential confounding factor of differential cell type proportions across heterogeneous bulk samples. We propose a novel unsupervised deconvolution method for inferring cell type composition from bulk Hi-C data, the Two-step Hi-c UNsupervised DEconvolution appRoach (THUNDER). We conducted extensive simulations to test THUNDER based on combining two published single-cell Hi-C (scHi-C) datasets. THUNDER more accurately estimates the underlying cell type proportions compared to reference-free methods (e.g., TOAST, and NMF) and is more robust than reference-dependent methods (e.g. MuSiC). We further demonstrate the practical utility of THUNDER to estimate cell type proportions and identify cell-type-specific interactions in Hi-C data from adult human cortex tissue samples. THUNDER will be a useful tool in adjusting for varying cell type composition in population samples, facilitating valid and more powerful downstream analysis such as differential chromatin organization studies. Additionally, THUNDER estimated contact profiles provide a useful exploratory framework to investigate cell-type-specificity of the chromatin interactome while experimental data is still rare.     

A unified framework for inferring the multi-scale organization of chromatin domains from Hi-C

Chromosomes are giant chain molecules organized into an ensemble of three-dimensional structures characterized with its genomic state and the corresponding biological functions. Despite the strong cell-to-cell heterogeneity, the cell-type specific pattern demonstrated in high-throughput chromosome conformation capture (Hi-C) data hints at a valuable link between structure and function, which makes inference of chromatin domains (CDs) from the pattern of Hi-C a central problem in genome research. Here we present a unified method for analyzing Hi-C data to determine spatial organization of CDs over multiple genomic scales. By applying statistical physics-based clustering analysis to a polymer physics model of the chromosome, our method identifies the CDs that best represent the global pattern of correlation manifested in Hi-C. The multi-scale intra-chromosomal structures compared across different cell types uncover the principles underlying the multi-scale organization of chromatin chain: (i) Sub-TADs, TADs, and meta-TADs constitute a robust hierarchical structure. (ii) The assemblies of compartments and TAD-based domains are governed by different organizational principles. (iii) Sub-TADs are the common building blocks of chromosome architecture. Our physically principled interpretation and analysis of Hi-C not only offer an accurate and quantitative view of multi-scale chromatin organization but also help decipher its connections with genome function.     

Dynamic analyses of the expression of the HISTONE::YFP fusion protein in arabidopsis show that syncytial endosperm is divided in mitotic domains

During early seed development, nuclear divisions in the endosperm are not followed by cell division, leading to the development of a syncytium. The simple organization of the Arabidopsis endosperm provides a model in which to study the regulation of the cell cycle in relation to development. To monitor nuclear divisions, we constructed a HISTONE 2B::YELLOW FLUORESCENT PROTEIN gene fusion (H2B::YFP). To validate its use as a vital marker for chromatin in plants, H2B::YFP was expressed constitutively in Arabidopsis. This enabled the observation of mitoses in living root meristems. H2B::YFP was expressed specifically in Arabidopsis syncytial endosperm by using GAL4 transactivation. Monitoring mitotic activity in living syncytial endosperm showed that the syncytium was organized into three domains in which nuclei divide simultaneously with a specific time course. Each mitotic domain has a distinct spatiotemporal pattern of mitotic CYCLIN B1;1 accumulation. The polar spatial organization of the three mitotic domains suggests interactions between developmental mechanisms and the regulation of the cell cycle.     

The 4D Nucleome: Genome Compartmentalization in an Evolutionary Context

4D nucleome research aims to understand the impact of nuclear organization in space and time on nuclear functions, such as gene expression patterns, chromatin replication, and the maintenance of genome integrity. In this review we describe evidence that the origin of 4D genome compartmentalization can be traced back to the prokaryotic world. In cell nuclei of animals and plants chromosomes occupy distinct territories, built up from ~1 Mb chromatin domains, which in turn are composed of smaller chromatin subdomains and also form larger chromatin domain clusters. Microscopic evidence for this higher order chromatin landscape was strengthened by chromosome conformation capture studies, in particular Hi-C. This approach demonstrated ~1 Mb sized, topologically associating domains in mammalian cell nuclei separated by boundaries. Mutations, which destroy boundaries, can result in developmental disorders and cancer. Nucleosomes appeared first as tetramers in the Archaea kingdom and later evolved to octamers built up each from two H2A, two H2B, two H3, and two H4 proteins. Notably, nucleosomes were lost during the evolution of the Dinoflagellata phylum. Dinoflagellate chromosomes remain condensed during the entire cell cycle, but their chromosome architecture differs radically from the architecture of other eukaryotes. In summary, the conservation of fundamental features of higher order chromatin arrangements throughout the evolution of metazoan animals suggests the existence of conserved, but still unknown mechanism(s) controlling this architecture. Notwithstanding this conservation, a comparison of metazoans and protists also demonstrates species-specific structural and functional features of nuclear organization.     

Minimalistic 3D chromatin models: Sparse interactions in single cells drive the chromatin fold and form many-body units

Computational three-dimensional chromatin modeling has helped uncover principles of genome organization. Here, we discuss methods for modeling three-dimensional chromatin structures, with focus on a minimalistic polymer model which inverts population Hi-C into single-cell conformations. Utilizing only basic physical properties, this model reveals that a few specific Hi-C interactions can fold chromatin into conformations consistent with single-cell imaging, Dip-C, and FISH measurements. Aggregated single-cell chromatin conformations also reproduce Hi-C frequencies. This approach allows quantification of structural heterogeneity and discovery of many-body interaction units and has revealed additional insights, including (1) topologically associating domains as a byproduct of folding driven by specific interactions, (2) cell subpopulations with different structural scaffolds are developmental stage dependent, and (3) the functional landscape of many-body units within enhancer-rich regions. We also discuss these findings in relation to the genome structure–function relationship.     

Higher-order folding of heterochromatin: protein bridges span the nucleosome arrays.

In interphase eukaryotic nuclei, chromatin is divided into two morphologically distinct types known as heterochromatin and euchromatin. It has been long suggested that the two types of chromatin differ at the level of higher-order folding. Recent studies have revealed the features of chromatin 3D architecture that distinguish the higher-order folding of repressed and active chromatin and have identified chromosomal proteins and their modifications associated with these structural transitions. This review discusses the molecular and structural determinants of chromatin higher-order folding in relation to mechanism(s) of heterochromatin formation and genetic silencing during cell differentiation and tissue development.     

HiCImpute: A Bayesian hierarchical model for identifying structural zeros and enhancing single cell Hi-C data

Single cell Hi-C techniques enable one to study cell to cell variability in chromatin interactions. However, single cell Hi-C (scHi-C) data suffer severely from sparsity, that is, the existence of excess zeros due to insufficient sequencing depth. Complicating the matter further is the fact that not all zeros are created equal: some are due to loci truly not interacting because of the underlying biological mechanism (structural zeros); others are indeed due to insufficient sequencing depth (sampling zeros or dropouts), especially for loci that interact infrequently. Differentiating between structural zeros and dropouts is important since correct inference would improve downstream analyses such as clustering and discovery of subtypes. Nevertheless, distinguishing between these two types of zeros has received little attention in the single cell Hi-C literature, where the issue of sparsity has been addressed mainly as a data quality improvement problem. To fill this gap, in this paper, we propose HiCImpute, a Bayesian hierarchical model that goes beyond data quality improvement by also identifying observed zeros that are in fact structural zeros. HiCImpute takes spatial dependencies of scHi-C 2D data structure into account while also borrowing information from similar single cells and bulk data, when such are available. Through an extensive set of analyses of synthetic and real data, we demonstrate the ability of HiCImpute for identifying structural zeros with high sensitivity, and for accurate imputation of dropout values. Downstream analyses using data improved from HiCImpute yielded much more accurate clustering of cell types compared to using observed data or data improved by several comparison methods. Most significantly, HiCImpute-improved data have led to the identification of subtypes within each of the excitatory neuronal cells of L4 and L5 in the prefrontal cortex.     

Microrheology for Hi-C Data Reveals the Spectrum of the Dynamic 3D Genome Organization

The one-dimensional information of genomic DNA is hierarchically packed inside the eukaryotic cell nucleus and organized in a three-dimensional (3D) space. Genome-wide chromosome conformation capture (Hi-C) methods have uncovered the 3D genome organization and revealed multiscale chromatin domains of compartments and topologically associating domains (TADs). Moreover, single-nucleosome live-cell imaging experiments have revealed the dynamic organization of chromatin domains caused by stochastic thermal fluctuations. However, the mechanism underlying the dynamic regulation of such hierarchical and structural chromatin units within the microscale thermal medium remains unclear. Microrheology is a way to measure dynamic viscoelastic properties coupling between thermal microenvironment and mechanical response. Here, we propose a new, to our knowledge, microrheology for Hi-C data to analyze the dynamic compliance property as a measure of rigidness and flexibility of genomic regions along with the time evolution. Our method allows the conversion of an Hi-C matrix into the spectrum of the dynamic rheological property along the genomic coordinate of a single chromosome. To demonstrate the power of the technique, we analyzed Hi-C data during the neural differentiation of mouse embryonic stem cells. We found that TAD boundaries behave as more rigid nodes than the intra-TAD regions. The spectrum clearly shows the dynamic viscoelasticity of chromatin domain formation at different timescales. Furthermore, we characterized the appearance of synchronous and liquid-like intercompartment interactions in differentiated cells. Together, our microrheology data derived from Hi-C data provide physical insights into the dynamics of the 3D genome organization.     

Predicting chromatin organization using histone marks

Genome-wide mapping of three dimensional chromatin organization is an important yet technically challenging task. To aid experimental effort and to understand the determinants of long-range chromatin interactions, we have developed a computational model integrating Hi-C and histone mark ChIP-seq data to predict two important features of chromatin organization: chromatin interaction hubs and topologically associated domain (TAD) boundaries. Our model accurately and robustly predicts these features across datasets and cell types. Cell-type specific histone mark information is required for prediction of chromatin interaction hubs but not for TAD boundaries. Our predictions provide a useful guide for the exploration of chromatin organization.

Electronic supplementary material

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

Suppression of liquid–liquid phase separation by 1,6-hexanediol partially compromises the 3D genome organization in living cells

Liquid–liquid phase separation (LLPS) contributes to the spatial and functional segregation of molecular processes within the cell nucleus. However, the role played by LLPS in chromatin folding in living cells remains unclear. Here, using stochastic optical reconstruction microscopy (STORM) and Hi-C techniques, we studied the effects of 1,6-hexanediol (1,6-HD)-mediated LLPS disruption/modulation on higher-order chromatin organization in living cells. We found that 1,6-HD treatment caused the enlargement of nucleosome clutches and their more uniform distribution in the nuclear space. At a megabase-scale, chromatin underwent moderate but irreversible perturbations that resulted in the partial mixing of A and B compartments. The removal of 1,6-HD from the culture medium did not allow chromatin to acquire initial configurations, and resulted in more compact repressed chromatin than in untreated cells. 1,6-HD treatment also weakened enhancer-promoter interactions and TAD insulation but did not considerably affect CTCF-dependent loops. Our results suggest that 1,6-HD-sensitive LLPS plays a limited role in chromatin spatial organization by constraining its folding patterns and facilitating compartmentalization at different levels.     

Comparison of Hi-C results using in-solution versus in-nucleus ligation

BackgroundChromosome conformation capture and various derivative methods such as 4C, 5C and Hi-C have emerged as standard tools to analyze the three-dimensional organization of the genome in the nucleus. These methods employ ligation of diluted cross-linked chromatin complexes, intended to favor proximity-dependent, intra-complex ligation. During development of single-cell Hi-C, we devised an alternative Hi-C protocol with ligation in preserved nuclei rather than in solution. Here we directly compare Hi-C methods employing in-nucleus ligation with the standard in-solution ligation.ResultsWe show in-nucleus ligation results in consistently lower levels of inter-chromosomal contacts. Through chromatin mixing experiments we show that a significantly large fraction of inter-chromosomal contacts are the result of spurious ligation events formed during in-solution ligation. In-nucleus ligation significantly reduces this source of experimental noise, and results in improved reproducibility between replicates. We also find that in-nucleus ligation eliminates restriction fragment length bias found with in-solution ligation. These improvements result in greater reproducibility of long-range intra-chromosomal and inter-chromosomal contacts, as well as enhanced detection of structural features such as topologically associated domain boundaries.ConclusionsWe conclude that in-nucleus ligation captures chromatin interactions more consistently over a wider range of distances, and significantly reduces both experimental noise and bias. In-nucleus ligation creates higher quality Hi-C libraries while simplifying the experimental procedure. We suggest that the entire range of 3C applications are likely to show similar benefits from in-nucleus ligation.

Electronic supplementary material

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

Advances in higher-order chromatin architecture: the move towards 4D genome

In eukaryotes, the genome is hierarchically packed inside the nucleus, which facilitates physical contact between cis-regulatory elements (CREs), such as enhancers and promoters. Accumulating evidence highlights the critical role of higher-order chromatin structure in precise regulation of spatiotemporal gene expression under diverse biological contexts including lineage commitment and cell activation by external stimulus. Genomics and imaging-based technologies, such as Hi-C and DNA fluorescence in situ hybridization (FISH), have revealed the key principles of genome folding, while newly developed tools focus on improvement in resolution, throughput and modality at single-cell and population levels, and challenge the knowledge obtained through conventional approaches. In this review, we discuss recent advances in our understanding of principles of higher-order chromosome conformation and technologies to investigate 4D chromatin interactions.     

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.