The three-dimensional genome: regulating gene expression during pluripotency and development

BINDER achieves accurate identification of hierarchical TADs by comprehensively characterizing consensus TAD boundaries

As crucial chromatin structures, hierarchical TADs play important roles in epigenetic organization, transcriptional activity, gene regulation, and cell differentiation. Currently, it remains a highly challenging task to accurately identify hierarchical TADs in a computational manner. The key bottleneck for existing TAD callers lies in the difficulty in the prediction of precise TAD boundaries. We solve this problem by introducing a novel algorithm, called BINDER, which conducts a boundary consensus approach, and then precisely locate hierarchical TAD boundaries by developing a multifaceted boundary characterization strategy. In comparison with other leading TAD callers, BINDER shows significant improvement in identifying hierarchical TADs and exhibits the strongest robustness with ultrasparse data, which supports the importance of boundary identification in calling hierarchical TADs. Applying BINDER to experimental data and mouse hematopoietic cases, we find that the hierarchical TADs identified by BINDER show strong biological relevance in their epigenetic organization, transcriptional activity, DNA motifs, and coregulation during cellular differentiation. BINDER discovers differences in the enrichment of two specific transcription factors, CHD1 and CHD2, at TAD boundaries with different hierarchies. It also observes variations in the gene expression of TADs with different hierarchies during cellular differentiation.

Genome structural dynamics: insights from Gaussian network analysis of Hi-C data

Characterization of the spatiotemporal properties of the chromatin is essential to gaining insights into the physical bases of gene co-expression, transcriptional regulation and epigenetic modifications. The Gaussian network model (GNM) has proven in recent work to serve as a useful tool for modeling chromatin structural dynamics, using as input high-throughput chromosome conformation capture data. We focus here on the exploration of the collective dynamics of chromosomal structures at hierarchical levels of resolution, from single gene loci to topologically associating domains or entire chromosomes. The GNM permits us to identify long-range interactions between gene loci, shedding light on the role of cross-correlations between distal regions of the chromosomes in regulating gene expression. Notably, GNM analysis performed across diverse cell lines highlights the conservation of the global/cooperative movements of the chromatin across different types of cells. Variations driven by localized couplings between genomic loci, on the other hand, underlie cell differentiation, underscoring the significance of the four-dimensional properties of the genome in defining cellular identity. Finally, we demonstrate the close relation between the cell type-dependent mobility profiles of gene loci and their gene expression patterns, providing a clear demonstration of the role of chromosomal 4D features in defining cell-specific differential expression of genes.

RAD21 deficiency drives corneal to scleral differentiation fate switching via upregulating WNT9B

The cornea and sclera are distinct adjacent tissues, yet their stromal cells originate from common neural crest cells (NCCs). Sclerocornea is a disease characterized by an indistinguishable boundary between the cornea and sclera. Previously, we identified a RAD21 mutation in a sclerocornea pedigree. Here, we investigated the impacts of RAD21 on NCC activities during eye development. RAD21 deficiency caused upregulation of PCDHGC3. Both RAD21 knockdown and PCDHGC3 upregulation disrupted the migration of NCCs. Transcriptome analysis indicated that WNT9B had 190.9-fold higher expression in scleral stroma than in corneal stroma. WNT9B was also significantly upregulated by both RAD21 knockdown and PCDHGC3 overexpression, and knock down of WNT9B rescued the differentiation and migration of NCCs with RAD21 deficiency. Consistently, overexpressing wnt9b in Xenopus tropicalis led to ocular developmental abnormalities. In summary, WNT9B is a determinant factor during NCC differentiation into corneal keratocytes or scleral stromal cells and is affected by RAD21 expression.

Multiplexed chromatin imaging reveals predominantly pairwise long-range coordination between Drosophila Polycomb genes

Polycomb (Pc) group proteins are transcriptional regulators with key roles in development, cell identity, and differentiation. Pc-bound chromatin regions form repressive domains that interact in 3D to assemble repressive nuclear compartments. Here, we use multiplexed chromatin imaging to investigate whether Pc compartments involve the clustering of multiple Pc domains during Drosophila development. Notably, 3D proximity between Pc targets is rare and involves predominantly pairwise interactions. These 3D proximities are particularly enhanced in segments where Pc genes are co-repressed. In addition, segment-specific expression of Hox Pc targets leads to their spatial segregation from Pc-repressed genes. Finally, non-Hox Pc targets are more proximal in regions where they are co-expressed. These results indicate that long-range Pc interactions are temporally and spatially regulated during differentiation and development but do not induce frequent clustering of multiple distant Pc genes.

Analysis of long-range chromatin contacts, compartments and looping between mouse embryonic stem cells, lens epithelium and lens fibers

BACKGROUND

Nuclear organization of interphase chromosomes involves individual chromosome territories, "open" and "closed" chromatin compartments, topologically associated domains (TADs) and chromatin loops. The DNA- and RNA-binding transcription factor CTCF together with the cohesin complex serve as major organizers of chromatin architecture. Cellular differentiation is driven by temporally and spatially coordinated gene expression that requires chromatin changes of individual loci of various complexities. Lens differentiation represents an advantageous system to probe transcriptional mechanisms underlying tissue-specific gene expression including high transcriptional outputs of individual crystallin genes until the mature lens fiber cells degrade their nuclei.

RESULTS

Chromatin organization between mouse embryonic stem (ES) cells, newborn (P0.5) lens epithelium and fiber cells were analyzed using Hi-C. Localization of CTCF in both lens chromatins was determined by ChIP-seq and compared with ES cells. Quantitative analyses show major differences between number and size of TADs and chromatin loop size between these three cell types. In depth analyses show similarities between lens samples exemplified by overlaps between compartments A and B. Lens epithelium-specific CTCF peaks are found in mostly methylated genomic regions while lens fiber-specific and shared peaks occur mostly within unmethylated DNA regions. Major differences in TADs and loops are illustrated at the ~ 500 kb Pax6 locus, encoding the critical lens regulatory transcription factor and within a larger ~ 15 Mb WAGR locus, containing Pax6 and other loci linked to human congenital diseases. Lens and ES cell Hi-C data (TADs and loops) together with ATAC-seq, CTCF, H3K27ac, H3K27me3 and ENCODE cis-regulatory sites are shown in detail for the Pax6, Sox1 and Hif1a loci, multiple crystallin genes and other important loci required for lens morphogenesis. The majority of crystallin loci are marked by unexpectedly high CTCF-binding across their transcribed regions.

CONCLUSIONS

Our study has generated the first data on 3-dimensional (3D) nuclear organization in lens epithelium and lens fibers and directly compared these data with ES cells. These findings generate novel insights into lens-specific transcriptional gene control, open new research avenues to study transcriptional condensates in lens fiber cells, and enable studies of non-coding genetic variants linked to cataract and other lens and ocular abnormalities.

Noncoding mutations cause super-enhancer retargeting resulting in protein synthesis dysregulation during B cell lymphoma progression

