The genome-wide dynamics of the binding of Ldb1 complexes during erythroid differentiation

An intricate regulatory circuit between FLI1 and GATA1/GATA2/LDB1/ERG dictates erythroid vs. megakaryocytic differentiation

During hematopoiesis, megakaryocytic erythroid progenitors (MEPs) differentiate into megakaryocytic or erythroid lineages in response to specific transcriptional factors, yet the regulatory mechanism remains to be elucidated. Using the MEP‑like cell line HEL western blotting, RT‑qPCR, lentivirus‑mediated downregulation, flow cytometry as well as chromatin immunoprecipitation (ChIp) assay demonstrated that the E26 transformation‑specific (ETS) transcription factor friend leukemia integration factor 1 (Fli‑1) inhibits erythroid differentiation. The present study using these methods showed that while FLI1‑mediated downregulation of GATA binding protein 1 (GATA1) suppresses erythropoiesis, its direct transcriptional induction of GATA2 promotes megakaryocytic differentiation. GATA1 is also involved in megakaryocytic differentiation through regulation of GATA2. By contrast to FLI1, the ETS member erythroblast transformation‑specific‑related gene (ERG) negatively controls GATA2 and its overexpression through exogenous transfection blocks megakaryocytic differentiation. In addition, FLI1 regulates expression of LIM Domain Binding 1 (LDB1) during erythroid and megakaryocytic commitment, whereas shRNA‑mediated depletion of LDB1 downregulates FLI1 and GATA2 but increases GATA1 expression. In agreement, LDB1 ablation using shRNA lentivirus expression blocks megakaryocytic differentiation and modestly suppresses erythroid maturation. These results suggested that a certain threshold level of LDB1 expression enables FLI1 to block erythroid differentiation. Overall, FLI1 controlled the commitment of MEP to either erythroid or megakaryocytic lineage through an intricate regulation of GATA1/GATA2, LDB1 and ERG, exposing multiple targets for cell fate commitment and therapeutic intervention.

The ETO2 transcriptional cofactor maintains acute leukemia by driving a MYB/EP300-dependent stemness program

Transcriptional cofactors of the ETO family are recurrent fusion partners in acute leukemia. We characterized the ETO2 regulome by integrating transcriptomic and chromatin binding analyses in human erythroleukemia xenografts and controlled ETO2 depletion models. We demonstrate that beyond its well-established repressive activity, ETO2 directly activates transcription of MYB, among other genes. The ETO2-activated signature is associated with a poorer prognosis in erythroleukemia but also in other acute myeloid and lymphoid leukemia subtypes. Mechanistically, ETO2 colocalizes with EP300 and MYB at enhancers supporting the existence of an ETO2/MYB feedforward transcription activation loop (e.g., on MYB itself). Both small-molecule and PROTAC-mediated inhibition of EP300 acetyltransferases strongly reduced ETO2 protein, chromatin binding, and ETO2-activated transcripts. Taken together, our data show that ETO2 positively enforces a leukemia maintenance program that is mediated in part by the MYB transcription factor and that relies on acetyltransferase cofactors to stabilize ETO2 scaffolding activity.

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.

Efficient genome editing in erythroid cells unveils novel MYB target genes and regulatory functions

Targeted genome editing holds great promise in biology. However, efficient genome modification, including gene knock-in (KI), remains an unattained goal in multiple cell types and loci due to poor transfection efficiencies and low target genes expression, impeding the positive selection of recombined cells. Here, we describe a genome editing approach to achieve efficient gene targeting using hard to transfect erythroid cell lines. We demonstrate robust fluorescent protein KI efficiency in low expressed transcription factor (TF) genes (e.g., Myb or Zeb1). We further show the ability to target two independent loci in individual cells, exemplified by MYB-GFP and NuMA-Cherry double KI, allowing multicolor labeling of regulatory factors at physiological endogenous levels. Our KI tagging approach allowed us to perform genome-wide TF analysis at increased signal-to-noise ratios, and highlighted previously unidentified MYB target genes and pathways. Overall, we establish a versatile CRISPR-Cas9-based platform, offering attractive opportunities for the dissection of the erythroid differentiation process.

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.

The NFIA-ETO2 fusion blocks erythroid maturation and induces pure erythroid leukemia in cooperation with mutant TP53

The NFIA-ETO2 fusion is the product of a t(1;16)(p31;q24) chromosomal translocation, so far, exclusively found in pediatric patients with pure erythroid leukemia (PEL). To address the role for the pathogenesis of the disease, we facilitated the expression of the NFIA-ETO2 fusion in murine erythroblasts (EBs). We observed that NFIA-ETO2 significantly increased proliferation and impaired erythroid differentiation of murine erythroleukemia cells and of primary fetal liver-derived EBs. However, NFIA-ETO2-expressing EBs acquired neither aberrant in vitro clonogenic activity nor disease-inducing potential upon transplantation into irradiated syngenic mice. In contrast, in the presence of 1 of the most prevalent erythroleukemia-associated mutations, TP53R248Q, expression of NFIA-ETO2 resulted in aberrant clonogenic activity and induced a fully penetrant transplantable PEL-like disease in mice. Molecular studies support that NFIA-ETO2 interferes with erythroid differentiation by preferentially binding and repressing erythroid genes that contain NFI binding sites and/or are decorated by ETO2, resulting in a activity shift from GATA- to ETS-motif-containing target genes. In contrast, TP53R248Q does not affect erythroid differentiation but provides self-renewal and survival potential, mostly via downregulation of known TP53 targets. Collectively, our work indicates that NFIA-ETO2 initiates PEL by suppressing gene expression programs of terminal erythroid differentiation and cooperates with TP53 mutation to induce erythroleukemia.

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.

Heterochromatin rewiring and domain disruption-mediated chromatin compaction during erythropoiesis

Mammalian erythropoiesis involves progressive chromatin compaction and subsequent enucleation in terminal differentiation, but the mechanisms underlying the three-dimensional chromatin reorganization remain obscure. Here, we systematically analyze the higher-order chromatin in purified populations of primary human erythroblasts. Our results reveal that heterochromatin regions undergo substantial compression, with H3K9me3 markers relocalizing to the nuclear periphery and forming a significant number of long-range interactions, and that ~58% of the topologically associating domain (TAD) boundaries are disrupted, while certain TADs enriched for markers of the active transcription state and erythroid master regulators, GATA1 and KLF1, are selectively maintained during terminal erythropoiesis. Finally, we demonstrate that GATA1 is involved in safeguarding selected essential chromatin domains during terminal erythropoiesis. Our study therefore delineates the molecular characteristics of a development-driven chromatin compaction process, which reveals transcription competence as a key indicator of the selected domain maintenance to ensure appropriate gene expression during the extreme compaction of chromatin.

Multi-modular structure of the gene regulatory network for specification and commitment of murine T cells

