GATA-1 mediates auto-regulation of Gfi-1B transcription in K562 cells

Machine-guided design of synthetic cell type-specific cis-regulatory elements

Cis-regulatory elements (CREs) control gene expression, orchestrating tissue identity, developmental timing, and stimulus responses, which collectively define the thousands of unique cell types in the body. While there is great potential for strategically incorporating CREs in therapeutic or biotechnology applications that require tissue specificity, there is no guarantee that an optimal CRE for an intended purpose has arisen naturally through evolution. Here, we present a platform to engineer and validate synthetic CREs capable of driving gene expression with programmed cell type specificity. We leverage innovations in deep neural network modeling of CRE activity across three cell types, efficient in silico optimization, and massively parallel reporter assays (MPRAs) to design and empirically test thousands of CREs. Through in vitro and in vivo validation, we show that synthetic sequences outperform natural sequences from the human genome in driving cell type-specific expression. Synthetic sequences leverage unique sequence syntax to promote activity in the on-target cell type and simultaneously reduce activity in off-target cells. Together, we provide a generalizable framework to prospectively engineer CREs and demonstrate the required literacy to write regulatory code that is fit-for-purpose in vivo across vertebrates.

Static electromagnetic field and recombinant human fibroblasts encoding miR-451 and miR-16 increased cell trans-differentiation to CD71+ and CD235a+ erythroid like progenitor

Introduction

Ex vivo blood production is an urgent need of most countries, and creating production protocols can save the lives of many patients. Despite the recent advances in blood production in ex vivo conditions, its high-scale production is not yet possible, and requires further studies. Therefore, by transfecting fibroblast cells with miR-16, and miR-451 genes, as well as applying low frequency electromagnetic fields (ELF-EMF) treatment, we tried to increase the differentiation of these cells into CD71+ and CD235a+ erythroid like progenitors.

Methods

After preparation, and cultivation of human dermal transgenic fibroblast cells, they were transfected by Plenti3-hsa-miR451, Plenti3-hsa-miR16 and Plenti3-backbone inserted into E. coli Stbl4 genome. Then, transgenic fibroblast cells were treated with 10mT ELF-EMF every day for 20 minutes for 7 days. Using a flow cytometer, the expressions of CD71, and CD235a were studied in these cells, and the expressions of genes involved in hematopoiesis were studied using the RT-PCR technique.

Results

The results indicated an increase in the differentiation of fibroblast cells treated with 10mT ELF-EMF to erythroid like progenitors. Furthermore, the percentage of CD71+ and CD235a+ cells was the highest in irradiated cells encoding miR-16 and miR-451, which indicates their differentiation into erythroid like progenitors. Also, in the transgenic cells treated with ELF-EMF, an increase in the expressions of α-chain, β-chain, γ-chain and GATA1 genes was observed, which indicates the potential of these cells for hematopoiesis. However, there was no significant difference in the expression of CD34 and CD38 genes in these cell lines.

Conclusion

Both ELF-EMF and upregulations of miR-16 and miR-451 lead to improved differentiation of fibroblast cells into erythroid like progenitors.

LSD1 inhibition modulates transcription factor networks in myeloid malignancies

Acute Myeloid Leukemia (AML) is a type of cancer of the blood system that is characterized by an accumulation of immature hematopoietic cells in the bone marrow and blood. Its pathogenesis is characterized by an increase in self-renewal and block in differentiation in hematopoietic stem and progenitor cells. Underlying its pathogenesis is the acquisition of mutations in these cells. As there are many different mutations found in AML that can occur in different combinations the disease is very heterogeneous. There has been some progress in the treatment of AML through the introduction of targeted therapies and a broader application of the stem cell transplantation in its treatment. However, many mutations found in AML are still lacking defined interventions. These are in particular mutations and dysregulation in important myeloid transcription factors and epigenetic regulators that also play a crucial role in normal hematopoietic differentiation. While a direct targeting of the partial loss-of-function or change in function observed in these factors is very difficult to imagine, recent data suggests that the inhibition of LSD1, an important epigenetic regulator, can modulate interactions in the network of myeloid transcription factors and restore differentiation in AML. Interestingly, the impact of LSD1 inhibition in this regard is quite different between normal and malignant hematopoiesis. The effect of LSD1 inhibition involves transcription factors that directly interact with LSD1 such as GFI1 and GFI1B, but also transcription factors that bind to enhancers that are modulated by LSD1 such as PU.1 and C/EBPα as well as transcription factors that are regulated downstream of LSD1 such as IRF8. In this review, we are summarizing the current literature on the impact of LSD1 modulation in normal and malignant hematopoietic cells and the current knowledge how the involved transcription factor networks are altered. We are also exploring how these modulation of transcription factors play into the rational selection of combination partners with LSD1 inhibitors, which is an intense area of clinical investigation.

TF-Prioritizer: a Java pipeline to prioritize condition-specific transcription factors

BACKGROUND

Eukaryotic gene expression is controlled by cis-regulatory elements (CREs), including promoters and enhancers, which are bound by transcription factors (TFs). Differential expression of TFs and their binding affinity at putative CREs determine tissue- and developmental-specific transcriptional activity. Consolidating genomic datasets can offer further insights into the accessibility of CREs, TF activity, and, thus, gene regulation. However, the integration and analysis of multimodal datasets are hampered by considerable technical challenges. While methods for highlighting differential TF activity from combined chromatin state data (e.g., chromatin immunoprecipitation [ChIP], ATAC, or DNase sequencing) and RNA sequencing data exist, they do not offer convenient usability, have limited support for large-scale data processing, and provide only minimal functionality for visually interpreting results.

RESULTS

We developed TF-Prioritizer, an automated pipeline that prioritizes condition-specific TFs from multimodal data and generates an interactive web report. We demonstrated its potential by identifying known TFs along with their target genes, as well as previously unreported TFs active in lactating mouse mammary glands. Additionally, we studied a variety of ENCODE datasets for cell lines K562 and MCF-7, including 12 histone modification ChIP sequencing as well as ATAC and DNase sequencing datasets, where we observe and discuss assay-specific differences.

CONCLUSION

TF-Prioritizer accepts ATAC, DNase, or ChIP sequencing and RNA sequencing data as input and identifies TFs with differential activity, thus offering an understanding of genome-wide gene regulation, potential pathogenesis, and therapeutic targets in biomedical research.