Whole-genome sequencing of longitudinal tumor pairs representing transformation of follicular lymphoma to high-grade B cell lymphoma with MYC and BCL2 rearrangements (double-hit lymphoma) identified coding and noncoding genomic alterations acquired during lymphoma progression. Many of these transformation-associated alterations recurrently and focally occur at topologically associating domain resident regulatory DNA elements, including H3K4me3 promoter marks located within H3K27ac super-enhancer clusters in B cell non-Hodgkin lymphoma. One region found to undergo recurrent alteration upon transformation overlaps a super-enhancer affecting the expression of the PAX5/ZCCHC7 gene pair. ZCCHC7 encodes a subunit of the Trf4/5-Air1/2-Mtr4 polyadenylation-like complex and demonstrated copy number gain, chromosomal translocation and enhancer retargeting-mediated transcriptional upregulation upon lymphoma transformation. Consequently, lymphoma cells demonstrate nucleolar dysregulation via altered noncoding 5.8S ribosomal RNA processing. We find that a noncoding mutation acquired during lymphoma progression affects noncoding rRNA processing, thereby rewiring protein synthesis leading to oncogenic changes in the lymphoma proteome.

PRR14 organizes H3K9me3-modified heterochromatin at the nuclear lamina

The eukaryotic genome is organized in three dimensions within the nucleus. Transcriptionally active chromatin is spatially separated from silent heterochromatin, a large fraction of which is located at the nuclear periphery. However, the mechanisms by which chromatin is localized at the nuclear periphery remain poorly understood. Here we demonstrate that Proline Rich 14 (PRR14) protein organizes H3K9me3-modified heterochromatin at the nuclear lamina. We show that PRR14 dynamically associates with both the nuclear lamina and heterochromatin, and is able to reorganize heterochromatin in the nucleus of interphase cells independent of mitosis. We characterize two functional HP1-binding sites within PRR14 that contribute to its association with heterochromatin. We also demonstrate that PPR14 forms an anchoring surface for heterochromatin at the nuclear lamina where it interacts dynamically with HP1-associated chromatin. Our study proposes a model of dynamic heterochromatin organization at the nuclear lamina via the PRR14 tethering protein.

The secret life of chromatin tethers

The nuclear envelope plays an essential role in organizing the genome inside of the nucleus. The inner nuclear membrane is coated with a meshwork of filamentous lamin proteins that provide a surface to organize a variety of cellular processes. A subset of nuclear lamina- and membrane-associated proteins functions as anchors to hold transcriptionally silent heterochromatin at the nuclear periphery. While most chromatin tethers are integral membrane proteins, a limited number are lamina-bound. One example is the mammalian proline-rich 14 (PRR14) protein. PRR14 is a recently characterized protein with unique function that is different from other known chromatin tethers. Here, we review our current understanding of PRR14 structure and function in organizing heterochromatin at the nuclear periphery.

Combined loss of CDH1 and downstream regulatory sequences drive early-onset diffuse gastric cancer and increase penetrance of hereditary diffuse gastric cancer

BACKGROUND

Germline CDH1 pathogenic or likely pathogenic variants cause hereditary diffuse gastric cancer (HDGC). Once a genetic cause is identified, stomachs' and breasts' surveillance and/or prophylactic surgery is offered to asymptomatic CDH1 carriers, which is life-saving. Herein, we characterized an inherited mechanism responsible for extremely early-onset gastric cancer and atypical HDGC high penetrance.

METHODS

Whole-exome sequencing (WES) re-analysis was performed in an unsolved HDGC family. Accessible chromatin and CDH1 promoter interactors were evaluated in normal stomach by ATAC-seq and 4C-seq, and functional analysis was performed using CRISPR-Cas9, RNA-seq and pathway analysis.

RESULTS

We identified a germline heterozygous 23 Kb CDH1-TANGO6 deletion in a family with eight diffuse gastric cancers, six before age 30. Atypical HDGC high penetrance and young cancer-onset argued towards a role for the deleted region downstream of CDH1, which we proved to present accessible chromatin, and CDH1 promoter interactors in normal stomach. CRISPR-Cas9 edited cells mimicking the CDH1-TANGO6 deletion display the strongest CDH1 mRNA downregulation, more impacted adhesion-associated, type-I interferon immune-associated and oncogenic signalling pathways, compared to wild-type or CDH1-deleted cells. This finding solved an 18-year family odyssey and engaged carrier family members in a cancer prevention pathway of care.

CONCLUSION

In this work, we demonstrated that regulatory elements lying down-stream of CDH1 are part of a chromatin network that control CDH1 expression and influence cell transcriptome and associated signalling pathways, likely explaining high disease penetrance and very young cancer-onset. This study highlights the importance of incorporating scientific-technological updates and clinical guidelines in routine diagnosis, given their impact in timely genetic diagnosis and disease prevention.

Remodeling of the 3D chromatin architecture in the marine microalga Nannochloropsis oceanica during lipid accumulation

BACKGROUND

Genomic three-dimensional (3D) spatial organization plays a key role in shaping gene expression and associated chromatin modification, and it is highly sensitive to environmental stress conditions. In microalgae, exposure to nitrogen stress can drive lipid accumulation, yet the associated functional alterations in the spatial organization of the microalgal genome have yet to be effectively characterized.

RESULTS

Accordingly, the present study employed RNA-seq, Hi-C, and ChIP-seq approaches to explore the relationship between 3D chromosomal architecture and gene expression during lipid accumulation in the marine microalga Nannochloropsis oceanica in response to nitrogen deprivation (ND). These analyses revealed that ND resulted in various changes in chromosomal organization, including A/B compartment transitions, topologically associating domain (TAD) shifts, and the disruption of short-range interactions. Significantly higher levels of gene expression were evident in A compartments and TAD boundary regions relative to B compartments and TAD interior regions, consistent with observed histone modification enrichment in these areas. ND-induced differentially expressed genes (DEGs) were notably enriched in altered TAD-associated regions and regions exhibiting differential genomic contact. These DEGs were subjected to Gene Ontology (GO) term analyses that indicated they were enriched in the 'fatty acid metabolism', 'response to stress', 'carbon fixation' and 'photosynthesis' functional categories, in line with the ND treatment conditions used to conduct this study. These data indicate that Nannochloropsis cells exhibit a clear association between chromatin organization and transcriptional activity under nitrogen stress conditions. Pronounced and extensive histone modifications were evident in response to ND. Observed changes in chromatin architecture were linked to shifts in histone modifications and gene expression.

CONCLUSIONS

Overall, the reprogramming of many lipid metabolism-associated genes was evident under nitrogen stress conditions with respect to both histone modifications and chromosomal organization. Together these results revealed that higher-order chromatin architecture represents a new layer that can guide efforts to understand the transcriptional regulation of lipid metabolism in nitrogen-deprived microalgae.

RNA polymerase II pausing regulates chromatin organization in erythrocytes

Chicken erythrocytes are nucleated cells often considered to be transcriptionally inactive, although the epigenetic changes and chromatin remodeling that would mediate transcriptional repression and the extent of gene silencing during avian terminal erythroid differentiation are not fully understood. Here, we characterize the changes in gene expression, chromatin accessibility, genome organization and chromatin nuclear disposition during the terminal stages of erythropoiesis in chicken and uncover complex chromatin reorganization at different genomic scales. We observe a robust decrease in transcription in erythrocytes, but a set of genes maintains their expression, including genes involved in RNA polymerase II (Pol II) promoter-proximal pausing. Erythrocytes exhibit a reoriented nuclear architecture, with accessible chromatin positioned towards the nuclear periphery together with the paused RNA Pol II. In erythrocytes, chromatin domains are partially lost genome-wide, except at minidomains retained around paused promoters. Our results suggest that promoter-proximal pausing of RNA Pol II contributes to the transcriptional regulation of the erythroid genome and highlight the role of RNA polymerase in the maintenance of local chromatin organization.

