Dynamics of the epigenetic landscape during erythroid differentiation after GATA1 restoration

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.

Identification of mammalian transcription factors that bind to inaccessible chromatin

Transcription factors (TFs) are proteins that affect gene expression by binding to regulatory regions of DNA in a sequence specific manner. The binding of TFs to DNA is controlled by many factors, including the DNA sequence, concentration of TF, chromatin accessibility and co-factors. Here, we systematically investigated the binding mechanism of hundreds of TFs by analysing ChIP-seq data with our explainable statistical model, ChIPanalyser. This tool uses as inputs the DNA sequence binding motif; the capacity to distinguish between strong and weak binding sites; the concentration of TF; and chromatin accessibility. We found that approximately one third of TFs are predicted to bind the genome in a DNA accessibility independent fashion, which includes TFs that can open the chromatin, their co-factors and TFs with similar motifs. Our model predicted this to be the case when the TF binds to its strongest binding regions in the genome, and only a small number of TFs have the capacity to bind dense chromatin at their weakest binding regions, such as CTCF, USF2 and CEBPB. Our study demonstrated that the binding of hundreds of human and mouse TFs is predicted by ChIPanalyser with high accuracy and showed that many TFs can bind dense chromatin.

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.

Gene silencing dynamics are modulated by transiently active regulatory elements

Transcriptional enhancers have been extensively characterized, but cis-regulatory elements involved in acute gene repression have received less attention. Transcription factor GATA1 promotes erythroid differentiation by activating and repressing distinct gene sets. Here, we study the mechanism by which GATA1 silences the proliferative gene Kit during murine erythroid cell maturation and define stages from initial loss of activation to heterochromatinization. We find that GATA1 inactivates a potent upstream enhancer but concomitantly creates a discrete intronic regulatory region marked by H3K27ac, short noncoding RNAs, and de novo chromatin looping. This enhancer-like element forms transiently and serves to delay Kit silencing. The element is ultimately erased via the FOG1/NuRD deacetylase complex, as revealed by the study of a disease-associated GATA1 variant. Hence, regulatory sites can be self-limiting by dynamic co-factor usage. Genome-wide analyses across cell types and species uncover transiently active elements at numerous genes during repression, suggesting that modulation of silencing kinetics is widespread.

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.

Epigenetic and Transcriptional Control of Erythropoiesis

Erythropoiesis is a process of enormous magnitude, with the average person generating two to three million red cells every second. Erythroid progenitors start as large cells with large nuclei, and over the course of three to four cell divisions they undergo a dramatic decrease in cell size accompanied by profound nuclear condensation, which culminates in enucleation. As maturing erythroblasts are undergoing these dramatic phenotypic changes, they accumulate hemoglobin and express high levels of other erythroid-specific genes, while silencing much of the non-erythroid transcriptome. These phenotypic and gene expression changes are associated with distinct changes in the chromatin landscape, and require close coordination between transcription factors and epigenetic regulators, as well as precise regulation of RNA polymerase II activity. Disruption of these processes are associated with inherited anemias and myelodysplastic syndromes. Here, we review the epigenetic mechanisms that govern terminal erythroid maturation, and their role in human disease.

Regulation of RNA polymerase II activity is essential for terminal erythroid maturation

The terminal maturation of human erythroblasts requires significant changes in gene expression in the context of dramatic nuclear condensation. Defects in this process are associated with inherited anemias and myelodysplastic syndromes. The progressively dense appearance of the condensing nucleus in maturing erythroblasts led to the assumption that heterochromatin accumulation underlies this process, but despite extensive study, the precise mechanisms underlying this essential biologic process remain elusive. To delineate the epigenetic changes associated with the terminal maturation of human erythroblasts, we performed mass spectrometry of histone posttranslational modifications combined with chromatin immunoprecipitation coupled with high-throughput sequencing, Assay for Transposase Accessible Chromatin, and RNA sequencing. Our studies revealed that the terminal maturation of human erythroblasts is associated with a dramatic decline in histone marks associated with active transcription elongation, without accumulation of heterochromatin. Chromatin structure and gene expression were instead correlated with dynamic changes in occupancy of elongation competent RNA polymerase II, suggesting that terminal erythroid maturation is controlled largely at the level of transcription. We further demonstrate that RNA polymerase II "pausing" is highly correlated with transcriptional repression, with elongation competent RNA polymerase II becoming a scare resource in late-stage erythroblasts, allocated to erythroid-specific genes. Functional studies confirmed an essential role for maturation stage-specific regulation of RNA polymerase II activity during erythroid maturation and demonstrate a critical role for HEXIM1 in the regulation of gene expression and RNA polymerase II activity in maturing erythroblasts. Taken together, our findings reveal important insights into the mechanisms that regulate terminal erythroid maturation and provide a novel paradigm for understanding normal and perturbed erythropoiesis.

AIAP: A Quality Control and Integrative Analysis Package to Improve ATAC-seq Data Analysis

Assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) is a technique widely used to investigate genome-wide chromatin accessibility. The recently published Omni-ATAC-seq protocol substantially improves the signal/noise ratio and reduces the input cell number. High-quality data are critical to ensure accurate analysis. Several tools have been developed for assessing sequencing quality and insertion size distribution for ATAC-seq data; however, key quality control (QC) metrics have not yet been established to accurately determine the quality of ATAC-seq data. Here, we optimized the analysis strategy for ATAC-seq and defined a series of QC metrics for ATAC-seq data, including reads under peak ratio (RUPr), background (BG), promoter enrichment (ProEn), subsampling enrichment (SubEn), and other measurements. We incorporated these QC tests into our recently developed ATAC-seq Integrative Analysis Package (AIAP) to provide a complete ATAC-seq analysis system, including quality assurance, improved peak calling, and downstream differential analysis. We demonstrated a significant improvement of sensitivity (20%–60%) in both peak calling and differential analysis by processing paired-end ATAC-seq datasets using AIAP. AIAP is compiled into Docker/Singularity, and it can be executed by one command line to generate a comprehensive QC report. We used ENCODE ATAC-seq data to benchmark and generate QC recommendations, and developed qATACViewer for the user-friendly interaction with the QC report. The software, source code, and documentation of AIAP are freely available at https://github.com/Zhang-lab/ATAC-seq_QC_analysis.