The scaffolding function of LSD1/KDM1A reinforces a negative feedback loop to repress stem cell gene expression during primitive hematopoiesis

Lsd1/Kdm1a functions both as a histone demethylase enzyme and as a scaffold for assembling chromatin modifier and transcription factor complexes to regulate gene expression. The relative contributions of Lsd1's demethylase and scaffolding functions during embryogenesis are not known. Here, we analyze two independent zebrafish lsd1/kdm1a mutant lines and show Lsd1 is required to repress primitive hematopoietic stem cell gene expression. Lsd1 rescue constructs containing point mutations that selectively abrogate its demethylase or scaffolding capacity demonstrate the scaffolding function of Lsd1, not its demethylase activity, is required for repression of gene expression in vivo. Lsd1's SNAG-binding domain mediates its scaffolding function and reinforces a negative feedback loop to repress the expression of SNAG-domain-containing genes during embryogenesis, including gfi1 and snai1/2. Our findings reveal a model in which the SNAG-binding and scaffolding function of Lsd1, and its associated negative feedback loop, provide transient and reversible regulation of gene expression during hematopoietic development.

GFI1B acts as a metabolic regulator in hematopoiesis and acute myeloid leukemia

Recent studies highlighted the role of transcription factors in metabolic regulation during hematopoiesis and leukemia development. GFI1B is a transcriptional repressor that plays a critical role in hematopoiesis, and its expression is negatively related to the prognosis of acute myeloid leukemia (AML) patients. We earlier reported a change in the metabolic state of hematopoietic stem cells upon Gfi1b deletion. Here we explored the role of Gfi1b in metabolism reprogramming during hematopoiesis and leukemogenesis. We demonstrated that Gfi1b deletion remarkably activated mitochondrial respiration and altered energy metabolism dependence toward oxidative phosphorylation (OXPHOS). Mitochondrial substrate dependency was shifted from glucose to fatty acids upon Gfi1b deletion via upregulating fatty acid oxidation (FAO). On a molecular level, Gfi1b epigenetically regulated multiple FAO-related genes. Moreover, we observed that metabolic phenotypes evolved as cells progressed from preleukemia to leukemia, and the correlation between Gfi1b expression level and metabolic phenotype was affected by genetic variations in AML cells. FAO or OXPHOS inhibition significantly impeded leukemia progression of Gfi1b-KO MLL/AF9 cells. Finally, we showed that Gfi1b-deficient AML cells were more sensitive to metformin as well as drugs implicated in OXPHOS and FAO inhibition, opening new potential therapeutic strategies.

Characterization of a genomic region 8 kb downstream of GFI1B associated with myeloproliferative neoplasms

A genomic locus 8 kb downstream of the transcription factor GFI1B (Growth Factor Independence 1B) predisposes to clonal hematopoiesis and myeloproliferative neoplasms. One of the most significantly associated polymorphisms in this region is rs621940-G. GFI1B auto-represses GFI1B, and altered GFI1B expression contributes to myeloid neoplasms. We studied whether rs621940-G affects GFI1B expression and growth of immature cells. GFI1B ChIP-seq showed clear binding to the rs621940 locus. Preferential binding of various hematopoietic transcription factors to either the rs621940-C or -G allele was observed, but GFI1B showed no preference. In gene reporter assays the rs621940 region inhibited GFI1B promoter activity with the G-allele having less suppressive effects compared to the C-allele. However, CRISPR-Cas9 mediated deletion of the locus in K562 cells did not alter GFI1B expression nor auto-repression. In healthy peripheral blood mononuclear cells GFI1B expression did not differ consistently between the rs621940 alleles. Long range and targeted deep sequencing did not detect consistent effects of rs621940-G on allelic GFI1B expression either. Finally, we observed that myeloid colony formation was not significantly affected by either rs621940 allele in 193 healthy donors. Together, these findings show no evidence that rs621940 or its locus affect GFI1B expression, auto-repression or growth of immature myeloid cells.

Motif-Raptor: a cell type-specific and transcription factor centric approach for post-GWAS prioritization of causal regulators

MOTIVATION

Genome-wide association studies (GWASs) have identified thousands of common trait-associated genetic variants but interpretation of their function remains challenging. These genetic variants can overlap the binding sites of transcription factors (TFs) and therefore could alter gene expression. However, we currently lack a systematic understanding on how this mechanism contributes to phenotype.

RESULTS

We present Motif-Raptor, a TF-centric computational tool that integrates sequence-based predictive models, chromatin accessibility, gene expression datasets and GWAS summary statistics to systematically investigate how TF function is affected by genetic variants. Given trait-associated non-coding variants, Motif-Raptor can recover relevant cell types and critical TFs to drive hypotheses regarding their mechanism of action. We tested Motif-Raptor on complex traits such as rheumatoid arthritis and red blood cell count and demonstrated its ability to prioritize relevant cell types, potential regulatory TFs and non-coding SNPs which have been previously characterized and validated.

AVAILABILITY AND IMPLEMENTATION

Motif-Raptor is freely available as a Python package at: https://github.com/pinellolab/MotifRaptor.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The transcription factors GFI1 and GFI1B as modulators of the innate and acquired immune response

GFI1 and GFI1B are small nuclear proteins of 45 and 37kDa, respectively, that have a simple two-domain structure: The first consists of a group of six c-terminal C₂H₂ zinc finger motifs that are almost identical in sequence and bind to very similar, specific DNA sites. The second is an N-terminal 20 amino acid SNAG domain that can bind to the pocket of the histone demethylase KDM1A (LSD1) near its active site. When bound to DNA, both proteins act as bridging factors that bring LSD1 and associated proteins into the vicinity of methylated substrates, in particular histone H3 or TP53. GFI1 can also bring methyl transferases such as PRMT1 together with its substrates that include the DNA repair proteins MRE11 and 53BP1, thereby enabling their methylation and activation. While GFI1B is expressed almost exclusively in the erythroid and megakaryocytic lineage, GFI1 has clear biological roles in the development and differentiation of lymphoid and myeloid immune cells. GFI1 is required for lymphoid/myeloid and monocyte/granulocyte lineage decision as well as the correct nuclear interpretation of a number of important immune-signaling pathways that are initiated by NOTCH1, interleukins such as IL2, IL4, IL5 or IL7, by the pre TCR or -BCR receptors during early lymphoid differentiation or by T and B cell receptors during activation of lymphoid cells. Myeloid cells also depend on GFI1 at both stages of early differentiation as well as later stages in the process of activation of macrophages through Toll-like receptors in response to pathogen-associated molecular patterns. The knowledge gathered on these factors over the last decades puts GFI1 and GFI1B at the center of many biological processes that are critical for both the innate and acquired immune system.

