Chromatin architecture reorganization during stem cell differentiation

SATB2 organizes the 3D genome architecture of cognition in cortical neurons

The DNA-binding protein SATB2 is genetically linked to human intelligence. We studied its influence on the three-dimensional (3D) epigenome by mapping chromatin interactions and accessibility in control versus SATB2-deficient cortical neurons. We find that SATB2 affects the chromatin looping between enhancers and promoters of neuronal-activity-regulated genes, thus influencing their expression. It also alters A/B compartments, topologically associating domains, and frequently interacting regions. Genes linked to SATB2-dependent 3D genome changes are implicated in highly specialized neuronal functions and contribute to cognitive ability and risk for neuropsychiatric and neurodevelopmental disorders. Non-coding DNA regions with a SATB2-dependent structure are enriched for common variants associated with educational attainment, intelligence, and schizophrenia. Our data establish SATB2 as a cell-type-specific 3D genome modulator, which operates both independently and in cooperation with CCCTC-binding factor (CTCF) to set up the chromatin landscape of pyramidal neurons for cognitive processes.

Methylation of histone H3 lysine 36 is a barrier for therapeutic interventions of head and neck squamous cell carcinoma

Approximately 20% of head and neck squamous cell carcinomas (HNSCCs) exhibit reduced methylation on lysine 36 of histone H3 (H3K36me) due to mutations in histone methylase NSD1 or a lysine-to-methionine mutation in histone H3 (H3K36M). Whether such alterations of H3K36me can be exploited for therapeutic interventions is still unknown. Here, we show that HNSCC models expressing H3K36M can be divided into two groups: those that display aberrant accumulation of H3K27me3 and those that maintain steady levels of H3K27me3. The former group exhibits reduced proliferation, genome instability, and heightened sensitivity to genotoxic agents like PARP1/2 inhibitors. Conversely, H3K36M HNSCC models with constant H3K27me3 levels lack these characteristics unless H3K27me3 is elevated by DNA hypomethylating agents or inhibiting H3K27me3 demethylases KDM6A/B. Mechanistically, H3K36M reduces H3K36me by directly impeding the activities of the histone methyltransferase NSD3 and the histone demethylase LSD2. Notably, aberrant H3K27me3 levels induced by H3K36M expression are not a bona fide epigenetic mark because they require continuous expression of H3K36M to be inherited. Moreover, increased sensitivity to PARP1/2 inhibitors in H3K36M HNSCC models depends solely on elevated H3K27me3 levels and diminishing BRCA1- and FANCD2-dependent DNA repair. Finally, a PARP1/2 inhibitor alone reduces tumor burden in a H3K36M HNSCC xenograft model with elevated H3K27me3, whereas in a model with consistent H3K27me3, a combination of PARP1/2 inhibitors and agents that up-regulate H3K27me3 proves to be successful. These findings underscore the crucial balance between H3K36 and H3K27 methylation in maintaining genome instability, offering new therapeutic options for patients with H3K36me-deficient tumors.

Baf-mediated transcriptional regulation of teashirt is essential for the development of neural progenitor cell lineages

Accumulating evidence hints heterochromatin anchoring to the inner nuclear membrane as an upstream regulatory process of gene expression. Given that the formation of neural progenitor cell lineages and the subsequent maintenance of postmitotic neuronal cell identity critically rely on transcriptional regulation, it seems possible that the development of neuronal cells is influenced by cell type-specific and/or context-dependent programmed regulation of heterochromatin anchoring. Here, we explored this possibility by genetically disrupting the evolutionarily conserved barrier-to-autointegration factor (Baf) in the Drosophila nervous system. Through single-cell RNA sequencing, we demonstrated that Baf knockdown induces prominent transcriptomic changes, particularly in type I neuroblasts. Among the differentially expressed genes, our genetic analyses identified teashirt (tsh), a transcription factor that interacts with beta-catenin, to be closely associated with Baf knockdown-induced phenotypes that were suppressed by the overexpression of tsh or beta-catenin. We also found that Baf and tsh colocalized in a region adjacent to heterochromatin in type I NBs. Notably, the subnuclear localization pattern remained unchanged when one of these two proteins was knocked down, indicating that both proteins contribute to the anchoring of heterochromatin to the inner nuclear membrane. Overall, this study reveals that the Baf-mediated transcriptional regulation of teashirt is a novel molecular mechanism that regulates the development of neural progenitor cell lineages.

Glioblastoma evolution and heterogeneity from a 3D whole-tumor perspective

Treatment failure for the lethal brain tumor glioblastoma (GBM) is attributed to intratumoral heterogeneity and tumor evolution. We utilized 3D neuronavigation during surgical resection to acquire samples representing the whole tumor mapped by 3D spatial coordinates. Integrative tissue and single-cell analysis revealed sources of genomic, epigenomic, and microenvironmental intratumoral heterogeneity and their spatial patterning. By distinguishing tumor-wide molecular features from those with regional specificity, we inferred GBM evolutionary trajectories from neurodevelopmental lineage origins and initiating events such as chromothripsis to emergence of genetic subclones and spatially restricted activation of differential tumor and microenvironmental programs in the core, periphery, and contrast-enhancing regions. Our work depicts GBM evolution and heterogeneity from a 3D whole-tumor perspective, highlights potential therapeutic targets that might circumvent heterogeneity-related failures, and establishes an interactive platform enabling 360° visualization and analysis of 3D spatial patterns for user-selected genes, programs, and other features across whole GBM tumors.

DiffDomain enables identification of structurally reorganized topologically associating domains

Topologically associating domains (TADs) are critical structural units in three-dimensional genome organization of mammalian genome. Dynamic reorganizations of TADs between health and disease states are associated with essential genome functions. However, computational methods for identifying reorganized TADs are still in the early stages of development. Here, we present DiffDomain, an algorithm leveraging high-dimensional random matrix theory to identify structurally reorganized TADs using high-throughput chromosome conformation capture (Hi-C) contact maps. Method comparison using multiple real Hi-C datasets reveals that DiffDomain outperforms alternative methods for false positive rates, true positive rates, and identifying a new subtype of reorganized TADs. Applying DiffDomain to Hi-C data from different cell types and disease states demonstrates its biological relevance. Identified reorganized TADs are associated with structural variations and epigenomic changes such as changes in CTCF binding sites. By applying to a single-cell Hi-C data from mouse neuronal development, DiffDomain can identify reorganized TADs between cell types with reasonable reproducibility using pseudo-bulk Hi-C data from as few as 100 cells per condition. Moreover, DiffDomain reveals differential cell-to-population variability and heterogeneous cell-to-cell variability in TADs. Therefore, DiffDomain is a statistically sound method for better comparative analysis of TADs using both Hi-C and single-cell Hi-C data.

