Chromatin architecture reorganization during stem cell differentiation

The Three-Dimensional Organization of Mammalian Genomes

Animal development depends on not only the linear genome sequence that embeds millions of cis-regulatory elements, but also the three-dimensional (3D) chromatin architecture that orchestrates the interplay between cis-regulatory elements and their target genes. Compared to our knowledge of the cis-regulatory sequences, the understanding of the 3D genome organization in human and other eukaryotes is still limited. Recent advances in technologies to map the 3D genome architecture have greatly accelerated the pace of discovery. Here, we review emerging concepts of chromatin organization in mammalian cells, discuss the dynamics of chromatin conformation during development, and highlight important roles for chromatin organization in cancer and other human diseases.

Between form and function: the complexity of genome folding

It has been known for over a century that chromatin is not randomly distributed within the nucleus. However, the question of how DNA is folded and the influence of such folding on nuclear processes remain topics of intensive current research. A longstanding, unanswered question is whether nuclear organization is simply a reflection of nuclear processes such as transcription and replication, or whether chromatin is folded by independent mechanisms and this per se encodes function? Evidence is emerging that both may be true. Here, using the α-globin gene cluster as an illustrative model, we provide an overview of the most recent insights into the layers of genome organization across different scales and how this relates to gene activity.

Chromosome topology guides the Drosophila Dosage Compensation Complex for target gene activation

X chromosome dosage compensation in Drosophila requires chromosome-wide coordination of gene activation. The male-specific lethal dosage compensation complex (DCC) identifies and binds to X-chromosomal high-affinity sites (HAS) from which it boosts transcription. A sub-class of HAS, PionX sites, represent first contacts on the X. Here, we explored the chromosomal interactions of representative PionX sites by high-resolution 4C and determined the global chromosome conformation by Hi-C in sex-sorted embryos. Male and female X chromosomes display similar nuclear architecture, concordant with clustered, constitutively active genes. PionX sites, like HAS, are evenly distributed in the active compartment and engage in short- and long-range interactions beyond compartment boundaries. Long-range, inter-domain interactions between DCC binding sites are stronger in males, suggesting that the complex refines chromatin organization. By de novo induction of DCC in female cells, we monitored the extent of activation surrounding PionX sites. This revealed a remarkable range of DCC action not only in linear proximity, but also at megabase distance if close in space, suggesting that DCC profits from pre-existing chromosome folding to activate genes.

Hierarchical role for transcription factors and chromatin structure in genome organization along adipogenesis

The three dimensional folding of mammalian genomes is cell type specific and difficult to alter suggesting that it is an important component of gene regulation. However, given the multitude of chromatin-associating factors, the mechanisms driving the colocalization of active chromosomal domains and the role of this organization in regulating the transcription program in adipocytes are not clear. Analysis of genome-wide chromosomal associations revealed cell type-specific spatial clustering of adipogenic genes in 3T3-L1 cells. Time course analysis demonstrated that the adipogenic 'hub', sampled by PPARγ and Lpin1, undergoes orchestrated reorganization during adipogenesis. Coupling the dynamics of genome architecture with multiple chromatin datasets indicated that among all the transcription factors (TFs) tested, RXR is central to genome reorganization at the beginning of adipogenesis. Interestingly, at the end of differentiation, the adipogenic hub was shifted to an H3K27me3-repressive environment in conjunction with attenuation of gene transcription. We propose a stage-specific hierarchy for the activity of TFs contributing to the establishment of an adipogenic genome architecture that brings together the adipogenic genetic program. In addition, the repositioning of this network in a H3K27me3-rich environment at the end of differentiation may contribute to the stabilization of gene transcription levels and reduce the developmental plasticity of these specialized cells.

DATABASE