Multifaceted Actions of GFI1 and GFI1B in Hematopoietic Stem Cell Self-Renewal and Lineage Commitment

Growth factor independence 1 (GFI1) and the closely related protein GFI1B are small nuclear proteins that act as DNA binding transcriptional repressors. Both recognize the same consensus DNA binding motif via their C-terminal zinc finger domains and regulate the expression of their target genes by recruiting chromatin modifiers such as histone deacetylases (HDACs) and demethylases (LSD1) by using an N-terminal SNAG domain that comprises only 20 amino acids. The only region that is different between both proteins is the region that separates the zinc finger domains and the SNAG domain. Both proteins are co-expressed in hematopoietic stem cells (HSCs) and, to some extent, in multipotent progenitors (MPPs), but expression is specified as soon as early progenitors and show signs of lineage bias. While expression of GFI1 is maintained in lymphoid primed multipotent progenitors (LMPPs) that have the potential to differentiate into both myeloid and lymphoid cells, GFI1B expression is no longer detectable in these cells. By contrast, GFI1 expression is lost in megakaryocyte precursors (MKPs) and in megakaryocyte-erythrocyte progenitors (MEPs), which maintain a high level of GFI1B expression. Consequently, GFI1 drives myeloid and lymphoid differentiation and GFI1B drives the development of megakaryocytes, platelets, and erythrocytes. How such complementary cell type- and lineage-specific functions of GFI1 and GFI1B are maintained is still an unresolved question in particular since they share an almost identical structure and very similar biochemical modes of actions. The cell type-specific accessibility of GFI1/1B binding sites may explain the fact that very similar transcription factors can be responsible for very different transcriptional programming. An additional explanation comes from recent data showing that both proteins may have additional non-transcriptional functions. GFI1 interacts with a number of proteins involved in DNA repair and lack of GFI1 renders HSCs highly susceptible to DNA damage-induced death and restricts their proliferation. In contrast, GFI1B binds to proteins of the beta-catenin/Wnt signaling pathway and lack of GFI1B leads to an expansion of HSCs and MKPs, illustrating the different impact that GFI1 or GFI1B has on HSCs. In addition, GFI1 and GFI1B are required for endothelial cells to become the first blood cells during early murine development and are among those transcription factors needed to convert adult endothelial cells or fibroblasts into HSCs. This role of GFI1 and GFI1B bears high significance for the ongoing effort to generate hematopoietic stem and progenitor cells de novo for the autologous treatment of blood disorders such as leukemia and lymphoma.

Learning common and specific patterns from data of multiple interrelated biological scenarios with matrix factorization

High-throughput biological technologies (e.g. ChIP-seq, RNA-seq and single-cell RNA-seq) rapidly accelerate the accumulation of genome-wide omics data in diverse interrelated biological scenarios (e.g. cells, tissues and conditions). Integration and differential analysis are two common paradigms for exploring and analyzing such data. However, current integrative methods usually ignore the differential part, and typical differential analysis methods either fail to identify combinatorial patterns of difference or require matched dimensions of the data. Here, we propose a flexible framework CSMF to combine them into one paradigm to simultaneously reveal Common and Specific patterns via Matrix Factorization from data generated under interrelated biological scenarios. We demonstrate the effectiveness of CSMF with four representative applications including pairwise ChIP-seq data describing the chromatin modification map between K562 and Huvec cell lines; pairwise RNA-seq data representing the expression profiles of two different cancers; RNA-seq data of three breast cancer subtypes; and single-cell RNA-seq data of human embryonic stem cell differentiation at six time points. Extensive analysis yields novel insights into hidden combinatorial patterns in these multi-modal data. Results demonstrate that CSMF is a powerful tool to uncover common and specific patterns with significant biological implications from data of interrelated biological scenarios.

Growth Factor Independence 1B-Mediated Transcriptional Repression and Lineage Allocation Require Lysine-Specific Demethylase 1-Dependent Recruitment of the BHC Complex

Growth factor independence 1B (GFI1B) coordinates assembly of transcriptional repressor complexes comprised of corepressors and histone-modifying enzymes to control gene expression programs governing lineage allocation in hematopoiesis. Enforced expression of GFI1B in K562 erythroleukemia cells favors erythroid over megakaryocytic differentiation, providing a platform to define molecular determinants of binary fate decisions triggered by GFI1B. We deployed proteome-wide proximity labeling to identify factors whose inclusion in GFI1B complexes depends upon GFI1B's obligate effector, lysine-specific demethylase 1 (LSD1). We show that GFI1B preferentially recruits core and putative elements of the BRAF-histone deacetylase (HDAC) (BHC) chromatin-remodeling complex (LSD1, RCOR1, HMG20A, HMG20B, HDAC1, HDAC2, PHF21A, GSE1, ZMYM2, and ZNF217) in an LSD1-dependent manner to control acquisition of erythroid traits by K562 cells. Among these elements, depletion of both HMG20A and HMG20B or of GSE1 blocks GFI1B-mediated erythroid differentiation, phenocopying impaired differentiation brought on by LSD1 depletion or disruption of GFI1B-LSD1 binding. These findings demonstrate the central role of the GFI1B-LSD1 interaction as a determinant of BHC complex recruitment to enable cell fate decisions driven by GFI1B.

Gfi1aa and Gfi1b set the pace for primitive erythroblast differentiation from hemangioblasts in the zebrafish embryo