Cancer Evolution: A Multifaceted Affair

Cancer cells adapt and survive through the acquisition and selection of molecular modifications. This process defines cancer evolution. Building on a theoretical framework based on heritable genetic changes has provided insights into the mechanisms supporting cancer evolution. However, cancer hallmarks also emerge via heritable nongenetic mechanisms, including epigenetic and chromatin topological changes, and interactions between tumor cells and the tumor microenvironment. Recent findings on tumor evolutionary mechanisms draw a multifaceted picture where heterogeneous forces interact and influence each other while shaping tumor progression. A comprehensive characterization of the cancer evolutionary toolkit is required to improve personalized medicine and biomarker discovery.

SIGNIFICANCE

Tumor evolution is fueled by multiple enabling mechanisms. Importantly, genetic instability, epigenetic reprogramming, and interactions with the tumor microenvironment are neither alternative nor independent evolutionary mechanisms. As demonstrated by findings highlighted in this perspective, experimental and theoretical approaches must account for multiple evolutionary mechanisms and their interactions to ultimately understand, predict, and steer tumor evolution.

The transcriptional activator Klf5 recruits p300-mediated H3K27ac for maintaining trophoblast stem cell pluripotency

The effective proliferation and differentiation of trophoblast stem cells (TSCs) is indispensable for the development of the placenta, which is the key to maintaining normal fetal growth during pregnancy. Kruppel-like factor 5 (Klf5) is implicated in the activation of pluripotency gene expression in embryonic stem cells (ESCs), yet its function in TSCs is poorly understood. Here, we showed that Klf5 knockdown resulted in the downregulation of core TSC-specific genes, consequently causing rapid differentiation of TSCs. Consistently, Klf5-depleted embryos lost the ability to establish TSCs in vitro. At the molecular level, Klf5 preferentially occupied the proximal promoter regions and maintained an open chromatin architecture of key TSC-specific genes. Deprivation of Klf5 impaired the enrichment of p300, a major histone acetyl transferase of H3 lysine 27 acetylation (H3K27ac), and further reduced the occupancy of H3K27ac at promoter regions, leading to decreased transcriptional activity of TSC pluripotency genes. Thus, our findings highlight a novel mechanism of Klf5 in regulating the self-renewal and differentiation of TSCs and provide a reference for understanding placental development and improving pregnancy rates.

Transcriptional reprogramming at the intersection of the heat shock response and proteostasis

Cellular homeostasis is constantly challenged by a myriad of extrinsic and intrinsic stressors. To mitigate the stress-induced damage, cells activate transient survival programs. The heat shock response (HSR) is an evolutionarily well-conserved survival program that is activated in response to proteotoxic stress. The HSR encompasses a dual regulation of transcription, characterized by rapid activation of genes encoding molecular chaperones and concomitant global attenuation of non-chaperone genes. Recent genome-wide approaches have delineated the molecular depth of stress-induced transcriptional reprogramming. The dramatic rewiring of gene and enhancer networks is driven by key transcription factors, including heat shock factors (HSFs), that together with chromatin-modifying enzymes remodel the 3D chromatin architecture, determining the selection of either gene activation or repression. Here, we highlight the current advancements of molecular mechanisms driving transcriptional reprogramming during acute heat stress. We also discuss the emerging implications of HSF-mediated stress signaling in the context of physiological and pathological conditions.

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.

Computationally defined and in vitro validated putative genomic safe harbour loci for transgene expression in human cells

Selection of the target site is an inherent question for any project aiming for directed transgene integration. Genomic safe harbour (GSH) loci have been proposed as safe sites in the human genome for transgene integration. Although several sites have been characterised for transgene integration in the literature, most of these do not meet criteria set out for a GSH and the limited set that do have not been characterised extensively. Here, we conducted a computational analysis using publicly available data to identify 25 unique putative GSH loci that reside in active chromosomal compartments. We validated stable transgene expression and minimal disruption of the native transcriptome in three GSH sites in vitro using human embryonic stem cells (hESCs) and their differentiated progeny. Furthermore, for easy targeted transgene expression, we have engineered constitutive landing pad expression constructs into the three validated GSH in hESCs.

MotifHub: Detection of trans-acting DNA motif group with probabilistic modeling algorithm

BACKGROUND

Trans-acting factors are of special importance in transcription regulation, which is a group of proteins that can directly or indirectly recognize or bind to the 8-12 bp core sequence of cis-acting elements and regulate the transcription efficiency of target genes. The progressive development in high-throughput chromatin capture technology (e.g., Hi-C) enables the identification of chromatin-interacting sequence groups where trans-acting DNA motif groups can be discovered. The problem difficulty lies in the combinatorial nature of DNA sequence pattern matching and its underlying sequence pattern search space.

METHOD

Here, we propose to develop MotifHub for trans-acting DNA motif group discovery on grouped sequences. Specifically, the main approach is to develop probabilistic modeling for accommodating the stochastic nature of DNA motif patterns.

RESULTS

Based on the modeling, we develop global sampling techniques based on EM and Gibbs sampling to address the global optimization challenge for model fitting with latent variables. The results reflect that our proposed approaches demonstrate promising performance with linear time complexities.

CONCLUSION

MotifHub is a novel algorithm considering the identification of both DNA co-binding motif groups and trans-acting TFs. Our study paves the way for identifying hub TFs of stem cell development (OCT4 and SOX2) and determining potential therapeutic targets of prostate cancer (FOXA1 and MYC). To ensure scientific reproducibility and long-term impact, its matrix-algebra-optimized source code is released at http://bioinfo.cs.cityu.edu.hk/MotifHub.

Chromatin remodeling analysis reveals the RdDM pathway responds to low-phosphorus stress in maize

Chromatin in eukaryotes folds into a complex three-dimensional (3D) structure that is essential for controlling gene expression and cellular function and is dynamically regulated in biological processes. Studies on plant phosphorus signaling have concentrated on single genes and gene interactions. It is critical to expand the existing signaling pathway in terms of its 3D structure. In this study, low-Pi treatment led to greater chromatin volume. Furthermore, low-Pi stress increased the insulation score and the number of TAD-like domains, but the effects on the A/B compartment were not obvious. The methylation levels of target sites (hereafter as RdDM levels) peaked at specific TAD-like boundaries, whereas RdDM peak levels at conserved TAD-like boundaries shifted and decreased sharply. The distribution pattern of RdDM sites originating from the Helitron transposons matched that of genome-wide RdDM sites near TAD-like boundaries. RdDM pathway genes were upregulated in the middle or early stages and downregulated in the later stages under low-Pi conditions. The RdDM pathway mutant ddm1a showed increased tolerance to low-Pi stress, with shortened and thickened roots contributing to higher Pi uptake from the shallow soil layer. ChIP-seq results revealed that ZmDDM1A could bind to Pi- and root development-related genes. Strong associations were found between interacting genes in significantly different chromatin-interaction regions and root traits. These findings not only expand the mechanisms by which plants respond to low-Pi stress through the RdDM pathway but also offer a crucial framework for the analysis of biological issues using 3D genomics.

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.

