Dynamic Control of Enhancer Repertoires Drives Lineage and Stage-Specific Transcription during Hematopoiesis

Endogenous small molecule effectors in GATA transcription factor mechanisms governing biological and pathological processes

Transcriptional mechanisms establish and maintain complex genetic and protein networks to control cell state transitions. The hematopoietic transcription factor GATA1 is a master regulator of erythropoiesis and megakaryopoiesis, and human GATA1 genetic variants cause anemia and megakaryoblastic leukemia. Multiomic analyses revealed that GATA1 controls expression of transporters and metabolic enzymes that dictate intracellular levels of endogenous small molecules, including heme, metal ions, and sphingolipids. Besides its canonical function as a hemoglobin component, heme facilitates or antagonizes GATA1 function to regulate erythropoiesis via mechanisms dependent or independent of the heme-binding transcription factor BTB domain and CNC homology 1 (BACH1). GATA1 regulates the expression of genes encoding heme biosynthetic enzymes and BACH1. GATA1 maintains homeostasis of bioactive ceramides during erythroid differentiation by regulating genes encoding sphingolipid metabolic enzymes. Disrupting ceramide homeostasis impairs critical cytokine signaling and is detrimental to erythroid cells. During erythroid maturation, GATA1 induces a zinc transporter switch that favors export versus import, thus dictating the intracellular zinc level, erythroblast survival, and differentiation. In aggregate, these studies support an emerging paradigm in which GATA factor-dependent transcriptional mechanisms control the intracellular levels of endogenous small molecules and small molecule-dependent feedback loops that serve as vital effectors of transcription factor activity, genome function, and cell state transitions.

Triterpenes from Ganoderma lucidum inhibit hepatocellular carcinoma by regulating enhancer-associated lncRNA in vivo

ETHNOPHARMACOLOGICAL RELEVANCE

Ganoderma lucidum (G. lucidum) has been widely used as adjuvant of anti-tumor therapy for variety tumors. The bioactive ingredients of G. lucidum mainly include triterpenes, such as Ganoderic acid A, Ganoderic acid B, Ganoderenic acid A, Ganoderenic acid B, Ganoderenic acid D, and Ganoderic acid X. However, the effects and underlying mechanisms of G. lucidum are often challenging in hepatocellular carcinoma (HCC) treatment.

AIM OF THE STUDY

To explore the potential role and mechanism of enhancer-associated lncRNAs (en-lncRNAs) in G. lucidum treated HCC through the in vivo and in vitro experiments.

MATERIALS AND METHODS

Hepa1-6-bearing C57 BL/6 mice model were established to evaluate the therapeutic efficacy of G. lucidum treated HCC. Ki67 and TUNEL staining were used to detect the tumor cell proliferation and apoptosis in vivo. The Mouse lncRNA 4*180K array was implemented to identify the differentially expressed (DE) lncRNAs and mRNAs of G. lucidum treated tumor mice. The constructed lncRNA-mRNA co-expression network and bioinformatics analysis were used to selected core en-lncRNAs and its neighboring genes. The UPLC-MS method was used to identify the triterpenes of G. lucidum, and the in vitro experiments were used to verify which triterpene monomers regulated en-lncRNAs in tumor cells. Finally, a stable knockdown/overexpression cell lines were used to confirm the relationship between en-lncRNA and neighboring gene.

RESULTS

Ki67 and TUNEL staining demonstrated G. lucidum significantly inhibited tumor growth, suppressed cell proliferation and induced apoptosis in vivo. Transcriptomic analysis revealed the existence of 126 DE lncRNAs high correlated with 454 co-expressed mRNAs in G. lucidum treated tumor mice. Based on lncRNA-mRNA network and qRT-PCR validation, 6 core lncRNAs were selected and considered high correlated with G. lucidum treatment. Bioinformatics analysis revealed FR036820 and FR121302 might act as enhancers, and qRT-PCR results suggested FR121302 might enhance Popdc2 mRNA level in HCC. Furthermore, 6 main triterpene monomers of G. lucidum were identified by UPLC-MS method, and in vitro experiments showed FR121302 and Popdc2 were significantly suppressed by Ganoderenic acid A and Ganoderenic acid B, respectively. The knock/overexpression results demonstrated that FR121302 activating and enhancing Popdc2 expression levels, and Ganoderenic acid A and Ganoderenic acid B dramatically suppressed FR121302 and decreased Popdc2 level in Hepa1-6 cells.

CONCLUSIONS

Enhancer-associated lncRNA plays a crucial role as an enhancer during hepatocarcinogenesis, and triterpenes of G. lucidum significantly inhibited tumor cell proliferation and induced apoptosis by regulating en-lncRNAs. Our study demonstrated Ganoderenic acid A and Ganoderenic acid B suppressed en-lncRNA FR121302 may be one of the critical strategies of G. lucidum inhibit hepatocellular carcinoma growth.

Concurrent loss of LKB1 and KEAP1 enhances SHMT-mediated antioxidant defence in KRAS-mutant lung cancer

Non-small-cell lung cancer (NSCLC) with concurrent mutations in KRAS and the tumour suppressor LKB1 (KL NSCLC) is refractory to most therapies and has one of the worst predicted outcomes. Here we describe a KL-induced metabolic vulnerability associated with serine-glycine-one-carbon (SGOC) metabolism. Using RNA-seq and metabolomics data from human NSCLC, we uncovered that LKB1 loss enhanced SGOC metabolism via serine hydroxymethyltransferase (SHMT). LKB1 loss, in collaboration with KEAP1 loss, activated SHMT through inactivation of the salt-induced kinase (SIK)-NRF2 axis and satisfied the increased demand for one-carbon units necessary for antioxidant defence. Chemical and genetic SHMT suppression increased cellular sensitivity to oxidative stress and cell death. Further, the SHMT inhibitor enhanced the in vivo therapeutic efficacy of paclitaxel (first-line NSCLC therapy inducing oxidative stress) in KEAP1-mutant KL tumours. The data reveal how this highly aggressive molecular subtype of NSCLC fulfills their metabolic requirements and provides insight into therapeutic strategies.

Epigenetic dissection of human blood group genes reveals regulatory elements and detailed characteristics of KEL and four other loci

BACKGROUND

Blood typing is essential for safe transfusions and is performed serologically or genetically. Genotyping predominantly focuses on coding regions, but non-coding variants may affect gene regulation, as demonstrated in the ABO, FY and XG systems. To uncover regulatory loci, we expanded a recently developed bioinformatics pipeline for discovery of non-coding variants by including additional epigenetic datasets.

METHODS

Multiple datasets including ChIP-seq with erythroid transcription factors (TFs), histone modifications (H3K27ac, H3K4me1), and chromatin accessibility (ATAC-seq) were analyzed. Candidate regulatory regions were investigated for activity (luciferase assays) and TF binding (electrophoretic mobility shift assay, EMSA, and mass spectrometry, MS).

RESULTS

In total, 814 potential regulatory sites in 47 blood-group-related genes were identified where one or more erythroid TFs bound. Enhancer candidates in CR1, EMP3, ABCB6, and ABCC4 indicated by ATAC-seq, histone markers, and co-occupancy of 4 TFs (GATA1/KLF1/RUNX1/NFE2) were investigated but only CR1 and ABCC4 showed increased transcription. Co-occupancy of GATA1 and KLF1 was observed in the KEL promoter, previously reported to contain GATA1 and Sp1 sites. TF binding energy scores decreased when three naturally occurring variants were introduced into GATA1 and KLF1 motifs. Two of three GATA1 sites and the KLF1 site were confirmed functionally. EMSA and MS demonstrated increased GATA1 and KLF1 binding to the wild-type compared to variant motifs.