T cells develop from multipotent progenitors by a gradual process dependent on intrathymic Notch signaling and coupled with extensive proliferation. The stages leading them to T-cell lineage commitment are well characterized by single-cell and bulk RNA analyses of sorted populations and by direct measurements of precursor-product relationships. This process depends not only on Notch signaling but also on multiple transcription factors, some associated with stemness and multipotency, some with alternative lineages, and others associated with T-cell fate. These factors interact in opposing or semi-independent T cell gene regulatory network (GRN) subcircuits that are increasingly well defined. A newly comprehensive picture of this network has emerged. Importantly, because key factors in the GRN can bind to markedly different genomic sites at one stage than they do at other stages, the genes they significantly regulate are also stage-specific. Global transcriptome analyses of perturbations have revealed an underlying modular structure to the T-cell commitment GRN, separating decisions to lose "stem-ness" from decisions to block alternative fates. Finally, the updated network sheds light on the intimate relationship between the T-cell program, which depends on the thymus, and the innate lymphoid cell (ILC) program, which does not.

Elucidation of the Role of FAM210B in Mitochondrial Metabolism and Erythropoiesis

Mitochondria play essential and specific roles during erythroid differentiation. Recently, FAM210B, encoding a mitochondrial inner membrane protein, has been identified as a novel target of GATA-1, as well as an erythropoietin-inducible gene. While FAM210B protein is involved in regulate mitochondrial metabolism and heme biosynthesis, its detailed function remains unknown. Here, we generated both knockout and knockdown of endogenous FAM210B in human induced pluripotent stem-derived erythroid progenitor (HiDEP) cells using CRISPR/Cas9 methodology. Intriguingly, erythroid differentiation was more pronounced in the FAM210B-depleted cells, and this resulted in increased frequency of orthochromatic erythroblasts and decreased frequencies of basophilic/polychromatic erythroblasts. Comprehensive metabolite analysis and functional analysis indicated that oxygen consumption rates and the NAD (NAD⁺)/NADH ratio were significantly decreased, while lactate production was significantly increased in FAM210B deletion HiDEP cells, indicating involvement of FAM210B in mitochondrial energy metabolism in erythroblasts. Finally, we purified FAM210B-interacting protein from K562 cells that stably expressed His/biotin-tagged FAM210B. Mass spectrometry analysis of the His/biotin-purified material indicated interactions with multiple subunits of mitochondrial ATP synthases, such as subunit alpha (ATP5A) and beta (ATP5B). Our results suggested that FAM210B contributes prominently to erythroid differentiation by regulating mitochondrial energy metabolism. Our results provide insights into the pathophysiology of dysregulated hematopoiesis.

Reinforced erythroid differentiation inhibits leukemogenic potential of t(8;21) leukemia

Oncoprotein AML1-ETO (AE) derived from t(8;21)(q22;q22) translocation is typically present in a portion of French-American-British-M2 subtype of acute myeloid leukemia (AML). Although these patients have relatively favorable prognoses, substantial numbers of them would relapse after conventional therapy. Here, we explored whether reinforcing the endogenous differentiation potential of t(8;21) AML cells would diminish the associated malignancy. In doing so, we noticed an expansion of immature erythroid blasts featured in both AML1-ETO9a (AE9a) and AE plus c-KIT (N822K) (AK) murine leukemic models. Interestingly, in the AE9a murine model, a spontaneous step-wise erythroid differentiation path, as characterized by the differential expression of CD43/c-Kit and the upregulation of several key erythroid transcription factors (TFs), accompanied the decline or loss of leukemia-initiating potential. Notably, overexpression of one of the key erythroid TFs, Ldb1, potently disrupted the repopulation of AE9a leukemic cells in vivo, suggesting a new promising intervention strategy of t(8;21) AML through enforcing their erythroid differentiation.

MTG8 interacts with LHX6 to specify cortical interneuron subtype identity

Cortical interneurons originating in the embryonic medial ganglionic eminence (MGE) diverge into a range of different subtypes found in the adult mouse cerebral cortex. The mechanisms underlying this divergence and the timing when subtype identity is set up remain unclear. We identify the highly conserved transcriptional co-factor MTG8 as being pivotal in the development of a large subset of MGE cortical interneurons that co-expresses Somatostatin (SST) and Neuropeptide Y (NPY). MTG8 interacts with the pan-MGE transcription factor LHX6 and together the two factors are sufficient to promote expression of critical cortical interneuron subtype identity genes. The SST-NPY cortical interneuron fate is initiated early, well before interneurons migrate into the cortex, demonstrating an early onset specification program. Our findings suggest that transcriptional co-factors and modifiers of generic lineage specification programs may hold the key to the emergence of cortical interneuron heterogeneity from the embryonic telencephalic germinal zones.

There is a large diversity of inhibitory interneurons in the mammalian cerebral cortex. How this emerges during embryogenesis remains unclear. Here, the authors identify MTG8 as a co-factor of LHX6 and a new regulator of cortical interneuron development.

Enhancer looping protein LDB1 regulates hepatocyte gene expression by cooperating with liver transcription factors

Enhancers establish proximity with distant target genes to regulate temporospatial gene expression and specify cell identity. Lim domain binding protein 1 (LDB1) is a conserved and widely expressed protein that functions as an enhancer looping factor. Previous studies in erythroid cells and neuronal cells showed that LDB1 forms protein complexes with different transcription factors to regulate cell-specific gene expression. Here, we show that LDB1 regulates expression of liver genes by occupying enhancer elements and cooperating with hepatic transcription factors HNF4A, FOXA1, TCF7 and GATA4. Using the glucose transporter SLC2A2 gene, encoding GLUT2, as an example, we find that LDB1 regulates gene expression by mediating enhancer-promoter interactions. In vivo, we find that LDB1 deficiency in primary mouse hepatocytes dysregulates metabolic gene expression and changes the enhancer landscape. Conditional deletion of LDB1 in adult mouse liver induces glucose intolerance. However, Ldb1 knockout hepatocytes show improved liver pathology under high-fat diet conditions associated with increased expression of genes related to liver fatty acid metabolic processes. Thus, LDB1 is linked to liver metabolic functions under normal and obesogenic conditions.

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.

Current insights into the role of Fli-1 in hematopoiesis and malignant transformation

Fli-1, a member of the ETS family of transcription factors, was discovered in 1991 through retroviral insertional mutagenesis as a driver of mouse erythroleukemias. In the past 30 years, nearly 2000 papers have defined its biology and impact on normal development and cancer. In the hematopoietic system, Fli-1 controls self-renewal of stem cells and their differentiation into diverse mature blood cells. Fli-1 also controls endothelial survival and vasculogenesis, and high and low levels of Fli-1 are implicated in the auto-immune diseases systemic lupus erythematosus and systemic sclerosis, respectively. In addition, aberrant Fli-1 expression is observed in, and is essential for, the growth of multiple hematological malignancies and solid cancers. Here, we review the historical context and latest research on Fli-1, focusing on its role in hematopoiesis, immune response, and malignant transformation. The importance of identifying Fli-1 modulators (both agonists and antagonists) and their potential clinical applications is discussed.

Application of Weighted Gene Coexpression Network Analysis to Identify Key Modules and Hub Genes in Systemic Juvenile Idiopathic Arthritis