The transcriptional repressors Gfi1(a) and Gfi1b are epigenetic regulators with unique and overlapping roles in hematopoiesis. In different contexts, Gfi1 and Gfi1b restrict or promote cell proliferation, prevent apoptosis, influence cell fate decisions, and are essential for terminal differentiation. Here, we show in primitive red blood cells (prRBCs) that they can also set the pace for cellular differentiation. In zebrafish, prRBCs express 2 of 3 zebrafish Gfi1/1b paralogs, Gfi1aa and Gfi1b. The recently identified zebrafish gfi1aa gene trap allele qmc551 drives erythroid green fluorescent protein (GFP) instead of Gfi1aa expression, yet homozygous carriers have normal prRBCs. prRBCs display a maturation defect only after splice morpholino-mediated knockdown of Gfi1b in gfi1aa qmc551 homozygous embryos. To study the transcriptome of the Gfi1aa/1b double-depleted cells, we performed an RNA-Seq experiment on GFP-positive prRBCs sorted from 20-hour-old embryos that were heterozygous or homozygous for gfi1aa qmc551 , as well as wt or morphant for gfi1b We subsequently confirmed and extended these data in whole-mount in situ hybridization experiments on newly generated single- and double-mutant embryos. Combined, the data showed that in the absence of Gfi1aa, the synchronously developing prRBCs were delayed in activating late erythroid differentiation, as they struggled to suppress early erythroid and endothelial transcription programs. The latter highlighted the bipotent nature of the progenitors from which prRBCs arise. In the absence of Gfi1aa, Gfi1b promoted erythroid differentiation as stepwise loss of wt gfi1b copies progressively delayed Gfi1aa-depleted prRBCs even further, showing that Gfi1aa and Gfi1b together set the pace for prRBC differentiation from hemangioblasts.

Genetic Modifiers of the Breast Tumor Microenvironment

Multiple nonmalignant cell types in the tumor microenvironment (TME) impact breast cancer risk, metastasis, and response to therapy, yet most heritable mechanisms that influence TME cell function and breast cancer outcomes are largely unknown. Breast cancer risk is ∼30% heritable and >170 genetic loci have been associated with breast cancer traits. However, the majority of candidate genes have poorly defined mechanistic roles in breast cancer biology. Research indicates that breast cancer risk modifiers directly impact cancer cells, yet it is equally plausible that some modifier alleles impact the nonmalignant TME. The objective of this review is to examine the list of current breast cancer candidate genes that may modify breast cancer risk and outcome through the TME.

SMIM1 variants rs1175550 and rs143702418 independently modulate Vel blood group antigen expression

The Vel blood group antigen is expressed on the red blood cells of most individuals. Recently, we described that homozygosity for inactivating mutations in SMIM1 defines the rare Vel-negative phenotype. Still, Vel-positive individuals show great variability in Vel antigen expression, creating a risk for Vel blood typing errors and transfusion reactions. We fine-mapped the regulatory region located in SMIM1 intron 2 in Swedish blood donors, and observed a strong correlation between expression and rs1175550 as well as with a previously unreported tri-nucleotide insertion (rs143702418; C > CGCA). While the two variants are tightly linked in Caucasians, we separated their effects in African Americans, and found that rs1175550G and to a lesser extent rs143702418C independently increase SMIM1 and Vel antigen expression. Gel shift and luciferase assays indicate that both variants are transcriptionally active, and we identified binding of the transcription factor TAL1 as a potential mediator of the increased expression associated with rs1175550G. Our results provide insight into the regulatory logic of Vel antigen expression, and extend the set of markers for genetic Vel blood group typing.

A Novel LSD1 Inhibitor T-3775440 Disrupts GFI1B-Containing Complex Leading to Transdifferentiation and Impaired Growth of AML Cells

Dysregulation of lysine (K)-specific demethylase 1A (LSD1), also known as KDM1A, has been implicated in the development of various cancers, including leukemia. Here, we describe the antileukemic activity and mechanism of action of T-3775440, a novel irreversible LSD1 inhibitor. Cell growth analysis of leukemia cell lines revealed that acute erythroid leukemia (AEL) and acute megakaryoblastic leukemia cells (AMKL) were highly sensitive to this compound. T-3775440 treatment enforced transdifferentiation of erythroid/megakaryocytic lineages into granulomonocytic-like lineage cells. Mechanistically, T-3775440 disrupted the interaction between LSD1 and growth factor-independent 1B (GFI1B), a transcription factor critical for the differentiation processes of erythroid and megakaryocytic lineage cells. Knockdown of LSD1 and GFI1B recapitulated T-3775440-induced transdifferentiation and cell growth suppression, highlighting the significance of LSD1-GFI1B axis inhibition with regard to the anti-AML effects of T-3775440. Moreover, T-3775440 exhibited significant antitumor efficacy in AEL and AMKL xenograft models. Our findings provide a rationale for evaluating LSD1 inhibitors as potential treatments and indicate a novel mechanism of action against AML, particularly AEL and AMKL. Mol Cancer Ther; 16(2); 273-84. ©2016 AACR.

New insights into the regulation by RUNX1 and GFI1(s) proteins of the endothelial to hematopoietic transition generating primordial hematopoietic cells

The first hematopoietic cells are generated very early in ontogeny to support the growth of the embryo and to provide the foundation to the adult hematopoietic system. There is a considerable therapeutic interest in understanding how these first blood cells are generated in order to try to reproduce this process in vitro. This would allow generating blood products, or hematopoietic cell populations from embryonic stem (ES) cells, induced pluripotent stem cells or through directed reprogramming. Recent studies have clearly established that the first hematopoietic cells originate from a hemogenic endothelium (HE) through an endothelial to hematopoietic transition (EHT). The molecular mechanisms underlining this transition remain largely unknown with the exception that the transcription factor RUNX1 is critical for this process. In this Extra Views report, we discuss our recent studies demonstrating that the transcriptional repressors GFI1 and GFI1B have a critical role in the EHT. We established that these RUNX1 transcriptional targets are actively implicated in the downregulation of the endothelial program and the loss of endothelial identity during the formation of the first blood cells. In addition, our results suggest that GFI1 expression provides an ideal novel marker to identify, isolate and study the HE cell population.

ConBind: motif-aware cross-species alignment for the identification of functional transcription factor binding sites

Eukaryotic gene expression is regulated by transcription factors (TFs) binding to promoter as well as distal enhancers. TFs recognize short, but specific binding sites (TFBSs) that are located within the promoter and enhancer regions. Functionally relevant TFBSs are often highly conserved during evolution leaving a strong phylogenetic signal. While multiple sequence alignment (MSA) is a potent tool to detect the phylogenetic signal, the current MSA implementations are optimized to align the maximum number of identical nucleotides. This approach might result in the omission of conserved motifs that contain interchangeable nucleotides such as the ETS motif (IUPAC code: GGAW). Here, we introduce ConBind, a novel method to enhance alignment of short motifs, even if their mutual sequence similarity is only partial. ConBind improves the identification of conserved TFBSs by improving the alignment accuracy of TFBS families within orthologous DNA sequences. Functional validation of the Gfi1b + 13 enhancer reveals that ConBind identifies additional functionally important ETS binding sites that were missed by all other tested alignment tools. In addition to the analysis of known regulatory regions, our web tool is useful for the analysis of TFBSs on so far unknown DNA regions identified through ChIP-sequencing.

