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

High mobility group box-1 protein-mediated class II major histocompatibility complex transactivator superenhancers are critical for dendritic cell-trained immunity in acute-to-chronic progression of allograft rejection

Chronic allograft rejection is mainly mediated by indirect recognition. Dendritic cells (DCs), as the major antigen-presenting cells in indirect recognition, exhibit an enhanced antigen-presenting ability in chronic rejection, but the specific mechanism is still unclear. Here, we found that pretreatment with high mobility group box-1 protein (HMGB1) in vivo can induce trained immunity in DCs. These trained DCs demonstrated an enhanced ability to present alloantigen, accelerating allograft rejection in a CTLA4-Ig-induced chronic rejection model by upregulating the expression of major histocompatibility complex (MHC)-II and class II major histocompatibility complex transactivator (CIITA) molecules. Mechanistically, we found that HMGB1 promoted the formation of superenhancers (SEs) of CIITA, epigenetically reprogramming DCs and promoting trained immunity. The SE inhibitor JQ1 reduced the expression of CIITA and MHC-II in DCs, thereby delaying the occurrence of chronic rejection. Interestingly, we identified HMGB1 as a specific inducer of SE formation in a newly named SEa region of CIITA. Targeted knockout of the CIITA's SEa region inhibited HMGB1-induced trained immunity in DCs. Taken together, our data confirm that HMGB1 can induce the formation of the SEs of CIITA, promote trained immunity in DCs, and accelerate allograft rejection, thus offering a new potential target for the treatment of chronic rejection.

Cis-Regulation of the CFTR Gene in Pancreatic Cells

Genome organization is essential for precise spatial and temporal gene expression and relies on interactions between promoters and distal cis-regulatory elements (CREs), which constitute ~8% of the human genome. For the cystic fibrosis transmembrane conductance regulator (CFTR) gene, tissue-specific expression, especially in the pancreas, remains poorly understood. Unraveling its regulation could clarify the clinical heterogeneity observed in cystic fibrosis and CFTR-related disorders. To understand the role of 3D chromatin architecture in establishing tissue-specific expression of the CFTR gene, we mapped chromatin interactions and epigenomic regulation in Capan-1 pancreatic cells. Candidate CREs are validated by luciferase reporter assay and CRISPR knock-out. We identified active CREs not only around the CFTR gene but also outside the topologically associating domain (TAD). We demonstrate the involvement of multiple CREs upstream and downstream of the CFTR gene and reveal a cooperative effect of the -44 kb, -35 kb, +15.6 kb, and +37.7 kb regions, which share common predicted transcription factor (TF) motifs. We also extend our analysis to compare 3D chromatin conformation in intestinal and pancreatic cells, providing valuable insights into the tissue specificity of CREs in regulating CFTR gene expression.

The transcription factor PITX1 cooperates with super-enhancers to regulate the expression of DUSP4 and inhibit pyroptosis in pulmonary artery smooth muscle cells

BACKGROUND

Pulmonary hypertension (PH) is a highly fatal pathophysiological syndrome. The group 1 pulmonary arterial hypertension (PAH) is characterized by acute pulmonary vasoconstriction and chronic vascular remodeling caused by hyperplasia and hypertrophy of pulmonary artery smooth muscle cells (PASMCs) and chronic inflammation. Pyroptosis is an inflammatory mode of cell death that is regulated by super-enhancers (SEs) and occurs in the setting of tumors and cardiovascular diseases. However, whether SEs are involved in the pathological process of pyroptosis in PAH and the specific mechanism involved remain unclear.

METHODS

Here, we identified the SE target gene DUSP4 via ChIP-seq with an anti-H3K27ac antibody, and bioinformatics predictions revealed that the transcription factor PITX1 can bind to the promoter and SE sequences of DUSP4. The AAV5 vector was used to deliver shRNAs targeting PITX1 and DUSP4 to PASMCs.

RESULTS

PITX1 overexpression reversed the increase in right ventricular systolic pressure and pulmonary vascular remodeling, restored the PAAT/PAVTI ratio in hypoxic pulmonary hypertension (HPH, Group 3 PH) and SuHx PAH (Group 1 PAH) mice, and suppressed pyroptosis in pulmonary vascular cells. However, knockdown of DUSP4 counteracted the effects of PITX1 overexpression. Similar results were obtained in cultured PASMCs. In addition, treatment with the SE inhibitors JQ1 and iBET decreased the transcription of DUSP4 and increased the expression of hypoxia-induced pyroptosis proteins in PASMCs.

CONCLUSION

We confirmed that PITX1 can promote DUSP4 expression by binding to the DUSP4 promoter and SE to reduce pyroptosis in hypoxic PASMCs, providing new insights into the role of SEs and pyroptosis in pulmonary vascular remodeling and a theoretical basis for the treatment of PAH and related diseases.

Deciphering enhancers of hearing loss genes for efficient and targeted gene therapy of hereditary deafness

Hereditary hearing loss accounts for about 60% of congenital deafness. Although adeno-associated virus (AAV)-mediated gene therapy shows substantial potential for treating genetic hearing impairments, there remain significant concerns regarding the specificity and safety of AAV vectors. The sophisticated nature of the cochlea further complicates the challenge of precisely targeting gene delivery. Here, we introduced an AAV-reporter-based in vivo transcriptional enhancer reconstruction (ARBITER) workflow, enabling efficient and reliable dissection of enhancers. With ARBITER, we successfully demonstrated that the conserved non-coding elements (CNEs) within the gene locus collaboratively regulate the expression of Slc26a5, which was further validated using knockout mouse models. We also assessed the potential of identified enhancers to treat hereditary hearing loss by conducting gene therapy in Slc26a5 mutant mice. Based on the original Slc26a5 enhancer with limited efficiency, we engineered a highly efficient and outer hair cell (OHC)-specific enhancer, B8, which successfully restored hearing of Slc26a5 knockout mice.

Systemic genome-epigenome analysis captures a lineage-specific super-enhancer for MYB in gastrointestinal adenocarcinoma

Gastrointestinal adenocarcinoma is a major cancer type for the digestive system, ranking as the top cause of cancer-related deaths worldwide. While there has been extensive research on mutations in protein-coding regions, the knowledge of the landscape of its non-coding regulatory elements is still insufficient. Combining the analysis of active enhancer profiles and genomic structural variation, we discovered and validated a lineage-specific super-enhancer for MYB in gastrointestinal adenocarcinoma. This super-enhancer is composed of a predominant enhancer e4 and several additional enhancers, whose transcriptional activity is regulated by the direct binding of HNF4A and MYB itself. Suppression of the super-enhancer downregulated the expression of MYB, inhibited downstream Notch signaling and prevented the development of gastrointestinal adenocarcinoma both in vitro and in vivo. Our study uncovers a mechanism driven by non-coding variations that regulate MYB expression in a lineage-specific manner, offering new insights into the carcinogenic mechanism and potential therapeutic strategies for gastrointestinal adenocarcinoma.