All sequence data reported in this paper have been deposited at GEO (http://www.ncbi.nlm.nih.gov/geo/) (GSE92475).

Genome-Wide Association Study to Find Modifiers for Tetralogy of Fallot in the 22q11.2 Deletion Syndrome Identifies Variants in the GPR98 Locus on 5q14.3

BACKGROUND

The 22q11.2 deletion syndrome (22q11.2DS; DiGeorge syndrome/velocardiofacial syndrome) occurs in 1 of 4000 live births, and 60% to 70% of affected individuals have congenital heart disease, ranging from mild to severe. In our cohort of 1472 subjects with 22q11.2DS, a total of 62% (n=906) have congenital heart disease and 36% (n=326) of these have tetralogy of Fallot (TOF), comprising the largest subset of severe congenital heart disease in the cohort.

METHODS AND RESULTS

To identify common genetic variants associated with TOF in individuals with 22q11.2DS, we performed a genome-wide association study using Affymetrix 6.0 array and imputed genotype data. In our cohort, TOF was significantly associated with a genotyped single-nucleotide polymorphism (rs12519770, P=2.98×10-8) in an intron of the adhesion GPR98 (G-protein-coupled receptor V1) gene on chromosome 5q14.3. There was also suggestive evidence of association between TOF and several additional single-nucleotide polymorphisms in this region. Some genome-wide significant loci in introns or noncoding regions could affect regulation of genes nearby or at a distance. On the basis of this possibility, we examined existing Hi-C chromatin conformation data to identify genes that might be under shared transcriptional regulation within the region on 5q14.3. There are 6 genes in a topologically associated domain of chromatin with GPR98, including MEF2C (Myocyte-specific enhancer factor 2C). MEF2C is the only gene that is known to affect heart development in mammals and might be of interest with respect to 22q11.2DS.

CONCLUSIONS

In conclusion, common variants may contribute to TOF in 22q11.2DS and may function in cardiac outflow tract development.

Enhancer connectome in primary human cells identifies target genes of disease-associated DNA elements

The challenge of linking intergenic mutations to target genes has limited molecular understanding of human diseases. Here we show that H3K27ac HiChIP generates high-resolution contact maps of active enhancers and target genes in rare primary human T cell subtypes and coronary artery smooth muscle cells. Differentiation of naive T cells into T helper 17 cells or regulatory T cells creates subtype-specific enhancer-promoter interactions, specifically at regions of shared DNA accessibility. These data provide a principled means of assigning molecular functions to autoimmune and cardiovascular disease risk variants, linking hundreds of noncoding variants to putative gene targets. Target genes identified with HiChIP are further supported by CRISPR interference and activation at linked enhancers, by the presence of expression quantitative trait loci, and by allele-specific enhancer loops in patient-derived primary cells. The majority of disease-associated enhancers contact genes beyond the nearest gene in the linear genome, leading to a fourfold increase in the number of potential target genes for autoimmune and cardiovascular diseases.

Topological organization and dynamic regulation of human tRNA genes during macrophage differentiation

Background

The human genome is hierarchically organized into local and long-range structures that help shape cell-type-specific transcription patterns. Transfer RNA (tRNA) genes (tDNAs), which are transcribed by RNA polymerase III (RNAPIII) and encode RNA molecules responsible for translation, are dispersed throughout the genome and, in many cases, linearly organized into genomic clusters with other tDNAs. Whether the location and three-dimensional organization of tDNAs contribute to the activity of these genes has remained difficult to address, due in part to unique challenges related to tRNA sequencing. We therefore devised integrated tDNA expression profiling, a method that combines RNAPIII mapping with biotin-capture of nascent tRNAs. We apply this method to the study of dynamic tRNA gene regulation during macrophage development and further integrate these data with high-resolution maps of 3D chromatin structure.

Results

Integrated tDNA expression profiling reveals domain-level and loop-based organization of tRNA gene transcription during cellular differentiation. tRNA genes connected by DNA loops, which are proximal to CTCF binding sites and expressed at elevated levels compared to non-loop tDNAs, change coordinately with tDNAs and protein-coding genes at distal ends of interactions mapped by in situ Hi-C. We find that downregulated tRNA genes are specifically marked by enhanced promoter-proximal binding of MAF1, a transcriptional repressor of RNAPIII activity, altogether revealing multiple levels of tDNA regulation during cellular differentiation.

Conclusions

We present evidence of both local and coordinated long-range regulation of human tDNA expression, suggesting the location and organization of tRNA genes contribute to dynamic tDNA activity during macrophage development.

Electronic supplementary material

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

Genome Instability and γH2AX

γH2AX has emerged in the last 20 years as a central player in the DDR (DNA damage response), with specificity for DSBs (double-strand breaks). Upon the generation of DSBs, γ-phosphorylation extends along megabase-long domains in chromatin, both sides of the damage. The significance of this mechanism is of great importance; it depicts a biological amplification mechanism where one DSB induces the γ-phosphorylation of thousands of H2AX molecules along megabaselong domains of chromatin, that are adjusted to the sites of DSBs. A sequential recruitment of signal transduction factors that interact to each other and become activated to further amplify the signal that will travel to the cytoplasm take place on the γ-phosphorylated chromatin. γ-phosphorylation is an early event in the DSB damage response, induced in all phases of the cell cycle, and participates in both DSB repair pathways, the HR (homologous recombination) and NHEJ (non-homologous end joining). Today, numerous studies support the notion that γH2AX functions as a guardian of the genome by preventing misrepaired DSB that increase the mutation load of the cells and may further lead to genome instability and carcinogenesis.

Comparative analysis of three-dimensional chromosomal architecture identifies a novel fetal hemoglobin regulatory element

In this study, Huang et al. compared the chromosomal architectures of fetal and adult human erythroblasts and found that, globally, chromatin structures and compartments A/B are highly similar at both developmental stages. Their results uncover a new critical regulatory region as a potential target for therapeutic genome editing for hemoglobinopathies and highlight the power of chromosome conformation analysis in discovering new cis control elements.

Chromatin structure is tightly intertwined with transcription regulation. Here we compared the chromosomal architectures of fetal and adult human erythroblasts and found that, globally, chromatin structures and compartments A/B are highly similar at both developmental stages. At a finer scale, we detected distinct folding patterns at the developmentally controlled β-globin locus. Specifically, new fetal stage-specific contacts were uncovered between a region separating the fetal (γ) and adult (δ and β) globin genes (encompassing the HBBP1 and BGLT3 noncoding genes) and two distal chromosomal sites (HS5 and 3′HS1) that flank the locus. In contrast, in adult cells, the HBBP1–BGLT3 region contacts the embryonic ε-globin gene, physically separating the fetal globin genes from the enhancer (locus control region [LCR]). Deletion of the HBBP1 region in adult cells alters contact landscapes in ways more closely resembling those of fetal cells, including increased LCR–γ-globin contacts. These changes are accompanied by strong increases in γ-globin transcription. Notably, the effects of HBBP1 removal on chromatin architecture and gene expression closely mimic those of deleting the fetal globin repressor BCL11A, implicating BCL11A in the function of the HBBP1 region. Our results uncover a new critical regulatory region as a potential target for therapeutic genome editing for hemoglobinopathies and highlight the power of chromosome conformation analysis in discovering new cis control elements.

IC-Finder: inferring robustly the hierarchical organization of chromatin folding

The spatial organization of the genome plays a crucial role in the regulation of gene expression. Recent experimental techniques like Hi-C have emphasized the segmentation of genomes into interaction compartments that constitute conserved functional domains participating in the maintenance of a proper cell identity. Here, we propose a novel method, IC-Finder, to identify interaction compartments (IC) from experimental Hi-C maps. IC-Finder is based on a hierarchical clustering approach that we adapted to account for the polymeric nature of chromatin. Based on a benchmark of realistic in silico Hi-C maps, we show that IC-Finder is one of the best methods in terms of reliability and is the most efficient numerically. IC-Finder proposes two original options: a probabilistic description of the inferred compartments and the possibility to explore the various hierarchies of chromatin organization. Applying the method to experimental data in fly and human, we show how the predicted segmentation may depend on the normalization scheme and how 3D compartmentalization is tightly associated with epigenomic information. IC-Finder provides a robust and generic ‘all-in-one’ tool to uncover the general principles of 3D chromatin folding and their influence on gene regulation. The software is available at http://membres-timc.imag.fr/Daniel.Jost/DJ-TIMC/Software.html.

Tracking the evolution of 3D gene organization demonstrates its connection to phenotypic divergence

It has recently been shown that the organization of genes in eukaryotic genomes, and specifically in 3D, is strongly related to gene expression and function and partially conserved between organisms. However, previous studies of 3D genomic organization analyzed each organism independently from others. Here, we propose an approach for unified inter-organismal analysis of gene organization based on a network representation of Hi-C data. We define and detect four classes of spatially co-evolving orthologous modules (SCOMs), i.e. gene families that co-evolve in their 3D organization, based on patterns of divergence and conservation of distances. We demonstrate our methodology on Hi-C data from Saccharomyces cerevisiae and Schizosaccharomyces pombe, and identify, among others, modules relating to RNA splicing machinery and chromatin silencing by small RNA which are central to S. pombe's lifestyle. Our results emphasize the importance of 3D genomic organization in eukaryotes and suggest that the evolutionary mechanisms that shape gene organization affect the organism fitness and phenotypes. The proposed algorithms can be utilized in future studies of genome evolution and comparative analysis of spatial genomic organization in different tissues, conditions and single cells.

Dynamic chromatin accessibility modeled by Markov process of randomly-moving molecules in the 3D genome

Chromatin three-dimensional (3D) structure plays critical roles in gene expression regulation by influencing locus interactions and accessibility of chromatin regions. Here we propose a Markov process model to derive a chromosomal equilibrium distribution of randomly-moving molecules as a functional consequence of spatially organized genome 3D structures. The model calculates steady-state distributions (SSD) from Hi-C data as quantitative measures of each chromatin region's dynamic accessibility for transcription factors and histone modification enzymes. Different from other Hi-C derived features such as compartment A/B and interaction hubs, or traditional methods measuring chromatin accessibility such as DNase-seq and FAIRE-seq, SSD considers both chromatin–chromatin and protein–chromatin interactions. Through our model, we find that SSD could capture the chromosomal equilibrium distributions of activation histone modifications and transcription factors. Compared with compartment A/B, SSD has higher correlations with the binding of these histone modifications and transcription factors. In addition, we find that genes located in high SSD regions tend to be expressed at higher level. Furthermore, we track the change of genome organization during stem cell differentiation, and propose a two-stage model to explain the dynamic change of SSD and gene expression during differentiation, where chromatin organization genes first gain chromatin accessibility and are expressed before lineage-specific genes do. We conclude that SSD is a novel and better measure of dynamic chromatin activity and accessibility.

Advances in Genomic Profiling and Analysis of 3D Chromatin Structure and Interaction

Recent sequence-based profiling technologies such as high-throughput sequencing to detect fragment nucleotide sequence (Hi-C) and chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) have revolutionized the field of three-dimensional (3D) chromatin architecture. It is now recognized that human genome functions as folded 3D chromatin units and looping paradigm is the basic principle of gene regulation. To better interpret the 3D data dramatically accumulating in past five years and to gain deep biological insights, huge efforts have been made in developing novel quantitative analysis methods. However, the full understanding of genome regulation requires thorough knowledge in both genomic technologies and their related data analyses. We summarize the recent advances in genomic technologies in identifying the 3D chromatin structure and interaction, and illustrate the quantitative analysis methods to infer functional domains and chromatin interactions, and further elucidate the emerging single-cell Hi-C technique and its computational analysis, and finally discuss the future directions such as advances of 3D chromatin techniques in diseases.

