The 3D genome in transcriptional regulation and pluripotency

Rewiring of SINE-MIR enhancer topology and Esrrb modulation in expanded and naive pluripotency

BACKGROUND

The interplay between 3D genomic structure and transposable elements (TE) in regulating cell state-specific gene expression program is largely unknown. Here, we explore the utilization of TE-derived enhancers in naïve and expanded pluripotent states by integrative analysis of genome-wide Hi-C-defined enhancer interactions, H3K27ac HiChIP profiling and CRISPR-guided TE proteomics landscape.

RESULTS

We find that short interspersed nuclear elements (SINEs) are the more involved TEs in the active chromatin and 3D genome architecture. In particular, mammalian-wide interspersed repeat (MIR), a SINE family member, is highly associated with naïve-specific genomic interactions compared to the expanded state. Primarily, in the naïve pluripotent state, MIR enhancer is co-opted by ESRRB for naïve-specific gene expression program. This ESRRB and MIR enhancer interaction is crucial for the formation of loops that build a network of enhancers and super-enhancers regulating pluripotency genes. We demonstrate that loss of a ESRRB-bound MIR enhancer impairs self-renewal. We also find that MIR is co-bound by structural protein complex, ESRRB-YY1, in the naïve pluripotent state.

CONCLUSIONS

Altogether, our study highlights the topological regulation of ESRRB on MIR in the naïve potency state.

The establishment of nuclear organization in mouse embryos is orchestrated by multiple epigenetic pathways

The folding of the genome in the 3D nuclear space is fundamental for regulating all DNA-related processes. The association of the genome with the nuclear lamina into lamina-associated domains (LADs) represents the earliest feature of nuclear organization during development. Here, we performed a gain-of-function screen in mouse embryos to obtain mechanistic insights. We find that perturbations impacting histone H3 modifications, heterochromatin, and histone content are crucial for the establishment of nuclear architecture in zygotes and/or 2-cell-stage embryos. Notably, some perturbations exerted differential effects on zygotes versus 2-cell-stage embryos. Moreover, embryos with disrupted LADs can rebuild nuclear architecture at the 2-cell stage, indicating that the initial establishment of LADs in zygotes might be dispensable for early development. Our findings provide valuable insights into the functional interplay between chromatin and structural components of the nucleus that guide genome-lamina interactions during the earliest developmental stages.

Identification of ZNF384 as a regulator of epigenome in leukemia

Leukemia is a complex hematologic cancer driven by genetic and epigenetic changes that impact gene expression. Understanding these molecular mechanisms is essential for improving leukemia diagnosis and prognosis. This study examines the role of the zinc finger protein ZNF384 in the epigenome and its influence on gene regulation in leukemia. We analyzed next-generation sequencing data from The Encyclopedia of DNA Elements (ENCODE), integrating datasets such as chromatin immunoprecipitation sequencing (ChIP-seq) of ZNF384 and regulatory histone marks, RNA sequencing (RNA-seq), and Hi-C data from K562 and GM12878 cells. Additionally, we used RNA-seq from K562 ZNF384 knock-down (KD) cells generated via CRISPR interference (CRISPRi) to validate our findings. This enabled us to explore the chromatin interaction patterns of ZNF384 and its regulatory impact. Our results demonstrate that ZNF384 associates with promoters and enhancers in K562 and GM12878 cells, facilitating increased transcription levels. We also found ZNF384 enriched at topologically associating domain (TAD) boundaries and chromatin loops, suggesting a role in three-dimensional (3D) chromatin organization. Furthermore, we identified a significant binding of ZNF384 at SINE-Alu elements in both K562 and GM12878 cells. In summary, this study highlights the regulatory role of ZNF384 in the leukemia epigenome and its impact on gene expression. Understanding the oncogenic implications of ZNF384 may improve leukemia diagnosis and prognosis.

Tracing the Shared Foundations of Gene Expression and Chromatin Structure

The three-dimensional organization of chromatin into topologically associating domains (TADs) may impact gene regulation by bringing distant genes into contact. However, many questions about TADs' function and their influence on transcription remain unresolved due to technical limitations in defining TAD boundaries and measuring the direct effect that TADs have on gene expression. Here, we develop consensus TAD maps for human and mouse with a novel "bag-of-genes" approach for defining the gene composition within TADs. This approach enables new functional interpretations of TADs by providing a way to capture species-level differences in chromatin organization. We also leverage a generative AI foundation model computed from 33 million transcriptomes to define contextual similarity, an embedding-based metric that is more powerful than co-expression at representing functional gene relationships. Our analytical framework directly leads to testable hypotheses about chromatin organization across cellular states. We find that TADs play an active role in facilitating gene co-regulation, possibly through a mechanism involving transcriptional condensates. We also discover that the TAD-linked enhancement of transcriptional context is strongest in early developmental stages and systematically declines with aging. Investigation of cancer cells show distinct patterns of TAD usage that shift with chemotherapy treatment, suggesting specific roles for TAD-mediated regulation in cellular development and plasticity. Finally, we develop "TAD signatures" to improve statistical analysis of single-cell transcriptomic data sets in predicting cancer cell-line drug response. These findings reshape our understanding of cellular plasticity in development and disease, indicating that chromatin organization acts through probabilistic mechanisms rather than deterministic rules.

Software availability

https://singhlab.net/tadmap.

DconnLoop: a deep learning model for predicting chromatin loops based on multi-source data integration

BACKGROUND

Chromatin loops are critical for the three-dimensional organization of the genome and gene regulation. Accurate identification of chromatin loops is essential for understanding the regulatory mechanisms in disease. However, current mainstream detection methods rely primarily on single-source data, such as Hi-C, which limits these methods' ability to capture the diverse features of chromatin loop structures. In contrast, multi-source data integration and deep learning approaches, though not yet widely applied, hold significant potential.

RESULTS

In this study, we developed a method called DconnLoop to integrate Hi-C, ChIP-seq, and ATAC-seq data to predict chromatin loops. This method achieves feature extraction and fusion of multi-source data by integrating residual mechanisms, directional connectivity excitation modules, and interactive feature space decoders. Finally, we apply density estimation and density clustering to the genome-wide prediction results to identify more representative loops. The code is available from https://github.com/kuikui-C/DconnLoop .

CONCLUSIONS

The results demonstrate that DconnLoop outperforms existing methods in both precision and recall. In various experiments, including Aggregate Peak Analysis and peak enrichment comparisons, DconnLoop consistently shows advantages. Extensive ablation studies and validation across different sequencing depths further confirm DconnLoop's robustness and generalizability.

Evolutionary divergence in CTCF-mediated chromatin topology drives transcriptional innovation in humans

Chromatin topology can impact gene regulation, but how evolutionary divergence in chromatin topology has shaped gene regulatory landscapes for distinctive human traits remains poorly understood. CTCF sites determine chromatin topology by forming domains and loops. Here, we show evolutionary divergence in CTCF-mediated chromatin topology at the domain and loop scales during primate evolution, elucidating distinct mechanisms for shaping regulatory landscapes. Human-specific divergent domains lead to a broad rewiring of transcriptional landscapes. Divergent CTCF loops concord with species-specific enhancer activity, influencing enhancer connectivity to target genes in a concordant yet constrained manner. Under this concordant mechanism, we establish the role of human-specific CTCF loops in shaping transcriptional isoform diversity, with functional implications for disease susceptibility. Furthermore, we validate the function of these human-specific CTCF loops using human forebrain organoids. This study advances our understanding of genetic evolution from the perspective of genome architecture.