A somatic mutation of GFI1B identified in leukemia alters cell fate via a SPI1 (PU.1) centered genetic regulatory network

We identify a mutation (D262N) in the erythroid-affiliated transcriptional repressor GFI1B, in an acute myeloid leukemia (AML) patient with antecedent myelodysplastic syndrome (MDS). The GFI1B-D262N mutant functionally antagonizes the transcriptional activity of wild-type GFI1B. GFI1B-D262N promoted myelomonocytic versus erythroid output from primary human hematopoietic precursors and enhanced cell survival of both normal and MDS derived precursors. Re-analysis of AML transcriptome data identifies a distinct group of patients in whom expression of wild-type GFI1B and SPI1 (PU.1) have an inverse pattern. In delineating this GFI1B-SPI1 relationship we show that (i) SPI1 is a direct target of GFI1B, (ii) expression of GFI1B-D262N produces elevated expression of SPI1, and (iii) SPI1-knockdown restores balanced lineage output from GFI1B-D262N-expressing precursors. These results table the SPI1-GFI1B transcriptional network as an important regulatory axis in AML as well as in the development of erythroid versus myelomonocytic cell fate.

Growth factor independence 1b (gfi1b) is important for the maturation of erythroid cells and the regulation of embryonic globin expression

Growth factor independence 1b (GFI1B) is a DNA binding repressor of transcription with vital functions in hematopoiesis. Gfi1b-null embryos die at midgestation very likely due to defects in erythro- and megakaryopoiesis. To analyze the full functionality of Gfi1b, we used conditionally deficient mice that harbor floxed Gfi1b alleles and inducible (Mx-Cre, Cre-ERT) or erythroid specific (EpoR-Cre) Cre expressing transgenes. In contrast to the germline knockout, EpoR-Cre mediated erythroid specific ablation of Gfi1b allows full gestation, but causes perinatal lethality with very few mice surviving to adulthood. Both the embryonic deletion of Gfi1b by EpoR-Cre and the deletion in adult mice by Mx-Cre or Cre-ERT leads to reduced numbers of erythroid precursors, perturbed and delayed erythroid maturation, anemia and extramedullary erythropoiesis. Global expression analyses showed that the Hba-x, Hbb-bh1 and Hbb-y embryonic globin genes were upregulated in Gfi1b deficient TER119⁺ fetal liver cells over the gestation period from day 12.5–17.5 p.c. and an increased level of Hbb-bh1 and Hbb-y embryonic globin gene expression was even maintained in adult Gfi1b deficient mice. While the expression of Bcl11a, a regulator of embryonic globin expression was not affected by Gfi1b deficiency, the expression of Gata1 was reduced and the expression of Sox6, also involved in globin switch, was almost entirely lost when Gfi1b was absent. These findings establish Gfi1b as a regulator of embryonic globin expression and embryonic and adult erythroid maturation.

Comparative annotation of functional regions in the human genome using epigenomic data

Epigenetic regulation is dynamic and cell-type dependent. The recently available epigenomic data in multiple cell types provide an unprecedented opportunity for a comparative study of epigenetic landscape. We developed a machine-learning method called ChroModule to annotate the epigenetic states in eight ENCyclopedia Of DNA Elements cell types. The trained model successfully captured the characteristic histone-modification patterns associated with regulatory elements, such as promoters and enhancers, and showed superior performance on identifying enhancers compared with the state-of-art methods. In addition, given the fixed number of epigenetic states in the model, ChroModule allows straightforward illustration of epigenetic variability in multiple cell types. Using this feature, we found that invariable and variable epigenetic states across cell types correspond to housekeeping functions and stimulus response, respectively. Especially, we observed that enhancers, but not the other regulatory elements, dictate cell specificity, as similar cell types share common enhancers, and cell-type-specific enhancers are often bound by transcription factors playing critical roles in that cell type. More interestingly, we found some genomic regions are dormant in cell type but primed to become active in other cell types. These observations highlight the usefulness of ChroModule in comparative analysis and interpretation of multiple epigenomes.

Predicting cell-type-specific gene expression from regions of open chromatin

Complex patterns of cell-type-specific gene expression are thought to be achieved by combinatorial binding of transcription factors (TFs) to sequence elements in regulatory regions. Predicting cell-type-specific expression in mammals has been hindered by the oftentimes unknown location of distal regulatory regions. To alleviate this bottleneck, we used DNase-seq data from 19 diverse human cell types to identify proximal and distal regulatory elements at genome-wide scale. Matched expression data allowed us to separate genes into classes of cell-type-specific up-regulated, down-regulated, and constitutively expressed genes. CG dinucleotide content and DNA accessibility in the promoters of these three classes of genes displayed substantial differences, highlighting the importance of including these aspects in modeling gene expression. We associated DNase I hypersensitive sites (DHSs) with genes, and trained classifiers for different expression patterns. TF sequence motif matches in DHSs provided a strong performance improvement in predicting gene expression over the typical baseline approach of using proximal promoter sequences. In particular, we achieved competitive performance when discriminating up-regulated genes from different cell types or genes up- and down-regulated under the same conditions. We identified previously known and new candidate cell-type-specific regulators. The models generated testable predictions of activating or repressive functions of regulators. DNase I footprints for these regulators were indicative of their direct binding to DNA. In summary, we successfully used information of open chromatin obtained by a single assay, DNase-seq, to address the problem of predicting cell-type-specific gene expression in mammalian organisms directly from regulatory sequence.

PhosphoSiteAnalyzer: a bioinformatic platform for deciphering phospho proteomes using kinase predictions retrieved from NetworKIN