3D Chromatin Architecture Re-Wiring at the CDH3/CDH1 Loci Contributes to E-Cadherin to P-Cadherin Expression Switch in Gastric Cancer

Cadherins are cell-cell adhesion molecules, fundamental for cell architecture and polarity. E-cadherin to P-cadherin switch can rescue adherens junctions in epithelial tumours. Herein, we disclose a mechanism for E-cadherin to P-cadherin switch in gastric cancers. CDH1 and CDH3 mRNA expression was obtained from 42 gastric tumours' RNA-seq data. CRISPR-Cas9 was used to knock out CDH1 and a putative regulatory element. CDH1-depleted and parental cells were submitted to proteomics and enrichment GO terms analysis; ATAC-seq/4C-seq with a CDH1 promoter viewpoint to assess chromatin accessibility and conformation; and RT-PCR/flow cytometry to assess CDH1/E-cadherin and CDH3/P-cadherin expression. In 42% of gastric tumours analysed, CDH1 to CDH3 switch was observed. CDH1 knockout triggered CDH1/E-cadherin complete loss and CDH3/P-cadherin expression increase at plasma membrane. This switch, likely rescuing adherens junctions, increased cell migration/proliferation, commonly observed in aggressive tumours. E- to P-cadherin switch accompanied increased CDH1 promoter interactions with CDH3-eQTL, absent in normal stomach and parental cells. CDH3-eQTL deletion promotes CDH3/CDH1 reduced expression. These data provide evidence that loss of CDH1/E-cadherin expression alters the CDH3 locus chromatin conformation, allowing a CDH1 promoter interaction with a CDH3-eQTL, and promoting CDH3/P-cadherin expression. These data highlight a novel mechanism triggering E- to P-cadherin switch in gastric cancer.

DNA architectural protein CTCF facilitates subset-specific chromatin interactions to limit the formation of memory CD8+ T cells

Although the importance of genome organization for transcriptional regulation of cell-fate decisions and function is clear, the changes in chromatin architecture and how these impact effector and memory CD8+ T cell differentiation remain unknown. Using Hi-C, we studied how genome configuration is integrated with CD8+ T cell differentiation during infection and investigated the role of CTCF, a key chromatin remodeler, in modulating CD8+ T cell fates through CTCF knockdown approaches and perturbation of specific CTCF-binding sites. We observed subset-specific changes in chromatin organization and CTCF binding and revealed that weak-affinity CTCF binding promotes terminal differentiation of CD8+ T cells through the regulation of transcriptional programs. Further, patients with de novo CTCF mutations had reduced expression of the terminal-effector genes in peripheral blood lymphocytes. Therefore, in addition to establishing genome architecture, CTCF regulates effector CD8+ T cell heterogeneity through altering interactions that regulate the transcription factor landscape and transcriptome.

Chromatin accessibility differences between alpha, beta, and delta cells identifies common and cell type-specific enhancers

BACKGROUND

High throughput sequencing has enabled the interrogation of the transcriptomic landscape of glucagon-secreting alpha cells, insulin-secreting beta cells, and somatostatin-secreting delta cells. These approaches have furthered our understanding of expression patterns that define healthy or diseased islet cell types and helped explicate some of the intricacies between major islet cell crosstalk and glucose regulation. All three endocrine cell types derive from a common pancreatic progenitor, yet alpha and beta cells have partially opposing functions, and delta cells modulate and control insulin and glucagon release. While gene expression signatures that define and maintain cellular identity have been widely explored, the underlying epigenetic components are incompletely characterized and understood. However, chromatin accessibility and remodeling is a dynamic attribute that plays a critical role to determine and maintain cellular identity.

RESULTS

Here, we compare and contrast the chromatin landscape between mouse alpha, beta, and delta cells using ATAC-Seq to evaluate the significant differences in chromatin accessibility. The similarities and differences in chromatin accessibility between these related islet endocrine cells help define their fate in support of their distinct functional roles. We identify patterns that suggest that both alpha and delta cells are poised, but repressed, from becoming beta-like. We also identify patterns in differentially enriched chromatin that have transcription factor motifs preferentially associated with different regions of the genome. Finally, we not only confirm and visualize previously discovered common endocrine- and cell specific- enhancer regions across differentially enriched chromatin, but identify novel regions as well. We compiled our chromatin accessibility data in a freely accessible database of common endocrine- and cell specific-enhancer regions that can be navigated with minimal bioinformatics expertise.

CONCLUSIONS

Both alpha and delta cells appear poised, but repressed, from becoming beta cells in murine pancreatic islets. These data broadly support earlier findings on the plasticity in identity of non-beta cells under certain circumstances. Furthermore, differential chromatin accessibility shows preferentially enriched distal-intergenic regions in beta cells, when compared to either alpha or delta cells.

Emerging RUNX2-Mediated Gene Regulatory Mechanisms Consisting of Multi-Layered Regulatory Networks in Skeletal Development

Skeletal development is tightly coordinated by chondrocytes and osteoblasts, which are derived from skeletal progenitors, and distinct cell-type gene regulatory programs underlie the specification and differentiation of cells. Runt-related transcription factor 2 (Runx2) is essential to chondrocyte hypertrophy and osteoblast differentiation. Genetic studies have revealed the biological functions of Runx2 and its involvement in skeletal genetic diseases. Meanwhile, molecular biology has provided a framework for our understanding of RUNX2-mediated transactivation at a limited number of cis-regulatory elements. Furthermore, studies using next-generation sequencing (NGS) have provided information on RUNX2-mediated gene regulation at the genome level and novel insights into the multiple layers of gene regulatory mechanisms, including the modes of action of RUNX2, chromatin accessibility, the concept of pioneer factors and phase separation, and three-dimensional chromatin organization. In this review, I summarize the emerging RUNX2-mediated regulatory mechanism from a multi-layer perspective and discuss future perspectives for applications in the treatment of skeletal diseases.

Regulation of organelle size and organization during development

During early embryogenesis, as cells divide in the developing embryo, the size of intracellular organelles generally decreases to scale with the decrease in overall cell size. Organelle size scaling is thought to be important to establish and maintain proper cellular function, and defective scaling may lead to impaired development and disease. However, how the cell regulates organelle size and organization are largely unanswered questions. In this review, we summarize the process of size scaling at both the cell and organelle levels and discuss recently discovered mechanisms that regulate this process during early embryogenesis. In addition, we describe how some recently developed techniques and Xenopus as an animal model can be used to investigate the underlying mechanisms of size regulation and to uncover the significance of proper organelle size scaling and organization.

Combined method of whole mount and block-face imaging: Acquisition of 3D data of gene expression pattern from conventional in situ hybridization

Visualization of spatiotemporal expression of a gene of interest is a fundamental technique for analyzing the involvements of genes in organ development. In situ hybridization (ISH) is one of the most popular methods for visualizing gene expression. When conventional ISH is performed on sections or whole-mount specimens, the gene expression pattern is represented in 2-dimensional (2D) microscopic images or in the surface view of the specimen. To obtain 3-dimensional (3D) data of gene expression from conventional ISH, the "serial section method" has traditionally been employed. However, this method requires an extensive amount of time and labor because it requires researchers to collect a tremendous number of sections, label all sections by ISH, and image them before 3D reconstruction. Here, we proposed a rapid and low-cost 3D imaging method that can create 3D gene expression patterns from conventional ISH-labeled specimens. Our method consists of a combination of whole-mount ISH and Correlative Microscopy and Blockface imaging (CoMBI). The whole-mount ISH-labeled specimens were sliced using a microtome or cryostat, and all block-faces were imaged and used to reconstruct 3D images by CoMBI. The 3D data acquired using our method showed sufficient quality to analyze the morphology and gene expression patterns in the developing mouse heart. In addition, 2D microscopic images of the sections can be obtained when needed. Correlating 2D microscopic images and 3D data can help annotate gene expression patterns and understand the anatomy of developing organs. These results indicated that our method can be useful in the field of developmental biology.

