Mechanisms of Interplay between Transcription Factors and the 3D Genome

Identifying the DNA methylation preference of transcription factors using ProtBERT and SVM

Transcription factors (TFs) can affect gene expression by binding to certain specific DNA sequences. This binding process of TFs may be modulated by DNA methylation. A subset of TFs that serve as methylation readers preferentially binds to certain methylated DNA and is defined as TFPM. The identification of TFPMs enhances our understanding of DNA methylation's role in gene regulation. However, their experimental identification is resource-demanding. In this study, we propose a novel two-step computational approach to classify TFs and TFPMs. First, we employed a fine-tuned ProtBERT model to differentiate between the classes of TFs and non-TFs. Second, we combined the Reduced Amino Acid Category (RAAC) with K-mer and SVM to predict the potential of TFs to bind to methylated DNA. Comparative experiments demonstrate that our proposed methods outperform all existing approaches and emphasize the efficiency of our computational framework in classifying TFs and TFPMs. Cross-species validation on an independent mouse dataset further demonstrates the generalizability of our proposed framework In addition, we conducted predictions on all human transcription factors and found that most of the top 20 proteins belong to the Krueppel C2H2-type Zinc-finger family. So far, some studies have demonstrated a partial correlation between this family and DNA methylation and confirmed the preference of some of its members, thereby showing the robustness of our approach.

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

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

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.

Regulation of 3D genome organization during T cell activation

Within the three-dimensional (3D) nuclear space, the genome organizes into a series of orderly structures that impose important influences on gene regulation. T lymphocytes, crucial players in adaptive immune responses, undergo intricate transcriptional remodeling upon activation, leading to differentiation into specific effector and memory T cell subsets. Recent evidence suggests that T cell activation is accompanied by dynamic changes in genome architecture at multiple levels, providing a unique biological context to explore the functional relevance and molecular mechanisms of 3D genome organization. Here, we summarize recent advances that link the reorganization of genome architecture to the remodeling of transcriptional programs and conversion of cell fates during T cell activation and differentiation. We further discuss how various chromatin architecture regulators, including CCCTC-binding factor and several transcription factors, collectively modulate the genome architecture during this process.

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.

Maf-family bZIP transcription factor NRL interacts with RNA-binding proteins and R-loops in retinal photoreceptors

RNA-binding proteins (RBPs) perform diverse functions including the regulation of chromatin dynamics and the coupling of transcription with RNA processing. However, our understanding of their actions in mammalian neurons remains limited. Using affinity purification, yeast-two-hybrid and proximity ligation assays, we identified interactions of multiple RBPs with neural retina leucine (NRL) zipper, a Maf-family transcription factor critical for retinal rod photoreceptor development and function. In addition to splicing, many NRL-interacting RBPs are associated with R-loops, which form during transcription and increase during photoreceptor maturation. Focusing on DHX9 RNA helicase, we demonstrate that its expression is modulated by NRL and that the NRL-DHX9 interaction is positively influenced by R-loops. ssDRIP-Seq analysis reveals both stranded and unstranded R-loops at distinct genomic elements, characterized by active and inactive epigenetic signatures and enriched at neuronal genes. NRL binds to both types of R-loops, suggesting an epigenetically independent function. Our findings suggest additional functions of NRL during transcription and highlight complex interactions among transcription factors, RBPs, and R-loops in regulating photoreceptor gene expression in the mammalian retina.

Molecular mechanisms altering cell identity in cancer

Intrinsic and extrinsic factors influence cancer cell identity throughout its lifespan. During tumor progression and metastasis formation, cancer cells are exposed to different environmental stimuli, resulting in a stepwise cellular reprogramming. Similar stepwise changes of cell identity have been shown as a major consequence of cancer treatment, as cells are exposed to extracellular stress that can result in the establishment of subpopulations exhibiting different epigenetic and transcriptional patterns, indicating a rapid adaptation mechanism of cellular identity by extrinsic stress factors. Both mechanisms, tumor progression-mediated changes and therapy response, rely on signaling pathways affecting the epigenetic and subsequent transcriptional landscape, which equip the cells with mechanisms for survival and tumor progression. These non-genetic alterations are propagated to the daughter cells, indicating a need for successful information propagation and transfer to the daughter generations, thereby allowing for a stepwise adaptation to environmental cues. However, the exact mechanisms how these cell identity changes are occurring, which context-specific mechanisms are behind and how this can be exploited for future therapeutic interventions is not yet fully understood and exploited. In this review, we discuss the current knowledge on cell identity maintenance mechanisms intra- and intergenerational in development and disease and how these mechanisms are altered in cancer. We will as well address how cancer treatment might target these properties.

Pioneer factor GATA6 promotes colorectal cancer through 3D genome regulation

Colorectal cancer (CRC) is one of the most lethal and prevalent malignancies. While the overexpression of pioneer factor GATA6 in CRC has been linked with metastasis, its role in genome-wide gene expression dysregulation remains unclear. Through studies of primary human CRC tissues and analysis of the TCGA data, we found that GATA6 preferentially binds at CRC-specific active enhancers, with enrichment at enhancer-promoter loop anchors. GATA6 protein also physically interacts with CTCF, suggesting its critical role in 3D genome organization. The ablation of GATA6 through AID and CRISPR systems severely impaired cancer cell clonogenicity and proliferation. Mechanistically, GATA6 knockout induced global loss of CRC-specific open chromatins and extensive alterations of critical enhancer-promoter interactions for CRC oncogenes. Last, we showed that GATA6 knockout greatly reduced tumor growth and improved survival in mice. Together, we revealed a previously unidentified mechanism by which GATA6 contributes to the pathogenesis of colorectal cancer.

A guide to studying 3D genome structure and dynamics in the kidney

The human genome is tightly packed into the 3D environment of the cell nucleus. Rapidly evolving and sophisticated methods of mapping 3D genome architecture have shed light on fundamental principles of genome organization and gene regulation. The genome is physically organized on different scales, from individual genes to entire chromosomes. Nuclear landmarks such as the nuclear envelope and nucleoli have important roles in compartmentalizing the genome within the nucleus. Genome activity (for example, gene transcription) is also functionally partitioned within this 3D organization. Rather than being static, the 3D organization of the genome is tightly regulated over various time scales. These dynamic changes in genome structure over time represent the fourth dimension of the genome. Innovative methods have been used to map the dynamic regulation of genome structure during important cellular processes including organism development, responses to stimuli, cell division and senescence. Furthermore, disruptions to the 4D genome have been linked to various diseases, including of the kidney. As tools and approaches to studying the 4D genome become more readily available, future studies that apply these methods to study kidney biology will provide insights into kidney function in health and disease.

RUNX2 Phase Separation Mediates Long-Range Regulation Between Osteoporosis-Susceptibility Variant and XCR1 to Promote Osteoblast Differentiation

