Genome-wide analysis of simultaneous GATA1/2, RUNX1, FLI1, and SCL binding in megakaryocytes identifies hematopoietic regulators

BRD9 regulates normal human hematopoietic stem cell function and lineage differentiation

Bromodomain containing protein 9 (BRD9), a member of the non-canonical BRG1/BRM-associated factor (ncBAF) chromatin remodeling complex, has been implicated as a synthetic lethal target in AML but its function in normal human hematopoiesis is unknown. In hematopoietic stem and progenitor cells (HSPC) genomic or chemical inhibition of BRD9 led to a proliferative disadvantage and loss of stem cells in vitro. Human HSPCs with reduced BRD9 protein levels produced lower numbers of immature mixed multipotent GEMM colonies in semi-solid media. In lineage-promoting culture conditions, cells with reduced BRD9 levels failed to differentiate into the megakaryocytic lineage and showed delayed differentiation into erythroid cells but enhanced terminal myeloid differentiation. HSPCs with BRD9 knock down (KD) had reduced long-term multilineage engraftment in a xenotransplantation assay. An increased number of downregulated genes in RNAseq analysis after BRD9 KD coupled with a gain in chromatin accessibility at the promoters of several repressive transcription factors (TF) suggest that BRD9 functions in the maintenance of active transcription during HSC differentiation. In particular, the hematopoietic master regulator GATA1 was identified as one of the core TFs regulating the gene networks modulated by BRD9 loss in HSPCs. BRD9 inhibition reduced a GATA1-luciferase reporter signal, further suggesting a role for BRD9 in regulating GATA1 activity. BRD9 is therefore an additional example of epigenetic regulation of human hematopoiesis.

RUNX1 Isoforms Regulate RUNX1 and Target-Genes Differentially in Platelets-Megakaryocytes: Association with Clinical Cardiovascular Events

Background

Hematopoietic transcription factor RUNX1 is expressed from proximal P2 and distal P1 promoter to yield isoforms RUNX1 B and C, respectively. The roles of these isoforms in RUNX1 autoregulation and downstream-gene regulation in megakaryocytes and platelets are unknown.

Objectives

To understand the regulation of RUNX1 and its target genes by RUNX1 isoforms.

Methods

We performed studies on RUNX1 isoforms in megakaryocytic HEL cells and HeLa cells (lack endogenous RUNX1), in platelets from 85 healthy volunteers administered aspirin or ticagrelor, and on the association of RUNX1 target genes with acute events in 587 patients with cardiovascular disease (CVD).

Results

In chromatin immunoprecipitation and luciferase promoter assays, RUNX1 isoforms B and C bound and regulated P1 and P2 promoters. In HeLa cells RUNX1B decreased and RUNX1C increased P1 and P2 activities, respectively. In HEL cells, RUNX1B overexpression decreased RUNX1C and RUNX1A expression; RUNX1C increased RUNX1B and RUNX1A. RUNX1B and RUNX1C regulated target genes (MYL9, F13A1, PCTP, PDE5A and others) differentially in HEL cells. In platelets RUNX1B transcripts (by RNAseq) correlated negatively with RUNX1C and RUNX1A; RUNX1C correlated positively with RUNX1A. RUNX1B correlated positively with F13A1, PCTP, PDE5A, RAB1B, and others, and negatively with MYL9. In our previous studies, RUNX1C transcripts in whole blood were protective against acute events in CVD patients. We found that higher expression of RUNX1 targets F13A1 and RAB31 associated with acute events.

Conclusions

RUNX1 isoforms B and C autoregulate RUNX1 and regulate downstream genes in a differential manner and this associates with acute events in CVD.

Hematopoietic stem cell division is governed by distinct RUNX1 binding partners

A defined number of hematopoietic stem cell (HSC) clones are born during development and expand to form the pool of adult stem cells. An intricate balance between self-renewal and differentiation of these HSCs supports hematopoiesis for life. HSC fate is determined by complex transcription factor networks that drive cell-type specific gene programs. The transcription factor RUNX1 is required for definitive hematopoiesis, and mutations in Runx1 have been shown to reduce clonal diversity. The RUNX1 cofactor, CBFý, stabilizes RUNX1 binding to DNA, and disruption of their interaction alters downstream gene expression. Chemical screening for modulators of Runx1 and HSC expansion in zebrafish led us to identify a new mechanism for the RUNX1 inhibitor, Ro5-3335. We found that Ro5-3335 increased HSC divisions in zebrafish, and animals transplanted with Ro5-3335 treated cells had enhanced chimerism compared to untreated cells. Using human CD34+ cells, we show that Ro5-3335 remodels the RUNX1 transcription complex by binding to ELF1, independent of CBFý. This allows specific expression of cell cycle and hematopoietic genes that enhance HSC self-renewal and prevent differentiation. Furthermore, we provide the first evidence to show that it is possible to pharmacologically increase the number of stem cell clones in vivo , revealing a previously unknown mechanism for enhancing clonal diversity. Our studies have revealed a mechanism by which binding partners of RUNX1 determine cell fate, with ELF transcription factors guiding cell division. This information could lead to treatments that enhance clonal diversity for blood diseases.

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

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

Circadian regulation of stereotypic chromatin conformations at enhancers

Cooperation between the circadian transcription factor (TF) CLOCK:BMAL1 and other TFs at cis-regulatory elements (CREs) is critical to daily rhythms of transcription. Yet, the modalities of this cooperation are unclear. Here, we analyzed the co-binding of multiple TFs on single DNA molecules in mouse liver using single molecule footprinting (SMF). We found that SMF reads clustered in stereotypic chromatin states that reflect distinguishable organization of TFs and nucleosomes, and that were remarkably conserved between all samples. DNA protection at CLOCK:BMAL1 binding motif (E-box) varied between CREs, from E-boxes being solely bound by CLOCK:BMAL1 to situations where other TFs competed with CLOCK:BMAL1 for E-box binding. SMF also uncovered CLOCK:BMAL1 cooperative binding at E-boxes separated by 250 bp, which structurally altered the CLOCK:BMAL1-DNA interface. Importantly, we discovered multiple nucleosomes with E-boxes at entry/exit sites that were removed upon CLOCK:BMAL1 DNA binding, thereby promoting the formation of open chromatin states that facilitate DNA binding of other TFs and that were associated with rhythmic transcription. These results demonstrate the utility of SMF for studying how CLOCK:BMAL1 and other TFs regulate stereotypical chromatin states at CREs to promote transcription.

An update on recent studies of extracellular vesicles and their role in hypercoagulability in thalassemia (Review)

Thromboembolic events are a significant clinical concern in thalassemia and hemoglobinopathies, highlighting the need for new strategies to treat and detect these specific hematologic complications. In recent years, extracellular vesicles (EVs) have garnered interest due to their role in cell-to-cell communication, including angiogenesis, immune responses and coagulation activation. Their multifaceted role depends on the cellular origin and cargo, making them potential diagnostic biomarkers and therapeutic agents. The present review highlights recent advances in understanding the involvement of EVs in hypercoagulability in thalassemia, the characterization of circulating EVs and the potential for using EVs as predictive biomarkers. β-Thalassemia intermedia exhibits a high incidence of thromboembolic events, contributing to significant morbidity and mortality. Advanced technologies have enabled the profiling and characterization of circulating EVs in patients with β-thalassemia through various techniques, including flow cytometry, proteomic studies, reverse transcription-quantitative PCR, transmission electron microscopy, nanoparticle tracking analysis and western blot analysis. Microparticles from splenectomized β-thalassemia/hemoglobin E patients induce platelet activation and aggregation, potentially contributing to thrombus formation. The abundance of these microparticles, primarily released from platelets and damaged red cells, may have a role in thromboembolic events and other clinical complications in thalassemia. This suggests a promising future for EVs as diagnostic and predictive biomarkers in thalassemia management.

Genomic landscape of patients with germline RUNX1 variants and familial platelet disorder with myeloid malignancy

ABSTRACT