Coordinated changes in gene expression kinetics underlie both mouse and human erythroid maturation

BACKGROUND

Single-cell technologies are transforming biomedical research, including the recent demonstration that unspliced pre-mRNA present in single-cell RNA-Seq permits prediction of future expression states. Here we apply this RNA velocity concept to an extended timecourse dataset covering mouse gastrulation and early organogenesis.

RESULTS

Intriguingly, RNA velocity correctly identifies epiblast cells as the starting point, but several trajectory predictions at later stages are inconsistent with both real-time ordering and existing knowledge. The most striking discrepancy concerns red blood cell maturation, with velocity-inferred trajectories opposing the true differentiation path. Investigating the underlying causes reveals a group of genes with a coordinated step-change in transcription, thus violating the assumptions behind current velocity analysis suites, which do not accommodate time-dependent changes in expression dynamics. Using scRNA-Seq analysis of chimeric mouse embryos lacking the major erythroid regulator Gata1, we show that genes with the step-changes in expression dynamics during erythroid differentiation fail to be upregulated in the mutant cells, thus underscoring the coordination of modulating transcription rate along a differentiation trajectory. In addition to the expected block in erythroid maturation, the Gata1-chimera dataset reveals induction of PU.1 and expansion of megakaryocyte progenitors. Finally, we show that erythropoiesis in human fetal liver is similarly characterized by a coordinated step-change in gene expression.

CONCLUSIONS

By identifying a limitation of the current velocity framework coupled with in vivo analysis of mutant cells, we reveal a coordinated step-change in gene expression kinetics during erythropoiesis, with likely implications for many other differentiation processes.

RNA-regulatory exosome complex confers cellular survival to promote erythropoiesis

Cellular differentiation requires vast remodeling of transcriptomes, and therefore machinery mediating remodeling controls differentiation. Relative to transcriptional mechanisms governing differentiation, post-transcriptional processes are less well understood. As an important post-transcriptional determinant of transcriptomes, the RNA exosome complex (EC) mediates processing and/or degradation of select RNAs. During erythropoiesis, the erythroid transcription factor GATA1 represses EC subunit genes. Depleting EC structural subunits prior to GATA1-mediated repression is deleterious to erythroid progenitor cells. To assess the importance of the EC catalytic subunits Dis3 and Exosc10 in this dynamic process, we asked if these subunits function non-redundantly to control erythropoiesis. Dis3 or Exosc10 depletion in primary murine hematopoietic progenitor cells reduced erythroid progenitors and their progeny, while sparing myeloid cells. Dis3 loss severely compromised erythroid progenitor and erythroblast survival, rendered erythroblasts hypersensitive to apoptosis-inducing stimuli and induced γ-H2AX, indicative of DNA double-stranded breaks. Dis3 loss-of-function phenotypes were more severe than those caused by Exosc10 depletion. We innovated a genetic rescue system to compare human Dis3 with multiple myeloma-associated Dis3 mutants S447R and R750K, and only wild type Dis3 was competent to rescue progenitors. Thus, Dis3 establishes a disease mutation-sensitive, cell type-specific survival mechanism to enable a differentiation program.

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.

GATA factor-regulated solute carrier ensemble reveals a nucleoside transporter-dependent differentiation mechanism

Developmental-regulatory networks often include large gene families encoding mechanistically-related proteins like G-protein-coupled receptors, zinc finger transcription factors and solute carrier (SLC) transporters. In principle, a common mechanism may confer expression of multiple members integral to a developmental process, or diverse mechanisms may be deployed. Using genetic complementation and enhancer-mutant systems, we analyzed the 456 member SLC family that establishes the small molecule constitution of cells. This analysis identified SLC gene cohorts regulated by GATA1 and/or GATA2 during erythroid differentiation. As >50 SLC genes shared GATA factor regulation, a common mechanism established multiple members of this family. These genes included Slc29a1 encoding an equilibrative nucleoside transporter (Slc29a1/ENT1) that utilizes adenosine as a preferred substrate. Slc29a1 promoted erythroblast survival and differentiation ex vivo. Targeted ablation of murine Slc29a1 in erythroblasts attenuated erythropoiesis and erythrocyte regeneration in response to acute anemia. Our results reveal a GATA factor-regulated SLC ensemble, with a nucleoside transporter component that promotes erythropoiesis and prevents anemia, and establish a mechanistic link between GATA factor and adenosine mechanisms. We propose that integration of the GATA factor-adenosine circuit with other components of the GATA factor-regulated SLC ensemble establishes the small molecule repertoire required for progenitor cells to efficiently generate erythrocytes.

GATA transcription factors endow blood stem and progenitor cells with activities to produce progeny that transport oxygen to protect cells and tissues, evade pathogens and control physiological processes. GATA factors regulate hundreds of genes, and the actions of these genes mediate important biological functions. While the genes have been documented, many questions remain regarding how the “network” components mediate biological functions. The networks include members of large gene families, and the relationships between the regulation and function of individual family members is not well understood. Analyzing datasets from genetic complementation and enhancer mutant systems revealed that GATA factors regulate an ensemble of membrane transporters termed solute carrier proteins (SLCs), which dictate the small molecule composition of cells. Genetic analyses with Slc29a1, which transports adenosine, revealed its function to promote erythrocyte development, and Slc29a1 attenuated anemia in a mouse model. This study revealed the importance of SLC transporters in GATA factor networks. We propose that the GATA factor-adenosine circuit integrates with other SLCs to establish/maintain the small molecule constitution of progenitor cells as a new mechanism to control blood cell development.

Advances in Nutritional Epigenetics-A Fresh Perspective for an Old Idea. Lessons Learned, Limitations, and Future Directions