GWASs have identified many loci associated with osteoporosis, but the underlying genetic regulatory mechanisms and the potential drug target need to be explored. Here, a new regulatory mechanism is found that a GWAS intergenic SNP (rs4683184) functions as an enhancer to influence the binding affinity of transcription factor RUNX2, whose phase separation can mediate the long-range chromatin interaction between enhancer and target gene XCR1 (a member of the GPCR family), leading to changes of XCR1 expression and osteoblast differentiation. Bone-targeting AAV of Xcr1 can improve bone formation in osteoporosis mice, suggesting that XCR1 can be a new susceptibility gene for osteoporosis. This study is the first to link non-coding SNP with phase separation, providing a new insight into long-range chromatin regulation mechanisms with susceptibility to complex diseases, and finding a potential target for the development of osteoporosis drugs and corresponding translational research.

Interphase chromatin biophysics and mechanics: new perspectives and open questions

The physical organization and properties of chromatin within the interphase nucleus are intimately linked to a wide range of functional DNA-based processes. In this context, interphase chromatin mechanics - that is, how chromatin, physically, responds to forces - is gaining increasing attention. Recent methodological advances for probing the force-response of chromatin in cellulo open new avenues for research, as well as new questions. This review discusses emerging views from these approaches and others, including recent in vitro single-molecule studies of cohesin and condensin motor activities, providing insights into physical and material aspects of chromatin, its plasticity in the context of functional processes, its nonequilibrium or 'active matter' properties, and the importance of factors such as chromatin fiber tension and stiffness. This growing field offers exciting opportunities to better understand the interplay between interphase chromosome structure, dynamics, mechanics, and functions.

MYBL2 promotes proliferation of clear cell renal cell carcinoma by regulating TOP2A and activating AKT/mTOR signaling pathway

Renal cell carcinoma (RCC) is one of the most common malignancies in the urinary system, and clear cell renal cell carcinoma (ccRCC) is the most common subtype. MYBL2 has been reported to be overexpressed in various tumors and associated with poor prognosis in patients, but its biological role in ccRCC remains unclear. In this study, we investigated the mRNA and protein expression levels of MYBL2 in ccRCC samples and evaluated the prognostic value of MYBL2 using TCGA dataset. In vitro functional assays were performed using CCK-8, EdU, colony formation, cell scratch, and transwell assays, as well as in vivo tumorigenesis assays to investigate the biological functions of MYBL2 in ccRCC. Additionally, gene set enrichment analysis (GSEA) was used to explore the downstream pathways of MYBL2, which were further validated. Finally, we predicted the target genes of MYBL2 using bioinformatics and validated them using ChIP and dual-luciferase reporter gene assays. MYBL2 expression was significantly higher in ccRCC than in adjacent normal tissues and was associated with poor prognosis. MYBL2 expression was positively correlated with the pathological tumor grade and clinical TNM stage of ccRCC patients. Knockdown of MYBL2 significantly inhibited the proliferation of renal cancer cells in vitro and in vivo, and knockdown of MYBL2 could inhibit cell invasion and migration, while overexpression of MYBL2 had the opposite effect. GSEA revealed that MYBL2 was associated with the mTOR signaling pathway and cell cycle pathway, which was confirmed by our study. Finally, we found that TOP2A was a target gene of MYBL2, and MYBL2 could bind to the TOP2A promoter to regulate its transcriptional activity, promoting the proliferation of clear cell renal cell carcinoma cells. MYBL2 emerges as a highly expressed factor that significantly correlates with adverse patient prognosis in ccRCC. Mechanistically, MYBL2 transcriptionally upregulates TOP2A, thereby modulating the proliferation of ccRCC cells. Furthermore, MYBL2 activates the mTOR signaling pathway, a critical node in the progression of ccRCC. Collectively, these findings position MYBL2 as a promising candidate for both a biological marker and a therapeutic target in the management of ccRCC.

Dose-dependent sensitivity of human 3D chromatin to a heart disease-linked transcription factor

Dosage-sensitive transcription factors (TFs) underlie altered gene regulation in human developmental disorders, and cell-type specific gene regulation is linked to the reorganization of 3D chromatin during cellular differentiation. Here, we show dose-dependent regulation of chromatin organization by the congenital heart disease (CHD)-linked, lineage-restricted TF TBX5 in human cardiomyocyte differentiation. Genome organization, including compartments, topologically associated domains, and chromatin loops, are sensitive to reduced TBX5 dosage in a human model of CHD, with variations in response across individual cells. Regions normally bound by TBX5 are especially sensitive, while co-occupancy with CTCF partially protects TBX5-bound TAD boundaries and loop anchors. These results highlight the importance of lineage-restricted TF dosage in cell-type specific 3D chromatin dynamics, suggesting a new mechanism for TF-dependent disease.

Maf-family bZIP transcription factor NRL interacts with RNA-binding proteins and R-loops in retinal photoreceptors

RNA-binding proteins (RBPs) perform diverse functions including the regulation of chromatin dynamics and the coupling of transcription with RNA processing. However, our understanding of their actions in mammalian neurons remains limited. Using affinity purification, yeast-two-hybrid and proximity ligation assays, we identified interactions of multiple RBPs with NRL, a Maf-family bZIP transcription factor critical for retinal rod photoreceptor development and function. In addition to splicing, many NRL-interacting RBPs are associated with R-loops, which form during transcription and increase during photoreceptor maturation. Focusing on DHX9 RNA helicase, we demonstrate that its expression is modulated by NRL and that the NRL-DHX9 interaction is positively influenced by R-loops. ssDRIP-Seq analysis reveals both stranded and unstranded R-loops at distinct genomic elements, characterized by active and inactive epigenetic signatures and enriched at neuronal genes. NRL binds to both types of R-loops, suggesting an epigenetically independent function. Our findings suggest additional functions of NRL during transcription and highlight complex interactions among transcription factors, RBPs, and R-loops in regulating photoreceptor gene expression in the mammalian retina.

ZNF143 is a transcriptional regulator of nuclear-encoded mitochondrial genes that acts independently of looping and CTCF

Gene expression is orchestrated by transcription factors, which function within the context of a three-dimensional genome. Zinc-finger protein 143 (ZNF143/ZFP143) is a transcription factor that has been implicated in both gene activation and chromatin looping. To study the direct consequences of ZNF143/ZFP143 loss, we generated a ZNF143/ZFP143 depletion system in mouse embryonic stem cells. Our results show that ZNF143/ZFP143 degradation has no effect on chromatin looping. Systematic analysis of ZNF143/ZFP143 occupancy data revealed that a commonly used antibody cross-reacts with CTCF, leading to its incorrect association with chromatin loops. Nevertheless, ZNF143/ZFP143 specifically activates nuclear-encoded mitochondrial genes, and its loss leads to severe mitochondrial dysfunction. Using an in vitro embryo model, we find that ZNF143/ZFP143 is an essential regulator of organismal development. Our results establish ZNF143/ZFP143 as a conserved transcriptional regulator of cell proliferation and differentiation by safeguarding mitochondrial activity.