Deciphering single-cell genomic architecture: insights into cellular heterogeneity and regulatory dynamics

BACKGROUND

The genomic architecture of eukaryotes exhibits dynamic spatial and temporal changes, enabling cellular processes critical for maintaining viability and functional diversity. Recent advances in sequencing technologies have facilitated the dissection of genomic architecture and functional activity at single-cell resolution, moving beyond the averaged signals typically derived from bulk cell analyses.

MAIN BODY

The advent of single-cell genomics and epigenomics has yielded transformative insights into cellular heterogeneity, behavior, and biological complexity with unparalleled genomic resolution and reproducibility. This review summarizes recent progress in the characterization of genomic architecture at the single-cell level, emphasizing the impact of structural variation and chromatin organization on gene regulatory networks and cellular identity.

CONCLUSION

Future directions in single-cell genomics and high-resolution epigenomic methodologies are explored, focusing on emerging challenges and potential impacts on the understanding of cellular states, regulatory dynamics, and the intricate mechanisms driving cellular function and diversity. Future perspectives on the challenges and potential implications of single-cell genomics, along with high-resolution genomic and epigenomic technologies for understanding cellular states and regulatory dynamics, are also discussed.

ChromoGen: Diffusion model predicts single-cell chromatin conformations

Breakthroughs in high-throughput sequencing and microscopic imaging technologies have revealed that chromatin structures vary considerably between cells of the same type. However, a thorough characterization of this heterogeneity remains elusive due to the labor-intensive and time-consuming nature of these experiments. To address these challenges, we introduce ChromoGen, a generative model based on state-of-the-art artificial intelligence techniques that efficiently predicts three-dimensional, single-cell chromatin conformations de novo with both region and cell type specificity. These generated conformations accurately reproduce experimental results at both the single-cell and population levels. Moreover, ChromoGen successfully transfers to cell types excluded from the training data using just DNA sequence and widely available DNase-seq data, thus providing access to chromatin structures in myriad cell types. These achievements come at a remarkably low computational cost. Therefore, ChromoGen enables the systematic investigation of single-cell chromatin organization, its heterogeneity, and its relationship to sequencing data, all while remaining economical.

Chromatin phase separation and nuclear shape fluctuations are correlated in a polymer model of the nucleus

Abnormal cell nuclear shapes are hallmarks of diseases, including progeria, muscular dystrophy, and many cancers. Experiments have shown that disruption of heterochromatin and increases in euchromatin lead to nuclear deformations, such as blebs and ruptures. However, the physical mechanisms through which chromatin governs nuclear shape are poorly understood. To investigate how heterochromatin and euchromatin might govern nuclear morphology, we studied chromatin microphase separation in a composite coarse-grained polymer and elastic shell simulation model. By varying chromatin density, heterochromatin composition, and heterochromatin-lamina interactions, we show how the chromatin phase organization may perturb nuclear shape. Increasing chromatin density stabilizes the lamina against large fluctuations. However, increasing heterochromatin levels or heterochromatin-lamina interactions enhances nuclear shape fluctuations by a "wetting"-like interaction. In contrast, fluctuations are insensitive to heterochromatin's internal structure. Our simulations suggest that peripheral heterochromatin accumulation could perturb nuclear morphology, while nuclear shape stabilization likely occurs through mechanisms other than chromatin microphase organization.

Droplet Hi-C enables scalable, single-cell profiling of chromatin architecture in heterogeneous tissues

Current methods for analyzing chromatin architecture are not readily scalable to heterogeneous tissues. Here we introduce Droplet Hi-C, which uses a commercial microfluidic device for high-throughput, single-cell chromatin conformation profiling in droplets. Using Droplet Hi-C, we mapped the chromatin architecture of the mouse cortex and analyzed gene regulatory programs in major cortical cell types. In addition, we used this technique to detect copy number variations, structural variations and extrachromosomal DNA in human glioblastoma, colorectal and blood cancer cells, revealing clonal dynamics and other oncogenic events during treatment. We refined the technique to allow joint profiling of chromatin architecture and transcriptome in single cells, facilitating exploration of the links between chromatin architecture and gene expression in both normal tissues and tumors. Thus, Droplet Hi-C both addresses critical gaps in chromatin analysis of heterogeneous tissues and enhances understanding of gene regulation.

Polymer Physics Models Reveal Structural Folding Features of Single-Molecule Gene Chromatin Conformations

Here, we employ polymer physics models of chromatin to investigate the 3D folding of a 2 Mb wide genomic region encompassing the human LTN1 gene, a crucial DNA locus involved in key cellular functions. Through extensive Molecular Dynamics simulations, we reconstruct in silico the ensemble of single-molecule LTN1 3D structures, which we benchmark against recent in situ Hi-C 2.0 data. The model-derived single molecules are then used to predict structural folding features at the single-cell level, providing testable predictions for super-resolution microscopy experiments.

OpenNucleome for high-resolution nuclear structural and dynamical modeling

The intricate structural organization of the human nucleus is fundamental to cellular function and gene regulation. Recent advancements in experimental techniques, including high-throughput sequencing and microscopy, have provided valuable insights into nuclear organization. Computational modeling has played significant roles in interpreting experimental observations by reconstructing high-resolution structural ensembles and uncovering organization principles. However, the absence of standardized modeling tools poses challenges for furthering nuclear investigations. We present OpenNucleome-an open-source software designed for conducting GPU-accelerated molecular dynamics simulations of the human nucleus. OpenNucleome offers particle-based representations of chromosomes at a resolution of 100 KB, encompassing nuclear lamina, nucleoli, and speckles. This software furnishes highly accurate structural models of nuclear architecture, affording the means for dynamic simulations of condensate formation, fusion, and exploration of non-equilibrium effects. We applied OpenNucleome to uncover the mechanisms driving the emergence of 'fixed points' within the nucleus-signifying genomic loci robustly anchored in proximity to specific nuclear bodies for functional purposes. This anchoring remains resilient even amidst significant fluctuations in chromosome radial positions and nuclear shapes within individual cells. Our findings lend support to a nuclear zoning model that elucidates genome functionality. We anticipate OpenNucleome to serve as a valuable tool for nuclear investigations, streamlining mechanistic explorations and enhancing the interpretation of experimental observations.

Assessing next-generation sequencing-based computational methods for predicting transcriptional regulators with query gene sets

This article provides an in-depth review of computational methods for predicting transcriptional regulators (TRs) with query gene sets. Identification of TRs is of utmost importance in many biological applications, including but not limited to elucidating biological development mechanisms, identifying key disease genes, and predicting therapeutic targets. Various computational methods based on next-generation sequencing (NGS) data have been developed in the past decade, yet no systematic evaluation of NGS-based methods has been offered. We classified these methods into two categories based on shared characteristics, namely library-based and region-based methods. We further conducted benchmark studies to evaluate the accuracy, sensitivity, coverage, and usability of NGS-based methods with molecular experimental datasets. Results show that BART, ChIP-Atlas, and Lisa have relatively better performance. Besides, we point out the limitations of NGS-based methods and explore potential directions for further improvement.

