Dynamics of the 4D genome during in vivo lineage specification and differentiation

CTCF depletion decouples enhancer-mediated gene activation from chromatin hub formation

Enhancers and promoters interact in three-dimensional (3D) chromatin structures to regulate gene expression. Here we characterize the mechanisms that drive the formation and function of these structures in a lymphoid-to-myeloid transdifferentiation system. Based on analyses at base pair resolution, we demonstrate a close correlation between binding of regulatory proteins, formation of chromatin interactions and gene expression. Multi-way interaction analyses and computational modeling show that tissue-specific gene loci are organized into chromatin hubs, characterized by cooperative interactions between multiple enhancers, promoters and CTCF-binding sites. While depletion of CTCF strongly impairs the formation of these chromatin hubs, the effects of CTCF depletion on gene expression are modest and can be explained by rewired enhancer-promoter interactions. These findings demonstrate a role for enhancer-promoter interactions in gene regulation that is independent of cooperative interactions in chromatin hubs. Together, these results contribute to our understanding of the structure-function relationship of the genome during cellular differentiation.

Active regulatory elements recruit cohesin to establish cell specific chromatin domains

As the 3D structure of the genome is analysed at ever increasing resolution it is clear that there is considerable variation in the 3D chromatin architecture across different cell types. It has been proposed that this may, in part, be due to increased recruitment of cohesin to activated cis-elements (enhancers and promoters) leading to cell-type specific loop extrusion underlying the formation of new sub-TADs. Here we show that cohesin correlates well with the presence of active enhancers and that this varies in an allele-specific manner with the presence or absence of polymorphic enhancers which vary from one individual to another. Using the alpha globin cluster as a model, we show that when all enhancers are removed, peaks of cohesin disappear from these regions and the erythroid specific sub-TAD is no longer formed. Re-insertion of the major alpha globin enhancer (R2) is associated with re-establishment of recruitment and increased interactions. In complementary experiments insertion of the R2 enhancer element into a "neutral" region of the genome recruits cohesin, induces transcription and creates a new large (75 kb) erythroid-specific domain. Together these findings support the proposal that active enhancers recruit cohesin, stimulate loop extrusion and promote the formation of cell specific sub-TADs.

Ancient genomic linkage of α-globin and Nprl3 couples metabolism with erythropoiesis

Red blood cell development from erythroid progenitors requires profound reshaping of metabolism and gene expression. How these transcriptional and metabolic alterations are coupled is unclear. Nprl3 (an inhibitor of mTORC1) has remained in synteny with the α-globin genes for >500 million years, and harbours most of the a-globin enhancers. However, whether Nprl3 serves an erythroid role is unknown. We found that while haematopoietic progenitors require basal Nprl3 expression, erythroid Nprl3 expression is further boosted by the α-globin enhancers. This lineage-specific upregulation is required for sufficient erythropoiesis. Loss of Nprl3 affects erythroblast metabolism via elevating mTORC1 signalling, suppressing autophagy and disrupting glycolysis. Broadly consistent with these murine findings, human NPRL3-knockout erythroid progenitors produce fewer enucleated cells and demonstrate dysregulated mTORC1 signalling in response to nutrient availability and erythropoietin. Therefore, we propose that the anciently conserved linkage of NprI3, α-globin and their associated enhancers has coupled metabolic and developmental control of erythropoiesis.

Post-mitotic transcriptional activation and 3D regulatory interactions show locus- and differentiation-specific sensitivity to cohesin depletion

Prior studies showed that structural loops collapse upon acute cohesin depletion, while regulatory enhancer-promoter (E-P) loops largely persist, consistent with minimal transcriptional changes. However, these studies, conducted in asynchronous cells, could not resolve whether cohesin is required for the establishment of regulatory interactions and transcriptional activation during cell division or cell state transitions. To address this gap, we degraded RAD21, a core cohesin subunit, in naïve mouse embryonic stem cells (ESCs) transitioning from mitosis to G1 either in self-renewal condition or during differentiation toward formative pluripotency. Although most structural loops failed to be re-established without cohesin, about 35% of regulatory loops reformed at normal or higher frequencies. Cohesin-independent loops showed characteristics of strong active enhancers and promoters and a significant association with H3K27ac mitotic bookmarks. However, inhibition of CBP/p300 during mitotic exit did not impact these cohesin-independent interactions, suggesting the presence of complex compensatory mechanisms. At the transcriptional level, cohesin depletion induced only minor changes, supporting that post-mitotic transcriptional reactivation is largely independent of cohesin. The few genes with impaired reactivation were directly bound by RAD21 at their promoters, engaged in many structural loops, and located within strongly insulated TADs with low gene density. Importantly, degrading cohesin during the M-to-G1 transition in the presence of EpiLC differentiation signals revealed a larger group of susceptible genes, including key signature genes and transcription factors. Impaired activation of these genes was partly due to the failure to establish de novo EpiLC-specific interactions in the absence of cohesin. These experiments revealed locus-specific and context-specific dependencies between cohesin, E-P interactions, and transcription.

The α-globin super-enhancer acts in an orientation-dependent manner

Individual enhancers are defined as short genomic regulatory elements, bound by transcription factors, and able to activate cell-specific gene expression at a distance, in an orientation-independent manner. Within mammalian genomes, enhancer-like elements may be found individually or within clusters referred to as locus control regions or super-enhancers (SEs). While these behave similarly to individual enhancers with respect to cell specificity, distribution and distance, their orientation-dependence has not been formally tested. Here, using the α-globin locus as a model, we show that while an individual enhancer works in an orientation-independent manner, the direction of activity of a SE changes with its orientation. When the SE is inverted within its normal chromosomal context, expression of its normal targets, the α-globin genes, is severely reduced and the normally silent genes lying upstream of the α-globin locus are upregulated. These findings add to our understanding of enhancer-promoter specificity that precisely activate transcription.

Transient promoter interactions modulate developmental gene activation

Transcriptional induction coincides with the formation of various chromatin topologies. Strong evidence supports that gene activation is accompanied by a general increase in promoter-enhancer interactions. However, it remains unclear how these topological changes are coordinated across time and space during transcriptional activation. Here, we combine chromatin conformation capture with transcription and chromatin profiling during an embryonic stem cell (ESC) differentiation time course to determine how 3D genome restructuring is related to transcriptional transitions. This approach allows us to identify distinct topological alterations that are associated with the magnitude of transcriptional induction. We detect transiently formed interactions and demonstrate by genetic deletions that associated distal regulatory elements (DREs), as well as appropriate formation and disruption of these interactions, can contribute to the transcriptional induction of linked genes. Together, our study links topological dynamics to the magnitude of transcriptional induction and detects an uncharacterized type of transcriptionally important DREs.