Systemic juvenile idiopathic arthritis (sJIA) is a severe autoinflammatory disorder with a still not clearly defined molecular mechanism. To better understand the disease, we used scattered datasets from public domains and performed a weighted gene coexpression network analysis (WGCNA) to identify key modules and hub genes underlying sJIA pathogenesis. Two gene expression datasets, GSE7753 and GSE13501, were used to construct the WGCNA. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were applied to the genes and hub genes in the sJIA modules. Cytoscape was used to screen and visualize the hub genes. We further compared the hub genes with the genome-wide association study (GWAS) genes and used a consensus WGCNA to verify that our conclusions were conservative and reproducible across multiple independent datasets. A total of 5,414 genes were obtained for WGCNA, from which highly correlated genes were divided into 17 modules. The red module demonstrated the highest correlation with the sJIA module (r = 0.8, p = 3e⁻²⁹), whereas the green-yellow module was found to be closely related to the non-sJIA module (r = 0.62, p = 1e⁻¹⁴). Functional enrichment analysis demonstrated that the red module was mostly enriched in the activation of immune responses, infection, nucleosomes, and erythrocytes, and the green-yellow module was mostly enriched in immune responses and inflammation. Additionally, the hub genes in the red module were highly enriched in erythrocyte differentiation, including ALAS2, AHSP, TRIM10, TRIM58, and KLF1. The hub genes from the green-yellow module were mainly associated with immune responses, as exemplified by the genes KLRB1, KLRF1, CD160, and KIRs. We identified sJIA-related modules and several hub genes that might be associated with the development of sJIA. Particularly, the modules may help understand the mechanisms of sJIA, and the hub genes may become biomarkers and therapeutic targets of sJIA in the future.

Mechanisms of Binding Specificity among bHLH Transcription Factors

The transcriptome of every cell is orchestrated by the complex network of interaction between transcription factors (TFs) and their binding sites on DNA. Disruption of this network can result in many forms of organism malfunction but also can be the substrate of positive natural selection. However, understanding the specific determinants of each of these individual TF-DNA interactions is a challenging task as it requires integrating the multiple possible mechanisms by which a given TF ends up interacting with a specific genomic region. These mechanisms include DNA motif preferences, which can be determined by nucleotide sequence but also by DNA’s shape; post-translational modifications of the TF, such as phosphorylation; and dimerization partners and co-factors, which can mediate multiple forms of direct or indirect cooperative binding. Binding can also be affected by epigenetic modifications of putative target regions, including DNA methylation and nucleosome occupancy. In this review, we describe how all these mechanisms have a role and crosstalk in one specific family of TFs, the basic helix-loop-helix (bHLH), with a very conserved DNA binding domain and a similar DNA preferred motif, the E-box. Here, we compile and discuss a rich catalog of strategies used by bHLH to acquire TF-specific genome-wide landscapes of binding sites.

A Positive Regulatory Feedback Loop between EKLF/KLF1 and TAL1/SCL Sustaining the Erythropoiesis

The erythroid Krüppel-like factor EKLF/KLF1 is a hematopoietic transcription factor binding to the CACCC DNA motif and participating in the regulation of erythroid differentiation. With combined use of microarray-based gene expression profiling and the promoter-based ChIP-chip assay of E14.5 fetal liver cells from wild type (WT) and EKLF-knockout (Eklf^−/−) mouse embryos, we identified the pathways and direct target genes activated or repressed by EKLF. This genome-wide study together with the molecular/cellular analysis of the mouse erythroleukemic cells (MEL) indicate that among the downstream direct target genes of EKLF is Tal1/Scl. Tal1/Scl encodes another DNA-binding hematopoietic transcription factor TAL1/SCL, known to be an Eklf activator and essential for definitive erythroid differentiation. Further identification of the authentic Tal gene promoter in combination with the in vivo genomic footprinting approach and DNA reporter assay demonstrate that EKLF activates the Tal gene through binding to a specific CACCC motif located in its promoter. These data establish the existence of a previously unknow positive regulatory feedback loop between two DNA-binding hematopoietic transcription factors, which sustains mammalian erythropoiesis.

A gain-of-function single nucleotide variant creates a new promoter which acts as an orientation-dependent enhancer-blocker

Many single nucleotide variants (SNVs) associated with human traits and genetic diseases are thought to alter the activity of existing regulatory elements. Some SNVs may also create entirely new regulatory elements which change gene expression, but the mechanism by which they do so is largely unknown. Here we show that a single base change in an otherwise unremarkable region of the human α-globin cluster creates an entirely new promoter and an associated unidirectional transcript. This SNV downregulates α-globin expression causing α-thalassaemia. Of note, the new promoter lying between the α-globin genes and their associated super-enhancer disrupts their interaction in an orientation-dependent manner. Together these observations show how both the order and orientation of the fundamental elements of the genome determine patterns of gene expression and support the concept that active genes may act to disrupt enhancer-promoter interactions in mammals as in Drosophila. Finally, these findings should prompt others to fully evaluate SNVs lying outside of known regulatory elements as causing changes in gene expression by creating new regulatory elements.

The role of promoters as potential insulator elements has been largely unexplored in mammals. Here the authors show that a single nucleotide variant in the α-globin locus forms a new promoter and acts as an orientation-dependent enhancer-blocking insulator element.

ATF4 Regulates MYB to Increase γ-Globin in Response to Loss of β-Globin

β-Hemoglobinopathies can trigger rapid production of red blood cells in a process known as stress erythropoiesis. Cellular stress prompts differentiating erythroid precursors to express high levels of fetal γ-globin. However, the mechanisms underlying γ-globin production during cellular stress are still poorly defined. Here, we use CRISPR-Cas genome editing to model the stress caused by reduced levels of adult β-globin. We find that decreased β-globin is sufficient to induce robust re-expression of γ-globin, and RNA sequencing (RNA-seq) of differentiating isogenic erythroid precursors implicates ATF4 as a causal regulator of this response. ATF4 binds within the HBS1L-MYB intergenic enhancer and regulates expression of MYB, a known γ-globin regulator. Overall, the reduction of ATF4 upon β-globin knockout decreases the levels of MYB and BCL11A. Identification of ATF4 as a key regulator of globin compensation adds mechanistic insight to the poorly understood phenomenon of stress-induced globin compensation and could inform strategies to treat hemoglobinopathies.

Reduction of RUNX1 transcription factor activity by a CBFA2T3-mimicking peptide: application to B cell precursor acute lymphoblastic leukemia

Background

B Cell Precursor Acute Lymphoblastic Leukemia (BCP-ALL) is the most common pediatric cancer. Identifying key players involved in proliferation of BCP-ALL cells is crucial to propose new therapeutic targets. Runt Related Transcription Factor 1 (RUNX1) and Core-Binding Factor Runt Domain Alpha Subunit 2 Translocated To 3 (CBFA2T3, ETO2, MTG16) are master regulators of hematopoiesis and are implicated in leukemia.

Methods

We worked with BCP-ALL mononuclear bone marrow patients’ cells and BCP-ALL cell lines, and performed Chromatin Immunoprecipitations followed by Sequencing (ChIP-Seq), co-immunoprecipitations (co-IP), proximity ligation assays (PLA), luciferase reporter assays and mouse xenograft models.

Results

We demonstrated that CBFA2T3 transcript levels correlate with RUNX1 expression in the pediatric t(12;21) ETV6-RUNX1 BCP-ALL. By ChIP-Seq in BCP-ALL patients’ cells and cell lines, we found that RUNX1 is recruited on its promoter and on an enhancer of CBFA2T3 located − 2 kb upstream CBFA2T3 promoter and that, subsequently, the transcription factor RUNX1 drives both RUNX1 and CBFA2T3 expression. We demonstrated that, mechanistically, RUNX1 and CBFA2T3 can be part of the same complex allowing CBFA2T3 to strongly potentiate the activity of the transcription factor RUNX1. Finally, we characterized a CBFA2T3-mimicking peptide that inhibits the interaction between RUNX1 and CBFA2T3, abrogating the activity of this transcription complex and reducing BCP-ALL lymphoblast proliferation.