Static and Dynamic DNA Loops form AP-1-Bound Activation Hubs during Macrophage Development

The three-dimensional arrangement of the human genome comprises a complex network of structural and regulatory chromatin loops important for coordinating changes in transcription during human development. To better understand the mechanisms underlying context-specific 3D chromatin structure and transcription during cellular differentiation, we generated comprehensive in situ Hi-C maps of DNA loops in human monocytes and differentiated macrophages. We demonstrate that dynamic looping events are regulatory rather than structural in nature and uncover widespread coordination of dynamic enhancer activity at preformed and acquired DNA loops. Enhancer-bound loop formation and enhancer activation of preformed loops together form multi-loop activation hubs at key macrophage genes. Activation hubs connect 3.4 enhancers per promoter and exhibit a strong enrichment for activator protein 1 (AP-1)-binding events, suggesting that multi-loop activation hubs involving cell-type-specific transcription factors represent an important class of regulatory chromatin structures for the spatiotemporal control of transcription.

A Scalable Bayesian Method for Integrating Functional Information in Genome-wide Association Studies

Genome-wide association studies (GWASs) have identified many complex loci. However, most loci reside in noncoding regions and have unknown biological functions. Integrative analysis that incorporates known functional information into GWASs can help elucidate the underlying biological mechanisms and prioritize important functional variants. Hence, we develop a flexible Bayesian variable selection model with efficient computational techniques for such integrative analysis. Different from previous approaches, our method models the effect-size distribution and probability of causality for variants with different annotations and jointly models genome-wide variants to account for linkage disequilibrium (LD), thus prioritizing associations based on the quantification of the annotations and allowing for multiple associated variants per locus. Our method dramatically improves both computational speed and posterior sampling convergence by taking advantage of the block-wise LD structures in human genomes. In simulations, our method accurately quantifies the functional enrichment and performs more powerfully for prioritizing the true associations than alternative methods, where the power gain is especially apparent when multiple associated variants in LD reside in the same locus. We applied our method to an in-depth GWAS of age-related macular degeneration with 33,976 individuals and 9,857,286 variants. We find the strongest enrichment for causality among non-synonymous variants (54× more likely to be causal, 1.4× larger effect sizes) and variants in transcription, repressed Polycomb, and enhancer regions, as well as identify five additional candidate loci beyond the 32 known AMD risk loci. In conclusion, our method is shown to efficiently integrate functional information in GWASs, helping identify functional associated-variants and underlying biology.

Dynamic Rewiring of Promoter-Anchored Chromatin Loops during Adipocyte Differentiation

Interactions between transcriptional promoters and their distal regulatory elements play an important role in transcriptional regulation; however, the extent to which these interactions are subject to rapid modulations in response to signals is unknown. Here, we use promoter capture Hi-C to demonstrate a rapid reorganization of promoter-anchored chromatin loops within 4 hr after inducing differentiation of 3T3-L1 preadipocytes. The establishment of new promoter-enhancer loops is tightly coupled to activation of poised (histone H3 lysine 4 mono- and dimethylated) enhancers, as evidenced by the acquisition of histone H3 lysine 27 acetylation and the binding of MED1, SMC1, and P300 proteins to these regions, as well as to activation of target genes. Intriguingly, formation of loops connecting activated enhancers and promoters is also associated with extensive recruitment of corepressors such as NCoR and HDACs, indicating that this class of coregulators may play a previously unrecognized role during enhancer activation.

Topologically associating domains are ancient features that coincide with Metazoan clusters of extreme noncoding conservation

Developmental genes in metazoan genomes are surrounded by dense clusters of conserved noncoding elements (CNEs). CNEs exhibit unexplained extreme levels of sequence conservation, with many acting as developmental long-range enhancers. Clusters of CNEs define the span of regulatory inputs for many important developmental regulators and have been described previously as genomic regulatory blocks (GRBs). Their function and distribution around important regulatory genes raises the question of how they relate to 3D conformation of these loci. Here, we show that clusters of CNEs strongly coincide with topological organisation, predicting the boundaries of hundreds of topologically associating domains (TADs) in human and Drosophila. The set of TADs that are associated with high levels of noncoding conservation exhibit distinct properties compared to TADs devoid of extreme noncoding conservation. The close correspondence between extreme noncoding conservation and TADs suggests that these TADs are ancient, revealing a regulatory architecture conserved over hundreds of millions of years.

Metazoan genomes contain many clusters of conserved noncoding elements. Here, the authors provide evidence that these clusters coincide with distinct topologically associating domains in humans and Drosophila, revealing a conserved regulatory genomic architecture.

Polycomb Repressive Complex 1 Generates Discrete Compacted Domains that Change during Differentiation

Master regulatory genes require stable silencing by the polycomb group (PcG) to prevent misexpression during differentiation and development. Some PcG proteins covalently modify histones, which contributes to heritable repression. The role for other effects on chromatin structure is less understood. We characterized the organization of PcG target genes in ESCs and neural progenitors using 5C and super-resolution microscopy. The genomic loci of repressed PcG targets formed discrete, small (20-140 Kb) domains of tight interaction that corresponded to locations bound by canonical polycomb repressive complex 1 (PRC1). These domains changed during differentiation as PRC1 binding changed. Their formation depended upon the Polyhomeotic component of canonical PRC1 and occurred independently of PRC1-catalyzed ubiquitylation. PRC1 domains differ from topologically associating domains in size and boundary characteristics. These domains have the potential to play a key role in transmitting epigenetic silencing of PcG targets by linking PRC1 to formation of a repressive higher-order structure.