Nutritional epigenetics is a rapidly expanding field of research, and the natural modulation of the genome is a non-invasive, sustainable, and personalized alternative to gene-editing for chronic disease management. Genetic differences and epigenetic inflexibility resulting in abnormal gene expression, differential or aberrant methylation patterns account for the vast majority of diseases. The expanding understanding of biological evolution and the environmental influence on epigenetics and natural selection requires relearning of once thought to be well-understood concepts. This research explores the potential for natural modulation by the less understood epigenetic modifications such as ubiquitination, nitrosylation, glycosylation, phosphorylation, and serotonylation concluding that the under-appreciated acetylation and mitochondrial dependant downstream epigenetic post-translational modifications may be the pinnacle of the epigenomic hierarchy, essential for optimal health, including sustainable cellular energy production. With an emphasis on lessons learned, this conceptional exploration provides a fresh perspective on methylation, demonstrating how increases in environmental methane drive an evolutionary down regulation of endogenous methyl groups synthesis and demonstrates how epigenetic mechanisms are cell-specific, making supplementation with methyl cofactors throughout differentiation unpredictable. Interference with the epigenomic hierarchy may result in epigenetic inflexibility, symptom relief and disease concomitantly and may be responsible for the increased incidence of neurological disease such as autism spectrum disorder.

Genetic screens reveal CCDC115 as a modulator of erythroid iron and heme trafficking

Transferrin-bound iron (TBI), the physiological circulating iron form, is acquired by cells through the transferrin receptor (TfR1) by endocytosis. In erythroid cells, most of the acquired iron is incorporated into heme in the mitochondria. Cellular trafficking of heme is indispensable for erythropoiesis and many other essential biological processes. Comprehensive elucidation of molecular pathways governing and regulating cellular iron acquisition and heme trafficking is required to better understand physiological and pathological processes affecting erythropoiesis. Here, we report the first genome-wide clustered regularly interspaced short palindromic repeats (CRISPR) screens in human erythroid cells to identify determinants of iron and heme uptake, as well as heme-mediated erythroid differentiation. We identified several candidate modulators of TBI acquisition including TfR1, indicating that our approach effectively revealed players mechanistically relevant to the process. Interestingly, components of the endocytic pathway were also revealed as potential determinants of transferrin acquisition. We deciphered a role for the vacuolar-type H+ - ATPase (V- ATPase) assembly factor coiled-coil domain containing 115 (CCDC115) in TBI uptake and validated this role in CCDC115 deficient K562 cells. Our screen in hemin-treated cells revealed perturbations leading to cellular adaptation to heme, including those corresponding to trafficking mechanisms and transcription factors potentiating erythroid differentiation. Pathway analysis indicated that endocytosis and vesicle acidification are key processes for heme trafficking in erythroid precursors. Furthermore, we provided evidence that CCDC115, which we identified as required for TBI uptake, is also involved in cellular heme distribution. This work demonstrates a previously unappreciated common intersection in trafficking of transferrin iron and heme in the endocytic pathway of erythroid cells.

Nuclear interacting SET domain protein 1 inactivation impairs GATA1-regulated erythroid differentiation and causes erythroleukemia

The nuclear receptor binding SET domain protein 1 (NSD1) is recurrently mutated in human cancers including acute leukemia. We show that NSD1 knockdown alters erythroid clonogenic growth of human CD34⁺ hematopoietic cells. Ablation of Nsd1 in the hematopoietic system of mice induces a transplantable erythroleukemia. In vitro differentiation of Nsd1^−/− erythroblasts is majorly impaired despite abundant expression of GATA1, the transcriptional master regulator of erythropoiesis, and associated with an impaired activation of GATA1-induced targets. Retroviral expression of wildtype NSD1, but not a catalytically-inactive NSD1^N1918Q SET-domain mutant induces terminal maturation of Nsd1^−/− erythroblasts. Despite similar GATA1 protein levels, exogenous NSD1 but not NSD^N1918Q significantly increases the occupancy of GATA1 at target genes and their expression. Notably, exogenous NSD1 reduces the association of GATA1 with the co-repressor SKI, and knockdown of SKI induces differentiation of Nsd1^−/− erythroblasts. Collectively, we identify the NSD1 methyltransferase as a regulator of GATA1-controlled erythroid differentiation and leukemogenesis.

Loss of function mutations of NSD1 occur in blood cancers. Here, the authors report that NSD1 loss blocks erythroid differentiation which leads to an erythroleukemia-like disease in mice by impairing GATA1-induced target gene activation.

A Bayesian data fusion based approach for learning genome-wide transcriptional regulatory networks

BACKGROUND

Reverse engineering of transcriptional regulatory networks (TRN) from genomics data has always represented a computational challenge in System Biology. The major issue is modeling the complex crosstalk among transcription factors (TFs) and their target genes, with a method able to handle both the high number of interacting variables and the noise in the available heterogeneous experimental sources of information.

RESULTS

In this work, we propose a data fusion approach that exploits the integration of complementary omics-data as prior knowledge within a Bayesian framework, in order to learn and model large-scale transcriptional networks. We develop a hybrid structure-learning algorithm able to jointly combine TFs ChIP-Sequencing data and gene expression compendia to reconstruct TRNs in a genome-wide perspective. Applying our method to high-throughput data, we verified its ability to deal with the complexity of a genomic TRN, providing a snapshot of the synergistic TFs regulatory activity. Given the noisy nature of data-driven prior knowledge, which potentially contains incorrect information, we also tested the method's robustness to false priors on a benchmark dataset, comparing the proposed approach to other regulatory network reconstruction algorithms. We demonstrated the effectiveness of our framework by evaluating structural commonalities of our learned genomic network with other existing networks inferred by different DNA binding information-based methods.

CONCLUSIONS

This Bayesian omics-data fusion based methodology allows to gain a genome-wide picture of the transcriptional interplay, helping to unravel key hierarchical transcriptional interactions, which could be subsequently investigated, and it represents a promising learning approach suitable for multi-layered genomic data integration, given its robustness to noisy sources and its tailored framework for handling high dimensional data.

Protection from DNA re-methylation by transcription factors in primordial germ cells and pre-implantation embryos can explain trans-generational epigenetic inheritance

Background

A growing body of evidence suggests that certain epiphenotypes can be passed across generations via both the male and female germlines of mammals. These observations have been difficult to explain owing to a global loss of the majority of known epigenetic marks present in parental chromosomes during primordial germ cell development and after fertilization.

