Integrated genome-scale analysis of the transcriptional regulatory landscape in a blood stem/progenitor cell model

Yin Yang 1 regulates cohesin complex protein SMC3 in mouse hematopoietic stem cells

ABSTRACT

Yin Yang 1 (YY1) and structural maintenance of chromosomes 3 (SMC3) are 2 critical chromatin structural factors that mediate long-distance enhancer-promoter interactions and promote developmentally regulated changes in chromatin architecture in hematopoietic stem/progenitor cells (HSPCs). Although YY1 has critical functions in promoting hematopoietic stem cell (HSC) self-renewal and maintaining HSC quiescence, SMC3 is required for proper myeloid lineage differentiation. However, many questions remain unanswered regarding how YY1 and SMC3 interact with each other and affect hematopoiesis. We found that YY1 physically interacts with SMC3 and cooccupies with SMC3 at a large cohort of promoters genome wide, and YY1 deficiency deregulates the genetic network governing cell metabolism. YY1 occupies the Smc3 promoter and represses SMC3 expression in HSPCs. Although deletion of 1 Smc3 allele partially restores HSC numbers and quiescence in YY1 knockout mice, Yy1-/-Smc3+/- HSCs fail to reconstitute blood after bone marrow transplant. YY1 regulates HSC metabolic pathways and maintains proper intracellular reactive oxygen species levels in HSCs, and this regulation is independent of the YY1-SMC3 axis. Our results establish a distinct YY1-SMC3 axis and its impact on HSC quiescence and metabolism.

Reprogramming of 3D genome structure underlying HSPC development in zebrafish

BACKGROUND

Development of hematopoietic stem and progenitor cells (HSPC) is a multi-staged complex process that conserved between zebrafish and mammals. Understanding the mechanism underlying HSPC development is a holy grail of hematopoietic biology, which is helpful for HSPC clinical application. Chromatin conformation plays important roles in transcriptional regulation and cell fate decision; however, its dynamic and role in HSPC development is poorly investigated.

METHODS

We performed chromatin structure and multi-omics dissection across different stages of HSPC developmental trajectory in zebrafish for the first time, including Hi-C, RNA-seq, ATAC-seq, H3K4me3 and H3K27ac ChIP-seq.

RESULTS

The chromatin organization of zebrafish HSPC resemble mammalian cells with similar hierarchical structure. We revealed the multi-scale reorganization of chromatin structure and its influence on transcriptional regulation and transition of cell fate during HSPC development. Nascent HSPC is featured by loose conformation with obscure structure at all layers. Notably, PU.1 was identified as a potential factor mediating formation of promoter-involved loops and regulating gene expression of HSPC.

CONCLUSIONS

Our results provided a global view of chromatin structure dynamics associated with development of zebrafish HSPC and discovered key transcription factors involved in HSPC chromatin interactions, which will provide new insights into the epigenetic regulatory mechanisms underlying vertebrate HSPC fate decision.

An NFIX-mediated regulatory network governs the balance of hematopoietic stem and progenitor cells during hematopoiesis

The transcription factor (TF) nuclear factor I-X (NFIX) is a positive regulator of hematopoietic stem and progenitor cell (HSPC) transplantation. Nfix-deficient HSPCs exhibit a severe loss of repopulating activity, increased apoptosis, and a loss of colony-forming potential. However, the underlying mechanism remains elusive. Here, we performed cellular indexing of transcriptomes and epitopes by high-throughput sequencing (CITE-seq) on Nfix-deficient HSPCs and observed a loss of long-term hematopoietic stem cells and an accumulation of megakaryocyte and myelo-erythroid progenitors. The genome-wide binding profile of NFIX in primitive murine hematopoietic cells revealed its colocalization with other hematopoietic TFs, such as PU.1. We confirmed the physical interaction between NFIX and PU.1 and demonstrated that the 2 TFs co-occupy super-enhancers and regulate genes implicated in cellular respiration and hematopoietic differentiation. In addition, we provide evidence suggesting that the absence of NFIX negatively affects PU.1 binding at some genomic loci. Our data support a model in which NFIX collaborates with PU.1 at super-enhancers to promote the differentiation and homeostatic balance of hematopoietic progenitors.

Unveiling the complexity of transcription factor networks in hematopoietic stem cells: implications for cell therapy and hematological malignancies

The functionality and longevity of hematopoietic tissue is ensured by a tightly controlled balance between self-renewal, quiescence, and differentiation of hematopoietic stem cells (HSCs) into the many different blood lineages. Cell fate determination in HSCs is influenced by signals from extrinsic factors (e.g., cytokines, irradiation, reactive oxygen species, O2 concentration) that are translated and integrated by intrinsic factors such as Transcription Factors (TFs) to establish specific gene regulatory programs. TFs also play a central role in the establishment and/or maintenance of hematological malignancies, highlighting the need to understand their functions in multiple contexts. TFs bind to specific DNA sequences and interact with each other to form transcriptional complexes that directly or indirectly control the expression of multiple genes. Over the past decades, significant research efforts have unraveled molecular programs that control HSC function. This, in turn, led to the identification of more than 50 TF proteins that influence HSC fate. However, much remains to be learned about how these proteins interact to form molecular networks in combination with cofactors (e.g. epigenetics factors) and how they control differentiation, expansion, and maintenance of cellular identity. Understanding these processes is critical for future applications particularly in the field of cell therapy, as this would allow for manipulation of cell fate and induction of expansion, differentiation, or reprogramming of HSCs using specific cocktails of TFs. Here, we review recent findings that have unraveled the complexity of molecular networks controlled by TFs in HSCs and point towards possible applications to obtain functional HSCs ex vivo for therapeutic purposes including hematological malignancies. Furthermore, we discuss the challenges and prospects for the derivation and expansion of functional adult HSCs in the near future.