Conclusions

Altogether, our findings reveal a novel and important activation loop between the transcription regulator CBFA2T3 and the transcription factor RUNX1 that promotes BCP-ALL proliferation, supporting the development of an innovative therapeutic approach based on the NHR2 subdomain of CBFA2T3 protein.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13045-021-01051-z.

Interplay between FLI-1 and the LDB1 complex in murine erythroleukemia cells and during megakaryopoiesis

Transcription factors are key players in a broad range of cellular processes such as cell-fate decision. Understanding how they act to control these processes is of critical importance for therapy purposes. FLI-1 controls several hematopoietic lineage differentiation including megakaryopoiesis and erythropoiesis. Its aberrant expression is often observed in cancer and is associated with poor prognosis. We showed that FLI-1 interacts with the LDB1 complex, which also plays critical roles in erythropoiesis and megakaryopoiesis. In this study, we aimed to unravel how FLI-1 and the LDB1 complex act together in murine erythroleukemia cells and in megakaryocyte. Combining omics techniques, we show that FLI-1 enables the recruitment of the LDB1 complex to regulatory sequences of megakaryocytic genes and to enhancers. We show as well for the first time that FLI-1 is able to modulate the 3D chromatin organization by promoting chromatin looping between enhancers and promoters most likely through the LDB1 complex.

•

FLI-1 is important for the recruitment of the LDB1 complex

•

FLI-1 is important for chromatin looping

•

FLI-1 and the LDB1 complex co-regulate megakaryopoiesis

TATA and paused promoters active in differentiated tissues have distinct expression characteristics

Core promoter types differ in the extent to which RNA polymerase II (Pol II) pauses after initiation, but how this affects their tissue-specific gene expression characteristics is not well understood. While promoters with Pol II pausing elements are active throughout development, TATA promoters are highly active in differentiated tissues. We therefore used a genomics approach on late-stage Drosophila embryos to analyze the properties of promoter types. Using tissue-specific Pol II ChIP-seq, we found that paused promoters have high levels of paused Pol II throughout the embryo, even in tissues where the gene is not expressed, while TATA promoters only show Pol II occupancy when the gene is active. The promoter types are associated with different chromatin accessibility in ATAC-seq data and have different expression characteristics in single-cell RNA-seq data. The two promoter types may therefore be optimized for different properties: paused promoters show more consistent expression when active, while TATA promoters have lower background expression when inactive. We propose that tissue-specific genes have evolved to use two different strategies for their differential expression across tissues.

When basic science reaches into rational therapeutic design: from historical to novel leads for the treatment of β-globinopathies

PURPOSE OF REVIEW

β-hemoglobinopathies, such as β-Thalassemias (β-Thal) and sickle cell disease (SCD) are among the most common inherited genetic disorders in humans worldwide. These disorders are characterized by a quantitative (β-Thal) or qualitative (SCD) defects in adult hemoglobin production, leading to anemia, ineffective erythropoiesis and severe secondary complications. Reactivation of the fetal globin genes (γ-globin), making-up fetal hemoglobin (HbF), which are normally silenced in adults, represents a major strategy to ameliorate anemia and disease severity.

RECENT FINDINGS

Following the identification of the first 'switching factors' for the reactivation of fetal globin gene expression more than 10 years ago, a multitude of novel leads have recently been uncovered.

SUMMARY

Recent findings provided invaluable functional insights into the genetic and molecular networks controlling globin genes expression, revealing that complex repression systems evolved in erythroid cells to maintain HbF silencing in adults. This review summarizes these unique and exciting discoveries of the regulatory factors controlling the globin switch. New insights and novel leads for therapeutic strategies based on the pharmacological induction of HbF are discussed. This represents a major breakthrough for rational drug design in the treatment of β-Thal and SCD.

The transcription factor TAL1 and miR-17-92 create a regulatory loop in hematopoiesis

A network of gene regulatory factors such as transcription factors and microRNAs establish and maintain gene expression patterns during hematopoiesis. In this network, transcription factors regulate each other and are involved in regulatory loops with microRNAs. The microRNA cluster miR-17-92 is located within the MIR17HG gene and encodes six mature microRNAs. It is important for hematopoietic differentiation and plays a central role in malignant disease. However, the transcription factors downstream of miR-17-92 are largely elusive and the transcriptional regulation of miR-17-92 is not fully understood. Here we show that miR-17-92 forms a regulatory loop with the transcription factor TAL1. The miR-17-92 cluster inhibits expression of TAL1 and indirectly leads to decreased stability of the TAL1 transcriptional complex. We found that TAL1 and its heterodimerization partner E47 regulate miR-17-92 transcriptionally. Furthermore, miR-17-92 negatively influences erythroid differentiation, a process that depends on gene activation by the TAL1 complex. Our data give example of how transcription factor activity is fine-tuned during normal hematopoiesis. We postulate that disturbance of the regulatory loop between TAL1 and the miR-17-92 cluster could be an important step in cancer development and progression.

Comparison of the PU.1 transcriptional regulome and interactome in human and mouse inflammatory dendritic cells

Dendritic cells (DCs) are key immune modulators and are able to mount immune responses or tolerance. DC differentiation and activation imply a plethora of molecular and cellular responses, including transcriptional changes. PU.1 is a highly expressed transcription factor in DCs and coordinates relevant aspects of DC biology. Due to their role as immune regulators, DCs pose as a promising immunotherapy tool. However, some of their functional features, such as survival, activation, or migration, are compromised due to the limitations to simulate in vitro the physiologic DC differentiation process. A better knowledge of transcriptional programs would allow the identification of potential targets for manipulation with the aim of obtaining "qualified" DCs for immunotherapy purposes. Most of the current knowledge regarding DC biology derives from studies using mouse models, which not always find a parallel in human. In the present study, we dissect the PU.1 transcriptional regulome and interactome in mouse and human DCs, in the steady state or LPS activated. The PU.1 transcriptional regulome was identified by performing PU.1 chromatin immunoprecipitation followed by high-throughput sequencing and pairing these data with RNAsequencing data. The PU.1 interactome was identified by performing PU.1 immunoprecipitation followed by mass spectrometry analysis. Our results portray PU.1 as a pivotal factor that plays an important role in the regulation of genes required for proper DC activation and function, and assures the repression of nonlineage genes. The interspecies differences between human and mouse DCs are surprisingly substantial, highlighting the need to study the biology of human DCs.

Embryonic erythropoiesis and hemoglobin switching require transcriptional repressor ETO2 to modulate chromatin organization