DISCUSSION

This combined bioinformatics and experimental approach revealed multiple candidate regulatory regions and predicted TF co-occupancy sites. The KEL promoter was characterized in detail, indicating that two adjacent GATA1 and KLF1 motifs are most crucial for transcription.

Bidirectional epigenetic editing reveals hierarchies in gene regulation

CRISPR perturbation methods are limited in their ability to study non-coding elements and genetic interactions. In this study, we developed a system for bidirectional epigenetic editing, called CRISPRai, in which we apply activating (CRISPRa) and repressive (CRISPRi) perturbations to two loci simultaneously in the same cell. We developed CRISPRai Perturb-seq by coupling dual perturbation gRNA detection with single-cell RNA sequencing, enabling study of pooled perturbations in a mixed single-cell population. We applied this platform to study the genetic interaction between two hematopoietic lineage transcription factors, SPI1 and GATA1, and discovered novel characteristics of their co-regulation on downstream target genes, including differences in SPI1 and GATA1 occupancy at genes that are regulated through different modes. We also studied the regulatory landscape of IL2 (interleukin-2) in Jurkat T cells, primary T cells and chimeric antigen receptor (CAR) T cells and elucidated mechanisms of enhancer-mediated IL2 gene regulation. CRISPRai facilitates investigation of context-specific genetic interactions, provides new insights into gene regulation and will enable exploration of non-coding disease-associated variants.

Interspecies regulatory landscapes and elements revealed by novel joint systematic integration of human and mouse blood cell epigenomes

Knowledge of locations and activities of cis-regulatory elements (CREs) is needed to decipher basic mechanisms of gene regulation and to understand the impact of genetic variants on complex traits. Previous studies identified candidate CREs (cCREs) using epigenetic features in one species, making comparisons difficult between species. In contrast, we conducted an interspecies study defining epigenetic states and identifying cCREs in blood cell types to generate regulatory maps that are comparable between species, using integrative modeling of eight epigenetic features jointly in human and mouse in our Validated Systematic Integration (VISION) Project. The resulting catalogs of cCREs are useful resources for further studies of gene regulation in blood cells, indicated by high overlap with known functional elements and strong enrichment for human genetic variants associated with blood cell phenotypes. The contribution of each epigenetic state in cCREs to gene regulation, inferred from a multivariate regression, was used to estimate epigenetic state Regulatory Potential (esRP) scores for each cCRE in each cell type, which were used to categorize dynamic changes in cCREs. Groups of cCREs displaying similar patterns of regulatory activity in human and mouse cell types, obtained by joint clustering on esRP scores, harbored distinctive transcription factor binding motifs that were similar between species. An interspecies comparison of cCREs revealed both conserved and species-specific patterns of epigenetic evolution. Finally, we showed that comparisons of the epigenetic landscape between species can reveal elements with similar roles in regulation, even in the absence of genomic sequence alignment.

A novel cis-regulatory element regulates αD and αA-globin gene expression in chicken erythroid cells

Background

Cis-regulatory elements (CREs) play crucial roles in regulating gene expression during erythroid cell differentiation. Genome-wide erythroid-specific CREs have not been characterized in chicken erythroid cells, which is an organism model used to study epigenetic regulation during erythropoiesis.

Methods

Analysis of public genome-wide accessibility (ATAC-seq) maps, along with transcription factor (TF) motif analysis, CTCF, and RNA Pol II occupancy, as well as transcriptome analysis in fibroblasts and erythroid HD3 cells, were used to characterize erythroid-specific CREs. An α-globin CRE was identified, and its regulatory activity was validated in vitro and in vivo by luciferase activity and genome-editing assays in HD3 cells, respectively. Additionally, circular chromosome conformation capture (UMI-4C) assays were used to distinguish its role in structuring the α-globin domain in erythroid chicken cells.

Results

Erythroid-specific CREs displayed occupancy by erythroid TF binding motifs, CTCF, and RNA Pol II, as well as an association with genes involved in hematopoiesis and cell differentiation. An α-globin CRE, referred to as CRE-2, was identified as exhibiting enhancer activity over αD and αA genes in vitro and in vivo. Induction of terminal erythroid differentiation showed that α-globin CRE-2 is required for the induction of αD and αA. Analysis of TF binding motifs at α-globin CRE-2 shows apparent regulation mediated by GATA-1, YY1, and CTCF binding.

Conclusion

Our findings demonstrate that cell-specific CREs constitute a key mechanism that contributes to the fine-tuning gene regulation of erythroid cell differentiation and provide insights into the annotation and characterization of CREs in chicken cells.

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.

Three-Dimensional Chromosomal Landscape Revealing miR-146a Dysfunctional Enhancer in Lupus and Establishing a CRISPR-Mediated Approach to Inhibit the Interferon Pathway

OBJECTIVE

The diminished expression of microRNA-146a (miR-146a) in systemic lupus erythematosus (SLE) contributes to the aberrant activation of the interferon pathway. Despite its significance, the underlying mechanism driving this reduced expression remains elusive. Considering the integral role of enhancers in steering gene expression, our study seeks to pinpoint the SLE-affected enhancers responsible for modulating miR-146a expression. Additionally, we aim to elucidate the mechanisms by which these enhancers influence the contribution of miR-146a to the activation of the interferon pathway.

METHODS

Circular chromosome conformation capture sequencing and epigenomic profiles were applied to identify candidate enhancers of miR-146a. CRISPR activation was performed to screen functional enhancers. Differential analysis of chromatin accessibility was used to identify SLE-dysregulated enhancers, and the mechanism underlying enhancer dysfunction was investigated by analyzing transcription factor binding. The therapeutic value of a lupus-related enhancer was further evaluated by targeting it in the peripheral blood mononuclear cells (PBMCs) of patients with SLE through a CRISPR activation approach.

RESULTS

We identified shared and cell-specific enhancers of miR-146a in distinct immune cells. An enhancer 32.5 kb downstream of miR-146a possesses less accessibility in SLE, and its chromatin openness was negatively correlated with SLE disease activity. Moreover, CCAAT/enhancer binding protein α, a down-regulated transcription factor in patients with SLE, binds to the 32.5-kb enhancer and induces the epigenomic change of this locus. Furthermore, CRISPR-based activation of this enhancer in SLE PBMCs could inhibit the activity of interferon pathway.

CONCLUSION

Our work defines a promising target for SLE intervention. We adopted integrative approaches to define cell-specific and functional enhancers of the SLE critical gene and investigated the mechanism underlying its dysregulation mediated by a lupus-related enhancer.

Mechanisms controlling cellular and systemic iron homeostasis