Familial platelet disorder with associated myeloid malignancies (FPDMM) is caused by germline RUNX1 mutations and characterized by thrombocytopenia and increased risk of hematologic malignancies. We recently launched a longitudinal natural history study for patients with FPDMM. Among 27 families with research genomic data by the end of 2021, 26 different germline RUNX1 variants were detected. Besides missense mutations enriched in Runt homology domain and loss-of-function mutations distributed throughout the gene, splice-region mutations and large deletions were detected in 6 and 7 families, respectively. In 25 of 51 (49%) patients without hematologic malignancy, somatic mutations were detected in at least 1 of the clonal hematopoiesis of indeterminate potential (CHIP) genes or acute myeloid leukemia (AML) driver genes. BCOR was the most frequently mutated gene (in 9 patients), and multiple BCOR mutations were identified in 4 patients. Mutations in 6 other CHIP- or AML-driver genes (TET2, DNMT3A, KRAS, LRP1B, IDH1, and KMT2C) were also found in ≥2 patients without hematologic malignancy. Moreover, 3 unrelated patients (1 with myeloid malignancy) carried somatic mutations in NFE2, which regulates erythroid and megakaryocytic differentiation. Sequential sequencing data from 19 patients demonstrated dynamic changes of somatic mutations over time, and stable clones were more frequently found in older adult patients. In summary, there are diverse types of germline RUNX1 mutations and high frequency of somatic mutations related to clonal hematopoiesis in patients with FPDMM. Monitoring changes in somatic mutations and clinical manifestations prospectively may reveal mechanisms for malignant progression and inform clinical management. This trial was registered at www.clinicaltrials.gov as #NCT03854318.

Linking GATA2 to myeloid dysplasia and complex cytogenetics in adult myelodysplastic neoplasm and acute myeloid leukemia

ABSTRACT

GATA binding protein 2 (GATA2) is a conserved zinc finger transcription factor that regulates the emergence and maintenance of complex genetic programs driving development and function of hematopoietic stem and progenitor cells (HSPCs). Patients born with monoallelic GATA2 mutations develop myelodysplastic neoplasm (MDS) and acute myeloid leukemia (AML), whereas acquired GATA2 mutations are reported in 3% to 5% of sporadic AML cases. The mechanisms by which aberrant GATA2 activity promotes MDS and AML are incompletely understood. Efforts to understand GATA2 in basic biology and disease will be facilitated by the development of broadly efficacious antibodies recognizing physiologic levels of GATA2 in diverse tissue types and assays. Here, we purified a polyclonal anti-GATA2 antibody and generated multiple highly specific anti-GATA2 monoclonal antibodies, optimized them for immunohistochemistry on patient bone marrow bioosy samples, and analyzed GATA2 expression in adults with healthy bone marrow, MDS, and acute leukemia. In healthy bone marrow, GATA2 was detected in mast cells, subsets of CD34+ HSPCs, E-cadherin-positive erythroid progenitors, and megakaryocytes. In MDS, GATA2 expression correlates with bone marrow blast percentage, positively correlates with myeloid dysplasia and complex cytogenetics, and is a nonindependent negative predictor of overall survival. In acute leukemia, the percent of GATA2+ blasts closely associates with myeloid lineage, whereas a subset of lymphoblastic and undifferentiated leukemias with myeloid features also express GATA2. However, the percent of GATA2+ blasts in AML is highly variable. Elevated GATA2 expression in AML blasts correlates with peripheral neutropenia and complex AML cytogenetics but, unlike in MDS, does not predict survival.

GATA1 in Normal and Pathologic Megakaryopoiesis and Platelet Development

GATA1 is a highly conserved hematopoietic transcription factor (TF), essential for normal erythropoiesis and megakaryopoiesis, that encodes a full-length, predominant isoform and an amino (N) terminus-truncated isoform GATA1s. It is consistently expressed throughout megakaryocyte development and interacts with its target genes either independently or in association with binding partners such as FOG1 (friend of GATA1). While the N-terminus and zinc finger have classically been demonstrated to be necessary for the normal regulation of platelet-specific genes, murine models, cell-line studies, and human case reports indicate that the carboxy-terminal activation domain and zinc finger also play key roles in precisely controlling megakaryocyte growth, proliferation, and maturation. Murine models have shown that disruptions to GATA1 increase the proliferation of immature megakaryocytes with abnormal architecture and impaired terminal differentiation into platelets. In humans, germline GATA1 mutations result in variable cytopenias, including macrothrombocytopenia with abnormal platelet aggregation and excessive bleeding tendencies, while acquired GATA1s mutations in individuals with trisomy 21 (T21) result in transient abnormal myelopoiesis (TAM) and myeloid leukemia of Down syndrome (ML-DS) arising from a megakaryocyte-erythroid progenitor (MEP). Taken together, GATA1 plays a key role in regulating megakaryocyte differentiation, maturation, and proliferative capacity. As sequencing and proteomic technologies expand, additional GATA1 mutations and regulatory mechanisms contributing to human diseases of megakaryocytes and platelets are likely to be revealed.

ABCC4 impacts megakaryopoiesis and protects megakaryocytes against 6-mercaptopurine induced cytotoxicity

The role of ABCC4, an ATP-binding cassette transporter, in the process of platelet formation, megakaryopoiesis, is unknown. Here, we show that ABCC4 is highly expressed in megakaryocytes (MKs). Mining of public genomic data (ATAC-seq and genome wide chromatin interactions, Hi-C) revealed that key megakaryopoiesis transcription factors (TFs) interacted with ABCC4 regulatory elements and likely accounted for high ABCC4 expression in MKs. Importantly these genomic interactions for ABCC4 ranked higher than for genes with known roles in megakaryopoiesis suggesting a role for ABCC4 in megakaryopoiesis. We then demonstrate that ABCC4 is required for optimal platelet formation as in vitro differentiation of fetal liver derived MKs from Abcc4-/- mice exhibited impaired proplatelet formation and polyploidization, features required for optimal megakaryopoiesis. Likewise, a human megakaryoblastic cell line, MEG-01 showed that acute ABCC4 inhibition markedly suppressed key processes in megakaryopoiesis and that these effects were related to reduced cAMP export and enhanced dissociation of a negative regulator of megakaryopoiesis, protein kinase A (PKA) from ABCC4. PKA activity concomitantly increased after ABCC4 inhibition which was coupled with significantly reduced GATA-1 expression, a TF needed for optimal megakaryopoiesis. Further, ABCC4 protected MKs from 6-mercaptopurine (6-MP) as Abcc4-/- mice show a profound reduction in MKs after 6-MP treatment. In total, our studies show that ABCC4 not only protects the MKs but is also required for maximal platelet production from MKs, suggesting modulation of ABCC4 function might be a potential therapeutic strategy to regulate platelet production.

Transcription Factor TAL1 in Erythropoiesis

Lineage-specific transcription factors (TFs) regulate differentiation of hematopoietic stem cells (HSCs). They are decisive for the establishment and maintenance of lineage-specific gene expression programs during hematopoiesis. For this they create a regulatory network between TFs, epigenetic cofactors, and microRNAs. They activate cell-type specific genes and repress competing gene expression programs. Disturbance of this process leads to impaired lineage fidelity and diseases of the blood system. The TF T-cell acute leukemia 1 (TAL1) is central for erythroid differentiation and contributes to the formation of distinct gene regulatory complexes in progenitor cells and erythroid cells. A TAL1/E47 heterodimer binds to DNA with the TFs GATA-binding factor 1 and 2 (GATA1/2), the cofactors LIM domain only 1 and 2 (LMO1/2), and LIM domain-binding protein 1 (LDB1) to form a core TAL1 complex. Furthermore, cell-type-dependent interactions of TAL1 with other TFs such as with runt-related transcription factor 1 (RUNX1) and Kruppel-like factor 1 (KLF1) are established. Moreover, TAL1 activity is regulated by the formation of TAL1 isoforms, posttranslational modifications (PTMs), and microRNAs. Here, we describe the function of TAL1 in normal hematopoiesis with a focus on erythropoiesis.

Wnt/β-Catenin-Signaling Modulates Megakaryopoiesis at the Megakaryocyte-Erythrocyte Progenitor Stage in the Hematopoietic System