Phosphoproteomic experiments are routinely conducted in laboratories worldwide, and because of the fast development of mass spectrometric techniques and efficient phosphopeptide enrichment methods, researchers frequently end up having lists with tens of thousands of phosphorylation sites for further interrogation. To answer biologically relevant questions from these complex data sets, it becomes essential to apply computational, statistical, and predictive analytical methods. Here we provide an advanced bioinformatic platform termed "PhosphoSiteAnalyzer" to explore large phosphoproteomic data sets that have been subjected to kinase prediction using the previously published NetworKIN algorithm. NetworKIN applies sophisticated linear motif analysis and contextual network modeling to obtain kinase-substrate associations with high accuracy and sensitivity. PhosphoSiteAnalyzer provides an algorithm to retrieve kinase predictions from the public NetworKIN webpage in a semiautomated way and applies hereafter advanced statistics to facilitate a user-tailored in-depth analysis of the phosphoproteomic data sets. The interface of the software provides a high degree of analytical flexibility and is designed to be intuitive for most users. PhosphoSiteAnalyzer is a freeware program available at http://phosphosite.sourceforge.net .

Gfi1b negatively regulates Rag expression directly and via the repression of FoxO1

Gfi1b negatively regulates Rag expression through direct binding to the Rag locus and through inhibition of Foxo1; mice lacking both Gfi1b and Gfi1 exhibit a block in B cell development.

Precise regulation of Rag (recombination-activating gene) expression is crucial to prevent genomic instability caused by the generation of Rag-mediated DNA breaks. Although mechanisms of Rag activation have been well characterized, the mechanism by which Rag expression is down-regulated in early B cell development has not been fully elucidated. Using a complementary DNA library screen, we identified the transcriptional repressor Gfi1b as negative regulator of the Rag locus. Expression of Gfi1b causes repression of Rag1 and Rag2 in cell lines and primary mouse cells. Conversely, Gfi1b-deficient cell lines exhibit increased Rag expression, double-strand breaks and recombination, and cell cycle defects. In primary cells, transcription of Gfi1b inversely correlates with Rag transcription, and simultaneous inactivation of Gfi1 and Gfi1b leads to an increase in Rag transcription early in B cell development. In addition, deletion of Gfi1 and Gfi1b in vivo results in a severe block in B cell development. Gfi1b orchestrates Rag repression via a dual mechanism. Direct binding of Gfi1b to a site 5′ of the B cell–specific Erag enhancer results in epigenetic changes in the Rag locus, whereas indirect inhibition is achieved through repression of the trans-activator Foxo1. Together, our experiments show that Gfi family members are essential for normal B cell development and play an important role in modulating expression of the V(D)J recombinase.

Human promoter mutations unveil Oct-1 and GATA-1 opposite action on Gfi1b regulation

Growth factor-independence 1b (Gfi1b) is a zinc finger transcription factor essential for erythroid and megakaryocytic development. To better understand Gfi1b regulation and to know the implication of the level of expression of this gene in human pathology, we have searched for promoter punctual sequence variations in 214 patients with different hematological diseases. We found two previously unknown congenital mutations at evolutionary conserved GATA and octamer-binding (Oct) transcription factor sites. The Oct site mutation was also found in five relatives of the patient. The GATA motif mutation reduced promoter activity by 50% in vitro, while homozygous patients with the octamer site mutation showed a four-to-five times increase of Gfi1b RNA in platelets. Electrophoretic mobility shift analyses demonstrated that different protein complexes bind to both sites and that binding is reduced by the mutations. Finally, we found that GATA-1 and Oct-1 are the main components of each complex. This study provides evidences of a new mechanism for Gfi1b repression. This is also the first report of Gfi1b mutations with a functional implication; further investigation and follow-up will clarify the involvement of these mutations in hematological disease.

Gfi-1B promoter remains associated with active chromatin marks throughout erythroid differentiation of human primary progenitor cells

Growth Factor Independent-1B (Gfi-1B) is a transcriptional repressor that plays critical roles in the control of erythropoiesis and megakaryopoiesis. Gfi-1B expression was described to be repressed by an autoregulatory feedback control loop. Here, we show that Gfi-1 transcription is positively regulated early after induction of erythroid differentiation and remains highly active to late erythroblasts. Using chromatin immunoprecipitation assays in CD34⁺ cells from human cord blood, we found that Gfi-1 and GATA-2 in immature progenitors and then Gfi-1B and GATA-1 in erythroblasts are bound to the Gfi-1B promoter as well as to the promoter of c-myc, a known Gfi-1B target gene. Surprisingly, this Gfi-1/GATA-2–Gfi-1B/GATA-1 switch observed at erythroblast stages is associated to an increase in the Gfi-1B transcription whereas it triggers repression of c-myc transcription. Accordingly, analysis of chromatin modification patterns shows that HDAC, CoREST, and LSD1 are recruited to the c-myc promoter leading to appearance of repressive chromatin marks. In contrast, the Gfi-1B promoter remains associated with a transcriptionally active chromatin configuration as highlighted by an increase in histone H3 acetylation and concomitant release of the LSD1 and CoREST corepressors. The repressive function of Gfi-1B therefore depends on the nature of the proteins recruited to the target gene promoters and on chromatin modifications. We conclude that Gfi-1B behaves as a lineage-affiliated gene with an open chromatin configuration in multipotent progenitors and sustained activation as cells progress throughout erythroid differentiation.

High-mobility group protein HMGB2 regulates human erythroid differentiation through trans-activation of GFI1B transcription

Gfi-1B is a transcriptional repressor that is crucial for erythroid differentiation: inactivation of the GFI1B gene in mice leads to embryonic death due to failure to produce differentiated red cells. Accordingly, GFI1B expression is tightly regulated during erythropoiesis, but the mechanisms involved in such regulation remain partially understood. We here identify HMGB2, a high-mobility group HMG protein, as a key regulator of GFI1B transcription. HMGB2 binds to the GFI1B promoter in vivo and up-regulates its trans-activation most likely by enhancing the binding of Oct-1 and, to a lesser extent, of GATA-1 and NF-Y to the GFI1B promoter. HMGB2 expression increases during erythroid differentiation concomitantly to the increase of GfI1B transcription. Importantly, knockdown of HMGB2 in immature hematopoietic progenitor cells leads to decreased Gfi-1B expression and impairs their erythroid differentiation. We propose that HMGB2 potentiates GATA-1-dependent transcription of GFI1B by Oct-1 and thereby controls erythroid differentiation.

Erythropoiesis: model systems, molecular regulators, and developmental programs