Stage-specific dynamic reorganization of genome topology shapes transcriptional neighborhoods in developing human retinal organoids

We have generated a high-resolution Hi-C map of developing human retinal organoids to elucidate spatiotemporal dynamics of genomic architecture and its relationship with gene expression patterns. We demonstrate progressive stage-specific alterations in DNA topology and correlate these changes with transcription of cell-type-restricted gene markers during retinal differentiation. Temporal Hi-C reveals a shift toward A compartment for protein-coding genes and B compartment for non-coding RNAs, displaying high and low expression, respectively. Notably, retina-enriched genes are clustered near lost boundaries of topologically associated domains (TADs), and higher-order assemblages (i.e., TAD cliques) localize in active chromatin regions with binding sites for eye-field transcription factors. These genes gain chromatin contacts at their transcription start site as organoid differentiation proceeds. Our study provides a global view of chromatin architecture dynamics associated with diversification of cell types during retinal development and serves as a foundational resource for in-depth functional investigations of retinal developmental traits.

3D genome remodeling and homologous pairing during meiotic prophase of mouse oogenesis and spermatogenesis

During meiosis, the chromatin and transcriptome undergo prominent switches. Although recent studies have explored the genome reorganization during spermatogenesis, the chromatin remodeling in oogenesis and characteristics of homologous pairing remain largely elusive. We comprehensively compared chromatin structures and transcriptomes at successive substages of meiotic prophase in both female and male mice using low-input high-through chromosome conformation capture (Hi-C) and RNA sequencing (RNA-seq). Compartments and topologically associating domains (TADs) gradually disappeared and slowly recovered in both sexes. We found that homologs adopted different sex-conserved pairing strategies prior to and after the leptotene-to-zygotene transition, changing from long interspersed nuclear element (LINE)-enriched compartments B to short interspersed nuclear element (SINE)-enriched compartments A. We complemented marker genes and predicted the sex-specific meiotic sterile genes for each substage. This study provides valuable insights into the similarities and distinctions between sexes in chromosome architecture, homologous pairing, and transcriptome during meiotic prophase of both oogenesis and spermatogenesis.

Balanced chromosomal rearrangements implicate YIPF5 and SPATC1L in non-obstructive oligoasthenozoospermia and oligozoospermia and of a derivative chromosome 22 in recurrent miscarriage

Naturally occurring balanced, unbalanced, and complex chromosomal rearrangements have been reported to cause pathogenic genomic or genetic variants leading to infertility and recurrent miscarriage. Therefore, balanced chromosomal rearrangements were used as genomic signposts for identification of candidate genes or genomic loci associated with male infertility due to defects of spermatogenesis, or with recurrent miscarriage. In three male probands, structural chromosomal variants and copy number variants were identified at nucleotide resolution by long-insert genome sequencing approaches and Sanger sequencing. The pathogenic potential of these and affected candidate genes was assessed based on convergent genomic and genotype-phenotype correlation data. Identification of balanced chromosomal rearrangement breakpoints and interpretation in the context of their genomic background of structural and copy number variants led us to conclude that the infertility due to oligoasthenozoospermia and oligozoospermia is most likely associated with a position effect on YIPF5 and SPATC1L, respectively. In a third proband with intellectual disability and recurrent miscarriage, disruption of CAMK2B causing autosomal dominant, intellectual developmental disorder 54 and increased meiotic segregation during gametogenesis of a der(22) are responsible for the reported phenotype. Our data further support the existence of loci at 5q23 and 21q22.3 for these spermatogenesis defects and highlight the importance of the naturally occurring balanced chromosomal rearrangements for assessment of the pathogenic mechanisms. Furthermore, we show comorbidities due to the same balanced chromosomal rearrangement caused by different pathogenic mechanisms.

Maternal dominance contributes to subgenome differentiation in allopolyploid fishes

Teleost fishes, which are the largest and most diverse group of living vertebrates, have a rich history of ancient and recent polyploidy. Previous studies of allotetraploid common carp and goldfish (cyprinids) reported a dominant subgenome, which is more expressed and exhibits biased gene retention. However, the underlying mechanisms contributing to observed 'subgenome dominance' remains poorly understood. Here we report high-quality genomes of twenty-one cyprinids to investigate the origin and subsequent subgenome evolution patterns following three independent allopolyploidy events. We identify the closest extant relatives of the diploid progenitor species, investigate genetic and epigenetic differences among subgenomes, and conclude that observed subgenome dominance patterns are likely due to a combination of maternal dominance and transposable element densities in each polyploid. These findings provide an important foundation to understanding subgenome dominance patterns observed in teleost fishes, and ultimately the role of polyploidy in contributing to evolutionary innovations.

CRISPR-Cas12a-integrated transgenes in genomic safe harbors retain high expression in human hematopoietic iPSC-derived lineages and primary cells

Discovery of genomic safe harbor sites (SHSs) is fundamental for multiple transgene integrations, such as reporter genes, chimeric antigen receptors (CARs), and safety switches, which are required for safe cell products for regenerative cell therapies and immunotherapies. Here we identified and characterized potential SHS in human cells. Using the CRISPR-MAD7 system, we integrated transgenes at these sites in induced pluripotent stem cells (iPSCs), primary T and natural killer (NK) cells, and Jurkat cell line, and demonstrated efficient and stable expression at these loci. Subsequently, we validated the differentiation potential of engineered iPSC toward CD34+ hematopoietic stem and progenitor cells (HSPCs), lymphoid progenitor cells (LPCs), and NK cells and showed that transgene expression was perpetuated in these lineages. Finally, we demonstrated that engineered iPSC-derived NK cells retained expression of a non-virally integrated anti-CD19 CAR, suggesting that several of the investigated SHSs can be used to engineer cells for adoptive immunotherapies.

Reconstructing diploid 3D chromatin structures from single cell Hi-C data with a polymer-based approach