Dynamic 3D genome reorganization during senescence: defining cell states through chromatin

Cellular senescence, a cell state characterized by growth arrest and insensitivity to growth stimulatory hormones, is accompanied by a massive change in chromatin organization. Senescence can be induced by a range of physiological signals and pathological stresses and was originally thought to be an irreversible state, implicated in normal development, wound healing, tumor suppression and aging. Recently cellular senescence was shown to be reversible in some cases, with exit being triggered by the modulation of the cell's transcriptional program by the four Yamanaka factors, the suppression of p53 or H3K9me3, PDK1, and/or depletion of AP-1. Coincident with senescence reversal are changes in chromatin organization, most notably the loss of senescence-associated heterochromatin foci (SAHF) found in oncogene-induced senescence. In addition to fixed-cell imaging, chromatin conformation capture and multi-omics have been used to examine chromatin reorganization at different spatial resolutions during senescence. They identify determinants of SAHF formation and other key features that differentiate distinct types of senescence. Not surprisingly, multiple factors, including the time of induction, the type of stress experienced, and the type of cell involved, influence the global reorganization of chromatin in senescence. Here we discuss how changes in the three-dimensional organization of the genome contribute to the regulation of transcription at different stages of senescence. In particular, the distinct contributions of heterochromatin- and lamina-mediated interactions, changes in gene expression, and other cellular control mechanisms are discussed. We propose that high-resolution temporal and spatial analyses of the chromatin landscape during senescence will identify early markers of the different senescence states to help guide clinical diagnosis.

Epigenetic remodeling and 3D chromatin reorganization governed by NKX2-1 drive neuroendocrine prostate cancer

A significant number of castration-resistant prostate cancer (CRPC) evolve into a neuroendocrine (NE) subtype termed NEPC, leading to resistance to androgen receptor (AR) pathway inhibitors and poor clinical outcomes. Through Hi-C analyses of a panel of patient-derived xenograft tumors, here we report drastically different 3D chromatin architectures between NEPC and CRPC samples. Such chromatin re-organization was faithfully recapitulated in vitro on isogenic cells undergoing NE transformation (NET). Mechanistically, neural transcription factor (TF) NKX2-1 is selectively and highly expressed in NEPC tumors and is indispensable for NET across various models. NKX2-1 preferentially binds to gene promoters, but it interacts with chromatin-pioneering factors such as FOXA2 at enhancer elements through chromatin looping, further strengthening FOXA2 binding at NE enhancers. Conversely, FOXA2 mediates regional DNA demethylation, attributing to NE enhancer priming and inducing NKX2-1 expression, forming a feed-forward loop. Single-cell multiome analyses of isogenic cells over time-course NET cells identify individual cells amid luminal-to-NE transformation, exhibiting intermediate epigenetic and transcriptome states. Lastly, NKX2-1/FOXA2 interacts with, and recruits CBP/p300 proteins to activate NE enhancers, and pharmacological inhibitors of CBP/p300 effectively blunted NE gene expression and abolished NEPC tumor growth. Thus, our study reports a hierarchical network of TFs governed by NKX2-1 in regulating the 2D and 3D chromatin re-organization during NET and uncovers a promising therapeutic approach to eradicate NEPC.

Elevated EGR1 Binding at Enhancers in Excitatory Neurons Correlates with Neuronal Subtype-Specific Epigenetic Regulation

Brain development and neuronal cell specification are accompanied with epigenetic changes to achieve diverse gene expression regulation. Interacting with cell-type specific epigenetic marks, transcription factors bind to different sets of cis-regulatory elements in different types of cells. Currently, it remains largely unclear how cell-type specific gene regulation is achieved for neurons. In this study, we generated epigenetic maps to perform comparative histone modification analysis between excitatory and inhibitory neurons. We found that neuronal cell-type specific histone modifications are enriched in super enhancer regions containing abundant EGR1 motifs. Further CUT&RUN data validated that more EGR1 binding sites can be detected in excitatory neurons and primarily located in enhancers. Integrative analysis revealed that EGR1 binding is strongly correlated with various epigenetic markers for open chromatin regions and associated with distinct gene pathways with neuronal subtype-specific functions. In inhibitory neurons, the majority of genomic regions hosting EGR1 binding sites become accessible at early embryonic stages. In contrast, the super enhancers in excitatory neurons hosting EGR1 binding sites gained their accessibility during postnatal stages. This study highlights the significance of transcription factor binding to enhancer regions, which may play a crucial role in establishing cell-type specific gene regulation in neurons.

Multiscale 3D genome rewiring during PTF1A-mediated somatic cell reprogramming into neural stem cells

The genome is intricately folded into chromatin compartments, topologically associating domains (TADs) and loops unique to each cell type. How this higher-order genome organization regulates cell fate transition remains elusive. Here we show how a single non-neural progenitor transcription factor, PTF1A, reorchestrates the 3D genome during fibroblast transdifferentiation into neural stem cells (NSCs). Multiomics analyses integrating Hi-C data, PTF1A and CTCF DNA-binding profiles, H3K27ac modification, and gene expression, demonstrate that PTF1A binds to subTAD boundaries subsequently associated with elevated CTCF binding and enhanced boundary insulation, and reorganizes chromatin loops, leading to gene expression changes that drive transdifferentiation into NSCs. Moreover, PTF1A activates enhancers and super-enhancers near low-insulation boundaries and modulates H3K27ac deposition, promoting cell fate transitions. Together, our data implicate an involvement of 3D genome in transcriptional and cell fate alterations, and highlight an essential role for PTF1A in gene expression control and multiscale 3D genome remodeling during cell reprogramming.

Molecular mechanisms of MAZ targeting up-regulation of NDUFS3 expression to promote malignant progression in melanoma

Myc-associated Zinc-finger Protein (MAZ) has been implicated in the malignant progression of various tumors. However, its expression and functional relationship of MAZ in melanoma have not been previously investigated. This study confirms elevated expression of MAZ in melanoma, correlating with poor patient prognosis. Furthermore, our findings demonstrate that MAZ enhances melanoma progression by promoting proliferation, migration and invasion. It is worth noting that we found that MAZ can target and regulate the transcription of NADH dehydrogenase [ubiquinone] iron-sulfur protein 3 (NDUFS3), a core subunit of mitochondrial complex I, to enhance mitochondrial metabolism and thus promote malignant progression of melanoma. Predictive modeling indicates that the co-expression of MAZ and NDUFS3 could serve as a potential prognostic marker for melanoma patients.

The MVPs (masterful versatile players): Chromatin factors as pivotal mediators between 3D genome organization and the response to environment