Activation of HIV-1 proviruses increases downstream chromatin accessibility

It is unclear how the activation of HIV-1 transcription affects chromatin structure. We interrogated chromatin organization both genome-wide and nearby HIV-1 integration sites using Hi-C and ATAC-seq. In conjunction, we analyzed the transcription of the HIV-1 genome and neighboring genes. We found that long-range chromatin contacts did not differ significantly between uninfected cells and those harboring an integrated HIV-1 genome, whether the HIV-1 genome was actively transcribed or inactive. Instead, the activation of HIV-1 transcription changes chromatin accessibility immediately downstream of the provirus, demonstrating that HIV-1 can alter local cellular chromatin structure. Finally, we examined HIV-1 and neighboring host gene transcripts with long-read sequencing and found populations of chimeric RNAs both virus-to-host and host-to-virus. Thus, multiomics profiling revealed that the activation of HIV-1 transcription led to local changes in chromatin organization and altered the expression of neighboring host genes.

A Tremendous Reorganization Journey for the 3D Chromatin Structure from Gametes to Embryos

The 3D chromatin structure within the nucleus is important for gene expression regulation and correct developmental programs. Recently, the rapid development of low-input chromatin conformation capture technologies has made it possible to study 3D chromatin structures in gametes, zygotes and early embryos in a variety of species, including flies, vertebrates and mammals. There are distinct 3D chromatin structures within the male and female gametes. Following the fertilization of male and female gametes, fertilized eggs undergo drastic epigenetic reprogramming at multi levels, including the 3D chromatin structure, to convert the terminally differentiated gamete state into the totipotent state, which can give rise to an individual. However, to what extent the 3D chromatin structure reorganization is evolutionarily conserved and what the underlying mechanisms are for the tremendous reorganization in early embryos remain elusive. Here, we review the latest findings on the 3D chromatin structure reorganization during embryogenesis, and discuss the convergent and divergent reprogramming patterns and key molecular mechanisms for the 3D chromatin structure reorganization from gametes to embryos in different species. These findings shed light on how the 3D chromatin structure reorganization contribute to embryo development in different species. The findings also indicate the role of the 3D chromatin structure on the acquisition of totipotent developmental potential.

Repression and 3D-restructuring resolves regulatory conflicts in evolutionarily rearranged genomes

Regulatory landscapes drive complex developmental gene expression, but it remains unclear how their integrity is maintained when incorporating novel genes and functions during evolution. Here, we investigated how a placental mammal-specific gene, Zfp42, emerged in an ancient vertebrate topologically associated domain (TAD) without adopting or disrupting the conserved expression of its gene, Fat1. In ESCs, physical TAD partitioning separates Zfp42 and Fat1 with distinct local enhancers that drive their independent expression. This separation is driven by chromatin activity and not CTCF/cohesin. In contrast, in embryonic limbs, inactive Zfp42 shares Fat1's intact TAD without responding to active Fat1 enhancers. However, neither Fat1 enhancer-incompatibility nor nuclear envelope-attachment account for Zfp42's unresponsiveness. Rather, Zfp42's promoter is rendered inert to enhancers by context-dependent DNA methylation. Thus, diverse mechanisms enabled the integration of independent Zfp42 regulation in the Fat1 locus. Critically, such regulatory complexity appears common in evolution as, genome wide, most TADs contain multiple independently expressed genes.

Three-dimensional genome organization in immune cell fate and function

Immune cell development and activation demand the precise and coordinated control of transcriptional programmes. Three-dimensional (3D) organization of the genome has emerged as an important regulator of chromatin state, transcriptional activity and cell identity by facilitating or impeding long-range genomic interactions among regulatory elements and genes. Chromatin folding thus enables cell type-specific and stimulus-specific transcriptional responses to extracellular signals, which are essential for the control of immune cell fate, for inflammatory responses and for generating a diverse repertoire of antigen receptor specificities. Here, we review recent findings connecting 3D genome organization to the control of immune cell differentiation and function, and discuss how alterations in genome folding may lead to immune dysfunction and malignancy.

Multi-layered transcriptional control of cranial neural crest development

The neural crest (NC) is an emblematic population of embryonic stem-like cells with remarkable migratory ability. These distinctive attributes have inspired the curiosity of developmental biologists for over 150 years, however only recently the regulatory mechanisms controlling the complex features of the NC have started to become elucidated at genomic scales. Regulatory control of NC development is achieved through combinatorial transcription factor binding and recruitment of associated transcriptional complexes to distal cis-regulatory elements. Together, they regulate when, where and to what extent transcriptional programmes are actively deployed, ultimately shaping ontogenetic processes. Here, we discuss how transcriptional networks control NC ontogeny, with a special emphasis on the molecular mechanisms underlying specification of the cephalic NC. We also cover emerging properties of transcriptional regulation revealed in diverse developmental systems, such as the role of three-dimensional conformation of chromatin, and how they are involved in the regulation of NC ontogeny. Finally, we highlight how advances in deciphering the NC transcriptional network have afforded new insights into the molecular basis of human diseases.

Identification of enhancers responsible for the coordinated expression of myosin heavy chain isoforms in skeletal muscle

Background

Skeletal muscles consist of fibers of differing contractility and metabolic properties, which are primarily determined by the content of myosin heavy chain (MYH) isoforms (MYH7, MYH2, MYH1, and MYH4). The regulation of Myh genes transcription depends on three-dimensional chromatin conformation interaction, but the mechanistic details remain to be determined.

Results

In this study, we characterized the interaction profiles of Myh genes using 4C-seq (circular chromosome conformation capture coupled to high-throughput sequencing). The interaction profile of Myh genes changed between fast quadriceps and slow soleus muscles. Combining chromatin immunoprecipitation-sequencing (ChIP-seq) and transposase accessible chromatin with high-throughput sequencing (ATAC-seq), we found that a 38 kb intergenic region interacting simultaneously with fast Myh genes promoters controlled the coordinated expression of fast Myh genes. We also identified four active enhancers of Myh7, and revealed that binding of MYOG and MYOD increased the activity of Myh7 enhancers.

Conclusions

This study provides new insight into the chromatin interactions that regulate Myh genes expression.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08737-9.

Regulatory Architecture of the RCA Gene Cluster Captures an Intragenic TAD Boundary, CTCF-Mediated Chromatin Looping and a Long-Range Intergenic Enhancer