Results

By integrating previously published BS-seq, DNase-seq, ATAC-seq, and RNA-seq data collected during multiple stages of primordial germ cell and pre-implantation development, we find that the methylation status of the majority of CpGs genome-wide is restored after global de-methylation, despite the fact that global CpG methylation drops to 10% in primordial germ cells and 20% in the inner cell mass of the blastocyst. We estimate the proportion of such CpGs with preserved methylation status to be 78%. Further, we find that CpGs at sites bound by transcription factors during the global re-methylation phases of germline and embryonic development remain hypomethylated across all developmental stages observed. On the other hand, CpGs at sites not bound by transcription factors during the global re-methylation phase have high methylation levels prior to global de-methylation, become de-methylated during global de-methylation, and then become re-methylated.

Conclusions

The results suggest that transcription factors can act as carriers of epigenetic information during germ cell and pre-implantation development by ensuring that the methylation status of CpGs is maintained. These findings provide the basis for a mechanistic description of trans-generational inheritance of epigenetic information in mammals.

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

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.

Cross Talk between One-Carbon Metabolism, Eph Signaling, and Histone Methylation Promotes Neural Stem Cell Differentiation

Metabolic pathways, once seen as a mere consequence of cell states, have emerged as active players in dictating different cellular events such as proliferation, self-renewal, and differentiation. Several studies have reported a role for folate-dependent one-carbon (1C) metabolism in stem cells; however, its exact mode of action and how it interacts with other cues are largely unknown. Here, we report a link between the Eph:ephrin cell-cell communication pathway and 1C metabolism in controlling neural stem cell differentiation. Transcriptional and functional analyses following ephrin stimulation revealed alterations in folate metabolism-related genes and enzymatic activity. In vitro and in vivo data indicate that Eph-B forward signaling alters the methylation state of H3K4 by regulating 1C metabolism and locks neural stem cell in a differentiation-ready state. Our study highlights a functional link between cell-cell communication, metabolism, and epigenomic remodeling in the control of stem cell self-renewal.

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.

Long noncoding RNA PCED1B-AS1 promotes erythroid differentiation coordinating with GATA1 and chromatin remodeling

Erythropoiesis is a complex and sophisticated multi-stage process regulated by a variety of factors, including the transcription factor GATA1 and non-coding RNA. GATA1 is regarded as an essential transcriptional regulator promoting transcription of erythroid-specific genes—such as long non-coding RNAs (lncRNA). Here, we comprehensively screened lncRNAs that were potentially regulated by GATA1 in erythroid cells. We identified a novel lncRNA—PCED1B-AS1—and verified its role in promoting erythroid differentiation of K562 erythroid cells. We also predicted a model in which PCED1B-AS1 participates in erythroid differentiation via dynamic chromatin remodeling involving GATA1. The relationship between lncRNA and chromatin in the process of erythroid differentiation remains to be revealed, and in our study we have carried out preliminary explorations.

EpiAlignment: alignment with both DNA sequence and epigenomic data

Comparative epigenomics, which subjects both epigenome and genome to interspecies comparison, has become a powerful approach to reveal regulatory features of the genome. Thus elucidated regulatory features surpass the information derived from comparison of genomic sequences alone. Here, we present EpiAlignment, a web-based tool to align genomic regions with both DNA sequence and epigenomic data. EpiAlignment takes DNA sequence and epigenomic profiles derived by ChIP-seq from two species as input data, and outputs the best semi-global alignments. These alignments are based on EpiAlignment scores, computed by a dynamic programming algorithm that accounts for both sequence alignment and epigenome similarity. For timely response, the EpiAlignment web server automatically initiates up to 140 computing threads depending on the size of user input data. For users’ convenience, we have pre-compiled the comparable human and mouse epigenome datasets in matched cell types and tissues from the Roadmap Epigenomics and ENCODE consortia. Users can either upload their own data or select pre-compiled datasets as inputs for EpiAlignment analyses. Results are presented in graphical and tabular formats where the entries can be interactively expanded to visualize additional features of these aligned regions. EpiAlignment is available at https://epialign.ucsd.edu/.

Interrogating Histone Acetylation and BRD4 as Mitotic Bookmarks of Transcription

Global changes in chromatin organization and the cessation of transcription during mitosis are thought to challenge the resumption of appropriate transcription patterns after mitosis. The acetyl-lysine binding protein BRD4 has been previously suggested to function as a transcriptional "bookmark" on mitotic chromatin. Here, genome-wide location analysis of BRD4 in erythroid cells, combined with data normalization and peak characterization approaches, reveals that BRD4 widely occupies mitotic chromatin. However, removal of BRD4 from mitotic chromatin does not impair post-mitotic activation of transcription. Additionally, histone mass spectrometry reveals global preservation of most posttranslational modifications (PTMs) during mitosis. In particular, H3K14ac, H3K27ac, H3K122ac, and H4K16ac widely mark mitotic chromatin, especially at lineage-specific genes, and predict BRD4 mitotic binding genome wide. Therefore, BRD4 is likely not a mitotic bookmark but only a "passenger." Instead, mitotic histone acetylation patterns may constitute the actual bookmarks that restore lineage-specific transcription patterns after mitosis.

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.

Integrative view on how erythropoietin signaling controls transcription patterns in erythroid cells

PURPOSE OF REVIEW

Erythropoietin (EPO) is necessary and sufficient to trigger dynamic transcriptional patterns that drive the differentiation of erythroid precursor cells into mature, enucleated red cells. Because the molecular cloning and Food and Drug Administration approval for the therapeutic use of EPO over 30 years ago, a detailed understanding of how EPO works has advanced substantially. Yet, the precise epigenetic and transcriptional mechanisms by which EPO signaling controls erythroid expression patterns remains poorly understood. This review focuses on the current state of erythroid biology in regards to EPO signaling from human genetics and functional genomics perspectives.

RECENT FINDINGS