Functional associations of evolutionarily recent human genes exhibit sensitivity to the 3D genome landscape and disease

Genome organization is intricately tied to regulating genes and associated cell fate decisions. Here, we examine the positioning and functional significance of human genes, grouped by their lineage restriction level, within the 3D organization of the genome. We reveal that genes of different lineage restriction levels have distinct positioning relationships with both domains and loop anchors, and remarkably consistent relationships with boundaries across cell types. While the functional associations of each group of genes are primarily cell type-specific, associations of conserved genes maintain greater stability across 3D genomic features and disease than recently evolved genes. Furthermore, the expression of these genes across various tissues follows an evolutionary progression, such that RNA levels increase from young lineage restricted genes to ancient genes present in most species. Thus, the distinct relationships of gene evolutionary age, function, and positioning within 3D genomic features contribute to tissue-specific gene regulation in development and disease.

The characteristics of CTCF binding sequences contribute to enhancer blocking activity

While the elements encoding enhancers and promoters have been relatively well studied, the full spectrum of insulator elements which bind the CCCTC binding factor (CTCF), is relatively poorly characterized. This is partly due to the genomic context of CTCF sites greatly influencing their roles and activity. Here we have developed an experimental system to determine the ability of minimal, consistently sized, individual CTCF elements to interpose between enhancers and promoters and thereby reduce gene expression during differentiation. Importantly, each element is tested in the identical location thereby minimising the effect of genomic context. We found no correlation between the ability of CTCF elements to block enhancer-promoter activity with the degree of evolutionary conservation; their resemblance to the consensus core sequences; or the number of CTCF core motifs harboured in the element. Nevertheless, we have shown that the strongest enhancer-promoter blockers include a previously described bound element lying upstream of the CTCF core motif. In addition, we found other uncharacterised DNaseI footprints located close to the core motif that may affect function. We have developed an assay of CTCF sequences which will enable researchers to sub-classify individual CTCF elements in a uniform and unbiased way.

A moonlighting job for α-globin in blood vessels

ABSTRACT

Red blood cells express high levels of hemoglobin A tetramer (α2β2) to facilitate oxygen transport. Hemoglobin subunits and related proteins are also expressed at lower levels in other tissues across the animal kingdom. Physiological functions for most nonerythroid globins likely derive from their ability to catalyze reduction-oxidation (redox) reactions via electron transfer through heme-associated iron. An interesting example is illustrated by the recent discovery that α-globin without β-globin is expressed in some arteriolar endothelial cells (ECs). α-globin binds EC nitric oxide (NO) synthase (eNOS) and degrades its enzymatic product NO, a potent vasodilator. Thus, depletion of α-globin in ECs or inhibition of its association with eNOS causes arteriolar relaxation and lowering of blood pressure in mice. Some of these findings have been replicated in isolated human blood vessels, and genetic studies are tractable in populations in which α-thalassemia alleles are prevalent. Two small studies identified associations between loss of α-globin genes in humans and NO-regulated vascular responses elicited by local hypoxia-induced blood flow or thermal stimulation. In a few larger population-based studies, no associations were detected between loss of α-globin genes and blood pressure, ischemic stroke, or pulmonary hypertension. In contrast, a significant positive association between α-globin gene copy number and kidney disease was detected in an African American cohort. Further studies are required to define comprehensively the expression of α-globin in different vascular beds and ascertain their overall impact on normal and pathological vascular physiology.

scGrapHiC: deep learning-based graph deconvolution for Hi-C using single cell gene expression

SUMMARY

Single-cell Hi-C (scHi-C) protocol helps identify cell-type-specific chromatin interactions and sheds light on cell differentiation and disease progression. Despite providing crucial insights, scHi-C data is often underutilized due to the high cost and the complexity of the experimental protocol. We present a deep learning framework, scGrapHiC, that predicts pseudo-bulk scHi-C contact maps using pseudo-bulk scRNA-seq data. Specifically, scGrapHiC performs graph deconvolution to extract genome-wide single-cell interactions from a bulk Hi-C contact map using scRNA-seq as a guiding signal. Our evaluations show that scGrapHiC, trained on seven cell-type co-assay datasets, outperforms typical sequence encoder approaches. For example, scGrapHiC achieves a substantial improvement of 23.2% in recovering cell-type-specific Topologically Associating Domains over the baselines. It also generalizes to unseen embryo and brain tissue samples. scGrapHiC is a novel method to generate cell-type-specific scHi-C contact maps using widely available genomic signals that enables the study of cell-type-specific chromatin interactions.

AVAILABILITY AND IMPLEMENTATION

The GitHub link: https://github.com/rsinghlab/scGrapHiC contains the source code of scGrapHiC and associated scripts to preprocess publicly available datasets to produce the results and visualizations we have discuss in this manuscript.

FACT maintains chromatin architecture and thereby stimulates RNA polymerase II pausing during transcription in vivo

Facilitates chromatin transcription (FACT) is a histone chaperone that supports transcription through chromatin in vitro, but its functional roles in vivo remain unclear. Here, we analyze the in vivo functions of FACT with the use of multi-omics analysis after rapid FACT depletion from human cells. We show that FACT depletion destabilizes chromatin and leads to transcriptional defects, including defective promoter-proximal pausing and elongation, and increased premature termination of RNA polymerase II. Unexpectedly, our analysis revealed that promoter-proximal pausing depends not only on the negative elongation factor (NELF) but also on the +1 nucleosome, which is maintained by FACT.

Suv39h-catalyzed H3K9me3 is critical for euchromatic genome organization and the maintenance of gene transcription

H3K9me3-dependent heterochromatin is critical for the silencing of repeat-rich pericentromeric regions and also has key roles in repressing lineage-inappropriate protein-coding genes in differentiation and development. Here, we investigate the molecular consequences of heterochromatin loss in cells deficient in both SUV39H1 and SUV39H2 (Suv39DKO), the major mammalian histone methyltransferase enzymes that catalyze heterochromatic H3K9me3 deposition. We reveal a paradoxical repression of protein-coding genes in Suv39DKO cells, with these differentially expressed genes principally in euchromatic (Tn5-accessible, H3K4me3- and H3K27ac-marked) rather than heterochromatic (H3K9me3-marked) or polycomb (H3K27me3-marked) regions. Examination of the three-dimensional (3D) nucleome reveals that transcriptomic dysregulation occurs in euchromatic regions close to the nuclear periphery in 3D space. Moreover, this transcriptomic dysregulation is highly correlated with altered 3D genome organization in Suv39DKO cells. Together, our results suggest that the nuclear lamina-tethering of Suv39-dependent H3K9me3 domains provides an essential scaffold to support euchromatic genome organization and the maintenance of gene transcription for healthy cellular function.

Potential roles of inter-chromosomal interactions in cell fate determination