The underlying mechanism of transcriptional co-repressor ETO2 during early erythropoiesis and hemoglobin switching is unclear. We find that absence of ETO2 in mice interferes with down-regulation of PU.1 and GATA2 in the fetal liver, impeding a key step required for commitment to erythroid maturation. In human β-globin transgenic Eto2 null mice and in human CD34+ erythroid progenitor cells with reduced ETO2, loss of ETO2 results in ineffective silencing of embryonic/fetal globin gene expression, impeding hemoglobin switching during erythroid differentiation. ETO2 occupancy genome-wide occurs virtually exclusively at LDB1-complex binding sites in enhancers and ETO2 loss leads to increased enhancer activity and expression of target genes. ETO2 recruits the NuRD nucleosome remodeling and deacetylation complex to regulate histone acetylation and nucleosome occupancy in the β-globin locus control region and γ-globin gene. Loss of ETO2 elevates LDB1, MED1 and Pol II in the locus and facilitates fetal γ-globin/LCR looping and γ-globin transcription. Absence of the ETO2 hydrophobic heptad repeat region impairs ETO2-NuRD interaction and function in antagonizing γ-globin/LCR looping. Our results reveal a pivotal role for ETO2 in erythropoiesis and globin gene switching through its repressive role in the LDB1 complex, affecting the transcription factor and epigenetic environment and ultimately restructuring chromatin organization.

Human erythroleukemia genetics and transcriptomes identify master transcription factors as functional disease drivers

Acute erythroleukemia (AEL or acute myeloid leukemia [AML]-M6) is a rare but aggressive hematologic malignancy. Previous studies showed that AEL leukemic cells often carry complex karyotypes and mutations in known AML-associated oncogenes. To better define the underlying molecular mechanisms driving the erythroid phenotype, we studied a series of 33 AEL samples representing 3 genetic AEL subgroups including TP53-mutated, epigenetic regulator-mutated (eg, DNMT3A, TET2, or IDH2), and undefined cases with low mutational burden. We established an erythroid vs myeloid transcriptome-based space in which, independently of the molecular subgroup, the majority of the AEL samples exhibited a unique mapping different from both non-M6 AML and myelodysplastic syndrome samples. Notably, >25% of AEL patients, including in the genetically undefined subgroup, showed aberrant expression of key transcriptional regulators, including SKI, ERG, and ETO2. Ectopic expression of these factors in murine erythroid progenitors blocked in vitro erythroid differentiation and led to immortalization associated with decreased chromatin accessibility at GATA1-binding sites and functional interference with GATA1 activity. In vivo models showed development of lethal erythroid, mixed erythroid/myeloid, or other malignancies depending on the cell population in which AEL-associated alterations were expressed. Collectively, our data indicate that AEL is a molecularly heterogeneous disease with an erythroid identity that results in part from the aberrant activity of key erythroid transcription factors in hematopoietic stem or progenitor cells.

LDB1 Enforces Stability on Direct and Indirect Oncoprotein Partners in Leukemia

The LMO2/LDB1 macromolecular complex is critical in hematopoietic stem and progenitor cell specification and in the development of acute leukemia. This complex is comprised of core subunits of LMO2 and LDB1 as well as single-stranded DNA-binding protein (SSBP) cofactors and DNA-binding basic helix-loop-helix (bHLH) and GATA transcription factors. We analyzed the steady-state abundance and kinetic stability of LMO2 and its partners via Halo protein tagging in conjunction with variant proteins deficient in binding their respective direct protein partners. We discovered a hierarchy of protein stabilities (with half-lives in descending order) as follows: LDB1 > SSBP > LMO2 > TAL1. Importantly, LDB1 is a remarkably stable protein that confers enhanced stability upon direct and indirect partners, thereby nucleating the formation of the multisubunit protein complex. The data imply that free subunits are more rapidly degraded than those incorporated within the LMO2/LDB1 complex. Our studies provided significant insights into LMO2/LDB1 macromolecular protein complex assembly and stability, which has implications for understanding its role in blood cell formation and for therapeutically targeting this complex in human leukemias.

LDB1 and the SWI/SNF complex participate in both transcriptional activation and repression by Caenorhabditis elegans BLIMP1/PRDM1

Transcription factors of the BLIMP1/PRDM1 family are important regulators of development. BLIMP1/PRDM1 can both activate and repress gene expression, however, the mechanism of activation is not well understood. Therefore, we looked for factors involved in gene activation by C. elegans BLMP-1, the ortholog of BLIMP1/PRDM1. BLMP-1 activates the expression of bed-3, a gene involved in vulval development. By screening nuclear proteins that function in vulval development, we identified two proteins (LDB-1 and HAM-3) required for BLMP-1 dependent bed-3 expression. LDB-1 is the sole C. elegans member of the LIM Binding Protein (LDB) family, whereas HAM-3 is an accessory subunit of the SWI/SNF complex (ortholog of human SMARCD3/BAF60C). A core SWI/SNF subunit SWSN-1 (ortholog of human SMARCC1/BAF155) is also involved. We found that LDB-1 and HAM-3 bind to BLMP-1, suggesting that BLMP-1 recruits LDB-1 and the SWI/SNF complex to activate bed-3 expression. Interestingly, LDB-1 and HAM-3 are involved in both transcriptional activation and repression. In particular, BLMP-1, LDB-1 and HAM-3 co-regulate a set of hypodermal genes including bed-3 (activated), col-124 (activated) and lin-29 (repressed). On the other hand, LDB-1 and HAM-3 are not required for activation or repression of some genes regulated by BLMP-1 (e.g. T09D3.8, nas-10). We also found that human LDB1, SMARCD3/BAF60C and SMARCC1/BAF155 all physically interact with human BLIMP1/PRDM1 in vitro and are closely associated with BLIMP1/PRDM1 in vivo. Taken together, these results identify LDB1 and SWI/SNF as likely conserved cofactors of BLIMP1/PRDM1, which participate in activation and repression of a subset of BLIMP1/PRDM1-regulated genes.

A non-linear reverse-engineering method for inferring genetic regulatory networks

Hematopoiesis is a highly complex developmental process that produces various types of blood cells. This process is regulated by different genetic networks that control the proliferation, differentiation, and maturation of hematopoietic stem cells (HSCs). Although substantial progress has been made for understanding hematopoiesis, the detailed regulatory mechanisms for the fate determination of HSCs are still unraveled. In this study, we propose a novel approach to infer the detailed regulatory mechanisms. This work is designed to develop a mathematical framework that is able to realize nonlinear gene expression dynamics accurately. In particular, we intended to investigate the effect of possible protein heterodimers and/or synergistic effect in genetic regulation. This approach includes the Extended Forward Search Algorithm to infer network structure (top-down approach) and a non-linear mathematical model to infer dynamical property (bottom-up approach). Based on the published experimental data, we study two regulatory networks of 11 genes for regulating the erythrocyte differentiation pathway and the neutrophil differentiation pathway. The proposed algorithm is first applied to predict the network topologies among 11 genes and 55 non-linear terms which may be for heterodimers and/or synergistic effect. Then, the unknown model parameters are estimated by fitting simulations to the expression data of two different differentiation pathways. In addition, the edge deletion test is conducted to remove possible insignificant regulations from the inferred networks. Furthermore, the robustness property of the mathematical model is employed as an additional criterion to choose better network reconstruction results. Our simulation results successfully realized experimental data for two different differentiation pathways, which suggests that the proposed approach is an effective method to infer the topological structure and dynamic property of genetic regulations.

Regulation of GATA1 levels in erythropoiesis