The goal of this review is to provide an integrative view of the gene regulatory underpinnings for erythroid expression patterns that are dynamically shaped during erythroid differentiation. Here, we highlight vignettes connecting recent insights into a genome-wide association study linking an EPO mutation to anemia, a study linking EPO-signaling to signal transducer and activator of transcription 5 (STAT5) chromatin occupancy and enhancers, and studies that examine the molecular mechanisms driving topological chromatin organization in erythroid cells.

SUMMARY

The genetic, epigenetic, and gene regulatory mechanisms underlying how hormone signal transduction influences erythroid gene expression remains only partly understood. A detailed understanding of these molecular pathways and how they intersect with one another will provide the basis for novel strategies to treat anemia and potentially other hematological diseases. As new regulators and signal transducers of EPO-signaling continue to emerge, new clinically relevant targets may be identified that improve the specificity and effectiveness of EPO therapy.

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.

Whole-transcriptome splicing profiling of E7.5 mouse primary germ layers reveals frequent alternative promoter usage during mouse early embryogenesis

Alternative splicing (AS) and alternative promoter (AP) usage expand the repertories of mammalian transcriptome profiles and thus diversify gene functions. However, our knowledge about the extent and functions of AS and AP usage in mouse early embryogenesis remains elusive. Here, by performing whole-transcriptome splicing profiling with high-throughput next generation sequencing, we report that AS extensively occurs in embryonic day (E) 7.5 mouse primary germ layers, and may be involved in multiple developmental processes. In addition, numerous RNA splicing factors are differentially expressed and alternatively spliced across the three germ layers, implying the potential importance of AS machinery in shaping early embryogenesis. Notably, AP usage is remarkably frequent at this stage, accounting for more than one quarter (430/1,648) of the total significantly different AS events. Genes generating the 430 AP events participate in numerous biological processes, and include important regulators essential for mouse early embryogenesis, suggesting that AP usage is widely used and might be relevant to mouse germ layer specification. Our data underline the potential significance of AP usage in mouse gastrulation, providing a rich data source and opening another dimension for understanding the regulatory mechanisms of mammalian early development.

Summary: This study seeks to capture the alternative splicing landscape during mouse gastrulation, underlining the potential importance of alternative promoter usage in mammalian early embryogenesis.

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.

Direct targets of pSTAT5 signalling in erythropoiesis

Erythropoietin (EPO) acts through the dimeric erythropoietin receptor to stimulate proliferation, survival, differentiation and enucleation of erythroid progenitor cells. We undertook two complimentary approaches to find EPO-dependent pSTAT5 target genes in murine erythroid cells: RNA-seq of newly transcribed (4sU-labelled) RNA, and ChIP-seq for pSTAT5 30 minutes after EPO stimulation. We found 302 pSTAT5-occupied sites: ~15% of these reside in promoters while the rest reside within intronic enhancers or intergenic regions, some >100kb from the nearest TSS. The majority of pSTAT5 peaks contain a central palindromic GAS element, TTCYXRGAA. There was significant enrichment for GATA motifs and CACCC-box motifs within the neighbourhood of pSTAT5-bound peaks, and GATA1 and/or KLF1 co-occupancy at many sites. Using 4sU-RNA-seq we determined the EPO-induced transcriptome and validated differentially expressed genes using dynamic CAGE data and qRT-PCR. We identified known direct pSTAT5 target genes such as Bcl2l1, Pim1 and Cish, and many new targets likely to be involved in driving erythroid cell differentiation including those involved in mRNA splicing (Rbm25), epigenetic regulation (Suv420h2), and EpoR turnover (Clint1/EpsinR). Some of these new EpoR-JAK2-pSTAT5 target genes could be used as biomarkers for monitoring disease activity in polycythaemia vera, and for monitoring responses to JAK inhibitors.

A hyperactive transcriptional state marks genome reactivation at the mitosis-G1 transition

Hsiung et al. tracked Pol II occupancy genome-wide in mammalian cells progressing from mitosis through late G1. During the earliest rounds of transcription at the mitosis–G1 transition, ∼50% of active genes and distal enhancers exhibit a spike in transcription, exceeding levels observed later in G1 phase. The transcriptional spike occurs heterogeneously and propagates to cell-to-cell differences in mature mRNA expression.

During mitosis, RNA polymerase II (Pol II) and many transcription factors dissociate from chromatin, and transcription ceases globally. Transcription is known to restart in bulk by telophase, but whether de novo transcription at the mitosis–G1 transition is in any way distinct from later in interphase remains unknown. We tracked Pol II occupancy genome-wide in mammalian cells progressing from mitosis through late G1. Unexpectedly, during the earliest rounds of transcription at the mitosis–G1 transition, ∼50% of active genes and distal enhancers exhibit a spike in transcription, exceeding levels observed later in G1 phase. Enhancer–promoter chromatin contacts are depleted during mitosis and restored rapidly upon G1 entry but do not spike. Of the chromatin-associated features examined, histone H3 Lys27 acetylation levels at individual loci in mitosis best predict the mitosis–G1 transcriptional spike. Single-molecule RNA imaging supports that the mitosis–G1 transcriptional spike can constitute the maximum transcriptional activity per DNA copy throughout the cell division cycle. The transcriptional spike occurs heterogeneously and propagates to cell-to-cell differences in mature mRNA expression. Our results raise the possibility that passage through the mitosis–G1 transition might predispose cells to diverge in gene expression states.

Epo reprograms the epigenome of erythroid cells

The hormone erythropoietin (Epo) is required for erythropoiesis, yet its molecular mechanism of action remains poorly understood, particularly with respect to chromatin dynamics. To investigate how Epo modulates the erythroid epigenome, we performed epigenetic profiling using an ex vivo murine cell system that undergoes synchronous erythroid maturation in response to Epo stimulation. Our findings define the repertoire of Epo-modulated enhancers, illuminating a new facet of Epo signaling. First, a large number of enhancers rapidly responded to Epo stimulation, revealing a cis-regulatory network of Epo-responsive enhancers. In contrast, most of the other identified enhancers remained in an active acetylated state during Epo signaling, suggesting that most erythroid enhancers are established at an earlier precursor stage. Second, we identified several hundred super-enhancers that were linked to key erythroid genes, such as Tal1, Bcl11a, and Mir144/451. Third, experimental and computational validation revealed that many predicted enhancer regions were occupied by TAL1 and enriched with DNA-binding motifs for GATA1, KLF1, TAL1/E-box, and STAT5. Additionally, many of these cis-regulatory regions were conserved evolutionarily and displayed correlated enhancer:promoter acetylation. Together, these findings define a cis-regulatory enhancer network for Epo signaling during erythropoiesis, and provide the framework for future studies involving the interplay of epigenetics and Epo signaling.