Histone H2A deubiquitinases in the transcriptional programs of development and hematopoiesis: a consolidated analysis

Monoubiquitinated lysine 119 of histone H2A (H2AK119ub) is a highly abundant epigenetic mark, associated with gene repression and deposited on chromatin by the polycomb repressor complex 1 (PRC1), which is an essential regulator of diverse transcriptional programs in mammalian development and tissue homeostasis. While multiple deubiquitinases (DUBs) with catalytic activity for H2AK119ub (H2A-DUBs) have been identified, we lack systematic analyses of their roles and cross-talk in transcriptional regulation. Here, we address H2A-DUB functions in epigenetic regulation of mammalian development and tissue maintenance by conducting a meta-analysis of 248 genomics datasets from 32 independent studies, focusing on the mouse model and covering embryonic stem cells (ESCs), hematopoietic, and immune cell lineages. This covers all the publicly available datasets that map genomic H2A-DUB binding and H2AK119ub distributions (ChIP-Seq), and all datasets assessing dysregulation in gene expression in the relevant H2A-DUB knockout models (RNA-Seq). Many accessory datasets for PRC1-2 and DUB-interacting proteins are also analyzed and interpreted, as well as further data assessing chromatin accessibility (ATAC-Seq) and transcriptional activity (RNA-seq). We report co-localization in the binding of H2A-DUBs BAP1, USP16, and to a lesser extent others that is conserved across different cell-types, and also the enrichment of antagonistic PRC1-2 protein complexes at the same genomic locations. Such conserved sites enriched for the H2A-DUBs and PRC1-2 are proximal to transcriptionally active genes that engage in housekeeping cellular functions. Nevertheless, they exhibit H2AK119ub levels significantly above the genomic average that can undergo further increase with H2A-DUB knockout. This indicates a cooperation between H2A-DUBs and PRC1-2 in the modulation of housekeeping transcriptional programs, conserved across many cell types, likely operating through their antagonistic effects on H2AK119ub and the regulation of local H2AK119ub turnover. Our study further highlights existing knowledge gaps and discusses important directions for future work.

A genome-wide relay of signalling-responsive enhancers drives hematopoietic specification

Developmental control of gene expression critically depends on distal cis-regulatory elements including enhancers which interact with promoters to activate gene expression. To date no global experiments have been conducted that identify their cell type and cell stage-specific activity within one developmental pathway and in a chromatin context. Here, we describe a high-throughput method that identifies thousands of differentially active cis-elements able to stimulate a minimal promoter at five stages of hematopoietic progenitor development from embryonic stem (ES) cells, which can be adapted to any ES cell derived cell type. We show that blood cell-specific gene expression is controlled by the concerted action of thousands of differentiation stage-specific sets of cis-elements which respond to cytokine signals terminating at signalling responsive transcription factors. Our work provides an important resource for studies of hematopoietic specification and highlights the mechanisms of how and where extrinsic signals program a cell type-specific chromatin landscape driving hematopoietic differentiation.

Dynamic Runx1 chromatin boundaries affect gene expression in hematopoietic development

The transcription factor RUNX1 is a critical regulator of developmental hematopoiesis and is frequently disrupted in leukemia. Runx1 is a large, complex gene that is expressed from two alternative promoters under the spatiotemporal control of multiple hematopoietic enhancers. To dissect the dynamic regulation of Runx1 in hematopoietic development, we analyzed its three-dimensional chromatin conformation in mouse embryonic stem cell (ESC) differentiation cultures. Runx1 resides in a 1.1 Mb topologically associating domain (TAD) demarcated by convergent CTCF motifs. As ESCs differentiate to mesoderm, chromatin accessibility, Runx1 enhancer-promoter (E-P) interactions, and CTCF-CTCF interactions increase in the TAD, along with initiation of Runx1 expression from the P2 promoter. Differentiation to hematopoietic progenitor cells is associated with the formation of tissue-specific sub-TADs over Runx1, a shift in E-P interactions, P1 promoter demethylation, and robust expression from both Runx1 promoters. Deletion of promoter-proximal CTCF sites at the sub-TAD boundaries has no obvious effects on E-P interactions but leads to partial loss of domain structure, mildly affects gene expression, and delays hematopoietic development. Together, our analysis of gene regulation at a large multi-promoter developmental gene reveals that dynamic sub-TAD chromatin boundaries play a role in establishing TAD structure and coordinated gene expression.

CDK19 regulates the proliferation of hematopoietic stem cells and acute myeloid leukemia cells by suppressing p53-mediated transcription of p21

The cell cycle progression of hematopoietic stem cells (HSCs) and acute myeloid leukemia (AML) cells is precisely controlled by multiple regulatory factors. However, the underlying mechanisms are not fully understood. Here, we find that cyclin-dependent kinase 19 (CDK19), not its paralogue CDK8, is relatively enriched in mouse HSCs, and its expression is more significantly increased than CDK8 after proliferative stresses. Furthermore, SenexinB (a CDK8/19 inhibitor) treatment impairs the proliferation and self-renewal ability of HSCs. Moreover, overexpression of CDK19 promotes HSC function better than CDK8 overexpression. Using CDK19 knockout mice, we observe that CDK19^-/- HSCs exhibit similar phenotypes to those of cells treated with SenexinB. Interestingly, the p53 signaling pathway is significantly activated in HSCs lacking CDK19 expression. Further investigations show that CDK19 can interact with p53 to inhibit p53-mediated transcription of p21 in HSCs and treatment with a specific p53 inhibitor (PFTβ) partially rescues the defects of CDK19-null HSCs. Importantly, SenexinB treatment markedly inhibits the proliferation of AML cells. Collectively, our findings indicate that CDK19 is involved in regulating HSC and AML cell proliferation via the p53-p21 pathway, revealing a new mechanism underlying cell cycle regulation in normal and malignant hematopoietic cells.