GATA1 is considered as the "master" transcription factor in erythropoiesis. It regulates at the transcriptional level all aspects of erythroid maturation and function, as revealed by gene knockout studies in mice and by genome-wide occupancies in erythroid cells. The GATA1 protein contains two zinc finger domains and an N-terminal transactivation domain. GATA1 translation results in the production of the full-length protein and of a shorter variant (GATA1s) lacking the N-terminal transactivation domain, which is functionally deficient in supporting erythropoiesis. GATA1 protein abundance is highly regulated in erythroid cells at different levels, including transcription, mRNA translation, posttranslational modifications, and protein degradation, in a differentiation-stage-specific manner. Maintaining high GATA1 protein levels is essential in the early stages of erythroid maturation, whereas downregulating GATA1 protein levels is a necessary step in terminal erythroid differentiation. The importance of maintaining proper GATA1 protein homeostasis in erythropoiesis is demonstrated by the fact that both GATA1 loss and its overexpression result in lethal anemia. Importantly, alterations in any of those GATA1 regulatory checkpoints have been recognized as an important cause of hematological disorders such as dyserythropoiesis (with or without thrombocytopenia), β-thalassemia, Diamond-Blackfan anemia, myelodysplasia, or leukemia. In this review, we provide an overview of the multilevel regulation of GATA1 protein homeostasis in erythropoiesis and of its deregulation in hematological disease.

The dynamic emergence of GATA1 complexes identified in in vitro embryonic stem cell differentiation and in vivo mouse fetal liver

GATA1 is an essential transcriptional regulator of myeloid hematopoietic differentiation towards red blood cells. During erythroid differentiation, GATA1 forms different complexes with other transcription factors such as LDB1, TAL1, E2A and LMO2 ("the LDB1 complex") or with FOG1. The functions of GATA1 complexes have been studied extensively in definitive erythroid differentiation; however, the temporal and spatial formation of these complexes during erythroid development is unknown. We applied proximity ligation assay (PLA) to detect, localize and quantify individual interactions during embryonic stem cell differentiation and in mouse fetal liver (FL) tissue. We show that GATA1/LDB1 interactions appear before the proerythroblast stage and increase in a subset of the CD71⁺/TER119^- cells to activate the terminal erythroid differentiation program in 12.5 day FL. Using Ldb1 and Gata1 knockdown FL cells, we studied the functional contribution of the GATA1/LDB1 complex during differentiation. This shows that the active LDB1 complex appears quite late at the proerythroblast stage of differentiation and confirms the power of PLA in studying the dynamic interaction of proteins in cell differentiation at the single cell level. We provide dynamic insight into the temporal and spatial formation of the GATA1 and LDB1 transcription factor complexes during hematopoietic development and differentiation.

GATA factor transcriptional activity: Insights from genome-wide binding profiles

The members of the GATA family of transcription factors have homologous zinc fingers and bind to similar sequence motifs. Recent advances in genome-wide technologies and the integration of bioinformatics data have led to a better understanding of how GATA factors regulate gene expression; GATA-factor-induced transcriptional and epigenetic changes have now been analyzed at unprecedented levels of detail. Here, we review the results of genome-wide studies of GATA factor occupancy in human and murine cell lines and primary cells (as determined by chromatin immunoprecipitation sequencing), and then discuss the molecular mechanisms underlying the mediation of transcriptional and epigenetic regulation by GATA factors.

Mouse Tryptase Gene Expression is Coordinately Regulated by GATA1 and GATA2 in Bone Marrow-Derived Mast Cells

Mast cell tryptases have crucial roles in allergic and inflammatory diseases. The mouse tryptase genes represent a cluster of loci on chromosome 16p3.3. While their functional studies have been extensively performed, transcriptional regulation of tryptase genes is poorly understood. In this study, we examined the molecular basis of the tryptase gene expression in bone marrow-derived mast cells (BMMCs) of C57BL/6 mice and in MEDMC-BRC6 mast cells. The expression of the Tpsb2 and Tpsg1 genes, which reside at the 3′-end of the tryptase locus, is significantly decreased by the reduction of the GATA transcription factors GATA1 or GATA2. Chromatin immunoprecipitation assays have shown that the GATA factors bind at multiple regions within the locus, including 1.0 and 72.8 kb upstream of the Tpsb2 gene, and that GATA1 and GATA2 facilitate each other’s DNA binding activity to these regions. Deletion of the −72.8 kb region by genome editing significantly reduced the Tpsb2 and Tpsg1 mRNA levels in MEDMC-BRC6 cells. Furthermore, binding of CTCF and the cohesin subunit Rad21 was found upstream of the −72.8 kb region and was significantly reduced in the absence of GATA1. These results suggest that mouse tryptase gene expression is coordinately regulated by GATA1 and GATA2 in BMMCs.

Mbd2-CP2c loop drives adult-type globin gene expression and definitive erythropoiesis

During hematopoiesis, red blood cells originate from the hematopoietic stem cell reservoir. Although the regulation of erythropoiesis and globin expression has been intensively investigated, the underlining mechanisms are not fully understood, including the interplay between transcription factors and epigenetic factors. Here, we uncover that the Mbd2-free NuRD chromatin remodeling complex potentiates erythroid differentiation of proerythroblasts via managing functions of the CP2c complexes. We found that both Mbd2 and Mbd3 expression is downregulated during differentiation of MEL cells in vitro and in normal erythropoiesis in mouse bone marrow, and Mbd2 downregulation is crucial for erythropoiesis. In uninduced MEL cells, the Mbd2-NuRD complex is recruited to the promoter via Gata1/Fog1, and, via direct binding through p66α, it acts as a transcriptional inhibitor of the CP2c complexes, preventing their DNA binding and promoting degradation of the CP2c family proteins to suppress globin gene expression. Conversely, during erythropoiesis in vitro and in vivo, the Mbd2-free NuRD does not dissociate from the chromatin and acts as a transcriptional coactivator aiding the recruitment of the CP2c complexes to chromatin, and thereby leading to the induction of the active hemoglobin synthesis and erythroid differentiation. Our study highlights the regulation of erythroid differentiation by the Mbd2-CP2c loop.

DNA·RNA triple helix formation can function as a cis-acting regulatory mechanism at the human β-globin locus

We have identified regulatory mechanisms in which an RNA transcript forms a DNA duplex·RNA triple helix with a gene or one of its regulatory elements, suggesting potential auto-regulatory mechanisms in vivo. We describe an interaction at the human β-globin locus, in which an RNA segment embedded in the second intron of the β-globin gene forms a DNA·RNA triplex with the HS2 sequence within the β-globin locus control region, a major regulator of globin expression. We show in human K562 cells that the triplex is stable in vivo. Its formation causes displacement from HS2 of major transcription factors and RNA Polymerase II, and consequently in loss of factors and polymerase that bind to the human ε- and γ-globin promoters, which are activated by HS2 in K562 cells. This results in reduced expression of these genes. These effects are observed when a small length of triplex-forming RNA is introduced into cells, or when a full-length intron-containing human β-globin transcript is expressed. Related results are obtained in human umbilical cord blood-derived erythroid progenitor-2 cells, in which β-globin expression is similarly affected by triplex formation. These results suggest a model in which RNAs conforming to the strict sequence rules for DNA·RNA triplex formation may participate in feedback regulation of genes in cis.

O-GlcNAc homeostasis contributes to cell fate decisions during hematopoiesis