From Nucleosomes to Compartments: Physicochemical Interactions Underlying Chromatin Organization

Chromatin organization plays a critical role in cellular function by regulating access to genetic information. However, understanding chromatin folding is challenging due to its complex, multiscale nature. Significant progress has been made in studying in vitro systems, uncovering the structure of individual nucleosomes and their arrays, and elucidating the role of physicochemical forces in stabilizing these structures. Additionally, remarkable advancements have been achieved in characterizing chromatin organization in vivo, particularly at the whole-chromosome level, revealing important features such as chromatin loops, topologically associating domains, and nuclear compartments. However, bridging the gap between in vitro and in vivo studies remains challenging. The resemblance between in vitro and in vivo chromatin conformations and the relevance of internucleosomal interactions for chromatin folding in vivo are subjects of debate. This article reviews experimental and computational studies conducted at various length scales, highlighting the significance of intrinsic interactions between nucleosomes and their roles in chromatin folding in vivo.

Yin Yang 1 regulates cohesin complex protein SMC3 in mouse hematopoietic stem cells

ABSTRACT

Yin Yang 1 (YY1) and structural maintenance of chromosomes 3 (SMC3) are 2 critical chromatin structural factors that mediate long-distance enhancer-promoter interactions and promote developmentally regulated changes in chromatin architecture in hematopoietic stem/progenitor cells (HSPCs). Although YY1 has critical functions in promoting hematopoietic stem cell (HSC) self-renewal and maintaining HSC quiescence, SMC3 is required for proper myeloid lineage differentiation. However, many questions remain unanswered regarding how YY1 and SMC3 interact with each other and affect hematopoiesis. We found that YY1 physically interacts with SMC3 and cooccupies with SMC3 at a large cohort of promoters genome wide, and YY1 deficiency deregulates the genetic network governing cell metabolism. YY1 occupies the Smc3 promoter and represses SMC3 expression in HSPCs. Although deletion of 1 Smc3 allele partially restores HSC numbers and quiescence in YY1 knockout mice, Yy1-/-Smc3+/- HSCs fail to reconstitute blood after bone marrow transplant. YY1 regulates HSC metabolic pathways and maintains proper intracellular reactive oxygen species levels in HSCs, and this regulation is independent of the YY1-SMC3 axis. Our results establish a distinct YY1-SMC3 axis and its impact on HSC quiescence and metabolism.

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.

Multiscale 3D genome organization underlies duck fatty liver with no adipose inflammation or serious injury

Non-alcoholic fatty liver disease (NAFLD) is the most common form of chronic liver disease. Little is known about how gene expression and chromatin structure are regulated in NAFLD due to lack of suitable model. Ducks naturally develop fatty liver similar to serious human non-alcoholic fatty liver (NAFL) without adipose inflammation and liver fibrosis, thus serves as a good model for investigating molecular mechanisms of adipose metabolism and anti-inflammation. Here, we constructed a NAFLD model without adipose inflammation and liver fibrosis in ducks. By performing dynamic pathological and transcriptomic analyses, we identified critical genes involving in regulation of the NF-κB and MHCII signaling, which usually lead to adipose inflammation and liver fibrosis. We further generated dynamic three-dimensional chromatin maps during liver fatty formation and recovery. This showed that ducks enlarged hepatocyte cell nuclei to reduce inter-chromosomal interaction, decompress chromatin structure, and alter strength of intra-TAD and loop interactions during fatty liver formation. These changes partially contributed to the tight control the NF-κB and the MHCII signaling. Our analysis uncovers duck chromatin reorganization might be advantageous to maintain liver regenerative capacity and reduce adipose inflammation. These findings shed light on new strategies for NAFLD control.

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.

CD-Loop: a chromatin loop detection method based on the diffusion model

Motivation

In recent years, there have been significant advances in various chromatin conformation capture techniques, and annotating the topological structure from Hi-C contact maps has become crucial for studying the three-dimensional structure of chromosomes. However, the structure and function of chromatin loops are highly dynamic and diverse, influenced by multiple factors. Therefore, obtaining the three-dimensional structure of the genome remains a challenging task. Among many chromatin loop prediction methods, it is difficult to fully extract features from the contact map and make accurate predictions at low sequencing depths.

Results

In this study, we put forward a deep learning framework based on the diffusion model called CD-Loop for predicting accurate chromatin loops. First, by pre-training the input data, we obtain prior probabilities for predicting the classification of the Hi-C contact map. Then, by combining the denoising process based on the diffusion model and the prior probability obtained by pre-training, candidate loops were predicted from the input Hi-C contact map. Finally, CD-Loop uses a density-based clustering algorithm to cluster the candidate chromatin loops and predict the final chromatin loops. We compared CD-Loop with the currently popular methods, such as Peakachu, Chromosight, and Mustache, and found that in different cell types, species, and sequencing depths, CD-Loop outperforms other methods in loop annotation. We conclude that CD-Loop can accurately predict chromatin loops and reveal cell-type specificity. The code is available at https://github.com/wangyang199897/CD-Loop.

H3K4me1 facilitates promoter-enhancer interactions and gene activation during embryonic stem cell differentiation

Histone H3 lysine 4 mono-methylation (H3K4me1) marks poised or active enhancers. KMT2C (MLL3) and KMT2D (MLL4) catalyze H3K4me1, but their histone methyltransferase activities are largely dispensable for transcription during early embryogenesis in mammals. To better understand the role of H3K4me1 in enhancer function, we analyze dynamic enhancer-promoter (E-P) interactions and gene expression during neural differentiation of the mouse embryonic stem cells. We found that KMT2C/D catalytic activities were only required for H3K4me1 and E-P contacts at a subset of candidate enhancers, induced upon neural differentiation. By contrast, a majority of enhancers retained H3K4me1 in KMT2C/D catalytic mutant cells. Surprisingly, H3K4me1 signals at these KMT2C/D-independent sites were reduced after acute depletion of KMT2B, resulting in aggravated transcriptional defects. Our observations therefore implicate KMT2B in the catalysis of H3K4me1 at enhancers and provide additional support for an active role of H3K4me1 in enhancer-promoter interactions and transcription in mammalian cells.

Droplet Hi-C for Fast and Scalable Profiling of Chromatin Architecture in Single Cells

Comprehensive analysis of chromatin architecture is crucial for understanding the gene regulatory programs during development and in disease pathogenesis, yet current methods often inadequately address the unique challenges presented by analysis of heterogeneous tissue samples. Here, we introduce Droplet Hi-C, which employs a commercial microfluidic device for high-throughput, single-cell chromatin conformation profiling in droplets. Using Droplet Hi-C, we mapped the chromatin architecture at single-cell resolution from the mouse cortex and analyzed gene regulatory programs in major cortical cell types. Additionally, we used this technique to detect copy number variation (CNV), structural variations (SVs) and extrachromosomal DNA (ecDNA) in cancer cells, revealing clonal dynamics and other oncogenic events during treatment. We further refined this technique to allow for joint profiling of chromatin architecture and transcriptome in single cells, facilitating a more comprehensive exploration of the links between chromatin architecture and gene expression in both normal tissues and tumors. Thus, Droplet Hi-C not only addresses critical gaps in chromatin analysis of heterogeneous tissues but also emerges as a versatile tool enhancing our understanding of gene regulation in health and disease.