3D genome organization links non-coding disease-associated variants to genes

Genome sequencing has revealed over 300 million genetic variations in human populations. Over 90% of variants are single nucleotide polymorphisms (SNPs), the remainder include short deletions or insertions, and small numbers of structural variants. Hundreds of thousands of these variants have been associated with specific phenotypic traits and diseases through genome wide association studies which link significant differences in variant frequencies with specific phenotypes among large groups of individuals. Only 5% of disease-associated SNPs are located in gene coding sequences, with the potential to disrupt gene expression or alter of the function of encoded proteins. The remaining 95% of disease-associated SNPs are located in non-coding DNA sequences which make up 98% of the genome. The role of non-coding, disease-associated SNPs, many of which are located at considerable distances from any gene, was at first a mystery until the discovery that gene promoters regularly interact with distal regulatory elements to control gene expression. Disease-associated SNPs are enriched at the millions of gene regulatory elements that are dispersed throughout the non-coding sequences of the genome, suggesting they function as gene regulation variants. Assigning specific regulatory elements to the genes they control is not straightforward since they can be millions of base pairs apart. In this review we describe how understanding 3D genome organization can identify specific interactions between gene promoters and distal regulatory elements and how 3D genomics can link disease-associated SNPs to their target genes. Understanding which gene or genes contribute to a specific disease is the first step in designing rational therapeutic interventions.

p53-dependent induction of P2X7 on hematopoietic stem and progenitor cells regulates hematopoietic response to genotoxic stress

Stem and progenitor cells are the main mediators of tissue renewal and repair, both under homeostatic conditions and in response to physiological stress and injury. Hematopoietic system is responsible for the regeneration of blood and immune cells and is maintained by bone marrow-resident hematopoietic stem and progenitor cells (HSPCs). Hematopoietic system is particularly susceptible to injury in response to genotoxic stress, resulting in the risk of bone marrow failure and secondary malignancies in cancer patients undergoing radiotherapy. Here we analyze the in vivo transcriptional response of HSPCs to genotoxic stress in a mouse whole-body irradiation model and, together with p53 ChIP-Seq and studies in p53-knockout (p53KO) mice, characterize the p53-dependent and p53-independent branches of this transcriptional response. Our work demonstrates the p53-independent induction of inflammatory transcriptional signatures in HSPCs in response to genotoxic stress and identifies multiple novel p53-target genes induced in HSPCs in response to whole-body irradiation. In particular, we establish the direct p53-mediated induction of P2X7 expression on HSCs and HSPCs in response to genotoxic stress. We further demonstrate the role of P2X7 in hematopoietic response to acute genotoxic stress, with P2X7 deficiency significantly extending mouse survival in irradiation-induced hematopoietic failure. We also demonstrate the role of P2X7 in the context of long-term HSC regenerative fitness following sublethal irradiation. Overall our studies provide important insights into the mechanisms of HSC response to genotoxic stress and further suggest P2X7 as a target for pharmacological modulation of HSC fitness and hematopoietic response to genotoxic injury.

Transcriptional Regulation of RUNX1: An Informatics Analysis

The RUNX1/AML1 gene encodes a developmental transcription factor that is an important regulator of haematopoiesis in vertebrates. Genetic disruptions to the RUNX1 gene are frequently associated with acute myeloid leukaemia. Gene regulatory elements (REs), such as enhancers located in non-coding DNA, are likely to be important for Runx1 transcription. Non-coding elements that modulate Runx1 expression have been investigated over several decades, but how and when these REs function remains poorly understood. Here we used bioinformatic methods and functional data to characterise the regulatory landscape of vertebrate Runx1. We identified REs that are conserved between human and mouse, many of which produce enhancer RNAs in diverse tissues. Genome-wide association studies detected single nucleotide polymorphisms in REs, some of which correlate with gene expression quantitative trait loci in tissues in which the RE is active. Our analyses also suggest that REs can be variant in haematological malignancies. In summary, our analysis identifies features of the RUNX1 regulatory landscape that are likely to be important for the regulation of this gene in normal and malignant haematopoiesis.

Loss of MYSM1 inhibits the oncogenic activity of cMYC in B cell lymphoma

MYSM1 is a chromatin‐binding protein, widely investigated for its functions in haematopoiesis in human and mouse; however, its role in haematologic malignancies remains unexplored. Here, we investigate the cross‐talk between MYSM1 and oncogenic cMYC in the transcriptional regulation of genes encoding ribosomal proteins, and the implications of these mechanisms for cMYC‐driven carcinogenesis. We demonstrate that in cMYC‐driven B cell lymphoma in mouse models, MYSM1‐loss represses ribosomal protein gene expression and protein synthesis. Importantly, the loss of MYSM1 also strongly inhibits cMYC oncogenic activity and protects against B cell lymphoma onset and progression in the mouse models. This advances the understanding of the molecular and transcriptional mechanisms of lymphomagenesis, and suggests MYSM1 as a possible drug target for cMYC‐driven malignancies.

Effects of sheared chromatin length on ChIP-seq quality and sensitivity