Mammalian genomic DNA is packed in a small nucleus, and its folding and organization in the nucleus are critical for gene regulation and cell fate determination. In interphase, chromosomes are compartmentalized into certain nuclear spaces and territories that are considered incompatible with each other. The regulation of gene expression is influenced by the epigenetic characteristics of topologically associated domains and A/B compartments within chromosomes (intrachromosomal). Previously, interactions among chromosomes detected via chromosome conformation capture-based methods were considered noise or artificial errors. However, recent studies based on newly developed ligation-independent methods have shown that inter-chromosomal interactions play important roles in gene regulation. This review summarizes the recent understanding of spatial genomic organization in mammalian interphase nuclei and discusses the potential mechanisms that determine cell identity. In addition, this review highlights the potential role of inter-chromosomal interactions in early mouse development.

The probability of chromatin to be at the nuclear lamina has no systematic effect on its transcription level in fruit flies

BACKGROUND

Multiple studies have demonstrated a negative correlation between gene expression and positioning of genes at the nuclear envelope (NE) lined by nuclear lamina, but the exact relationship remains unclear, especially in light of the highly stochastic, transient nature of the gene association with the NE.

RESULTS

In this paper, we ask whether there is a causal, systematic, genome-wide relationship between the expression levels of the groups of genes in topologically associating domains (TADs) of Drosophila nuclei and the probabilities of TADs to be found at the NE. To investigate the nature of this possible relationship, we combine a coarse-grained dynamic model of the entire Drosophila nucleus with genome-wide gene expression data; we analyze the TAD averaged transcription levels of genes against the probabilities of individual TADs to be in contact with the NE in the control and lamins-depleted nuclei. Our findings demonstrate that, within the statistical error margin, the stochastic positioning of Drosophila melanogaster TADs at the NE does not, by itself, systematically affect the mean level of gene expression in these TADs, while the expected negative correlation is confirmed. The correlation is weak and disappears completely for TADs not containing lamina-associated domains (LADs) or TADs containing LADs, considered separately. Verifiable hypotheses regarding the underlying mechanism for the presence of the correlation without causality are discussed. These include the possibility that the epigenetic marks and affinity to the NE of a TAD are determined by various non-mutually exclusive mechanisms and remain relatively stable during interphase.

CONCLUSIONS

At the level of TADs, the probability of chromatin being in contact with the nuclear envelope has no systematic, causal effect on the transcription level in Drosophila. The conclusion is reached by combining model-derived time-evolution of TAD locations within the nucleus with their experimental gene expression levels.

Pas de deux: the coordinated coupling of erythroid differentiation with the cell cycle

PURPOSE OF REVIEW

Recent work reveals that cell cycle duration and structure are remodeled in lock-step with distinct stages of erythroid differentiation. These cell cycle features have regulatory roles in differentiation, beyond the generic function of increasing cell number.

RECENT FINDINGS

Developmental progression through the early erythroid progenitor stage (known as colony-forming-erythroid, or 'CFU-e') is characterized by gradual shortening of G1 phase of the cycle. This process culminates in a key transcriptional switch to erythroid terminal differentiation (ETD) that is synchronized with, and dependent on, S phase progression. Further, the CFU-e/ETD switch takes place during an unusually short S phase, part of an exceptionally short cell cycle that is characterized by globally fast replication fork speeds. Cell cycle and S phase speed can alter developmental events during erythroid differentiation, through pathways that are targeted by glucocorticoid and erythropoietin signaling during the erythroid stress response.

SUMMARY

There is close inter-dependence between cell cycle structure and duration, S phase and replication fork speeds, and erythroid differentiation stage. Further, modulation of cell cycle structure and speed cycle impacts developmental progression and cell fate decisions during erythroid differentiation. These pathways may offer novel mechanistic insights and potential therapeutic targets.

Rewiring of the epigenome and chromatin architecture by exogenously induced retinoic acid signaling during zebrafish embryonic development

Retinoic acid (RA) is the ligand of RA receptors (RARs), transcription factors that bind to RA response elements. RA signaling is required for multiple processes during embryonic development, including body axis extension, hindbrain antero-posterior patterning and forelimb bud initiation. Although some RA target genes have been identified, little is known about the genome-wide effects of RA signaling during in vivo embryonic development. Here, we stimulate the RA pathway by treating zebrafish embryos with all-trans-RA (atRA) and use a combination of RNA-seq, ATAC-seq, ChIP-seq and HiChIP to gain insight into the molecular mechanisms by which exogenously induced RA signaling controls gene expression. We find that RA signaling is involved in anterior/posterior patterning, central nervous system development, and the transition from pluripotency to differentiation. AtRA treatment also alters chromatin accessibility during early development and promotes chromatin binding of RARαa and the RA targets Hoxb1b, Meis2b and Sox3, which cooperate in central nervous system development. Finally, we show that exogenous RA induces a rewiring of chromatin architecture, with alterations in chromatin 3D interactions involving target genes. Altogether, our findings identify genome-wide targets of RA signaling and provide a molecular mechanism by which developmental signaling pathways regulate target gene expression by altering chromatin topology.

Physical and functional interaction among Irf8 enhancers during dendritic cell differentiation

The production of type 1 conventional dendritic cells (cDC1s) requires high expression of the transcription factor IRF8. Three enhancers at the Irf8 3' region function in a differentiation stage-specific manner. However, whether and how these enhancers interact physically and functionally remains unclear. Here, we show that the Irf8 3' enhancers directly interact with each other and contact the Irf8 gene body during cDC1 differentiation. The +56 kb enhancer, which functions from multipotent progenitor stages, activates the other 3' enhancers through an IRF8-dependent transcription factor program, that is, in trans. Then, the +32 kb enhancer, which operates in cDC1-committed cells, reversely acts in cis on the other 3' enhancers to maintain the high expression of Irf8. Indeed, mice with compound heterozygous deletion of the +56 and +32 kb enhancers are unable to generate cDC1s. These results illustrate how multiple enhancers cooperate to induce a lineage-determining transcription factor gene during cell differentiation.

Contributions of transcriptional noise to leukaemia evolution: KAT2A as a case-study