The Regulators of Complement Activation (RCA) gene cluster comprises several tandemly arranged genes with shared functions within the immune system. RCA members, such as complement receptor 2 (CR2), are well-established susceptibility genes in complex autoimmune diseases. Altered expression of RCA genes has been demonstrated at both the functional and genetic level, but the mechanisms underlying their regulation are not fully characterised. We aimed to investigate the structural organisation of the RCA gene cluster to identify key regulatory elements that influence the expression of CR2 and other genes in this immunomodulatory region. Using 4C, we captured extensive CTCF-mediated chromatin looping across the RCA gene cluster in B cells and showed these were organised into two topologically associated domains (TADs). Interestingly, an inter-TAD boundary was located within the CR1 gene at a well-characterised segmental duplication. Additionally, we mapped numerous gene-gene and gene-enhancer interactions across the region, revealing extensive co-regulation. Importantly, we identified an intergenic enhancer and functionally demonstrated this element upregulates two RCA members (CR2 and CD55) in B cells. We have uncovered novel, long-range mechanisms whereby autoimmune disease susceptibility may be influenced by genetic variants, thus highlighting the important contribution of chromatin topology to gene regulation and complex genetic disease.

3D or Not 3D: Shaping the Genome during Development

One of the most fundamental questions in developmental biology is how one fertilized cell can give rise to a fully mature organism and how gene regulation governs this process. Precise spatiotemporal gene expression is required for development and is believed to be achieved through a complex interplay of sequence-specific information, epigenetic modifications, trans-acting factors, and chromatin folding. Here we review the role of chromatin folding during development, the mechanisms governing 3D genome organization, and how it is established in the embryo. Furthermore, we discuss recent advances and debated questions regarding the contribution of the 3D genome to gene regulation during organogenesis. Finally, we describe the mechanisms that can reshape the 3D genome, including disease-causing structural variations and the emerging view that transposable elements contribute to chromatin organization.

Persistent spectral simplicial complex-based machine learning for chromosomal structural analysis in cellular differentiation

The three-dimensional (3D) chromosomal structure plays an essential role in all DNA-templated processes, including gene transcription, DNA replication and other cellular processes. Although developing chromosome conformation capture (3C) methods, such as Hi-C, which can generate chromosomal contact data characterized genome-wide chromosomal structural properties, understanding 3D genomic nature-based on Hi-C data remains lacking. Here, we propose a persistent spectral simplicial complex (PerSpectSC) model to describe Hi-C data for the first time. Specifically, a filtration process is introduced to generate a series of nested simplicial complexes at different scales. For each of these simplicial complexes, its spectral information can be calculated from the corresponding Hodge Laplacian matrix. PerSpectSC model describes the persistence and variation of the spectral information of the nested simplicial complexes during the filtration process. Different from all previous models, our PerSpectSC-based features provide a quantitative global-scale characterization of chromosome structures and topology. Our descriptors can successfully classify cell types and also cellular differentiation stages for all the 24 types of chromosomes simultaneously. In particular, persistent minimum best characterizes cell types and Dim (1) persistent multiplicity best characterizes cellular differentiation. These results demonstrate the great potential of our PerSpectSC-based models in polymeric data analysis.

Enhancer-Promoter Communication: It's Not Just About Contact

Cis-regulatory elements such as enhancers can be located even a million base pairs away from their cognate promoter and yet modulate gene transcription. Indeed, the 3D organisation of chromatin enables the establishment of long-range enhancer-promoter communication. The observation of long-range enhancer-promoter chromatin loops at active genes originally led to a model in which enhancers and promoters form physical contacts between each other to control transcription. Yet, recent microscopy data has challenged this prevailing activity-by-contact model of enhancer-promoter communication in transcriptional activation. Live single-cell imaging approaches do not systematically reveal a correlation between enhancer-proximity and transcriptional activation. We therefore discuss the need to move from a static to a dynamic view of enhancer-promoter relationships. We highlight recent studies that not only reveal considerable chromatin movement in specific cell types, but suggest links between chromatin compaction, chromatin movement and transcription. We describe the interplay between enhancer-promoter proximity within the context of biomolecular condensates and the need to understand how condensate microenvironments influence the chromatin binding kinetics of proteins that bind at cis-regulatory elements to activate transcription. Finally, given the complex multi-scale interplay between regulatory proteins, enhancer-promoter proximity and movement, we propose the need to integrate information from complementary single-cell next-generation sequencing and live-cell imaging approaches to derive unified 3D theoretical models of enhancer-promoter communication that are ultimately predictive of transcriptional output and cell fate. In time, improved models will shed light on how tissues grow and diseases emerge.

Cis-Regulatory Control of Mammalian Sex Determination

Sex determination is the process by which an initial bipotential gonad adopts either a testicular or ovarian cell fate. The inability to properly complete this process leads to a group of developmental disorders classified as disorders of sex development (DSD). To date, dozens of genes were shown to play roles in mammalian sex determination, and mutations in these genes can cause DSD in humans or gonadal sex reversal/dysfunction in mice. However, exome sequencing currently provides genetic diagnosis for only less than half of DSD patients. This points towards a major role for the non-coding genome during sex determination. In this review, we highlight recent advances in our understanding of non-coding, cis-acting gene regulatory elements and discuss how they may control transcriptional programmes that underpin sex determination in the context of the 3-dimensional folding of chromatin. As a paradigm, we focus on the Sox9 gene, a prominent pro-male factor and one of the most extensively studied genes in gonadal cell fate determination.

Genetically engineered birds; pre-CRISPR and CRISPR era

Generating biopharmaceuticals in genetically engineered bioreactors continues to reign supreme. Hence, genetically engineered birds have attracted considerable attention from the biopharmaceutical industry. Fairly recent genome engineering methods have made genome manipulation an easy and affordable task. In this review, we first provide a broad overview of the approaches and main impediments ahead of generating efficient and reliable genetically engineered birds, and various factors that affect the fate of a transgene. This section provides an essential background for the rest of the review, in which we discuss and compare different genome manipulation methods in the pre-CRISPR and CRISPR era in the field of avian genome engineering.

BRD4 orchestrates genome folding to promote neural crest differentiation

Higher-order chromatin structure regulates gene expression, and mutations in proteins mediating genome folding underlie developmental disorders known as cohesinopathies. However, the relationship between three-dimensional genome organization and embryonic development remains unclear. Here we define a role for bromodomain-containing protein 4 (BRD4) in genome folding, and leverage it to understand the importance of genome folding in neural crest progenitor differentiation. Brd4 deletion in neural crest results in cohesinopathy-like phenotypes. BRD4 interacts with NIPBL, a cohesin agonist, and BRD4 depletion or loss of the BRD4-NIPBL interaction reduces NIPBL occupancy, suggesting that BRD4 stabilizes NIPBL on chromatin. Chromatin interaction mapping and imaging experiments demonstrate that BRD4 depletion results in compromised genome folding and loop extrusion. Finally, mutation of individual BRD4 amino acids that mediate an interaction with NIPBL impedes neural crest differentiation into smooth muscle. Remarkably, loss of WAPL, a cohesin antagonist, rescues attenuated smooth muscle differentiation resulting from BRD4 loss. Collectively, our data reveal that BRD4 choreographs genome folding and illustrates the relevance of balancing cohesin activity for progenitor differentiation.

Global chromatin relabeling accompanies spatial inversion of chromatin in rod photoreceptors