The bone marrow (BM) hematopoietic system (HS) gives rise to blood cells originating from hematopoietic stem cells (HSCs), including megakaryocytes (MKs) and red blood cells (erythrocytes; RBCs). Many steps of the cell-fate decision remain to be elucidated, being important for cancer treatment. To explore the role of Wnt/β-catenin for MK and RBC differentiation, we activated β-catenin signaling in platelet-derived growth factor b (Pdgfb)-expressing cells of the HS using a Cre-lox approach (Ctnnb1BM-GOF). FACS analysis revealed that Pdgfb is mainly expressed by megakaryocytic progenitors (MKPs), MKs and platelets. Recombination resulted in a lethal phenotype in mutants (Ctnnb1BM-GOFwt/fl, Ctnnb1BM-GOFfl/fl) 3 weeks after tamoxifen injection, showing an increase in MKs in the BM and spleen, but no pronounced anemia despite reduced erythrocyte counts. BM transplantation (BMT) of Ctnnb1BM-GOF BM into lethally irradiated wildtype recipients (BMT-Ctnnb1BM-GOF) confirmed the megakaryocytic, but not the lethal phenotype. CFU-MK assays in vitro with BM cells of Ctnnb1BM-GOF mice supported MK skewing at the expense of erythroid colonies. Molecularly, the runt-related transcription factor 1 (RUNX1) mRNA, known to suppress erythropoiesis, was upregulated in Ctnnb1BM-GOF BM cells. In conclusion, β-catenin activation plays a key role in cell-fate decision favoring MK development at the expense of erythroid production.

Tumor cell-released kynurenine biases MEP differentiation into megakaryocytes in individuals with cancer by activating AhR-RUNX1

Tumor-derived factors are thought to regulate thrombocytosis and erythrocytopenia in individuals with cancer; however, such factors have not yet been identified. Here we show that tumor cell-released kynurenine (Kyn) biases megakaryocytic-erythroid progenitor cell (MEP) differentiation into megakaryocytes in individuals with cancer by activating the aryl hydrocarbon receptor-Runt-related transcription factor 1 (AhR-RUNX1) axis. During tumor growth, large amounts of Kyn from tumor cells are released into the periphery, where they are taken up by MEPs via the transporter SLC7A8. In the cytosol, Kyn binds to and activates AhR, leading to its translocation into the nucleus where AhR transactivates RUNX1, thus regulating MEP differentiation into megakaryocytes. In addition, activated AhR upregulates SLC7A8 in MEPs to induce positive feedback. Importantly, Kyn-AhR-RUNX1-regulated MEP differentiation was demonstrated in both humanized mice and individuals with cancer, providing potential strategies for the prevention of thrombocytosis and erythrocytopenia.

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.

Transcription factor genetics and biology in predisposition to bone marrow failure and hematological malignancy

Transcription factors (TFs) play a critical role as key mediators of a multitude of developmental pathways, with highly regulated and tightly organized networks crucial for determining both the timing and pattern of tissue development. TFs can act as master regulators of both primitive and definitive hematopoiesis, tightly controlling the behavior of hematopoietic stem and progenitor cells (HSPCs). These networks control the functional regulation of HSPCs including self-renewal, proliferation, and differentiation dynamics, which are essential to normal hematopoiesis. Defining the key players and dynamics of these hematopoietic transcriptional networks is essential to understanding both normal hematopoiesis and how genetic aberrations in TFs and their networks can predispose to hematopoietic disease including bone marrow failure (BMF) and hematological malignancy (HM). Despite their multifaceted and complex involvement in hematological development, advances in genetic screening along with elegant multi-omics and model system studies are shedding light on how hematopoietic TFs interact and network to achieve normal cell fates and their role in disease etiology. This review focuses on TFs which predispose to BMF and HM, identifies potential novel candidate predisposing TF genes, and examines putative biological mechanisms leading to these phenotypes. A better understanding of the genetics and molecular biology of hematopoietic TFs, as well as identifying novel genes and genetic variants predisposing to BMF and HM, will accelerate the development of preventative strategies, improve clinical management and counseling, and help define targeted treatments for these diseases.

PIK3R3 is a candidate regulator of platelet count in people of Bangladeshi ancestry

Background

Blood platelets are mediators of atherothrombotic disease and are regulated by complex sets of genes. Association studies in European ancestry populations have already detected informative platelet regulatory loci. Studies in other ancestries can potentially reveal new associations because of different allele frequencies, linkage structures, and variant effects.

Objectives

To reveal new regulatory genes for platelet count (PLT).

Methods

Genome-wide association studies (GWAS) were performed in 20,218 Bangladeshi and 9198 Pakistani individuals from the Genes & Health study. Loci significantly associated with PLT underwent fine-mapping to identify candidate genes.

Results

Of 1588 significantly associated variants (P < 5 × 10-8) at 20 loci in the Bangladeshi analysis, most replicated findings in prior transancestry GWAS and in the Pakistani analysis. However, the Bangladeshi locus defined by rs946528 (chr1:46019890) did not associate with PLT in the Pakistani analysis but was in the same linkage disequilibrium block (r2 ≥ 0.5) as PLT-associated variants in prior East Asian GWAS. The single independent association signal was refined to a 95% credible set of 343 variants spanning 8 coding genes. Functional annotation, mapping to megakaryocyte regulatory regions, and colocalization with blood expression quantitative trait loci identified the likely mediator of the PLT phenotype to be PIK3R3 encoding a regulator of phosphoinositol 3-kinase (PI3K).

Conclusion

Abnormal PI3K activity in the vessel wall is already implicated in the pathogenesis of atherothrombosis. Our identification of a new association between PIK3R3 and PLT provides further mechanistic insights into the contribution of the PI3K pathway to platelet biology.

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

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

Zebrafish ELL-associated factors Eaf1/2 modulate erythropoiesis via regulating gata1a expression and WNT signaling to facilitate hypoxia tolerance

EAF1 and EAF2, the eleven-nineteen lysine-rich leukemia (ELL)-associated factors which can assemble to the super elongation complex (AFF1/4, AF9/ENL, ELL, and P-TEFb), are reported to participate in RNA polymerase II to actively regulate a variety of biological processes, including leukemia and embryogenesis, but whether and how EAF1/2 function in hematopoietic system related hypoxia tolerance during embryogenesis remains unclear. Here, we unveiled that deletion of EAF1/2 (eaf1^-/- and eaf2^-/-) caused reduction in hypoxia tolerance in zebrafish, leading to reduced erythropoiesis during hematopoietic processes. Meanwhile, eaf1^-/- and eaf2^-/- mutants showed significant reduction in the expression of key transcriptional regulators scl, lmo2, and gata1a in erythropoiesis at both 24 h post fertilization (hpf) and 72 hpf, with gata1a downregulated while scl and lmo2 upregulated at 14 hpf. Mechanistically, eaf1^-/- and eaf2^-/- mutants exhibited significant changes in the expression of epigenetic modified histones, with a significant increase in the binding enrichment of modified histone H3K27me3 in gata1a promoter rather than scl and lmo2 promoters. Additionally, eaf1^-/- and eaf2^-/- mutants exhibited a dynamic expression of canonical WNT/β-catenin signaling during erythropoiesis, with significant reduction in p-β-Catenin level and in the binding enrichment of both scl and lmo2 promoters with the WNT transcriptional factor TCF4 at 24 hpf. These findings demonstrate an important role of Eaf1/2 in erythropoiesis in zebrafish and may have shed some light on regeneration medicine for anemia and related diseases and on molecular basis for fish economic or productive traits, such as growth, disease resistance, hypoxia tolerance, and so on.

In Vitro Human Haematopoietic Stem Cell Expansion and Differentiation

The haematopoietic system plays an essential role in our health and survival. It is comprised of a range of mature blood and immune cell types, including oxygen-carrying erythrocytes, platelet-producing megakaryocytes and infection-fighting myeloid and lymphoid cells. Self-renewing multipotent haematopoietic stem cells (HSCs) and a range of intermediate haematopoietic progenitor cell types differentiate into these mature cell types to continuously support haematopoietic system homeostasis throughout life. This process of haematopoiesis is tightly regulated in vivo and primarily takes place in the bone marrow. Over the years, a range of in vitro culture systems have been developed, either to expand haematopoietic stem and progenitor cells or to differentiate them into the various haematopoietic lineages, based on the use of recombinant cytokines, co-culture systems and/or small molecules. These approaches provide important tractable models to study human haematopoiesis in vitro. Additionally, haematopoietic cell culture systems are being developed and clinical tested as a source of cell products for transplantation and transfusion medicine. This review discusses the in vitro culture protocols for human HSC expansion and differentiation, and summarises the key factors involved in these biological processes.