In mammals, hundreds of proteins use iron in a multitude of cellular functions, including vital processes such as mitochondrial respiration, gene regulation and DNA synthesis or repair. Highly orchestrated regulatory systems control cellular and systemic iron fluxes ensuring sufficient iron delivery to target proteins is maintained, while limiting its potentially deleterious effects in iron-mediated oxidative cell damage and ferroptosis. In this Review, we discuss how cells acquire, traffick and export iron and how stored iron is mobilized for iron-sulfur cluster and haem biogenesis. Furthermore, we describe how these cellular processes are fine-tuned by the combination of various sensory and regulatory systems, such as the iron-regulatory protein (IRP)-iron-responsive element (IRE) network, the nuclear receptor co-activator 4 (NCOA4)-mediated ferritinophagy pathway, the prolyl hydroxylase domain (PHD)-hypoxia-inducible factor (HIF) axis or the nuclear factor erythroid 2-related factor 2 (NRF2) regulatory hub. We further describe how these pathways interact with systemic iron homeostasis control through the hepcidin-ferroportin axis to ensure appropriate iron fluxes. This knowledge is key for the identification of novel therapeutic opportunities to prevent diseases of cellular and/or systemic iron mismanagement.

Genomic landscape of patients with germline RUNX1 variants and familial platelet disorder with myeloid malignancy

ABSTRACT

Familial platelet disorder with associated myeloid malignancies (FPDMM) is caused by germline RUNX1 mutations and characterized by thrombocytopenia and increased risk of hematologic malignancies. We recently launched a longitudinal natural history study for patients with FPDMM. Among 27 families with research genomic data by the end of 2021, 26 different germline RUNX1 variants were detected. Besides missense mutations enriched in Runt homology domain and loss-of-function mutations distributed throughout the gene, splice-region mutations and large deletions were detected in 6 and 7 families, respectively. In 25 of 51 (49%) patients without hematologic malignancy, somatic mutations were detected in at least 1 of the clonal hematopoiesis of indeterminate potential (CHIP) genes or acute myeloid leukemia (AML) driver genes. BCOR was the most frequently mutated gene (in 9 patients), and multiple BCOR mutations were identified in 4 patients. Mutations in 6 other CHIP- or AML-driver genes (TET2, DNMT3A, KRAS, LRP1B, IDH1, and KMT2C) were also found in ≥2 patients without hematologic malignancy. Moreover, 3 unrelated patients (1 with myeloid malignancy) carried somatic mutations in NFE2, which regulates erythroid and megakaryocytic differentiation. Sequential sequencing data from 19 patients demonstrated dynamic changes of somatic mutations over time, and stable clones were more frequently found in older adult patients. In summary, there are diverse types of germline RUNX1 mutations and high frequency of somatic mutations related to clonal hematopoiesis in patients with FPDMM. Monitoring changes in somatic mutations and clinical manifestations prospectively may reveal mechanisms for malignant progression and inform clinical management. This trial was registered at www.clinicaltrials.gov as #NCT03854318.

MYB: A Key Transcription Factor in the Hematopoietic System Subject to Many Levels of Control

MYB is a master regulator and pioneer factor highly expressed in hematopoietic progenitor cells (HPCs) where it contributes to the reprogramming processes operating during hematopoietic development. MYB plays a complex role being involved in several lineages of the hematopoietic system. At the molecular level, the MYB gene is subject to intricate regulation at many levels through several enhancer and promoter elements, through transcriptional elongation control, as well as post-transcriptional regulation. The protein is modulated by post-translational modifications (PTMs) such as SUMOylation restricting the expression of its downstream targets. Together with a range of interaction partners, cooperating transcription factors (TFs) and epigenetic regulators, MYB orchestrates a fine-tuned symphony of genes expressed during various stages of haematopoiesis. At the same time, the complex MYB system is vulnerable, being a target for unbalanced control and cancer development.

The chromatin signatures of enhancers and their dynamic regulation

Enhancers are cis-regulatory elements that can stimulate gene expression from distance, and drive precise spatiotemporal gene expression profiles during development. Functional enhancers display specific features including an open chromatin conformation, Histone H3 lysine 27 acetylation, Histone H3 lysine 4 mono-methylation enrichment, and enhancer RNAs production. These features are modified upon developmental cues which impacts their activity. In this review, we describe the current state of knowledge about enhancer functions and the diverse chromatin signatures found on enhancers. We also discuss the dynamic changes of enhancer chromatin signatures, and their impact on lineage specific gene expression profiles, during development or cellular differentiation.

Biomanufacturing of glycosylated antibodies: Challenges, solutions, and future prospects

Traditionally, recombinant protein production has been done in several expression hosts of bacteria, fungi, and majorly CHO (Chinese Hamster Ovary) cells; few have high production costs and are susceptible to harmful toxin contamination. Green algae have the potential to produce recombinant proteins in a more sustainable manner. Microalgal diversity leads to offer excellent opportunities to produce glycosylated antibodies. An antibody with humanized glycans plays a crucial role in cellular communication that works to regulate cells and molecules, to control disease, and to stimulate immunity. Therefore, it becomes necessary to understand the role of abiotic factors (light, temperature, pH, etc.) in the production of bioactive molecules and molecular mechanisms of product synthesis from microalgae which would lead to harnessing the potential of algal bio-refinery. However, the potential of microalgae as the source of bio-refinery has been less explored. In the present review, omics approaches for microalgal engineering, methods of humanized glycoproteins production focusing majorly on N-glycosylation pathways, light-based regulation of glycosylation machinery, and production of antibodies with humanized glycans in microalgae with a major emphasis on modulation of post-translation machinery of microalgae which might play a role in better understanding of microalgal potential as a source for antibody production along with future perspectives.

Spatial promoter-enhancer hubs in cancer: organization, regulation, and function

Transcriptional dysregulation is a hallmark of cancer and can be driven by altered enhancer landscapes. Recent studies in genome organization have revealed that multiple enhancers and promoters can spatially coalesce to form dynamic topological assemblies, known as promoter-enhancer hubs, which strongly correlate with elevated gene expression. In this review, we discuss the structure and complexity of promoter-enhancer hubs recently identified in multiple cancer types. We further discuss underlying mechanisms driving dysregulation of promoter-enhancer hubs and speculate on their functional role in pathogenesis. Understanding the role of promoter-enhancer hubs in transcriptional dysregulation can provide insight into new therapeutic approaches to target these complex features of genome organization.

Iron Oxide Nanoparticles with and without Cobalt Functionalization Provoke Changes in the Transcription Profile via Epigenetic Modulation of Enhancer Activity

Despite the progress in the field of nanotoxicology, much about the cellular mechanisms that mediate the adverse effects of nanoparticles (NPs) and, in particular, the possible role of epigenetics in nanotoxicity, remains to be clarified. Therefore, we studied the changes occurring in the genome-wide distribution of H3K27ac, H3K4me1, H3K9me2, and H3K27me3 histone modifications and compared them with the transcriptome after exposing NIH3T3 cells to iron-based magnetic NPs (i.e., Fe2O3 and Fe2O3@Co NPs). We found that the transcription response is mainly due to changes in the genomic distribution of H3K27ac that can modulate the activity of enhancers. We propose that alteration of the epigenetic landscape is a key mechanism in defining the gene expression program changes resulting in nanotoxicity. With this approach, it is possible to construct a data set of genomic regions that could be useful for defining toxicity in a manner that is more comprehensive than what is possible with the present toxicology assays.

Understanding blood development and leukemia using sequencing-based technologies and human cell systems