Transcriptional noise is proposed to participate in cell fate changes, but contributions to mammalian cell differentiation systems, including cancer, remain associative. Cancer evolution is driven by genetic variability, with modulatory or contributory participation of epigenetic variants. Accumulation of epigenetic variants enhances transcriptional noise, which can facilitate cancer cell fate transitions. Acute myeloid leukaemia (AML) is an aggressive cancer with strong epigenetic dependencies, characterized by blocked differentiation. It constitutes an attractive model to probe links between transcriptional noise and malignant cell fate regulation. Gcn5/KAT2A is a classical epigenetic transcriptional noise regulator. Its loss increases transcriptional noise and modifies cell fates in stem and AML cells. By reviewing the analysis of KAT2A-depleted pre-leukaemia and leukaemia models, I discuss that the net result of transcriptional noise is diversification of cell fates secondary to alternative transcriptional programmes. Cellular diversification can enable or hinder AML progression, respectively, by differentiation of cell types responsive to mutations, or by maladaptation of leukaemia stem cells. KAT2A-dependent noise-responsive genes participate in ribosome biogenesis and KAT2A loss destabilizes translational activity. I discuss putative contributions of perturbed translation to AML biology, and propose KAT2A loss as a model for mechanistic integration of transcriptional and translational control of noise and fate decisions. This article is part of a discussion meeting issue 'Causes and consequences of stochastic processes in development and disease'.

Enhancer-promoter interactions become more instructive in the transition from cell-fate specification to tissue differentiation

To regulate expression, enhancers must come in proximity to their target gene. However, the relationship between the timing of enhancer-promoter (E-P) proximity and activity remains unclear, with examples of uncoupled, anticorrelated and correlated interactions. To assess this, we selected 600 characterized enhancers or promoters with tissue-specific activity in Drosophila embryos and performed Capture-C in FACS-purified myogenic or neurogenic cells during specification and tissue differentiation. This enabled direct comparison between E-P proximity and activity transitioning from OFF-to-ON and ON-to-OFF states across developmental conditions. This showed remarkably similar E-P topologies between specified muscle and neuronal cells, which are uncoupled from activity. During tissue differentiation, many new distal interactions emerge where changes in E-P proximity reflect changes in activity. The mode of E-P regulation therefore appears to change as embryogenesis proceeds, from largely permissive topologies during cell-fate specification to more instructive regulation during terminal tissue differentiation, when E-P proximity is coupled to activation.

Chromatin organization of muscle stem cell

The proper functioning of skeletal muscles is essential throughout life. A crucial crosstalk between the environment and several cellular mechanisms allows striated muscles to perform successfully. Notably, the skeletal muscle tissue reacts to an injury producing a completely functioning tissue. The muscle's robust regenerative capacity relies on the fine coordination between muscle stem cells (MuSCs or "satellite cells") and their specific microenvironment that dictates stem cells' activation, differentiation, and self-renewal. Critical for the muscle stem cell pool is a fine regulation of chromatin organization and gene expression. Acquiring a lineage-specific 3D genome architecture constitutes a crucial modulator of muscle stem cell function during development, in the adult stage, in physiological and pathological conditions. The context-dependent relationship between genome structure, such as accessibility and chromatin compartmentalization, and their functional effects will be analysed considering the improved 3D epigenome knowledge, underlining the intimate liaison between environmental encounters and epigenetics.

Shaping gene expression and its evolution by chromatin architecture and enhancer activity

Transcriptional regulation plays a pivotal role in orchestrating the intricate genetic programs governing embryonic development. The expression of developmental genes relies on the combined activity of several cis-regulatory elements (CREs), such as enhancers and silencers, which can be located at long linear distances from the genes that they regulate and that interact with them through establishment of chromatin loops. Mutations affecting their activity or interaction with their target genes can lead to developmental disorders and are thought to have importantly contributed to the evolution of the animal body plan. The income of next-generation-sequencing approaches has allowed identifying over a million of sequences with putative regulatory potential in the human genome. Characterizing their function and establishing gene-CREs maps is essential to decode the logic governing developmental gene expression and is one of the major challenges of the post-genomic era. Chromatin 3D organization plays an essential role in determining how CREs specifically contact their target genes while avoiding deleterious off-target interactions. Our understanding of these aspects has greatly advanced with the income of chromatin conformation capture techniques and fluorescence microscopy approaches to visualize the organization of DNA elements in the nucleus. Here we will summarize relevant aspects of how the interplay between CRE activity and chromatin 3D organization regulates developmental gene expression and how it relates to pathological conditions and the evolution of animal body plan.

Super-enhancers include classical enhancers and facilitators to fully activate gene expression

Super-enhancers are compound regulatory elements that control expression of key cell identity genes. They recruit high levels of tissue-specific transcription factors and co-activators such as the Mediator complex and contact target gene promoters with high frequency. Most super-enhancers contain multiple constituent regulatory elements, but it is unclear whether these elements have distinct roles in activating target gene expression. Here, by rebuilding the endogenous multipartite α-globin super-enhancer, we show that it contains bioinformatically equivalent but functionally distinct element types: classical enhancers and facilitator elements. Facilitators have no intrinsic enhancer activity, yet in their absence, classical enhancers are unable to fully upregulate their target genes. Without facilitators, classical enhancers exhibit reduced Mediator recruitment, enhancer RNA transcription, and enhancer-promoter interactions. Facilitators are interchangeable but display functional hierarchy based on their position within a multipartite enhancer. Facilitators thus play an important role in potentiating the activity of classical enhancers and ensuring robust activation of target genes.

Targeting the epigenome to reinvigorate T cells for cancer immunotherapy

Cancer immunotherapy using immune-checkpoint inhibitors (ICIs) has revolutionized the field of cancer treatment; however, ICI efficacy is constrained by progressive dysfunction of CD8+ tumor-infiltrating lymphocytes (TILs), which is termed T cell exhaustion. This process is driven by diverse extrinsic factors across heterogeneous tumor immune microenvironment (TIME). Simultaneously, tumorigenesis entails robust reshaping of the epigenetic landscape, potentially instigating T cell exhaustion. In this review, we summarize the epigenetic mechanisms governing tumor microenvironmental cues leading to T cell exhaustion, and discuss therapeutic potential of targeting epigenetic regulators for immunotherapies. Finally, we outline conceptual and technical advances in developing potential treatment paradigms involving immunostimulatory agents and epigenetic therapies.

Intra- and interchromosomal contact mapping reveals the Igh locus has extensive conformational heterogeneity and interacts with B-lineage genes

To produce a diverse antibody repertoire, immunoglobulin heavy-chain (Igh) loci undergo large-scale alterations in structure to facilitate juxtaposition and recombination of spatially separated variable (VH), diversity (DH), and joining (JH) genes. These chromosomal alterations are poorly understood. Uncovering their patterns shows how chromosome dynamics underpins antibody diversity. Using tiled Capture Hi-C, we produce a comprehensive map of chromatin interactions throughout the 2.8-Mb Igh locus in progenitor B cells. We find that the Igh locus folds into semi-rigid subdomains and undergoes flexible looping of the VH genes to its 3' end, reconciling two views of locus organization. Deconvolution of single Igh locus conformations using polymer simulations identifies thousands of different structures. This heterogeneity may underpin the diversity of V(D)J recombination events. All three immunoglobulin loci also participate in a highly specific, developmentally regulated network of interchromosomal interactions with genes encoding B cell-lineage factors. This suggests a model of interchromosomal coordination of B cell development.