OpenNucleome for high resolution nuclear structural and dynamical modeling

The intricate structural organization of the human nucleus is fundamental to cellular function and gene regulation. Recent advancements in experimental techniques, including high-throughput sequencing and microscopy, have provided valuable insights into nuclear organization. Computational modeling has played significant roles in interpreting experimental observations by reconstructing high-resolution structural ensembles and uncovering organization principles. However, the absence of standardized modeling tools poses challenges for furthering nuclear investigations. We present OpenNucleome-an open-source software designed for conducting GPU-accelerated molecular dynamics simulations of the human nucleus. OpenNucleome offers particle-based representations of chromosomes at a resolution of 100 KB, encompassing nuclear lamina, nucleoli, and speckles. This software furnishes highly accurate structural models of nuclear architecture, affording the means for dynamic simulations of condensate formation, fusion, and exploration of non-equilibrium effects. We applied OpenNucleome to uncover the mechanisms driving the emergence of "fixed points" within the nucleus-signifying genomic loci robustly anchored in proximity to specific nuclear bodies for functional purposes. This anchoring remains resilient even amidst significant fluctuations in chromosome radial positions and nuclear shapes within individual cells. Our findings lend support to a nuclear zoning model that elucidates genome functionality. We anticipate OpenNucleome to serve as a valuable tool for nuclear investigations, streamlining mechanistic explorations and enhancing the interpretation of experimental observations.

Assessing NGS-based computational methods for predicting transcriptional regulators with query gene sets

This article provides an in-depth review of computational methods for predicting transcriptional regulators with query gene sets. Identification of transcriptional regulators is of utmost importance in many biological applications, including but not limited to elucidating biological development mechanisms, identifying key disease genes, and predicting therapeutic targets. Various computational methods based on next-generation sequencing (NGS) data have been developed in the past decade, yet no systematic evaluation of NGS-based methods has been offered. We classified these methods into two categories based on shared characteristics, namely library-based and region-based methods. We further conducted benchmark studies to evaluate the accuracy, sensitivity, coverage, and usability of NGS-based methods with molecular experimental datasets. Results show that BART, ChIP-Atlas, and Lisa have relatively better performance. Besides, we point out the limitations of NGS-based methods and explore potential directions for further improvement.

Key points

An introduction to available computational methods for predicting functional TRs from a query gene set.A detailed walk-through along with practical concerns and limitations.A systematic benchmark of NGS-based methods in terms of accuracy, sensitivity, coverage, and usability, using 570 TR perturbation-derived gene sets.NGS-based methods outperform motif-based methods. Among NGS methods, those utilizing larger databases and adopting region-centric approaches demonstrate favorable performance. BART, ChIP-Atlas, and Lisa are recommended as these methods have overall better performance in evaluated scenarios.

Giant pandas in captivity undergo short-term adaptation in nerve-related pathways

BACKGROUND

Behaviors in captive animals, including changes in appetite, activity level, and social interaction, are often seen as adaptive responses. However, these behaviors may become progressively maladaptive, leading to stress, anxiety, depression, and other negative reactions in animals.

RESULTS

In this study, we investigated the whole-genome sequencing data of 39 giant panda individuals, including 11 in captivity and 28 in the wild. To eliminate the mountain range effect and focus on the factor of captivity only, we first performed a principal component analysis. We then enumerated the 21,474,180 combinations of wild giant pandas (11 chosen from 28) and calculated their distances from the 11 captive individuals. The 11 wild individuals with the closest distances were used for the subsequent analysis. The linkage disequilibrium (LD) patterns demonstrated that the population was almost eliminated. We identified 505 robust selected genomic regions harboring at least one SNP, and the absolute frequency difference was greater than 0.6 between the two populations. GO analysis revealed that genes in these regions were mainly involved in nerve-related pathways. Furthermore, we identified 22 GO terms for which the selection strength significantly differed between the two populations, and there were 10 nerve-related pathways among them. Genes in the differentially abundant regions were involved in nerve-related pathways, indicating that giant pandas in captivity underwent minor genomic selection. Additionally, we investigated the relationship between genetic variation and chromatin conformation structures. We found that nucleotide diversity (θπ) in the captive population was correlated with chromatin conformation structures, which included A/B compartments, topologically associated domains (TADs) and TAD-cliques. For each GO term, we then compared the expression level of genes regulated by the above four factors (AB index, TAD intactness, TAD clique and PEI) with the corresponding genomic background. The retained 10 GO terms were all coordinately regulated by the four factors, and three of them were associated with nerve-related pathways.

CONCLUSIONS

This study revealed that giant pandas in captivity undergo short-term adaptation in nerve-related pathways. Furthermore, it provides new insights into the molecular mechanism of gene expression regulation under short-term adaptation to environmental change.

Chromatin phase separation and nuclear shape fluctuations are correlated in a polymer model of the nucleus

Abnormalities in the shapes of mammalian cell nuclei are hallmarks of a variety of diseases, including progeria, muscular dystrophy, and various cancers. Experiments have shown that there is a causal relationship between chromatin organization and nuclear morphology. Decreases in heterochromatin levels, perturbations to heterochromatin organization, and increases in euchromatin levels all lead to misshapen nuclei, which exhibit deformations, such as nuclear blebs and nuclear ruptures. However, the polymer physical mechanisms of how chromatin governs nuclear shape and integrity are poorly understood. To investigate how heterochromatin and euchromatin, which are thought to microphase separate in vivo , govern nuclear morphology, we implemented a composite coarse-grained polymer and elastic shell model. By varying chromatin volume fraction (density), heterochromatin levels and structure, and heterochromatin-lamina interactions, we show how the spatial organization of chromatin polymer phases within the nucleus could perturb nuclear shape in some scenarios. Increasing the volume fraction of chromatin in the cell nucleus stabilizes the nuclear lamina against large fluctuations. However, surprisingly, we find that increasing heterochromatin levels or heterochromatin-lamina interactions enhances nuclear shape fluctuations in our simulations by a "wetting"-like interaction. In contrast, shape fluctuations are largely insensitive to the internal structure of the heterochromatin, such as the presence or absence of chromatin-chromatin crosslinks. Therefore, our simulations suggest that heterochromatin accumulation at the nuclear periphery could perturb nuclear morphology in a nucleus or nuclear region that is sufficiently soft, while stabilization of the nucleus via heterochromatin likely occurs through mechanisms other than chromatin microphase organization.

Hi-BDiSCO: folding 3D mesoscale genome structures from Hi-C data using brownian dynamics