5C analysis of the Epidermal Differentiation Complex locus reveals distinct chromatin interaction networks between gene-rich and gene-poor TADs in skin epithelial cells

Mammalian genomes contain several dozens of large (>0.5 Mbp) lineage-specific gene loci harbouring functionally related genes. However, spatial chromatin folding, organization of the enhancer-promoter networks and their relevance to Topologically Associating Domains (TADs) in these loci remain poorly understood. TADs are principle units of the genome folding and represents the DNA regions within which DNA interacts more frequently and less frequently across the TAD boundary. Here, we used Chromatin Conformation Capture Carbon Copy (5C) technology to characterize spatial chromatin interaction network in the 3.1 Mb Epidermal Differentiation Complex (EDC) locus harbouring 61 functionally related genes that show lineage-specific activation during terminal keratinocyte differentiation in the epidermis. 5C data validated by 3D-FISH demonstrate that the EDC locus is organized into several TADs showing distinct lineage-specific chromatin interaction networks based on their transcription activity and the gene-rich or gene-poor status. Correlation of the 5C results with genome-wide studies for enhancer-specific histone modifications (H3K4me1 and H3K27ac) revealed that the majority of spatial chromatin interactions that involves the gene-rich TADs at the EDC locus in keratinocytes include both intra- and inter-TAD interaction networks, connecting gene promoters and enhancers. Compared to thymocytes in which the EDC locus is mostly transcriptionally inactive, these interactions were found to be keratinocyte-specific. In keratinocytes, the promoter-enhancer anchoring regions in the gene-rich transcriptionally active TADs are enriched for the binding of chromatin architectural proteins CTCF, Rad21 and chromatin remodeler Brg1. In contrast to gene-rich TADs, gene-poor TADs show preferential spatial contacts with each other, do not contain active enhancers and show decreased binding of CTCF, Rad21 and Brg1 in keratinocytes. Thus, spatial interactions between gene promoters and enhancers at the multi-TAD EDC locus in skin epithelial cells are cell type-specific and involve extensive contacts within TADs as well as between different gene-rich TADs, forming the framework for lineage-specific transcription.

Blood monocyte transcriptome and epigenome analyses reveal loci associated with human atherosclerosis

Little is known regarding the epigenetic basis of atherosclerosis. Here we present the CD14+ blood monocyte transcriptome and epigenome signatures associated with human atherosclerosis. The transcriptome signature includes transcription coactivator, ARID5B, which is known to form a chromatin derepressor complex with a histone H3K9Me2-specific demethylase and promote adipogenesis and smooth muscle development. ARID5B CpG (cg25953130) methylation is inversely associated with both ARID5B expression and atherosclerosis, consistent with this CpG residing in an ARID5B enhancer region, based on chromatin capture and histone marks data. Mediation analysis supports assumptions that ARID5B expression mediates effects of cg25953130 methylation and several cardiovascular disease risk factors on atherosclerotic burden. In lipopolysaccharide-stimulated human THP1 monocytes, ARID5B knockdown reduced expression of genes involved in atherosclerosis-related inflammatory and lipid metabolism pathways, and inhibited cell migration and phagocytosis. These data suggest that ARID5B expression, possibly regulated by an epigenetically controlled enhancer, promotes atherosclerosis by dysregulating immunometabolism towards a chronic inflammatory phenotype.The molecular mechanisms mediating the impact of environmental factors in atherosclerosis are unclear. Here, the authors examine CD14+ blood monocyte's transcriptome and epigenome signatures to find differential methylation and expression of ARID5B to be associated with human atherosclerosis.

Genome organization: connecting the developmental origins of disease and genetic variation

An adverse early life environment can increase the risk of metabolic and other disorders later in life. Genetic variation can modify an individual's susceptibility to these environmental challenges. These gene by environment interactions are important, but difficult, to dissect. The nucleus is the primary organelle where environmental responses impact directly on the genetic variants within the genome, resulting in changes to the biology of the genome and ultimately the phenotype. Understanding genome biology requires the integration of the linear DNA sequence, epigenetic modifications and nuclear proteins that are present within the nucleus. The interactions between these layers of information may be captured in the emergent spatial genome organization. As such genome organization represents a key research area for decoding the role of genetic variation in the Developmental Origins of Health and Disease.

Gene functioning and storage within a folded genome

In mammals, genomic DNA that is roughly 2 m long is folded to fit the size of the cell nucleus that has a diameter of about 10 μm. The folding of genomic DNA is mediated via assembly of DNA-protein complex, chromatin. In addition to the reduction of genomic DNA linear dimensions, the assembly of chromatin allows to discriminate and to mark active (transcribed) and repressed (non-transcribed) genes. Consequently, epigenetic regulation of gene expression occurs at the level of DNA packaging in chromatin. Taking into account the increasing attention of scientific community toward epigenetic systems of gene regulation, it is very important to understand how DNA folding in chromatin is related to gene activity. For many years the hierarchical model of DNA folding was the most popular. It was assumed that nucleosome fiber (10-nm fiber) is folded into 30-nm fiber and further on into chromatin loops attached to a nuclear/chromosome scaffold. Recent studies have demonstrated that there is much less regularity in chromatin folding within the cell nucleus. The very existence of 30-nm chromatin fibers in living cells was questioned. On the other hand, it was found that chromosomes are partitioned into self-interacting spatial domains that restrict the area of enhancers action. Thus, TADs can be considered as structural-functional domains of the chromosomes. Here we discuss the modern view of DNA packaging within the cell nucleus in relation to the regulation of gene expression. Special attention is paid to the possible mechanisms of the chromatin fiber self-assembly into TADs. We discuss the model postulating that partitioning of the chromosome into TADs is determined by the distribution of active and inactive chromatin segments along the chromosome.

This article was specially invited by the editors and represents work by leading researchers.

Topological Domains, Metagenes, and the Emergence of Pleiotropic Regulations at Hox Loci

Hox gene clusters of jaw vertebrates display a tight genomic organization, which has no equivalent in any other bilateria genomes sequenced thus far. It was previously argued that such a topological consolidation toward a condensed, metagenic structure was due to the accumulation of long-range regulations flanking Hox loci, a phenomenon made possible by the successive genome duplications that occurred at the root of the vertebrate lineage, similar to gene neofunctionalization but applied to a coordinated multigenic system. Here, we propose that the emergence of such large vertebrate regulatory landscapes containing a range of global enhancers was greatly facilitated by the presence of topologically associating domains (TADs). These chromatin domains, mostly constitutive, may have been used as genomic niches where novel regulations could evolve due to both the preexistence of a structural backbone poised to integrate novel regulatory inputs, and a highly adaptive transcriptional readout. We propose a scenario for the coevolution of such TADs and the emergence of pleiotropy at ancestral vertebrate Hox loci.

Chromatin Domains: The Unit of Chromosome Organization