Three-dimensional genome structure and function

Linear DNA undergoes a series of compression and folding events, forming various three-dimensional (3D) structural units in mammalian cells, including chromosomal territory, compartment, topologically associating domain, and chromatin loop. These structures play crucial roles in regulating gene expression, cell differentiation, and disease progression. Deciphering the principles underlying 3D genome folding and the molecular mechanisms governing cell fate determination remains a challenge. With advancements in high-throughput sequencing and imaging techniques, the hierarchical organization and functional roles of higher-order chromatin structures have been gradually illuminated. This review systematically discussed the structural hierarchy of the 3D genome, the effects and mechanisms of cis-regulatory elements interaction in the 3D genome for regulating spatiotemporally specific gene expression, the roles and mechanisms of dynamic changes in 3D chromatin conformation during embryonic development, and the pathological mechanisms of diseases such as congenital developmental abnormalities and cancer, which are attributed to alterations in 3D genome organization and aberrations in key structural proteins. Finally, prospects were made for the research about 3D genome structure, function, and genetic intervention, and the roles in disease development, prevention, and treatment, which may offer some clues for precise diagnosis and treatment of related diseases.

Multidimensional profiling reveals GATA1-modulated stage-specific chromatin states and functional associations during human erythropoiesis

Mammalian erythroid development can be divided into three stages: hematopoietic stem and progenitor cell (HSPC), erythroid progenitor (Ery-Pro), and erythroid precursor (Ery-Pre). However, the mechanisms by which the 3D genome changes to establish the stage-specific transcription programs that are critical for erythropoiesis remain unclear. Here, we analyze the chromatin landscape at multiple levels in defined populations from primary human erythroid culture. While compartments and topologically associating domains remain largely unchanged, ∼50% of H3K27Ac-marked enhancers are dynamic in HSPC versus Ery-Pre. The enhancer anchors of enhancer-promoter loops are enriched for occupancy of respective stage-specific transcription factors (TFs), indicating these TFs orchestrate the enhancer connectome rewiring. The master TF of erythropoiesis, GATA1, is found to occupy most erythroid gene promoters at the Ery-Pro stage, and mediate conspicuous local rewiring through acquiring binding at the distal regions in Ery-Pre, promoting productive erythroid transcription output. Knocking out GATA1 binding sites precisely abrogates local rewiring and corresponding gene expression. Interestingly, knocking down GATA1 can transiently revert the cell state to an earlier stage and prolong the window of progenitor state. This study reveals mechanistic insights underlying chromatin rearrangements during development by integrating multidimensional chromatin landscape analyses to associate with transcription output and cellular states.

The Mediator complex regulates enhancer-promoter interactions

Enhancer-mediated gene activation generally requires physical proximity between enhancers and their target gene promoters. However, the molecular mechanisms by which interactions between enhancers and promoters are formed are not well understood. Here, we investigate the function of the Mediator complex in the regulation of enhancer-promoter interactions, by combining rapid protein depletion and high-resolution MNase-based chromosome conformation capture approaches. We show that depletion of Mediator leads to reduced enhancer-promoter interaction frequencies, which are associated with a strong decrease in gene expression. In addition, we find increased interactions between CTCF-binding sites upon Mediator depletion. These changes in chromatin architecture are associated with a redistribution of the Cohesin complex on chromatin and a reduction in Cohesin occupancy at enhancers. Together, our results indicate that the Mediator and Cohesin complexes contribute to enhancer-promoter interactions and provide insights into the molecular mechanisms by which communication between enhancers and promoters is regulated.

Linking genome structures to functions by simultaneous single-cell Hi-C and RNA-seq

Much progress has been made recently in single-cell chromosome conformation capture technologies. However, a method that allows simultaneous profiling of chromatin architecture and gene expression has not been reported. Here, we developed an assay named "Hi-C and RNA-seq employed simultaneously" (HiRES) and performed it on thousands of single cells from developing mouse embryos. Single-cell three-dimensional genome structures, despite being heavily determined by the cell cycle and developmental stages, gradually diverged in a cell type-specific manner as development progressed. By comparing the pseudotemporal dynamics of chromatin interactions with gene expression, we found a widespread chromatin rewiring that occurred before transcription activation. Our results demonstrate that the establishment of specific chromatin interactions is tightly related to transcriptional control and cell functions during lineage specification.

Region Capture Micro-C reveals coalescence of enhancers and promoters into nested microcompartments

Although enhancers are central regulators of mammalian gene expression, the mechanisms underlying enhancer-promoter (E-P) interactions remain unclear. Chromosome conformation capture (3C) methods effectively capture large-scale three-dimensional (3D) genome structure but struggle to achieve the depth necessary to resolve fine-scale E-P interactions. Here, we develop Region Capture Micro-C (RCMC) by combining micrococcal nuclease (MNase)-based 3C with a tiling region-capture approach and generate the deepest 3D genome maps reported with only modest sequencing. By applying RCMC in mouse embryonic stem cells and reaching the genome-wide equivalent of ~317 billion unique contacts, RCMC reveals previously unresolvable patterns of highly nested and focal 3D interactions, which we term microcompartments. Microcompartments frequently connect enhancers and promoters, and although loss of loop extrusion and inhibition of transcription disrupts some microcompartments, most are largely unaffected. We therefore propose that many E-P interactions form through a compartmentalization mechanism, which may partially explain why acute cohesin depletion only modestly affects global gene expression.

Determining chromatin architecture with Micro Capture-C

Micro Capture-C (MCC) is a chromatin conformation capture (3C) method for visualizing reproducible three-dimensional contacts of specified regions of the genome at base pair resolution. These methods are an established family of techniques that use proximity ligation to assay the topology of chromatin. MCC can generate data at substantially higher resolution than previous techniques through multiple refinements of the 3C method. Using a sequence agnostic nuclease, the maintenance of cellular integrity and full sequencing of the ligation junctions, MCC achieves subnucleosomal levels of resolution, which can be used to reveal transcription factor binding sites analogous to DNAse I footprinting. Gene dense regions, close-range enhancer-promoter contacts, individual enhancers within super-enhancers and multiple other types of loci or regulatory regions that were previously challenging to assay with conventional 3C techniques, are readily observed using MCC. MCC requires training in common molecular biology techniques and bioinformatics to perform the experiment and analyze the data. The protocol can be expected to be completed in a 3 week timeframe for experienced molecular biologists.

Live imaging reveals chromatin compaction transitions and dynamic transcriptional bursting during stem cell differentiation in vivo