The structure and dynamics of the eukaryotic genome are intimately linked to gene regulation and transcriptional activity. Many chromosome conformation capture experiments like Hi-C have been developed to detect genome-wide contact frequencies and quantify loop/compartment structures for different cellular contexts and time-dependent processes. However, a full understanding of these events requires explicit descriptions of representative chromatin and chromosome configurations. With the exponentially growing amount of data from Hi-C experiments, many methods for deriving 3D structures from contact frequency data have been developed. Yet, most reconstruction methods use polymer models with low resolution to predict overall genome structure. Here we present a Brownian Dynamics (BD) approach termed Hi-BDiSCO for producing 3D genome structures from Hi-C and Micro-C data using our mesoscale-resolution chromatin model based on the Discrete Surface Charge Optimization (DiSCO) model. Our approach integrates reconstruction with chromatin simulations at nucleosome resolution with appropriate biophysical parameters. Following a description of our protocol, we present applications to the NXN, HOXC, HOXA and Fbn2 mouse genes ranging in size from 50 to 100 kb. Such nucleosome-resolution genome structures pave the way for pursuing many biomedical applications related to the epigenomic regulation of chromatin and control of human disease.

Cell-type differential targeting of SETDB1 prevents aberrant CTCF binding, chromatin looping, and cis-regulatory interactions

SETDB1 is an essential histone methyltransferase that deposits histone H3 lysine 9 trimethylation (H3K9me3) to transcriptionally repress genes and repetitive elements. The function of differential H3K9me3 enrichment between cell-types remains unclear. Here, we demonstrate mutual exclusivity of H3K9me3 and CTCF across mouse tissues from different developmental timepoints. We analyze SETDB1 depleted cells and discover that H3K9me3 prevents aberrant CTCF binding independently of DNA methylation and H3K9me2. Such sites are enriched with SINE B2 retrotransposons. Moreover, analysis of higher-order genome architecture reveals that large chromatin structures including topologically associated domains and subnuclear compartments, remain intact in SETDB1 depleted cells. However, chromatin loops and local 3D interactions are disrupted, leading to transcriptional changes by modifying pre-existing chromatin landscapes. Specific genes with altered expression show differential interactions with dysregulated cis-regulatory elements. Collectively, we find that cell-type specific targets of SETDB1 maintain cellular identities by modulating CTCF binding, which shape nuclear architecture and transcriptomic networks.

3D Enhancer-promoter networks provide predictive features for gene expression and coregulation in early embryonic lineages

Mammalian embryogenesis commences with two pivotal and binary cell fate decisions that give rise to three essential lineages: the trophectoderm, the epiblast and the primitive endoderm. Although key signaling pathways and transcription factors that control these early embryonic decisions have been identified, the non-coding regulatory elements through which transcriptional regulators enact these fates remain understudied. Here, we characterize, at a genome-wide scale, enhancer activity and 3D connectivity in embryo-derived stem cell lines that represent each of the early developmental fates. We observe extensive enhancer remodeling and fine-scale 3D chromatin rewiring among the three lineages, which strongly associate with transcriptional changes, although distinct groups of genes are irresponsive to topological changes. In each lineage, a high degree of connectivity, or 'hubness', positively correlates with levels of gene expression and enriches for cell-type specific and essential genes. Genes within 3D hubs also show a significantly stronger probability of coregulation across lineages compared to genes in linear proximity or within the same contact domains. By incorporating 3D chromatin features, we build a predictive model for transcriptional regulation (3D-HiChAT) that outperforms models using only 1D promoter or proximal variables to predict levels and cell-type specificity of gene expression. Using 3D-HiChAT, we identify, in silico, candidate functional enhancers and hubs in each cell lineage, and with CRISPRi experiments, we validate several enhancers that control gene expression in their respective lineages. Our study identifies 3D regulatory hubs associated with the earliest mammalian lineages and describes their relationship to gene expression and cell identity, providing a framework to comprehensively understand lineage-specific transcriptional behaviors.

Phase Separation as a Driver of Stem Cell Organization and Function during Development

A properly organized subcellular composition is essential to cell function. The canonical organizing principle within eukaryotic cells involves membrane-bound organelles; yet, such structures do not fully explain cellular complexity. Furthermore, discrete non-membrane-bound structures have been known for over a century. Liquid-liquid phase separation (LLPS) has emerged as a ubiquitous mode of cellular organization without the need for formal lipid membranes, with an ever-expanding and diverse list of cellular functions that appear to be regulated by this process. In comparison to traditional organelles, LLPS can occur across wider spatial and temporal scales and involves more distinct protein and RNA complexes. In this review, we discuss the impacts of LLPS on the organization of stem cells and their function during development. Specifically, the roles of LLPS in developmental signaling pathways, chromatin organization, and gene expression will be detailed, as well as its impacts on essential processes of asymmetric cell division. We will also discuss how the dynamic and regulated nature of LLPS may afford stem cells an adaptable mode of organization throughout the developmental time to control cell fate. Finally, we will discuss how aberrant LLPS in these processes may contribute to developmental defects and disease.

Perturb-tracing enables high-content screening of multiscale 3D genome regulators

Three-dimensional (3D) genome organization becomes altered during development, aging, and disease1-23, but the factors regulating chromatin topology are incompletely understood and currently no technology can efficiently screen for new regulators of multiscale chromatin organization. Here, we developed an image-based high-content screening platform (Perturb-tracing) that combines pooled CRISPR screen, a new cellular barcode readout method (BARC-FISH), and chromatin tracing. We performed a loss-of-function screen in human cells, and visualized alterations to their genome organization from 13,000 imaging target-perturbation combinations, alongside perturbation-paired barcode readout in the same single cells. Using 1.4 million 3D positions along chromosome traces, we discovered tens of new regulators of chromatin folding at different length scales, ranging from chromatin domains and compartments to chromosome territory. A subset of the regulators exhibited 3D genome effects associated with loop-extrusion and A-B compartmentalization mechanisms, while others were largely unrelated to these known 3D genome mechanisms. We found that the ATP-dependent helicase CHD7, the loss of which causes the congenital neural crest syndrome CHARGE24 and a chromatin remodeler previously shown to promote local chromatin openness25-27, counter-intuitively compacts chromatin over long range in different genomic contexts and cell backgrounds including neural crest cells, and globally represses gene expression. The DNA compaction effect of CHD7 is independent of its chromatin remodeling activity and does not require other protein partners. Finally, we identified new regulators of nuclear architectures and found a functional link between chromatin compaction and nuclear shape. Altogether, our method enables scalable, high-content identification of chromatin and nuclear topology regulators that will stimulate new insights into the 3D genome functions, such as global gene and nuclear regulation, in health and disease.

Quantifying the large-scale chromosome structural dynamics during the mitosis-to-G1 phase transition of cell cycle

Cell cycle is known to be regulated by the underlying gene network. Chromosomes, which serve as the scaffold for gene expressions, undergo significant structural reorganizations during mitosis. Understanding the mechanism of the cell cycle from the chromosome structural perspective remains a grand challenge. In this study, we applied an integrated theoretical approach to investigate large-scale chromosome structural dynamics during the mitosis-to-G1 phase transition. We observed that the chromosome structural expansion and adaptation of the structural asphericity do not occur synchronously and attributed this behaviour to the unique unloading sequence of the two types of condensins. Furthermore, we observed that the coherent motions between the chromosomal loci are primarily enhanced within the topologically associating domains (TADs) as cells progress to the G1 phase, suggesting that TADs can be considered as both structural and dynamical units for organizing the three-dimensional chromosome. Our analysis also reveals that the quantified pathways of chromosome structural reorganization during the mitosis-to-G1 phase transition exhibit high stochasticity at the single-cell level and show nonlinear behaviours in changing TADs and contacts formed at the long-range regions. Our findings offer valuable insights into large-scale chromosome structural dynamics after mitosis.