Detailed understanding of the 3D structure of chromatin is a key ingredient to investigate a variety of processes inside the cell. Since direct methods to experimentally ascertain these structures lack the desired spatial fidelity, computational inference methods based on single cell Hi-C data have gained significant interest. Here, we develop a progressive simulation protocol to iteratively improve the resolution of predicted interphase structures by maximum-likelihood association of ambiguous Hi-C contacts using lower-resolution predictions. Compared to state-of-the-art methods, our procedure is not limited to haploid cell data and allows us to reach a resolution of up to 5,000 base pairs per bead. High resolution chromatin models grant access to a multitude of structural phenomena. Exemplarily, we verify the formation of chromosome territories and holes near aggregated chromocenters as well as the inversion of the CpG content for rod photoreceptor cells.

Cell-type-specific 3D-genome organization and transcription regulation in the brain

3D organization of the genome plays a critical role in regulating gene expression. However, it remains unclear how chromatin organization differs among different cell types in the brain. Here we used genome-scale DNA and RNA imaging to investigate 3D-genome organization in transcriptionally distinct cell types in the primary motor cortex of the mouse brain. We uncovered a wide spectrum of differences in the nuclear architecture and 3D-genome organization among different cell types, ranging from the physical size of the cell nucleus to the active-inactive chromatin compartmentalization and radial positioning of chromatin loci within the nucleus. These cell-type-dependent variations in nuclear architecture and chromatin organization exhibited strong correlation with both total transcriptional activity of the cell and transcriptional regulation of cell-type-specific marker genes. Moreover, we found that the methylated-DNA-binding protein MeCP2 regulates transcription in a divergent manner, depending on the nuclear radial positions of chromatin loci, through modulating active-inactive chromatin compartmentalization.

'Enhancing' skeletal muscle and stem cells in three-dimensions: genome regulation of skeletal muscle in development and disease

The noncoding genome imparts important regulatory control over gene expression. In particular, gene enhancers represent a critical layer of control that integrates developmental and differentiation signals outside the cell into transcriptional outputs inside the cell. Recently, there has been an explosion in genomic techniques to probe enhancer control, function, and regulation. How enhancers are regulated and integrate signals in stem cell development and differentiation is largely an open question. In this review, we focus on the role gene enhancers play in muscle stem cell specification, differentiation, and progression. We pay specific attention toward the identification of muscle-specific enhancers, the binding of transcription factors to these enhancers, and how enhancers communicate to their target genes via three-dimensional looping.

Transcription induces context-dependent remodeling of chromatin architecture during differentiation

Metazoan chromosomes are organized into discrete spatial domains (TADs), believed to contribute to the regulation of transcriptional programs. Despite extensive correlation between domain organization and gene activity, a direct mechanistic link is unclear, with perturbation studies often showing little effect. To follow chromatin architecture changes during development, we used Capture Hi-C to interrogate the domains around key differentially expressed genes during mouse thymocyte maturation, uncovering specific remodeling events. Notably, one TAD boundary was broadened to accommodate RNA polymerase elongation past the border, and subdomains were formed around some activated genes without changes in CTCF binding. The ectopic induction of some genes was sufficient to recapitulate domain formation in embryonic stem cells, providing strong evidence that transcription can directly remodel chromatin structure. These results suggest that transcriptional processes drive complex chromosome folding patterns that can be important in certain genomic contexts.

The Information Theory of Aging

Information storage and retrieval is essential for all life. In biology, information is primarily stored in two distinct ways: the genome, comprising nucleic acids, acts as a foundational blueprint and the epigenome, consisting of chemical modifications to DNA and histone proteins, regulates gene expression patterns and endows cells with specific identities and functions. Unlike the stable, digital nature of genetic information, epigenetic information is stored in a digital-analog format, susceptible to alterations induced by diverse environmental signals and cellular damage. The Information Theory of Aging (ITOA) states that the aging process is driven by the progressive loss of youthful epigenetic information, the retrieval of which via epigenetic reprogramming can improve the function of damaged and aged tissues by catalyzing age reversal.

Enhancers within the Ig V Gene Region Orchestrate Chromatin Topology and Regulate V Gene Rearrangement Frequency to Shape the B Cell Receptor Repertoire Specificities

Effective Ab-mediated responses depend on a highly diverse Ab repertoire with the ability to bind a wide range of epitopes in disease-causing agents. The generation of this repertoire depends on the somatic recombination of the variable (V), diversity (D), and joining (J) genes in the Ig loci of developing B cells. It has been known for some time that individual V, D, and J gene segments rearrange at different frequencies, but the mechanisms behind this unequal V gene usage have not been well understood. However, recent work has revealed that newly described enhancers scattered throughout the V gene-containing portion of the Ig loci regulate the V gene recombination frequency in a regional manner. Deletion of three of these enhancers revealed that these elements exert many layers of control during V(D)J recombination, including long-range chromatin interactions, epigenetic milieu, chromatin accessibility, and compartmentalization.

Large-scale data-driven and physics-based models offer insights into the relationships among the structures, dynamics, and functions of chromosomes

The organized three-dimensional chromosome architecture in the cell nucleus provides scaffolding for precise regulation of gene expression. When the cell changes its identity in the cell-fate decision-making process, extensive rearrangements of chromosome structures occur accompanied by large-scale adaptations of gene expression, underscoring the importance of chromosome dynamics in shaping genome function. Over the last two decades, rapid development of experimental methods has provided unprecedented data to characterize the hierarchical structures and dynamic properties of chromosomes. In parallel, these enormous data offer valuable opportunities for developing quantitative computational models. Here, we review a variety of large-scale polymer models developed to investigate the structures and dynamics of chromosomes. Different from the underlying modeling strategies, these approaches can be classified into data-driven ('top-down') and physics-based ('bottom-up') categories. We discuss their contributions to offering valuable insights into the relationships among the structures, dynamics, and functions of chromosomes and propose the perspective of developing data integration approaches from different experimental technologies and multidisciplinary theoretical/simulation methods combined with different modeling strategies.

Cellular reprogramming is driven by widespread rewiring of promoter-enhancer interactions

BACKGROUND

Long-range interactions between promoters and cis-regulatory elements, such as enhancers, play critical roles in gene regulation. However, the role of three-dimensional (3D) chromatin structure in orchestrating changes in transcriptional regulation during direct cell reprogramming is not fully understood.

RESULTS