Our current understanding of human hematopoiesis has undergone significant transformation throughout the years, challenging conventional views. The evolution of high-throughput technologies has enabled the accumulation of diverse data types, offering new avenues for investigating key regulatory processes in blood cell production and disease. In this review, we will explore the opportunities presented by these advancements for unraveling the molecular mechanisms underlying normal and abnormal hematopoiesis. Specifically, we will focus on the importance of enhancer-associated regulatory networks and highlight the crucial role of enhancer-derived transcription regulation. Additionally, we will discuss the unprecedented power of single-cell methods and the progression in using in vitro human blood differentiation system, in particular induced pluripotent stem cell models, in dissecting hematopoietic processes. Furthermore, we will explore the potential of ever more nuanced patient profiling to allow precision medicine approaches. Ultimately, we advocate for a multiparameter, regulatory network-based approach for providing a more holistic understanding of normal hematopoiesis and blood disorders.

Role of the pioneer transcription factor GATA2 in health and disease

The transcription factor GATA2 is involved in human diseases ranging from hematopoietic disorders, to cancer, to infectious diseases. GATA2 is one of six GATA-family transcription factors that act as pioneering transcription factors which facilitate the opening of heterochromatin and the subsequent binding of other transcription factors to induce gene expression from previously inaccessible regions of the genome. Although GATA2 is essential for hematopoiesis and lymphangiogenesis, it is also expressed in other tissues such as the lung, prostate gland, gastrointestinal tract, central nervous system, placenta, fetal liver, and fetal heart. Gene or transcriptional abnormalities of GATA2 causes or predisposes patients to several diseases including the hematological cancers acute myeloid leukemia and acute lymphoblastic leukemia, the primary immunodeficiency MonoMAC syndrome, and to cancers of the lung, prostate, uterus, kidney, breast, gastric tract, and ovaries. Recent data has also linked GATA2 expression and mutations to responses to infectious diseases including SARS-CoV-2 and Pneumocystis carinii pneumonia, and to inflammatory disorders such as atherosclerosis. In this article we review the role of GATA2 in the etiology and progression of these various diseases.

Insight into the regulatory mechanism of dynamic chromatin 3D interactions during cardiomyocyte differentiation in human

Cardiogenesis is an extremely complicated process involved with DNA regulatory elements, and trans factors regulate gene expression pattern spatiotemporally. Enhancers, as the well-known DNA elements, activate target gene expression by transcription factors (TFs) occupied to organize dynamic three-dimensional (3D) interactions, which when affected or interrupted might cause heart defects or diseases. In this study, we integrated transcriptome, 3D genome, and regulatome to reorganize the global 3D genome in cardiomyogenesis, showing a gradually decreased trend of both chromatin interactions and topological associating domains (TADs) during cardiomyocyte differentiation. And almost all of the chromatin interactions occurred within the same or between adjacent TADs involved with enhancers, indicating that dynamical rewiring of enhancer-related chromatin interactions in the continuous expansive TADs is closely correlated to cardiogenesis. Moreover, we found stage-specific interactions activate stage-specific expression to be involved within corresponding biological functions, and the stage-specific combined regulations of enhancers and binding TFs form connected networks to control stage-specific expression and biological processes, which promote cardiomyocyte differentiation. Finally, we identified markers based on regulatory networks, which might drive cardiac development. This study demonstrates the power of enhancer interactome combined with active TFs to reveal insights into transcriptional regulatory networks during cardiomyogenesis.

Dictys: dynamic gene regulatory network dissects developmental continuum with single-cell multiomics

Gene regulatory networks (GRNs) are key determinants of cell function and identity and are dynamically rewired during development and disease. Despite decades of advancement, challenges remain in GRN inference, including dynamic rewiring, causal inference, feedback loop modeling and context specificity. To address these challenges, we develop Dictys, a dynamic GRN inference and analysis method that leverages multiomic single-cell assays of chromatin accessibility and gene expression, context-specific transcription factor footprinting, stochastic process network and efficient probabilistic modeling of single-cell RNA-sequencing read counts. Dictys improves GRN reconstruction accuracy and reproducibility and enables the inference and comparative analysis of context-specific and dynamic GRNs across developmental contexts. Dictys' network analyses recover unique insights in human blood and mouse skin development with cell-type-specific and dynamic GRNs. Its dynamic network visualizations enable time-resolved discovery and investigation of developmental driver transcription factors and their regulated targets. Dictys is available as a free, open-source and user-friendly Python package.

Differentiation block in acute myeloid leukemia regulated by intronic sequences of FTO

Iroquois transcription factor gene IRX3 is highly expressed in 20-30% of acute myeloid leukemia (AML) and contributes to the pathognomonic differentiation block. Intron 8 FTO sequences ∼220kB downstream of IRX3 exhibit histone acetylation, DNA methylation, and contacts with the IRX3 promoter, which correlate with IRX3 expression. Deletion of these intronic elements confirms a role in positively regulating IRX3. RNAseq revealed long non-coding (lnc) transcripts arising from this locus. FTO-lncAML knockdown (KD) induced differentiation of AML cells, loss of clonogenic activity, and reduced FTO intron 8:IRX3 promoter contacts. While both FTO-lncAML KD and IRX3 KD induced differentiation, FTO-lncAML but not IRX3 KD led to HOXA downregulation suggesting transcript activity in trans. FTO-lncAMLhigh AML samples expressed higher levels of HOXA and lower levels of differentiation genes. Thus, a regulatory module in FTO intron 8 consisting of clustered enhancer elements and a long non-coding RNA is active in human AML, impeding myeloid differentiation.

Multidimensional profiling reveals GATA1-modulated stage-specific chromatin states and functional associations during human erythropoiesis

Mammalian erythroid development can be divided into three stages: hematopoietic stem and progenitor cell (HSPC), erythroid progenitor (Ery-Pro), and erythroid precursor (Ery-Pre). However, the mechanisms by which the 3D genome changes to establish the stage-specific transcription programs that are critical for erythropoiesis remain unclear. Here, we analyze the chromatin landscape at multiple levels in defined populations from primary human erythroid culture. While compartments and topologically associating domains remain largely unchanged, ∼50% of H3K27Ac-marked enhancers are dynamic in HSPC versus Ery-Pre. The enhancer anchors of enhancer-promoter loops are enriched for occupancy of respective stage-specific transcription factors (TFs), indicating these TFs orchestrate the enhancer connectome rewiring. The master TF of erythropoiesis, GATA1, is found to occupy most erythroid gene promoters at the Ery-Pro stage, and mediate conspicuous local rewiring through acquiring binding at the distal regions in Ery-Pre, promoting productive erythroid transcription output. Knocking out GATA1 binding sites precisely abrogates local rewiring and corresponding gene expression. Interestingly, knocking down GATA1 can transiently revert the cell state to an earlier stage and prolong the window of progenitor state. This study reveals mechanistic insights underlying chromatin rearrangements during development by integrating multidimensional chromatin landscape analyses to associate with transcription output and cellular states.

Erythroid lineage chromatin accessibility maps facilitate identification and validation of NFIX as a fetal hemoglobin repressor

Human genetics has validated de-repression of fetal gamma globin (HBG) in adult erythroblasts as a powerful therapeutic paradigm in diseases involving defective adult beta globin (HBB)¹. To identify factors involved in the switch from HBG to HBB expression, we performed Assay for Transposase Accessible Chromatin with high-throughput sequencing (ATAC-seq)² on sorted erythroid lineage cells derived from bone marrow (BM) or cord blood (CB), representing adult and fetal states, respectively. BM to CB cell ATAC-seq profile comparisons revealed genome-wide enrichment of NFI DNA binding motifs and increased NFIX promoter chromatin accessibility, suggesting that NFIX may repress HBG. NFIX knockdown in BM cells increased HBG mRNA and fetal hemoglobin (HbF) protein levels, coincident with increased chromatin accessibility and decreased DNA methylation at the HBG promoter. Conversely, overexpression of NFIX in CB cells reduced HbF levels. Identification and validation of NFIX as a new target for HbF activation has implications in the development of therapeutics for hemoglobinopathies.