Active regulatory elements recruit cohesin to establish cell specific chromatin domains

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

Principles of long-range gene regulation

Transcription from gene promoters occurs in specific spatiotemporal patterns in multicellular organisms, controlled by genomic regulatory elements. The communication between a regulatory element and a promoter requires a certain degree of physical proximity between them; hence, most gene regulation occurs locally in the genome. However, recent discoveries have revealed long-range gene regulation strategies that enhance interactions between regulatory elements and promoters by overcoming the distances between them in the linear genome. These new findings challenge the traditional view of how gene expression patterns are controlled. This review examines long-range gene regulation strategies recently reported in Drosophila and mammals, offering insights into their mechanisms and evolution.

An unusual transfusion-dependent hemoglobin H disease caused by a novel complex inverted duplication involving the α-globin regulatory elements and α-thalassemia--SEA deletion

Four multi-species conserved sequences (MCSs) are important enhancers which affect α-globin expression. Deletions of MCS can cause α-thalassemia. So far, duplication of MCS has never been reported to account for thalassemia. In this study, an unusual transfusion-dependent case of hemoglobin H disease was identified by whole-genome sequencing, optical genome mapping and longer PCR with special primers, which was caused by a familial 96,620-bp inverted duplication (from MCS-R1 to MCS-R4), inserted between chr16:199348 and 199349 (GRCh37/hg19) within MCSs. The duplication segment included an inverted repeat sequence from chr16:102712 to176193 and one direct repeat sequence from chr16:176208 to 199348. The associated α-thalassemia trait was confirmed to result from disrupted topological chromatin domains using ATAC-seq and the dual‑luciferase reporter assay system. This case presents a new mechanism of α-thalassemia, and may aid our understanding of the effects of enhancers on gene expression and the differential contribution of the four enhancer elements in the human a-globin locus.

Synergistic and antagonistic activities of IRF8 and FOS enhancer pairs during an immune-cell fate switch

Cell fate instructive genes tend to be regulated by large clusters of enhancers. Whether and how individual enhancers within such clusters cooperate in regulating gene expression is poorly understood. We have previously developed a computational method, SEGCOND, which identifies hubs that we termed Putative Transcriptional Condensates (PTCs), consisting of enhancer clusters and associated target genes. Here, we use SEGCOND to identify PTCs in a CEBPA-induced B-cell-to-macrophage transdifferentiation system. We find that PTCs are enriched for highly expressed, lineage-restricted genes and associate with BRD4, a component of transcriptional condensates. Further, we performed single and combinatorial deletions of enhancers within two PTCs active during induced transdifferentiation, harboring IRF8 and FOS. Two enhancers within the IRF8 PTC were found to provide a backup mechanism when combined, safeguarding IRF8 expression and efficient transdifferentiation. Unexpectedly, two individual enhancers within the FOS PTC antagonize each other on day 1 of transdifferentiation, delaying the conversion of B-cells into macrophages and reducing FOS expression, while on day 7, they cooperate to increase FOS levels induced cells. Our results reveal complex, differentiation-stage-specific interactions between individual enhancers within enhancer clusters.

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

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

Potential target within the tumor microenvironment - MT1-MMP

Matrix metalloproteinases are integral to the modification of the tumor microenvironment and facilitate tumor progression by degrading the extracellular matrix, releasing cytokines, and influencing the recruitment of immune cells. Among the matrix metalloproteinases, membrane-type matrix metalloproteinase 1 (MT1-MMP/MMP14) is the first identified membrane-type MMP and acts as an essential proteolytic enzyme that enables tumor infiltration and metastatic progression. Given the pivotal role of MT1-MMP in tumor progression and the correlation between its overexpression in tumors and unfavorable prognoses across multiple cancer types, a comprehensive understanding of the potential functional mechanisms of MT1-MMP is essential. This knowledge will aid in the advancement of diverse anti-tumor therapies aimed at targeting MT1-MMP. Although contemporary research has highlighted the considerable potential of MT1-MMP in targeted cancer therapy, studies pertaining to its application in cell therapy remain relatively limited. In this review, we delineate the structural characteristics and regulatory mechanisms of MT1-MMP expression, as well as its biological significance in tumorigenesis. Finally, we discussed the current status and prospects of anti-tumor therapies targeting MT1-MMP.

In vivo deletion of a GWAS-identified Myb distal enhancer acts on Myb expression, globin switching, and clinical erythroid parameters in β-thalassemia

Genome-wide association studies (GWAS) have identified numerous genetic variants linked to human diseases, mostly located in non-coding regions of the genome, particularly in putative enhancers. However, functional assessment of the non-coding GWAS variants has progressed at slow pace, since the functions of the vast majority of genomic enhancers have not been defined, impeding interpretation of disease-susceptibility variants. The HBS1L-MYB intergenic region harbors multiple SNPs associated with clinical erythroid parameters, including fetal hemoglobin levels, a feature impacting disease severity of beta-hemoglobinopathies such as sickle cell anemia and beta-thalassemia. HBS1L-MYB variants cluster in the vicinity of several MYB enhancers, altering MYB expression and globin switching. We and others have highlighted the conserved human MYB - 84kb enhancer, known as the - 81kb enhancer in the mouse, as likely candidate linked to these traits. We report here the generation of a Myb - 81kb enhancer knock-out mouse model, and shed light for the first time on its impact on steady state erythropoiesis and in beta-thalassemia in vivo.

Long-range enhancer-controlled genes are hypersensitive to regulatory factor perturbations

Cell-type-specific gene activation is regulated by enhancers, sometimes located at large genomic distances from target gene promoters. Whether distal enhancers require specific factors to orchestrate gene regulation remains unclear. Here, we used enhancer distance-controlled reporter screens to find candidate factors. We depleted them and employed activity-by-contact predictions to genome-wide classify genes based on enhancer distance. Predicted distal enhancers typically control tissue-restricted genes and often are strong enhancers. We find cohesin, but also mediator, most specifically required for long-range activation, with cohesin repressing short-range gene activation and prioritizing distal over proximal HBB genes competing for shared enhancers. Long-range controlled genes are also most sensitive to perturbations of other regulatory proteins and to BET inhibitor JQ1, this being more a consequence of their distinct enhancer features than distance. Our work predicts that lengthening of intervening sequences can help limit the expression of target genes to specialized cells with optimal trans-factor environments.

Bioinformatics analysis of electroacupuncture treatment for ischemic stroke: exploring transcriptional regulatory mechanisms mediated by super-enhancers

Background

Ischemic stroke is a leading cause of disability and mortality, imposing substantial physical, emotional, and economic burdens on patients and society. This study aimed to explore the regulatory effects of super-enhancers (SEs) on gene expression in the context of ischemic stroke and their potential transcriptional regulatory mechanisms.