Here, we performed integrated analyses of chromosomal architecture, epigenetics, and gene expression using Hi-C, promoter Capture Hi-C (PCHi-C), ChIP-seq, and RNA-seq during trans-differentiation of Pre-B cells into macrophages with a β-estradiol inducible C/EBPαER transgene. Within 1h of β-estradiol induction, C/EBPα translocated from the cytoplasm to the nucleus, binding to thousands of promoters and putative regulatory elements, resulting in the downregulation of Pre-B cell-specific genes and induction of macrophage-specific genes. Hi-C results were remarkably consistent throughout trans-differentiation, revealing only a small number of TAD boundary location changes, and A/B compartment switches despite significant changes in the expression of thousands of genes. PCHi-C revealed widespread changes in promoter-anchored loops with decreased interactions in parallel with decreased gene expression, and new and increased promoter-anchored interactions in parallel with increased expression of macrophage-specific genes.

CONCLUSIONS

Overall, our data demonstrate that C/EBPα-induced trans-differentiation involves few changes in genome architecture at the level of TADs and A/B compartments, in contrast with widespread reorganization of thousands of promoter-anchored loops in association with changes in gene expression and cell identity.

Uncovering oligodendrocyte enhancers that control Cnp expression

Oligodendrocytes (OLs) produce myelin sheaths around axons in the central nervous system (CNS). Myelin accelerates the propagation of action potentials along axons and supports the integrity of axons. Impaired myelination has been linked to neurological and neuropsychiatric disorders. As a major component of CNS myelin, 2',3'-cyclic nucleotide 3'-phosphodiesterase (CNP) plays an indispensable role in the axon-supportive function of myelin. Notably, this function requires a high-level expression of CNP in OLs, as evidenced by downregulated expression of CNP in mental disorders and animal models. Little is known about how CNP expression is regulated in OLs. Especially, OL enhancers that govern CNP remain elusive. We have recently developed a powerful method that links OL enhancers to target genes in a principled manner. Here, we applied it to Cnp, uncovering two OL enhancers for it (termed Cnp-E1 and Cnp-E2). Epigenome editing analysis revealed that Cnp-E1 and Cnp-E2 are dedicated to Cnp. ATAC-seq and ChIP-seq data show that Cnp-E1 and Cnp-E2 are conserved OL-specific enhancers. Single cell multi-omics data that jointly profile gene expression and chromatin accessibility suggest that Cnp-E2 plays an important role in Cnp expression in the early stage of OL differentiation while Cnp-E1 sustains it in mature OLs.

The Molecular Impacts of Retrotransposons in Development and Diseases

Retrotransposons are invasive genetic elements that constitute substantial portions of mammalian genomes. They have the potential to influence nearby gene expression through their cis-regulatory sequences, reverse transcription machinery, and the ability to mold higher-order chromatin structures. Due to their multifaceted functions, it is crucial for host fitness to maintain strict regulation of these parasitic sequences to ensure proper growth and development. This review explores how subsets of retrotransposons have undergone evolutionary exaptation to enhance the complexity of mammalian genomes. It also highlights the significance of regulating these elements, drawing on recent studies conducted in human and murine systems.

Correlative single molecule lattice light sheet imaging reveals the dynamic relationship between nucleosomes and the local chromatin environment

In the nucleus, biological processes are driven by proteins that diffuse through and bind to a meshwork of nucleic acid polymers. To better understand this interplay, we developed an imaging platform to simultaneously visualize single protein dynamics together with the local chromatin environment in live cells. Together with super-resolution imaging, new fluorescent probes, and biophysical modeling, we demonstrated that nucleosomes display differential diffusion and packing arrangements as chromatin density increases whereas the viscoelastic properties and accessibility of the interchromatin space remain constant. Perturbing nuclear functions impacted nucleosome diffusive properties in a manner that was dependent on local chromatin density and supportive of a model wherein transcription locally stabilizes nucleosomes while simultaneously allowing for the free exchange of nuclear proteins. Our results reveal that nuclear heterogeneity arises from both active and passive process and highlights the need to account for different organizational principals when modeling different chromatin environments.

Lineage specific 3D genome structure in the adult human brain and neurodevelopmental changes in the chromatin interactome

The human brain is a complex organ comprised of distinct cell types, and the contribution of the 3D genome to lineage specific gene expression remains poorly understood. To decipher cell type specific genome architecture, and characterize fine scale changes in the chromatin interactome across neural development, we compared the 3D genome of the human fetal cortical plate to that of neurons and glia isolated from the adult prefrontal cortex. We found that neurons have weaker genome compartmentalization compared to glia, but stronger TADs, which emerge during fetal development. Furthermore, relative to glia, the neuronal genome shifts more strongly towards repressive compartments. Neurons have differential TAD boundaries that are proximal to active promoters involved in neurodevelopmental processes. CRISPRi on CNTNAP2 in hIPSC-derived neurons reveals that transcriptional inactivation correlates with loss of insulation at the differential boundary. Finally, re-wiring of chromatin loops during neural development is associated with transcriptional and functional changes. Importantly, differential loops in the fetal cortex are associated with autism GWAS loci, suggesting a neuropsychiatric disease mechanism affecting the chromatin interactome. Furthermore, neural development involves gaining enhancer-promoter loops that upregulate genes that control synaptic activity. Altogether, our study provides multi-scale insights on the 3D genome in the human brain.

Statistical learning quantifies transposable element-mediated cis-regulation

BACKGROUND

Transposable elements (TEs) have colonized the genomes of most metazoans, and many TE-embedded sequences function as cis-regulatory elements (CREs) for genes involved in a wide range of biological processes from early embryogenesis to innate immune responses. Because of their repetitive nature, TEs have the potential to form CRE platforms enabling the coordinated and genome-wide regulation of protein-coding genes by only a handful of trans-acting transcription factors (TFs).

RESULTS

Here, we directly test this hypothesis through mathematical modeling and demonstrate that differences in expression at protein-coding genes alone are sufficient to estimate the magnitude and significance of TE-contributed cis-regulatory activities, even in contexts where TE-derived transcription fails to do so. We leverage hundreds of overexpression experiments and estimate that, overall, gene expression is influenced by TE-embedded CREs situated within approximately 500 kb of promoters. Focusing on the cis-regulatory potential of TEs within the gene regulatory network of human embryonic stem cells, we find that pluripotency-specific and evolutionarily young TE subfamilies can be reactivated by TFs involved in post-implantation embryogenesis. Finally, we show that TE subfamilies can be split into truly regulatorily active versus inactive fractions based on additional information such as matched epigenomic data, observing that TF binding may better predict TE cis-regulatory activity than differences in histone marks.

CONCLUSION

Our results suggest that TE-embedded CREs contribute to gene regulation during and beyond gastrulation. On a methodological level, we provide a statistical tool that infers TE-dependent cis-regulation from RNA-seq data alone, thus facilitating the study of TEs in the next-generation sequencing era.

Genomic characterization of IDH-mutant astrocytoma progression to grade 4 in the treatment setting