Multiple lineage-specific epigenetic landscapes at the antigen receptor loci

Antigen receptors (AgRs) expressed on B and T cells provide the adaptive immune system with ability to detect numerous foreign antigens. Epigenetic features of B cell receptor (BCR) and T cell receptor (TCR) genes were previously studied in lymphocytes, but little is known about their epigenetic features in other cells. Here, we explored histone modifications and transcription markers at the BCR and TCR loci in lymphocytes (pro-B, DP T cells, and mature CD4+ T cells), compared to embryonic stem (ES) cells and neurons. In B cells, the BCR loci exhibited active histone modifications and transcriptional markers indicative of active loci. Similar results were observed at the TCR loci in T cells. All loci were largely inactive in neurons. Surprisingly, in ES cells all AgR loci displayed a high degree of active histone modifications and markers of active transcription. Locations of these active histone modifications in ES cells were largely distinct from those in pro-B cells, and co-localized at numerous binding locations for transcription factors Oct4, Sox2, and Nanog. ES and pro-B cells also showed distinct binding patterns for the ubiquitous transcription factor YY1 and chromatin remodeler Brg1. On the contrary, there were many overlapping CCCTC-binding factor (CTCF) binding patterns when comparing ES cells, pro-B cells, and neurons. Our study identifies epigenetic features in ES cells and lymphocytes that may be related to ES cell pluripotency and lymphocyte tissue-specific activation at the AgR loci.

3D organization of regulatory elements for transcriptional regulation in Arabidopsis

BACKGROUND

Although spatial organization of compartments and topologically associating domains at large scale is relatively well studied, the spatial organization of regulatory elements at fine scale is poorly understood in plants.

RESULTS

Here we perform high-resolution chromatin interaction analysis using paired-end tag sequencing approach. We map chromatin interactions tethered with RNA polymerase II and associated with heterochromatic, transcriptionally active, and Polycomb-repressive histone modifications in Arabidopsis. Analysis of the regulatory repertoire shows that distal active cis-regulatory elements are linked to their target genes through long-range chromatin interactions with increased expression of the target genes, while poised cis-regulatory elements are linked to their target genes through long-range chromatin interactions with depressed expression of the target genes. Furthermore, we demonstrate that transcription factor MYC2 is critical for chromatin spatial organization, and propose that MYC2 occupancy and MYC2-mediated chromatin interactions coordinately facilitate transcription within the framework of 3D chromatin architecture. Analysis of functionally related gene-defined chromatin connectivity networks reveals that genes implicated in flowering-time control are functionally compartmentalized into separate subdomains via their spatial activity in the leaf or shoot apical meristem, linking active mark- or Polycomb-repressive mark-associated chromatin conformation to coordinated gene expression.

CONCLUSION

The results reveal that the regulation of gene transcription in Arabidopsis is not only by linear juxtaposition, but also by long-range chromatin interactions. Our study uncovers the fine scale genome organization of Arabidopsis and the potential roles of such organization in orchestrating transcription and development.

Crosstalk between RNA m6A and DNA methylation regulates transposable element chromatin activation and cell fate in human pluripotent stem cells

Transposable elements (TEs) are parasitic DNA sequences accounting for over half of the human genome. Tight control of the repression and activation states of TEs is critical for genome integrity, development, immunity and diseases, including cancer. However, precisely how this regulation is achieved remains unclear. Here we develop a targeted proteomic proximity labeling approach to capture TE-associated proteins in human embryonic stem cells (hESCs). We find that the RNA N6-methyladenosine (m6A) reader, YTHDC2, occupies genomic loci of the primate-specific TE, LTR7/HERV-H, specifically through its interaction with m6A-modified HERV-H RNAs. Unexpectedly, YTHDC2 recruits the DNA 5-methylcytosine (5mC)-demethylase, TET1, to remove 5mC from LTR7/HERV-H and prevent epigenetic silencing. Functionally, the YTHDC2/LTR7 axis inhibits neural differentiation of hESCs. Our results reveal both an underappreciated crosstalk between RNA m6A and DNA 5mC, the most abundant regulatory modifications of RNA and DNA in eukaryotes, and the fact that in hESCs this interplay controls TE activity and cell fate.

Cell-type-specific prediction of 3D chromatin organization enables high-throughput in silico genetic screening

Investigating how chromatin organization determines cell-type-specific gene expression remains challenging. Experimental methods for measuring three-dimensional chromatin organization, such as Hi-C, are costly and have technical limitations, restricting their broad application particularly in high-throughput genetic perturbations. We present C.Origami, a multimodal deep neural network that performs de novo prediction of cell-type-specific chromatin organization using DNA sequence and two cell-type-specific genomic features-CTCF binding and chromatin accessibility. C.Origami enables in silico experiments to examine the impact of genetic changes on chromatin interactions. We further developed an in silico genetic screening approach to assess how individual DNA elements may contribute to chromatin organization and to identify putative cell-type-specific trans-acting regulators that collectively determine chromatin architecture. Applying this approach to leukemia cells and normal T cells, we demonstrate that cell-type-specific in silico genetic screening, enabled by C.Origami, can be used to systematically discover novel chromatin regulation circuits in both normal and disease-related biological systems.

Systematic mapping and modeling of 3D enhancer-promoter interactions in early mouse embryonic lineages reveal regulatory principles that determine the levels and cell-type specificity of gene expression

Mammalian embryogenesis commences with two pivotal and binary cell fate decisions that give rise to three essential lineages, the trophectoderm (TE), the epiblast (EPI) and the primitive endoderm (PrE). Although key signaling pathways and transcription factors that control these early embryonic decisions have been identified, the non-coding regulatory elements via which transcriptional regulators enact these fates remain understudied. To address this gap, we have characterized, at a genome-wide scale, enhancer activity and 3D connectivity in embryo-derived stem cell lines that represent each of the early developmental fates. We observed extensive enhancer remodeling and fine-scale 3D chromatin rewiring among the three lineages, which strongly associate with transcriptional changes, although there are distinct groups of genes that are irresponsive to topological changes. In each lineage, a high degree of connectivity or "hubness" positively correlates with levels of gene expression and enriches for cell-type specific and essential genes. Genes within 3D hubs also show a significantly stronger probability of coregulation across lineages, compared to genes in linear proximity or within the same contact domains. By incorporating 3D chromatin features, we build a novel predictive model for transcriptional regulation (3D-HiChAT), which outperformed models that use only 1D promoter or proximal variables in predicting levels and cell-type specificity of gene expression. Using 3D-HiChAT, we performed genome-wide in silico perturbations to nominate candidate functional enhancers and hubs in each cell lineage, and with CRISPRi experiments we validated several novel enhancers that control expression of one or more genes in their respective lineages. Our study comprehensively identifies 3D regulatory hubs associated with the earliest mammalian lineages and describes their relationship to gene expression and cell identity, providing a framework to understand lineage-specific transcriptional behaviors.

Systematic characterization of regulatory variants of blood pressure genes