Methods

Super-enhancers were identified via H3K27ac chromatin immunoprecipitation sequencing (ChIP-seq) and ROSE software. RNA-sequencing (RNA-seq) was employed to screen for differentially expressed genes. A comparative analysis of ChIP-seq and RNA-seq data initially identified SE target genes, followed by further screening of key core differentially expressed SE target genes via the random forest method. The identified core SE target genes were initially validated through immunofluorescence and immunoblotting techniques. Additionally, potential core transcriptional regulatory circuits were preliminarily screened via the Coltron algorithm.

Results

We identified SE-associated genes in the ischemic stroke model and electroacupuncture-treated groups, revealing 41 genes uniquely regulated by SEs in the electroacupuncture group compared with 367 in the model group. Enrichment analyses revealed that pathways involved in axon guidance, regulation of lipolysis in adipocytes and sphingolipid signaling pathway were significantly enriched in the SE target genes, suggesting that these pathways may be involved in the therapeutic effects of electroacupuncture. Notably, HDAC7 emerged as a key SE-driven gene; its expression was significantly reduced following electroacupuncture treatment, indicating its potential as a therapeutic target. Protein expression analyses confirmed elevated levels of HDAC7 in the model group, which were reduced by electroacupuncture intervention (p < 0.05). Furthermore, core transcriptional regulatory circuitries involving SOX8, FOXK1, and KLF13 were identified, highlighting their roles in the modulation of SE-mediated gene regulation by acupuncture in the ischemic stroke context.

Conclusion

Overall, our findings provide novel insights into the molecular mechanisms by which acupuncture may treat ischemic stroke, identifying key SE target genes and transcriptional circuits as promising targets for future therapeutic strategies. Further research is warranted to validate these findings in clinical settings and explore the translational potential of acupuncture in ischemic stroke treatment.

Two unrelated distal genes activated by a shared enhancer benefit from localizing inside the same small topological domain

Enhancers are tissue-specific regulatory DNA elements that can activate transcription of genes over distance. Their target genes most often are located in the same contact domain-chromosomal entities formed by cohesin DNA loop extrusion and typically flanked by CTCF-bound boundaries. Enhancers shared by multiple unrelated genes are underexplored but may be more common than anticipated. Here, we analyzed the interplay between an enhancer and two distal functionally unrelated genes residing at opposite domain boundaries. The enhancer strongly activated their expression and supported their frequent interactions. Cohesin structured the domain and supported their transcription, but the genes did not rely on each other's transcription or show gene competition. Deleting either domain boundary not only extended the contact domain but led to reduced contacts within the original domain and reduction in the expression of both genes. Conversely, by isolating either gene with the enhancer in shorter domains, through insertion of new CTCF boundaries, intradomain contact frequencies increased, and the gene isolated with the enhancer was upregulated. Collectively, this shows that an enhancer can independently activate unrelated distal genes and that long-range gene regulation benefits from operating in small contact domains.

Engineered CRISPR-Cas12a for higher-order combinatorial chromatin perturbations

Multiplexed genetic perturbations are critical for testing functional interactions among coding or non-coding genetic elements. Compared to double-stranded DNA cutting, repressive chromatin formation using CRISPR interference (CRISPRi) avoids genotoxicity and is more effective for perturbing non-coding regulatory elements in pooled assays. However, current CRISPRi pooled screening approaches are limited to targeting one to three genomic sites per cell. We engineer an Acidaminococcus Cas12a (AsCas12a) variant, multiplexed transcriptional interference AsCas12a (multiAsCas12a), that incorporates R1226A, a mutation that stabilizes the ribonucleoprotein-DNA complex via DNA nicking. The multiAsCas12a-KRAB fusion improves CRISPRi activity over DNase-dead AsCas12a-KRAB fusions, often rescuing the activities of lentivirally delivered CRISPR RNAs (crRNA) that are inactive when used with the latter. multiAsCas12a-KRAB supports CRISPRi using 6-plex crRNA arrays in high-throughput pooled screens. Using multiAsCas12a-KRAB, we discover enhancer elements and dissect the combinatorial function of cis-regulatory elements in human cells. These results instantiate a group testing framework for efficiently surveying numerous combinations of chromatin perturbations for biological discovery and engineering.

Iterative improvement of deep learning models using synthetic regulatory genomics

Generative deep learning models can accurately reconstruct genome-wide epigenetic tracks from the reference genome sequence alone. But it is unclear what predictive power they have on sequence diverging from the reference, such as disease- and trait-associated variants or engineered sequences. Recent work has applied synthetic regulatory genomics to characterized dozens of deletions, inversions, and rearrangements of DNase I hypersensitive sites (DHSs). Here, we use the state-of-the-art model Enformer to predict DNA accessibility across these engineered sequences when delivered at their endogenous loci. At high level, we observe a good correlation between accessibility predicted by Enformer and experimentally measured values. But model performance was best for sequences that more resembled the reference, such as single deletions or combinations of multiple DHSs. Predictive power was poorer for rearrangements affecting DHS order or orientation. We use these data to fine-tune Enformer, yielding significant reduction in prediction error. We show that this fine-tuning retains strong predictive performance for other tracks. Our results show that current deep learning models perform poorly when presented with novel sequence diverging in certain critical features from their training set. Thus an iterative approach incorporating profiling of synthetic constructs can improve model generalizability, and ultimately enable functional classification of regulatory variants identified by population studies.

Master transcription factor binding sites constitute the core of early replication control elements

Eukaryotic genomes replicate in a defined temporal order called the replication timing (RT) program. RT is developmentally regulated with potential to drive cell fate transitions, but mechanisms controlling RT remain elusive. We previously identified "Early Replication Control Elements" (ERCEs) necessary for early RT, domain-wide transcription, 3D chromatin architecture and compartmentalization in mouse embryonic stem cells (mESCs) but, deletions identifying ERCEs were large and encompassed many putative regulatory elements. Here, we show that ERCEs are compound elements whose RT activity can largely be accounted for by multiple sites of diverse master transcription factor binding (subERCEs), distinguished from other such sites by their long-range interactions. While deletion of subERCEs had large effects on both transcription and RT, deleting transcription start sites eliminated nearly all transcription with moderate effects on RT. Our results suggest a model in which subERCEs respond to diverse master transcription factors by functioning both as transcription enhancers and as elements that organize chromatin domains structurally and support early RT, potentially providing a feed-forward loop to drive robust epigenomic change during cell fate transitions.

Gene editing without ex vivo culture evades genotoxicity in human hematopoietic stem cells