Chromatin immunoprecipitation followed by massively parallel, high throughput sequencing (ChIP-seq) is the method of choice for genome-wide identification of DNA segments bound by specific transcription factors or in chromatin with particular histone modifications. However, the quality of ChIP-seq datasets varies widely, with a substantial fraction being of intermediate to poor quality. Thus, it is important to discern and control the factors that contribute to variation in ChIP-seq. In this study, we focused on sonication, a user-controlled variable, to produce sheared chromatin. We systematically varied the amount of shearing of fixed chromatin from a mouse erythroid cell line, carefully measuring the distribution of resultant fragment lengths prior to ChIP-seq. This systematic study was complemented with a retrospective analysis of additional experiments. We found that the level of sonication had a pronounced impact on the quality of ChIP-seq signals. Over-sonication consistently reduced quality, while the impact of under-sonication differed among transcription factors, with no impact on sites bound by CTCF but frequently leading to the loss of sites occupied by TAL1 or bound by POL2. The bound sites not observed in low quality datasets were inferred to be a mix of both direct and indirect binding. We leveraged these findings to produce a set of CTCF ChIP-seq datasets in rare, primary hematopoietic progenitor cells. Our observation that the amount of chromatin sonication is a key variable in success of ChIP-seq experiments indicates that monitoring the level of sonication can improve ChIP-seq quality and reproducibility and facilitate ChIP-seq in rare cell types.

Kmt2c mutations enhance HSC self-renewal capacity and convey a selective advantage after chemotherapy

The myeloid tumor suppressor KMT2C is recurrently deleted in myelodysplastic syndrome (MDS) and acute myeloid leukemia (AML), particularly therapy-related MDS/AML (t-MDS/t-AML), as part of larger chromosome 7 deletions. Here, we show that KMT2C deletions convey a selective advantage to hematopoietic stem cells (HSCs) after chemotherapy treatment that may precipitate t-MDS/t-AML. Kmt2c deletions markedly enhance murine HSC self-renewal capacity without altering proliferation rates. Haploid Kmt2c deletions convey a selective advantage only when HSCs are driven into cycle by a strong proliferative stimulus, such as chemotherapy. Cycling Kmt2c-deficient HSCs fail to differentiate appropriately, particularly in response to interleukin-1. Kmt2c deletions mitigate histone methylation/acetylation changes that accrue as HSCs cycle after chemotherapy, and they impair enhancer recruitment during HSC differentiation. These findings help explain why Kmt2c deletions are more common in t-MDS/t-AML than in de novo AML or clonal hematopoiesis: they selectively protect cycling HSCs from differentiation without inducing HSC proliferation themselves.

The autism-related protein CHD8 contributes to the stemness and differentiation of mouse hematopoietic stem cells

Chromodomain helicase DNA-binding protein 8 (CHD8) is an ATP-dependent chromatin-remodeling factor that is encoded by the most frequently mutated gene in individuals with autism spectrum disorder. CHD8 is expressed not only in neural tissues but also in many other organs; however, its functions are largely unknown. Here, we show that CHD8 is highly expressed in and maintains the stemness of hematopoietic stem cells (HSCs). Conditional deletion of Chd8 specifically in mouse bone marrow induces cell cycle arrest, apoptosis, and a differentiation block in HSCs in association with upregulation of the expression of p53 target genes. A colony formation assay and bone marrow transplantation reveal that CHD8 deficiency also compromises the stemness of HSCs. Furthermore, additional ablation of p53 rescues the impaired stem cell function and differentiation block of CHD8-deficient HSCs. Our results thus suggest that the CHD8-p53 axis plays a key role in regulation of the stemness and differentiation of HSCs.

Interactions between lineage-associated transcription factors govern haematopoietic progenitor states

Recent advances in molecular profiling provide descriptive datasets of complex differentiation landscapes including the haematopoietic system, but the molecular mechanisms defining progenitor states and lineage choice remain ill-defined. Here, we employed a cellular model of murine multipotent haematopoietic progenitors (Hoxb8-FL) to knock out 39 transcription factors (TFs) followed by RNA-Seq analysis, to functionally define a regulatory network of 16,992 regulator/target gene links. Focussed analysis of the subnetworks regulated by the B-lymphoid TF Ebf1 and T-lymphoid TF Gata3 revealed a surprising role in common activation of an early myeloid programme. Moreover, Gata3-mediated repression of Pax5 emerges as a mechanism to prevent precocious B-lymphoid differentiation, while Hox-mediated activation of Meis1 suppresses myeloid differentiation. To aid interpretation of large transcriptomics datasets, we also report a new method that visualises likely transitions that a progenitor will undergo following regulatory network perturbations. Taken together, this study reveals how molecular network wiring helps to establish a multipotent progenitor state, with experimental approaches and analysis tools applicable to dissecting a broad range of both normal and perturbed cellular differentiation landscapes.

Single-Cell Analysis of Neonatal HSC Ontogeny Reveals Gradual and Uncoordinated Transcriptional Reprogramming that Begins before Birth

Fetal and adult hematopoietic stem cells (HSCs) have distinct proliferation rates, lineage biases, gene expression profiles, and gene dependencies. Although these differences are widely recognized, it is not clear how the transition from fetal to adult identity is coordinated. Here we show that murine HSCs and committed hematopoietic progenitor cells (HPCs) undergo a gradual, rather than precipitous, transition from fetal to adult transcriptional states. The transition begins prior to birth and is punctuated by a late prenatal spike in type I interferon signaling that promotes perinatal HPC expansion and sensitizes progenitors to the leukemogenic FLT3ITD mutation. Most other changes in gene expression and enhancer activation are imprecisely timed and poorly coordinated. Thus, heterochronic enhancer elements, and their associated transcripts, are activated independently of one another rather than as part of a robust network. This simplifies the regulatory programs that guide neonatal HSC/HPC ontogeny, but it creates heterogeneity within these populations.

Lysine acetyltransferase Tip60 is required for hematopoietic stem cell maintenance