In recent years, the study of genome dynamics has become a prominent research field due to its influence on understanding the control of gene expression. The study of 3D genome organization has unveiled multiple mechanisms in orchestrating chromosome folding. Growing evidence reveals that these mechanisms are not only important for genome organization, but play a pivotal role in enabling plants to adapt to environmental stimuli. In this review, we provide an overview of the current knowledge concerning epigenetic factors and regulatory elements driving 3D genome dynamics and their responses to external stimuli. We discuss the most recent findings, previous evidence, and explore their implications for future research.

ETS1 and RBPJ transcriptionally regulate METTL14 to suppress TGF-β1-induced epithelial-mesenchymal transition in human bronchial epithelial cells

Asthma is a chronic respiratory disease characterized by airway inflammation and remodeling. Epithelial-mesenchymal transition (EMT) of bronchial epithelial cells is considered to be a crucial player in asthma. Methyltransferase-like 14 (METTL14), an RNA methyltransferase, is implicated in multiple pathological processes, including EMT, cell proliferation and migration. However, the role of METTL14 in asthma remains uncertain. This research aimed to explore the biological functions of METTL14 in asthma and its underlying upstream mechanisms. METTL14 expression was down-regulated in asthmatic from three GEO datasets (GSE104468, GSE165934, and GSE74986). Consistent with this trend, METTL14 was decreased in the lung tissues of OVA-induced asthmatic mice and transforming growth factor-β1 (TGF-β1)-stimulated human bronchial epithelial cells (Beas-2B) in this study. Overexpression of METTL14 caused reduction in mesenchymal markers (FN1, N-cad, Col-1 and α-SMA) in TGF-β1-treated cells, but caused increase in epithelial markers (E-cad), thus inhibiting EMT. Also, METTL14 suppressed the proliferation and migration ability of TGF-β1-treated Beas-2B cells. Two transcription factors, ETS1 and RBPJ, could both bind to the promoter region of METTL14 and drive its expression. Elevating METTL14 expression could reversed EMT, cell proliferation and migration promoted by ETS1 or RBPJ deficiency. These results indicate that the ETS1/METTL14 and RBPJ/METTL14 transcription axes exhibit anti-EMT, anti-proliferation and anti-migration functions in TGF-β1-induced bronchial epithelial cells, implying that METTL14 may be considered an alternative candidate target for the treatment of asthma.

YY1-controlled regulatory connectivity and transcription are influenced by the cell cycle

Few transcription factors have been examined for their direct roles in physically connecting enhancers and promoters. Here acute degradation of Yin Yang 1 (YY1) in erythroid cells revealed its requirement for the maintenance of numerous enhancer-promoter loops, but not compartments or domains. Despite its reported ability to interact with cohesin, the formation of YY1-dependent enhancer-promoter loops does not involve stalling of cohesin-mediated loop extrusion. Integrating mitosis-to-G1-phase dynamics, we observed partial retention of YY1 on mitotic chromatin, predominantly at gene promoters, followed by rapid rebinding during mitotic exit, coinciding with enhancer-promoter loop establishment. YY1 degradation during the mitosis-to-G1-phase interval revealed a set of enhancer-promoter loops that require YY1 for establishment during G1-phase entry but not for maintenance in interphase, suggesting that cell cycle stage influences YY1's architectural function. Thus, as revealed here for YY1, chromatin architectural functions of transcription factors can vary in their interplay with CTCF and cohesin as well as by cell cycle stage.

NUP98 oncofusions in myeloid malignancies: An update on molecular mechanisms and therapeutic opportunities

Acute myeloid leukemia (AML) is an aggressive hematological malignancy with a heterogeneous molecular landscape. In the pediatric context, the NUP98 gene is a frequent target of chromosomal rearrangements that are linked to poor prognosis and unfavorable treatment outcomes in different AML subtypes. The translocations fuse NUP98 to a diverse array of partner genes, resulting in fusion proteins with novel functions. NUP98 fusion oncoproteins induce aberrant biomolecular condensation, abnormal gene expression programs, and re-wired protein interactions which ultimately cause alterations in the cell cycle and changes in cellular structures, all of which contribute to leukemia development. The extent of these effects is steered by the functional domains of the fusion partners and the influence of concomitant somatic mutations. In this review, we discuss the complex characteristics of NUP98 fusion proteins and potential novel therapeutic approaches for NUP98 fusion-driven AML.

The relationship between nanoscale genome organization and gene expression in mouse embryonic stem cells during pluripotency transition

During early development, gene expression is tightly regulated. However, how genome organization controls gene expression during the transition from naïve embryonic stem cells to epiblast stem cells is still poorly understood. Using single-molecule microscopy approaches to reach nanoscale resolution, we show that genome remodeling affects gene transcription during pluripotency transition. Specifically, after exit from the naïve pluripotency state, chromatin becomes less compacted, and the OCT4 transcription factor has lower mobility and is more bound to its cognate sites. In epiblast cells, the active transcription hallmark, H3K9ac, decreases within the Oct4 locus, correlating with reduced accessibility of OCT4 and, in turn, with reduced expression of Oct4 nascent RNAs. Despite the high variability in the distances between active pluripotency genes, distances between Nodal and Oct4 decrease during epiblast specification. In particular, highly expressed Oct4 alleles are closer to nuclear speckles during all stages of the pluripotency transition, while only a distinct group of highly expressed Nodal alleles are in close proximity to Oct4 when associated with a nuclear speckle in epiblast cells. Overall, our results provide new insights into the role of the spatiotemporal genome remodeling during mouse pluripotency transition and its correlation with the expression of key pluripotency genes.

Systematic dissection of sequence features affecting binding specificity of a pioneer factor reveals binding synergy between FOXA1 and AP-1

Despite the unique ability of pioneer factors (PFs) to target nucleosomal sites in closed chromatin, they only bind a small fraction of their genomic motifs. The underlying mechanism of this selectivity is not well understood. Here, we design a high-throughput assay called chromatin immunoprecipitation with integrated synthetic oligonucleotides (ChIP-ISO) to systematically dissect sequence features affecting the binding specificity of a classic PF, FOXA1, in human A549 cells. Combining ChIP-ISO with in vitro and neural network analyses, we find that (1) FOXA1 binding is strongly affected by co-binding transcription factors (TFs) AP-1 and CEBPB; (2) FOXA1 and AP-1 show binding cooperativity in vitro; (3) FOXA1's binding is determined more by local sequences than chromatin context, including eu-/heterochromatin; and (4) AP-1 is partially responsible for differential binding of FOXA1 in different cell types. Our study presents a framework for elucidating genetic rules underlying PF binding specificity and reveals a mechanism for context-specific regulation of its binding.

Dynamic Changes in Histone Modifications Are Associated with Differential Chromatin Interactions