Super-enhancer-associated gene CAPG promotes AML progression

Acute myeloid leukemia is the most common acute leukemia in adults, the barrier of refractory and drug resistance has yet to be conquered in the clinical. Abnormal gene expression and epigenetic changes play an important role in pathogenesis and treatment. A super-enhancer is an epigenetic modifier that promotes pro-tumor genes and drug resistance by activating oncogene transcription. Multi-omics integrative analysis identifies the super-enhancer-associated gene CAPG and its high expression level was correlated with poor prognosis in AML. CAPG is a cytoskeleton protein but has an unclear function in AML. Here we show the molecular function of CAPG in regulating NF-κB signaling pathway by proteomic and epigenomic analysis. Knockdown of Capg in the AML murine model resulted in exhausted AML cells and prolonged survival of AML mice. In conclusion, SEs-associated gene CAPG can contributes to AML progression through NF-κB.

3D enhancer-promoter interactions and multi-connected hubs: Organizational principles and functional roles

The spatiotemporal control of gene expression is dependent on the activity of cis-acting regulatory sequences, called enhancers, which regulate target genes over variable genomic distances and, often, by skipping intermediate promoters, suggesting mechanisms that control enhancer-promoter communication. Recent genomics and imaging technologies have revealed highly complex enhancer-promoter interaction networks, whereas advanced functional studies have started interrogating the forces behind the physical and functional communication among multiple enhancers and promoters. In this review, we first summarize our current understanding of the factors involved in enhancer-promoter communication, with a particular focus on recent papers that have revealed new layers of complexities to old questions. In the second part of the review, we focus on a subset of highly connected enhancer-promoter "hubs" and discuss their potential functions in signal integration and gene regulation, as well as the putative factors that might determine their dynamics and assembly.

High-throughput methods for the analysis of transcription factors and chromatin modifications: Low input, single cell and spatial genomic technologies

Genome-wide analysis of transcription factors and epigenomic features is instrumental to shed light on DNA-templated regulatory processes such as transcription, cellular differentiation or to monitor cellular responses to environmental cues. Two decades of technological developments have led to a rich set of approaches progressively pushing the limits of epigenetic profiling towards single cells. More recently, disruptive technologies using innovative biochemistry came into play. Assays such as CUT&RUN, CUT&Tag and variations thereof show considerable potential to survey multiple TFs or histone modifications in parallel from a single experiment and in native conditions. These are in the path to become the dominant assays for genome-wide analysis of TFs and chromatin modifications in bulk, single-cell, and spatial genomic applications. The principles together with pros and cons are discussed.

Experimental Validation and Prediction of Super-Enhancers: Advances and Challenges

Super-enhancers (SEs) are cis-regulatory elements of the human genome that have been widely discussed since the discovery and origin of the term. Super-enhancers have been shown to be strongly associated with the expression of genes crucial for cell differentiation, cell stability maintenance, and tumorigenesis. Our goal was to systematize research studies dedicated to the investigation of structure and functions of super-enhancers as well as to define further perspectives of the field in various applications, such as drug development and clinical use. We overviewed the fundamental studies which provided experimental data on various pathologies and their associations with particular super-enhancers. The analysis of mainstream approaches for SE search and prediction allowed us to accumulate existing data and propose directions for further algorithmic improvements of SEs' reliability levels and efficiency. Thus, here we provide the description of the most robust algorithms such as ROSE, imPROSE, and DEEPSEN and suggest their further use for various research and development tasks. The most promising research direction, which is based on topic and number of published studies, are cancer-associated super-enhancers and prospective SE-targeted therapy strategies, most of which are discussed in this review.

Comparative genomics analyses reveal sequence determinants underlying interspecies variations in injury-responsive enhancers

BACKGROUND

Injury induces profound transcriptional remodeling events, which could lead to only wound healing, partial tissue repair, or perfect regeneration in different species. Injury-responsive enhancers (IREs) are cis-regulatory elements activated in response to injury signals, and have been demonstrated to promote tissue regeneration in some organisms such as zebrafish and flies. However, the functional significances of IREs in mammals remain elusive. Moreover, whether the transcriptional responses elicited by IREs upon injury are conserved or specialized in different species, and what sequence features may underlie the functional variations of IREs have not been elucidated.

RESULTS

We identified a set of IREs that are activated in both regenerative and non-regenerative neonatal mouse hearts upon myocardial ischemia-induced damage by integrative epigenomic and transcriptomic analyses. Motif enrichment analysis showed that AP-1 and ETS transcription factor binding motifs are significantly enriched in both zebrafish and mouse IREs. However, the IRE-associated genes vary considerably between the two species. We further found that the IRE-related sequences in zebrafish and mice diverge greatly, with the loss of IRE inducibility accompanied by a reduction in AP-1 and ETS motif frequencies. The functional turnover of IREs between zebrafish and mice is correlated with changes in transcriptional responses of the IRE-associated genes upon injury. Using mouse cardiomyocytes as a model, we demonstrated that the reduction in AP-1 or ETS motif frequency attenuates the activation of IREs in response to hypoxia-induced damage.

CONCLUSIONS

By performing comparative genomics analyses on IREs, we demonstrated that inter-species variations in AP-1 and ETS motifs may play an important role in defining the functions of enhancers during injury response. Our findings provide important insights for understanding the molecular mechanisms of transcriptional remodeling in response to injury across species.

Light-activated macromolecular phase separation modulates transcription by reconfiguring chromatin interactions

Biomolecular condensates participate in the regulation of gene transcription, yet the relationship between nuclear condensation and transcriptional activation remains elusive. Here, we devised a biotinylated CRISPR-dCas9-based optogenetic method, light-activated macromolecular phase separation (LAMPS), to enable inducible formation, affinity purification, and multiomic dissection of nuclear condensates at the targeted genomic loci. LAMPS-induced condensation at enhancers and promoters activates endogenous gene transcription by chromatin reconfiguration, causing increased chromatin accessibility and de novo formation of long-range chromosomal loops. Proteomic profiling of light-induced condensates by dCas9-mediated affinity purification uncovers multivalent interaction-dependent remodeling of macromolecular composition, resulting in the selective enrichment of transcriptional coactivators and chromatin structure proteins. Our findings support a model whereby the formation of nuclear condensates at native genomic loci reconfigures chromatin architecture and multiprotein assemblies to modulate gene transcription. Hence, LAMPS facilitates mechanistic interrogation of the relationship between nuclear condensation, genome structure, and gene transcription in living cells.

An intellectual disability-related MED23 mutation dysregulates gene expression by altering chromatin conformation and enhancer activities