The pioneer factor OCT4 requires the chromatin remodeller BRG1 to support gene regulatory element function in mouse embryonic stem cells

Pioneer transcription factors recognise and bind their target sequences in inaccessible chromatin to establish new transcriptional networks throughout development and cellular reprogramming. During this process, pioneer factors establish an accessible chromatin state to facilitate additional transcription factor binding, yet it remains unclear how different pioneer factors achieve this. Here, we discover that the pluripotency-associated pioneer factor OCT4 binds chromatin to shape accessibility, transcription factor co-binding, and regulatory element function in mouse embryonic stem cells. Chromatin accessibility at OCT4-bound sites requires the chromatin remodeller BRG1, which is recruited to these sites by OCT4 to support additional transcription factor binding and expression of the pluripotency-associated transcriptome. Furthermore, the requirement for BRG1 in shaping OCT4 binding reflects how these target sites are used during cellular reprogramming and early mouse development. Together this reveals a distinct requirement for a chromatin remodeller in promoting the activity of the pioneer factor OCT4 and regulating the pluripotency network.

DOI:http://dx.doi.org/10.7554/eLife.22631.001

All cells in your body contain the same genetic information in the form of genes encoded within DNA. Yet, cells use this information in different ways so that the activities of individual genes within that DNA can vary from cell to cell. This allows identical cells to become different to each other and to adapt to changing circumstances.

A group of proteins called transcription factors control the activity of certain genes by binding to specific sites on DNA. However, this isn’t a straightforward process because DNA in human and other animal cells is usually associated with structures called nucleosomes that can block access to the DNA. Pioneer transcription factors, such as OCT4, are a specific group of transcription factors that can attach to DNA in spite of the nucleosomes, but it’s not clear how this is possible. Once pioneer transcription factors attach to DNA they can help other transcription factors to bind alongside them.

King et al. studied OCT4 in stem cells from mouse embryos to investigate how it is able to act as a pioneer transcription factor and control gene activity. The experiments show that several other transcription factors lose the ability to bind to DNA when OCT4 is absent. This leads to widespread changes in gene activity in the cells, which seems to be due to other transcription factors being unable to get past the nucleosomes to attach to the DNA.

Further experiments showed that OCT4 needs a protein called BRG1 in order to act as a pioneer transcription factor. BRG1 is an enzyme that is able to move and remove (remodel) nucleosomes attached to DNA, suggesting that normal transcription factor binding requires this activity. The next challenge is to investigate whether BRG1, or similar enzymes, are also needed by other pioneer transcription factors that are required for normal gene activity and cell identity. This will be important because many enzymes that remodel nucleosomes are disrupted in human diseases like cancer where cells lose their normal identity.

DOI:http://dx.doi.org/10.7554/eLife.22631.002

Identification of recurrent combinatorial patterns of chromatin modifications at promoters across various tissue types

Background

Identification and analysis of recurrent combinatorial patterns of multiple chromatin modifications provide invaluable information for understanding epigenetic regulations. Furthermore, as more data becomes available, it is computationally expensive and unnecessary to study combinatorial patterns of all modifications.

Methods

A novel framework is proposed to investigate recurrent combinatorial patterns of a subset of quantitatively selected chromatin modifications. The framework is based on heirarchical clustering and selects subsets of chromatin modifications that form distinct recurrent patterns at regulatory regions. The identified recurrent combinatorial patterns can be further utilized to discover novel regulatory regions. Data is in the form of genome wide maps of histone acetylations, methylations, and histone variant of human skeletal muscular and B-lymphocyte cells both derived from the ENCODE project.

Results

A case study conducted at promoter regions is presented: four out of twelve chromatin modifications were selected, eight different promoter states were identified and the identified patterns of active promoters were further utilized to discover novel promoter regions. Several previously un-annotated promoters were discovered, further investigations confirm their promoter functions.

Conclusions

This framework is approproiately general and could lead to better understanding of epigenetic regulations by discovering previously unknown regulatory regions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1346-5) contains supplementary material, which is available to authorized users.

Widespread and dynamic translational control of red blood cell development

Cell development requires tight yet dynamic control of protein production. Here, we use parallel RNA and ribosome profiling to study translational regulatory dynamics during murine terminal erythropoiesis. Our results uncover pervasive translational control of protein synthesis, with widespread alternative translation initiation and termination, robust discrimination of long noncoding from micropeptide-encoding RNAs, and dynamic use of upstream open reading frames. Further, we identify hundreds of messenger RNAs (mRNAs) whose translation efficiency is dynamically controlled during erythropoiesis and that enrich for target sites of RNA-binding proteins that are specific to hematopoietic cells, thus unraveling potential regulators of erythroid translational programs. A major such program involves enhanced decoding of specific mRNAs that are depleted in terminally differentiating/enucleating cells with decreasing transcriptional capacity. We find that RBM38, an erythroid-specific RNA-binding protein previously implicated in splicing, interacts with the general translation initiation factor eIF4G and promotes translation of a subset of these irreplaceable mRNAs. Inhibition of RBM38 compromises translation in erythroblasts and impairs their maturation, highlighting a key function for this protein during erythropoiesis. These findings thus reveal critical roles for dynamic translational control in supporting specialized mammalian cell formation.

Distinct gene expression program dynamics during erythropoiesis from human induced pluripotent stem cells compared with adult and cord blood progenitors

BACKGROUND