Gene editing the BCL11A erythroid enhancer is a validated approach to fetal hemoglobin (HbF) induction for β-hemoglobinopathy therapy, though heterogeneity in edit allele distribution and HbF response may impact its safety and efficacy. Here, we compare combined CRISPR-Cas9 editing of the BCL11A +58 and +55 enhancers with leading gene modification approaches under clinical investigation. Dual targeting of the BCL11A +58 and +55 enhancers with 3xNLS-SpCas9 and two single guide RNAs (sgRNAs) resulted in superior HbF induction, including in sickle cell disease (SCD) patient xenografts, attributable to simultaneous disruption of core half E-box/GATA motifs at both enhancers. Unintended on-target outcomes of double-strand break (DSB) repair in hematopoietic stem and progenitor cells (HSPCs), such as long deletions and centromere-distal chromosome fragment loss, are a byproduct of cellular proliferation stimulated by ex vivo culture. Editing quiescent HSPCs bypasses long deletion and micronuclei formation and preserves efficient on-target editing and engraftment function.

Signaling Pathways Driving MSC Osteogenesis: Mechanisms, Regulation, and Translational Applications

Mesenchymal stem cells (MSCs) are crucial for skeletal development, homeostasis, and repair, primarily through their differentiation into osteoblasts and other skeletal lineage cells. Key signaling pathways, including Wnt, TGF-β/BMP, PTH, Hedgehog, and IGF, act as critical regulators of MSC osteogenesis, playing pivotal roles in maintaining bone homeostasis and facilitating regeneration. These pathways interact in distinct ways at various stages of bone development, mineralization, and remodeling. This review provides an overview of the molecular mechanisms by which these pathways regulate MSC osteogenesis, their influence on bone tissue formation, and their implications in bone diseases and therapeutic strategies. Additionally, we explore the potential applications of these pathways in bone tissue engineering, with a particular focus on promoting the use of MSCs as seed cells for bone defect repair. Ultimately, this review aims to highlight potential avenues for advancing bone biology research, treating bone disorders, and enhancing regenerative medicine.

Genetic factors mediating long-range enhancer-promoter communication in mammalian development

Enhancers are remotely located noncoding DNA sequences that regulate gene expression in response to developmental, homeostatic, and environmental cues. Canonical short-range enhancers located <50 kb from their cognate promoters function by binding transcription factors, coactivators, and chromatin modifiers. In this review, we discuss recent evidence that medium-range (50-400 kb) and long-range (>400 kb) enhancers rely on additional mechanisms, including cohesin, CCCTC-binding factor, and high-affinity protein-protein interactions. These mechanisms are crucial for establishing the physical proximity and interaction between enhancers and their target promoters over extended genomic distances and ensuring robust gene activation during mammalian development. Future studies will be critical to unravel their prevalence and evolutionary significance across various genomic loci, cell types, and species.

Super-enhancer Activates Master Transcription Factor NR3C1 Expression and Promotes 5-FU Resistance in Gastric Cancer

Poor response to 5-fluorouracil (5-FU) remains an obstacle in the treatment of gastric cancer (GC). Super enhancers (SEs) are crucial for determining tumor cell survival under drug pressure. SE landscapes related to 5-FU-resistance are mapped to GC using chromatin immunoprecipitation-sequencing (ChIP-Seq). SiRNA transcription factors (TFs) screen determines master TF Nuclear Receptor Subfamily 3 Group C Member 1 (NR3C1) activated by SE. High NR3C1 expression driven by SE correlated with 5-FU resistance in patient-derived organoids (PDOs). Phase separation formed by NR3C1 is observed using fluorescence recovery after photobleaching (FRAP). NR3C1 protein and Mediator promoted SE-related gene transcription via phase separation. SEs and NR3C1 co-binding patterns are explored using Cleavage Under Targets and Tagmentation (CUT&Tag) sequencing. 5-FU-related genes driven by NR3C1 are identified using epigenetic reader inhibitor JQ1 and NR3C1 specific inhibitor Cort108297. NR3C1 knockdown increases 5-FU sensitivity and alters the SE landscape through enhancer reprogramming, reducing downstream 5-FU-related target genes. JQ1 and Cort108297 both improve 5-FU efficacy in PDOs and patient-derived xenografts (PDXs) by destroying SEs or inhibiting NR3C1. In conclusion, SE-driven NR3C1 promotes 5-FU resistance in GC. SE destruction and NR3C1 inhibition lead to enhancer reconstruction and reduce 5-FU-related gene transcription, providing alternative therapeutic strategies for improving 5-FU sensitivity.

Recent advances in genome mining and synthetic biology for discovery and biosynthesis of natural products

Natural products have long served as critical raw materials in chemical and pharmaceutical manufacturing, primarily which can provide superior scaffolds or intermediates for drug discovery and development. Over the last century, natural products have contributed to more than a third of therapeutic drug production. However, traditional methods of producing drugs from natural products have become less efficient and more expensive over the past few decades. The combined utilization of genome mining and synthetic biology based on genome sequencing, bioinformatics tools, big data analytics, genetic engineering, metabolic engineering, and systems biology promises to counter this trend. Here, we reviewed recent (2020-2023) examples of genome mining and synthetic biology used to resolve challenges in the production of natural products, such as less variety, poor efficiency, and low yield. Additionally, the emerging efficient tools, design principles, and building strategies of synthetic biology and its application prospects in NPs synthesis have also been discussed.

Exploring the interplay between enhancer-promoter interactions and transcription

Enhancers in metazoan genomes are known to activate their target genes across both short and long genomic distances. Recent advances in chromosome conformation capture assays and single-cell imaging have shed light on the underlying chromatin contacts and dynamics. Yet the relationship between 3D physical enhancer-promoter (E-P) interactions and transcriptional activation remains unresolved. In this brief review, we discuss recent studies exploring this relationship across scales: from developmental stages to the minutes surrounding transcriptional activation and from the tissue level to single-allele subcellular dynamics. We discuss how seemingly contradictory observations might be reconciled and contribute to a refined causal relationship between E-P interactions and transcription, with mutual influences.

Biomolecular condensates in immune cell fate

Fate decisions during immune cell development require temporally precise changes in gene expression. Evidence suggests that the dynamic modulation of these changes is associated with the formation of diverse, membrane-less nucleoprotein assemblies that are termed biomolecular condensates. These condensates are thought to orchestrate fate-determining transcriptional and post-transcriptional processes by locally and transiently concentrating DNA or RNA molecules alongside their regulatory proteins. Findings have established a link between condensate formation and the gene regulatory networks that ensure the proper development of immune cells. Conversely, condensate dysregulation has been linked to impaired immune cell fates, including ageing and malignant transformation. This Review explores the putative mechanistic links between condensate assembly and the gene regulatory frameworks that govern normal and pathological development in the immune system.

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

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

Improved ChIP Sequencing for H3K27ac Profiling and Super-Enhancer Analysis Assisted by Fluorescence-Activated Sorting of Formalin-Fixed Paraffin-Embedded Tissues