Transcriptional Mediator controls diverse gene programs for various developmental and pathological processes. The human Mediator MED23/R617Q mutation was reported in a familial intellectual disability (ID) disorder, although the underlying mechanisms remain poorly understood. Constructed by gene editing, the Med23/R617Q knock-in mutant mice exhibited embryonic lethality due to the largely reduced Med23/R617Q protein level, but the R617Q mutation in HEK293T cells didn't change its expression and incorporation into Mediator Complex. RNA-seq revealed that MED23/R617Q mutation disturbed gene expression, related to neural development, learning and memory. Specifically, R617Q mutation reduced the MED23-dependent activities of ELK1 and E1A, but in contrast, upregulated the MAPK/ELK1-driven early immediate genes (IEGs) JUN and FOS. ChIP-seq and Hi-C revealed that the MED23 R617Q mutation reprogramed a subset of enhancers and local chromatin interactions, which correlated well with the corresponding gene expression. Importantly, the enhancers and chromatin interactions surrounding IEGs were unchanged by the R617Q mutation, but DACH1, an upstream repressor of IEGs, showed reduced enhancer-promoter interactions and decreased expression in mutant cells, thus relieving its inhibition to the intellectual-related IEGs. Overall, unraveling the MED23-DACH1-IEG axis provides a mechanistic explanation for the effects of the MED23/R617Q mutation on gene dysregulation and inherited ID.

eRNA profiling uncovers the enhancer landscape of oesophageal adenocarcinoma and reveals new deregulated pathways

Cancer is driven by both genetic and epigenetic changes that impact on gene expression profiles and the resulting tumourigenic phenotype. Enhancers are transcriptional regulatory elements that are key to our understanding of how this rewiring of gene expression is achieved in cancer cells. Here, we have harnessed the power of RNA-seq data from hundreds of patients with oesophageal adenocarcinoma (OAC) or its precursor state Barrett's oesophagus coupled with open chromatin maps to identify potential enhancer RNAs and their associated enhancer regions in this cancer. We identify ~1000 OAC-specific enhancers and use these data to uncover new cellular pathways that are operational in OAC. Among these are enhancers for JUP, MYBL2, and CCNE1, and we show that their activity is required for cancer cell viability. We also demonstrate the clinical utility of our dataset for identifying disease stage and patient prognosis. Our data therefore identify an important set of regulatory elements that enhance our molecular understanding of OAC and point to potential new therapeutic directions.

The aberrant epigenome of DNMT3B-mutated ICF1 patient iPSCs is amenable to correction, with the exception of a subset of regions with H3K4me3- and/or CTCF-based epigenetic memory

Bi-allelic hypomorphic mutations in DNMT3B disrupt DNA methyltransferase activity and lead to immunodeficiency, centromeric instability, facial anomalies syndrome, type 1 (ICF1). Although several ICF1 phenotypes have been linked to abnormally hypomethylated repetitive regions, the unique genomic regions responsible for the remaining disease phenotypes remain largely uncharacterized. Here we explored two ICF1 patient-derived induced pluripotent stem cells (iPSCs) and their CRISPR-Cas9-corrected clones to determine whether DNMT3B correction can globally overcome DNA methylation defects and related changes in the epigenome. Hypomethylated regions throughout the genome are highly comparable between ICF1 iPSCs carrying different DNMT3B variants, and significantly overlap with those in ICF1 patient peripheral blood and lymphoblastoid cell lines. These regions include large CpG island domains, as well as promoters and enhancers of several lineage-specific genes, in particular immune-related, suggesting that they are premarked during early development. CRISPR-corrected ICF1 iPSCs reveal that the majority of phenotype-related hypomethylated regions reacquire normal DNA methylation levels following editing. However, at the most severely hypomethylated regions in ICF1 iPSCs, which also display the highest increases in H3K4me3 levels and/or abnormal CTCF binding, the epigenetic memory persists, and hypomethylation remains uncorrected. Overall, we demonstrate that restoring the catalytic activity of DNMT3B can reverse the majority of the aberrant ICF1 epigenome. However, a small fraction of the genome is resilient to this rescue, highlighting the challenge of reverting disease states that are due to genome-wide epigenetic perturbations. Uncovering the basis for the persistent epigenetic memory will promote the development of strategies to overcome this obstacle.

Complexity of enhancer networks predicts cell identity and disease genes revealed by single-cell multi-omics analysis

Many enhancers exist as clusters in the genome and control cell identity and disease genes; however, the underlying mechanism remains largely unknown. Here, we introduce an algorithm, eNet, to build enhancer networks by integrating single-cell chromatin accessibility and gene expression profiles. The complexity of enhancer networks is assessed by two metrics: the number of enhancers and the frequency of predicted enhancer interactions (PEIs) based on chromatin co-accessibility. We apply eNet algorithm to a human blood dataset and find cell identity and disease genes tend to be regulated by complex enhancer networks. The network hub enhancers (enhancers with frequent PEIs) are the most functionally important. Compared with super-enhancers, enhancer networks show better performance in predicting cell identity and disease genes. eNet is robust and widely applicable in various human or mouse tissues datasets. Thus, we propose a model of enhancer networks containing three modes: Simple, Multiple and Complex, which are distinguished by their complexity in regulating gene expression. Taken together, our work provides an unsupervised approach to simultaneously identify key cell identity and disease genes and explore the underlying regulatory relationships among enhancers in single cells.

Genomic Landscape of Patients with Germline RUNX1 Variants and Familial Platelet Disorder with Myeloid Malignancy

Germline RUNX1 mutations lead to familial platelet disorder with associated myeloid malignancies (FPDMM), which is characterized by thrombocytopenia and a life-long risk (35-45%) of hematological malignancies. We recently launched a longitudinal natural history study for patients with FPDMM at the NIH Clinical Center. Among 29 families with research genomic data, 28 different germline RUNX1 variants were detected. Besides missense mutations enriched in Runt homology domain and loss-of-function mutations distributed throughout the gene, splice-region mutations and large deletions were detected in 6 and 7 families, respectively. In 24 of 54 (44.4%) non-malignant patients, somatic mutations were detected in at least one of the clonal hematopoiesis of indeterminate potential (CHIP) genes or acute myeloid leukemia (AML) driver genes. BCOR was the most frequently mutated gene (in 9 patients), and multiple BCOR mutations were identified in 4 patients. Mutations in 7 other CHIP or AML driver genes ( DNMT3A, TET2, NRAS, SETBP1, SF3B1, KMT2C , and LRP1B ) were also found in more than one non-malignant patient. Moreover, three unrelated patients (one with myeloid malignancy) carried somatic mutations in NFE2 , which regulates erythroid and megakaryocytic differentiation. Sequential sequencing data from 19 patients demonstrated dynamic changes of somatic mutations over time, and stable clones were more frequently found in elderly patients. In summary, there are diverse types of germline RUNX1 mutations and high frequency of somatic mutations related to clonal hematopoiesis in patients with FPDMM. Monitoring dynamic changes of somatic mutations prospectively will benefit patients’ clinical management and reveal mechanisms for progression to myeloid malignancies.

Key Points

Comprehensive genomic profile of patients with FPDMM with germline RUNX1 mutations. Rising clonal hematopoiesis related secondary mutations that may lead to myeloid malignancies.

Single-cell chromatin accessibility identifies enhancer networks driving gene expression during spinal cord development in mouse

Spinal cord development is precisely orchestrated by spatiotemporal gene regulatory programs. However, the underlying epigenetic mechanisms remain largely elusive. Here, we profiled single-cell chromatin accessibility landscapes in mouse neural tubes spanning embryonic days 9.5-13.5. We identified neuronal-cell-cluster-specific cis-regulatory elements in neural progenitors and neurons. Furthermore, we applied a novel computational method, eNet, to build enhancer networks by integrating single-cell chromatin accessibility and gene expression data and identify the hub enhancers within enhancer networks. It was experimentally validated in vivo for Atoh1 that knockout of the hub enhancers, but not the non-hub enhancers, markedly decreased Atoh1 expression and reduced dp1/dI1 cells. Together, our work provides insights into the epigenetic regulation of spinal cord development and a proof-of-concept demonstration of enhancer networks as a general mechanism in transcriptional regulation.