Hematopoietic stem cells (HSCs) have the potential to replenish the blood system for the lifetime of the organism. Their 2 defining properties, self-renewal and differentiation, are tightly regulated by the epigenetic machineries. Using conditional gene-knockout models, we demonstrated a critical requirement of lysine acetyltransferase 5 (Kat5, also known as Tip60) for murine HSC maintenance in both the embryonic and adult stages, which depends on its acetyltransferase activity. Genome-wide chromatin and transcriptome profiling in murine hematopoietic stem and progenitor cells revealed that Tip60 colocalizes with c-Myc and that Tip60 deletion suppress the expression of Myc target genes, which are associated with critical biological processes for HSC maintenance, cell cycling, and DNA repair. Notably, acetylated H2A.Z (acH2A.Z) was enriched at the Tip60-bound active chromatin, and Tip60 deletion induced a robust reduction in the acH2A.Z/H2A.Z ratio. These results uncover a critical epigenetic regulatory layer for HSC maintenance, at least in part through Tip60-dependent H2A.Z acetylation to activate Myc target genes.

CHD7 and Runx1 interaction provides a braking mechanism for hematopoietic differentiation

Hematopoietic stem and progenitor cell (HSPC) formation and lineage differentiation involve gene expression programs orchestrated by transcription factors and epigenetic regulators. Genetic disruption of the chromatin remodeler chromodomain-helicase-DNA-binding protein 7 (CHD7) expanded phenotypic HSPCs, erythroid, and myeloid lineages in zebrafish and mouse embryos. CHD7 acts to suppress hematopoietic differentiation. Binding motifs for RUNX and other hematopoietic transcription factors are enriched at sites occupied by CHD7, and decreased RUNX1 occupancy correlated with loss of CHD7 localization. CHD7 physically interacts with RUNX1 and suppresses RUNX1-induced expansion of HSPCs during development through modulation of RUNX1 activity. Consequently, the RUNX1:CHD7 axis provides proper timing and function of HSPCs as they emerge during hematopoietic development or mature in adults, representing a distinct and evolutionarily conserved control mechanism to ensure accurate hematopoietic lineage differentiation.

MYSM1 maintains ribosomal protein gene expression in hematopoietic stem cells to prevent hematopoietic dysfunction

Ribosomopathies are congenital disorders caused by mutations in the genes encoding ribosomal and other functionally related proteins. They are characterized by anemia, other hematopoietic and developmental abnormalities, and p53 activation. Ribosome assembly requires coordinated expression of many ribosomal protein (RP) genes; however, the regulation of RP gene expression, especially in hematopoietic stem cells (HSCs), remains poorly understood. MYSM1 is a transcriptional regulator essential for HSC function and hematopoiesis. We established that HSC dysfunction in Mysm1 deficiency is driven by p53; however, the mechanisms of p53 activation remained unclear. Here, we describe the transcriptome of Mysm1-deficient mouse HSCs and identify MYSM1 genome-wide DNA binding sites. We establish a direct role for MYSM1 in RP gene expression and show a reduction in protein synthesis in Mysm1-/- HSCs. Loss of p53 in mice fully rescues Mysm1-/- anemia phenotype but not RP gene expression, indicating that RP gene dysregulation is a direct outcome of Mysm1 deficiency and an upstream mediator of Mysm1-/- phenotypes through p53 activation. We characterize a patient with a homozygous nonsense MYSM1 gene variant, and we demonstrate reduced protein synthesis and increased p53 levels in patient hematopoietic cells. Our work provides insights into the specialized mechanisms regulating RP gene expression in HSCs and establishes a common etiology of MYSM1 deficiency and ribosomopathy syndromes.

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.

Induction of developmental hematopoiesis mediated by transcription factors and the hematopoietic microenvironment

The induction of hematopoiesis in various cell types via transcription factor (TF) reprogramming has been demonstrated by several strategies. The eventual goal of these approaches is to generate a product for unmet needs in hematopoietic cell transplantation therapies. The most successful strategies hew closely to clues provided from developmental hematopoiesis in terms of factor expression and environmental cues. In this review, we aim to summarize the TFs that play important roles in developmental hematopoiesis primarily and to also touch on adult hematopoiesis. Several aspects of cellular and molecular biology coalesce in this process, with TFs and surrounding cellular signals playing a major role in the overall development of the hematopoietic lineage. We attempt to put these elements into the context of reprogramming and highlight their roles.

Joint inference and alignment of genome structures enables characterization of compartment-independent reorganization across cell types

BACKGROUND

Comparisons of Hi-C data sets between cell types and conditions have revealed differences in topologically associated domains (TADs) and A/B compartmentalization, which are correlated with differences in gene regulation. However, previous comparisons have focused on known forms of 3D organization while potentially neglecting other functionally relevant differences. We aimed to create a method to quantify all locus-specific differences between two Hi-C data sets.

RESULTS

We developed MultiMDS to jointly infer and align 3D chromosomal structures from two Hi-C data sets, thereby enabling a new way to comprehensively quantify relocalization of genomic loci between cell types. We demonstrate this approach by comparing Hi-C data across a variety of cell types. We consistently find relocalization of loci with minimal difference in A/B compartment score. For example, we identify compartment-independent relocalizations between GM12878 and K562 cells that involve loci displaying enhancer-associated histone marks in one cell type and polycomb-associated histone marks in the other.

CONCLUSIONS

MultiMDS is the first tool to identify all loci that relocalize between two Hi-C data sets. Our method can identify 3D localization differences that are correlated with cell-type-specific regulatory activities and which cannot be identified using other methods.

PRDM16s transforms megakaryocyte-erythroid progenitors into myeloid leukemia-initiating cells