Archived clinical formalin-fixed paraffin-embedded tissue (FFPE) is valuable for the study of tumor epigenetics. Although protocol of chromatin immunoprecipitation coupled with next generation sequencing (NGS) (ChIP-seq) using FFPE samples has been established, removal of interference signals from non-target cell components in the samples is still needed. In this study, the protocol of ChIP-seq with purified cells from FFPE lymphoid tissue of nodal T follicular helper cell lymphoma, angioimmunoblastic type (nTFHL-AI) after fluorescence-activated cell sorting (FACS) was established and optimized. Essential steps included single cell preparation, heat treatment enhancing antigen retrieval and labeling, cell sorting, chromatin shearing, ChIP and NGS. Through assistance of FACS, we successfully isolated tumor cells from FFPE lymph node samples of nTFHL-AI and profiled super-enhancers (SEs) mapping by enrichment of H3K27ac signals. The data indicated that the SEs mapping of the sorted cells was different from that of the entire unsorted tissue sample. The H3K27ac signals with cell lineage specificity from background cell components were successfully removed, and the remaining SEs mapping was more similar to T follicular helper cell in an unsupervised clustering analysis, rather than the primary tissue. In addition, we also evaluated the protocol using cultured pure cell lines, and the results indicated that the sequencing data obtained through this protocol had high fidelity and reproducibility. These results show that ChIP-seq for H3K27ac profiling and SEs mapping assisted by FACS with pathological FFPE tissue is available for research of histone modification. Precise epigenetic characteristics of the tumor cell can be described with this protocol.

Enhancer cooperativity can compensate for loss of activity over large genomic distances

Enhancers are short DNA sequences that activate their target promoter from a distance; however, increasing the genomic distance between the enhancer and the promoter decreases expression levels. Many genes are controlled by combinations of multiple enhancers, yet the interaction and cooperation of individual enhancer elements are not well understood. Here, we developed a synthetic platform in mouse embryonic stem cells that allows building complex regulatory landscapes from the bottom up. We tested the system by integrating individual enhancers at different distances and confirmed that the strength of an enhancer contributes to how strongly it is affected by increased genomic distance. Furthermore, synergy between two enhancer elements depends on the distance at which the two elements are integrated: introducing a weak enhancer between a strong enhancer and the promoter strongly increases reporter gene expression, allowing enhancers to activate from increased genomic distances.

Bridging the gap: How enhancers cooperate to regulate gene expression over large genomic distances

By building synthetic regulatory landscapes, Jensen et al.1 and Thomas et al.2 demonstrate in this issue of Molecular Cell that gene expression levels strongly depend on the genomic distance between enhancers and promoters and that enhancer cooperation can compensate for reduced enhancer activity over large genomic distances.

In vivo mechanism of the interaction between trimethylamine lyase expression and glycolytic pathways

Recent studies confirmed that host-gut microbiota interactions modulate disease-linked metabolite TMA production via TMA lyase. However, microbial enzyme production mechanisms remain unclear. In the present study, we investigated the impact of dietary and intervention factors on gut microbiota, microbial gene expression, and the interplay between TMA lyase and glycolytic pathways in mice. Using 16S rRNA gene sequencing, metagenomics, and metabolomics, the gut microbiota composition and microbial functional gene expression profiles related to TMA lyase and glycolytic enzymes were determined. The results revealed that distinct diets and intervention factors altered gut microbiota, gene expression, and metabolites linked to glycine metabolism and glycolysis. Notably, an arabinoxylan-rich diet suppressed genes linked to choline, glycine, glycolysis, and TMA lyase, favoring glycine utilization via pyruvate pathways. Glycolytic inhibitors amplified these effects, mainly inhibiting pyruvate kinase. Our findings underscored the crosstalk between TMA lyase and glycolytic pathways, regulating glycine levels, and suggested avenues for targeted interventions and personalized diets to curb choline TMA lyase production.

The perinucleolar compartment and the oncogenic super-enhancers are part of the same phase-separated structure filled with phosphatidylinositol 4,5-bisphosphate and long non-coding RNA HANR

The liquid-liquid phase separation in the cell nucleus regulates various processes such as gene regulation and transcription control, chromatin organization, and DNA repair. A plethora of proteins and RNAs contribute to the formation of biomolecular condensates and recently, several nuclear phosphoinositides were shown to be a part of these membrane-less complexes within the nucleus as well. Here we lipid-interacting RNA sequencing (LIPRNAseq) and confocal microscopy to uncover the RNA-binding capacity and localization of phosphatidylinositol 4,5 bisphosphate (PIP2). We discovered the consensus PIP2-binding AU-rich RNA motif and identified long non-coding RNA HANR (lncHANR) to colocalize with PIP2 in the proximity to the nucleolus in the perinucleolar compartment (PNC). Colocalization studies with different nuclear markers reveal that PIP2-HANR presence in the PNC correlates with oncogenic super-enhancers, and both PNC and oncogenic enhancers are part of the same structure. As lncHANR, PNC, and oncogenic super-enhancers are associated with cancer cell lines and tumors, we suggest that they can serve as interchangeable prognostic markers. Understanding of the interplay between lipid metabolism, and lncRNAs in subnuclear compartment phase separation can lead to future improvement in treatment strategies and personalized cancer management approaches.

The first intron and promoter of Arabidopsis DIACYLGLYCEROL ACYLTRANSFERASE 1 exert synergistic effects on pollen and embryo lipid accumulation

Accumulation of triacylglycerols (TAGs) is crucial during various stages of plant development. In Arabidopsis, two enzymes share overlapping functions to produce TAGs, namely acyl-CoA:diacylglycerol acyltransferase 1 (DGAT1) and phospholipid:diacylglycerol acyltransferase 1 (PDAT1). Loss of function of both genes in a dgat1-1/pdat1-2 double mutant is gametophyte lethal. However, the key regulatory elements controlling tissue-specific expression of either gene has not yet been identified. We transformed a dgat1-1/dgat1-1//PDAT1/pdat1-2 parent with transgenic constructs containing the Arabidopsis DGAT1 promoter fused to the AtDGAT1 open reading frame either with or without the first intron. Triple homozygous plants were obtained, however, in the absence of the DGAT1 first intron anthers fail to fill with pollen, seed yield is c. 10% of wild-type, seed oil content remains reduced (similar to dgat1-1/dgat1-1), and non-Mendelian segregation of the PDAT1/pdat1-2 locus occurs. Whereas plants expressing the AtDGAT1pro:AtDGAT1 transgene containing the first intron mostly recover phenotypes to wild-type. This study establishes that a combination of the promoter and first intron of AtDGAT1 provides the proper context for temporal and tissue-specific expression of AtDGAT1 in pollen. Furthermore, we discuss possible mechanisms of intron mediated regulation and how regulatory elements can be used as genetic tools to functionally replace TAG biosynthetic enzymes in Arabidopsis.

Transient promoter interactions modulate developmental gene activation

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

A mathematical model clarifies the ABC Score formula used in enhancer-gene prediction