The nuclear architecture of rod photoreceptor cells in nocturnal mammals is unlike that of other animal cells. Murine rod cells have an “inverted” chromatin organization with euchromatin at the nuclear periphery and heterochromatin packed in the center of the nucleus. In conventional nuclear architecture, euchromatin is mostly in the interior, and heterochromatin is largely at the nuclear periphery. We demonstrate that inverted nuclear architecture is achieved through global relabeling of the rod cell epigenome. During rod cell maturation, H3K9me2-labeled nuclear peripheral heterochromatin is relabeled with H3K9me3 and repositioned to the nuclear center, while transcriptionally active euchromatin is labeled with H3K9me2 and positioned at the nuclear periphery. Global chromatin relabeling is correlated with spatial rearrangement, suggesting a critical role for histone modifications, specifically H3K9 methylation, in nuclear architecture. These results reveal a dramatic example of genome-wide epigenetic relabeling of chromatin that accompanies altered nuclear architecture in a postnatal, postmitotic cell.

Chromatin Conformation in Development and Disease

Chromatin domains and loops are important elements of chromatin structure and dynamics, but much remains to be learned about their exact biological role and nature. Topological associated domains and functional loops are key to gene expression and hold the answer to many questions regarding developmental decisions and diseases. Here, we discuss new findings, which have linked chromatin conformation with development, differentiation and diseases and hypothesized on various models while integrating all recent findings on how chromatin architecture affects gene expression during development, evolution and disease.

Understanding 3D genome organization by multidisciplinary methods

Understanding how chromatin is folded in the nucleus is fundamental to understanding its function. Although 3D genome organization has been historically difficult to study owing to a lack of relevant methodologies, major technological breakthroughs in genome-wide mapping of chromatin contacts and advances in imaging technologies in the twenty-first century considerably improved our understanding of chromosome conformation and nuclear architecture. In this Review, we discuss methods of 3D genome organization analysis, including sequencing-based techniques, such as Hi-C and its derivatives, Micro-C, DamID and others; microscopy-based techniques, such as super-resolution imaging coupled with fluorescence in situ hybridization (FISH), multiplex FISH, in situ genome sequencing and live microscopy methods; and computational and modelling approaches. We describe the most commonly used techniques and their contribution to our current knowledge of nuclear architecture and, finally, we provide a perspective on up-and-coming methods that open possibilities for future major discoveries.

The nuclear periphery is a scaffold for tissue-specific enhancers

Nuclear architecture influences gene regulation and cell identity by controlling the three-dimensional organization of genes and their distal regulatory sequences, which may be far apart in linear space. The genome is functionally and spatially segregated in the eukaryotic nucleus with transcriptionally active regions in the nuclear interior separated from repressive regions, including those at the nuclear periphery. Here, we describe the identification of a novel type of nuclear peripheral chromatin domain that is enriched for tissue-specific transcriptional enhancers. Like other chromatin at the nuclear periphery, these regions are marked by H3K9me2. But unlike the nuclear peripheral Lamina-Associated Domains (LADs), these novel, enhancer-rich domains have limited Lamin B interaction. We therefore refer to them as H3K9me2-Only Domains (KODs). In mouse embryonic stem cells, KODs are found in Hi-C-defined A compartments and feature relatively accessible chromatin. KODs are characterized by low gene expression and enhancers located in these domains bear the histone marks of an inactive or poised state. These results indicate that KODs organize a subset of inactive, tissue-specific enhancers at the nuclear periphery. We hypothesize that KODs may play a role in facilitating and perhaps constraining the enhancer-promoter interactions underlying spatiotemporal regulation of gene expression programs in differentiation and development.

Mitochondrial DNA Methylation and Human Diseases

Epigenetic modifications of the nuclear genome, including DNA methylation, histone modifications and non-coding RNA post-transcriptional regulation, are increasingly being involved in the pathogenesis of several human diseases. Recent evidence suggests that also epigenetic modifications of the mitochondrial genome could contribute to the etiology of human diseases. In particular, altered methylation and hydroxymethylation levels of mitochondrial DNA (mtDNA) have been found in animal models and in human tissues from patients affected by cancer, obesity, diabetes and cardiovascular and neurodegenerative diseases. Moreover, environmental factors, as well as nuclear DNA genetic variants, have been found to impair mtDNA methylation patterns. Some authors failed to find DNA methylation marks in the mitochondrial genome, suggesting that it is unlikely that this epigenetic modification plays any role in the control of the mitochondrial function. On the other hand, several other studies successfully identified the presence of mtDNA methylation, particularly in the mitochondrial displacement loop (D-loop) region, relating it to changes in both mtDNA gene transcription and mitochondrial replication. Overall, investigations performed until now suggest that methylation and hydroxymethylation marks are present in the mtDNA genome, albeit at lower levels compared to those detectable in nuclear DNA, potentially contributing to the mitochondria impairment underlying several human diseases.

A TAD Skeptic: Is 3D Genome Topology Conserved?

The notion that topologically associating domains (TADs) are highly conserved across species is prevalent in the field of 3D genomics. However, what exactly is meant by 'highly conserved' and what are the actual comparative data that support this notion? To address these questions, we performed a historical review of the relevant literature and retraced numerous citation chains to reveal the primary data that were used as the basis for the widely accepted conclusion that TADs are highly conserved across evolution. A thorough review of the available evidence suggests the answer may be more complex than what is commonly presented.

Germ cell differentiation requires Tdrd7-dependent chromatin and transcriptome reprogramming marked by germ plasm relocalization

In many animal models, primordial germ cell (PGC) development depends on maternally deposited germ plasm, which prevents somatic cell fate. Here, we show that PGCs respond to regulatory information from the germ plasm in two distinct phases using two distinct mechanisms in zebrafish. We demonstrate that PGCs commence zygotic genome activation together with the somatic blastocysts with no demonstrable differences in transcriptional and chromatin opening. Unexpectedly, both PGC and somatic blastocysts activate germ-cell-specific genes, which are only stabilized in PGCs by cytoplasmic germ plasm determinants. Disaggregated perinuclear relocalization of germ plasm during PGC migration is regulated by the germ plasm determinant Tdrd7 and is coupled to dramatic divergence between PGC and somatic transcriptomes. This transcriptional divergence relies on PGC-specific cis-regulatory elements characterized by promoter-proximal distribution. We show that Tdrd7-dependent reconfiguration of chromatin accessibility is required for elaboration of PGC fate but not for PGC migration.

•

No evidence for transcriptional activation delay in zebrafish PGCs

•

Germ-plasm-associated post-transcriptional divergence during ZGA

•

Epigenetic reprogramming marks onset of PGC migration

•

Epigenetic reprogramming in PGCs relies on Tdrd7, coupled to germ plasm relocalization

Dynamics of genome architecture and chromatin function during human B cell differentiation and neoplastic transformation

To investigate the three-dimensional (3D) genome architecture across normal B cell differentiation and in neoplastic cells from different subtypes of chronic lymphocytic leukemia and mantle cell lymphoma patients, here we integrate in situ Hi-C and nine additional omics layers. Beyond conventional active (A) and inactive (B) compartments, we uncover a highly-dynamic intermediate compartment enriched in poised and polycomb-repressed chromatin. During B cell development, 28% of the compartments change, mostly involving a widespread chromatin activation from naive to germinal center B cells and a reversal to the naive state upon further maturation into memory B cells. B cell neoplasms are characterized by both entity and subtype-specific alterations in 3D genome organization, including large chromatin blocks spanning key disease-specific genes. This study indicates that 3D genome interactions are extensively modulated during normal B cell differentiation and that the genome of B cell neoplasias acquires a tumor-specific 3D genome architecture.

Global Genome Conformational Programming during Neuronal Development Is Associated with CTCF and Nuclear FGFR1-The Genome Archipelago Model