High blood pressure (BP) is the major risk factor for cardiovascular disease. Genome-wide association studies have identified genetic variants for BP, but functional insights into causality and related molecular mechanisms lag behind. We functionally characterize 4,608 genetic variants in linkage with 135 BP loci in vascular smooth muscle cells and cardiomyocytes by massively parallel reporter assays. High densities of regulatory variants at BP loci (i.e., ULK4, MAP4, CFDP1, PDE5A) indicate that multiple variants drive genetic association. Regulatory variants are enriched in repeats, alter cardiovascular-related transcription factor motifs, and spatially converge with genes controlling specific cardiovascular pathways. Using heuristic scoring, we define likely causal variants, and CRISPR prime editing finally determines causal variants for KCNK9, SFXN2, and PCGF6, which are candidates for developing high BP. Our systems-level approach provides a catalog of functionally relevant variants and their genomic architecture in two trait-relevant cell lines for a better understanding of BP gene regulation.

Ultrastructure and fractal property of chromosomes in close-to-native yeast nuclei visualized using X-ray laser diffraction

Genome compaction and activity in the nucleus depend on spatiotemporal changes in the organization of chromatins in chromosomes. However, the direct imaging of the chromosome structures in the nuclei has been difficult and challenging. Herein, we directly visualized the structure of chromosomes in frozen-hydrated nuclei of budding yeast in the interphase using X-ray laser diffraction. The reconstructed projection electron density maps revealed inhomogeneous distributions of chromosomes, such as a 300 nm assembly and fibrous substructures in the elliptic-circular shaped nuclei of approximately 800 nm. In addition, from the diffraction patterns, we confirmed the absence of regular arrangements of chromosomes and chromatins with 400-20 nm spacing, and demonstrated that chromosomes were composed of self-similarly assembled substructural domains with an average radius of gyration of 58 nm and smooth surfaces. Based on these analyses, we constructed putative models to discuss the organization of 16 chromosomes, carrying DNA of 4.1 mm in 800 nm ellipsoid of the nucleus at the interphase. We anticipate the structural parameters on the fractal property of chromosomes and the experimental images to be a starting point for constructing more sophisticated 3D structural models of the nucleus.

The circadian clock, aging and its implications in cancer

Circadian clock orchestrates the intergenic biochemical, physiological and behavioral changes to form an approximate 24h oscillation through the transcription-translation feedback loop (TTFL). Mechanistically, a heterodimer of transcriptional activator formed by BMAL1 and CLOCK, governs the expression of its transcriptional repressors, CRY, PER and REV-ERBα/β proteins, thereby controlling more than 50 % of protein encoding genes in human. There is also increasing evidence showing that tumor microenvironment can disrupt specific clock gene functions to facilitate tumorigenesis. Although there is great progress in understanding the molecular mechanisms of the circadian clock, aging and cancer, elucidating their complex relationships among these processes remains challenging. Herein, the optimization of the chronochemotherapy regimen has not been justified yet for treatment of cancer. Here, we discuss the hypothesis of relocalization of chromatin modifiers (RCM) along with function(s) of the circadian rhythm on aging and carcinogenesis. We will also introduce the function of the chromatin remodeling as a new avenue for rejuvenation of competent tissues to combat aging and cancer.

Protein arginine methyltransferase 5 (Prmt5) localizes to chromatin loop anchors and modulates expression of genes at TAD boundaries during early adipogenesis

Protein arginine methyltransferase 5 (Prmt5) is an essential regulator of embryonic development and adult progenitor cell functions. Prmt5 expression is mis-regulated in many cancers, and the development of Prmt5 inhibitors as cancer therapeutics is an active area of research. Prmt5 functions via effects on gene expression, splicing, DNA repair, and other critical cellular processes. We examined whether Prmt5 functions broadly as a genome-wide regulator of gene transcription and higher-order chromatin interactions during the initial stages of adipogenesis using ChIP-Seq, RNA-seq, and Hi-C using 3T3-L1 cells, a frequently utilized model for adipogenesis. We observed robust genome-wide Prmt5 chromatin-binding at the onset of differentiation. Prmt5 localized to transcriptionally active genomic regions, acting as both a positive and a negative regulator. A subset of Prmt5 binding sites co-localized with mediators of chromatin organization at chromatin loop anchors. Prmt5 knockdown decreased insulation strength at the boundaries of topologically associating domains (TADs) adjacent to sites with Prmt5 and CTCF co-localization. Genes overlapping such weakened TAD boundaries showed transcriptional dysregulation. This study identifies Prmt5 as a broad regulator of gene expression, including regulation of early adipogenic factors, and reveals an unappreciated requirement for Prmt5 in maintaining strong insulation at TAD boundaries and overall chromatin organization.

Dynamic enhancer interactome promotes senescence and aging

Gene expression programs are regulated by enhancers which act in a context-specific manner, and can reside at great distances from their target genes. Extensive three-dimensional (3D) genome reorganization occurs in senescence, but how enhancer interactomes are reconfigured during this process is just beginning to be understood. Here we generated high-resolution contact maps of active enhancers and their target genes, assessed chromatin accessibility, and established one-dimensional maps of various histone modifications and transcription factors to comprehensively understand the regulation of enhancer configuration during senescence. Hyper-connected enhancer communities/cliques formed around genes that are highly expressed and within essential gene pathways in each cell state. In addition, motif analysis indicates the involvement of specific transcription factors in hyper-connected regulatory elements in each condition; importantly, MafK, a bZIP family transcription factor, was upregulated in senescence, and reduced expression of MafK ameliorated the senescence phenotypes. Because the accumulation of senescent cells is a key feature of aging, we further investigated enhancer connectomes in the liver of young and aged mice. Hyper-connected enhancer communities were identified during aging, which regulate essential genes that maintain cell differentiation and homeostasis. These findings reveal that hyper-connected enhancer communities correlate with high gene expression in senescence and aging and provide potential hotspots for therapeutic intervention in aging and age-associated diseases.

A multi-omics integrative analysis based on CRISPR screens re-defines the pluripotency regulatory network in ESCs

A comprehensive and precise definition of the pluripotency gene regulatory network (PGRN) is crucial for clarifying the regulatory mechanisms in embryonic stem cells (ESCs). Here, after a CRISPR/Cas9-based functional genomics screen and integrative analysis with other functional genomes, transcriptomes, proteomes and epigenome data, an expanded pluripotency-associated gene set is obtained, and a new PGRN with nine sub-classes is constructed. By integrating the DNA binding, epigenetic modification, chromatin conformation, and RNA expression profiles, the PGRN is resolved to six functionally independent transcriptional modules (CORE, MYC, PAF, PRC, PCGF and TBX). Spatiotemporal transcriptomics reveal activated CORE/MYC/PAF module activity and repressed PRC/PCGF/TBX module activity in both mouse ESCs (mESCs) and pluripotent cells of early embryos. Moreover, this module activity pattern is found to be shared by human ESCs (hESCs) and cancers. Thus, our results provide novel insights into elucidating the molecular basis of ESC pluripotency.

Compartmentalization with nuclear landmarks yields random, yet precise, genome organization