How eukaryotic chromosomes fold inside the nucleus is an age-old question that remains unanswered today. Early biochemical and microscopic studies revealed the existence of chromatin domains and loops as a pervasive feature of interphase chromosomes, but the biological implications of such organizational features were obscure. Genome-wide analysis of pair-wise chromatin interactions using chromatin conformation capture (3C)-based techniques has shed new light on the organization of chromosomes in interphase nuclei. Particularly, the finding of cell-type invariant, evolutionarily conserved topologically associating domains (TADs) in a broad spectrum of cell types has provided a new molecular framework for the study of animal development and human diseases. Here, we review recent progress in characterization of such chromatin domains and delineation of mechanisms of their formation in animal cells.

Lineage-specific dynamic and pre-established enhancer-promoter contacts cooperate in terminal differentiation

Chromosome conformation is an important feature of metazoan gene regulation; however, enhancer-promoter contact remodeling during cellular differentiation remains poorly understood. To address this, genome-wide promoter capture Hi-C (CHi-C) was performed during epidermal differentiation. Two classes of enhancer-promoter contacts associated with differentiation-induced genes were identified. The first class ('gained') increased in contact strength during differentiation in concert with enhancer acquisition of the H3K27ac activation mark. The second class ('stable') were pre-established in undifferentiated cells, with enhancers constitutively marked by H3K27ac. The stable class was associated with the canonical conformation regulator cohesin, whereas the gained class was not, implying distinct mechanisms of contact formation and regulation. Analysis of stable enhancers identified a new, essential role for a constitutively expressed, lineage-restricted ETS-family transcription factor, EHF, in epidermal differentiation. Furthermore, neither class of contacts was observed in pluripotent cells, suggesting that lineage-specific chromatin structure is established in tissue progenitor cells and is further remodeled in terminal differentiation.

The interplay of epigenetic marks during stem cell differentiation and development

Chromatin, the template for epigenetic regulation, is a highly dynamic entity that is constantly reshaped during early development and differentiation. Epigenetic modification of chromatin provides the necessary plasticity for cells to respond to environmental and positional cues, and enables the maintenance of acquired information without changing the DNA sequence. The mechanisms involve, among others, chemical modifications of chromatin, changes in chromatin constituents and reconfiguration of chromatin interactions and 3D structure. New advances in genome-wide technologies have paved the way towards an integrative view of epigenome dynamics during cell state transitions, and recent findings in embryonic stem cells highlight how the interplay between different epigenetic layers reshapes the transcriptional landscape.

"Cat's Cradling" the 3D Genome by the Act of LncRNA Transcription

There is growing evidence that transcription and nuclear organization are tightly linked. Yet, whether transcription of thousands of long noncoding RNAs (lncRNAs) could play a role in this packaging process remains elusive. Although some lncRNAs have been found to have clear roles in nuclear architecture (e.g., FIRRE, NEAT1, XIST, and others), the vast majority remain poorly understood. In this Perspective, we highlight how the act of transcription can affect nuclear architecture. We synthesize several recent findings into a proposed model where the transcription of lncRNAs can serve as guide-posts for shaping genome organization. This model is similar to the game "cat's cradle," where the shape of a string is successively changed by opening up new sites for finger placement. Analogously, transcription of lncRNAs could serve as "grip holds" for nuclear proteins to pull the genome into new positions. This model could explain general lncRNA properties such as low abundance and tissue specificity. Overall, we propose a general framework for how the act of lncRNA transcription could play a role in organizing the 3D genome.

A critical assessment of topologically associating domain prediction tools

Topologically associating domains (TADs) have been proposed to be the basic unit of chromosome folding and have been shown to play key roles in genome organization and gene regulation. Several different tools are available for TAD prediction, but their properties have never been thoroughly assessed. In this manuscript, we compare the output of seven different TAD prediction tools on two published Hi-C data sets. TAD predictions varied greatly between tools in number, size distribution and other biological properties. Assessed against a manual annotation of TADs, individual TAD boundary predictions were found to be quite reliable, but their assembly into complete TAD structures was much less so. In addition, many tools were sensitive to sequencing depth and resolution of the interaction frequency matrix. This manuscript provides users and designers of TAD prediction tools with information that will help guide the choice of tools and the interpretation of their predictions.

Large scale genomic reorganization of topological domains at the HoxD locus

BACKGROUND

The transcriptional activation of HoxD genes during mammalian limb development involves dynamic interactions with two topologically associating domains (TADs) flanking the HoxD cluster. In particular, the activation of the most posterior HoxD genes in developing digits is controlled by regulatory elements located in the centromeric TAD (C-DOM) through long-range contacts.

RESULTS

To assess the structure-function relationships underlying such interactions, we measured compaction levels and TAD discreteness using a combination of chromosome conformation capture (4C-seq) and DNA FISH. We assessed the robustness of the TAD architecture by using a series of genomic deletions and inversions that impact the integrity of this chromatin domain and that remodel long-range contacts. We report multi-partite associations between HoxD genes and up to three enhancers. We find that the loss of native chromatin topology leads to the remodeling of TAD structure following distinct parameters.

CONCLUSIONS

Our results reveal that the recomposition of TAD architectures after large genomic re-arrangements is dependent on a boundary-selection mechanism in which CTCF mediates the gating of long-range contacts in combination with genomic distance and sequence specificity. Accordingly, the building of a recomposed TAD at this locus depends on distinct functional and constitutive parameters.

Identification of SNPs associated with muscle yield and quality traits using allelic-imbalance analyses of pooled RNA-Seq samples in rainbow trout

Background

Coding/functional SNPs change the biological function of a gene and, therefore, could serve as “large-effect” genetic markers. In this study, we used two bioinformatics pipelines, GATK and SAMtools, for discovering coding/functional SNPs with allelic-imbalances associated with total body weight, muscle yield, muscle fat content, shear force, and whiteness. Phenotypic data were collected for approximately 500 fish, representing 98 families (5 fish/family), from a growth-selected line, and the muscle transcriptome was sequenced from 22 families with divergent phenotypes (4 low- versus 4 high-ranked families per trait).

Results

GATK detected 59,112 putative SNPs; of these SNPs, 4798 showed allelic imbalances (>2.0 as an amplification and <0.5 as loss of heterozygosity). SAMtools detected 87,066 putative SNPs; and of them, 4962 had allelic imbalances between the low- and high-ranked families. Only 1829 SNPs with allelic imbalances were common between the two datasets, indicating significant differences in algorithms. The two datasets contained 7930 non-redundant SNPs of which 4439 mapped to 1498 protein-coding genes (with 6.4% non-synonymous SNPs) and 684 mapped to 295 lncRNAs. Validation of a subset of 92 SNPs revealed 1) 86.7–93.8% success rate in calling polymorphic SNPs and 2) 95.4% consistent matching between DNA and cDNA genotypes indicating a high rate of identifying SNPs with allelic imbalances. In addition, 4.64% SNPs revealed random monoallelic expression. Genome distribution of the SNPs with allelic imbalances exhibited high density for all five traits in several chromosomes, especially chromosome 9, 20 and 28. Most of the SNP-harboring genes were assigned to important growth-related metabolic pathways.