A Myb enhancer-guided analysis of basophil and mast cell differentiation

The transcription factor MYB is a crucial regulator of hematopoietic stem and progenitor cells. However, the nature of lineage-specific enhancer usage of the Myb gene is largely unknown. We identify the Myb -68 enhancer, a regulatory element which marks basophils and mast cells. Using the Myb -68 enhancer activity, we show a population of granulocyte-macrophage progenitors with higher potential to differentiate into basophils and mast cells. Single cell RNA-seq demonstrates the differentiation trajectory is continuous from progenitors to mature basophils in vivo, characterizes bone marrow cells with a gene signature of mast cells, and identifies LILRB4 as a surface marker of basophil maturation. Together, our study leads to a better understanding of how MYB expression is regulated in a lineage-associated manner, and also shows how a combination of lineage-related reporter mice and single-cell transcriptomics can overcome the rarity of target cells and enhance our understanding of gene expression programs that control cell differentiation in vivo.

Mitoferrin, Cellular and Mitochondrial Iron Homeostasis

Iron is essential for many cellular processes, but cellular iron homeostasis must be maintained to ensure the balance of cellular signaling processes and prevent disease. Iron transport in and out of the cell and cellular organelles is crucial in this regard. The transport of iron into the mitochondria is particularly important, as heme and the majority of iron-sulfur clusters are synthesized in this organelle. Iron is also required for the production of mitochondrial complexes that contain these iron-sulfur clusters and heme. As the principal iron importers in the mitochondria of human cells, the mitoferrins have emerged as critical regulators of cytosolic and mitochondrial iron homeostasis. Here, we review the discovery and structure of the mitoferrins, as well as the significance of these proteins in maintaining cytosolic and mitochondrial iron homeostasis for the prevention of cancer and many other diseases.

Super-enhancers conserved within placental mammals maintain stem cell pluripotency

Despite pluripotent stem cells sharing key transcription factors, their maintenance involves distinct genetic inputs. Emerging evidence suggests that super-enhancers (SEs) can function as master regulatory hubs to control cell identity and pluripotency in humans and mice. However, whether pluripotency-associated SEs share an evolutionary origin in mammals remains elusive. Here, we performed comprehensive comparative epigenomic and transcription factor binding analyses among pigs, humans, and mice to identify pluripotency-associated SEs. Like typical enhancers, SEs displayed rapid evolution in mammals. We showed that BRD4 is an essential and conserved activator for mammalian pluripotency-associated SEs. Comparative motif enrichment analysis revealed 30 shared transcription factor binding motifs among the three species. The majority of transcriptional factors that bind to identified motifs are known regulators associated with pluripotency. Further, we discovered three pluripotency-associated SEs (SE-SOX2, SE-PIM1, and SE-FGFR1) that displayed remarkable conservation in placental mammals and were sufficient to drive reporter gene expression in a pluripotency-dependent manner. Disruption of these conserved SEs through the CRISPR-Cas9 approach severely impaired stem cell pluripotency. Our study provides insights into the understanding of conserved regulatory mechanisms underlying the maintenance of pluripotency as well as species-specific modulation of the pluripotency-associated regulatory networks in mammals.

Nested epistasis enhancer networks for robust genome regulation

Mammalian genomes possess multiple enhancers spanning an ultralong distance (>megabases) to modulate important genes, yet it is unclear how these enhancers coordinate to achieve this task. Here, we combine multiplexed CRISPRi screening with machine learning to define quantitative enhancer-enhancer interactions. We find that the ultralong distance enhancer network possesses a nested multi-layer architecture that confers functional robustness of gene expression. Experimental characterization reveals that enhancer epistasis is maintained by three-dimensional chromosomal interactions and BRD4 condensation. Machine learning prediction of synergistic enhancers provides an effective strategy to identify non-coding variant pairs associated with pathogenic genes in diseases beyond Genome-Wide Association Studies (GWAS) analysis. Our work unveils nested epistasis enhancer networks, which can better explain enhancer functions within cells and in diseases.

HDAC1 and PRC2 mediate combinatorial control in SPI1/PU.1-dependent gene repression in murine erythroleukaemia

Although originally described as transcriptional activator, SPI1/PU.1, a major player in haematopoiesis whose alterations are associated with haematological malignancies, has the ability to repress transcription. Here, we investigated the mechanisms underlying gene repression in the erythroid lineage, in which SPI1 exerts an oncogenic function by blocking differentiation. We show that SPI1 represses genes by binding active enhancers that are located in intergenic or gene body regions. HDAC1 acts as a cooperative mediator of SPI1-induced transcriptional repression by deacetylating SPI1-bound enhancers in a subset of genes, including those involved in erythroid differentiation. Enhancer deacetylation impacts on promoter acetylation, chromatin accessibility and RNA pol II occupancy. In addition to the activities of HDAC1, polycomb repressive complex 2 (PRC2) reinforces gene repression by depositing H3K27me3 at promoter sequences when SPI1 is located at enhancer sequences. Moreover, our study identified a synergistic relationship between PRC2 and HDAC1 complexes in mediating the transcriptional repression activity of SPI1, ultimately inducing synergistic adverse effects on leukaemic cell survival. Our results highlight the importance of the mechanism underlying transcriptional repression in leukemic cells, involving complex functional connections between SPI1 and the epigenetic regulators PRC2 and HDAC1.

Erythroid Cell Research: 3D Chromatin, Transcription Factors and Beyond

Studies of the regulatory networks and signals controlling erythropoiesis have brought important insights in several research fields of biology and have been a rich source of discoveries with far-reaching implications beyond erythroid cells biology. The aim of this review is to highlight key recent discoveries and show how studies of erythroid cells bring forward novel concepts and refine current models related to genome and 3D chromatin organization, signaling and disease, with broad interest in life sciences.

Dual function NFI factors control fetal hemoglobin silencing in adult erythroid cells

The mechanisms by which the fetal-type β-globin-like genes HBG1 and HBG2 are silenced in adult erythroid precursor cells remain a fundamental question in human biology and have therapeutic relevance to sickle cell disease (SCD) and β-thalassemia. Here, we identify via a CRISPR-Cas9 genetic screen two members of the NFI transcription factor family – NFIA and NFIX – as HBG1/2 repressors. NFIA and NFIX are expressed at elevated levels in adult erythroid cells compared to fetal cells, and function cooperatively to repress HBG1/2 in cultured cells and in human-to-mouse xenotransplants. Genomic profiling, genome editing, and DNA binding assays demonstrate that the potent concerted activity of NFIA and NFIX is explained in part by their ability to stimulate the expression of BCL11A, a known silencer of the HBG1/2 genes, and in part by directly repressing the HBG1/2 genes. Thus, NFI factors emerge as versatile regulators of the fetal-to-adult switch in β-globin production.

Molecular and cellular mechanisms that regulate human erythropoiesis