Enhancers are discrete DNA elements that regulate the expression of eukaryotic genes. They are important not only for their regulatory function, but also as loci that are frequently associated with disease traits. Despite their significance, our conceptual understanding of how enhancers work remains limited. CRISPR-interference methods have recently provided the means to systematically screen for enhancers in cell culture, from which a formula for predicting whether an enhancer regulates a gene, the Activity-by-Contact (ABC) Score, has emerged and has been widely adopted. While useful as a binary classifier, it is less effective at predicting the quantitative effect of an enhancer on gene expression. It is also unclear how the algebraic form of the ABC Score arises from the underlying molecular mechanisms and what assumptions are needed for it to hold. Here, we use the graph-theoretic linear framework, previously introduced to analyze gene regulation, to formulate the default model, a mathematical model of how multiple enhancers independently regulate a gene. We show that the algebraic form of the ABC Score arises from this model. However, the default model assumptions also imply that enhancers act additively on steady-state gene expression. This is known to be false for certain genes and we show how modifying the assumptions can accommodate this discrepancy. Overall, our approach lays a rigorous, biophysical foundation for future studies of enhancer-gene regulation.

Engineering structural variants to interrogate genome function

Structural variation, such as deletions, duplications, inversions and complex rearrangements, can have profound effects on gene expression, genome stability, phenotypic diversity and disease susceptibility. Structural variants can encompass up to millions of bases and have the potential to rearrange substantial segments of the genome. They contribute considerably more to genetic diversity in human populations and have larger effects on phenotypic traits than point mutations. Until recently, our understanding of the effects of structural variants was driven mainly by studying naturally occurring variation. New genome-engineering tools capable of generating deletions, insertions, inversions and translocations, together with the discovery of new recombinases and advances in creating synthetic DNA constructs, now enable the design and generation of an extended range of structural variation. Here, we discuss these tools and examples of their application and highlight existing challenges that will need to be overcome to fully harness their potential.

DNA regulatory element cooperation and competition in transcription

Regulation of eukaryotic transcription is a complex process that enables precise temporal and spatial control of gene expression. Promoters, which are cis-regulatory elements (CREs) located proximal to the transcription start site (TSS), selectively integrate regulatory cues from distal CREs, or enhancers, and their associated transcriptional machinery. In this review, we discuss current knowledge regarding CRE cooperation and competition impacting gene expression, including features of enhancer-promoter, enhancer-enhancer, and promoter-promoter interplay. We also provide an overview of recent insights into the underlying molecular mechanisms that facilitate physical and functional interaction of regulatory elements, such as the involvement of enhancer RNAs and biomolecular condensates. [BMB Reports 2024; 57(12): 509-520].

SE-lncRNAs in Cancer: Classification, Subcellular Localisation, Function and Corresponding TFs

Emerging evidence highlights certain long noncoding RNAs (lncRNAs) transcribed from or interacting with super-enhancer (SE) regulatory elements. These lncRNAs, known as SE-lncRNAs, are strongly linked to cancer and regulate cancer progression through multiple interactions with downstream targets. The expression of SE-lncRNAs is controlled by various transcription factors (TFs), and dysregulation of these TFs can contribute to cancer development. In this review, we discuss the characteristics, classification and subcellular distribution of SE-lncRNAs and summarise the role of key TFs in the transcription and regulation of SE-lncRNAs. Moreover, we examine the distinct functions and potential mechanisms of SE-lncRNAs in cancer progression.

Entangled and non-modular enhancer sequences producing independent spatial activities

The modularity of transcriptional enhancers is central to our understanding of morphological evolution, allowing specific changes to a gene expression pattern component, without affecting others. Enhancer modularity refers to physically separated stretches of regulatory sequence producing discrete spatiotemporal transcriptional activity. This concept stems from assays that test the sufficiency of a DNA segment to drive spatial reporter expression resembling that of the corresponding gene. Focusing on spatial patterns, it overlooks quantitative aspects of gene expression, underestimating the regulatory sequence actually required to reach full endogenous expression levels. Here, we show that five regulatory activities of the gene yellow in Drosophila, classically described as modular, result from extensively overlapping sequences, with broadly distributed regulatory information. Nevertheless, the independent regulatory activities of these entangled enhancers appear to be nucleated by specific segments that we called enhancer cores. Our work calls for a reappraisal of enhancer definition and properties, as well as of the consequences on regulatory evolution.

Investigation of Polymorphisms Induced by the Solo Long Terminal Repeats (Solo-LTRs) in Porcine Endogenous Retroviruses (ERVs)

Homologous recombination events take place between the 5' and 3' long terminal repeats (LTRs) of ERVs, resulting in the generation of solo-LTR, which can cause solo-LTR-associated polymorphism across different genomes. In the current study, specific criteria were established for the filtration of solo-LTRs, resulting in an average of 5630 solo-LTRs being identified in 21 genomes. Subsequently, a protocol was developed for detecting solo-LTR polymorphisms in the pig genomes, resulting in the discovery of 927 predicted solo-LTR polymorphic sites. Following verification and filtration processes, 603 highly reliable solo-LTR polymorphic sites were retained, involving 446 solo-LTR presence sites (solo-LTR+) and 157 solo-LTR absence sites (solo-LTR-) relative to the reference genome. Intersection analysis with gene/functional regions revealed that 248 solo-LTR- sites and 23 solo-LTR+ sites overlapped with genes or were in the vicinity of genes or functional regions, impacting a diverse range of gene structures. Moreover, through the utilization of 156 solo-LTR polymorphic sites for population genetic analysis, it was observed that these solo-LTR loci effectively clustered various breeds together, aligning with expectations and underscoring their practical utility. This study successfully established a methodology for detecting solo-LTR polymorphic sites. By applying these methods, a total of 603 high-reliability solo-LTR polymorphic sites were pinpointed, with nearly half of them being linked to genes or functional regions.

Super-enhancer-driven ZFP36L1 promotes PD-L1 expression in infiltrative gastric cancer

Gastric cancer (GC) is a major cause of cancer-related mortality worldwide. Despite the widespread recognition of tumor immunotherapy in treating unresectable GC, challenges, including ineffective immunotherapy and drug resistance, persist. Therefore, understanding the regulatory mechanisms of PD-L1, particularly in the context of super-enhancers (SEs) and zinc finger protein 36 ring finger protein-like 1 (ZFP36L1) RNA-binding protein, is crucial. In this study, we performed H3K27ac Cleavage Under Targets and Tagmentation (CUT&Tag) sequencing, investigated the heterogeneity of SEs between two GC subtypes with differential growth patterns, and revealed the immune escape signatures driven by ZFP36L1-SE in infiltrative GC through SEs inhibitors treatment. The regulation of ZFP36L1 to PD-L1 was evaluated by quantitative PCR, western blot, flow cytometry, and immunohistochemistry. Furthermore, we explored its regulatory mechanisms using a combination of molecular biology techniques, including luciferase reporter assay, GST/RNA pull-down, chromatin immunoprecipitation (ChIP)/RIP experiments, and in vivo functional assays. We demonstrated that ZFP36L1, driven by an SE, enhances IFN-γ-induced PD-L1 expression, with SPI1 identified as the specific transcription factor binding to ZFP36L1-SE. Mechanistically, ZFP36L1 binds to the adenylate uridylate-rich element in the 3' untranslated region (3'UTR) of HDAC3 mRNA, exacerbating its mRNA decay, and thereby facilitating PD-L1 abnormal transcriptional activation. Collectively, our findings provide mechanistic insights into the role of the SPI1-ZFP36L1-HDAC3-PD-L1 signaling axis in orchestrating immune escape mechanisms in GC, thereby offering valuable insights into the potential targets for immune checkpoint therapy in GC management.