Conclusion

These results demonstrate utility of RNA-Seq in assessing phenotype-associated allelic imbalances in pooled RNA-Seq samples. The SNPs identified in this study were included in a new SNP-Chip design (available from Affymetrix) for genomic and genetic analyses in rainbow trout.

Electronic supplementary material

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

Computational Prediction of Position Effects of Apparently Balanced Human Chromosomal Rearrangements

Interpretation of variants of uncertain significance, especially chromosomal rearrangements in non-coding regions of the human genome, remains one of the biggest challenges in modern molecular diagnosis. To improve our understanding and interpretation of such variants, we used high-resolution three-dimensional chromosomal structural data and transcriptional regulatory information to predict position effects and their association with pathogenic phenotypes in 17 subjects with apparently balanced chromosomal abnormalities. We found that the rearrangements predict disruption of long-range chromatin interactions between several enhancers and genes whose annotated clinical features are strongly associated with the subjects' phenotypes. We confirm gene-expression changes for a couple of candidate genes to exemplify the utility of our analysis of position effect. These results highlight the important interplay between chromosomal structure and disease and demonstrate the need to utilize chromatin conformational data for the prediction of position effects in the clinical interpretation of non-coding chromosomal rearrangements.

3D Chromatin Structures of Mature Gametes and Structural Reprogramming during Mammalian Embryogenesis

High-order chromatin structure plays important roles in gene expression regulation. Knowledge of the dynamics of 3D chromatin structures during mammalian embryo development remains limited. We report the 3D chromatin architecture of mouse gametes and early embryos using an optimized Hi-C method with low-cell samples. We find that mature oocytes at the metaphase II stage do not have topologically associated domains (TADs). In sperm, extra-long-range interactions (>4 Mb) and interchromosomal interactions occur frequently. The high-order structures of both the paternal and maternal genomes in zygotes and two-cell embryos are obscure but are gradually re-established through development. The establishment of the TAD structure requires DNA replication but not zygotic genome activation. Furthermore, unmethylated CpGs are enriched in A compartment, and methylation levels are decreased to a greater extent in A compartment than in B compartment in embryos. In summary, the global reprogramming of chromatin architecture occurs during early mammalian development.

Local Genome Topology Can Exhibit an Incompletely Rewired 3D-Folding State during Somatic Cell Reprogramming

Pluripotent genomes are folded in a topological hierarchy that reorganizes during differentiation. The extent to which chromatin architecture is reconfigured during somatic cell reprogramming is poorly understood. Here we integrate fine-resolution architecture maps with epigenetic marks and gene expression in embryonic stem cells (ESCs), neural progenitor cells (NPCs), and NPC-derived induced pluripotent stem cells (iPSCs). We find that most pluripotency genes reconnect to target enhancers during reprogramming. Unexpectedly, some NPC interactions around pluripotency genes persist in our iPSC clone. Pluripotency genes engaged in both "fully-reprogrammed" and "persistent-NPC" interactions exhibit over/undershooting of target expression levels in iPSCs. Additionally, we identify a subset of "poorly reprogrammed" interactions that do not reconnect in iPSCs and display only partially recovered, ESC-specific CTCF occupancy. 2i/LIF can abrogate persistent-NPC interactions, recover poorly reprogrammed interactions, reinstate CTCF occupancy, and restore expression levels. Our results demonstrate that iPSC genomes can exhibit imperfectly rewired 3D-folding linked to inaccurately reprogrammed gene expression.

Activation of the alpha-globin gene expression correlates with dramatic upregulation of nearby non-globin genes and changes in local and large-scale chromatin spatial structure

BACKGROUND

In homeotherms, the alpha-globin gene clusters are located within permanently open genome regions enriched in housekeeping genes. Terminal erythroid differentiation results in dramatic upregulation of alpha-globin genes making their expression comparable to the rRNA transcriptional output. Little is known about the influence of the erythroid-specific alpha-globin gene transcription outburst on adjacent, widely expressed genes and large-scale chromatin organization. Here, we have analyzed the total transcription output, the overall chromatin contact profile, and CTCF binding within the 2.7 Mb segment of chicken chromosome 14 harboring the alpha-globin gene cluster in cultured lymphoid cells and cultured erythroid cells before and after induction of terminal erythroid differentiation.

RESULTS

We found that, similarly to mammalian genome, the chicken genomes is organized in TADs and compartments. Full activation of the alpha-globin gene transcription in differentiated erythroid cells is correlated with upregulation of several adjacent housekeeping genes and the emergence of abundant intergenic transcription. An extended chromosome region encompassing the alpha-globin cluster becomes significantly decompacted in differentiated erythroid cells, and depleted in CTCF binding and CTCF-anchored chromatin loops, while the sub-TAD harboring alpha-globin gene cluster and the upstream major regulatory element (MRE) becomes highly enriched with chromatin interactions as compared to lymphoid and proliferating erythroid cells. The alpha-globin gene domain and the neighboring loci reside within the A-like chromatin compartment in both lymphoid and erythroid cells and become further segregated from the upstream gene desert upon terminal erythroid differentiation.

CONCLUSIONS

Our findings demonstrate that the effects of tissue-specific transcription activation are not restricted to the host genomic locus but affect the overall chromatin structure and transcriptional output of the encompassing topologically associating domain.

Connected Gene Communities Underlie Transcriptional Changes in Cornelia de Lange Syndrome

Cornelia de Lange syndrome (CdLS) is a complex multisystem developmental disorder caused by mutations in cohesin subunits and regulators. While its precise molecular mechanisms are not well defined, they point toward a global deregulation of the transcriptional gene expression program. Cohesin is associated with the boundaries of chromosome domains and with enhancer and promoter regions connecting the three-dimensional genome organization with transcriptional regulation. Here, we show that connected gene communities, structures emerging from the interactions of noncoding regulatory elements and genes in the three-dimensional chromosomal space, provide a molecular explanation for the pathoetiology of CdLS associated with mutations in the cohesin-loading factor NIPBL and the cohesin subunit SMC1A NIPBL and cohesin are important constituents of connected gene communities that are centrally positioned at noncoding regulatory elements. Accordingly, genes deregulated in CdLS are positioned within reach of NIPBL- and cohesin-occupied regions through promoter-promoter interactions. Our findings suggest a dynamic model where NIPBL loads cohesin to connect genes in communities, offering an explanation for the gene expression deregulation in the CdLS.

Comparison of computational methods for Hi-C data analysis

Hi-C is a genome-wide sequencing technique used to investigate 3D chromatin conformation inside the nucleus. Computational methods are required to analyze Hi-C data and identify chromatin interactions and topologically associating domains (TADs) from genome-wide contact probability maps. We quantitatively compared the performance of 13 algorithms in their analyses of Hi-C data from six landmark studies and simulations. This comparison revealed differences in the performance of methods for chromatin interaction identification, but more comparable results for TAD detection between algorithms.