To enable effective oxygen transport, ∼200 billion red blood cells (RBCs) need to be produced every day in the bone marrow through the fine-tuned process of erythropoiesis. Erythropoiesis is regulated at multiple levels to ensure that defective RBC maturation or overproduction can be avoided. Here, we provide an overview of different layers of this control, ranging from cytokine signaling mechanisms that enable extrinsic regulation of RBC production to intrinsic transcriptional pathways necessary for effective erythropoiesis. Recent studies have also elucidated the importance of posttranscriptional regulation and highlighted additional gatekeeping mechanisms necessary for effective erythropoiesis. We additionally discuss the insights gained by studying human genetic variation affecting erythropoiesis and highlight the discovery of BCL11A as a regulator of hemoglobin switching through genetic studies. Finally, we provide an outlook of how our ability to measure multiple facets of this process at single-cell resolution, while accounting for the impact of human variation, will continue to refine our knowledge of erythropoiesis and how this process is perturbed in disease. As we learn more about this intricate and important process, additional opportunities to modulate erythropoiesis for therapeutic purposes will undoubtedly emerge.

Developmental maturation of the hematopoietic system controlled by a Lin28b-let-7-Cbx2 axis

Hematopoiesis changes over life to meet the demands of maturation and aging. Here, we find that the definitive hematopoietic stem and progenitor cell (HSPC) compartment is remodeled from gestation into adulthood, a process regulated by the heterochronic Lin28b/let-7 axis. Native fetal and neonatal HSPCs distribute with a pro-lymphoid/erythroid bias with a shift toward myeloid output in adulthood. By mining transcriptomic data comparing juvenile and adult HSPCs and reconstructing coordinately activated gene regulatory networks, we uncover the Polycomb repressor complex 1 (PRC1) component Cbx2 as an effector of Lin28b/let-7's control of hematopoietic maturation. We find that juvenile Cbx2-/- hematopoietic tissues show impairment of B-lymphopoiesis, a precocious adult-like myeloid bias, and that Cbx2/PRC1 regulates developmental timing of expression of key hematopoietic transcription factors. These findings define a mechanism of regulation of HSPC output via chromatin modification as a function of age with potential impact on age-biased pediatric and adult blood disorders.

Domain adaptive neural networks improve cross-species prediction of transcription factor binding

The intrinsic DNA sequence preferences and cell type–specific cooperative partners of transcription factors (TFs) are typically highly conserved. Hence, despite the rapid evolutionary turnover of individual TF binding sites, predictive sequence models of cell type–specific genomic occupancy of a TF in one species should generalize to closely matched cell types in a related species. To assess the viability of cross-species TF binding prediction, we train neural networks to discriminate ChIP-seq peak locations from genomic background and evaluate their performance within and across species. Cross-species predictive performance is consistently worse than within-species performance, which we show is caused in part by species-specific repeats. To account for this domain shift, we use an augmented network architecture to automatically discourage learning of training species–specific sequence features. This domain adaptation approach corrects for prediction errors on species-specific repeats and improves overall cross-species model performance. Our results show that cross-species TF binding prediction is feasible when models account for domain shifts driven by species-specific repeats.

Identification of a modular super-enhancer in murine retinal development

Super-enhancers are expansive regions of genomic DNA comprised of multiple putative enhancers that contribute to the dynamic gene expression patterns during development. This is particularly important in neurogenesis because many essential transcription factors have complex developmental stage- and cell-type specific expression patterns across the central nervous system. In the developing retina, Vsx2 is expressed in retinal progenitor cells and is maintained in differentiated bipolar neurons and Müller glia. A single super-enhancer controls this complex and dynamic pattern of expression. Here we show that deletion of one region disrupts retinal progenitor cell proliferation but does not affect cell fate specification. The deletion of another region has no effect on retinal progenitor cell proliferation but instead leads to a complete loss of bipolar neurons. This prototypical super-enhancer may serve as a model for dissecting the complex gene expression patterns for neurogenic transcription factors during development. Moreover, it provides a unique opportunity to alter expression of individual transcription factors in particular cell types at specific stages of development. This provides a deeper understanding of function that cannot be achieved with traditional knockout mouse approaches.

Analysis of long and short enhancers in melanoma cell states

Understanding how enhancers drive cell-type specificity and efficiently identifying them is essential for the development of innovative therapeutic strategies. In melanoma, the melanocytic (MEL) and the mesenchymal-like (MES) states present themselves with different responses to therapy, making the identification of specific enhancers highly relevant. Using massively parallel reporter assays (MPRAs) in a panel of patient-derived melanoma lines (MM lines), we set to identify and decipher melanoma enhancers by first focusing on regions with state-specific H3K27 acetylation close to differentially expressed genes. An in-depth evaluation of those regions was then pursued by investigating the activity of overlapping ATAC-seq peaks along with a full tiling of the acetylated regions with 190 bp sequences. Activity was observed in more than 60% of the selected regions, and we were able to precisely locate the active enhancers within ATAC-seq peaks. Comparison of sequence content with activity, using the deep learning model DeepMEL2, revealed that AP-1 alone is responsible for the MES enhancer activity. In contrast, SOX10 and MITF both influence MEL enhancer function with SOX10 being required to achieve high levels of activity. Overall, our MPRAs shed light on the relationship between long and short sequences in terms of their sequence content, enhancer activity, and specificity across melanoma cell states.

Alu insertion variants alter gene transcript levels

Alu are high copy number interspersed repeats that have accumulated near genes during primate and human evolution. They are a pervasive source of structural variation in modern humans. Impacts that Alu insertions may have on gene expression are not well understood, although some have been associated with expression quantitative trait loci (eQTLs). Here, we directly test regulatory effects of polymorphic Alu insertions in isolation of other variants on the same haplotype. To screen insertion variants for those with such effects, we used ectopic luciferase reporter assays and evaluated 110 Alu insertion variants, including more than 40 with a potential role in disease risk. We observed a continuum of effects with significant outliers that up- or down-regulate luciferase activity. Using a series of reporter constructs, which included genomic context surrounding the Alu, we can distinguish between instances in which the Alu disrupts another regulator and those in which the Alu introduces new regulatory sequence. We next focused on three polymorphic Alu loci associated with breast cancer that display significant effects in the reporter assay. We used CRISPR to modify the endogenous sequences, establishing cell lines varying in the Alu genotype. Our findings indicate that Alu genotype can alter expression of genes implicated in cancer risk, including PTHLH, RANBP9, and MYC These data show that commonly occurring polymorphic Alu elements can alter transcript levels and potentially contribute to disease risk.

Convergence of oncogenic cooperation at single-cell and single-gene levels drives leukemic transformation

Cancers develop from the accumulation of somatic mutations, yet it remains unclear how oncogenic lesions cooperate to drive cancer progression. Using a mouse model harboring NRasG12D and EZH2 mutations that recapitulates leukemic progression, we employ single-cell transcriptomic profiling to map cellular composition and gene expression alterations in healthy or diseased bone marrows during leukemogenesis. At cellular level, NRasG12D induces myeloid lineage-biased differentiation and EZH2-deficiency impairs myeloid cell maturation, whereas they cooperate to promote myeloid neoplasms with dysregulated transcriptional programs. At gene level, NRasG12D and EZH2-deficiency independently and synergistically deregulate gene expression. We integrate results from histopathology, leukemia repopulation, and leukemia-initiating cell assays to validate transcriptome-based cellular profiles. We use this resource to relate developmental hierarchies to leukemia phenotypes, evaluate oncogenic cooperation at single-cell and single-gene levels, and identify GEM as a regulator of leukemia-initiating cells. Our studies establish an integrative approach to deconvolute cancer evolution at single-cell resolution in vivo.