Stem cell differentiation requires dramatic changes in gene expression and global remodeling of chromatin architecture. How and when chromatin remodels relative to the transcriptional, behavioral, and morphological changes during differentiation remain unclear, particularly in an intact tissue context. Here, we develop a quantitative pipeline which leverages fluorescently-tagged histones and longitudinal imaging to track large-scale chromatin compaction changes within individual cells in a live mouse. Applying this pipeline to epidermal stem cells, we reveal that cell-to-cell chromatin compaction heterogeneity within the stem cell compartment emerges independent of cell cycle status, and instead is reflective of differentiation status. Chromatin compaction state gradually transitions over days as differentiating cells exit the stem cell compartment. Moreover, establishing live imaging of Keratin-10 (K10) nascent RNA, which marks the onset of stem cell differentiation, we find that Keratin-10 transcription is highly dynamic and largely precedes the global chromatin compaction changes associated with differentiation. Together, these analyses reveal that stem cell differentiation involves dynamic transcriptional states and gradual chromatin rearrangement.

Gene silencing dynamics are modulated by transiently active regulatory elements

Transcriptional enhancers have been extensively characterized, but cis-regulatory elements involved in acute gene repression have received less attention. Transcription factor GATA1 promotes erythroid differentiation by activating and repressing distinct gene sets. Here, we study the mechanism by which GATA1 silences the proliferative gene Kit during murine erythroid cell maturation and define stages from initial loss of activation to heterochromatinization. We find that GATA1 inactivates a potent upstream enhancer but concomitantly creates a discrete intronic regulatory region marked by H3K27ac, short noncoding RNAs, and de novo chromatin looping. This enhancer-like element forms transiently and serves to delay Kit silencing. The element is ultimately erased via the FOG1/NuRD deacetylase complex, as revealed by the study of a disease-associated GATA1 variant. Hence, regulatory sites can be self-limiting by dynamic co-factor usage. Genome-wide analyses across cell types and species uncover transiently active elements at numerous genes during repression, suggesting that modulation of silencing kinetics is widespread.

Chromosome conformation capture approaches to investigate 3D genome architecture in Ankylosing Spondylitis

Ankylosing Spondylitis (AS) is a chronic inflammatory arthritis of the spine exhibiting a strong genetic background. The mechanistic and functional understanding of the AS-associated genomic loci, identified with Genome Wide Association Studies (GWAS), remains challenging. Chromosome conformation capture (3C) and derivatives are recent techniques which are of great help in elucidating the spatial genome organization and of enormous support in uncover a mechanistic explanation for disease-associated genetic variants. The perturbation of three-dimensional (3D) genome hierarchy may lead to a plethora of human diseases, including rheumatological disorders. Here we illustrate the latest approaches and related findings on the field of genome organization, highlighting how the instability of 3D genome conformation may be among the causes of rheumatological disease phenotypes. We suggest a new perspective on the inclusive potential of a 3C approach to inform GWAS results in rheumatic diseases. 3D genome organization may ultimately lead to a more precise and comprehensive functional interpretation of AS association, which is the starting point for emerging and more specific therapies.

Low input capture Hi-C (liCHi-C) identifies promoter-enhancer interactions at high-resolution

Long-range interactions between regulatory elements and promoters are key in gene transcriptional control; however, their study requires large amounts of starting material, which is not compatible with clinical scenarios nor the study of rare cell populations. Here we introduce low input capture Hi-C (liCHi-C) as a cost-effective, flexible method to map and robustly compare promoter interactomes at high resolution. As proof of its broad applicability, we implement liCHi-C to study normal and malignant human hematopoietic hierarchy in clinical samples. We demonstrate that the dynamic promoter architecture identifies developmental trajectories and orchestrates transcriptional transitions during cell-state commitment. Moreover, liCHi-C enables the identification of disease-relevant cell types, genes and pathways potentially deregulated by non-coding alterations at distal regulatory elements. Finally, we show that liCHi-C can be harnessed to uncover genome-wide structural variants, resolve their breakpoints and infer their pathogenic effects. Collectively, our optimized liCHi-C method expands the study of 3D chromatin organization to unique, low-abundance cell populations, and offers an opportunity to uncover factors and regulatory networks involved in disease pathogenesis.

On-microscope staging of live cells reveals changes in the dynamics of transcriptional bursting during differentiation

Determining the mechanisms by which genes are switched on and off during development is a key aim of current biomedical research. Gene transcription has been widely observed to occur in a discontinuous fashion, with short bursts of activity interspersed with periods of inactivity. It is currently not known if or how this dynamic behaviour changes as mammalian cells differentiate. To investigate this, using an on-microscope analysis, we monitored mouse α-globin transcription in live cells throughout erythropoiesis. We find that changes in the overall levels of α-globin transcription are most closely associated with changes in the fraction of time a gene spends in the active transcriptional state. We identify differences in the patterns of transcriptional bursting throughout differentiation, with maximal transcriptional activity occurring in the mid-phase of differentiation. Early in differentiation, we observe increased fluctuation in transcriptional activity whereas at the peak of gene expression, in early erythroblasts, transcription is relatively stable. Later during differentiation as α-globin expression declines, we again observe more variability in transcription within individual cells. We propose that the observed changes in transcriptional behaviour may reflect changes in the stability of active transcriptional compartments as gene expression is regulated during differentiation.

Natural and Experimental Rewiring of Gene Regulatory Regions

The successful development and ongoing functioning of complex organisms depend on the faithful execution of the genetic code. A critical step in this process is the correct spatial and temporal expression of genes. The highly orchestrated transcription of genes is controlled primarily by cis-regulatory elements: promoters, enhancers, and insulators. The medical importance of this key biological process can be seen by the frequency with which mutations and inherited variants that alter cis-regulatory elements lead to monogenic and complex diseases and cancer. Here, we provide an overview of the methods available to characterize and perturb gene regulatory circuits. We then highlight mechanisms through which regulatory rewiring contributes to disease, and conclude with a perspective on how our understanding of gene regulation can be used to improve human health.

The 3D genome landscape: Diverse chromosomal interactions and their functional implications

Genome organization includes contacts both within a single chromosome and between distinct chromosomes. Thus, regulatory organization in the nucleus may include interplay of these two types of chromosomal interactions with genome activity. Emerging advances in omics and single-cell imaging technologies have allowed new insights into chromosomal contacts, including those of homologs and sister chromatids, and their significance to genome function. In this review, we highlight recent studies in this field and discuss their impact on understanding the principles of chromosome organization and associated functional implications in diverse cellular processes. Specifically, we describe the contributions of intra-chromosomal, inter-homolog, and inter-sister chromatid contacts to genome organization and gene expression.

The role of specialized cell cycles during erythroid lineage development: insights from single-cell RNA sequencing