Developmental and housekeeping transcriptional programs display distinct modes of enhancer-enhancer cooperativity in Drosophila

Genomic enhancers are key transcriptional regulators which, upon the binding of sequence-specific transcription factors, activate their cognate target promoters. Although enhancers have been extensively studied in isolation, a substantial number of genes have more than one simultaneously active enhancer, and it remains unclear how these cooperate to regulate transcription. Using Drosophila melanogaster S2 cells as a model, we assay the activities of more than a thousand individual enhancers and about a million enhancer pairs toward housekeeping and developmental core promoters with STARR-seq. We report that housekeeping and developmental enhancers show distinct modes of enhancer-enhancer cooperativity: while housekeeping enhancers are additive such that their combined activity mirrors the sum of their individual activities, developmental enhancers are super-additive and combine multiplicatively. Super-additivity between developmental enhancers is promiscuous and neither depends on the enhancers' endogenous genomic contexts nor on specific transcription factor motif signatures. However, it can be further boosted by Twist and Trl motifs and saturates for the highest levels of enhancer activity. These results have important implications for our understanding of gene regulation in complex multi-enhancer developmental loci and genomically clustered housekeeping genes, providing a rationale to interpret the transcriptional impact of non-coding mutations at different loci.

Context transcription factors establish cooperative environments and mediate enhancer communication

Many enhancers control gene expression by assembling regulatory factor clusters, also referred to as condensates. This process is vital for facilitating enhancer communication and establishing cellular identity. However, how DNA sequence and transcription factor (TF) binding instruct the formation of high regulatory factor environments remains poorly understood. Here we developed a new approach leveraging enhancer-centric chromatin accessibility quantitative trait loci (caQTLs) to nominate regulatory factor clusters genome-wide. By analyzing TF-binding signatures within the context of caQTLs and comparing episomal versus endogenous enhancer activities, we discovered a class of regulators, 'context-only' TFs, that amplify the activity of cell type-specific caQTL-binding TFs, that is, 'context-initiator' TFs. Similar to super-enhancers, enhancers enriched for context-only TF-binding sites display high coactivator binding and sensitivity to bromodomain-inhibiting molecules. We further show that binding sites for context-only and context-initiator TFs underlie enhancer coordination, providing a mechanistic rationale for how a loose TF syntax confers regulatory specificity.

ImAge quantitates aging and rejuvenation

For efficient, cost-effective and personalized healthcare, biomarkers that capture aspects of functional, biological aging, thus predicting disease risk and lifespan more accurately and reliably than chronological age, are essential. We developed an imaging-based chromatin and epigenetic age (ImAge) that captures intrinsic age-related trajectories of the spatial organization of chromatin and epigenetic marks in single nuclei, in mice. We show that such trajectories readily emerge as principal changes in each individual dataset without regression on chronological age, and that ImAge can be computed using several epigenetic marks and DNA labeling. We find that interventions known to affect biological aging induce corresponding effects on ImAge, including increased ImAge upon chemotherapy treatment and decreased ImAge upon caloric restriction and partial reprogramming by transient OSKM expression in liver and skeletal muscle. Further, ImAge readouts from chronologically identical mice inversely correlated with their locomotor activity, suggesting that ImAge may capture elements of biological and functional age. In sum, we developed ImAge, an imaging-based biomarker of aging with single-cell resolution rooted in the analysis of spatial organization of epigenetic marks.

Super-enhancer omics in stem cell

The hallmarks of stem cells, such as proliferation, self-renewal, development, differentiation, and regeneration, are critical to maintain stem cell identity which is sustained by genetic and epigenetic factors. Super-enhancers (SEs), which consist of clusters of active enhancers, play a central role in maintaining stemness hallmarks by specifically transcriptional model. The SE-navigated transcriptional complex, including SEs, non-coding RNAs, master transcriptional factors, Mediators and other co-activators, forms phase-separated condensates, which offers a toggle for directing diverse stem cell fate. With the burgeoning technologies of multiple-omics applied to examine different aspects of SE, we firstly raise the concept of "super-enhancer omics", inextricably linking to Pan-omics. In the review, we discuss the spatiotemporal organization and concepts of SEs, and describe links between SE-navigated transcriptional complex and stem cell features, such as stem cell identity, self-renewal, pluripotency, differentiation and development. We also elucidate the mechanism of stemness and oncogenic SEs modulating cancer stem cells via genomic and epigenetic alterations hijack in cancer stem cell. Additionally, we discuss the potential of targeting components of the SE complex using small molecule compounds, genome editing, and antisense oligonucleotides to treat SE-associated organ dysfunction and diseases, including cancer. This review also provides insights into the future of stem cell research through the paradigm of SEs.

Cell signaling and transcriptional regulation of osteoblast lineage commitment, differentiation, bone formation, and homeostasis

The initiation of osteogenesis primarily occurs as mesenchymal stem cells undergo differentiation into osteoblasts. This differentiation process plays a crucial role in bone formation and homeostasis and is regulated by two intricate processes: cell signal transduction and transcriptional gene expression. Various essential cell signaling pathways, including Wnt, BMP, TGF-β, Hedgehog, PTH, FGF, Ephrin, Notch, Hippo, and Piezo1/2, play a critical role in facilitating osteoblast differentiation, bone formation, and bone homeostasis. Key transcriptional factors in this differentiation process include Runx2, Cbfβ, Runx1, Osterix, ATF4, SATB2, and TAZ/YAP. Furthermore, a diverse array of epigenetic factors also plays critical roles in osteoblast differentiation, bone formation, and homeostasis at the transcriptional level. This review provides an overview of the latest developments and current comprehension concerning the pathways of cell signaling, regulation of hormones, and transcriptional regulation of genes involved in the commitment and differentiation of osteoblast lineage, as well as in bone formation and maintenance of homeostasis. The paper also reviews epigenetic regulation of osteoblast differentiation via mechanisms, such as histone and DNA modifications. Additionally, we summarize the latest developments in osteoblast biology spurred by recent advancements in various modern technologies and bioinformatics. By synthesizing these insights into a comprehensive understanding of osteoblast differentiation, this review provides further clarification of the mechanisms underlying osteoblast lineage commitment, differentiation, and bone formation, and highlights potential new therapeutic applications for the treatment of bone diseases.

