Erythroid GATA1 function revealed by genome-wide analysis of transcription factor occupancy, histone modifications, and mRNA expression

TFTG: A comprehensive database for human transcription factors and their targets

Transcription factors (TFs) are major contributors to gene transcription, especially in controlling cell-specific gene expression and disease occurrence and development. Uncovering the relationship between TFs and their target genes is critical to understanding the mechanism of action of TFs. With the development of high-throughput sequencing techniques, a large amount of TF-related data has accumulated, which can be used to identify their target genes. In this study, we developed TFTG (Transcription Factor and Target Genes) database (http://tf.liclab.net/TFTG), which aimed to provide a large number of available human TF-target gene resources by multiple strategies, besides performing a comprehensive functional and epigenetic annotations and regulatory analyses of TFs. We identified extensive available TF-target genes by collecting and processing TF-associated ChIP-seq datasets, perturbation RNA-seq datasets and motifs. We also obtained experimentally confirmed relationships between TF and target genes from available resources. Overall, the target genes of TFs were obtained through integrating the relevant data of various TFs as well as fourteen identification strategies. Meanwhile, TFTG was embedded with user-friendly search, analysis, browsing, downloading and visualization functions. TFTG is designed to be a convenient resource for exploring human TF-target gene regulations, which will be useful for most users in the TF and gene expression regulation research.

Gapped-kmer sequence modeling robustly identifies regulatory vocabularies and distal enhancers conserved between evolutionarily distant mammals

Gene regulatory elements drive complex biological phenomena and their mutations are associated with common human diseases. The impacts of human regulatory variants are often tested using model organisms such as mice. However, mapping human enhancers to conserved elements in mice remains a challenge, due to both rapid enhancer evolution and limitations of current computational methods. We analyze distal enhancers across 45 matched human/mouse cell/tissue pairs from a comprehensive dataset of DNase-seq experiments, and show that while cell-specific regulatory vocabulary is conserved, enhancers evolve more rapidly than promoters and CTCF binding sites. Enhancer conservation rates vary across cell types, in part explainable by tissue specific transposable element activity. We present an improved genome alignment algorithm using gapped-kmer features, called gkm-align, and make genome wide predictions for 1,401,803 orthologous regulatory elements. We show that gkm-align discovers 23,660 novel human/mouse conserved enhancers missed by previous algorithms, with strong evidence of conserved functional activity.

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

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

BIT: Bayesian Identification of Transcriptional Regulators

BIT is a novel Bayesian hierarchical model capable of predicting transcriptional regulators (TRs) from the input of user-provided epigenomic regions. TRs are critical molecules in transcriptional regulation. Many diseases and cancers are linked to the dysfunction of TRs. Knowing TRs in certain biological process can help find new biomarkers or therapeutic targets. Thus, BIT formulates a novel Bayesian hierarchical model with the Pólya-gamma data augmentation strategy. Based on collected ChIP-seq datasets, BIT can identify TRs responsible for the genome-wide binding pattern within the user-provided epigenomic regions. BIT has been validated by using a simulation study and three applications.

Functional dissection of complex and molecular trait variants at single nucleotide resolution

Identifying the causal variants and mechanisms that drive complex traits and diseases remains a core problem in human genetics. The majority of these variants have individually weak effects and lie in non-coding gene-regulatory elements where we lack a complete understanding of how single nucleotide alterations modulate transcriptional processes to affect human phenotypes. To address this, we measured the activity of 221,412 trait-associated variants that had been statistically fine-mapped using a Massively Parallel Reporter Assay (MPRA) in 5 diverse cell-types. We show that MPRA is able to discriminate between likely causal variants and controls, identifying 12,025 regulatory variants with high precision. Although the effects of these variants largely agree with orthogonal measures of function, only 69% can plausibly be explained by the disruption of a known transcription factor (TF) binding motif. We dissect the mechanisms of 136 variants using saturation mutagenesis and assign impacted TFs for 91% of variants without a clear canonical mechanism. Finally, we provide evidence that epistasis is prevalent for variants in close proximity and identify multiple functional variants on the same haplotype at a small, but important, subset of trait-associated loci. Overall, our study provides a systematic functional characterization of likely causal common variants underlying complex and molecular human traits, enabling new insights into the regulatory grammar underlying disease risk.

HEXIM1 is an essential transcription regulator during human erythropoiesis

ABSTRACT

Regulation of RNA polymerase II (RNAPII) activity is an essential process that governs gene expression; however, its contribution to the fundamental process of erythropoiesis remains unclear. hexamethylene bis-acetamide inducible 1 (HEXIM1) regulates RNAPII activity by controlling the location and activity of positive transcription factor β. We identified a key role for HEXIM1 in controlling erythroid gene expression and function, with overexpression of HEXIM1 promoting erythroid proliferation and fetal globin expression. HEXIM1 regulated erythroid proliferation by enforcing RNAPII pausing at cell cycle check point genes and increasing RNAPII occupancy at genes that promote cycle progression. Genome-wide profiling of HEXIM1 revealed that it was increased at both repressed and activated genes. Surprisingly, there were also genome-wide changes in the distribution of GATA-binding factor 1 (GATA1) and RNAPII. The most dramatic changes occurred at the β-globin loci, where there was loss of RNAPII and GATA1 at β-globin and gain of these factors at γ-globin. This resulted in increased expression of fetal globin, and BGLT3, a long noncoding RNA in the β-globin locus that regulates fetal globin expression. GATA1 was a key determinant of the ability of HEXIM1 to repress or activate gene expression. Genes that gained both HEXIM1 and GATA1 had increased RNAPII and increased gene expression, whereas genes that gained HEXIM1 but lost GATA1 had an increase in RNAPII pausing and decreased expression. Together, our findings reveal a central role for universal transcription machinery in regulating key aspects of erythropoiesis, including cell cycle progression and fetal gene expression, which could be exploited for therapeutic benefit.

Defining a cohort of anemia-activated cis elements reveals a mechanism promoting erythroid precursor function

Acute anemia elicits broad transcriptional changes in erythroid progenitors and precursors. We previously discovered a cis-regulatory transcriptional enhancer at the sterile alpha motif domain-14 enhancer locus (S14E), defined by a CANNTG-spacer-AGATAA composite motif and occupied by GATA1 and TAL1 transcription factors, is required for survival in severe anemia. However, S14E is only 1 of dozens of anemia-activated genes containing similar motifs. In a mouse model of acute anemia, we identified populations of expanding erythroid precursors, which increased expression of genes that contain S14E-like cis elements. We reveal that several S14E-like cis elements provide important transcriptional control of newly identified anemia-inducing genes, including the Ssx-2 interacting protein (Ssx2ip). Ssx2ip expression was determined to play an important role in erythroid progenitor/precursor cell activities, cell cycle regulation, and cell proliferation. Over a weeklong course of acute anemia recovery, we observed that erythroid gene activation mediated by S14E-like cis elements occurs during a phase coincident with low hematocrit and high progenitor activities, with distinct transcriptional programs activated at earlier and later time points. Our results define a genome-wide mechanism in which S14E-like enhancers control transcriptional responses during erythroid regeneration. These findings provide a framework to understand anemia-specific transcriptional mechanisms, ineffective erythropoiesis, anemia recovery, and phenotypic variability within human populations.

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.

The role of GATA switch in benzene metabolite hydroquinone inhibiting erythroid differentiation in K562 cells

The phenolic metabolite of benzene, hydroquinone (HQ), has potential risks for hematological disorders and hematotoxicity in humans. Previous studies have revealed that reactive oxygen species, DNA methylation, and histone acetylation participate in benzene metabolites inhibiting erythroid differentiation in hemin-induced K562 cells. GATA1 and GATA2 are crucial erythroid-specific transcription factors that exhibit dynamic expression patterns during erythroid differentiation. We investigated the role of GATA factors in HQ-inhibited erythroid differentiation in K562 cells. When K562 cells were induced with 40 μM hemin for 0-120 h, the mRNA and protein levels of GATA1 and GATA2 changed dynamically. After exposure to 40 μM HQ for 72 h, K562 cells were induced with 40 μM hemin for 48 h. HQ considerably reduced the percentage of hemin-induced Hb-positive cells, decreased the GATA1 mRNA, protein, and occupancy levels at α-globin and β-globin gene clusters, and increased the GATA2 mRNA and protein levels significantly. ChIP-seq analysis revealed that HQ reduced GATA1 occupancy, and increased GATA2 occupancy at most gene loci in hemin-induced K562 cells. And GATA1 and GATA2 might play essential roles in the erythroid differentiation protein interaction network. These results elucidate that HQ decreases GATA1 occupancy and increases GATA2 occupancy at the erythroid gene loci, thereby downregulating GATA1 and upregulating GATA2 expression, which in turn modulates the expression of erythroid genes and inhibits erythroid differentiation. This partially explains the mechanism of benzene hematotoxicity.

Toward a comprehensive catalog of regulatory elements

Regulatory elements are the genomic regions that interact with transcription factors to control cell-type-specific gene expression in different cellular environments. A precise and complete catalog of functional elements encoded by the human genome is key to understanding mammalian gene regulation. Here, we review the current state of regulatory element annotation. We first provide an overview of assays for characterizing functional elements, including genome, epigenome, transcriptome, three-dimensional chromatin interaction, and functional validation assays. We then discuss computational methods for defining regulatory elements, including peak-calling and other statistical modeling methods. Finally, we introduce several high-quality lists of regulatory element annotations and suggest potential future directions.

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.

Pathogenic human variant that dislocates GATA2 zinc fingers disrupts hematopoietic gene expression and signaling networks

Although certain human genetic variants are conspicuously loss of function, decoding the impact of many variants is challenging. Previously, we described a patient with leukemia predisposition syndrome (GATA2 deficiency) with a germline GATA2 variant that inserts 9 amino acids between the 2 zinc fingers (9aa-Ins). Here, we conducted mechanistic analyses using genomic technologies and a genetic rescue system with Gata2 enhancer-mutant hematopoietic progenitor cells to compare how GATA2 and 9aa-Ins function genome-wide. Despite nuclear localization, 9aa-Ins was severely defective in occupying and remodeling chromatin and regulating transcription. Variation of the inter-zinc finger spacer length revealed that insertions were more deleterious to activation than repression. GATA2 deficiency generated a lineage-diverting gene expression program and a hematopoiesis-disrupting signaling network in progenitors with reduced granulocyte-macrophage colony-stimulating factor (GM-CSF) and elevated IL-6 signaling. As insufficient GM-CSF signaling caused pulmonary alveolar proteinosis and excessive IL-6 signaling promoted bone marrow failure and GATA2 deficiency patient phenotypes, these results provide insight into mechanisms underlying GATA2-linked pathologies.

Hematopoietic/erythroid enhancers activate nearby target genes by extending histone H3K27ac and transcribing intergenic RNA

Enhancers activate gene transcription remotely, which requires tissue specific transcription factors binding to them. GATA1 and TAL1 are hematopoietic/erythroid-specific factors and often bind together to enhancers, activating target genes. Interestingly, we found that some hematopoietic/erythroid genes are transcribed in a GATA1-dependent but TAL1-independnet manner. They appear to have enhancers within a relatively short distance. In this study, we paired highly transcribed hematopoietic/erythroid genes with the nearest GATA1/TAL1-binding enhancers and analyzed these putative enhancer-gene pairs depending on distance between them. Enhancers located at various distances from genes in the pairs, which was not related to transcription level of the genes. However, genes with enhancers at short distances away tended to be transcriptionally unaffected by TAL1 depletion. Histone H3K27ac extended from the enhancers to target genes. The H3K27ac extension was maintained without TAL1, even though it disappeared owing to the loss of GATA1. Intergenic RNA was highly transcribed from the enhancers to nearby target genes, independent of TAL1. Taken together, TAL1-independent transcription of hematopoietic/erythroid genes appears to be promoted by enhancers present in a short distance. These enhancers are likely to activate nearby target genes by tracking the intervening regions.

Regulation of Cell Plasticity by Bromodomain and Extraterminal Domain (BET) Proteins: A New Perspective in Glioblastoma Therapy

BET proteins are a family of multifunctional epigenetic readers, mainly involved in transcriptional regulation through chromatin modelling. Transcriptome handling ability of BET proteins suggests a key role in the modulation of cell plasticity, both in fate decision and in lineage commitment during embryonic development and in pathogenic conditions, including cancerogenesis. Glioblastoma is the most aggressive form of glioma, characterized by a very poor prognosis despite the application of a multimodal therapy. Recently, new insights are emerging about the glioblastoma cellular origin, leading to the hypothesis that several putative mechanisms occur during gliomagenesis. Interestingly, epigenome dysregulation associated with loss of cellular identity and functions are emerging as crucial features of glioblastoma pathogenesis. Therefore, the emerging roles of BET protein in glioblastoma onco-biology and the compelling demand for more effective therapeutic strategies suggest that BET family members could be promising targets for translational breakthroughs in glioblastoma treatment. Primarily, “Reprogramming Therapy”, which is aimed at reverting the malignant phenotype, is now considered a promising strategy for GBM therapy.

BET proteins: Biological functions and therapeutic interventions

Bromodomain and extra-terminal (BET) family member proteins (BRD2, BRD3, BRD4 and BRDT) play a pivotal role in interpreting the epigenetic information of histone Kac modification, thus controlling gene expression, remodeling chromatin structures and avoid replicative stress-induced DNA damages. Abnormal activation of BET proteins is tightly correlated to various human diseases, including cancer. Therefore, BET bromodomain inhibitors (BBIs) were considered as promising therapeutics to treat BET-related diseases, raising >70 clinical trials in the past decades. Despite preliminary effects achieved, drug resistance and adverse events represent two major challenges for current BBIs development. In this review, we will introduce the biological functions of BET proteins in both physiological and pathological conditions; and summarize the progress in current BBI drug development. Moreover, we will also discuss the major challenges in the front of BET inhibitor development and provide rational strategies to overcome these obstacles.

RNA-regulatory exosome complex suppresses an apoptotic program to confer erythroid progenitor cell survival in vivo

The RNA-regulatory exosome complex (EC) posttranscriptionally and cotranscriptionally processes and degrades RNAs in a context-dependent manner. Although the EC functions in diverse cell types, its contributions to stem and progenitor cell development are not well understood. Previously, we demonstrated that the transcriptional regulator of erythrocyte development, GATA1, represses EC subunit genes, and the EC maintains erythroid progenitors in vitro. To determine if this mechanism operates in vivo, we used the hematopoietic-specific Vav1-Cre and "conditional by inversion" mouse system to ablate Exosc3, encoding an EC structural subunit. Although Exosc3C/C Cre+ embryos developed normally until embryonic day 14.5, Exosc3 ablation was embryonic lethal and severely reduced erythromyeloid progenitor activity. RNA sequencing analysis of Exosc3-ablated burst-forming unit-erythroid revealed elevated transcripts encoding multiple proapoptotic factors, and the mutant erythroid progenitors exhibited increased apoptosis. We propose that the EC controls an ensemble of apoptosis-regulatory RNAs, thereby promoting erythroid progenitor survival and developmental erythropoiesis in vivo.

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.

Beyond neurodegenerative diseases: α-synuclein in erythropoiesis

α-synuclein (α-syn) is a highly conserved and thermostable protein that is widely distributed in human brain. An intracellular aggregation of α-syn in dopaminergic neurons is the hallmark of a group of neurodegenerative diseases including Parkinson's disease. Interestingly, α-syn is also highly expressed in red blood cells and is considered as one of the most abundant proteins in red blood cells. Moreover, α-syn is thought to play a regulatory role during normal erythropoiesis. However, whether α-syn participates in the pathogenesis of erythroid diseases has not been reported. In this review, we discuss the protein structure of α-syn and the importance of α-syn in erythropoiesis.

Molecular and cellular mechanisms that regulate human erythropoiesis

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

HDAC1 is required for GATA-1 transcription activity, global chromatin occupancy and hematopoiesis

The activity of hematopoietic factor GATA-1 is modulated through p300/CBP-mediated acetylation and FOG-1 mediated indirect interaction with HDAC1/2 containing NuRD complex. Although GATA-1 acetylation is implicated in GATA-1 activation, the role of deacetylation is not studied. Here, we found that the FOG-1/NuRD does not deacetylate GATA-1. However, HDAC1/2 can directly bind and deacetylate GATA-1. Two arginine residues within the GATA-1 linker region mediates direct interaction with HDAC1. The arginine to alanine mutation (2RA) blocks GATA-1 deacetylation and fails to induce erythroid differentiation. Gene expression profiling and ChIP-seq analysis further demonstrate the importance of GATA-1 deacetylation for gene activation and chromatin recruitment. GATA-12RA knock-in (KI) mice suffer mild anemia and thrombocytopenia with accumulation of immature erythrocytes and megakaryocytes in bone marrow and spleen. Single cell RNA-seq analysis of Lin- cKit+ (LK) cells further reveal a profound change in cell subpopulations and signature gene expression patterns in HSC, myeloid progenitors, and erythroid/megakaryocyte clusters in KI mice. Thus, GATA-1 deacetylation and its interaction with HDAC1 modulates GATA-1 chromatin binding and transcriptional activity that control erythroid/megakaryocyte commitment and differentiation.

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.

Discovering How Heme Controls Genome Function Through Heme-omics

Protein ensembles control genome function by establishing, maintaining, and deconstructing cell-type-specific chromosomal landscapes. A plethora of small molecules orchestrate cellular functions and therefore may link physiological processes with genome biology. The metabolic enzyme and hemoglobin cofactor heme induces proteolysis of a transcriptional repressor, Bach1, and regulates gene expression post-transcriptionally. However, whether heme controls genome function broadly or through prescriptive actions is unclear. Using assay for transposase-accessible chromatin sequencing (ATAC-seq), we establish a heme-dependent chromatin atlas in wild-type and mutant erythroblasts lacking enhancers that confer normal heme synthesis. Amalgamating chromatin landscapes and transcriptomes in cells with sub-physiological heme and post-heme rescue reveals parallel Bach1-dependent and Bach1-independent mechanisms that target heme-sensing chromosomal hotspots. The hotspots harbor a DNA motif demarcating heme-regulated chromatin and genes encoding proteins not known to be heme regulated, including metabolic enzymes. The heme-omics analysis establishes how an essential biochemical cofactor controls genome function and cellular physiology.

FBXO11-mediated proteolysis of BAHD1 relieves PRC2-dependent transcriptional repression in erythropoiesis

The histone mark H3K27me3 and its reader/writer polycomb repressive complex 2 (PRC2) mediate widespread transcriptional repression in stem and progenitor cells. Mechanisms that regulate this activity are critical for hematopoietic development but are poorly understood. Here we show that the E3 ubiquitin ligase F-box only protein 11 (FBXO11) relieves PRC2-mediated repression during erythroid maturation by targeting its newly identified substrate bromo adjacent homology domain-containing 1 (BAHD1), an H3K27me3 reader that recruits transcriptional corepressors. Erythroblasts lacking FBXO11 are developmentally delayed, with reduced expression of maturation-associated genes, most of which harbor bivalent histone marks at their promoters. In FBXO11-/- erythroblasts, these gene promoters bind BAHD1 and fail to recruit the erythroid transcription factor GATA1. The BAHD1 complex interacts physically with PRC2, and depletion of either component restores FBXO11-deficient erythroid gene expression. Our studies identify BAHD1 as a novel effector of PRC2-mediated repression and reveal how a single E3 ubiquitin ligase eliminates PRC2 repression at many developmentally poised bivalent genes during erythropoiesis.

MEIRLOP: improving score-based motif enrichment by incorporating sequence bias covariates

BACKGROUND

Motif enrichment analysis (MEA) identifies over-represented transcription factor binding (TF) motifs in the DNA sequence of regulatory regions, enabling researchers to infer which transcription factors can regulate transcriptional response to a stimulus, or identify sequence features found near a target protein in a ChIP-seq experiment. Score-based MEA determines motifs enriched in regions exhibiting extreme differences in regulatory activity, but existing methods do not control for biases in GC content or dinucleotide composition. This lack of control for sequence bias, such as those often found in CpG islands, can obscure the enrichment of biologically relevant motifs.

RESULTS

We developed Motif Enrichment In Ranked Lists of Peaks (MEIRLOP), a novel MEA method that determines enrichment of TF binding motifs in a list of scored regulatory regions, while controlling for sequence bias. In this study, we compare MEIRLOP against other MEA methods in identifying binding motifs found enriched in differentially active regulatory regions after interferon-beta stimulus, finding that using logistic regression and covariates improves the ability to call enrichment of ISGF3 binding motifs from differential acetylation ChIP-seq data compared to other methods. Our method achieves similar or better performance compared to other methods when quantifying the enrichment of TF binding motifs from ENCODE TF ChIP-seq datasets. We also demonstrate how MEIRLOP is broadly applicable to the analysis of numerous types of NGS assays and experimental designs.

CONCLUSIONS

Our results demonstrate the importance of controlling for sequence bias when accurately identifying enriched DNA sequence motifs using score-based MEA. MEIRLOP is available for download from https://github.com/npdeloss/meirlop under the MIT license.

A Unique Epigenomic Landscape Defines Human Erythropoiesis

Mammalian erythropoiesis yields a highly specialized cell type, the mature erythrocyte, evolved to meet the organismal needs of increased oxygen-carrying capacity. To better understand the regulation of erythropoiesis, we performed genome-wide studies of chromatin accessibility, DNA methylation, and transcriptomics using a recently developed strategy to obtain highly purified populations of primary human erythroid cells. The integration of gene expression, DNA methylation, and chromatin state dynamics reveals that stage-specific gene regulation during erythropoiesis is a stepwise and hierarchical process involving many cis-regulatory elements. Erythroid-specific, nonpromoter sites of chromatin accessibility are linked to erythroid cell phenotypic variation and inherited disease. Comparative analyses of stage-specific chromatin accessibility indicate that there is limited early chromatin priming of erythroid genes during hematopoiesis. The epigenome of terminally differentiating erythroid cells defines a distinct subset of highly specialized cells that are vastly dissimilar from other hematopoietic and nonhematopoietic cell types. These epigenomic and transcriptome data are powerful tools to study human erythropoiesis.

Schulz et al. use genome-wide studies of chromatin accessibility, DNA methylation, and transcriptomes in primary human erythroid cells to reveal important characteristics of erythropoiesis. Chromatin accessibility of terminal erythroid differentiation is markedly dissimilar from other hematopoietic cell types. Epigenomic changes are linked to erythroid cell traits and disease genes.

miR-144/451 inhibits c-Myc to promote erythroid differentiation

Ablation of miR-144/451 disrupts homeostasis of erythropoiesis. Myc, a protooncogenic protein, is essential for erythroblast proliferation but commits rapid downregulation during erythroid maturation. How erythroblasts orchestrate maturation processes through coding and non-coding genes is largely unknown. In this study, we use miR-144/451 knockout mice as in vivo model, G1E, MEL erythroblast lines and erythroblasts from fresh mouse fetal livers as in vitro systems to demonstrate that targeted depletion of miR-144/451 blocks erythroid nuclear condensation and enucleation. This is due, at least in part, to the continued high expression of Myc in erythroblasts when miR-144/451 is absent. Specifically, miR-144/451 directly inhibits Myc in erythroblasts. Loss of miR-144/451 locus derepresses, and thus, increases the expression of Myc. Sustained high levels of Myc in miR-144/451-depleted erythroblasts blocks erythroid differentiation. Moreover, Myc reversely regulates the expression of miR-144/451, forming a positive miR-144/451-Myc feedback to ensure the complete shutoff of Myc during erythropoiesis. Given that erythroid-specific transcription factor GATA1 activates miR-144/451 and inactivates Myc, our findings indicate that GATA1-miR-144/451-Myc network safeguards normal erythroid differentiation. Our findings also demonstrate that disruption of the miR-144/451-Myc crosstalk causes anemia, suggesting that miR-144/451 might be a potential therapeutic target in red cell diseases.

Methylation of a CGATA element inhibits binding and regulation by GATA-1

Alterations in DNA methylation occur during development, but the mechanisms by which they influence gene expression remain uncertain. There are few examples where modification of a single CpG dinucleotide directly affects transcription factor binding and regulation of a target gene in vivo. Here, we show that the erythroid transcription factor GATA-1 — that typically binds T/AGATA sites — can also recognise CGATA elements, but only if the CpG dinucleotide is unmethylated. We focus on a single CGATA site in the c-Kit gene which progressively becomes unmethylated during haematopoiesis. We observe that methylation attenuates GATA-1 binding and gene regulation in cell lines. In mice, converting the CGATA element to a TGATA site that cannot be methylated leads to accumulation of megakaryocyte-erythroid progenitors. Thus, the CpG dinucleotide is essential for normal erythropoiesis and this study illustrates how a single methylated CpG can directly affect transcription factor binding and cellular regulation.

While DNA methylation is thought to play a regulatory role, there are few examples where modification of a single CpG dinucleotide directly affects transcription factor binding. Here the authors show that methylation of a single CGATA element within the c-Kit gene inhibits binding and regulation by erythroid transcription factor GATA-1, both in cells and in mice, suggesting that methylation at this site plays an essential role in erythropoiesis.

KnockTF: a comprehensive human gene expression profile database with knockdown/knockout of transcription factors

Transcription factors (TFs) and their target genes have important functions in human diseases and biological processes. Gene expression profile analysis before and after knockdown or knockout is one of the most important strategies for obtaining target genes of TFs and exploring TF functions. Human gene expression profile datasets with TF knockdown and knockout are accumulating rapidly. Based on the urgent need to comprehensively and effectively collect and process these data, we developed KnockTF (http://www.licpathway.net/KnockTF/index.html), a comprehensive human gene expression profile database of TF knockdown and knockout. KnockTF provides a number of resources for human gene expression profile datasets associated with TF knockdown and knockout and annotates TFs and their target genes in a tissue/cell type-specific manner. The current version of KnockTF has 570 manually curated RNA-seq and microarray datasets associated with 308 TFs disrupted by different knockdown and knockout techniques and across multiple tissue/cell types. KnockTF collects upstream pathway information of TFs and functional annotation results of downstream target genes. It provides details about TFs binding to promoters, super-enhancers and typical enhancers of target genes. KnockTF constructs a TF-differentially expressed gene network and performs network analyses for genes of interest. KnockTF will help elucidate TF-related functions and potential biological effects.

GATA Factor-Mediated Gene Regulation in Human Erythropoiesis

Erythroid commitment and differentiation are regulated by the coordinated action of a host of transcription factors, including GATA2 and GATA1. Here, we explored GATA-mediated transcriptional regulation through the integrative analysis of gene expression, chromatin modifications, and GATA factors' binding in human multipotent hematopoietic stem/progenitor cells, early erythroid progenitors, and late precursors. A progressive loss of H3K27 acetylation and a diminished usage of active enhancers and super-enhancers were observed during erythroid commitment and differentiation. GATA factors mediate transcriptional changes through a stage-specific interplay with regulatory elements: GATA1 binds different sets of regulatory elements in erythroid progenitors and precursors and controls the transcription of distinct genes during commitment and differentiation. Importantly, our results highlight a pivotal role of promoters in determining the transcriptional program activated upon erythroid differentiation. Finally, we demonstrated that GATA1 binding to a stage-specific super-enhancer sustains the expression of the KIT receptor in human erythroid progenitors.

•

GATA2/1 binding to regulatory regions and transcriptional changes during erythropoiesis

•

GATA1 sustains KIT expression in human erythroid progenitors

Multiplexed capture of spatial configuration and temporal dynamics of locus-specific 3D chromatin by biotinylated dCas9

The spatiotemporal control of 3D genome is fundamental for gene regulation, yet it remains challenging to profile high-resolution chromatin structure at cis-regulatory elements (CREs). Using C-terminally biotinylated dCas9, endogenous biotin ligases, and pooled sgRNAs, we describe the dCas9-based CAPTURE method for multiplexed analysis of locus-specific chromatin interactions. The redesigned system allows for quantitative analysis of the spatial configuration of a few to hundreds of enhancers or promoters in a single experiment, enabling comparisons across CREs within and between gene clusters. Multiplexed analyses of the spatiotemporal configuration of erythroid super-enhancers and promoter-centric interactions reveal organizational principles of genome structure and function.

An integrative view of the regulatory and transcriptional landscapes in mouse hematopoiesis

Thousands of epigenomic data sets have been generated in the past decade, but it is difficult for researchers to effectively use all the data relevant to their projects. Systematic integrative analysis can help meet this need, and the VISION project was established for validated systematic integration of epigenomic data in hematopoiesis. Here, we systematically integrated extensive data recording epigenetic features and transcriptomes from many sources, including individual laboratories and consortia, to produce a comprehensive view of the regulatory landscape of differentiating hematopoietic cell types in mouse. By using IDEAS as our integrative and discriminative epigenome annotation system, we identified and assigned epigenetic states simultaneously along chromosomes and across cell types, precisely and comprehensively. Combining nuclease accessibility and epigenetic states produced a set of more than 200,000 candidate cis-regulatory elements (cCREs) that efficiently capture enhancers and promoters. The transitions in epigenetic states of these cCREs across cell types provided insights into mechanisms of regulation, including decreases in numbers of active cCREs during differentiation of most lineages, transitions from poised to active or inactive states, and shifts in nuclease accessibility of CTCF-bound elements. Regression modeling of epigenetic states at cCREs and gene expression produced a versatile resource to improve selection of cCREs potentially regulating target genes. These resources are available from our VISION website to aid research in genomics and hematopoiesis.

Chromatin occupancy and epigenetic analysis reveal new insights into the function of the GATA1 N terminus in erythropoiesis

Mutations in GATA1, which lead to expression of the GATA1s isoform that lacks the GATA1 N terminus, are seen in patients with Diamond-Blackfan anemia (DBA). In our efforts to better understand the connection between GATA1s and DBA, we comprehensively studied erythropoiesis in Gata1s mice. Defects in yolks sac and fetal liver hematopoiesis included impaired terminal maturation and reduced numbers of erythroid progenitors. RNA-sequencing revealed that both erythroid and megakaryocytic gene expression patterns were altered by the loss of the N terminus, including aberrant upregulation of Gata2 and Runx1. Dysregulation of global H3K27 methylation was found in the erythroid progenitors upon loss of N terminus of GATA1. Chromatin-binding assays revealed that, despite similar occupancy of GATA1 and GATA1s, there was a striking reduction of H3K27me3 at regulatory elements of the Gata2 and Runx1 genes. Consistent with the observation that overexpression of GATA2 has been reported to impair erythropoiesis, we found that haploinsufficiency of Gata2 rescued the erythroid defects of Gata1s fetuses. Together, our integrated genomic analysis of transcriptomic and epigenetic signatures reveals that, Gata1 mice provide novel insights into the role of the N terminus of GATA1 in transcriptional regulation and red blood cell maturation which may potentially be useful for DBA patients.

Systematic integration of GATA transcription factors and epigenomes via IDEAS paints the regulatory landscape of hematopoietic cells

Members of the GATA family of transcription factors play key roles in the differentiation of specific cell lineages by regulating the expression of target genes. Three GATA factors play distinct roles in hematopoietic differentiation. In order to better understand how these GATA factors function to regulate genes throughout the genome, we are studying the epigenomic and transcriptional landscapes of hematopoietic cells in a model-driven, integrative fashion. We have formed the collaborative multi-lab VISION project to conduct ValIdated Systematic IntegratiON of epigenomic data in mouse and human hematopoiesis. The epigenomic data included nuclease accessibility in chromatin, CTCF occupancy, and histone H3 modifications for 20 cell types covering hematopoietic stem cells, multilineage progenitor cells, and mature cells across the blood cell lineages of mouse. The analysis used the Integrative and Discriminative Epigenome Annotation System (IDEAS), which learns all common combinations of features (epigenetic states) simultaneously in two dimensions-along chromosomes and across cell types. The result is a segmentation that effectively paints the regulatory landscape in readily interpretable views, revealing constitutively active or silent loci as well as the loci specifically induced or repressed in each stage and lineage. Nuclease accessible DNA segments in active chromatin states were designated candidate cis-regulatory elements in each cell type, providing one of the most comprehensive registries of candidate hematopoietic regulatory elements to date. Applications of VISION resources are illustrated for the regulation of genes encoding GATA1, GATA2, GATA3, and Ikaros. VISION resources are freely available from our website http://usevision.org.

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.

GATA1 mutations in red cell disorders

GATA1 is an essential regulator of erythroid cell gene expression and maturation. In its absence, erythroid progenitors are arrested in differentiation and undergo apoptosis. Much has been learned about GATA1 function through animal models, which include genetic knockouts as well as ones with decreased levels of expression. However, even greater insights have come from the finding that a number of rare red cell disorders, including Diamond-Blackfan anemia, are associated with GATA1 mutations. These mutations affect the amino-terminal zinc finger (N-ZF) and the amino-terminus of the protein, and in both cases can alter the DNA-binding activity, which is primarily conferred by the third functional domain, the carboxyl-terminal zinc finger (C-ZF). Here we discuss the role of GATA1 in erythropoiesis with an emphasis on the mutations found in human patients with red cell disorders.

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.

Hit and Run Transcriptional Repressors Are Difficult to Catch in the Act

Transcriptional silencing may not necessarily depend on the continuous residence of a sequence-specific repressor at a control element and may act via a "hit and run" mechanism. Due to limitations in assays that detect transcription factor (TF) binding, such as chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq), this phenomenon may be challenging to detect and therefore its prevalence may be underappreciated. To explore this possibility, erythroid gene promoters that are regulated directly by GATA1 in an inducible system are analyzed. It is found that many regulated genes are bound immediately after induction of GATA1 but the residency of GATA1 decreases over time, particularly at repressed genes. Furthermore, it is shown that the repressive mark H3K27me3 is seldom associated with bound repressors, whereas, in contrast, the active (H3K4me3) histone mark is overwhelmingly associated with TF binding. It is hypothesized that during cellular differentiation and development, certain genes are silenced by repressive TFs that subsequently vacate the region. Catching such repressor TFs in the act of silencing via assays such as ChIP-seq is thus a temporally challenging prospect. The use of inducible systems, epitope tags, and alternative techniques may provide opportunities for detecting elusive "hit and run" transcriptional silencing. Also see the video abstract here https://youtu.be/vgrsoP_HF3g.

Neutrophil GM-CSF signaling in inflammatory bowel disease patients is influenced by non-coding genetic variants

Neutrophil dysfunction and GM-CSF auto-antibodies are observed in pediatric and adult patients with Crohn’s disease (CD). We associated damaging coding variants with low GM-CSF induced STAT5 stimulation index (GMSI) in pediatric CD patients and implicated variation of neutrophil GM-CSF signaling in cell function and disease complications. Because many CD patients with low GMSI do not carry damaging coding mutations, we sought to test the hypothesis that non-coding variants contribute to this phenotype. We enrolled, performed whole genome sequencing, and measured the GMSI in 77 CD and ulcerative colitis (UC) patients (24 low and 53 normal GMSI). We identified 4 non-coding variants (rs3808851, rs10974787, rs10974788 and rs10974789) in RCL1 significantly associated with variation of GMSI level (p < 0.011). They were validated in two independent cohorts with: RNAseq data (n = 50) and blood eQTL dataset (n = 31,684). These variants are in LD and affect expression of JAK2 (p 0.005 to 0.013), RCL1 (p 8.17E-13 to 2.98E-11) and AK3 (p 2.00E-68 to 3.03E-55) genes. Additionally, they influence proteins involved in differentiation of gut epithelium, inflammation, and immune system regulation. In summary, our study outlines the contribution of non-coding variants in neutrophil GM-CSF signaling and the potential importance of RCL1 and AK3 in CD pathogenesis.

Insights from resolving protein-DNA interactions at near base-pair resolution

One of the central goals in molecular biology is to understand how cell-type-specific expression patterns arise through selective recruitment of RNA polymerase II (Pol II) to a subset of gene promoters. Pol II needs to be recruited to a precise genomic position at the proper time to produce messenger RNA from a DNA template. Ostensibly, transcription is a relatively simple cellular process; yet, experimentally measuring and then understanding the combinatorial possibilities of transcriptional regulators remain a daunting task. Since its introduction in 1985, chromatin immunoprecipitation (ChIP) has remained a key tool for investigating protein-DNA contacts in vivo. Over 30 years of intensive research using ChIP have provided numerous insights into mechanisms of gene regulation. As functional genomic technologies improve, they present new opportunities to address key biological questions. ChIP-exo is a refined version of ChIP-seq that significantly reduces background signal, while providing near base-pair mapping resolution for protein-DNA interactions. This review discusses the evolution of the ChIP assay over the years; the methodological differences between ChIP-seq, ChIP-exo and ChIP-nexus; and highlight new insights into epigenetic and transcriptional mechanisms that were uniquely enabled with the near base-pair resolution of ChIP-exo.

Transcription-dependent spreading of the Dal80 yeast GATA factor across the body of highly expressed genes

GATA transcription factors are highly conserved among eukaryotes and play roles in transcription of genes implicated in cancer progression and hematopoiesis. However, although their consensus binding sites have been well defined in vitro, the in vivo selectivity for recognition by GATA factors remains poorly characterized. Using ChIP-Seq, we identified the Dal80 GATA factor targets in yeast. Our data reveal Dal80 binding to a large set of promoters, sometimes independently of GATA sites, correlating with nitrogen- and/or Dal80-sensitive gene expression. Strikingly, Dal80 was also detected across the body of promoter-bound genes, correlating with high expression. Mechanistic single-gene experiments showed that Dal80 spreading across gene bodies requires active transcription. Consistently, Dal80 co-immunoprecipitated with the initiating and post-initiation forms of RNA Polymerase II. Our work suggests that GATA factors could play dual, synergistic roles during transcription initiation and post-initiation steps, promoting efficient remodeling of the gene expression program in response to environmental changes.

GATA transcription factors are highly conserved among eukaryotes and play key roles in cancer progression and hematopoiesis. In budding yeast, four GATA transcription factors are involved in the response to the quality of nitrogen supply. Here, we have determined the whole genome binding profile of the Dal80 GATA factor, and revealed that it also associates with the body of promoter-bound genes. The observation that intragenic spreading correlates with high expression levels and exquisite Dal80 sensitivity suggests that GATA factors could play other, unexpected roles at post-initiation stages in eukaryotes.

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.

A natural regulatory mutation in the proximal promoter elevates fetal globin expression by creating a de novo GATA1 site

β-hemoglobinopathies, such as sickle cell disease and β-thalassemia, result from mutations in the adult β-globin gene. Reactivating the developmentally silenced fetal γ-globin gene elevates fetal hemoglobin levels and ameliorates symptoms of β-hemoglobinopathies. The continued expression of fetal γ-globin into adulthood occurs naturally in a genetic condition termed hereditary persistence of fetal hemoglobin (HPFH). Point mutations in the fetal γ-globin proximal promoter can cause HPFH. The -113A>G HPFH mutation falls within the -115 cluster of HPFH mutations, a binding site for the fetal globin repressor BCL11A. We demonstrate that the -113A>G HPFH mutation, unlike other mutations in the cluster, does not disrupt BCL11A binding but rather creates a de novo binding site for the transcriptional activator GATA1. Introduction of the -113A>G HPFH mutation into erythroid cells using the clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated protein 9 (Cas9) system increases GATA1 binding and elevates fetal globin levels. These results reveal the mechanism by which the -113A>G HPFH mutation elevates fetal globin and demonstrate the sensitivity of the fetal globin promoter to point mutations that often disrupt repressor binding sites but here create a de novo site for an erythroid activator.

Genomics and transcriptomics of megakaryocytes and platelets: Implications for health and disease

The field of megakaryocyte and platelet biology has been transformed with the implementation of high throughput sequencing. The use of modern sequencing technologies has led to the discovery of causative mutations in congenital platelet disorders and has been a useful tool in uncovering many other mechanisms of altered platelet formation and function. Although the understanding of the presence of RNA in platelets is relatively novel, mRNA and miRNA expression profiles are being shown to play an increasingly important role in megakaryopoiesis and platelet function in normal physiology as well as in disease states. Understanding the genetic perturbations underlying platelet dysfunction provides insight into normal megakaryopoiesis and thrombopoiesis, as well as guiding the development of novel therapeutics.

PI3k/AKT signaling pathway: Erythropoiesis and beyond

Erythropoiesis is a multi-step process that involves the differentiation of hematopoietic stem cells into mature red blood cells (RBCs). This process is regulated by several signaling pathways, transcription factors and microRNAs (miRNAs). Many studies have shown that dysregulation of this process can lead to hematologic disorders. PI3K/AKT is one of the most important pathways that control many cellular processes including, cell division, autophagy, survival, and differentiation. In this review, we focus on the role of PI3K/AKT pathway in erythropoiesis and discuss the function of some of the most important genes, transcription factors, and miRNAs that regulate different stages of erythropoiesis which play roles in differentiation and maturation of RBCs, prevention of apoptosis, and autophagy induction. Understanding the role of the PI3K pathway in erythropoiesis may provide new insights into diagnosing erythrocyte disorders.

DiNAMO: highly sensitive DNA motif discovery in high-throughput sequencing data

Background

Discovering over-represented approximate motifs in DNA sequences is an essential part of bioinformatics. This topic has been studied extensively because of the increasing number of potential applications. However, it remains a difficult challenge, especially with the huge quantity of data generated by high throughput sequencing technologies. To overcome this problem, existing tools use greedy algorithms and probabilistic approaches to find motifs in reasonable time. Nevertheless these approaches lack sensitivity and have difficulties coping with rare and subtle motifs.

Results

We developed DiNAMO (for DNA MOtif), a new software based on an exhaustive and efficient algorithm for IUPAC motif discovery. We evaluated DiNAMO on synthetic and real datasets with two different applications, namely ChIP-seq peaks and Systematic Sequencing Error analysis. DiNAMO proves to compare favorably with other existing methods and is robust to noise.

Conclusions

We shown that DiNAMO software can serve as a tool to search for degenerate motifs in an exact manner using IUPAC models. DiNAMO can be used in scanning mode with sliding windows or in fixed position mode, which makes it suitable for numerous potential applications.

Availability

https://github.com/bonsai-team/DiNAMO.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2215-1) contains supplementary material, which is available to authorized users.

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.

GATA1 insufficiencies in primary myelofibrosis and other hematopoietic disorders: consequences for therapy

INTRODUCTION

GATA1, the founding member of a family of transcription factors, plays important roles in the development of hematopoietic cells of several lineages. Although loss of GATA1 has been known to impair hematopoiesis in animal models for nearly 25 years, the link between GATA1 defects and human blood diseases has only recently been realized. Areas covered: Here the current understanding of the functions of GATA1 in normal hematopoiesis and how it is altered in disease is reviewed. GATA1 is indispensable mainly for erythroid and megakaryocyte differentiation. In erythroid cells, GATA1 regulates early stages of differentiation, and its deficiency results in apoptosis. In megakaryocytes, GATA1 controls terminal maturation and its deficiency induces proliferation. GATA1 alterations are often found in diseases involving these two lineages, such as congenital erythroid and/or megakaryocyte deficiencies, including Diamond Blackfan Anemia (DBA), and acquired neoplasms, such as acute megakaryocytic leukemia (AMKL) and the myeloproliferative neoplasms (MPNs). Expert commentary: Since the first discovery of GATA1 mutations in AMKL, the number of diseases that are associated with impaired GATA1 function has increased to include DBA and MPNs. With respect to the latter, we are only just now appreciating the link between enhanced JAK/STAT signaling, GATA1 deficiency and disease pathogenesis.

Exploiting genetic variation to uncover rules of transcription factor binding and chromatin accessibility

Single-nucleotide variants that underlie phenotypic variation can affect chromatin occupancy of transcription factors (TFs). To delineate determinants of in vivo TF binding and chromatin accessibility, we introduce an approach that compares ChIP-seq and DNase-seq data sets from genetically divergent murine erythroid cell lines. The impact of discriminatory single-nucleotide variants on TF ChIP signal enables definition at single base resolution of in vivo binding characteristics of nuclear factors GATA1, TAL1, and CTCF. We further develop a facile complementary approach to more deeply test the requirements of critical nucleotide positions for TF binding by combining CRISPR-Cas9-mediated mutagenesis with ChIP and targeted deep sequencing. Finally, we extend our analytical pipeline to identify nearby contextual DNA elements that modulate chromatin binding by these three TFs, and to define sequences that impact kb-scale chromatin accessibility. Combined, our approaches reveal insights into the genetic basis of TF occupancy and their interplay with chromatin features.