During the development of mouse embryonic stem cells (ESC) to neuronal committed cells (NCC), coordinated changes in the expression of 2851 genes take place, mediated by the nuclear form of FGFR1. In this paper, widespread differences are demonstrated in the ESC and NCC inter- and intra-chromosomal interactions, chromatin looping, the formation of CTCF- and nFGFR1-linked Topologically Associating Domains (TADs) on a genome-wide scale and in exemplary HoxA-D loci. The analysis centered on HoxA cluster shows that blocking FGFR1 disrupts the loop formation. FGFR1 binding and genome locales are predictive of the genome interactions; likewise, chromatin interactions along with nFGFR1 binding are predictive of the genome function and correlate with genome regulatory attributes and gene expression. This study advances a topologically integrated genome archipelago model that undergoes structural transformations through the formation of nFGFR1-associated TADs. The makeover of the TAD islands serves to recruit distinct ontogenic programs during the development of the ESC to NCC.

Muscle progenitor specification and myogenic differentiation are associated with changes in chromatin topology

Using Hi-C, promoter-capture Hi-C (pCHi-C), and other genome-wide approaches in skeletal muscle progenitors that inducibly express a master transcription factor, Pax7, we systematically characterize at high-resolution the spatio-temporal re-organization of compartments and promoter-anchored interactions as a consequence of myogenic commitment and differentiation. We identify key promoter-enhancer interaction motifs, namely, cliques and networks, and interactions that are dependent on Pax7 binding. Remarkably, Pax7 binds to a majority of super-enhancers, and together with a cadre of interacting transcription factors, assembles feed-forward regulatory loops. During differentiation, epigenetic memory and persistent looping are maintained at a subset of Pax7 enhancers in the absence of Pax7. We also identify and functionally validate a previously uncharacterized Pax7-bound enhancer hub that regulates the essential myosin heavy chain cluster during skeletal muscle cell differentiation. Our studies lay the groundwork for understanding the role of Pax7 in orchestrating changes in the three-dimensional chromatin conformation in muscle progenitors.

Epigenetics of neuromuscular disorders

Neuromuscular disorders are a heterogeneous group of conditions affecting the neuromuscular system. The aim of this article is to review the major epigenetic findings in motor neuron diseases and major hereditary muscular dystrophies. DNA methylation changes are observed in both hereditary and sporadic forms, and combining DNA methylation analysis with mutational screening holds the potential for better diagnostic and prognostic accuracy. Novel, less toxic and more selective epigenetic drugs are designed and tested in animal and cell culture models of neuromuscular disorders, and non-coding RNAs are being investigated as either disease biomarkers or targets of therapeutic approaches to restore gene expression levels. Overall, neuromuscular disorder epigenetic biomarkers have a strong potential for clinical applications in the near future.

Phylogenic Determinants of Cardiovascular Frailty, Focus on Hemodynamics and Arterial Smooth Muscle Cells

The evolution of the circulatory system from invertebrates to mammals has involved the passage from an open system to a closed in-parallel system via a closed in-series system, accompanying the increasing complexity and efficiency of life’s biological functions. The archaic heart enables pulsatile motion waves of hemolymph in invertebrates, and the in-series circulation in fish occurs with only an endothelium, whereas mural smooth muscle cells appear later. The present review focuses on evolution of the circulatory system. In particular, we address how and why this evolution took place from a closed, flowing, longitudinal conductance at low pressure to a flowing, highly pressurized and bifurcating arterial compartment. However, although arterial pressure was the latest acquired hemodynamic variable, the general teleonomy of the evolution of species is the differentiation of individual organ function, supported by specific fueling allowing and favoring partial metabolic autonomy. This was achieved via the establishment of an active contractile tone in resistance arteries, which permitted the regulation of blood supply to specific organ activities via its localized function-dependent inhibition (active vasodilation). The global resistance to viscous blood flow is the peripheral increase in frictional forces caused by the tonic change in arterial and arteriolar radius, which backscatter as systemic arterial blood pressure. Consequently, the arterial pressure gradient from circulating blood to the adventitial interstitium generates the unidirectional outward radial advective conductance of plasma solutes across the wall of conductance arteries. This hemodynamic evolution was accompanied by important changes in arterial wall structure, supported by smooth muscle cell functional plasticity, including contractility, matrix synthesis and proliferation, endocytosis and phagocytosis, etc. These adaptive phenotypic shifts are due to epigenetic regulation, mainly related to mechanotransduction. These paradigms actively participate in cardio-arterial pathologies such as atheroma, valve disease, heart failure, aneurysms, hypertension, and physiological aging.

Identification and Characterization of a Transcribed Distal Enhancer Involved in Cardiac Kcnh2 Regulation

The human ether-a-go-go-related gene KCNH2 encodes the voltage-gated potassium channel underlying I_Kr, a current critical for the repolarization phase of the cardiac action potential. Mutations in KCNH2 that cause a reduction of the repolarizing current can result in cardiac arrhythmias associated with long-QT syndrome. Here, we investigate the regulation of KCNH2 and identify multiple active enhancers. A transcribed enhancer ∼85 kbp downstream of Kcnh2 physically contacts the promoters of two Kcnh2 isoforms in a cardiac-specific manner in vivo. Knockdown of its ncRNA transcript results in reduced expression of Kcnh2b and two neighboring mRNAs, Nos3 and Abcb8, in vitro. Genomic deletion of the enhancer, including the ncRNA transcription start site, from the mouse genome causes a modest downregulation of both Kcnh2a and Kcnh2b in the ventricles. These findings establish that the regulation of Kcnh2a and Kcnh2b is governed by a complex regulatory landscape that involves multiple partially redundantly acting enhancers.

Oplr16 serves as a novel chromatin factor to control stem cell fate by modulating pluripotency-specific chromosomal looping and TET2-mediated DNA demethylation

Formation of a pluripotency-specific chromatin network is a critical event in reprogramming somatic cells into pluripotent status. To characterize the regulatory components in this process, we used ‘chromatin RNA in situ reverse transcription sequencing’ (CRIST-seq) to profile RNA components that interact with the pluripotency master gene Oct4. Using this approach, we identified a novel nuclear lncRNA Oplr16 that was closely involved in the initiation of reprogramming. Oplr16 not only interacted with the Oct4 promoter and regulated its activity, but it was also specifically activated during reprogramming to pluripotency. Active expression of Oplr16 was required for optimal maintenance of pluripotency in embryonic stem cells. Oplr16 was also able to enhance reprogramming of fibroblasts into pluripotent cells. RNA reverse transcription-associated trap sequencing (RAT-seq) indicated that Oplr16 interacted with multiple target genes related to stem cell self-renewal. Of note, Oplr16 utilized its 3′-fragment to recruit the chromatin factor SMC1 to orchestrate pluripotency-specific intrachromosomal looping. After binding to the Oct4 promoter, Oplr16 recruited TET2 to induce DNA demethylation and activate Oct4 in fibroblasts, leading to enhanced reprogramming. These data suggest that Oplr16 may act as a pivotal chromatin factor to control stem cell fate by modulating chromatin architecture and DNA demethylation.

Evaluating Enhancer Function and Transcription