Human-induced pluripotent stem cells (hiPSCs) are a potentially invaluable resource for regenerative medicine, including the in vitro manufacture of blood products. HiPSC-derived red blood cells are an attractive therapeutic option in hematology, yet exhibit unexplained proliferation and enucleation defects that presently preclude such applications. We hypothesised that substantial differential regulation of gene expression during erythroid development accounts for these important differences between hiPSC-derived cells and those from adult or cord-blood progenitors. We thus cultured erythroblasts from each source for transcriptomic analysis to investigate differential gene expression underlying these functional defects.

RESULTS

Our high resolution transcriptional view of definitive erythropoiesis captures the regulation of genes relevant to cell-cycle control and confers statistical power to deploy novel bioinformatics methods. Whilst the dynamics of erythroid program elaboration from adult and cord blood progenitors were very similar, the emerging erythroid transcriptome in hiPSCs revealed radically different program elaboration compared to adult and cord blood cells. We explored the function of differentially expressed genes in hiPSC-specific clusters defined by our novel tunable clustering algorithms (SMART and Bi-CoPaM). HiPSCs show reduced expression of c-KIT and key erythroid transcription factors SOX6, MYB and BCL11A, strong HBZ-induction, and aberrant expression of genes involved in protein degradation, lysosomal clearance and cell-cycle regulation.

CONCLUSIONS

Together, these data suggest that hiPSC-derived cells may be specified to a primitive erythroid fate, and implies that definitive specification may more accurately reflect adult development. We have therefore identified, for the first time, distinct gene expression dynamics during erythroblast differentiation from hiPSCs which may cause reduced proliferation and enucleation of hiPSC-derived erythroid cells. The data suggest several mechanistic defects which may partially explain the observed aberrant erythroid differentiation from hiPSCs.

Distinct transcription factor complexes act on a permissive chromatin landscape to establish regionalized gene expression in CNS stem cells

Spatially distinct gene expression profiles in neural stem cells (NSCs) are a prerequisite to the formation of neuronal diversity, but how these arise from the regulatory interactions between chromatin accessibility and transcription factor activity has remained unclear. Here, we demonstrate that, despite their distinct gene expression profiles, NSCs of the mouse cortex and spinal cord share the majority of their DNase I hypersensitive sites (DHSs). Regardless of this similarity, domain-specific gene expression is highly correlated with the relative accessibility of associated DHSs, as determined by sequence read density. Notably, the binding pattern of the general NSC transcription factor SOX2 is also largely cell type specific and coincides with an enrichment of LHX2 motifs in the cortex and HOXA9 motifs in the spinal cord. Interestingly, in a zebrafish reporter gene system, these motifs were critical determinants of patterned gene expression along the rostral-caudal axis. Our findings establish a predictive model for patterned NSC gene expression, whereby domain-specific expression of LHX2 and HOX proteins act on their target motifs within commonly accessible cis-regulatory regions to specify SOX2 binding. In turn, this binding correlates strongly with these DHSs relative accessibility—a robust predictor of neighboring gene expression.

Unlinking an lncRNA from Its Associated cis Element

Long non-coding (lnc) RNAs can regulate gene expression and protein functions. However, the proportion of lncRNAs with biological activities among the thousands expressed in mammalian cells is controversial. We studied Lockd (lncRNA downstream of Cdkn1b), a 434-nt polyadenylated lncRNA originating 4 kb 3' to the Cdkn1b gene. Deletion of the 25-kb Lockd locus reduced Cdkn1b transcription by approximately 70% in an erythroid cell line. In contrast, homozygous insertion of a polyadenylation cassette 80 bp downstream of the Lockd transcription start site reduced the entire lncRNA transcript level by >90% with no effect on Cdkn1b transcription. The Lockd promoter contains a DNase-hypersensitive site, binds numerous transcription factors, and physically associates with the Cdkn1b promoter in chromosomal conformation capture studies. Therefore, the Lockd gene positively regulates Cdkn1b transcription through an enhancer-like cis element, whereas the lncRNA itself is dispensable, which may be the case for other lncRNAs.

Insight into GATA1 transcriptional activity through interrogation of cis elements disrupted in human erythroid disorders

Whole-exome sequencing has been incredibly successful in identifying causal genetic variants and has revealed a number of novel genes associated with blood and other diseases. One limitation of this approach is that it overlooks mutations in noncoding regulatory elements. Furthermore, the mechanisms by which mutations in transcriptionalcis-regulatory elements result in disease remain poorly understood. Here we used CRISPR/Cas9 genome editing to interrogate three such elements harboring mutations in human erythroid disorders, which in all cases are predicted to disrupt a canonical binding motif for the hematopoietic transcription factor GATA1. Deletions of as few as two to four nucleotides resulted in a substantial decrease (>80%) in target gene expression. Isolated deletions of the canonical GATA1 binding motif completely abrogated binding of the cofactor TAL1, which binds to a separate motif. Having verified the functionality of these three GATA1 motifs, we demonstrate strong evolutionary conservation of GATA1 motifs in regulatory elements proximal to other genes implicated in erythroid disorders, and show that targeted disruption of such elements results in altered gene expression. By modeling transcription factor binding patterns, we show that multiple transcription factors are associated with erythroid gene expression, and have created predictive maps modeling putative disruptions of their binding sites at key regulatory elements. Our study provides insight into GATA1 transcriptional activity and may prove a useful resource for investigating the pathogenicity of noncoding variants in human erythroid disorders.

The DEK Oncoprotein Is a Critical Component of the EKLF/KLF1 Enhancer in Erythroid Cells

Understanding how transcriptional regulators are themselves controlled is important in attaining a complete picture of the intracellular effects that follow signaling cascades during early development and cell-restricted differentiation. We have addressed this issue by focusing on the regulation of EKLF/KLF1, a zinc finger transcription factor that plays a necessary role in the global regulation of erythroid gene expression. Using biochemical affinity purification, we have identified the DEK oncoprotein as a critical factor that interacts with an essential upstream enhancer element of the EKLF promoter and exerts a positive effect on EKLF levels. This element also binds a core set of erythroid transcription factors, suggesting that DEK is part of a tissue-restricted enhanceosome that contains BMP4-dependent and -independent components. Together with local enrichment of properly coded histones and an open chromatin domain, optimal transcriptional activation of the EKLF locus can be established.

Genome-Wide Organization of GATA1 and TAL1 Determined at High Resolution