RUNX1 isoform disequilibrium promotes the development of trisomy 21-associated myeloid leukemia

Gain of chromosome 21 (Hsa21) is among the most frequent aneuploidies in leukemia. However, it remains unclear how partial or complete amplifications of Hsa21 promote leukemogenesis and why children with Down syndrome (DS) (ie, trisomy 21) are particularly at risk of leukemia development. Here, we propose that RUNX1 isoform disequilibrium with RUNX1A bias is key to DS-associated myeloid leukemia (ML-DS). Starting with Hsa21-focused CRISPR-CRISPR-associated protein 9 screens, we uncovered a strong and specific RUNX1 dependency in ML-DS cells. Expression of the RUNX1A isoform is elevated in patients with ML-DS, and mechanistic studies using murine ML-DS models and patient-derived xenografts revealed that excess RUNX1A synergizes with the pathognomonic Gata1s mutation during leukemogenesis by displacing RUNX1C from its endogenous binding sites and inducing oncogenic programs in complex with the MYC cofactor MAX. These effects were reversed by restoring the RUNX1A:RUNX1C equilibrium in patient-derived xenografts in vitro and in vivo. Moreover, pharmacological interference with MYC:MAX dimerization using MYCi361 exerted strong antileukemic effects. Thus, our study highlights the importance of alternative splicing in leukemogenesis, even on a background of aneuploidy, and paves the way for the development of specific and targeted therapies for ML-DS, as well as for other leukemias with Hsa21 aneuploidy or RUNX1 isoform disequilibrium.

Platelet transcriptome analysis in patients with germline RUNX1 mutations

BACKGROUND

Germline mutations in RUNX1 can cause a familial platelet disorder that may lead to acute myeloid leukemia, an autosomal dominant disorder characterized by moderate thrombocytopenia, platelet dysfunction, and a high risk of developing acute myeloid leukemia or myelodysplastic syndrome. Discerning the pathogenicity of novel RUNX1 variants is critical for patient management.

OBJECTIVES

To extend the characterization of RUNX1 variants and evaluate their effects by transcriptome analysis.

METHODS

Three unrelated patients with long-standing thrombocytopenia carrying heterozygous RUNX1 variants were included: P1, who is a subject with recent development of myelodysplastic syndrome, with c.802 C>T[p.Gln268∗] de novo; P2 with c.586A>G[p.Thr196Ala], a variant that segregates with thrombocytopenia and myeloid neoplasia in the family; and P3 with c.476A>G[p.Asn159Ser], which did not segregate with thrombocytopenia or neoplasia. Baseline platelet evaluations were performed. Ultrapure platelets were prepared for platelet transcriptome analysis.

RESULTS

In P1 and P2, but not in P3, transcriptome analysis confirmed aberrant expression of genes recognized as RUNX1 targets. Data allowed grouping patients by distinct gene expression profiles, which were partitioned with clinical parameters. Functional studies and platelet mRNA expression identified alterations in the actin cytoskeleton, downregulation of GFI1B, defective GPVI downstream signaling, and reduction of alpha granule proteins, such as thrombospondin-1, as features likely implicated in thrombocytopenia and platelet dysfunction.

CONCLUSION

Platelet phenotype, familial segregation, and platelet transcriptomics support the pathogenicity of RUNX1 variants p.Gln268∗ and p.Thr196Ala, but not p.Asn159Ser. This study is an additional proof of concept that platelet RNA analysis could be a tool to help classify pathogenic RUNX1 variants and identify novel RUNX1 targets.

Histone deacetylation-regulated cell surface Siglec-7 expression promoted megakaryocytic maturation and enhanced platelet-like particle release

BACKGROUND

Functioning as important hematologic cells for hemostasis, wound healing and immune defense platelets are produced before being released into the blood by cytoplasmic fragmentation at the end of the megakaryocyte (MK) differentiation, during which the involvement of both apoptosis and autophagy has been reported. Inhibitory sialic acid-binding immunoglobulin-like lectin-7 gene (Siglec-7) can be expressed on platelets and induce apoptosis on activation for uncharacterized function.

OBJECTIVE

We aimed to investigate the regulatory mechanism for Siglec-7 activation along MK differentiation and its physiologic role during the MK maturation and platelet formation.

METHODS

By using 2 well-established MK differentiation models (HEL and K562) and human primary CD34⁺ cell, we examined the upregulations of transcript and protein levels of Siglec-7 during MK differentiation, and the effect of Siglec-7 surface presence on MK differentiation and platelet-like particles (PLPs) release.

RESULTS

We show that both transcripts and surface Siglec-7 were elevated during MK differentiation, and the histone deacetylase 1 (HDAC1) acted as a negative regulator for Siglec-7 activation. By increasing Siglec-7 surface expression, we found that increased presence of Siglec-7 not only enhanced MK maturation but also the release of PLPs by activating caspase 3-dependent signaling, as evidenced in the observation of more CD41, polyploidy, and platelet factor 4 transcript formations.

CONCLUSION

In this study, we demonstrated that Siglec-7 activation was subjected to epigenetic regulation, and the resulting induced expression of surface Siglec-7 played an important regulatory role in promoting MK differentiation, maturation, and PLP formation.

Genomic Landscape of Patients with Germline RUNX1 Variants and Familial Platelet Disorder with Myeloid Malignancy

Germline RUNX1 mutations lead to familial platelet disorder with associated myeloid malignancies (FPDMM), which is characterized by thrombocytopenia and a life-long risk (35-45%) of hematological malignancies. We recently launched a longitudinal natural history study for patients with FPDMM at the NIH Clinical Center. Among 29 families with research genomic data, 28 different germline RUNX1 variants were detected. Besides missense mutations enriched in Runt homology domain and loss-of-function mutations distributed throughout the gene, splice-region mutations and large deletions were detected in 6 and 7 families, respectively. In 24 of 54 (44.4%) non-malignant patients, somatic mutations were detected in at least one of the clonal hematopoiesis of indeterminate potential (CHIP) genes or acute myeloid leukemia (AML) driver genes. BCOR was the most frequently mutated gene (in 9 patients), and multiple BCOR mutations were identified in 4 patients. Mutations in 7 other CHIP or AML driver genes ( DNMT3A, TET2, NRAS, SETBP1, SF3B1, KMT2C , and LRP1B ) were also found in more than one non-malignant patient. Moreover, three unrelated patients (one with myeloid malignancy) carried somatic mutations in NFE2 , which regulates erythroid and megakaryocytic differentiation. Sequential sequencing data from 19 patients demonstrated dynamic changes of somatic mutations over time, and stable clones were more frequently found in elderly patients. In summary, there are diverse types of germline RUNX1 mutations and high frequency of somatic mutations related to clonal hematopoiesis in patients with FPDMM. Monitoring dynamic changes of somatic mutations prospectively will benefit patients’ clinical management and reveal mechanisms for progression to myeloid malignancies.

Key Points

Comprehensive genomic profile of patients with FPDMM with germline RUNX1 mutations. Rising clonal hematopoiesis related secondary mutations that may lead to myeloid malignancies.

Weighted gene co-expression network analysis identifies key hub genes and pathways in acute myeloid leukemia