Early erythroid progenitors known as CFU-e undergo multiple self-renewal cell cycles. The CFU-e developmental stage ends with the onset of erythroid terminal differentiation (ETD). The transition from CFU-e to ETD is a critical cell fate decision that determines erythropoietic rate. Here we review recent insights into the regulation of this transition, garnered from flow cytometric and single-cell RNA sequencing studies. We find that the CFU-e/ETD transition is a rapid S phase-dependent transcriptional switch. It takes place during an S phase that is much shorter than in preceding or subsequent cycles, as a result of globally faster replication forks. Furthermore, it is preceded by cycles in which G1 becomes gradually shorter. These dramatic cell cycle and S phase remodeling events are directly linked to regulation of the CFU-e/ETD switch. Moreover, regulators of erythropoietic rate exert their effects by modulating cell cycle duration and S phase speed. Glucocorticoids increase erythropoietic rate by inducing the CDK inhibitor p57^KIP2, which slows replication forks, inhibiting the CFU-e/ETD switch. Conversely, erythropoietin promotes induction of ETD by shortening the cycle. S phase shortening was reported during cell fate decisions in non-erythroid lineages, suggesting a fundamentally new developmental role for cell cycle speed.

Promoter and enhancer RNAs regulate chromatin reorganization and activation of miR-10b/HOXD locus, and neoplastic transformation in glioma

miR-10b is silenced in normal neuroglial cells of the brain but commonly activated in glioma, where it assumes an essential tumor-promoting role. We demonstrate that the entire miR-10b-hosting HOXD locus is activated in glioma via the cis-acting mechanism involving 3D chromatin reorganization and CTCF-cohesin-mediated looping. This mechanism requires two interacting lncRNAs, HOXD-AS2 and LINC01116, one associated with HOXD3/HOXD4/miR-10b promoter and another with the remote enhancer. Knockdown of either lncRNA in glioma cells alters CTCF and cohesin binding, abolishes chromatin looping, inhibits the expression of all genes within HOXD locus, and leads to glioma cell death. Conversely, in cortical astrocytes, enhancer activation is sufficient for HOXD/miR-10b locus reorganization, gene derepression, and neoplastic cell transformation. LINC01116 RNA is essential for this process. Our results demonstrate the interplay of two lncRNAs in the chromatin folding and concordant regulation of miR-10b and multiple HOXD genes normally silenced in astrocytes and triggering the neoplastic glial transformation.

RASER-FISH: non-denaturing fluorescence in situ hybridization for preservation of three-dimensional interphase chromatin structure

DNA fluorescence in situ hybridization (FISH) has been a central technique in advancing our understanding of how chromatin is organized within the nucleus. With the increasing resolution offered by super-resolution microscopy, the optimal maintenance of chromatin structure within the nucleus is essential for accuracy in measurements and interpretation of data. However, standard 3D-FISH requires potentially destructive heat denaturation in the presence of chaotropic agents such as formamide to allow access to the DNA strands for labeled FISH probes. To avoid the need to heat-denature, we developed Resolution After Single-strand Exonuclease Resection (RASER)-FISH, which uses exonuclease digestion to generate single-stranded target DNA for efficient probe binding over a 2 d process. Furthermore, RASER-FISH is easily combined with immunostaining of nuclear proteins or the detection of RNAs. Here, we provide detailed procedures for RASER-FISH in mammalian cultured cells to detect single loci, chromatin tracks and topologically associating domains with conventional and super-resolution 3D structured illumination microscopy. Moreover, we provide a validation and characterization of our method, demonstrating excellent preservation of chromatin structure and nuclear integrity, together with improved hybridization efficiency, compared with classic 3D-FISH protocols.

Analysis of sub-kilobase chromatin topology reveals nano-scale regulatory interactions with variable dependence on cohesin and CTCF

Enhancers and promoters predominantly interact within large-scale topologically associating domains (TADs), which are formed by loop extrusion mediated by cohesin and CTCF. However, it is unclear whether complex chromatin structures exist at sub-kilobase-scale and to what extent fine-scale regulatory interactions depend on loop extrusion. To address these questions, we present an MNase-based chromosome conformation capture (3C) approach, which has enabled us to generate the most detailed local interaction data to date (20 bp resolution) and precisely investigate the effects of cohesin and CTCF depletion on chromatin architecture. Our data reveal that cis-regulatory elements have distinct internal nano-scale structures, within which local insulation is dependent on CTCF, but which are independent of cohesin. In contrast, we find that depletion of cohesin causes a subtle reduction in longer-range enhancer-promoter interactions and that CTCF depletion can cause rewiring of regulatory contacts. Together, our data show that loop extrusion is not essential for enhancer-promoter interactions, but contributes to their robustness and specificity and to precise regulation of gene expression.

Making connections: enhancers in cellular differentiation

Deciphering the process by which hundreds of distinct cell types emerge from a single zygote to form a complex multicellular organism remains one of the greatest challenges in biological research. Enhancers are known to be central to cell type-specific gene expression, yet many questions regarding how these genomic elements interact both temporally and spatially with other cis- and trans-acting factors to control transcriptional activity during differentiation and development remain unanswered. Here, we review our current understanding of the role of enhancers and their interactions in this context and highlight recent progress achieved with experimental methods of unprecedented resolution.

Epigenetic and Transcriptional Control of Erythropoiesis

Erythropoiesis is a process of enormous magnitude, with the average person generating two to three million red cells every second. Erythroid progenitors start as large cells with large nuclei, and over the course of three to four cell divisions they undergo a dramatic decrease in cell size accompanied by profound nuclear condensation, which culminates in enucleation. As maturing erythroblasts are undergoing these dramatic phenotypic changes, they accumulate hemoglobin and express high levels of other erythroid-specific genes, while silencing much of the non-erythroid transcriptome. These phenotypic and gene expression changes are associated with distinct changes in the chromatin landscape, and require close coordination between transcription factors and epigenetic regulators, as well as precise regulation of RNA polymerase II activity. Disruption of these processes are associated with inherited anemias and myelodysplastic syndromes. Here, we review the epigenetic mechanisms that govern terminal erythroid maturation, and their role in human disease.

HiCAR is a robust and sensitive method to analyze open-chromatin-associated genome organization

The long-range interactions of cis-regulatory elements (cREs) play a central role in gene regulation. cREs can be characterized as accessible chromatin sequences. However, it remains technically challenging to comprehensively identify their spatial interactions. Here, we report a new method HiCAR (Hi-C on accessible regulatory DNA), which utilizes Tn5 transposase and chromatin proximity ligation, for the analysis of open-chromatin-anchored interactions with low-input cells. By applying HiCAR in human embryonic stem cells and lymphoblastoid cells, we demonstrate that HiCAR identifies high-resolution chromatin contacts with an efficiency comparable with that of in situ Hi-C over all distance ranges. Interestingly, we found that the "poised" gene promoters exhibit silencer-like function to repress the expression of distal genes via promoter-promoter interactions. Lastly, we applied HiCAR to 30,000 primary human muscle stem cells and demonstrated that HiCAR is capable of analyzing chromatin accessibility and looping using low-input primary cells and clinical samples.