Cell-type- and condition-specific profiles of gene expression require coordination between protein-coding gene promoters and cis-regulatory sequences called enhancers. Enhancers can stimulate gene activity at great genomic distances from their targets, raising questions about how enhancers communicate with specific gene promoters and what molecular mechanisms underlie enhancer function. Characterization of enhancer loci has identified the molecular features of active enhancers that accompany the binding of transcription factors and local opening of chromatin. These characteristics include coactivator recruitment, histone modifications, and noncoding RNA transcription. However, it remains unclear which of these features functionally contribute to enhancer activity. Here, we discuss what is known about how enhancers regulate their target genes and how enhancers and promoters communicate. Further, we describe recent data demonstrating many similarities between enhancers and the gene promoters they control, and we highlight unanswered questions in the field, such as the potential roles of transcription at enhancers.

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

Mammalian gene expression patterns are controlled by regulatory elements, which interact within topologically associating domains (TADs). The relationship between activation of regulatory elements, formation of structural chromatin interactions and gene expression during development is unclear. Here, we present Tiled-C, a low-input chromosome conformation capture (3C) technique. We use this approach to study chromatin architecture at high spatial and temporal resolution through in vivo mouse erythroid differentiation. Integrated analysis of chromatin accessibility and single-cell expression data shows that regulatory elements gradually become accessible within pre-existing TADs during early differentiation. This is followed by structural re-organization within the TAD and formation of specific contacts between enhancers and promoters. Our high-resolution data show that these enhancer-promoter interactions are not established prior to gene expression, but formed gradually during differentiation, concomitant with progressive upregulation of gene activity. Together, these results provide new insight into the close, interdependent relationship between chromatin architecture and gene regulation during development.

The corrected gene proximity map for analyzing the 3D genome organization using Hi-C data

BACKGROUND

Genome-wide ligation-based assays such as Hi-C provide us with an unprecedented opportunity to investigate the spatial organization of the genome. Results of a typical Hi-C experiment are often summarized in a chromosomal contact map, a matrix whose elements reflect the co-location frequencies of genomic loci. To elucidate the complex structural and functional interactions between those genomic loci, networks offer a natural and powerful framework.

RESULTS

We propose a novel graph-theoretical framework, the Corrected Gene Proximity (CGP) map to study the effect of the 3D spatial organization of genes in transcriptional regulation. The starting point of the CGP map is a weighted network, the gene proximity map, whose weights are based on the contact frequencies between genes extracted from genome-wide Hi-C data. We derive a null model for the network based on the signal contributed by the 1D genomic distance and use it to "correct" the gene proximity for cell type 3D specific arrangements. The CGP map, therefore, provides a network framework for the 3D structure of the genome on a global scale. On human cell lines, we show that the CGP map can detect and quantify gene co-regulation and co-localization more effectively than the map obtained by raw contact frequencies. Analyzing the expression pattern of metabolic pathways of two hematopoietic cell lines, we find that the relative positioning of the genes, as captured and quantified by the CGP, is highly correlated with their expression change. We further show that the CGP map can be used to form an inter-chromosomal proximity map that allows large-scale abnormalities, such as chromosomal translocations, to be identified.

CONCLUSIONS

The Corrected Gene Proximity map is a map of the 3D structure of the genome on a global scale. It allows the simultaneous analysis of intra- and inter- chromosomal interactions and of gene co-regulation and co-localization more effectively than the map obtained by raw contact frequencies, thus revealing hidden associations between global spatial positioning and gene expression. The flexible graph-based formalism of the CGP map can be easily generalized to study any existing Hi-C datasets.

Three-dimensional genome architectural CCCTC-binding factor makes choice in duplicated enhancers at Pcdhα locus

During development, gene expression is spatiotemporally regulated by long-distance chromatin interactions between distal enhancers and target promoters. However, how specificity of the interactions between enhancers and promoters is achieved remains largely unknown. As there are far more enhancers than promoters in mammalian genomes, the complexities of enhancer choice during gene regulation remain obscure. CTCF, the CCCTC-binding factor that directionally binds to a vast range of genomic sites known as CBSs (CTCF-binding sites), mediates oriented chromatin looping between a substantial set of distal enhancers and target promoters. To investigate mechanisms by which CTCF engages in enhancer choice, we used CRISPR/Cas9-based DNA-fragment editing to duplicate CBS-containing enhancers and promoters in the native genomic locus of the clustered Pcdhα genes. We found that the promoter is regulated by the proximal one among duplicated enhancers and that this choice is dependent on CTCF-mediated directional enhancer-promoter looping. In addition, gene expression is unaltered upon the switch of enhancers. Moreover, after promoter duplication, only the proximal promoter is chosen by CTCF-mediated directional chromatin looping to contact with the distal enhancer. Finally, we demonstrated that both enhancer activation and chromatin looping with the promoter are essential for gene expression. These findings have important implications regarding the role of CTCF in specific interactions between enhancers and promoters as well as developmental regulation of gene expression by enhancer switching.

Differences in the intrinsic spatial dynamics of the chromatin contribute to cell differentiation

Advances in chromosome conformation capture techniques as well as computational characterization of genomic loci structural dynamics open new opportunities for exploring the mechanistic aspects of genome-scale differences across different cell types. We examined here the dynamic basis of variabilities between different cell types by investigating their chromatin mobility profiles inferred from Hi-C data using an elastic network model representation of the chromatin. Our comparative analysis of sixteen cell lines reveals close similarities between chromosomal dynamics across different cell lines on a global scale, but notable cell-specific variations emerge in the detailed spatial mobilities of genomic loci. Closer examination reveals that the differences in spatial dynamics mainly originate from the difference in the frequencies of their intrinsically accessible modes of motion. Thus, even though the chromosomes of different types of cells have access to similar modes of collective movements, not all modes are deployed by all cells, such that the effective mobilities and cross-correlations of genomic loci are cell-type-specific. Comparison with RNA-seq expression data reveals a strong overlap between highly expressed genes and those distinguished by high mobilities in the present study, in support of the role of the intrinsic spatial dynamics of chromatin as a determinant of cell differentiation.

Robust CTCF-Based Chromatin Architecture Underpins Epigenetic Changes in the Heart Failure Stress-Gene Response

BACKGROUND

The human genome folds in 3 dimensions to form thousands of chromatin loops inside the nucleus, encasing genes and cis-regulatory elements for accurate gene expression control. Physical tethers of loops are anchored by the DNA-binding protein CTCF and the cohesin ring complex. Because heart failure is characterized by hallmark gene expression changes, it was recently reported that substantial CTCF-related chromatin reorganization underpins the myocardial stress-gene response, paralleled by chromatin domain boundary changes observed in CTCF knockout.

METHODS

We undertook an independent and orthogonal analysis of chromatin organization with mouse pressure-overload model of myocardial stress (transverse aortic constriction) and cardiomyocyte-specific knockout of Ctcf. We also downloaded published data sets of similar cardiac mouse models and subjected them to independent reanalysis.

RESULTS

We found that the cardiomyocyte chromatin architecture remains broadly stable in transverse aortic constriction hearts, whereas Ctcf knockout resulted in ≈99% abolition of global chromatin loops. Disease gene expression changes correlated instead with differential histone H3K27-acetylation enrichment at their respective proximal and distal interacting genomic enhancers confined within these static chromatin structures. Moreover, coregulated genes were mapped out as interconnected gene sets on the basis of their multigene 3D interactions.

CONCLUSIONS

This work reveals a more stable genome-wide chromatin framework than previously described. Myocardial stress-gene transcription responds instead through H3K27-acetylation enhancer enrichment dynamics and gene networks of coregulation. Robust and intact CTCF looping is required for the induction of a rapid and accurate stress response.