Introduction: Acute myeloid leukemia (AML) is the most common type of leukemia in adults. However, there is a gap in understanding the molecular basis of the disease, partly because key genes associated with AML have not been extensively explored. In the current study, we aimed to identify genes that have strong association with AML based on a cross-species integrative approach. Methods: We used Weighted Gene Co-Expression Network Analysis (WGCNA) to identify co-expressed gene modules significantly correlated with human AML, and further selected the genes exhibiting a significant difference in expression between AML and healthy mouse. Protein-protein interactions, transcription factors, gene function, genetic regulation, and coding sequence variants were integrated to identify key hub genes in AML. Results: The cross-species approach identified a total of 412 genes associated with both human and mouse AML. Enrichment analysis confirmed an association of these genes with hematopoietic and immune-related functions, phenotypes, processes, and pathways. Further, the integrated analysis approach identified a set of important module genes including Nfe2, Trim27, Mef2c, Ets1, Tal1, Foxo1, and Gata1 in AML. Six of these genes (except ETS1) showed significant differential expression between human AML and healthy samples in an independent microarray dataset. All of these genes are known to be involved in immune/hematopoietic functions, and in transcriptional regulation. In addition, Nfe2, Trim27, Mef2c, and Ets1 harbor coding sequence variants, whereas Nfe2 and Trim27 are cis-regulated, making them attractive candidates for validation. Furthermore, subtype-specific analysis of the hub genes in human AML indicated high expression of NFE2 across all the subtypes (M0 through M7) and enriched expression of ETS1, LEF1, GATA1, and TAL1 in M6 and M7 subtypes. A significant correlation between methylation status and expression level was observed for most of these genes in AML patients. Conclusion: Findings from the current study highlight the importance of our cross-species approach in the identification of multiple key candidate genes in AML, which can be further studied to explore their detailed role in leukemia/AML.

FLI1 accelerates leukemogenesis through transcriptional regulation of pyruvate kinase-L/R and other glycolytic genes

In cancer cells, multiple oncogenes and tumor suppressors control glycolysis to sustain rapid proliferation. The ETS-related transcription factor Fli1 plays a critical role in the induction and progression of leukemia, yet, the underlying mechanism of this oncogenic event is still not fully understood. In this study, RNAseq analysis of FLI1-depleted human leukemic cells revealed transcriptional suppression of the PKLR gene and activation of multiple glycolytic genes, such as PKM1/2. Pharmacological inhibition of glycolysis by PKM2 inhibitor, Shikonin, significantly suppressed leukemic cell proliferation. FLI1 directly binds to the PKLR promoter, leading to the suppression of this inhibitor of glycolysis. In accordance, shRNA-mediated depletion of PKLR in leukemic HEL cells expressing high levels of FLI1 accelerated leukemia proliferation, pointing for the first time to its tumor suppressor function. PKLR knockdown also led to downregulation of the erythroid markers EPOR, HBA1, and HBA2 and suppression of erythroid differentiation. Interestingly, silencing of PKLR in HEL cells significantly increased FLI1 expression, which was associated with faster proliferation in culture. In FLI1-expressing leukemic cells, lower PKLR expression was associated with higher expression of PKM1 and PKM2, which promote aerobic glycolysis. Finally, injection of pyruvate, a known inhibitor of glycolysis, into leukemia mice significantly suppressed leukemogenesis. These results demonstrate that FLI1 promotes leukemia in part by inducing glycolysis, implicates PKLR in erythroid differentiation, and suggests that targeting glycolysis may be an attractive therapeutic strategy for cancers driven by FLI1 overexpression.

A Novel GATA1 Variant in the C-Terminal Zinc Finger Compared with the Platelet Phenotype of Patients with A Likely Pathogenic Variant in the N-Terminal Zinc Finger

The GATA1 transcription factor is essential for normal erythropoiesis and megakaryocytic differentiation. Germline GATA1 pathogenic variants in the N-terminal zinc finger (N-ZF) are typically associated with X-linked thrombocytopenia, platelet dysfunction, and dyserythropoietic anemia. A few variants in the C-terminal ZF (C-ZF) domain are described with normal platelet count but altered platelet function as the main characteristic. Independently performed molecular genetic analysis identified a novel hemizygous variant (c.865C>T, p.H289Y) in the C-ZF region of GATA1 in a German patient and in a Spanish patient. We characterized the bleeding and platelet phenotype of these patients and compared these findings with the parameters of two German siblings carrying the likely pathogenic variant p.D218N in the GATA1 N-ZF domain. The main difference was profound thrombocytopenia in the brothers carrying the p.D218N variant compared to a normal platelet count in patients carrying the p.H289Y variant; only the Spanish patient occasionally developed mild thrombocytopenia. A functional platelet defect affecting αIIbβ3 integrin activation and α-granule secretion was present in all patients. Additionally, mild anemia, anisocytosis, and poikilocytosis were observed in the patients with the C-ZF variant. Our data support the concept that GATA1 variants located in the different ZF regions can lead to clinically diverse manifestations.

Hirsutine, a novel megakaryopoiesis inducer, promotes thrombopoiesis via MEK/ERK/FOG1/TAL1 signaling

BACKGROUND

Thrombocytopenia (TP) remains a challenge in clinical hematology. TP may have serious consequences, such as recurrent skin and mucosal bleeding and increased risk of intracranial and internal organ hemorrhage. However, effective and safe therapeutic drugs for the long-term management of TP are still lacking.

PURPOSE

This study aimed to identify more effective active compounds for TP therapy.

METHODS

Liquid chromatography-mass spectrometry-nuclear magnetic resonance analysis was used to confirm the medicinal species and chemical structure of Hirsutine (HS). The proliferation of HS was examined by Cell Counting Kit (CCK-8) assay on cells lines. The effect of HS on megakaryocyte differentiation was analyzed by evaluating the expression of CD41, CD42b, and DNA ploidy via flow cytometry (FCM). The morphology of megakaryocytes and intermediate cells was observed using an optical microscope. K562 cells were then stained with Giemsa and benzidine. qRT-PCR was used to examine the mRNA expression of GATA-1, GATA-2, FOG-1, TAL-1, RUNX-1, NF-E2, and KLF-1 in K562 cells. Protein levels of the transcription factors were analyzed by western blotting. An MEK inhibitor was used to verify the relationship between the MEK/ERK signaling pathway and CD41/CD42b (FCM), FOG-1, and TAL-1. The Kunming thrombocytopenia mouse model was established by X-ray irradiation (4 Gy) and used to test HS activity and related hematopoietic organ index in vivo. Finally, computer simulations of molecular docking were used to predict the binding energies between HS-MEK and HS-ERK.

RESULTS

We preliminarily identified HS by screening a plant-sourced compound library for natural compounds with megakaryocytic differentiation and maturation (MKD/MKM)-promoting activity. We found that HS not only enhanced MKD/MKM of K562 and Meg01 cells, but also suppressed the decline of peripheral platelet levels in X-ray-induced myelosuppressive mice. In addition, HS promoted MKD via activation of MEK-ERK-FOG1/TAL1 signaling, which may be the key molecular mechanism of HS action in TP treatment. Molecular docking simulations further verified that HS targets the signaling protein MEK with high-affinity.

CONCLUSION

In this study, we report for the first time that hirsutine boosts MKD/MKM through the MEK/ERK/FOG1/TAL1 signaling pathway and thus represents a promising treatment option for TP.

Leukemia-initiating HSCs in chronic lymphocytic leukemia reveal clonal leukemogenesis and differential drug sensitivity

The existence of "leukemia-initiating cells" (LICs) in chronic lymphocytic leukemia (CLL) remains controversial due to the difficulty in isolating and identifying the tumor-initiating cells. Here, we demonstrate a microchannel electroporation (MEP) microarray that injects RNA-detecting probes into single live cells, allowing the imaging and characterization of heterogeneous LICs by intracellular RNA expression. Using limited-cell FACS sequencing (LC-FACSeq), we can detect and monitor rare live LICs during leukemogenesis and characterize their differential drug sensitivity. Disease-associated mutation accumulation in developing B lymphoid but not myeloid lineage in CLL patient hematopoietic stem cells (CLL-HSCs), and development of independent clonal CLL-like cells in murine patient-derived xenograft models, suggests the existence of CLL LICs. Furthermore, we identify differential protein ubiquitination and unfolding response signatures in GATA2^high CLL-HSCs that exhibit increased sensitivity to lenalidomide and resistance to fludarabine compared to GATA2^lowCLL-HSCs. These results highlight the existence of therapeutically targetable disease precursors in CLL.

Stepwise GATA1 and SMC3 mutations alter megakaryocyte differentiation in a Down syndrome leukemia model