As the progression of low-grade diffuse astrocytomas into grade 4 tumors significantly impacts patient prognosis, a better understanding of this process is of paramount importance for improved patient care. In this project, we analyzed matched IDH-mutant astrocytomas before and after progression to grade 4 from six patients (discovery cohort) with genome-wide sequencing, 21 additional patients with targeted sequencing, and 33 patients from Glioma Longitudinal AnalySiS cohort for validation. The Cancer Genome Atlas data from 595 diffuse gliomas provided supportive information. All patients in our discovery cohort received radiation, all but one underwent chemotherapy, and no patient received temozolomide (TMZ) before progression to grade 4 disease. One case in the discovery cohort exhibited a hypermutation signature associated with the inactivation of the MSH2 and DNMT3A genes. In other patients, the number of chromosomal rearrangements and deletions increased in grade 4 tumors. The cell cycle checkpoint gene CDKN2A, or less frequently RB1, was most commonly inactivated after receiving both chemo- and radiotherapy when compared to other treatment groups. Concomitant activating PDGFRA/MET alterations were detected in tumors that acquired a homozygous CDKN2A deletion. NRG3 gene was significantly downregulated and recurrently altered in progressed tumors. Its decreased expression was associated with poorer overall survival in both univariate and multivariate analysis. We also detected progression-related alterations in RAD51B and other DNA repair pathway genes associated with the promotion of error-prone DNA repair, potentially facilitating tumor progression. In our retrospective analysis of patient treatment and survival timelines (n = 75), the combination of postoperative radiation and chemotherapy (mainly TMZ) outperformed radiation, especially in the grade 3 tumor cohort, in which it was typically given after primary surgery. Our results provide further insight into the contribution of treatment and genetic alterations in cell cycle, growth factor signaling, and DNA repair-related genes to tumor evolution and progression.

Multi-Scale Imaging of the Dynamic Organization of Chromatin

Chromatin is now regarded as a heterogeneous and dynamic structure occupying a non-random position within the cell nucleus, where it plays a key role in regulating various functions of the genome. This current view of chromatin has emerged thanks to high spatiotemporal resolution imaging, among other new technologies developed in the last decade. In addition to challenging early assumptions of chromatin being regular and static, high spatiotemporal resolution imaging made it possible to visualize and characterize different chromatin structures such as clutches, domains and compartments. More specifically, super-resolution microscopy facilitates the study of different cellular processes at a nucleosome scale, providing a multi-scale view of chromatin behavior within the nucleus in different environments. In this review, we describe recent imaging techniques to study the dynamic organization of chromatin at high spatiotemporal resolution. We also discuss recent findings, elucidated by these techniques, on the chromatin landscape during different cellular processes, with an emphasis on the DNA damage response.

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.

Epstein-Barr virus and host cell 3D genome organization

The human genome is organized in an extremely complexed yet ordered way within the nucleus. Genome organization plays a critical role in the regulation of gene expression. Viruses manipulate the host machinery to influence host genome organization to favor their survival and promote disease development. Epstein-Barr virus (EBV) is a common human virus, whose infection is associated with various diseases, including infectious mononucleosis, cancer, and autoimmune disorders. This review summarizes our current knowledge of how EBV uses different strategies to control the cellular 3D genome organization to affect cell gene expression to transform normal cells into lymphoblasts.

International Union of Basic and Clinical Pharmacology CXIII: Nuclear Receptor Superfamily-Update 2023

The NR superfamily comprises 48 transcription factors in humans that control a plethora of gene network programs involved in a wide range of physiologic processes. This review will summarize and discuss recent progress in NR biology and drug development derived from integrating various approaches, including biophysical techniques, structural studies, and translational investigation. We also highlight how defective NR signaling results in various diseases and disorders and how NRs can be targeted for therapeutic intervention via modulation via binding to synthetic lipophilic ligands. Furthermore, we also review recent studies that improved our understanding of NR structure and signaling. SIGNIFICANCE STATEMENT: Nuclear receptors (NRs) are ligand-regulated transcription factors that are critical regulators of myriad physiological processes. NRs serve as receptors for an array of drugs, and in this review, we provide an update on recent research into the roles of these drug targets.

Gene regulatory network inference in the era of single-cell multi-omics

The interplay between chromatin, transcription factors and genes generates complex regulatory circuits that can be represented as gene regulatory networks (GRNs). The study of GRNs is useful to understand how cellular identity is established, maintained and disrupted in disease. GRNs can be inferred from experimental data - historically, bulk omics data - and/or from the literature. The advent of single-cell multi-omics technologies has led to the development of novel computational methods that leverage genomic, transcriptomic and chromatin accessibility information to infer GRNs at an unprecedented resolution. Here, we review the key principles of inferring GRNs that encompass transcription factor-gene interactions from transcriptomics and chromatin accessibility data. We focus on the comparison and classification of methods that use single-cell multimodal data. We highlight challenges in GRN inference, in particular with respect to benchmarking, and potential further developments using additional data modalities.

Active maintenance of CD8+ T cell naivety through regulation of global genome architecture

The differentiation of naive CD8+ T lymphocytes into cytotoxic effector and memory CTL results in large-scale changes in transcriptional and phenotypic profiles. Little is known about how large-scale changes in genome organization underpin these transcriptional programs. We use Hi-C to map changes in the spatial organization of long-range genome contacts within naive, effector, and memory virus-specific CD8+ T cells. We observe that the architecture of the naive CD8+ T cell genome is distinct from effector and memory genome configurations, with extensive changes within discrete functional chromatin domains associated with effector/memory differentiation. Deletion of BACH2, or to a lesser extent, reducing SATB1 DNA binding, within naive CD8+ T cells results in a chromatin architecture more reminiscent of effector/memory states. This suggests that key transcription factors within naive CD8+ T cells act to restrain T cell differentiation by actively enforcing a unique naive chromatin state.

Genome-wide transcription factor-binding maps reveal cell-specific changes in the regulatory architecture of human HSPCs