Eukaryotic genomes are organized into chromatin domains through long-range chromatin interactions which are mediated by the binding of architectural proteins, such as CTCF and cohesin, and histone modifications. Based on the published Hi-C and ChIP-seq datasets in human monocyte-derived macrophages, we identified 206 and 127 differential chromatin interactions (DCIs) that were not located within transcription readthrough regions in influenza A virus- and interferon β-treated cells, respectively, and found that the binding positions of CTCF and RAD21 within more than half of the DCI sites did not change. However, five histone modifications, H3K4me3, H3K27ac, H3K36me3, H3K9me3, and H3K27me3, showed significantly more dramatic changes than CTCF and RAD21 within the DCI sites. For H3K4me3, H3K27ac, H3K36me3, and H3K27me3, significantly more dramatic changes were observed outside than within the DCI sites. We further applied a motif scanning approach to discover proteins that might correlate with changes in histone modifications and chromatin interactions and found that PRDM9, ZNF384, and STAT2 frequently bound to DNA sequences corresponding to 1 kb genomic intervals with gains or losses of a histone modification within the DCI sites. This study explores the dynamic regulation of chromatin interactions and extends the current knowledge of the relationship between histone modifications and chromatin interactions.

How Transcription Factor Clusters Shape the Transcriptional Landscape

In eukaryotic cells, gene transcription typically occurs in discrete periods of promoter activity, interspersed with intervals of inactivity. This pattern deviates from simple stochastic events and warrants a closer examination of the molecular interactions that activate the promoter. Recent studies have identified transcription factor (TF) clusters as key precursors to transcriptional bursting. Often, these TF clusters form at chromatin segments that are physically distant from the promoter, making changes in chromatin conformation crucial for promoter-TF cluster interactions. In this review, I explore the formation and constituents of TF clusters, examining how the dynamic interplay between chromatin architecture and TF clustering influences transcriptional bursting. Additionally, I discuss techniques for visualizing TF clusters and provide an outlook on understanding the remaining gaps in this field.

A continuum of zinc finger transcription factor retention on native chromatin underlies dynamic genome organization

Transcription factor (TF) residence on chromatin translates into quantitative transcriptional or structural outcomes on genome. Commonly used formaldehyde crosslinking fixes TF-DNA interactions cumulatively and compromises the measured occupancy level. Here we mapped the occupancy level of global or individual zinc finger TFs like CTCF and MAZ, in the form of highly resolved footprints, on native chromatin. By incorporating reinforcing perturbation conditions, we established S-score, a quantitative metric to proxy the continuum of CTCF or MAZ retention across different motifs on native chromatin. The native chromatin-retained CTCF sites harbor sequence features within CTCF motifs better explained by S-score than the metrics obtained from other crosslinking or native assays. CTCF retention on native chromatin correlates with local SUMOylation level, and anti-correlates with transcriptional activity. The S-score successfully delineates the otherwise-masked differential stability of chromatin structures mediated by CTCF, or by MAZ independent of CTCF. Overall, our study established a paradigm continuum of TF retention across binding sites on native chromatin, explaining the dynamic genome organization.

SETDB1 deletion causes DNA demethylation and upregulation of multiple zinc-finger genes

BACKGROUND

SETDB1 (SET domain bifurcated-1) is a histone H3-lysine 9 (H3K9)-specific methyltransferase that mediates heterochromatin formation and repression of target genes. Despite the assumed functional link between DNA methylation and SETDB1-mediated H3K9 trimethylations, several studies have shown that SETDB1 operates autonomously of DNA methylation in a region- and cell-specific manner. This study analyzes SETDB1-null HAP1 cells through a linked methylome and transcriptome analysis, intending to explore genes controlled by SETDB1-involved DNA methylation.

METHODS AND RESULTS

We investigated SETDB1-mediated regulation of DNA methylation and gene transcription in human HAP1 cells using reduced-representation bisulfite sequencing (RRBS) and RNA sequencing. While two-thirds of differentially methylated CpGs (DMCs) in genic regions were hypomethylated in SETDB1-null cells, we detected a plethora of C2H2-type zinc-finger protein genes (C2H2-ZFP, 223 of 749) among the DMC-associated genes. Most C2H2-ZFPs with DMCs in their promoters were found hypomethylated in SETDB1-KO cells, while other non-ZFP genes with promoter DMCs were not. These C2H2-ZFPs with DMCs in their promoters were significantly upregulated in SETDB1-KO cells. Similarly, C2H2-ZFP genes were upregulated in SETDB1-null 293T cells, suggesting that SETDB1's function in ZFP gene repression is widespread. There are several C2H2-ZFP gene clusters on chromosome 19, which were selectively hypomethylated in SETDB1-KO cells.

CONCLUSIONS

SETDB1 collectively and specifically represses a substantial fraction of the C2H2-ZFP gene family. Through the en-bloc silencing of a set of ZFP genes, SETDB1 may help establish a panel of ZFP proteins that are expressed cell-type specifically and thereby can serve as signature proteins for cellular identity.

Transcriptome analysis reveals the effect of cold storage time on the expression of genes related to oxidative metabolism in Chinese black truffle

Chinese black truffle (Tuber indicum) is a hypogenous fungus of great value due to its distinctive aroma. In this study, both transcriptome and physicochemical analyses were performed to investigate the changes of nutrients and gene expression in truffle fruiting bodies during cold storage. The results of physicochemical analysis revealed the active metabolism of fruiting bodies in cold storage, showing the decreased contents of protein and soluble sugar, the variations in both polyphenol oxidase activity and total phenol content, and the detrimental effect of reactive oxygen species production caused by heavy metals (cadmium and lead) in truffles. Transcriptome analysis identified a total of 139,489 unigenes. Down-regulated expression of genes encoding the catalase-like domain-containing protein (katE), glutaredoxin protein (GRX), a copper/zinc superoxide dismutase (Sod_Cu), and aspartate aminotransferase (AAT) affected the degradation metabolism of intracellular oxides. Ribulose-5-phosphate-3-epimerase (RPE) was a key enzyme in response to oxidative stress in truffle cells through the pentose phosphate pathway (PPP). A total of 51,612 simple sequence repeats were identified, providing valuable resources for further genetic diversity analysis, molecular breeding, and genetic map-ping in T. indicum. Transcription factors GAL4 and SUF4-like protein were involved in glucose metabolism and histone methylation processes, respectively. Our study provided a fundamental characterization of the physicochemical and molecular variations in T. indicum during the cold storage at 4°C, providing strong experimental evidence to support the improvement of storage quality of T. indicum.

Integrated enhancer regulatory network by enhancer-promoter looping in gastric cancer