Human erythropoiesis is a complex multistep developmental process that begins at the level of pluripotent hematopoietic stem cells (HSCs) at bone marrow microenvironment (HSCs niche) and terminates with the production of erythrocytes (RBCs). This review covers the basic and contemporary aspects of erythropoiesis. These include the: (a) cell-lineage restricted pathways of differentiation originated from HSCs and going downward toward the blood cell development; (b) model systems employed to study erythropoiesis in culture (erythroleukemia cell lines and embryonic stem cells) and in vivo (knockout animals: avian, mice, zebrafish, and xenopus); (c) key regulators of erythropoiesis (iron, hypoxia, stress, and growth factors); (d) signaling pathways operating at hematopoietic stem cell niche for homeostatic regulation of self renewal (SCF/c-kit receptor, Wnt, Notch, and Hox) and for erythroid differentiation (HIF and EpoR). Furthermore, this review presents the mechanisms through which transcriptional factors (GATA-1, FOG-1, TAL-1/SCL/MO2/Ldb1/E2A, EKLF, Gfi-1b, and BCL11A) and miRNAs regulate gene pattern expression during erythroid differentiation. New insights regarding the transcriptional regulation of alpha- and beta-globin gene clusters were also presented. Emphasis was also given on (i) the developmental program of erythropoiesis, which consists of commitment to terminal erythroid maturation and hemoglobin production, (two closely coordinated events of erythropoieis) and (ii) the capacity of human embryonic and umbilical cord blood (UCB) stem cells to differentiate and produce RBCs in culture with highly selective media. These most recent developments will eventually permit customized red blood cell production needed for transfusion.

GFI1B controls its own expression binding to multiple sites

BACKGROUND

Transcription factors play essential roles in both normal and malignant hematopoiesis. This is the case for the growth factor independent 1b (GFI1B) transcription factor, which is required for erythroid and megakaryocytic differentiation and over-expressed in leukemic patients and cell lines.

DESIGN AND METHODS

To investigate GFI1B regulation, we searched for multispecies conserved non-coding elements between GFI1B and neighboring genes. We used a formaldehyde-assisted isolation of regulatory elements (FAIRE) assay and DNase1 hypersensitivity to assess the chromatin conformation of these sites. Next, we analyzed transcription factor binding and histone modifications at the GFI1B locus including the conserved non-coding elements by a chromatin immunoprecipitation assay. Finally, we studied the interaction of the GFI1B promoter and the conserved non-coding elements with the chromatin conformation capture technique and used immunofluorescence to evaluate GFI1B levels in individual cells.

RESULTS

We localized several conserved non-coding elements containing multiple erythroid specific transcription factor binding sites at the GFI1B locus. In GFI1B-expressing cells a subset of these conserved non-coding elements and the promoter adopt a close spatial conformation, localize with open chromatin sites, harbor chromatin modifications associated with gene activation and bind multiple transcription factors and co-repressors. Conclusions Our findings indicate that GFI1B regulatory elements behave as activators and repressors. Different protein levels within a cell population suggest that cells must activate and repress GFI1B continuously to control its final level. These data are consistent with a model of GFI1B regulation in which GFI1B binds to its own promoter and to the conserved non-coding elements as its levels rise. This would attract repressor complexes that progressively down-regulate the gene. GFI1B expression would decrease until a stage at which the activating complexes predominate and expression increases.

Growth factor independent 1b (Gfi1b) and a new splice variant of Gfi1b are highly expressed in patients with acute and chronic leukemia

Gfi1b is a transcriptional repressor that is essential for erythroid cells and megakaryocytes, but is also expressed in hematopoietic stem cells and early myeloid progenitors. The chromosomal localization of the Gfi1b gene at 9q34 and its functional homology with the proto-oncogene Gfi1 were suggestive for a role of Gfi1b in malignant transformation and myeloid leukemia. We show here that the expression of Gfi1b is strongly elevated in CML and AML patients compared to normal healthy controls and that imatinib, a drug widely used to treat CML, further enhances Gfi1b expression in patients even after remission. Our data suggest that Gfi1b may be an important factor to establish or maintain myeloid leukemia and myeloproliferative diseases and that, high expression levels of Gfi1b might be associated with the emergence of Philadelphia chromosome negative myeloid malignancies after imatinib withdrawal or after the development of imatinib resistance.

Phosphorylation site mapping of endogenous proteins: a combined MS and bioinformatics approach

We present a novel approach that combines MALDI-TOF profile analysis and bioinformatics-based inclusion criteria to comprehensively predict phosphorylation sites on a single protein of interest from limiting sample. It is technologically difficult to unambiguously identify phosphorylated residues, as many physiologically important phosphorylation sites are of too low abundance in vivo to be unambiguously assigned by mass spectrometry. Conversely, phosphorylation site prediction algorithms, while increasingly accurate, nevertheless overestimate the number of phosphorylation sites. In this study, we show that MODICAS, an MS data management and analysis tool, can be effectively merged with the bioinformatics attributes of residue conservation and phosphosite prediction to generate a short list of putative phosphorylation sites that can be subsequently verified by additional methodologies such as phosphospecific antibodies or mutational analysis. Therefore, the combination of MODICAS driven MS data analysis with bioinformatics-based filtering represents a substantial increase in the ability to putatively identify physiologically relevant phosphosites from limited starting material.

A HaemAtlas: characterizing gene expression in differentiated human blood cells

Hematopoiesis is a carefully controlled process that is regulated by complex networks of transcription factors that are, in part, controlled by signals resulting from ligand binding to cell-surface receptors. To further understand hematopoiesis, we have compared gene expression profiles of human erythroblasts, megakaryocytes, B cells, cytotoxic and helper T cells, natural killer cells, granulocytes, and monocytes using whole genome microarrays. A bioinformatics analysis of these data was performed focusing on transcription factors, immunoglobulin superfamily members, and lineage-specific transcripts. We observed that the numbers of lineage-specific genes varies by 2 orders of magnitude, ranging from 5 for cytotoxic T cells to 878 for granulocytes. In addition, we have identified novel coexpression patterns for key transcription factors involved in hematopoiesis (eg, GATA3-GFI1 and GATA2-KLF1). This study represents the most comprehensive analysis of gene expression in hematopoietic cells to date and has identified genes that play key roles in lineage commitment and cell function. The data, which are freely accessible, will be invaluable for future studies on hematopoiesis and the role of specific genes and will also aid the understanding of the recent genome-wide association studies.

Regulatory factors for initial T lymphocyte lineage specification