Erythroid development and differentiation from multiprogenitor cells into red blood cells requires precise transcriptional regulation. Key erythroid transcription factors, GATA1 and TAL1, cooperate, along with other proteins, to regulate many aspects of this process. How GATA1 and TAL1 are juxtaposed along the DNA and their cognate DNA binding site across the mouse genome remains unclear. We applied high-resolution ChIP-exo (chromatin immunoprecipitation followed by 5′-to-3′ exonuclease treatment and then massively parallel DNA sequencing) to GATA1 and TAL1 to study their positional organization across the mouse genome during GATA1-dependent maturation. Two complementary methods, MultiGPS and peak pairing, were used to determine high-confidence binding locations by ChIP-exo. We identified ∼10,000 GATA1 and ∼15,000 TAL1 locations, which were essentially confirmed by ChIP-seq (chromatin immunoprecipitation followed by massively parallel DNA sequencing). Of these, ∼4,000 locations were bound by both GATA1 and TAL1. About three-quarters of them were tightly linked to a partial E-box located 7 or 8 bp upstream of a WGATAA motif. Both TAL1 and GATA1 generated distinct characteristic ChIP-exo peaks around WGATAA motifs that reflect their positional arrangement within a complex. We show that TAL1 and GATA1 form a precisely organized complex at a compound motif consisting of a TG 7 or 8 bp upstream of a WGATAA motif across thousands of genomic locations.

Heterogeneity of epigenetic changes at ischemia/reperfusion- and endotoxin-induced acute kidney injury genes

Aberrant gene expression is a molecular hallmark of acute kidney injury (AKI). As epigenetic processes control gene expression in a cell- and environment-defined manner, understanding the epigenetic pathways that regulate genes altered by AKI may open vital new insights into the complexities of disease pathogenesis and identify possible therapeutic targets. Here we used matrix chromatin immunoprecipitation and integrative analysis to study 20 key permissive and repressive epigenetic histone marks at transcriptionally induced Tnf, Ngal, Kim-1, and Icam-1 genes in mouse models of AKI; unilateral renal ischemia/reperfusion, lipopolysaccharide (LPS), and their synergistically injurious combination. Results revealed unexpected heterogeneity of transcriptional and epigenetic responses. Tnf and Ngal were transcriptionally upregulated in response to both treatments individually, and to combination treatment. Kim-1 was induced by ischemia/reperfusion and Icam-1 by LPS only. Epigenetic alterations at these genes exhibited distinct time-dependent changes that shared some similarities, such as reduction in repressive histone modifications, and also had major ischemia/reperfusion versus endotoxin differences. Thus, diversity of changes at AKI genes in response to different insults indicates involvement of several epigenetic pathways. This could be exploited pharmacologically through rational-drug design to alter the course and improve clinical outcomes of this syndrome.

Hematopoietic Signaling Mechanism Revealed from a Stem/Progenitor Cell Cistrome

Thousands of cis-elements in genomes are predicted to have vital functions. Although conservation, activity in surrogate assays, polymorphisms, and disease mutations provide functional clues, deletion from endogenous loci constitutes the gold-standard test. A GATA-2-binding, Gata2 intronic cis-element (+9.5) required for hematopoietic stem cell genesis in mice is mutated in a human immunodeficiency syndrome. Because +9.5 is the only cis-element known to mediate stem cell genesis, we devised a strategy to identify functionally comparable enhancers ("+9.5-like") genome-wide. Gene editing revealed +9.5-like activity to mediate GATA-2 occupancy, chromatin opening, and transcriptional activation. A +9.5-like element resided in Samd14, which encodes a protein of unknown function. Samd14 increased hematopoietic progenitor levels/activity and promoted signaling by a pathway vital for hematopoietic stem/progenitor cell regulation (stem cell factor/c-Kit), and c-Kit rescued Samd14 loss-of-function phenotypes. Thus, the hematopoietic stem/progenitor cell cistrome revealed a mediator of a signaling pathway that has broad importance for stem/progenitor cell biology.

Dynamics of GATA1 binding and expression response in a GATA1-induced erythroid differentiation system

During the maturation phase of mammalian erythroid differentiation, highly proliferative cells committed to the erythroid lineage undergo dramatic changes in morphology and function to produce circulating, enucleated erythrocytes. These changes are caused by equally dramatic alterations in gene expression, which in turn are driven by changes in the abundance and binding patterns of transcription factors such as GATA1. We have studied the dynamics of GATA1 binding by ChIP-seq and the global expression responses by RNA-seq in a GATA1-dependent mouse cell line model for erythroid maturation, in both cases examining seven progressive stages during differentiation. Analyses of these data should provide insights both into mechanisms of regulation (early versus late targets) and the consequences in cell physiology (e.g., distinctive categories of genes regulated at progressive stages of differentiation). The data are deposited in the Gene Expression Omnibus, series GSE36029, GSE40522, GSE49847, and GSE51338.

The role of chromatin dynamics in immune cell development

In immune cells, as in all mammalian cells, nuclear DNA is wrapped around histones to form nucleosomes. The positioning and modifications of nucleosomes throughout the genome defines the chromatin state of the cell and has a large impact on gene regulation. Chromatin state is dynamic throughout immune cell development and activation. High-throughput open chromatin assays, such as DNase-seq, can be used to find regulatory element across the genome and, when combined with histone modifications, can specify their function. During hematopoiesis, distal regulatory elements, known as enhancers, are established by pioneer factors that alter chromatin state. Some of these enhancers are lost, some are gained, and some are maintained as a memory of the cell's developmental origin. The enhancer landscape is unique to the cell lineage-with different enhancers regulating the same promoter-and determines the mechanism of cell type-specific activation after exposure to stimuli. Histone modification and promoter architecture govern the diverse responses to stimulation. Furthermore, chromatin dynamics may explain the high plasticity of certain tissue-resident immune cell types. Future epigenomic research will depend on the development of more efficient experiments and better methods to associate enhancers with genes. The ultimate goal of mapping genome-wide chromatin state throughout the hematopoietic tree will help illuminate the mechanisms behind immune cell development and function.