Enhancer cis-regulatory elements play critical roles in gene regulation at many stages of cell growth. Enhancers in cancer cells also regulate the transcription of oncogenes. In this study, we performed a comprehensive analysis of long-range chromatin interactions, histone modifications, chromatin accessibility and expression in two gastric cancer (GC) cell lines compared to normal gastric epithelial cells. We found that GC-specific enhancers marked by histone modifications can activate a population of genes, including some oncogenes, by interacting with their proximal promoters. In addition, motif analysis of enhancer-promoter interacting enhancers showed that GC-specific transcription factors are enriched. Among them, we found that MYB is crucial for GC cell growth and activated by the enhancer with an enhancer-promoter loop and TCF7 upregulation. Clinical GC samples showed epigenetic activation of enhancers at the MYB locus and significant upregulation of TCF7 and MYB, regardless of molecular GC subtype and clinicopathological factors. Single-cell RNA sequencing of gastric mucosa with intestinal metaplasia showed high expression of TCF7 and MYB in intestinal stem cells. When we inactivated the loop-forming enhancer at the MYB locus using CRISPR interference (dCas9-KRAB), GC cell growth was significantly inhibited. In conclusion, we identified MYB as an oncogene activated by a loop-forming enhancer and contributing to GC cell growth.

Deciphering the dynamic code: DNA recognition by transcription factors in the ever-changing genome

Transcription factors (TFs) intricately navigate the vast genomic landscape to locate and bind specific DNA sequences for the regulation of gene expression programs. These interactions occur within a dynamic cellular environment, where both DNA and TF proteins experience continual chemical and structural perturbations, including epigenetic modifications, DNA damage, mechanical stress, and post-translational modifications (PTMs). While many of these factors impact TF-DNA binding interactions, understanding their effects remains challenging and incomplete. This review explores the existing literature on these dynamic changes and their potential impact on TF-DNA interactions.

Predicting RNA polymerase II transcriptional elongation pausing and associated histone code

RNA Polymerase II (Pol II) transcriptional elongation pausing is an integral part of the dynamic regulation of gene transcription in the genome of metazoans. It plays a pivotal role in many vital biological processes and disease progression. However, experimentally measuring genome-wide Pol II pausing is technically challenging and the precise governing mechanism underlying this process is not fully understood. Here, we develop RP3 (RNA Polymerase II Pausing Prediction), a network regularized logistic regression machine learning method, to predict Pol II pausing events by integrating genome sequence, histone modification, gene expression, chromatin accessibility, and protein-protein interaction data. RP3 can accurately predict Pol II pausing in diverse cellular contexts and unveil the transcription factors that are associated with the Pol II pausing machinery. Furthermore, we utilize a forward feature selection framework to systematically identify the combination of histone modification signals associated with Pol II pausing. RP3 is freely available at https://github.com/AMSSwanglab/RP3.

The single-cell chromatin landscape in gonadal cell lineage specification

Gonad development includes sex determination and divergent maturation of the testes and ovaries. Recent advances in measuring gene expression in single cells are providing new insights into this complex process. However, the underlying epigenetic regulatory mechanisms remain unclear. Here, we profiled chromatin accessibility in mouse gonadal cells of both sexes from embryonic day 11.5 to 14.5 using single-cell assay for transposase accessible chromatin by sequencing (scATAC-seq). Our results showed that individual cell types can be inferred by the chromatin landscape, and that cells can be temporally ordered along developmental trajectories. Integrative analysis of transcriptomic and chromatin-accessibility maps identified multiple putative regulatory elements proximal to key gonadal genes Nr5a1, Sox9 and Wt1. We also uncover cell type-specific regulatory factors underlying cell type specification. Overall, our results provide a better understanding of the epigenetic landscape associated with the progressive restriction of cell fates in the gonad.

TEAD2 initiates ground-state pluripotency by mediating chromatin looping

The transition of mouse embryonic stem cells (ESCs) between serum/LIF and 2i(MEK and GSK3 kinase inhibitor)/LIF culture conditions serves as a valuable model for exploring the mechanisms underlying ground and confused pluripotent states. Regulatory networks comprising core and ancillary pluripotency factors drive the gene expression programs defining stable naïve pluripotency. In our study, we systematically screened factors essential for ESC pluripotency, identifying TEAD2 as an ancillary factor maintaining ground-state pluripotency in 2i/LIF ESCs and facilitating the transition from serum/LIF to 2i/LIF ESCs. TEAD2 exhibits increased binding to chromatin in 2i/LIF ESCs, targeting active chromatin regions to regulate the expression of 2i-specific genes. In addition, TEAD2 facilitates the expression of 2i-specific genes by mediating enhancer-promoter interactions during the serum/LIF to 2i/LIF transition. Notably, deletion of Tead2 results in reduction of a specific set of enhancer-promoter interactions without significantly affecting binding of chromatin architecture proteins, CCCTC-binding factor (CTCF), and Yin Yang 1 (YY1). In summary, our findings highlight a novel prominent role of TEAD2 in orchestrating higher-order chromatin structures of 2i-specific genes to sustain ground-state pluripotency.

A fine-scale Arabidopsis chromatin landscape reveals chromatin conformation-associated transcriptional dynamics

Plants, as sessile organisms, deploy transcriptional dynamics for adapting to extreme growth conditions such as cold stress. Emerging evidence suggests that chromatin architecture contributes to transcriptional regulation. However, the relationship between chromatin architectural dynamics and transcriptional reprogramming in response to cold stress remains unclear. Here, we apply a chemical-crosslinking assisted proximity capture (CAP-C) method to elucidate the fine-scale chromatin landscape, revealing chromatin interactions within gene bodies closely associated with RNA polymerase II (Pol II) densities across initiation, pausing, and termination sites. We observe dynamic changes in chromatin interactions alongside Pol II activity alterations during cold stress, suggesting local chromatin dynamics may regulate Pol II activity. Notably, cold stress does not affect large-scale chromatin conformations. We further identify a comprehensive promoter-promoter interaction (PPI) network across the genome, potentially facilitating co-regulation of gene expression in response to cold stress. Our study deepens the understanding of chromatin conformation-associated gene regulation in plant response to cold.

Integrative analysis of transcriptomic and epigenomic data reveals distinct patterns for developmental and housekeeping gene regulation

BACKGROUND

Regulation of transcription is central to the emergence of new cell types during development, and it often involves activation of genes via proximal and distal regulatory regions. The activity of regulatory elements is determined by transcription factors (TFs) and epigenetic marks, but despite extensive mapping of such patterns, the extraction of regulatory principles remains challenging.

RESULTS

Here we study differentially and similarly expressed genes along with their associated epigenomic profiles, chromatin accessibility and DNA methylation, during lineage specification at gastrulation in mice. Comparison of the three lineages allows us to identify genomic and epigenomic features that distinguish the two classes of genes. We show that differentially expressed genes are primarily regulated by distal elements, while similarly expressed genes are controlled by proximal housekeeping regulatory programs. Differentially expressed genes are relatively isolated within topologically associated domains, while similarly expressed genes tend to be located in gene clusters. Transcription of differentially expressed genes is associated with differentially open chromatin at distal elements including enhancers, while that of similarly expressed genes is associated with ubiquitously accessible chromatin at promoters.