Oncogenic mutations confer on cells the ability to propagate indefinitely, but whether oncogenes alter the cell fate of these cells is unknown. Here, we show that the transcriptional regulator PRDM16s causes oncogenic fate conversion by transforming cells fated to form platelets and erythrocytes into myeloid leukemia stem cells (LSCs). Prdm16s expression in megakaryocyte-erythroid progenitors (MEPs), which normally lack the potential to generate granulomonocytic cells, caused AML by converting MEPs into LSCs. Prdm16s blocked megakaryocytic/erythroid potential by interacting with super enhancers and activating myeloid master regulators, including PU.1. A CRISPR dropout screen confirmed that PU.1 is required for Prdm16s-induced leukemia. Ablating PU.1 attenuated leukemogenesis and reinstated the megakaryocytic/erythroid potential of leukemic MEPs in mouse models and human AML with PRDM16 rearrangement. Thus, oncogenic PRDM16 s expression gives MEPs an LSC fate by activating myeloid gene regulatory networks.

Promoter Capture Hi-C: High-resolution, Genome-wide Profiling of Promoter Interactions

The three-dimensional organization of the genome is linked to its function. For example, regulatory elements such as transcriptional enhancers control the spatio-temporal expression of their target genes through physical contact, often bridging considerable (in some cases hundreds of kilobases) genomic distances and bypassing nearby genes. The human genome harbors an estimated one million enhancers, the vast majority of which have unknown gene targets. Assigning distal regulatory regions to their target genes is thus crucial to understand gene expression control. We developed Promoter Capture Hi-C (PCHi-C) to enable the genome-wide detection of distal promoter-interacting regions (PIRs), for all promoters in a single experiment. In PCHi-C, highly complex Hi-C libraries are specifically enriched for promoter sequences through in-solution hybrid selection with thousands of biotinylated RNA baits complementary to the ends of all promoter-containing restriction fragments. The aim is to then pull-down promoter sequences and their frequent interaction partners such as enhancers and other potential regulatory elements. After high-throughput paired-end sequencing, a statistical test is applied to each promoter-ligated restriction fragment to identify significant PIRs at the restriction fragment level. We have used PCHi-C to generate an atlas of long-range promoter interactions in dozens of human and mouse cell types. These promoter interactome maps have contributed to a greater understanding of mammalian gene expression control by assigning putative regulatory regions to their target genes and revealing preferential spatial promoter-promoter interaction networks. This information also has high relevance to understanding human genetic disease and the identification of potential disease genes, by linking non-coding disease-associated sequence variants in or near control sequences to their target genes.

UTX-mediated enhancer and chromatin remodeling suppresses myeloid leukemogenesis through noncatalytic inverse regulation of ETS and GATA programs

The histone H3 Lys27-specific demethylase UTX (or KDM6A) is targeted by loss-of-function mutations in multiple cancers. Here, we demonstrate that UTX suppresses myeloid leukemogenesis through noncatalytic functions, a property shared with its catalytically inactive Y-chromosome paralog, UTY (or KDM6C). In keeping with this, we demonstrate concomitant loss/mutation of KDM6A (UTX) and UTY in multiple human cancers. Mechanistically, global genomic profiling showed only minor changes in H3K27me3 but significant and bidirectional alterations in H3K27ac and chromatin accessibility; a predominant loss of H3K4me1 modifications; alterations in ETS and GATA-factor binding; and altered gene expression after Utx loss. By integrating proteomic and genomic analyses, we link these changes to UTX regulation of ATP-dependent chromatin remodeling, coordination of the COMPASS complex and enhanced pioneering activity of ETS factors during evolution to AML. Collectively, our findings identify a dual role for UTX in suppressing acute myeloid leukemia via repression of oncogenic ETS and upregulation of tumor-suppressive GATA programs.

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.

Thrombopoietin signaling to chromatin elicits rapid and pervasive epigenome remodeling within poised chromatin architectures

Thrombopoietin (TPO) is a critical cytokine regulating hematopoietic stem cell maintenance and differentiation into the megakaryocytic lineage. However, the transcriptional and chromatin dynamics elicited by TPO signaling are poorly understood. Here, we study the immediate early transcriptional and cis-regulatory responses to TPO in hematopoietic stem/progenitor cells (HSPCs) and use this paradigm of cytokine signaling to chromatin to dissect the relation between cis- regulatory activity and chromatin architecture. We show that TPO profoundly alters the transcriptome of HSPCs, with key hematopoietic regulators being transcriptionally repressed within 30 minutes of TPO. By examining cis-regulatory dynamics and chromatin architectures, we demonstrate that these changes are accompanied by rapid and extensive epigenome remodeling of cis-regulatory landscapes that is spatially coordinated within topologically associating domains (TADs). Moreover, TPO-responsive enhancers are spatially clustered and engage in preferential homotypic intra- and inter-TAD interactions that are largely refractory to TPO signaling. By further examining the link between cis-regulatory dynamics and chromatin looping, we show that rapid modulation of cis-regulatory activity is largely independent of chromatin looping dynamics. Finally, we show that, although activated and repressed cis-regulatory elements share remarkably similar DNA sequence compositions, transcription factor binding patterns accurately predict rapid cis-regulatory responses to TPO.

A CRISPR Dropout Screen Identifies Genetic Vulnerabilities and Therapeutic Targets in Acute Myeloid Leukemia

