Single Cell Transcriptome Analysis of GATA2 Deficiency in Hematopoiesis Modeled with Induced Pluripotent Stem Cells
Authors: Guitart Francis , Jung Moonjung , Cordes Stefan , Yu Shiqin , Zou Jizhong , Hong So Gun , Dang Vinh , Espinoza Diego A , Hsu Amy , Holland Steven M. , Townsley Danielle M. , Dunbar Cynthia E. , Winkler Thomas 
Published: October 2019

Abstract
GATA2 deficiency is a rare, inherited or sporadic genetic disorder characterized by variable onset of a pleomorphic constellation of immune, hematologic and lymphatic abnormalities linked to heterozygous mutations in the Gata2 gene. Patients develop monocyte, B cell, NK cell and dendritic cell deficiencies resulting in vulnerabilities to unusual infections. Patients with GATA2 deficiency also frequently progress to bone marrow failure, myelodysplastic syndrome and/or acute myelogenous leukemia. GATA proteins are transcription factors with central roles in early embryonic development and lineage specification. GATA2 is a master regulator of hematopoiesis, implicated in the initial generation and maintenance of hematopoietic stem cells (HSC).
Murine models recapitulate the human phenotype incompletely: GATA2 heterozygous knockout mice do not manifest loss of monocyte, B cells or NK cells; however, serial repopulation assays show decreased engraftment potential. Direct studies of primary HSC from patients with GATA2 deficiency are challenging due to the generally hypocellular marrow. We hypothesized that human pluripotent stem cells, particularly patient-specific iPSC, could be used to study potential developmental defects in GATA2 deficiency, overcoming a lack of primary HSC. In order to gain insights into the impact of human GATA2 deficiency on hematopoietic differentiation, we compared the single cell transcriptomes of HSPC differentiated from (i) iPSCs from a patient with GATA2 deficiency due to a mutation p.R337X (c.1009C>T) (ii) isogenic iPSCs created via homology-directed repair of Gata2 p.R337X, (iii) iPSCs from a healthy control and (iv) isogenic Gata2 heterozygous mutant iPSCs with monoallelic frameshift mutations in the second zinc finger domain.
Mesodermal and hematopoietic differentiation was performed under feeder-free, defined media conditions. At day 0, iPSCs were plated in mesodermal induction media containing VEGF, SCF, Activin A and Y27632 in STEMdiff APEL media. Mesodermal induction was continued until day 4, when embryoid bodies were cultured in hematopoietic specification media with SCF, FLT3L, IL3, IL6, G-CSF and BMP4 until day 16, when CD34+CD45+ iPSC-derived hematopoietic stem and progenitor cells (iHSC) were enumerated and sorted by fluorescence-activated cell sorting.
Single cell RNA-seq was performed using the 10XGenomics Chromium platform and primary analysis via CellRanger. Scater was used to filter outlier cells. Seurat as used to compute multiple manifold alignment and differential gene expression. Cell classification, pseudotemporal ordering and branch point analysis were performed with monocle. URD was used to calculate confirmatory diffusion maps and pseudotemporal ordering.
We analyzed 7,855 iHSPC (2952 from GATA2-deficient patient, 241 isogenic iHSPCs after repair of Gata2 mutation, 2,605 from a healthy volunteer and 2,057 from isogenic heterozygous Gata2 knockout iHSPCs) after filtering of outliers. We computed multiple manifold alignment to mitigate batch effects. Differential gene expression across Gata2 mutation status found that 42 out of 102 (42%) target genes of GATA2 (c.f. TRANSFAC database of curated transcription factor targets) were differentially expressed with adjusted p-values less than 0.05. Semi-supervised classification of cell-types and pseudotemporal ordering via monocle revealed two branch points, consistent with developmental branchings at the level of CLP and CMP multipotent progenitors. The numbers of cells along each branch was found to be statistically different (χ2=30.07, p-value = 3e-7) with the biggest differences noted in the lymphoid branch (state 4). Differential gene expression in this branch revealed a differential up-regulation of Notch1, CD69 and FKBPs and differential down-regulation of CD14.
In conclusion, iPSC/iHSPC differentiation models combined with single cell transcriptome analysis may be a valuable tool to identify pathways responsible for impaired hematopoietic/lymphatic development in GATA2 deficiency.
Figure. Figure.

Disclosures
Dunbar: National Institute of Health: Research Funding. Winkler:National Institute of Health: Research Funding.