Acute megakaryoblastic leukemia of Down syndrome (DS-AMKL) is a model of clonal evolution from a preleukemic transient myeloproliferative disorder requiring both a trisomy 21 (T21) and a GATA1s mutation to a leukemia driven by additional driver mutations. We modeled the megakaryocyte differentiation defect through stepwise gene editing of GATA1s, SMC3^+/–, and MPL^W515K, providing 20 different T21 or disomy 21 (D21) induced pluripotent stem cell (iPSC) clones. GATA1s profoundly reshaped iPSC-derived hematopoietic architecture with gradual myeloid-to-megakaryocyte shift and megakaryocyte differentiation alteration upon addition of SMC3 and MPL mutations. Transcriptional, chromatin accessibility, and GATA1-binding data showed alteration of essential megakaryocyte differentiation genes, including NFE2 downregulation that was associated with loss of GATA1s binding and functionally involved in megakaryocyte differentiation blockage. T21 enhanced the proliferative phenotype, reproducing the cellular and molecular abnormalities of DS-AMKL. Our study provides an array of human cell–based models revealing individual contributions of different mutations to DS-AMKL differentiation blockage, a major determinant of leukemic progression.

Preliminary Interpretations of Epigenetic Profiling of Cord Blood in Preeclampsia

Preeclampsia (PE) is characterized by new-onset hypertension after 20 weeks of pregnancy and results in high maternal and fetal mortality worldwide. It has been reported that PE is associated with abnormalities in the umbilical cord and cord blood. However, previous studies were focused primarily on the transcriptomics level, while the underlying gene regulatory landscapes are still unclear. Thus, we performed the Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-seq) using the umbilical cord blood samples collected from a patient with superimposed PE and three healthy donors to uncover the chromatin accessibility changes attributed to PE. We have identified genes associated with immunomodulation and hypoxia response that have higher chromatin accessibility close to their transcription start sites. Motif analysis indicated that the GATA family transcription factor binding was enriched in PE and may play an essential regulatory role in the disease progression. Overall, our findings provide an overview of gene regulatory programs and the corresponding downstream pathways associated with PE that may influence the placenta function and fetal growth.

E Protein Transcription Factors as Suppressors of T Lymphocyte Acute Lymphoblastic Leukemia

T Lymphocyte Acute Lymphoblastic Leukemia (ALL) is an aggressive disease arising from transformation of T lymphocytes during their development. The mutation spectrum of T-ALL has revealed critical regulators of the growth and differentiation of normal and leukemic T lymphocytes. Approximately, 60% of T-ALLs show aberrant expression of the hematopoietic stem cell-associated helix-loop-helix transcription factors TAL1 and LYL1. TAL1 and LYL1 function in multiprotein complexes that regulate gene expression in T-ALL but they also antagonize the function of the E protein homodimers that are critical regulators of T cell development. Mice lacking E2A, or ectopically expressing TAL1, LYL1, or other inhibitors of E protein function in T cell progenitors, also succumb to an aggressive T-ALL-like disease highlighting that E proteins promote T cell development and suppress leukemogenesis. In this review, we discuss the role of E2A in T cell development and how alterations in E protein function underlie leukemogenesis. We focus on the role of TAL1 and LYL1 and the genes that are dysregulated in E2a^-/- T cell progenitors that contribute to human T-ALL. These studies reveal novel mechanisms of transformation and provide insights into potential therapeutic targets for intervention in this disease.

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

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

The megakaryocytic transcription factor ARID3A suppresses leukemia pathogenesis

Given the plasticity of hematopoietic stem and progenitor cells, multiple routes of differentiation must be blocked in the the pathogenesis of acute myeloid leukemia, the molecular basis of which is incompletely understood. We report that posttranscriptional repression of the transcription factor ARID3A by miR-125b is a key event in the pathogenesis of acute megakaryoblastic leukemia (AMKL). AMKL is frequently associated with trisomy 21 and GATA1 mutations (GATA1s), and children with Down syndrome are at a high risk of developing the disease. The results of our study showed that chromosome 21-encoded miR-125b synergizes with Gata1s to drive leukemogenesis in this context. Leveraging forward and reverse genetics, we uncovered Arid3a as the main miR-125b target behind this synergy. We demonstrated that, during normal hematopoiesis, this transcription factor promotes megakaryocytic differentiation in concert with GATA1 and mediates TGFβ-induced apoptosis and cell cycle arrest in complex with SMAD2/3. Although Gata1s mutations perturb erythroid differentiation and induce hyperproliferation of megakaryocytic progenitors, intact ARID3A expression assures their megakaryocytic differentiation and growth restriction. Upon knockdown, these tumor suppressive functions are revoked, causing a blockade of dual megakaryocytic/erythroid differentiation and subsequently of AMKL. Inversely, restoring ARID3A expression relieves the arrest of megakaryocytic differentiation in AMKL patient-derived xenografts. This work illustrates how mutations in lineage-determining transcription factors and perturbation of posttranscriptional gene regulation can interact to block multiple routes of hematopoietic differentiation and cause leukemia. In AMKL, surmounting this differentiation blockade through restoration of the tumor suppressor ARID3A represents a promising strategy for treating this lethal pediatric disease.

Domain adaptive neural networks improve cross-species prediction of transcription factor binding

The intrinsic DNA sequence preferences and cell type–specific cooperative partners of transcription factors (TFs) are typically highly conserved. Hence, despite the rapid evolutionary turnover of individual TF binding sites, predictive sequence models of cell type–specific genomic occupancy of a TF in one species should generalize to closely matched cell types in a related species. To assess the viability of cross-species TF binding prediction, we train neural networks to discriminate ChIP-seq peak locations from genomic background and evaluate their performance within and across species. Cross-species predictive performance is consistently worse than within-species performance, which we show is caused in part by species-specific repeats. To account for this domain shift, we use an augmented network architecture to automatically discourage learning of training species–specific sequence features. This domain adaptation approach corrects for prediction errors on species-specific repeats and improves overall cross-species model performance. Our results show that cross-species TF binding prediction is feasible when models account for domain shifts driven by species-specific repeats.

Sorting nexin 24 is required for α-granule biogenesis and cargo delivery in megakaryocytes

Germline defects affecting the DNA-binding domain of the transcription factor FLI1 are associated with a bleeding disorder that is characterized by the presence of large, fused α-granules in platelets. We investigated whether the genes showing abnormal expression in FLI1-deficient platelets could be involved in platelet α-granule biogenesis by undertaking transcriptome analysis of control platelets and platelets harboring a DNA-binding variant of FLI1. Our analysis identified 2,276 transcripts that were differentially expressed in FLI1-deficient platelets. Functional annotation clustering of the coding transcripts revealed significant enrichment for gene annotations relating to protein transport, and identified Sorting nexin 24 (SNX24) as a candidate for further investigation. Using an induced pluripotent stem cell-derived megakaryocyte model, SNX24 expression was found to be increased during the early stages of megakaryocyte differentiation and downregulated during proplatelet formation, indicating tight regulatory control during megakaryopoiesis. CRISPR-Cas9 mediated knockout (KO) of SNX24 led to decreased expression of immature megakaryocyte markers, CD41 and CD61, and increased expression of the mature megakaryocyte marker CD42b (P=0.0001), without affecting megakaryocyte polyploidisation, or proplatelet formation. Electron microscopic analysis revealed an increase in empty membrane-bound organelles in SNX24 KO megakaryocytes, a reduction in α-granules and an absence of immature and mature multivesicular bodies, consistent with a defect in the intermediate stage of α-granule maturation. Co-localization studies showed that SNX24 associates with each compartment of α-granule maturation. Reduced expression of CD62P and VWF was observed in SNX24 KO megakaryocytes. We conclude that SNX24 is required for α-granule biogenesis and intracellular trafficking of α-granule cargo within megakaryocytes.

G protein-coupled receptor kinase 5 regulates thrombin signaling in platelets via PAR-1

Platelet reactivity via the PAR-1 thrombin receptor is mediated by GRK5.

Platelet GRK5 is associated with thrombus formation in humans and mice.