Acute myeloid leukemia (AML) is an aggressive cancer with a poor prognosis, for which mainstream treatments have not changed for decades. To identify additional therapeutic targets in AML, we optimize a genome-wide clustered regularly interspaced short palindromic repeats (CRISPR) screening platform and use it to identify genetic vulnerabilities in AML cells. We identify 492 AML-specific cell-essential genes, including several established therapeutic targets such as DOT1L, BCL2, and MEN1, and many other genes including clinically actionable candidates. We validate selected genes using genetic and pharmacological inhibition, and chose KAT2A as a candidate for downstream study. KAT2A inhibition demonstrated anti-AML activity by inducing myeloid differentiation and apoptosis, and suppressed the growth of primary human AMLs of diverse genotypes while sparing normal hemopoietic stem-progenitor cells. Our results propose that KAT2A inhibition should be investigated as a therapeutic strategy in AML and provide a large number of genetic vulnerabilities of this leukemia that can be pursued in downstream studies.

•

Optimized CRISPR platform for identification of genome-wide genetic vulnerabilities

•

Catalog of genetic vulnerabilities in acute myeloid leukemia cell lines

•

KAT2A inhibition induces myeloid differentiation and apoptosis

•

KAT2A inhibition arrests the growth of primary AML cells, but not of normal progenitors

Mammalian Transcription Factor Networks: Recent Advances in Interrogating Biological Complexity

Transcription factor (TF) networks are a key determinant of cell fate decisions in mammalian development and adult tissue homeostasis and are frequently corrupted in disease. However, our inability to experimentally resolve and interrogate the complexity of mammalian TF networks has hampered the progress in this field. Recent technological advances, in particular large-scale genome-wide approaches, single-cell methodologies, live-cell imaging, and genome editing, are emerging as important technologies in TF network biology. Several recent studies even suggest a need to re-evaluate established models of mammalian TF networks. Here, we provide a brief overview of current and emerging methods to define mammalian TF networks. We also discuss how these emerging technologies facilitate new ways to interrogate complex TF networks, consider the current open questions in the field, and comment on potential future directions and biomedical applications.

A DNA Contact Map for the Mouse Runx1 Gene Identifies Novel Haematopoietic Enhancers

The transcription factor Runx1 is essential for definitive haematopoiesis, and the RUNX1 gene is frequently translocated or mutated in leukaemia. Runx1 is transcribed from two promoters, P1 and P2, to give rise to different protein isoforms. Although the expression of Runx1 must be tightly regulated for normal blood development, the mechanisms that regulate Runx1 isoform expression during haematopoiesis remain poorly understood. Gene regulatory elements located in non-coding DNA are likely to be important for Runx1 transcription. Here we use circular chromosome conformation capture sequencing to identify DNA interactions with the P1 and P2 promoters of Runx1, and the previously identified +24 enhancer, in the mouse multipotent haematopoietic progenitor cell line HPC-7. The active promoter, P1, interacts with nine non-coding regions that are occupied by transcription factors within a 1 Mb topologically associated domain. Eight of nine regions function as blood-specific enhancers in zebrafish, of which two were previously shown to harbour blood-specific enhancer activity in mice. Interestingly, the +24 enhancer interacted with multiple distant regions on chromosome 16, suggesting it may regulate the expression of additional genes. The Runx1 DNA contact map identifies connections with multiple novel and known haematopoietic enhancers that are likely to be involved in regulating Runx1 expression in haematopoietic progenitor cells.

Inference of cell type specific regulatory networks on mammalian lineages

Transcriptional regulatory networks are at the core of establishing cell type specific gene expression programs. In mammalian systems, such regulatory networks are determined by multiple levels of regulation, including by transcription factors, chromatin environment, and three-dimensional organization of the genome. Recent efforts to measure diverse regulatory genomic datasets across multiple cell types and tissues offer unprecedented opportunities to examine the context-specificity and dynamics of regulatory networks at a greater resolution and scale than before. In parallel, numerous computational approaches to analyze these data have emerged that serve as important tools for understanding mammalian cell type specific regulation. In this article, we review recent computational approaches to predict the expression and sequence-based regulators of a gene's expression level and examine long-range gene regulation. We highlight promising approaches, insights gained, and open challenges that need to be overcome to build a comprehensive picture of cell type specific transcriptional regulatory networks.

SCL/TAL1: a multifaceted regulator from blood development to disease

SCL/TAL1 (stem cell leukemia/T-cell acute lymphoblastic leukemia [T-ALL] 1) is an essential transcription factor in normal and malignant hematopoiesis. It is required for specification of the blood program during development, adult hematopoietic stem cell survival and quiescence, and terminal maturation of select blood lineages. Following ectopic expression, SCL contributes to oncogenesis in T-ALL. Remarkably, SCL’s activities are all mediated through nucleation of a core quaternary protein complex (SCL:E-protein:LMO1/2 [LIM domain only 1 or 2]:LDB1 [LIM domain-binding protein 1]) and dynamic recruitment of conserved combinatorial associations of additional regulators in a lineage- and stage-specific context. The finely tuned control of SCL’s regulatory functions (lineage priming, activation, and repression of gene expression programs) provides insight into fundamental developmental and transcriptional mechanisms, and highlights mechanistic parallels between normal and oncogenic processes. Importantly, recent discoveries are paving the way to the development of innovative therapeutic opportunities in SCL+ T-ALL.

Single-Cell Analysis Identifies Distinct Stages of Human Endothelial-to-Hematopoietic Transition

During development, hematopoietic cells originate from endothelium in a process known as endothelial-to-hematopoietic transition (EHT). To study human EHT, we coupled flow cytometry and single-cell transcriptional analyses of human pluripotent stem cell-derived CD34+ cells. The resulting transcriptional hierarchy showed a continuum of endothelial and hematopoietic signatures. At the interface of these two signatures, a unique group of cells displayed both an endothelial signature and high levels of key hematopoietic stem cell-associated genes. This interphase group was validated via sort and subculture as an immediate precursor to hematopoietic cells. Differential expression analyses further divided this population into subgroups, which, upon subculture, showed distinct hematopoietic lineage differentiation potentials. We therefore propose that immediate precursors to hematopoietic cells already have their hematopoietic lineage restrictions defined prior to complete downregulation of the endothelial signature. These findings increase our understanding of the processes of de novo hematopoietic cell generation in the human developmental context.