Constrained release of lamina-associated enhancers and genes from the nuclear envelope during T-cell activation facilitates their association in chromosome compartments

The 3D organization of the genome changes concomitantly with expression changes during hematopoiesis and immune activation. Studies have focused either on lamina-associated domains (LADs) or on topologically associated domains (TADs), defined by preferential local chromatin interactions, and chromosome compartments, defined as higher-order interactions between TADs sharing functionally similar states. However, few studies have investigated how these affect one another. To address this, we mapped LADs using Lamin B1-DamID during Jurkat T-cell activation, finding significant genome reorganization at the nuclear periphery dominated by release of loci frequently important for T-cell function. To assess how these changes at the nuclear periphery influence wider genome organization, our DamID data sets were contrasted with TADs and compartments. Features of specific repositioning events were then tested by fluorescence in situ hybridization during T-cell activation. First, considerable overlap between TADs and LADs was observed with the TAD repositioning as a unit. Second, A1 and A2 subcompartments are segregated in 3D space through differences in proximity to LADs along chromosomes. Third, genes and a putative enhancer in LADs that were released from the periphery during T-cell activation became preferentially associated with A2 subcompartments and were constrained to the relative proximity of the lamina. Thus, lamina associations influence internal nuclear organization, and changes in LADs during T-cell activation may provide an important additional mode of gene regulation.

Engineering cell identity: establishing new gene regulatory and chromatin landscapes

Cellular reprogramming can be achieved by ectopically expressing transcription factors that directly convert one differentiated cell type into another, bypassing embryonic states. A number of different cell types have been generated by such 'direct lineage reprogramming' methods, but their practical utility has been limited because, in most protocols, the resulting populations are often partially differentiated or incompletely specified. Here, we review mechanisms of lineage reprogramming by pioneer transcription factors, a unique class of transcriptional regulators that has the capacity to engage with silent chromatin to activate target gene regulatory networks. We assess the possible barriers to successful reprogramming in the context of higher-order chromatin landscape, considering how the mechanistic relationship between nuclear organization and cell identity will be crucial to unlocking the full potential of cell fate engineering.

Chromatin Architecture Emerges during Zygotic Genome Activation Independent of Transcription

Chromatin architecture is fundamental in regulating gene expression. To investigate when spatial genome organization is first established during development, we examined chromatin conformation during Drosophila embryogenesis and observed the emergence of chromatin architecture within a tight time window that coincides with the onset of transcription activation in the zygote. Prior to zygotic genome activation, the genome is mostly unstructured. Early expressed genes serve as nucleation sites for topologically associating domain (TAD) boundaries. Activation of gene expression coincides with the establishment of TADs throughout the genome and co-localization of housekeeping gene clusters, which remain stable in subsequent stages of development. However, the appearance of TAD boundaries is independent of transcription and requires the transcription factor Zelda for locus-specific TAD boundary insulation. These results offer insight into when spatial organization of the genome emerges and identify a key factor that helps trigger this architecture.

Chromatin Control of Developmental Dynamics and Plasticity

Chromatin structure is intimately connected with gene expression and cell identity. Here we review recent advances in the field and discuss how establishment of cell identity during development is accompanied by large-scale remodeling of the epigenetic landscape and how this remodeling drives and supports lineage specification and maintenance. We discuss maternal control of the early embryonic epigenetic landscape, selective usage of enhancer clusters via 3D chromatin contacts leading to activation of transcription factor networks, and conserved regulation of developmental pathways by specific DNA demethylation of key regulatory regions. Together, these processes establish an epigenetic framework regulating different phases of embryonic development.

Regulatory landscape fusion in rhabdomyosarcoma through interactions between the PAX3 promoter and FOXO1 regulatory elements

Background

The organisation of vertebrate genomes into topologically associating domains (TADs) is believed to facilitate the regulation of the genes located within them. A remaining question is whether TAD organisation is achieved through the interactions of the regulatory elements within them or if these interactions are favoured by the pre-existence of TADs. If the latter is true, the fusion of two independent TADs should result in the rewiring of the transcriptional landscape and the generation of ectopic contacts.

Results

We show that interactions within the PAX3 and FOXO1 domains are restricted to their respective TADs in normal conditions, while in a patient-derived alveolar rhabdomyosarcoma cell line, harbouring the diagnostic t(2;13)(q35;q14) translocation that brings together the PAX3 and FOXO1 genes, the PAX3 promoter interacts ectopically with FOXO1 sequences. Using a combination of 4C-seq datasets, we have modelled the three-dimensional organisation of the fused landscape in alveolar rhabdomyosarcoma.

Conclusions

The chromosomal translocation that leads to alveolar rhabdomyosarcoma development generates a novel TAD that is likely to favour ectopic PAX3:FOXO1 oncogene activation in non-PAX3 territories. Rhabdomyosarcomas may therefore arise from cells which do not normally express PAX3. The borders of this novel TAD correspond to the original 5'- and 3'- borders of the PAX3 and FOXO1 TADs, respectively, suggesting that TAD organisation precedes the formation of regulatory long-range interactions. Our results demonstrate that, upon translocation, novel regulatory landscapes are formed allowing new intra-TAD interactions between the original loci involved.

Electronic supplementary material

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

Revealing the vectors of cellular identity with single-cell genomics

Single-cell genomics has now made it possible to create a comprehensive atlas of human cells. At the same time, it has reopened definitions of a cell's identity and of the ways in which identity is regulated by the cell's molecular circuitry. Emerging computational analysis methods, especially in single-cell RNA sequencing (scRNA-seq), have already begun to reveal, in a data-driven way, the diverse simultaneous facets of a cell's identity, from discrete cell types to continuous dynamic transitions and spatial locations. These developments will eventually allow a cell to be represented as a superposition of 'basis vectors', each determining a different (but possibly dependent) aspect of cellular organization and function. However, computational methods must also overcome considerable challenges-from handling technical noise and data scale to forming new abstractions of biology. As the scale of single-cell experiments continues to increase, new computational approaches will be essential for constructing and characterizing a reference map of cell identities.

lncRNA-screen: an interactive platform for computationally screening long non-coding RNAs in large genomics datasets

Background

Long non-coding RNAs (lncRNAs) have emerged as a class of factors that are important for regulating development and cancer. Computational prediction of lncRNAs from ultra-deep RNA sequencing has been successful in identifying candidate lncRNAs. However, the complexity of handling and integrating different types of genomics data poses significant challenges to experimental laboratories that lack extensive genomics expertise.

Result

To address this issue, we have developed lncRNA-screen, a comprehensive pipeline for computationally screening putative lncRNA transcripts over large multimodal datasets. The main objective of this work is to facilitate the computational discovery of lncRNA candidates to be further examined by functional experiments. lncRNA-screen provides a fully automated easy-to-run pipeline which performs data download, RNA-seq alignment, assembly, quality assessment, transcript filtration, novel lncRNA identification, coding potential estimation, expression level quantification, histone mark enrichment profile integration, differential expression analysis, annotation with other type of segmented data (CNVs, SNPs, Hi-C, etc.) and visualization. Importantly, lncRNA-screen generates an interactive report summarizing all interesting lncRNA features including genome browser snapshots and lncRNA-mRNA interactions based on Hi-C data.