CONCLUSION

Based on these associations of (linearly) distal genes' transcription start sites (TSSs) and putative enhancers for developmental genes, our findings allow us to link putative enhancers to their target promoters and to infer lineage-specific repertoires of putative driver transcription factors, within which we define subgroups of pioneers and co-operators.

The influencers' era: how the environment shapes chromatin in 3D

Environment-epigenome interactions are emerging as contributors to disease risk and health outcomes. In fact, organisms outside of the laboratory are constantly exposed to environmental changes that can influence chromatin regulation at multiple levels, potentially impacting on genome function. In this review, we will summarize recent findings on how major external cues impact on 3D chromatin organization in different experimental systems. We will describe environment-induced 3D genome alterations ranging from chromatin accessibility to the spatial distribution of the genome and discuss their role in regulating gene expression.

ZFP541 and KCTD19 regulate chromatin organization and transcription programs for male meiotic progression

The successful progression of meiosis prophase I requires integrating information from the structural and molecular levels. In this study, we show that ZFP541 and KCTD19 work in the same genetic pathway to regulate the progression of male meiosis and thus fertility. The Zfp541 and/or Kctd19 knockout male mice show various structural and recombination defects including detached chromosome ends, aberrant localization of chromosome axis components and recombination proteins, and globally altered histone modifications. Further analyses on RNA-seq, ChIP-seq, and ATAC-seq data provide molecular evidence for the above defects and reveal that ZFP541/KCTD19 activates the expression of many genes by repressing several major transcription repressors. More importantly, we reveal an unexpected role of ZFP541/KCTD19 in directly modulating chromatin organization. These results suggest that ZFP541/KCTD19 simultaneously regulates the transcription cascade and chromatin organization to ensure the coordinated progression of multiple events at chromosome structural and biochemical levels during meiosis prophase I.

Differentiation is accompanied by a progressive loss in transcriptional memory

BACKGROUND

Cell differentiation requires the integration of two opposite processes, a stabilizing cellular memory, especially at the transcriptional scale, and a burst of gene expression variability which follows the differentiation induction. Therefore, the actual capacity of a cell to undergo phenotypic change during a differentiation process relies upon a modification in this balance which favors change-inducing gene expression variability. However, there are no experimental data providing insight on how fast the transcriptomes of identical cells would diverge on the scale of the very first two cell divisions during the differentiation process.

RESULTS

In order to quantitatively address this question, we developed different experimental methods to recover the transcriptomes of related cells, after one and two divisions, while preserving the information about their lineage at the scale of a single cell division. We analyzed the transcriptomes of related cells from two differentiation biological systems (human CD34+ cells and T2EC chicken primary erythrocytic progenitors) using two different single-cell transcriptomics technologies (scRT-qPCR and scRNA-seq).

CONCLUSIONS

We identified that the gene transcription profiles of differentiating sister cells are more similar to each other than to those of non-related cells of the same type, sharing the same environment and undergoing similar biological processes. More importantly, we observed greater discrepancies between differentiating sister cells than between self-renewing sister cells. Furthermore, a progressive increase in this divergence from first generation to second generation was observed when comparing differentiating cousin cells to self renewing cousin cells. Our results are in favor of a gradual erasure of transcriptional memory during the differentiation process.

SpecLoop predicts cell type-specific chromatin loop via transcription factor cooperation

Cell-type-Specific Chromatin Loops (CSCLs) are crucial for gene regulation and cell fate determination. However, the mechanisms governing their establishment remain elusive. Here, we present SpecLoop, a network regularization-based machine learning framework, to investigate the role of transcription factors (TFs) cooperation in CSCL formation. SpecLoop integrates multi-omics data, including gene expression, chromatin accessibility, sequence, protein-protein interaction, and TF binding motif data, to predict CSCLs and identify TF cooperations. Using high resolution Hi-C data as the gold standard, SpecLoop accurately predicts CSCL in GM12878, IMR90, HeLa-S3, K562, HUVEC, HMEC, and NHEK seven cell types, with the AUROC values ranging from 0.8645 to 0.9852 and AUPR values ranging from 0.8654 to 0.9734. Notably SpecLoop demonstrates improved accuracy in predicting long-distance CSCLs and identifies TF complexes with strong predictive ability. Our study systematically explores the TFs and TF pairs associated with CSCL through effective integration of diverse omics data. SpecLoop is freely available at https://github.com/AMSSwanglab/SpecLoop.

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.

Source cell-type epigenetic memory persists in induced pluripotent cells but is lost in subsequently derived germline cells

Introduction: Retention of source cell-type epigenetic memory may mitigate the potential for induced pluripotent stem cells (iPSCs) to fully achieve transitions in cell fate in vitro. While this may not preclude the use of iPSC-derived somatic cell types for therapeutic applications, it becomes a major concern impacting the potential use of iPSC-derived germline cell types for reproductive applications. The transition from a source somatic cell type to iPSCs and then on to germ-cell like cells (GCLCs) recapitulates two major epigenetic reprogramming events that normally occur during development in vivo-embryonic reprogramming in the epiblast and germline reprogramming in primordial germ cells (PGCs). We examined the extent of epigenetic and transcriptomic memory persisting first during the transition from differentiated source cell types to iPSCs, and then during the transition from iPSCs to PGC-like cells (PGCLCs). Methods: We derived iPSCs from four differentiated mouse cell types including two somatic and two germ cell types and tested the extent to which each resulting iPSC line resembled a) a validated ES cell reference line, and b) their respective source cell types, on the basis of genome-wide gene expression and DNA methylation patterns. We then induced each iPSC line to form PGCLCs, and assessed epigenomic and transcriptomic memory in each compared to endogenous PGCs/M-prospermatogonia. Results: In each iPSC line, we found residual gene expression and epigenetic programming patterns characteristic of the corresponding source differentiated cell type from which each was derived. However, upon deriving PGCLCs, we found very little evidence of lingering epigenetic or transcriptomic memory of the original source cell type. Discussion: This result indicates that derivation of iPSCs and then GCLCs from differentiated source cell types in vitro recapitulates the two-phase epigenetic reprogramming that normally occurs in vivo, and that, to a significant extent, germline cell types derived in vitro from pluripotent cells accurately recapitulate epigenetic programming and gene expression patterns corresponding to equivalent endogenous germ cell types, suggesting that they have the potential to form the basis of in vitro gametogenesis as a useful therapeutic strategy for treatment of infertility.

Deficiency of Transcription Factor Sp1 Contributes to Hypertrophic Cardiomyopathy

BACKGROUND

Hypertrophic cardiomyopathy (HCM) is the most prevalent monogenic heart disorder. However, the pathogenesis of HCM, especially its nongenetic mechanisms, remains largely unclear. Transcription factors are known to be involved in various biological processes including cell growth. We hypothesized that SP1 (specificity protein 1), the first purified TF in mammals, plays a role in the cardiomyocyte growth and cardiac hypertrophy of HCM.