The interindividual variation in the functional response of platelets to activation by agonists is heritable. Genome-wide association studies (GWASs) of quantitative measures of platelet function have identified fewer than 20 distinctly associated variants, some with unknown mechanisms. Here, we report GWASs of pathway-specific functional responses to agonism by adenosine 5′-diphosphate, a glycoprotein VI–specific collagen mimetic, and thrombin receptor-agonist peptides, each specific to 1 of the G protein–coupled receptors PAR-1 and PAR-4, in subsets of 1562 individuals. We identified an association (P = 2.75 × 10⁻⁴⁰) between a common intronic variant, rs10886430, in the G protein–coupled receptor kinase 5 gene (GRK5) and the sensitivity of platelets to activate through PAR-1. The variant resides in a megakaryocyte-specific enhancer that is bound by the transcription factors GATA1 and MEIS1. The minor allele (G) is associated with fewer GRK5 transcripts in platelets and the greater sensitivity of platelets to activate through PAR-1. We show that thrombin-mediated activation of human platelets causes binding of GRK5 to PAR-1 and that deletion of the mouse homolog Grk5 enhances thrombin-induced platelet activation sensitivity and increases platelet accumulation at the site of vascular injury. This corroborates evidence that the human G allele of rs10886430 is associated with a greater risk for cardiovascular disease. In summary, by combining the results of pathway-specific GWASs and expression quantitative trait locus studies in humans with the results from platelet function studies in Grk5^−/− mice, we obtain evidence that GRK5 regulates the human platelet response to thrombin via the PAR-1 pathway.

Megakaryopoiesis and Platelet Biology: Roles of Transcription Factors and Emerging Clinical Implications

Platelets play a critical role in hemostasis and thrombus formation. Platelets are small, anucleate, and short-lived blood cells that are produced by the large, polyploid, and hematopoietic stem cell (HSC)-derived megakaryocytes in bone marrow. Approximately 3000 platelets are released from one megakaryocyte, and thus, it is important to understand the physiologically relevant mechanism of development of mature megakaryocytes. Many genes, including several key transcription factors, have been shown to be crucial for platelet biogenesis. Mutations in these genes can perturb megakaryopoiesis or thrombopoiesis, resulting in thrombocytopenia. Metabolic changes owing to inflammation, ageing, or diseases such as cancer, in which platelets play crucial roles in disease development, can also affect platelet biogenesis. In this review, I describe the characteristics of platelets and megakaryocytes in terms of their differentiation processes. The role of several critical transcription factors have been discussed to better understand the changes in platelet biogenesis that occur during disease or ageing.

GATA1 pathogenic variants disrupt MYH10 silencing during megakaryopoiesis

BACKGROUND

GATA1 is an essential transcription factor for both polyploidization and megakaryocyte (MK) differentiation. The polyploidization defect observed in GATA1 variant carriers is not well understood.

OBJECTIVE

To extensively phenotype two pedigrees displaying different variants in the GATA1 gene and determine if GATA1 controls MYH10 expression levels, a key modulator of MK polyploidization.

METHOD

A total of 146 unrelated propositi with constitutional thrombocytopenia were screened on a multigene panel. We described the genotype-phenotype correlation in GATA1 variant carriers and investigated the effect of these novel variants on MYH10 transcription using luciferase constructs.

RESULTS

The clinical profile associated with the p.L268M variant localized in the C terminal zinc finger was unusual in that the patient displayed bleeding and severe platelet aggregation defects without early-onset thrombocytopenia. p.N206I localized in the N terminal zinc finger was associated, on the other hand, with severe thrombocytopenia (15G/L) in early life. High MYH10 levels were evidenced in platelets of GATA1 variant carriers. Analysis of MKs anti-GATA1 chromatin immunoprecipitation-sequencing data revealed two GATA1 binding sites, located in the 3' untranslated region and in intron 8 of the MYH10 gene. Luciferase reporter assays showed their respective role in the regulation of MYH10 gene expression. Both GATA1 variants significantly alter intron 8 driven MYH10 transcription.

CONCLUSION

The discovery of an association between MYH10 and GATA1 is a novel one. Overall, this study suggests that impaired MYH10 silencing via an intronic regulatory element is the most likely cause of GATA1-related polyploidization defect.

Specific proteome changes in platelets from individuals with GATA1-, GFI1B-, and RUNX1-linked bleeding disorders

Mutations in the transcription factors GATA binding factor 1 (GATA1), growth factor independence 1B (GFI1B), and Runt-related transcription factor 1 (RUNX1) cause familial platelet and bleeding disorders. Mutant platelets exhibit common abnormalities including an α-granule reduction resulting in a grayish appearance in blood smears. This suggests that similar pathways are deregulated by different transcription factor mutations. To identify common factors, full platelet proteomes from 11 individuals with mutant GATA1R216Q, GFI1BQ287*, RUNX1Q154Rfs, or RUNX1TD2-6 and 28 healthy controls were examined by label-free quantitative mass spectrometry. In total, 2875 platelet proteins were reliably quantified. Clustering analysis of more than 300 differentially expressed proteins revealed profound differences between cases and controls. Among cases, 44 of 143 significantly downregulated proteins were assigned to platelet function, hemostasis, and granule biology, in line with platelet dysfunction and bleedings. Remarkably, none of these proteins were significantly diminished in all affected cases. Similarly, no proteins were commonly overrepresented in all affected cases compared with controls. These data indicate that the studied transcription factor mutations alter platelet proteomes in distinct largely nonoverlapping manners. This work provides the quantitative landscape of proteins that affect platelet function when deregulated by mutated transcription factors in inherited bleeding disorders.

RUNX-1 haploinsufficiency causes a marked deficiency of megakaryocyte-biased hematopoietic progenitor cells

Patients with familial platelet disorder with a predisposition to myeloid malignancy (FPDMM) harbor germline monoallelic mutations in a key hematopoietic transcription factor, RUNX-1. Previous studies of FPDMM have focused on megakaryocyte (Mk) differentiation and platelet production and signaling. However, the effects of RUNX-1 haploinsufficiency on hematopoietic progenitor cells (HPCs) and subsequent megakaryopoiesis remains incomplete. We studied induced pluripotent stem cell (iPSC)-derived HPCs (iHPCs) and Mks (iMks) from both patient-derived lines and a wild-type (WT) line modified to be RUNX-1 haploinsufficient (RUNX-1+/-), each compared with their isogenic WT control. All RUNX-1+/- lines showed decreased iMk yield and depletion of an Mk-biased iHPC subpopulation. To investigate global and local gene expression changes underlying this iHPC shift, single-cell RNA sequencing was performed on sorted FPDMM and control iHPCs. We defined several cell subpopulations in the Mk-biased iHPCs. Analyses of gene sets upregulated in FPDMM iHPCs indicated enrichment for response to stress, regulation of signal transduction, and immune signaling-related gene sets. Immunoblot analyses in FPDMM iMks were consistent with these findings, but also identified augmented baseline c-Jun N-terminal kinase (JNK) phosphorylation, known to be activated by transforming growth factor-β1 (TGF-β1) and cellular stressors. These findings were confirmed in adult human CD34+-derived stem and progenitor cells (HSPCs) transduced with lentiviral RUNX1 short hairpin RNA to mimic RUNX-1+/-. In both iHPCs and CD34+-derived HSPCs, targeted inhibitors of JNK and TGF-β1 pathways corrected the megakaryopoietic defect. We propose that such intervention may correct the thrombocytopenia in patients with FPDMM.

Enhancer-gene rewiring in the pathogenesis of Quebec platelet disorder

Quebec platelet disorder (QPD) is an autosomal dominant bleeding disorder with a unique, platelet-dependent, gain-of-function defect in fibrinolysis, without systemic fibrinolysis. The hallmark feature of QPD is a >100-fold overexpression of PLAU, specifically in megakaryocytes. This overexpression leads to a >100-fold increase in platelet stores of urokinase plasminogen activator (PLAU/uPA); subsequent plasmin-mediated degradation of diverse α-granule proteins; and platelet-dependent, accelerated fibrinolysis. The causative mutation is a 78-kb tandem duplication of PLAU. How this duplication causes megakaryocyte-specific PLAU overexpression is unknown. To investigate the mechanism that causes QPD, we used epigenomic profiling, comparative genomics, and chromatin conformation capture approaches to study PLAU regulation in cultured megakaryocytes from participants with QPD and unaffected controls. QPD duplication led to ectopic interactions between PLAU and a conserved megakaryocyte enhancer found within the same topologically associating domain (TAD). Our results support a unique disease mechanism whereby the reorganization of sub-TAD genome architecture results in a dramatic, cell-type-specific blood disorder phenotype.