Conclusion

lncRNA-screen provides a comprehensive solution for lncRNA discovery and an intuitive interactive report for identifying promising lncRNA candidates. lncRNA-screen is available as open-source software on GitHub.

Electronic supplementary material

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

Physical properties of the chromosomes and implications for development

Remarkable progress has been made in understanding chromosome structures inside the cell nucleus. Recent advances in Hi-C technologies enable the detection of genome-wide chromatin interactions, providing insight into three-dimensional (3D) genome organization. Advancements in the spatial and temporal resolutions of imaging as well as in molecular biological techniques allow the tracking of specific chromosomal loci, improving our understanding of chromosome movements. From these data, we are beginning to understand how the intra-nuclear locations of chromatin loci and the 3D genome structure change during development and differentiation. This emerging field of genome structure and dynamics research requires an interdisciplinary approach including efficient collaborations between experimental biologists and physicists, informaticians, or engineers. Quantitative and mathematical analyses based on polymer physics are becoming increasingly important for processing and interpreting experimental data on 3D chromosome structures and dynamics. In this review, we aim to provide an overview of recent research on the physical aspects of chromosome structure and dynamics oriented for biologists. These studies have mainly focused on chromosomes at the cellular level, using unicellular organisms and cultured cells. However, physical parameters that change during development, such as nuclear size, may impact genome structure and dynamics. Here, we discuss how chromatin dynamics and genome structures in early embryos change during development, which we expect will be a hot topic in the field of chromatin dynamics in the near future. We hope this review helps developmental biologists to quantitatively investigate the physical natures of chromosomes in developmental biology research.

Coupling between chromosome intermingling and gene regulation during cellular differentiation

In this article, we summarize current findings for the emergence of biophysical properties such as nuclear stiffness, chromatin compaction, chromosome positioning, and chromosome intermingling during stem cell differentiation, which eventually correlated with the changes of gene expression profiles during cellular differentiation. An overview is first provided to link stem cell differentiation with alterations in nuclear architecture, chromatin compaction, along with nuclear and chromatin dynamics. Further, we highlight the recent biophysical and molecular approaches, imaging methods and computational developments in characterizing transcription-related chromosome organization especially chromosome intermingling and nano-scale chromosomal contacts. Finally, the article ends with an outlook towards the emergence of a functional roadmap in setting up chromosome positioning and intermingling in a cell type specific manner during cellular differentiation.

YY1 and CTCF orchestrate a 3D chromatin looping switch during early neural lineage commitment

CTCF is an architectural protein with a critical role in connecting higher-order chromatin folding in pluripotent stem cells. Recent reports have suggested that CTCF binding is more dynamic during development than previously appreciated. Here, we set out to understand the extent to which shifts in genome-wide CTCF occupancy contribute to the 3D reconfiguration of fine-scale chromatin folding during early neural lineage commitment. Unexpectedly, we observe a sharp decrease in CTCF occupancy during the transition from naïve/primed pluripotency to multipotent primary neural progenitor cells (NPCs). Many pluripotency gene-enhancer interactions are anchored by CTCF, and its occupancy is lost in parallel with loop decommissioning during differentiation. Conversely, CTCF binding sites in NPCs are largely preexisting in pluripotent stem cells. Only a small number of CTCF sites arise de novo in NPCs. We identify another zinc finger protein, Yin Yang 1 (YY1), at the base of looping interactions between NPC-specific genes and enhancers. Putative NPC-specific enhancers exhibit strong YY1 signal when engaged in 3D contacts and negligible YY1 signal when not in loops. Moreover, siRNA knockdown of Yy1 specifically disrupts interactions between key NPC enhancers and their target genes. YY1-mediated interactions between NPC regulatory elements are often nested within constitutive loops anchored by CTCF. Together, our results support a model in which YY1 acts as an architectural protein to connect developmentally regulated looping interactions; the location of YY1-mediated interactions may be demarcated in development by a preexisting topological framework created by constitutive CTCF-mediated interactions.

The second decade of 3C technologies: detailed insights into nuclear organization

The relevance of three-dimensional (3D) genome organization for transcriptional regulation and thereby for cellular fate at large is now widely accepted. Our understanding of the fascinating architecture underlying this function is based on microscopy studies as well as the chromosome conformation capture (3C) methods, which entered the stage at the beginning of the millennium. The first decade of 3C methods rendered unprecedented insights into genome topology. Here, we provide an update of developments and discoveries made over the more recent years. As we discuss, established and newly developed experimental and computational methods enabled identification of novel, functionally important chromosome structures. Regulatory and architectural chromatin loops throughout the genome are being cataloged and compared between cell types, revealing tissue invariant and developmentally dynamic loops. Architectural proteins shaping the genome were disclosed, and their mode of action is being uncovered. We explain how more detailed insights into the 3D genome increase our understanding of transcriptional regulation in development and misregulation in disease. Finally, to help researchers in choosing the approach best tailored for their specific research question, we explain the differences and commonalities between the various 3C-derived methods.

The Cohesin Release Factor WAPL Restricts Chromatin Loop Extension

The spatial organization of chromosomes influences many nuclear processes including gene expression. The cohesin complex shapes the 3D genome by looping together CTCF sites along chromosomes. We show here that chromatin loop size can be increased and that the duration with which cohesin embraces DNA determines the degree to which loops are enlarged. Cohesin's DNA release factor WAPL restricts this loop extension and also prevents looping between incorrectly oriented CTCF sites. We reveal that the SCC2/SCC4 complex promotes the extension of chromatin loops and the formation of topologically associated domains (TADs). Our data support the model that cohesin structures chromosomes through the processive enlargement of loops and that TADs reflect polyclonal collections of loops in the making. Finally, we find that whereas cohesin promotes chromosomal looping, it rather limits nuclear compartmentalization. We conclude that the balanced activity of SCC2/SCC4 and WAPL enables cohesin to correctly structure chromosomes.

Intrinsic Molecular Processes: Impact on Mutagenesis

Mutations provide resources for genome evolution by generating genetic variability. In addition, mutations act as a driving force leading to disease pathogenesis, and thus have important implications for disease diagnosis, prognosis, and treatment. Understanding the mechanisms underlying how mutations occur is therefore of prime importance for elucidating evolutionary and pathogenic processes. Recent genomics studies have revealed that mutations occur non-randomly across the human genome. In particular, the distribution of mutations is highly associated with intrinsic molecular processes including transcription, chromatin organization, DNA replication timing, and DNA repair. Interplay between intrinsic processes and extrinsic mutagenic exposure may thus imprint a characteristic mutational landscape on tumors. We discuss the impact of intrinsic molecular processes on mutation acquisition in cancer.