The addition of a single β-d-GlcNAc sugar (O-GlcNAc) by O-GlcNAc-transferase (OGT) and O-GlcNAc removal by O-GlcNAcase (OGA) maintain homeostatic O-GlcNAc levels on cellular proteins. Changes in protein O-GlcNAcylation regulate cellular differentiation and cell fate decisions, but how these changes affect erythropoiesis, an essential process in blood cell formation, remains unclear. Here, we investigated the role of O-GlcNAcylation in erythropoiesis by using G1E-ER4 cells, which carry the erythroid-specific transcription factor GATA-binding protein 1 (GATA-1) fused to the estrogen receptor (GATA-1-ER) and therefore undergo erythropoiesis after β-estradiol (E2) addition. We observed that during G1E-ER4 differentiation, overall O-GlcNAc levels decrease, and physical interactions of GATA-1 with both OGT and OGA increase. RNA-Seq-based transcriptome analysis of G1E-ER4 cells differentiated in the presence of the OGA inhibitor Thiamet-G (TMG) revealed changes in expression of 433 GATA-1 target genes. ChIP results indicated that the TMG treatment decreases the occupancy of GATA-1, OGT, and OGA at the GATA-binding site of the lysosomal protein transmembrane 5 (Laptm5) gene promoter. TMG also reduced the expression of genes involved in differentiation of NB4 and HL60 human myeloid leukemia cells, suggesting that O-GlcNAcylation is involved in the regulation of hematopoietic differentiation. Sustained treatment of G1E-ER4 cells with TMG before differentiation reduced hemoglobin-positive cells and increased stem/progenitor cell surface markers. Our results show that alterations in O-GlcNAcylation disrupt transcriptional programs controlling erythropoietic lineage commitment, suggesting a role for O-GlcNAcylation in regulating hematopoietic cell fate.

Helix-loop-helix proteins and the advent of cellular diversity: 30 years of discovery

In this review from Murre, the evolution of HLH genes, the structures of HLH domains, and the elaborate activities of HLH proteins in multicellular life are discussed.

Helix–loop–helix (HLH) proteins are dimeric transcription factors that control lineage- and developmental-specific gene programs. Genes encoding for HLH proteins arose in unicellular organisms >600 million years ago and then duplicated and diversified from ancestral genes across the metazoan and plant kingdoms to establish multicellularity. Hundreds of HLH proteins have been identified with diverse functions in a wide variety of cell types. HLH proteins orchestrate lineage specification, commitment, self-renewal, proliferation, differentiation, and homing. HLH proteins also regulate circadian clocks, protect against hypoxic stress, promote antigen receptor locus assembly, and program transdifferentiation. HLH proteins deposit or erase epigenetic marks, activate noncoding transcription, and sequester chromatin remodelers across the chromatin landscape to dictate enhancer–promoter communication and somatic recombination. Here the evolution of HLH genes, the structures of HLH domains, and the elaborate activities of HLH proteins in multicellular life are discussed.

Temporal autoregulation during human PU.1 locus SubTAD formation

Epigenetic control of gene expression occurs within discrete spatial chromosomal units called topologically associating domains (TADs), but the exact spatial requirements of most genes are unknown; this is of particular interest for genes involved in cancer. We therefore applied high-resolution chromosomal conformation capture sequencing to map the three-dimensional (3D) organization of the human locus encoding the key myeloid transcription factor PU.1 in healthy monocytes and acute myeloid leukemia (AML) cells. We identified a dynamic ∼75-kb unit (SubTAD) as the genomic region in which spatial interactions between PU.1 gene regulatory elements occur during myeloid differentiation and are interrupted in AML. Within this SubTAD, proper initiation of the spatial chromosomal interactions requires PU.1 autoregulation and recruitment of the chromatin-adaptor protein LDB1 (LIM domain-binding protein 1). However, once these spatial interactions have occurred, LDB1 stabilizes them independently of PU.1 autoregulation. Thus, our data support that PU.1 autoregulates its expression in a "hit-and-run" manner by initiating stable chromosomal loops that result in a transcriptionally active chromatin architecture.

Divergent Hemogen genes of teleosts and mammals share conserved roles in erythropoiesis: analysis using transgenic and mutant zebrafish

Hemogen is a vertebrate transcription factor that performs important functions in erythropoiesis and testicular development and may contribute to neoplasia. Here we identify zebrafish Hemogen and show that it is considerably smaller (∼22 kDa) than its human ortholog (∼55 kDa), a striking difference that is explained by an underlying modular structure. We demonstrate that Hemogens are largely composed of 21-25 amino acid repeats, some of which may function as transactivation domains (TADs). Hemogen expression in embryonic and adult zebrafish is detected in hematopoietic, renal, neural and gonadal tissues. Using Tol2- and CRISPR/Cas9-generated transgenic zebrafish, we show that Hemogen expression is controlled by two Gata1-dependent regulatory sequences that act alone and together to control spatial and temporal expression during development. Partial depletion of Hemogen in embryos by morpholino knockdown reduces the number of erythrocytes in circulation. CRISPR/Cas9-generated zebrafish lines containing either a frameshift mutation or an in-frame deletion in a putative, C-terminal TAD display anemia and embryonic tail defects. This work expands our understanding of Hemogen and provides mutant zebrafish lines for future study of the mechanism of this important transcription factor.

Summary: Transgenic and mutant zebrafish lines were created to characterize the expression and functions of Hemogen, a transcription factor involved in the formation of red blood cells and other processes.

Strap associates with Csde1 and affects expression of select Csde1-bound transcripts

Erythropoiesis is regulated at many levels, including control of mRNA translation. Changing environmental conditions, such as hypoxia or the availability of nutrients and growth factors, require a rapid response enacted by the enhanced or repressed translation of existing transcripts. Cold shock domain protein e1 (Csde1/Unr) is an RNA-binding protein required for erythropoiesis and strongly upregulated in erythroblasts relative to other hematopoietic progenitors. The aim of this study is to identify the Csde1-containing protein complexes and investigate their role in post-transcriptional expression control of Csde1-bound transcripts. We show that Serine/Threonine kinase receptor-associated protein (Strap/Unrip), was the protein most strongly associated with Csde1 in erythroblasts. Strap is a WD40 protein involved in signaling and RNA splicing, but its role when associated with Csde1 is unknown. Reduced expression of Strap did not alter the pool of transcripts bound by Csde1. Instead, it altered the mRNA and/or protein expression of several Csde1-bound transcripts that encode for proteins essential for translational regulation during hypoxia, such as Hmbs, eIF4g3 and Pabpc4. Also affected by Strap knockdown were Vim, a Gata-1 target crucial for erythrocyte enucleation, and Elavl1, which stabilizes Gata-1 mRNA. The major cellular processes affected by both Csde1 and Strap were ribosome function and cell cycle control.

Necessity and Sufficiency of Ldb1 in the Generation, Differentiation and Maintenance of Non-photoreceptor Cell Types During Retinal Development