Hematopoietic stem and progenitor cells (HSPCs) rely on a complex interplay among transcription factors (TFs) to regulate differentiation into mature blood cells. A heptad of TFs (FLI1, ERG, GATA2, RUNX1, TAL1, LYL1, LMO2) bind regulatory elements in bulk CD34+ HSPCs. However, whether specific heptad-TF combinations have distinct roles in regulating hematopoietic differentiation remains unknown. We mapped genome-wide chromatin contacts (HiC, H3K27ac, HiChIP), chromatin modifications (H3K4me3, H3K27ac, H3K27me3) and 10 TF binding profiles (heptad, PU.1, CTCF, STAG2) in HSPC subsets (stem/multipotent progenitors plus common myeloid, granulocyte macrophage, and megakaryocyte erythrocyte progenitors) and found TF occupancy and enhancer-promoter interactions varied significantly across cell types and were associated with cell-type-specific gene expression. Distinct regulatory elements were enriched with specific heptad-TF combinations, including stem-cell-specific elements with ERG, and myeloid- and erythroid-specific elements with combinations of FLI1, RUNX1, GATA2, TAL1, LYL1, and LMO2. Furthermore, heptad-occupied regions in HSPCs were subsequently bound by lineage-defining TFs, including PU.1 and GATA1, suggesting that heptad factors may prime regulatory elements for use in mature cell types. We also found that enhancers with cell-type-specific heptad occupancy shared a common grammar with respect to TF binding motifs, suggesting that combinatorial binding of TF complexes was at least partially regulated by features encoded in DNA sequence motifs. Taken together, this study comprehensively characterizes the gene regulatory landscape in rare subpopulations of human HSPCs. The accompanying data sets should serve as a valuable resource for understanding adult hematopoiesis and a framework for analyzing aberrant regulatory networks in leukemic cells.

Single-cell DNA methylation and 3D genome architecture in the human brain

Delineating the gene-regulatory programs underlying complex cell types is fundamental for understanding brain function in health and disease. Here, we comprehensively examined human brain cell epigenomes by probing DNA methylation and chromatin conformation at single-cell resolution in 517 thousand cells (399 thousand neurons and 118 thousand non-neurons) from 46 regions of three adult male brains. We identified 188 cell types and characterized their molecular signatures. Integrative analyses revealed concordant changes in DNA methylation, chromatin accessibility, chromatin organization, and gene expression across cell types, cortical areas, and basal ganglia structures. We further developed single-cell methylation barcodes that reliably predict brain cell types using the methylation status of select genomic sites. This multimodal epigenomic brain cell atlas provides new insights into the complexity of cell-type-specific gene regulation in adult human brains.

Identifying quantitatively differential chromosomal compartmentalization changes and their biological significance from Hi-C data using DARIC

BACKGROUND

Chromosomal compartmentalization plays a critical role in maintaining proper transcriptional programs in cell differentiation and oncogenesis. However, currently the prevalent method for comparative analysis of compartmentalization landscapes between different cell types is limited to the qualitative switched compartments.

RESULTS

To identify genomic regions with quantitatively differential compartmentalization changes from genome-wide chromatin conformation data like Hi-C, we developed a computational framework named DARIC. DARIC includes three modules: compartmentalization quantification, normalization, and differential analysis. Comparing DARIC with the conventional compartment switching analysis reveals substantial regions characterized by quantitatively significant compartmentalization changes without switching. These changes are accompanied by changes in gene expression, chromatin accessibility, H3K27ac intensity, as well as the interactions with nuclear lamina proteins and nuclear positioning, highlighting the functional importance of such quantitative changes in gene regulation. We applied DARIC to dissect the quantitative compartmentalization changes during human cardiomyocyte differentiation and identified two distinct mechanisms for gene activation based on the association with compartmentalization changes. Using the quantitative compartmentalization measurement module from DARIC, we further dissected the compartment variability landscape in the human genome by analyzing a compendium of 32 Hi-C datasets from 4DN. We discovered an interesting correlation between compartmentalization variability and sub-compartments.

CONCLUSIONS

DARIC is a useful tool for analyzing quantitative compartmentalization changes and mining novel biological insights from increasing Hi-C data. Our results demonstrate the functional significance of quantitative compartmentalization changes in gene regulation, and provide new insights into the relationship between compartmentalization variability and sub-compartments in the human genome.

Functional characterization of Alzheimer's disease genetic variants in microglia

Candidate cis-regulatory elements (cCREs) in microglia demonstrate the most substantial enrichment for Alzheimer's disease (AD) heritability compared to other brain cell types. However, whether and how these genome-wide association studies (GWAS) variants contribute to AD remain elusive. Here we prioritize 308 previously unreported AD risk variants at 181 cCREs by integrating genetic information with microglia-specific 3D epigenome annotation. We further establish the link between functional variants and target genes by single-cell CRISPRi screening in microglia. In addition, we show that AD variants exhibit allelic imbalance on target gene expression. In particular, rs7922621 is the effective variant in controlling TSPAN14 expression among other nominated variants in the same cCRE and exerts multiple physiological effects including reduced cell surface ADAM10 and altered soluble TREM2 (sTREM2) shedding. Our work represents a systematic approach to prioritize and characterize AD-associated variants and provides a roadmap for advancing genetic association to experimentally validated cell-type-specific phenotypes and mechanisms.

The dynamic role of cohesin in maintaining human genome architecture

Recent advances in genomic and imaging techniques have revealed the complex manner of organizing billions of base pairs of DNA necessary for maintaining their functionality and ensuring the proper expression of genetic information. The SMC proteins and cohesin complex primarily contribute to forming higher-order chromatin structures, such as chromosomal territories, compartments, topologically associating domains (TADs) and chromatin loops anchored by CCCTC-binding factor (CTCF) protein or other genome organizers. Cohesin plays a fundamental role in chromatin organization, gene expression and regulation. This review aims to describe the current understanding of the dynamic nature of the cohesin-DNA complex and its dependence on cohesin for genome maintenance. We discuss the current 3C technique and numerous bioinformatics pipelines used to comprehend structural genomics and epigenetics focusing on the analysis of Cohesin-centred interactions. We also incorporate our present comprehension of Loop Extrusion (LE) and insights from stochastic modelling.

Rare variation in noncoding regions with evolutionary signatures contributes to autism spectrum disorder risk

Little is known about the role of noncoding regions in the etiology of autism spectrum disorder (ASD). We examined three classes of noncoding regions: Human Accelerated Regions (HARs), which show signatures of positive selection in humans; experimentally validated neural Vista Enhancers (VEs); and conserved regions predicted to act as neural enhancers (CNEs). Targeted and whole genome analysis of >16,600 samples and >4900 ASD probands revealed that likely recessive, rare, inherited variants in HARs, VEs, and CNEs substantially contribute to ASD risk in probands whose parents share ancestry, which enriches for recessive contributions, but modestly, if at all, in simplex family structures. We identified multiple patient variants in HARs near IL1RAPL1 and in a VE near SIM1 and showed that they change enhancer activity. Our results implicate both human-evolved and evolutionarily conserved noncoding regions in ASD risk and suggest potential mechanisms of how changes in regulatory regions can modulate social behavior.

Context-dependent perturbations in chromatin folding and the transcriptome by cohesin and related factors