The 3D organization of eukaryotic genomes plays an important role in genome function. While significant progress has been made in deciphering the folding mechanisms of individual chromosomes, the principles of the dynamic large-scale spatial arrangement of all chromosomes inside the nucleus are poorly understood. We use polymer simulations to model the diploid human genome compartmentalization relative to nuclear bodies such as nuclear lamina, nucleoli, and speckles. We show that a self-organization process based on a cophase separation between chromosomes and nuclear bodies can capture various features of genome organization, including the formation of chromosome territories, phase separation of A/B compartments, and the liquid property of nuclear bodies. The simulated 3D structures quantitatively reproduce both sequencing-based genomic mapping and imaging assays that probe chromatin interaction with nuclear bodies. Importantly, our model captures the heterogeneous distribution of chromosome positioning across cells while simultaneously producing well-defined distances between active chromatin and nuclear speckles. Such heterogeneity and preciseness of genome organization can coexist due to the nonspecificity of phase separation and the slow chromosome dynamics. Together, our work reveals that the cophase separation provides a robust mechanism for us to produce functionally important 3D contacts without requiring thermodynamic equilibration that can be difficult to achieve.

Introduction to stem cells

Stem cells have self-renewal capability and can proliferate and differentiate into a variety of functionally active cells that can serve in various tissues and organs. This review discusses the history, definition, and classification of stem cells. Human pluripotent stem cells (hPSCs) mainly include embryonic stem cells (hESCs) and induced pluripotent stem cells (hiPSCs). Embryonic stem cells are derived from the inner cell mass of the embryo. Induced pluripotent stem cells are derived from reprogramming somatic cells. Pluripotent stem cells have the ability to differentiate into cells derived from all three germ layers (endoderm, mesoderm, and ectoderm). Adult stem cells can be multipotent or unipotent and can produce tissue-specific terminally differentiated cells. Stem cells can be used in cell therapy to replace and regenerate damaged tissues or organs.

The gut microbiota links disease to human genome evolution

A large number of studies have established a causal relationship between the gut microbiota and human disease. In addition, the composition of the microbiota is substantially influenced by the human genome. Modern medical research has confirmed that the pathogenesis of various diseases is closely related to evolutionary events in the human genome. Specific regions of the human genome known as human accelerated regions (HARs) have evolved rapidly over several million years since humans diverged from a common ancestor with chimpanzees, and HARs have been found to be involved in some human-specific diseases. Furthermore, the HAR-regulated gut microbiota has undergone rapid changes during human evolution. We propose that the gut microbiota may serve as an important mediator linking diseases to human genome evolution.

Loss of epigenetic information as a cause of mammalian aging

All living things experience an increase in entropy, manifested as a loss of genetic and epigenetic information. In yeast, epigenetic information is lost over time due to the relocalization of chromatin-modifying proteins to DNA breaks, causing cells to lose their identity, a hallmark of yeast aging. Using a system called "ICE" (inducible changes to the epigenome), we find that the act of faithful DNA repair advances aging at physiological, cognitive, and molecular levels, including erosion of the epigenetic landscape, cellular exdifferentiation, senescence, and advancement of the DNA methylation clock, which can be reversed by OSK-mediated rejuvenation. These data are consistent with the information theory of aging, which states that a loss of epigenetic information is a reversible cause of aging.

Dynamic chromatin organization and regulatory interactions in human endothelial cell differentiation

Vascular endothelial cells are a mesoderm-derived lineage with many essential functions, including angiogenesis and coagulation. The gene-regulatory mechanisms underpinning endothelial specialization are largely unknown, as are the roles of chromatin organization in regulating endothelial cell transcription. To investigate the relationships between chromatin organization and gene expression, we induced endothelial cell differentiation from human pluripotent stem cells and performed Hi-C and RNA-sequencing assays at specific time points. Long-range intrachromosomal contacts increase over the course of differentiation, accompanied by widespread heteroeuchromatic compartment transitions that are tightly associated with transcription. Dynamic topologically associating domain boundaries strengthen and converge on an endothelial cell state, and function to regulate gene expression. Chromatin pairwise point interactions (DNA loops) increase in frequency during differentiation and are linked to the expression of genes essential to vascular biology. Chromatin dynamics guide transcription in endothelial cell development and promote the divergence of endothelial cells from cardiomyocytes.

A Comparison of Topologically Associating Domain Callers Based on Hi-C Data

Topologically associating domains (TADs) are local chromatin interaction domains, which have been shown to play an important role in gene expression regulation. TADs were originally discovered in the investigation of 3D genome organization based on High-throughput Chromosome Conformation Capture (Hi-C) data. Continuous considerable efforts have been dedicated to developing methods for detecting TADs from Hi-C data. Different computational methods for TADs identification vary in their assumptions and criteria in calling TADs. As a consequence, the TADs called by these methods differ in their similarities and biological features they are enriched in. In this work, we performed a systematic comparison of twenty-six TAD callers. We first compared the TADs and gaps between adjacent TADs across different methods, resolutions, and sequencing depths. We then assessed the quality of TADs and TAD boundaries according to three criteria: the decay of contact frequencies over the genomic distance, enrichment and depletion of regulatory elements around TAD boundaries, and reproducibility of TADs and TAD boundaries in replicate samples. Last, due to the lack of a gold standard of TADs, we also evaluated the performance of the methods on synthetic datasets. We discussed the key principles of TAD callers, and pinpointed current situation in the detection of TADs. We provide a concise, comprehensive, and systematic framework for evaluating the performance of TAD callers, and expect our work will provide useful guidance in choosing suitable approaches for the detection and evaluation of TADs.

Chondrocyte De-Differentiation: Biophysical Cues to Nuclear Alterations

Autologous chondrocyte implantation (ACI) is a cell therapy to repair cartilage defects. In ACI a biopsy is taken from a non-load bearing area of the knee and expanded in-vitro. The expansion process provides the benefit of generating a large number of cells required for implantation; however, during the expansion these cells de-differentiate and lose their chondrocyte phenotype. In this review we focus on examining the de-differentiation phenotype from a mechanobiology and biophysical perspective, highlighting some of the nuclear mechanics and chromatin changes in chondrocytes seen during the expansion process and how this relates to the gene expression profile. We propose that manipulating chondrocyte nuclear architecture and chromatin organization will highlight mechanisms that will help to preserve the chondrocyte phenotype.

Insights into the cell fate decision-making processes from chromosome structural reorganizations

The cell fate decision-making process, which provides the capability of a cell transition to a new cell type, involves the reorganizations of 3D genome structures. Currently, the high temporal resolution picture of how the chromosome structural rearrangements occur and further influence the gene activities during the cell-state transition is still challenging to acquire. Here, we study the chromosome structural reorganizations during the cell-state transitions among the pluripotent embryonic stem cell, the terminally differentiated normal cell, and the cancer cell using a nonequilibrium landscape-switching model implemented in the molecular dynamics simulation. We quantify the chromosome (de)compaction pathways during the cell-state transitions and find that the two pathways having the same destinations can merge prior to reaching the final states. The chromosomes at the merging states have similar structural geometries but can differ in long-range compartment segregation and spatial distribution of the chromosomal loci and genes, leading to cell-type-specific transition mechanisms. We identify the irreversible pathways of chromosome structural rearrangements during the forward and reverse transitions connecting the same pair of cell states, underscoring the critical roles of nonequilibrium dynamics in the cell-state transitions. Our results contribute to the understanding of the cell fate decision-making processes from the chromosome structural perspective.