PURPOSE OF REVIEW

Initiation of T lymphocyte development depends on balanced regulatory inputs from multiple essential transcription factors. This review highlights contributions of E2A, hematopoietic transcription factor PU.1, growth factor independence (Gfi)-1, T cell factor (TCF)-1, and Runx factors and their interactions with the Notch pathway to promote T cell development.

RECENT FINDINGS

E2A and Runx family factors have been implicated in establishing competent precursors in which Notch signaling can induce the T cell program. An early role was also indicated for PU.1. Later PU.1 activities are antagonistic to pro-T cell factors, however, including E proteins, Myb, Gfi-1, and TCF-1. Diversion to a non-T lineage can be promoted by PU.1, CCAAT/enhancer binding protein, or even GATA and TCF, but these diversion mechanisms are blocked by Notch signaling. An emergent gene network summarizes the cross-regulatory relationships among these factors.

SUMMARY

Entry into the T-cell pathway is controlled by a dynamic balance among essential regulatory factors that depend on Notch signaling not only to trigger initiation of the T-cell program but also to maintain the lineage fidelity of their collective action.

From hematopoietic stem cells to platelets

Megakaryocytopoiesis is the process that leads to the production of platelets. This process involves the commitment of multipotent hematopoietic stem cells toward megakaryocyte (MK) progenitors, the proliferation and differentiation of MK progenitors, the polyploidization of MK precursors and the maturation of MK. Mature MK produce platelets by cytoplasmic fragmentation occurring through a dynamic and regulated process, called proplatelet formation, and consisting of long pseudopodial elongations that break in the blood flow. Recent insights have demonstrated that the MK and erythroid lineages are tightly associated at both the cellular and molecular levels, especially in the transcription factors that regulate their differentiation programs. Megakaryocytopoiesis is regulated by two types of transcription factors, those regulating the differentiation process, such as GATA-1, and those regulating proplatelet formation, such as NF-E2. The humoral factor thrombopoietin (TPO) is the primary regulator of MK differentiation and platelet production through the stimulation of its receptor MPL. Numerous acquired or congenital pathologies of the MK lineage are now explained by molecular abnormalities in the activity of the transcription factors involved in megakaryocytopoiesis, in the Tpo or c-mpl genes, as well as in signaling molecules associated with MPL. The recent development of MPL agonists may provide efficient agents for the treatment of some thrombocytopenias.

GATA-1 and Gfi-1B interplay to regulate Bcl-xL transcription

The induction of Bcl-x(L) is critical for the survival of late proerythroblasts. The erythroid-specific transcriptional network that regulates Bcl-x(L) expression in erythropoiesis remains unclear. The activation of the central erythropoietic transcriptional factor, GATA-1, leads to the early, transient induction of a transcription repressor, Gfi-1B, followed by the late induction of Bcl-x(L) during erythroid maturation in G1ER cells. Chromatin immunoprecipitation assays demonstrated that a constant level of GATA-1 binds to the Bcl-x promoter throughout the entire induction period, while Gfi-1B is transiently associated with the promoter in the early phase. The sustained expression of Gfi-1B abolished GATA-1-induced Bcl-x(L) expression. Here, we present evidence that GATA-1 binds to the noncanonical GATT motif of the Bcl-x promoter for trans-activation. Gfi-1B expressed at increased levels is recruited to the Bcl-x promoter through its association with GATA-1, suppressing Bcl-x(L) transcription. Therefore, the down-regulation of Gfi-1B in the late phase of erythroid maturation is necessary for Bcl-x(L) induction. Furthermore, we show that the inhibition of Bcr-Abl kinase by treatment with imatinib caused the up-regulation of Gfi-1B in K562 cells, where Gfi-1B also cooperated with GATA-1 to repress Bcl-x(L) transcription. Gfi-1B knockdown by RNA interference diminished imatinib-induced apoptosis, while the overexpression of Gfi-1B sensitized K562 cells to arsenic-induced death. These findings illuminate the role of Gfi-1B in GATA-1-mediated transcription in the survival aspect of erythroid cells.

Ovol1 represses its own transcription by competing with transcription activator c-Myb and by recruiting histone deacetylase activity

Ovol1 belongs to a family of evolutionarily conserved zinc finger proteins that act downstream of key developmental signaling pathways such as Wnt and TGF-β/BMP. It plays important roles in epithelial and germ cell development, particularly by repressing c-Myc and Id2 genes and modulating the balance between proliferation and differentiation of progenitor cells. In this study, we show that Ovol1 negatively regulates its own expression by binding to and repressing the activity of its promoter. We further demonstrate that Ovol1 uses both passive and active repression mechanisms to auto-repress: (1) it antagonizes transcriptional activation of c-Myb, a known positive regulator of proliferation, by competing for DNA binding; (2) it recruits histone deacetylase activity to the promoter via an N-terminal SNAG repressor domain. At Ovol1 cognate sites in the endogenous Ovol1 promoter, c-Myb binding correlates with increased histone acetylation, whereas the expression of Ovol1 correlates with a displacement of c-Myb from the DNA and decreased histone acetylation. Collectively, our data suggest that Ovol1 restricts its own expression by counteracting c-Myb activation and histone acetylation of the Ovol1 promoter.

Growth factor independent 1B (Gfi1b) is an E2A target gene that modulates Gata3 in T-cell lymphomas

The E2A transcription factors are required for normal T lymphopoiesis and to prevent T-lymphocyte progenitor transformation. Ectopic expression of E2A proteins in E2A-deficient lymphomas results in growth arrest and apoptosis, indicating that these cells remain responsive to the targets of E2A. Here we identify the transcriptional repressor growth factor independent 1B (Gfi1b) as a target of E2A that promotes growth arrest and apoptosis in lymphomas. Gfi1b expression in primary T-lymphocyte progenitors is dependent on E2A and excess Gfi1b prevents the outgrowth of T lymphocyte progenitors in vitro. Gfi1b represses expression of Gata3, a transcription factor whose appropriate regulation is required for survival of lymphomas and T-lymphocyte progenitors. We also show that ectopic expression of Gata3 in lymphomas promotes expression of Gfi1b, indicating that these proteins may function in an autoregulatory loop that maintains appropriate levels of Gata3. Therefore, we propose that E2A proteins prevent lymphoma cell expansion, at least in part through regulation of Gfi1b and modulation of Gata3 expression.