In vivo CD8+ T cell CRISPR screening reveals control by Fli1 in infection and cancer

Improving effector activity of antigen-specific T cells is a major goal in cancer immunotherapy. Despite the identification of several effector T cell (TEFF)-driving transcription factors (TFs), the transcriptional coordination of TEFF biology remains poorly understood. We developed an in vivo T cell CRISPR screening platform and identified a key mechanism restraining TEFF biology through the ETS family TF, Fli1. Genetic deletion of Fli1 enhanced TEFF responses without compromising memory or exhaustion precursors. Fli1 restrained TEFF lineage differentiation by binding to cis-regulatory elements of effector-associated genes. Loss of Fli1 increased chromatin accessibility at ETS:RUNX motifs, allowing more efficient Runx3-driven TEFF biology. CD8+ T cells lacking Fli1 provided substantially better protection against multiple infections and tumors. These data indicate that Fli1 safeguards the developing CD8+ T cell transcriptional landscape from excessive ETS:RUNX-driven TEFF cell differentiation. Moreover, genetic deletion of Fli1 improves TEFF differentiation and protective immunity in infections and cancer.

CBFβ-SMMHC Affects Genome-wide Polycomb Repressive Complex 1 Activity in Acute Myeloid Leukemia

Mutations and deletions of polycomb repressive complex (PRC) components are increasingly recognized to affect tumor biology in a range of cancers. However, little is known about how genetic alterations of PRC-interacting molecules such as the core binding factor (CBF) complex influence polycomb activity. We report that the acute myeloid leukemia (AML)-associated CBFβ-SMMHC fusion oncoprotein physically interacts with the PRC1 complex and that these factors co-localize across the AML genome in an apparently PRC2-independent manner. Depletion of CBFβ-SMMHC caused substantial increases in genome-wide PRC1 binding and marked changes in the association between PRC1 and the CBF DNA-binding subunit RUNX1. PRC1 was more likely to be associated with actively transcribed genes in CBFβ-SMMHC-expressing cells. CBFβ-SMMHC depletion had heterogeneous effects on gene expression, including significant reductions in transcription of ribosomal loci occupied by PRC1. Our results provide evidence that CBFβ-SMMHC markedly and diversely affects polycomb recruitment and transcriptional regulation across the AML genome.

Molecular Co-occupancy Identifies Transcription Factor Binding Cooperativity In Vivo

Gene activation requires the cooperative activity of multiple transcription factors at cis-regulatory elements (CREs). Yet, most transcription factors have short residence time, questioning the requirement of their physical co-occupancy on DNA to achieve cooperativity. Here, we present a DNA footprinting method that detects individual molecular interactions of transcription factors and nucleosomes with DNA in vivo. We apply this strategy to quantify the simultaneous binding of multiple transcription factors on single DNA molecules at mouse CREs. Analysis of the binary occupancy patterns at thousands of motif combinations reveals that high DNA co-occupancy occurs for most types of transcription factors, in the absence of direct physical interaction, at sites of competition with nucleosomes. Perturbation of pairwise interactions demonstrates the function of molecular co-occupancy in binding cooperativity. Our results reveal the interactions regulating CREs at molecular resolution and identify DNA co-occupancy as a widespread cooperativity mechanism used by transcription factors to remodel chromatin.

Molecular Fingerprint and Developmental Regulation of the Tegmental GABAergic and Glutamatergic Neurons Derived from the Anterior Hindbrain

Tegmental nuclei in the ventral midbrain and anterior hindbrain control motivated behavior, mood, memory, and movement. These nuclei contain inhibitory GABAergic and excitatory glutamatergic neurons, whose molecular diversity and development remain largely unraveled. Many tegmental neurons originate in the embryonic ventral rhombomere 1 (r1), where GABAergic fate is regulated by the transcription factor (TF) Tal1. We used single-cell mRNA sequencing of the mouse ventral r1 to characterize the Tal1-dependent and independent neuronal precursors. We describe gene expression dynamics during bifurcation of the GABAergic and glutamatergic lineages and show how active Notch signaling promotes GABAergic fate selection in post-mitotic precursors. We identify GABAergic precursor subtypes that give rise to distinct tegmental nuclei and demonstrate that Sox14 and Zfpm2, two TFs downstream of Tal1, are necessary for the differentiation of specific tegmental GABAergic neurons. Our results provide a framework for understanding the development of cellular diversity in the tegmental nuclei.

A Platelet Function Modulator of Thrombin Activation Is Causally Linked to Cardiovascular Disease and Affects PAR4 Receptor Signaling

Dual antiplatelet therapy reduces ischemic events in cardiovascular disease, but it increases bleeding risk. Thrombin receptors PAR1 and PAR4 are drug targets, but the role of thrombin in platelet aggregation remains largely unexplored in large populations. We performed a genome-wide association study (GWAS) of platelet aggregation in response to full-length thrombin, followed by clinical association analyses, Mendelian randomization, and functional characterization including iPSC-derived megakaryocyte and platelet experiments. We identified a single sentinel variant in the GRK5 locus (rs10886430-G, p = 3.0 × 10^-42) associated with increased thrombin-induced platelet aggregation (β = 0.70, SE = 0.05). We show that disruption of platelet GRK5 expression by rs10886430-G is associated with enhanced platelet reactivity. The proposed mechanism of a GATA1-driven megakaryocyte enhancer is confirmed in allele-specific experiments. Utilizing further data, we demonstrate that the allelic effect is highly platelet- and thrombin-specific and not likely due to effects on thrombin levels. The variant is associated with increased risk of cardiovascular disease outcomes in UK BioBank, most strongly with pulmonary embolism. The variant associates with increased risk of stroke in the MEGASTROKE, UK BioBank, and FinnGen studies. Mendelian randomization analyses in independent samples support a causal role for rs10886430-G in increasing risk for stroke, pulmonary embolism, and venous thromboembolism through its effect on thrombin-induced platelet reactivity. We demonstrate that G protein-coupled receptor kinase 5 (GRK5) promotes platelet activation specifically via PAR4 receptor signaling. GRK5 inhibitors in development for the treatment of heart failure and cancer could have platelet off-target deleterious effects. Common variants in GRK5 may modify clinical outcomes with PAR4 inhibitors, and upregulation of GRK5 activity or signaling in platelets may have therapeutic benefits.

Novel Genetic Variations in Acute Myeloid Leukemia in Pakistani Population

Acute myeloid leukemia (AML) is a hematological malignancy characterized by clonal expansion of blast cells that exhibit great genetic heterogeneity. In this study, we describe the mutational landscape and its clinico-pathological significance in 26 myeloid neoplasm patients from a South Asian population (Pakistan) by using ultra-deep targeted next-generation DNA sequencing of 54 genes (∼5000×) and its subsequent bioinformatics analysis. The data analysis indicated novel non-silent somatic mutational events previously not reported in AML, including nine non-synonymous and one stop-gain mutations. Notably, two recurrent somatic non-synonymous mutations, i.e., STAG2 (causing p.L526F) and BCORL1 (p.A400V), were observed in three unrelated cases each. The BCOR was found to have three independent non-synonymous somatic mutations in three cases. Further, the SRSF2 with a protein truncating somatic mutation (p.Q88X) was observed for the first time in AML in this study. The prioritization of germline mutations with ClinVar, SIFT, Polyphen2, and Combined Annotation Dependent Depletion (CADD) highlighted 18 predicted deleterious/pathogenic mutations, including two recurrent deleterious mutations, i.e., a novel heterozygous non-synonymous SNV in GATA2 (p.T358P) and a frameshift insertion in NPM1 (p.L258fs), found in two unrelated cases each. The WT1 was observed with three independent potential detrimental germline mutations in three different cases. Collectively, non-silent somatic and/or germline mutations were observed in 23 (88.46%) of the cases (0.92 mutation per case). Furthermore, the pharmGKB database exploration showed a missense SNV rs1042522 in TP53, exhibiting decreased response to anti-cancer drugs, in 19 (73%) of the cases. This genomic profiling of AML provides deep insight into the disease pathophysiology. Identification of pharmacogenomics markers will help to adopt personalized approach for the management of AML patients in Pakistan.