Predicting gene expression state and prioritizing putative enhancers using 5hmC signal

BACKGROUND

Like its parent base 5-methylcytosine (5mC), 5-hydroxymethylcytosine (5hmC) is a direct epigenetic modification of cytosines in the context of CpG dinucleotides. 5hmC is the most abundant oxidized form of 5mC, generated through the action of TET dioxygenases at gene bodies of actively-transcribed genes and at active or lineage-specific enhancers. Although such enrichments are reported for 5hmC, to date, predictive models of gene expression state or putative regulatory regions for genes using 5hmC have not been developed.

RESULTS

Here, by using only 5hmC enrichment in genic regions and their vicinity, we develop neural network models that predict gene expression state across 49 cell types. We show that our deep neural network models distinguish high vs low expression state utilizing only 5hmC levels and these predictive models generalize to unseen cell types. Further, in order to leverage 5hmC signal in distal enhancers for expression prediction, we employ an Activity-by-Contact model and also develop a graph convolutional neural network model with both utilizing Hi-C data and 5hmC enrichment to prioritize enhancer-promoter links. These approaches identify known and novel putative enhancers for key genes in multiple immune cell subsets.

CONCLUSIONS

Our work highlights the importance of 5hmC in gene regulation through proximal and distal mechanisms and provides a framework to link it to genome function. With the recent advances in 6-letter DNA sequencing by short and long-read techniques, profiling of 5mC and 5hmC may be done routinely in the near future, hence, providing a broad range of applications for the methods developed here.

Multiplex profiling of developmental cis-regulatory elements with quantitative single-cell expression reporters

The inability to scalably and precisely measure the activity of developmental cis-regulatory elements (CREs) in multicellular systems is a bottleneck in genomics. Here we develop a dual RNA cassette that decouples the detection and quantification tasks inherent to multiplex single-cell reporter assays. The resulting measurement of reporter expression is accurate over multiple orders of magnitude, with a precision approaching the limit set by Poisson counting noise. Together with RNA barcode stabilization via circularization, these scalable single-cell quantitative expression reporters provide high-contrast readouts, analogous to classic in situ assays but entirely from sequencing. Screening >200 regions of accessible chromatin in a multicellular in vitro model of early mammalian development, we identify 13 (8 previously uncharacterized) autonomous and cell-type-specific developmental CREs. We further demonstrate that chimeric CRE pairs generate cognate two-cell-type activity profiles and assess gain- and loss-of-function multicellular expression phenotypes from CRE variants with perturbed transcription factor binding sites. Single-cell quantitative expression reporters can be applied in developmental and multicellular systems to quantitatively characterize native, perturbed and synthetic CREs at scale, with high sensitivity and at single-cell resolution.

Conserved Cis-Acting Range Extender Element Mediates Extreme Long-Range Enhancer Activity in Mammals

While most mammalian enhancers regulate their cognate promoters over moderate distances of tens of kilobases (kb), some enhancers act over distances in the megabase range. The sequence features enabling such extreme-distance enhancer-promoter interactions remain elusive. Here, we used in vivo enhancer replacement experiments in mice to show that short- and medium-range enhancers cannot initiate gene expression at extreme-distance range. We uncover a novel conserved cis-acting element, Range EXtender (REX), that confers extreme-distance regulatory activity and is located next to a long-range enhancer of Sall1. The REX element itself has no endogenous enhancer activity. However, addition of the REX to other short- and mid-range enhancers substantially increases their genomic interaction range. In the most extreme example observed, addition of the REX increased the range of an enhancer by an order of magnitude, from its native 71kb to 840kb. The REX element contains highly conserved [C/T]AATTA homeodomain motifs. These motifs are enriched around long-range limb enhancers genome-wide, including the ZRS, a benchmark long-range limb enhancer of Shh. Mutating the [C/T]AATTA motifs within the ZRS does not affect its limb-specific enhancer activity at short range, but selectively abolishes its long-range activity, resulting in severe limb reduction in knock-in mice. In summary, we identify a sequence signature globally associated with long-range enhancer-promoter interactions and describe a prototypical REX element that is necessary and sufficient to confer extreme-distance gene activation by remote enhancers.

Genomic context sensitizes regulatory elements to genetic disruption

Genomic context critically modulates regulatory function but is difficult to manipulate systematically. The murine insulin-like growth factor 2 (Igf2)/H19 locus is a paradigmatic model of enhancer selectivity, whereby CTCF occupancy at an imprinting control region directs downstream enhancers to activate either H19 or Igf2. We used synthetic regulatory genomics to repeatedly replace the native locus with 157-kb payloads, and we systematically dissected its architecture. Enhancer deletion and ectopic delivery revealed previously uncharacterized long-range regulatory dependencies at the native locus. Exchanging the H19 enhancer cluster with the Sox2 locus control region (LCR) showed that the H19 enhancers relied on their native surroundings while the Sox2 LCR functioned autonomously. Analysis of regulatory DNA actuation across cell types revealed that these enhancer clusters typify broader classes of context sensitivity genome wide. These results show that unexpected dependencies influence even well-studied loci, and our approach permits large-scale manipulation of complete loci to investigate the relationship between regulatory architecture and function.

Multiple genetic variants at the SLC30A8 locus affect local super-enhancer activity and influence pancreatic β-cell survival and function

Variants at the SLC30A8 locus are associated with type 2 diabetes (T2D) risk. The lead variant, rs13266634, encodes an amino acid change, Arg325Trp (R325W), at the C-terminus of the secretory granule-enriched zinc transporter, ZnT8. Although this protein-coding variant was previously thought to be the sole driver of T2D risk at this locus, recent studies have provided evidence for lowered expression of SLC30A8 mRNA in protective allele carriers. In the present study, we examined multiple variants that influence SLC30A8 allele-specific expression. Epigenomic mapping has previously identified an islet-selective enhancer cluster at the SLC30A8 locus, hosting multiple T2D risk and cASE associations, which is spatially associated with the SLC30A8 promoter and additional neighboring genes. Here, we show that deletion of variant-bearing enhancer regions using CRISPR-Cas9 in human-derived EndoC-βH3 cells lowers the expression of SLC30A8 and several neighboring genes and improves glucose-stimulated insulin secretion. While downregulation of SLC30A8 had no effect on beta cell survival, loss of UTP23, RAD21, or MED30 markedly reduced cell viability. Although eQTL or cASE analyses in human islets did not support the association between these additional genes and diabetes risk, the transcriptional regulator JQ1 lowered the expression of multiple genes at the SLC30A8 locus and enhanced stimulated insulin secretion.