During mammalian retinal development, the multipotent progenitors differentiate into all classes of retinal cells under the delicate control of transcriptional factors. The deficiency of a transcription cofactor, the LIM-domain binding protein Ldb1, has been shown to cause proliferation and developmental defects in multiple tissues including cardiovascular, hematopoietic, and nervous systems; however, it remains unclear whether and how it regulates retinal development. By expression profiling, RNA in situ hybridization and immunostaining, here we show that Ldb1 is expressed in the progenitors during early retinal development, but later its expression gradually shifts to non-photoreceptor cell types including bipolar, amacrine, horizontal, ganglion, and Müller glial cells. Retina-specific ablation of Ldb1 in mice resulted in microphthalmia, optic nerve hypoplasia, retinal thinning and detachment, and profound vision impairment as determined by electroretinography. In the mutant retina, there was precocious differentiation of amacrine and horizontal cells, indicating a requirement of Ldb1 in maintaining the retinal progenitor pool. Additionally, all non-photoreceptor cell types were greatly reduced which appeared to be caused by a generation defect and/or retinal degeneration via excessive cell apoptosis. Furthermore, we showed that misexpressed Ldb1 was sufficient to promote the generation of bipolar, amacrine, horizontal, ganglion, and Müller glial cells at the expense of photoreceptors. Together, these results demonstrate that Ldb1 is not only necessary but also sufficient for the development and/or maintenance of non-photoreceptor cell types, and implicate that the pleiotropic functions of Ldb1 during retinal development are context-dependent and determined by its interaction with diverse LIM-HD (LIM-homeodomain) and LMO (LIM domain-only) binding protein partners.

The Co-operation of RUNX1 with LDB1, CDK9 and BRD4 Drives Transcription Factor Complex Relocation During Haematopoietic Specification

Haematopoietic cells arise from endothelial cells within the dorsal aorta of the embryo via a process called the endothelial-haematopoietic transition (EHT). This process crucially depends on the transcription factor RUNX1 which rapidly activates the expression of genes essential for haematopoietic development. Using an inducible version of RUNX1 in a mouse embryonic stem cell differentiation model we showed that prior to the EHT, haematopoietic genes are primed by the binding of the transcription factor FLI1. Once expressed, RUNX1 relocates FLI1 towards its binding sites. However, the nature of the transcription factor assemblies recruited by RUNX1 to reshape the chromatin landscape and initiate mRNA synthesis are unclear. Here, we performed genome-wide analyses of RUNX1-dependent binding of factors associated with transcription elongation to address this question. We demonstrate that RUNX1 induction moves FLI1 from distal ETS/GATA sites to RUNX1/ETS sites and recruits the basal transcription factors CDK9, BRD4, the Mediator complex and the looping factor LDB1. Our study explains how the expression of a single transcription factor can drive rapid and replication independent transitions in cellular shape which are widely observed in development and disease.

Evolution of hemoglobin loci and their regulatory elements

Across the expanse of vertebrate evolution, each species produces multiple forms of hemoglobin in erythroid cells at appropriate times and in the proper amounts. The multiple hemoglobins are encoded in two globin gene clusters in almost all species. One globin gene cluster, linked to the gene NPRL3, is preserved in all vertebrates, including a gene cluster encoding the highly divergent globins from jawless vertebrates. This preservation of synteny may reflect the presence of a powerful enhancer of globin gene expression in the NPRL3 gene. Despite substantial divergence in noncoding DNA sequences among mammals, several epigenetic features of the globin gene regulatory regions are preserved across vertebrates. The preserved features include multiple DNase hypersensitive sites, at least one of which is an enhancer, and binding by key lineage-restricted transcription factors such as GATA1 and TAL1, which in turn recruit coactivators such as P300 that catalyze acetylation of histones. The maps of epigenetic features are strongly correlated with activity in gene regulation, and resources for accessing and visualizing such maps are readily available to the community of researchers and students.

Selective dissociation between LSD1 and GFI1B by a LSD1 inhibitor NCD38 induces the activation of ERG super-enhancer in erythroleukemia cells

Lysine-specific demethylase 1 (LSD1) is a histone modifier for transcriptional repression involved in the regulation of hematopoiesis. We previously reported that a LSD1 inhibitor NCD38 induces transdifferentiation from erythroid lineage to granulomonocytic lineage and exerts anti-leukemia effect through de-repression of the specific super-enhancers of hematopoietic regulators including ERG in a human erythroleukemia cell line, HEL. However, the mechanistic basis for this specificity of NCD38 has remained unclear. Herein, we report major partners associated with LSD1 and clarify the mechanism in HEL cells. Proteome analysis identified 54 candidate proteins associated with LSD1, including several transcription factors such as GFI1B and RUNX1 as well as BRAF-histone deacetylase complex (BHC) components such as CoREST, HDAC1, and HDAC2. NCD38 selectively disrupted the interaction of LSD1 with GFI1B but not with RUNX1, CoREST, HDAC1 and HDAC2. Erg was downregulated in murine erythroid progenitors with prominent upregulation of Gfi1b. NCD38 induced ERG and attenuated an erythroid marker CD235a in HEL while this attenuation was mimicked by the lentiviral overexpression of ERG. The ERG super-enhancer contained the conserved binding motif of GFI1B and was actually occupied by GFI1B. NCD38 dissociated LSD1 and CoREST but not GFI1B from the ERG super-enhancer. Collectively, the selective separation of LSD1 from GFI1B by NCD38 restores the ERG super-enhancer activation and consequently upregulates ERG expression, inducing the transdifferentiation linked to the anti-leukemia effect.

The LDB1 Complex Co-opts CTCF for Erythroid Lineage-Specific Long-Range Enhancer Interactions

Lineage-specific transcription factors are critical for long-range enhancer interactions, but direct or indirect contributions of architectural proteins such as CCCTC-binding factor (CTCF) to enhancer function remain less clear. The LDB1 complex mediates enhancer-gene interactions at the β-globin locus through LDB1 self-interaction. We find that an LDB1-bound enhancer upstream of carbonic anhydrase 2 (Car2) activates its expression by interacting directly with CTCF at the gene promoter. Both LDB1 and CTCF are required for enhancer-Car2 looping, and the domain of LDB1 contacted by CTCF is necessary to rescue Car2 transcription in LDB1-deficient cells. Genome-wide studies and CRISPR/Cas9 genome editing indicate that LDB1-CTCF enhancer looping underlies activation of a substantial fraction of erythroid genes. Our results provide a mechanism by which long-range interactions of architectural protein CTCF can be tailored to achieve a tissue-restricted pattern of chromatin loops and gene expression.

Dynamics of promoter bivalency and RNAP II pausing in mouse stem and differentiated cells

Background

Mammalian embryonic stem cells display a unique epigenetic and transcriptional state to facilitate pluripotency by maintaining lineage-specification genes in a poised state. Two epigenetic and transcription processes involved in maintaining poised state are bivalent chromatin, characterized by the simultaneous presence of activating and repressive histone methylation marks, and RNA polymerase II (RNAPII) promoter proximal pausing. However, the dynamics of histone modifications and RNAPII at promoters in diverse cellular contexts remains underexplored.

Results

We collected genome wide data for bivalent chromatin marks H3K4me3 and H3K27me3, and RNAPII (8WG16) occupancy together with expression profiling in eight different cell types, including ESCs, in mouse. The epigenetic and transcription profiles at promoters grouped in over thirty clusters with distinct functional identities and transcription control.

Conclusion

The clustering analysis identified distinct bivalent clusters where genes in one cluster retained bivalency across cell types while in the other were mostly cell type specific, but neither showed a high RNAPII pausing. We noted that RNAPII pausing is more associated with active genes than bivalent genes in a cell type, and was globally reduced in differentiated cell types compared to multipotent.

Electronic supplementary material

The online version of this article (10.1186/s12861-018-0163-7) contains supplementary material, which is available to authorized users.