METHODS

Cardiac-specific conditional knockout of Sp1 mice were constructed to investigate the role of SP1 in the heart. The echocardiography, histochemical experiment, and transmission electron microscope were performed to analyze the cardiac phenotypes of cardiac-specific conditional knockout of Sp1 mice. RNA sequencing, chromatin immunoprecipitation sequencing, and adeno-associated virus experiments in vivo were performed to explore the downstream molecules of SP1. To examine the therapeutic effect of SP1 on HCM, an SP1 overexpression vector was constructed and injected into the mutant allele of Myh6 R404Q/+ (Myh6 c. 1211C>T) HCM mice. The human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) from a patient with HCM were used to detect the potential therapeutic effects of SP1 in human HCM.

RESULTS

The cardiac-specific conditional knockout of Sp1 mice developed a typical HCM phenotype, displaying overt myocardial hypertrophy, interstitial fibrosis, and disordered myofilament. In addition, Sp1 knockdown dramatically increased the cell area of hiPSC-CMs and caused intracellular myofibrillar disorganization, which was similar to the hypertrophic cardiomyocytes of HCM. Mechanistically, Tuft1 was identified as the key target gene of SP1. The hypertrophic phenotypes induced by Sp1 knockdown in both hiPSC-CMs and mice could be rescued by TUFT1 (tuftelin 1) overexpression. Furthermore, SP1 overexpression suppressed the development of HCM in the mutant allele of Myh6 R404Q/+ mice and also reversed the hypertrophic phenotype of HCM hiPSC-CMs.

CONCLUSIONS

Our study demonstrates that SP1 deficiency leads to HCM. SP1 overexpression exhibits significant therapeutic effects on both HCM mice and HCM hiPSC-CMs, suggesting that SP1 could be a potential intervention target for HCM.

Epigenetic regulatory layers in the 3D nucleus

Nearly 7 decades have elapsed since Francis Crick introduced the central dogma of molecular biology, as part of his ideas on protein synthesis, setting the fundamental rules of sequence information transfer from DNA to RNAs and proteins. We have since learned that gene expression is finely tuned in time and space, due to the activities of RNAs and proteins on regulatory DNA elements, and through cell-type-specific three-dimensional conformations of the genome. Here, we review major advances in genome biology and discuss a set of ideas on gene regulation and highlight how various biomolecular assemblies lead to the formation of structural and regulatory features within the nucleus, with roles in transcriptional control. We conclude by suggesting further developments that will help capture the complex, dynamic, and often spatially restricted events that govern gene expression in mammalian cells.

Enhancer-driven 3D chromatin domain folding modulates transcription in human mammary tumor cells

The genome is organized in functional compartments and structural domains at the sub-megabase scale. How within these domains interactions between numerous cis-acting enhancers and promoters regulate transcription remains an open question. Here, we determined chromatin folding and composition over several hundred kb around estrogen-responsive genes in human breast cancer cell lines after hormone stimulation. Modeling of 5C data at 1.8 kb resolution was combined with quantitative 3D analysis of multicolor FISH measurements at 100 nm resolution and integrated with ChIP-seq data on transcription factor binding and histone modifications. We found that rapid estradiol induction of the progesterone gene expression occurs in the context of preexisting, cell type-specific chromosomal architectures encompassing the 90 kb progesterone gene coding region and an enhancer-spiked 5' 300 kb upstream genomic region. In response to estradiol, interactions between estrogen receptor α (ERα) bound regulatory elements are reinforced. Whereas initial enhancer-gene contacts coincide with RNA Pol 2 binding and transcription initiation, sustained hormone stimulation promotes ERα accumulation creating a regulatory hub stimulating transcript synthesis. In addition to implications for estrogen receptor signaling, we uncover that preestablished chromatin architectures efficiently regulate gene expression upon stimulation without the need for de novo extensive rewiring of long-range chromatin interactions.

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.

Polymer Modeling Reveals Interplay between Physical Properties of Chromosomal DNA and the Size and Distribution of Condensin-Based Chromatin Loops

Transient DNA loops occur throughout the genome due to thermal fluctuations of DNA and the function of SMC complex proteins such as condensin and cohesin. Transient crosslinking within and between chromosomes and loop extrusion by SMCs have profound effects on high-order chromatin organization and exhibit specificity in cell type, cell cycle stage, and cellular environment. SMC complexes anchor one end to DNA with the other extending some distance and retracting to form a loop. How cells regulate loop sizes and how loops distribute along chromatin are emerging questions. To understand loop size regulation, we employed bead-spring polymer chain models of chromatin and the activity of an SMC complex on chromatin. Our study shows that (1) the stiffness of the chromatin polymer chain, (2) the tensile stiffness of chromatin crosslinking complexes such as condensin, and (3) the strength of the internal or external tethering of chromatin chains cooperatively dictate the loop size distribution and compaction volume of induced chromatin domains. When strong DNA tethers are invoked, loop size distributions are tuned by condensin stiffness. When DNA tethers are released, loop size distributions are tuned by chromatin stiffness. In this three-way interaction, the presence and strength of tethering unexpectedly dictates chromatin conformation within a topological domain.

Interchromosomal Colocalization with Parental Genes Is Linked to the Function and Evolution of Mammalian Retrocopies

Retrocopies are gene duplicates arising from reverse transcription of mature mRNA transcripts and their insertion back into the genome. While long being regarded as processed pseudogenes, more and more functional retrocopies have been discovered. How the stripped-down retrocopies recover expression capability and become functional paralogs continually intrigues evolutionary biologists. Here, we investigated the function and evolution of retrocopies in the context of 3D genome organization. By mapping retrocopy-parent pairs onto sequencing-based and imaging-based chromatin contact maps in human and mouse cell lines and onto Hi-C interaction maps in 5 other mammals, we found that retrocopies and their parental genes show a higher-than-expected interchromosomal colocalization frequency. The spatial interactions between retrocopies and parental genes occur frequently at loci in active subcompartments and near nuclear speckles. Accordingly, colocalized retrocopies are more actively transcribed and translated and are more evolutionarily conserved than noncolocalized ones. The active transcription of colocalized retrocopies may result from their permissive epigenetic environment and shared regulatory elements with parental genes. Population genetic analysis of retroposed gene copy number variants in human populations revealed that retrocopy insertions are not entirely random in regard to interchromosomal interactions and that colocalized retroposed gene copy number variants are more likely to reach high frequencies, suggesting that both insertion bias and natural selection contribute to the colocalization of retrocopy-parent pairs. Further dissection implies that reduced selection efficacy, rather than positive selection, contributes to the elevated allele frequency of colocalized retroposed gene copy number variants. Overall, our results hint a role of interchromosomal colocalization in the "resurrection" of initially neutral retrocopies.