Dynamic Runx1 chromatin boundaries affect gene expression in hematopoietic development

The transcription factor RUNX1 is a critical regulator of developmental hematopoiesis and is frequently disrupted in leukemia. Runx1 is a large, complex gene that is expressed from two alternative promoters under the spatiotemporal control of multiple hematopoietic enhancers. To dissect the dynamic regulation of Runx1 in hematopoietic development, we analyzed its three-dimensional chromatin conformation in mouse embryonic stem cell (ESC) differentiation cultures. Runx1 resides in a 1.1 Mb topologically associating domain (TAD) demarcated by convergent CTCF motifs. As ESCs differentiate to mesoderm, chromatin accessibility, Runx1 enhancer-promoter (E-P) interactions, and CTCF-CTCF interactions increase in the TAD, along with initiation of Runx1 expression from the P2 promoter. Differentiation to hematopoietic progenitor cells is associated with the formation of tissue-specific sub-TADs over Runx1, a shift in E-P interactions, P1 promoter demethylation, and robust expression from both Runx1 promoters. Deletion of promoter-proximal CTCF sites at the sub-TAD boundaries has no obvious effects on E-P interactions but leads to partial loss of domain structure, mildly affects gene expression, and delays hematopoietic development. Together, our analysis of gene regulation at a large multi-promoter developmental gene reveals that dynamic sub-TAD chromatin boundaries play a role in establishing TAD structure and coordinated gene expression.

Capture-C: a modular and flexible approach for high-resolution chromosome conformation capture

Chromosome conformation capture (3C) methods measure the spatial proximity between DNA elements in the cell nucleus. Many methods have been developed to sample 3C material, including the Capture-C family of protocols. Capture-C methods use oligonucleotides to enrich for interactions of interest from sequencing-ready 3C libraries. This approach is modular and has been adapted and optimized to work for sampling of disperse DNA elements (NuTi Capture-C), including from low cell inputs (LI Capture-C), as well as to generate Hi-C like maps for specific regions of interest (Tiled-C) and to interrogate multiway interactions (Tri-C). We present the design, experimental protocol and analysis pipeline for NuTi Capture-C in addition to the variations for generation of LI Capture-C, Tiled-C and Tri-C data. The entire procedure can be performed in 3 weeks and requires standard molecular biology skills and equipment, access to a next-generation sequencing platform, and basic bioinformatic skills. Implemented with other sequencing technologies, these methods can be used to identify regulatory interactions and to compare the structural organization of the genome in different cell types and genetic models.

The stochastic nature of genome organization and function

Genomes have complex three-dimensional structures. High-resolution population-based biochemical studies over the last decade have painted a mostly static picture of the genome characterized by universal organizational features, such as chromatin domains and compartments. Yet, when analyzed at the single cell level, these architectural elements are highly variable. The heterogeneity in genome organization is in line with the inherent stochasticity of transcription that shows high variation between individual cells. We highlight recent findings on single-cell variability in genome organization and describe a framework for how the stochastic nature of chromatin organization may relate to transcription dynamics.

Contribution of 3D genome topological domains to genetic risk of cancers: a genome-wide computational study

Background

Genome-wide association studies have identified statistical associations between various diseases, including cancers, and a large number of single-nucleotide polymorphisms (SNPs). However, they provide no direct explanation of the mechanisms underlying the association. Based on the recent discovery that changes in three-dimensional genome organization may have functional consequences on gene regulation favoring diseases, we investigated systematically the genome-wide distribution of disease-associated SNPs with respect to a specific feature of 3D genome organization: topologically associating domains (TADs) and their borders.

Results

For each of 449 diseases, we tested whether the associated SNPs are present in TAD borders more often than observed by chance, where chance (i.e., the null model in statistical terms) corresponds to the same number of pointwise loci drawn at random either in the entire genome, or in the entire set of disease-associated SNPs listed in the GWAS catalog. Our analysis shows that a fraction of diseases displays such a preferential localization of their risk loci. Moreover, cancers are relatively more frequent among these diseases, and this predominance is generally enhanced when considering only intergenic SNPs. The structure of SNP-based diseasome networks confirms that localization of risk loci in TAD borders differs between cancers and non-cancer diseases. Furthermore, different TAD border enrichments are observed in embryonic stem cells and differentiated cells, consistent with changes in topological domains along embryogenesis and delineating their contribution to disease risk.

Conclusions

Our results suggest that, for certain diseases, part of the genetic risk lies in a local genetic variation affecting the genome partitioning in topologically insulated domains. Investigating this possible contribution to genetic risk is particularly relevant in cancers. This study thus opens a way of interpreting genome-wide association studies, by distinguishing two types of disease-associated SNPs: one with an effect on an individual gene, the other acting in interplay with 3D genome organization.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40246-022-00375-2.

Assessment of Multiway Interactions with Tri-C

Tri-C is a chromosome conformation capture (3C) approach that can efficiently identify multiway chromatin interactions with viewpoints of interest. As opposed to pair-wise interactions identified in methods such as Hi-C, 4C, and Capture-C, the detection of multiway interactions allows researchers to investigate how multiple cis-regulatory elements interact together in higher-order structures in single nuclei and address questions regarding structural cooperation between these elements. Here, we describe the procedure for designing and performing a Tri-C experiment.

Discrete regulatory modules instruct hematopoietic lineage commitment and differentiation

Lineage commitment and differentiation is driven by the concerted action of master transcriptional regulators at their target chromatin sites. Multiple efforts have characterized the key transcription factors (TFs) that determine the various hematopoietic lineages. However, the temporal interactions between individual TFs and their chromatin targets during differentiation and how these interactions dictate lineage commitment remains poorly understood. Here we perform dense, daily, temporal profiling of chromatin accessibility (DNase I-seq) and gene expression changes (total RNA-seq) along ex vivo human erythropoiesis to comprehensively define developmentally regulated DNase I hypersensitive sites (DHSs) and transcripts. We link both distal DHSs to their target gene promoters and individual TFs to their target DHSs, revealing that the regulatory landscape is organized in distinct sequential regulatory modules that regulate lineage restriction and maturation. Finally, direct comparison of transcriptional dynamics (bulk and single-cell) and lineage potential between erythropoiesis and megakaryopoiesis uncovers differential fate commitment dynamics between the two lineages as they exit the stem and progenitor stage. Collectively, these data provide insights into the temporally regulated synergy of the cis- and the trans-regulatory components underlying hematopoietic lineage commitment and differentiation.