Cohesin regulates gene expression through context-specific chromatin folding mechanisms such as enhancer-promoter looping and topologically associating domain (TAD) formation by cooperating with factors such as cohesin loaders and the insulation factor CTCF. We developed a computational workflow to explore how three-dimensional (3D) structure and gene expression are regulated collectively or individually by cohesin and related factors. The main component is CustardPy, by which multi-omics datasets are compared systematically. To validate our methodology, we generated 3D genome, transcriptome, and epigenome data before and after depletion of cohesin and related factors and compared the effects of depletion. We observed diverse effects on the 3D genome and transcriptome, and gene expression changes were correlated with the splitting of TADs caused by cohesin loss. We also observed variations in long-range interactions across TADs, which correlated with their epigenomic states. These computational tools and datasets will be valuable for 3D genome and epigenome studies.

ARID4B loss leads to activated STAT1-dependent interferon pathway in mouse embryonic stem cells and during meso/endodermal differentiation

Objective

Proper deactivation of the pluripotency network and activation of a lineage-specific gene expression program are critical for mouse embryonic stem cell (mESC) differentiation. This is achieved by the coordinated action of transcription and chromatin factors. Our previous work identified ARID4B as a critical chromatin factor for mesoderm and endoderm differentiation. As part of a histone deacetylase complex, ARID4B plays a role in transcriptional suppression of its direct targets. Here, we investigated the mechanism of ARID4B function in mESC differentiation by focusing on genes and pathways that are upregulated in its absence.

Material and Methods

We analyzed transcriptomic results of wild-type and arid4bΔ endoderm or mesoderm differentiated cells through integrative genomics viewer and ingenuity pathway analysis. We performed real-time quantitative polymerase chain reaction for selected genes. To understand pathway activation, we performed Western blot for candidate proteins during the time-course of differentiation. We also analyzed H3K4me3, H3K27me3 and H3K27Ac ChIP-seq results to understand changes in the chromatin environment.

Results

Interferon-related genes were activated in arid4bΔ mESCs and endoderm or mesoderm directed cells. Consistent with this, higher phosphorylated STAT1 levels were found in arid4bΔ mESCs while a related phosphorylated STAT3 was unchanged. Finally, we observed a significant increase in H3K4me3 around interferon-related distal gene regulatory regions with a combination of either upregulation of H3K27Ac level or downregulation of H3K27me3 level.

Conclusion

These results provide evidence that ARID4B is involved in the suppression of interferon-related genes in mESCs and during meso/endoderm differentiation through modulation, mainly of H3K4me3. This regulation might be important for successful mESC differentiation.

The 3D chromatin landscape of rhabdomyosarcoma

Rhabdomyosarcoma (RMS) is a pediatric soft tissue cancer with a lack of precision therapy options for patients. We hypothesized that with a general paucity of known mutations in RMS, chromatin structural driving mechanisms are essential for tumor proliferation. Thus, we carried out high-depth in situ Hi-C in representative cell lines and patient-derived xenografts (PDXs) to define chromatin architecture in each major RMS subtype. We report a comprehensive 3D chromatin structural analysis and characterization of fusion-positive (FP-RMS) and fusion-negative RMS (FN-RMS). We have generated spike-in in situ Hi-C chromatin interaction maps for the most common FP-RMS and FN-RMS cell lines and compared our data with PDX models. In our studies, we uncover common and distinct structural elements in large Mb-scale chromatin compartments, tumor-essential genes within variable topologically associating domains and unique patterns of structural variation. Our high-depth chromatin interactivity maps and comprehensive analyses provide context for gene regulatory events and reveal functional chromatin domains in RMS.

Simultaneous profiling of chromatin architecture and transcription in single cells

The three-dimensional structure of chromatin plays a crucial role in development and disease, both of which are associated with transcriptional changes. However, given the heterogeneity in single-cell chromatin architecture and transcription, the regulatory relationship between the three-dimensional chromatin structure and gene expression is difficult to explain based on bulk cell populations. Here we develop a single-cell, multimodal, omics method allowing the simultaneous detection of chromatin architecture and messenger RNA expression by sequencing (single-cell transcriptome sequencing (scCARE-seq)). Applying scCARE-seq to examine chromatin architecture and transcription from 2i to serum single mouse embryonic stem cells, we observe improved separation of cell clusters compared with single-cell chromatin conformation capture. In addition, after defining the cell-cycle phase of each cell through chromatin architecture extracted by scCARE-seq, we find that periodic changes in chromatin architecture occur in parallel with transcription during the cell cycle. These findings highlight the potential of scCARE-seq to facilitate comprehensive analyses that may boost our understanding of chromatin architecture and transcription in the same single cell.

Chromosome-level organization of the regulatory genome in the Drosophila nervous system

Previous studies have identified topologically associating domains (TADs) as basic units of genome organization. We present evidence of a previously unreported level of genome folding, where distant TAD pairs, megabases apart, interact to form meta-domains. Within meta-domains, gene promoters and structural intergenic elements present in distant TADs are specifically paired. The associated genes encode neuronal determinants, including those engaged in axonal guidance and adhesion. These long-range associations occur in a large fraction of neurons but support transcription in only a subset of neurons. Meta-domains are formed by diverse transcription factors that are able to pair over long and flexible distances. We present evidence that two such factors, GAF and CTCF, play direct roles in this process. The relative simplicity of higher-order meta-domain interactions in Drosophila, compared with those previously described in mammals, allowed the demonstration that genomes can fold into highly specialized cell-type-specific scaffolds that enable megabase-scale regulatory associations.

Extravillous trophoblast cell lineage development is associated with active remodeling of the chromatin landscape

The extravillous trophoblast cell lineage is a key feature of placentation and successful pregnancy. Knowledge of transcriptional regulation driving extravillous trophoblast cell development is limited. Here, we map the transcriptome and epigenome landscape as well as chromatin interactions of human trophoblast stem cells and their transition into extravillous trophoblast cells. We show that integrating chromatin accessibility, long-range chromatin interactions, transcriptomic, and transcription factor binding motif enrichment enables identification of transcription factors and regulatory mechanisms critical for extravillous trophoblast cell development. We elucidate functional roles for TFAP2C, SNAI1, and EPAS1 in the regulation of extravillous trophoblast cell development. EPAS1 is identified as an upstream regulator of key extravillous trophoblast cell transcription factors, including ASCL2 and SNAI1 and together with its target genes, is linked to pregnancy loss and birth weight. Collectively, we reveal activation of a dynamic regulatory network and provide a framework for understanding extravillous trophoblast cell specification in trophoblast cell lineage development and human placentation.