Expression and regulation of C/EBPα in normal myelopoiesis and in malignant transformation

One of the most studied transcription factors in hematopoiesis is the leucine zipper CCAAT-enhancer binding protein α (C/EBPα), which is mainly involved in cell fate decisions for myeloid differentiation. Its involvement in acute myeloid leukemia (AML) is diverse, with patients frequently exhibiting mutations, deregulation of gene expression, or alterations in the function of C/EBPα. In this review, we emphasize the importance of C/EBPα for neutrophil maturation, its role in myeloid priming of hematopoietic stem and progenitor cells, and its indispensable requirement for AML development. We discuss that mutations in the open reading frame of CEBPA lead to an altered C/EBPα function, affecting the expression of downstream genes and consequently deregulating myelopoiesis. The emerging transcriptional mechanisms of CEBPA are discussed based on recent studies. Novel insights on how these mechanisms may be deregulated by oncoproteins or mutations/variants in CEBPA enhancers are suggested in principal to reveal novel mechanisms of how CEBPA is deregulated at the transcriptional level.

Capturing genomic relationships that matter

There is a strong interrelationship within the cell nucleus between form and function of the genome. This connection is exhibited across multiple hierarchies, ranging from grand-scale positioning of chromosomes and their intersection with specific nuclear functional activities, the segregation of chromosome structure into distinct domains and long-range regulatory contacts that drive spatial and temporal expression patterns of genes. Fifteen years ago, the development of the chromosome conformation capture method placed the nature of specific, long-range regulatory interactions under scrutiny. However, its development and integration with next-generation sequencing technologies has greatly expanded the breadth and scope of what is detected. The sheer scale of data offered by these important advances has come with new and challenging bottlenecks that are both experimental and bioinformatical. Here, we discuss the recent and prospective development and implementation of new methodologies and analytical tools that are allowing an in-depth, yet focussed characterisation of genomic contacts that are associated with functional activities in the nucleus.

Comprehensive population-based genome sequencing provides insight into hematopoietic regulatory mechanisms

Genetic variants affecting hematopoiesis can influence commonly measured blood cell traits. To identify factors that affect hematopoiesis, we performed association studies for blood cell traits in the population-based Estonian Biobank using high-coverage whole-genome sequencing (WGS) in 2,284 samples and SNP genotyping in an additional 14,904 samples. Using up to 7,134 samples with available phenotype data, our analyses identified 17 associations across 14 blood cell traits. Integration of WGS-based fine-mapping and complementary epigenomic datasets provided evidence for causal mechanisms at several loci, including at a previously undiscovered basophil count-associated locus near the master hematopoietic transcription factor CEBPA The fine-mapped variant at this basophil count association near CEBPA overlapped an enhancer active in common myeloid progenitors and influenced its activity. In situ perturbation of this enhancer by CRISPR/Cas9 mutagenesis in hematopoietic stem and progenitor cells demonstrated that it is necessary for and specifically regulates CEBPA expression during basophil differentiation. We additionally identified basophil count-associated variation at another more pleiotropic myeloid enhancer near GATA2, highlighting regulatory mechanisms for ordered expression of master hematopoietic regulators during lineage specification. Our study illustrates how population-based genetic studies can provide key insights into poorly understood cell differentiation processes of considerable physiologic relevance.

A CRISPR Dropout Screen Identifies Genetic Vulnerabilities and Therapeutic Targets in Acute Myeloid Leukemia

Acute myeloid leukemia (AML) is an aggressive cancer with a poor prognosis, for which mainstream treatments have not changed for decades. To identify additional therapeutic targets in AML, we optimize a genome-wide clustered regularly interspaced short palindromic repeats (CRISPR) screening platform and use it to identify genetic vulnerabilities in AML cells. We identify 492 AML-specific cell-essential genes, including several established therapeutic targets such as DOT1L, BCL2, and MEN1, and many other genes including clinically actionable candidates. We validate selected genes using genetic and pharmacological inhibition, and chose KAT2A as a candidate for downstream study. KAT2A inhibition demonstrated anti-AML activity by inducing myeloid differentiation and apoptosis, and suppressed the growth of primary human AMLs of diverse genotypes while sparing normal hemopoietic stem-progenitor cells. Our results propose that KAT2A inhibition should be investigated as a therapeutic strategy in AML and provide a large number of genetic vulnerabilities of this leukemia that can be pursued in downstream studies.

Finding partners to play the music of regulation

In this issue of Blood, Wilson et al generate and analyze a treasure trove of epigenetic data, such as transcription factor occupancy, histone modifications, and chromatin interaction frequencies, genome-wide (ie, epigenomic data), in a cell line model of hematopoietic stem/progenitor cells (HSPCs). To appreciate the importance of these data, consider an analogy of gene expression being a song or symphony (transcripts) played by musicians (transcription factors and transcriptional machinery) reading the score encoded in the genome sequence. Previous studies revealed the positions of a few transcription factors across the genome, so we only knew about, for example, the violinists and oboists. No wonder we did not understand how the music was being generated (how expression was regulated). By mapping the sites of occupancy of many more transcription factors (now a total of 29), as well as positions of 4 histone modifications and DNase hypersensitive sites, Wilson et al reveal many more of the players and their partners. Furthermore, their data on 3-dimensional interaction frequencies of chromatin show how groups of musicians (protein complexes) come together in an orchestra to read the score and perform a symphony.