O-GlcNAc homeostasis contributes to cell fate decisions during hematopoiesis

The addition of a single β-d-GlcNAc sugar (O-GlcNAc) by O-GlcNAc-transferase (OGT) and O-GlcNAc removal by O-GlcNAcase (OGA) maintain homeostatic O-GlcNAc levels on cellular proteins. Changes in protein O-GlcNAcylation regulate cellular differentiation and cell fate decisions, but how these changes affect erythropoiesis, an essential process in blood cell formation, remains unclear. Here, we investigated the role of O-GlcNAcylation in erythropoiesis by using G1E-ER4 cells, which carry the erythroid-specific transcription factor GATA-binding protein 1 (GATA-1) fused to the estrogen receptor (GATA-1–ER) and therefore undergo erythropoiesis after β-estradiol (E2) addition. We observed that during G1E-ER4 differentiation, overall O-GlcNAc levels decrease, and physical interactions of GATA-1 with both OGT and OGA increase. RNA-Seq–based transcriptome analysis of G1E-ER4 cells differentiated in the presence of the OGA inhibitor Thiamet-G (TMG) revealed changes in expression of 433 GATA-1 target genes. ChIP results indicated that the TMG treatment decreases the occupancy of GATA-1, OGT, and OGA at the GATA-binding site of the lysosomal protein transmembrane 5 (Laptm5) gene promoter. TMG also reduced the expression of genes involved in differentiation of NB4 and HL60 human myeloid leukemia cells, suggesting that O-GlcNAcylation is involved in the regulation of hematopoietic differentiation. Sustained treatment of G1E-ER4 cells with TMG before differentiation reduced hemoglobin-positive cells and increased stem/progenitor cell surface markers. Our results show that alterations in O-GlcNAcylation disrupt transcriptional programs controlling erythropoietic lineage commitment, suggesting a role for O-GlcNAcylation in regulating hematopoietic cell fate.

Hematopoiesis begins early in embryonic development and continues throughout adulthood to produce and replenish the various cell types comprising the blood system. Distinct cellular lineages are derived from hematopoietic stem cells via the activation of cell-specific gene expression programs, spatiotemporally controlled by lineage-specific transcription factors that activate or inactivate downstream genes (1) The molecular mechanisms by which lineage-specific transcriptional programs are established and maintained in different hematopoietic cell types are broadly understood, but the role of post-translational modifications (PTMs) 5 in these transcriptional cascades is not well defined. PTMs have many fundamental biological functions and present an exciting new avenue to explore in the context of hematopoiesis. Understanding how extracellular signals and the cellular environment regulate hematopoiesis through PTMs has significant implications regarding the modulation of lineage-specific gene expression and offers the potential for identification and development of novel therapies to treat blood-related disorders.
Erythropoiesis is the process by which hematopoietic stem cells and progenitors proliferate and differentiate into mature erythrocytes. Erythroid differentiation is regulated by lineagespecific transcription factors that orchestrate proper cellular differentiation. GATA-1, a dual zinc finger DNA-binding protein that recognizes WGATAR motifs, is a key erythroid-specific transcription factor that coordinates erythropoietic gene programs during maturation (2)(3)(4)(5). At the onset of erythroid differentiation and as it proceeds, many GATA-1 target genes are activated or repressed (6 -9). The co-occupancy of DNA-binding sites by GATA-1 and other erythroid-specific or ubiquitous transcription factors, including coactivators and corepressors, determines the gene transcription status to some extent (7, 10 -12). Epigenetic modification of histones and variation of WGATAR-adjacent DNA-binding motifs play a role (7,8,(13)(14)(15), but these alone do not adequately explain how GATA-1-containing protein complexes are modulated to positively or negatively regulate gene expression.
Post-translational modification of transcription complexes during hematopoiesis might regulate lineage-specific activity by altering protein-protein interactions within these complexes or by modifying nearby chromatin structure (16 -19). Recently, we demonstrated that O-GlcNAcylation modulates GATA-1-mediated repression of the A ␥-globin promoter (20). O-GlcNAcylation is the modification of serine and threonine residues by ␤-D-GlcNAc (O-GlcNAc). This monomeric sugar moiety is added and removed by the O-GlcNAc-processing enzymes, O-GlcNAc-transferase (OGT) and O-GlcNAcase (OGA), respectively (21). Dynamic O-GlcNAc cycling is critical for the proper regulation of gene transcription. O-GlcNAc is one of the histone code marks (22,23), and it modulates other epigenetic changes on chromatin, such as DNA methylation (24,25), histone acetylation (26), and histone methylation (25,27). TET (ten-eleven translocation) proteins and host cell factor 1 colocalize with OGT at gene promoters associated with activating histone marks to promote target gene transcription (23)(24)(25)28). OGT also interacts with mSin3A or components of the polycomb repressive complexes 1 and 2 to mediate gene silencing by forming stable repressor complexes (26,27,29). In addition, OGT and OGA regulate preinitiation complex formation and RNA polymerase II activity (30 -32).
OGT and OGA interact with the GATA-1-FOG-1-Mi2␤ (NuRD) repressor complex at the A ␥-globin promoter, and O-GlcNAcylation of the NuRD chromatin-remodeling complex is necessary for GATA-1-FOG-1-Mi2␤ repressor formation (20). Other studies have demonstrated that O-GlcNAc cycling on chromatin modifiers is critical for modulation of gene transcription and normal development (33,34). Thus, O-GlcNAc might regulate hematopoietic cell differentiation and maturation. During embryonic stem cell differentiation into neurons, cellular O-GlcNAc levels decrease (35) as well as those in epidermal keratinocytes (36), cardiomyocytes (37), and myoblasts (38). O-GlcNAc also decreases in K562 erythroleukemia cells following treatments that induce terminal differentiation (20). However, the overall O-GlcNAc level increases during differentiation of mouse chondrocyte (39) and osteoblast cells (40). These data suggest the hypothesis that O-GlcNAc cycling contributes to the regulation of gene expression programs during hematopoietic cell differentiation.
To test this hypothesis, we used well-established erythroid and myeloid cell systems that can be induced to undergo terminal differentiation. For an erythroid model, we used the G1E-ER4 cell line, an erythroid progenitor cell line derived from murine GATA-1 Ϫ erythroid (G1E) cells that stably express GATA-1-ER, a fusion product of GATA-1 and human estrogen receptor ligand-binding domain (6). GATA-1 is the master regulatory transcription factor for erythropoiesis. When GATA-1 activity is restored following the addition of ␤-estradiol (E 2 ), G1E-ER4 cells undergo erythroid differentiation (6,41). For a myeloid cell model, we used the NB4 and HL60 myeloid leukemia cell lines, which differentiate into neutrophil-like cells following addition of all-trans-retinoic acid (42). In this work, we demonstrate that 1) erythropoietic cell differentiation resulted in a marked reduction of total O-GlcNAc levels; 2) interaction among GATA-1-ER, OGT, and OGA was enhanced after E 2 -induced differentiation; 3) perturbation of O-GlcNAc cycling affected 1,173 differentially expressed genes, including 47 erythroid-specific GATA-1 target genes; 4) OGA inhibition by Thiamet-G (TMG) led to an amplification of GATA-1-mediated activation or repression of GATA-1 target genes; 5) TMG treatment impaired myeloid cell differentiation into neutrophil-like cells; and 6) prolonged TMG treatment disrupted erythroid differentiation in G1E-ER4 cells. We conclude that O-GlcNAcylation modulates lineage-specific transcriptional cascades, thus contributing to the control of hematopoietic differentiation.

O-GlcNAc levels decrease during erythropoiesis
To determine whether O-GlcNAcylation is involved in erythroid differentiation, we treated GATA-1-ER cells with E 2 to trigger erythroid differentiation. Expression of GATA-1-ER increased as reported previously (43,44); Gata2 transcription was repressed by GATA-1, resulting in the decrease of GATA-2 protein level (15). Expression of friend of GATA-1 (FOG-1) also increased as reported previously (6) (Fig. 1a). Interestingly, the total cellular O-GlcNAc level was dramatically decreased after E 2 treatment. We found that the level continued to decrease during differentiation, reaching its lowest level at 30 h post-E 2 administration. The decreased O-GlcNAc level was associated with a reduction of OGT protein levels and an increase of OGA protein levels ( Fig. 1, a, c, and e). However, the large decrease in cellular O-GlcNAcylation after E 2 treatment suggests that the targeting of OGT and OGA to specific substrates might also be altered during differentiation. The O-GlcNAc level, as well as OGT and OGA expression, did not change over time in G1E negative control cells ( Fig. 1, a, b, and d), which, unlike G1E-ER4 cells, lack inducible GATA-1 function (Fig.  1a). These data indicate that the O-GlcNAc level decreases during erythropoiesis and that the decline was a direct result of GATA-1 activity.
We treated GATA-1-ER cells for 30 h with both E 2 and TMG, a highly selective inhibitor of OGA that blocks removal of the O-GlcNAc moiety by the enzyme, to determine whether the reduction in total cellular O-GlcNAcylation is required for proper erythroid differentiation (45)(46)(47). Cells treated with TMG maintain the O-GlcNAc moiety on target proteins, leading to increased global levels of O-GlcNAcylation and altered O-GlcNAc cycling. TMG was added at the onset of E 2 -initiated differentiation (acute TMG treatment). Western blotting analyses showed an increased O-GlcNAc level following TMG treatment. FOG-1, GATA-1-ER, and GATA-2 protein levels were unaffected by TMG treatment (Fig. 1f). Morphologic mat-O-GlcNAc homeostasis regulates hematopoiesis uration (Fig. 1g) and hemoglobin induction (Fig. 1h) were similar in E 2 -induced, Ϯacute TMG G1E-ER4 cells, and TMG did not affect ␤ major -globin gene transcription (Fig. 1k). We did not observe changes in morphology (Fig. 1i) or hemoglobin induction in control G1E cells treated with E 2 Ϯ acute TMG (Fig. 1j), indicating that these cells did not undergo erythroid differentiation. We conclude that there were no significant observed differences in cell morphology, hemoglobin induction, or ␤ major -globin gene transcription between E 2 -induced, Ϯacute TMG-treated G1E-ER4 cells within the first 30 h of induction.

O-GlcNAc homeostasis regulates hematopoiesis
However, short-term TMG treatment did affect the overall O-GlcNAc levels in G1E-ER4 cells, suggesting the potential to disrupt cellular pathways involved in erythroid differentiation.

Prolonged OGA inhibition impairs erythropoiesis
Previously, we demonstrated in cell lines and mouse tissue that sustained TMG treatment (14 days or longer) alters the metabolic demands of a cell and reprograms the transcriptome, leading to a variety of adaptive changes to TMG (48); therefore, we hypothesized that long-term TMG treatment would dramatically alter erythropoiesis compared with acute TMG treatment. G1E-ER4 cells were treated for 2 weeks with a daily dose of TMG and then differentiated with E 2 . To analyze how perturbation of O-GlcNAc cycling affects erythropoiesis, we performed flow cytometry on acute (TMG ϩ E 2 ) and long-term (2-week) TMG-treated G1E-ER4 cells. There was no significant difference after 30 h of E 2 induction between treatment groups (data not shown). However, after 72 h of E 2 induction, acute TMG treatment increased the expression of all cell surface markers (Fig. 2a). This increase was enhanced further in longterm TMG-treated G1E-ER4 cells (Fig. 2b). Overall O-GlcNAc levels and OGA expression increased during long-term TMG exposure, whereas OGT expression decreased (Fig. 2c). Longterm TMG exposure led to a slight decrease of GATA-1-ER pre-E 2 induction and an increase in FOG-1 expression pre-and post-E 2 induction (Fig. 2c). Sustained TMG treatment resulted in a 10% reduction of hemoglobin-positive cells compared with E 2 -only treatment (Fig. 2, d and e). Sustained TMG ϩ E 2 treatment reduced ␤ major -globin expression by 34.5% compared with E 2 treatment only (Fig. 2f). Laptm5 and Fndc5 gene expression increased 2.8-and 2.6-fold, respectively, with sustained TMG treatment and E 2 (Fig. 2, g and h). These data demonstrate that sustained TMG treatment slows G1E-ER4 differentiation and alters the expression of erythroid genes.

Myeloid to neutrophil-like differentiation is impaired with TMG treatment
Based on our G1E-ER4 cell data, we hypothesized that O-GlcNAcylation might be involved in the differentiation of other hematopoietic lineages. Human NB4 and HL60 myeloid leukemia cells differentiate into neutrophil-like cells following the addition of all-trans-retinoic acid (ATRA) (42). Therefore, we treated cells with ATRA for 48 h prior to harvesting and analysis. In both cell lines, we measured a decline in O-GlcNAc levels. The O-GlcNAcylation change in the NB4 cells affected a large number of O-GlcNAc-modified proteins, whereas in HL60 cells proteins at 76, 44, 37, and 12 kDa showed declines in O-GlcNAcylation (Fig. 3a). Similar to G1E-ER4 cells, OGT levels did not decline during differentiation, suggesting that OGT substrate-targeting proteins change after ATRA induction. Next, we treated NB4 cells with TMG concomitantly with ATRA induction and measured transcription of the cathepsin D (CTSD) and defensin ␣1 (DEFA1) genes, which increase during neutrophil-like differentiation (49,50). ATRA treatment led to a robust increase in CTSD (Fig. 3b) and DEFA1 (Fig. 3c) mRNA expression. TMG treatment decreased CTSD and DEFA1 gene expression following ATRA induction. These data phenocopy gene expression changes during G1E-ER4 cell dif-ferentiation and suggest that 1) O-GlcNAc levels decline as hematopoietic cells differentiate, 2) substrate targeting of OGT and OGA changes during differentiation, and 3) disruption in O-GlcNAc homeostasis impairs differentiation.

GATA-1 interacts with OGT and OGA in G1E-ER4 cells
Our previous study in human erythroleukemia K562 cells or murine chemical inducer of dimerization-dependent human ␤-globin locus yeast artificial chromosome bone marrow cells demonstrated that OGT and OGA interact with the NuRD repressor complex at the Ϫ566 GATA site of the A ␥-globin gene promoter and that the Mi2␤-OGT-OGA interactions change when ␥-globin transcription is altered (20). Thus, we assessed whether OGT and OGA interact with GATA-1-ER in G1E-ER4 cell lysates treated with TMG and/or E 2 via coimmunoprecipitation experiments using antibodies against O-GlcNAc (RL2), OGT, OGA, and GATA-1. We found that GATA-1-ER was poorly coimmunoprecipitated by O-GlcNAc, OGT, and OGA antibodies prior to E 2 treatment. Following E 2 treatment, O-GlcNAc, OGT, and OGA robustly coimmunoprecipitated GATA-1-ER ( Fig. 4, a-c). O-GlcNAc, OGT, and OGA coimmunoprecipitated FOG-1 before and after E 2 treatment, although E 2 treatment increased the amount of coimmunoprecipitated FOG-1 with a molecular weight shift suggesting an increase in FOG-1 phosphorylation (51, 52) ( Fig. 4, a-c). GATA-1 antibody coimmunoprecipitated OGT and OGA before and after E 2 treatment (Fig. 4d). These data demonstrate that the interaction among GATA-1-ER, OGT, and OGA, as well as the O-GlcNAcylation of GATA-1-ER-interacting proteins, increases during G1E-ER4 cell differentiation.

Inhibition of OGA with TMG changes erythroid gene transcription networks
Based on our previous work (20) and the current study, we have observed that OGT and OGA interact with GATA-1 in multiple erythroid cell lines; next, we tested whether disruptions of O-GlcNAc homeostasis altered the erythroid gene transcription network during G1E-ER4 differentiation. We performed next generation RNA-sequence (RNA-Seq) analysis on G1E-ER4 cells without treatment, with acute TMG or E 2 treatment only, and with E 2 ϩ acute TMG treatment for 30 h because this time point corresponds to previously published GATA-1 next generation chromatin immunoprecipitation (ChIP)-sequencing data (7). 8,271 transcripts (Table S1) were included in the RNA-Seq analysis from the different treatments, and the shared genes that were down-regulated (Fig. 5a) or up-regulated (Fig. 5b) are shown. Our RNA-Seq data suggest that inhibition of OGA by TMG changes transcription levels of numerous genes during erythropoiesis. Gene set enrichment analysis (GSEA) of the 8,271 transcripts was performed to characterize the biological processes altered by inhibition of OGA (53,54). The top 10 up-regulated (top panel) and down-regulated (bottom panel) biological processes are listed comparing acute TMG ϩ E 2 with E 2 treatment only (Fig. 5c). Examples of enrichment analysis plots are shown in Fig. 5d. Interestingly, biological processes of myeloid-leukocyte activation and inflammatory response were up-regulated. To further evaluate O-GlcNAc homeostasis regulates hematopoiesis Figure 2. Sustained OGA inhibition modulates erythropoietic gene expression. TMG was added to G1E-ER4 cells at the same time as E 2 induction. Treatment was either acute (72 h) or long-term (2 weeks). Flow cytometry analysis was performed on E 2 -only (data not shown), E 2 ϩ acute TMG (a), and long-term TMG treatment (2 weeks (2w)) ϩ E 2 (b). Whole cell lysates were prepared at 0 and 72 h for E 2 -only and long-term TMG treatment ϩ E 2 and subjected to immunoblotting (c). Cells were analyzed histologically by May-Grünwald Giemsa (MGG) staining (d) for cell morphology and benzidine staining (BZ) (e) for hemoglobin content. Transcription levels of the ␤ major -globin gene (Hbb-b1) (f), Laptm5 (g), and Fndc5 (h) were measured by RT-qPCR; ribosome protein L-7 was used as an internal control. Boxes represent cropped blots. All experiments were performed with at least three biological replicates. Error bars represent S.D. WB, Western blotting.

O-GlcNAc homeostasis regulates hematopoiesis
our RNA-Seq data, we compared the transcription profile of G1E-ER4 cells treated with E 2 ϩ acute TMG with that of E 2 only-treated cells (Fig. 6a). 1,173 transcripts (Table S2) were differentially expressed in E 2 ϩ acute TMG compared with E 2 treatment only; 530 genes were up-regulated, and 643 were down-regulated (Fig. 6, b and c), including 47 erythroid-specific genes (Table S3) (55). Together, these results indicate that transcriptional networks were changed when OGA was inhib-ited with TMG during erythropoiesis and that alteration of O-GlcNAc cycling in GATA-1-ER cells modulates terminal differentiation of erythroid gene programs.

O-GlcNAc cycling modulates GATA-1-ER target gene transcription during erythropoiesis
To evaluate how O-GlcNAc cycling modulates GATA-1-ER target gene transcription, we used an online database, Harmo-

O-GlcNAc homeostasis regulates hematopoiesis
nizome (56), to search for GATA-1 target genes present in G1E-ER4 cells after E 2 treatment. A total of 4,072 GATA-1 target genes were selected from the database by combining nine sets of mouse GATA-1-binding site profiles in G1E-ER4 cells from the ENCODE database (56, 57). 433 of these genes (Table S4) were differentially expressed in E 2 ϩ acute TMG-treated cells compared with E 2 only-treated cells (Fig. 6d). Of the 433 genes, 243 were up-regulated, and 190 were down-regulated (Fig. 6d). GSEA was performed on these 433 genes. Immune response and defense response biological processes, which are characteristic of leukocytes, were observed to be up-regulated in E 2 ϩ acute TMG-treated cells compared with E 2 only-treated cells (data not shown).
Next, four genes were selected for reverse transcriptionquantitative polymerase chain reaction (RT-qPCR) analysis to assess the quality of our RNA-Seq data. Laptm5, Fndc5, Parp14, and Mvb12a are all GATA-1 target genes that showed robust transcription level changes following E 2 ϩ acute TMG treatment compared with E 2 -only treatment in our RNA-Seq data analysis. Following E 2 -only treatment for 30 h, Laptm5, Fndc5, and Parp14 transcription levels were increased compared with control ( Fig. 7, a-c), whereas the Mvb12a transcription level decreased (Fig. 7d). When cells were treated with E 2 ϩ acute TMG, Laptm5, Fndc5, and Parp14 transcription levels were increased further compared with E 2 -only treatment (Fig. 7,  a-c), whereas Mvb12a transcription decreased further (Fig.

O-GlcNAc homeostasis regulates hematopoiesis
7d). These results were consistent with the -fold changes from our RNA-Seq data. We repeated identical measurements on the expression of these genes in G1E cells (negative control); no change was detected in the expression of these genes with TMG, E 2 , or E 2 ϩ acute TMG treatment (data not shown). These data confirm that subsets of GATA-1 target genes are affected by alterations in O-GlcNAc homeostasis during erythropoiesis.

Inhibition of OGA decreases the occupancy of the Laptm5 intronic GATA-binding site by GATA-1, O-GlcNAc, OGT, and OGA
Based on our RNA-Seq and GSEA analysis, we observed that inhibition of OGA by TMG altered a subset of GATA-1 target genes during erythropoiesis, especially those associated with immune response and defense response genetic pathways. Laptm5, a target of GATA-1-ER activation, is specifically

O-GlcNAc homeostasis regulates hematopoiesis
expressed in hematopoietic cells and is preferentially expressed in immune cells (58,59). LAPTM5 positively modulates inflammatory signaling pathways and cytokine secretion in macrophages (59) and is up-regulated in G1E-ER4 after E 2 treatment for 30 h. We performed ChIP assays at the GATA-binding site located in the first intron of Laptm5, ϳ8 kb downstream from the transcription start site (TSS) (7). Using antibodies against GATA-1, OGT, OGA, and O-GlcNAc to assess the occupancy of this site by GATA-1-ER, OGT, OGA, and O-GlcNAc, we observed in untreated control G1E-ER4 cells an enrichment in the occupancy of OGT, OGA, and an O-GlcNAcylation (Fig. 8,  b-d), suggesting that O-GlcNAcylation at the Laptm5 GATAbinding site might be essential for maintaining the basal level of Laptm5 transcription. Following GATA-1-ER restoration by E 2 , GATA-1-ER occupancy dramatically increased; however, E 2 ϩ acute TMG treatment reduced the amount of GATA-1-ER occupancy at this site (Fig. 8a, left panel). OGT occupancy did not change, and OGA occupancy increased after E 2 treatment compared with the untreated control (Fig. 8, b and c,  left panel). We hypothesize that this contributes to the decrease in overall O-GlcNAc levels at this GATA-binding site (Fig. 8d,  left panel).
Interestingly, the occupancy of OGT and OGA and the overall O-GlcNAc level at this GATA-binding site decreased after E 2 ϩ acute TMG treatment compared with E 2 -only treatment ( Fig. 8, b-d, left panel). A similar pattern was observed when acute TMG-only treatment was compared with the untreated control ( Fig. 8, b-d, left panel). To show that these changes were specific to the GATA-binding site, we used a non-GATA-1-binding DNA region 50 kb downstream of the GATA-binding site as a negative control. We did not observe any changes when comparing E 2 ϩ acute TMG treatment with E 2 -only treatment (Fig. 8, a-d, right panel). These data suggest that O-GlcNAc homeostasis plays a role in GATA-1 regulation of the Laptm5 gene and that acute TMG treatment alters the occupancy by GATA-1-ER, OGT, and OGA at this binding site.

Discussion
In this study, we utilized a well-established cell model of erythropoiesis, G1E-ER4 cells, to investigate the role of O-GlcNAc during GATA-1-mediated erythroid differentiation. After the restoration of GATA-1 activity, overall O-GlcNAc levels dramatically decreased with an increase of OGA and a decrease of OGT protein levels. The decrease in O-GlcNAcylation during differentiation was also seen in two myeloid leukemia cell lines after induction into more neutrophil-like cells. In addition, GATA-1-ER restoration promoted the interaction of GATA-1-ER with OGT and OGA. RNA-Seq analysis of G1E-ER4 cells treated with E 2 ϩ acute TMG dem-  O-GlcNAc homeostasis regulates hematopoiesis onstrated a change in transcription level of 1,173 genes compared with cells treated with E 2 only, including 433 GATA-1 target genes. The transcription level changes were orthogonally validated on selected genes by RT-qPCR. One of these genes, Laptm5, was subjected to further proof-of-principle analysis. ChIP data demonstrated that the occupancy of the Laptm5 GATA-binding site by GATA-1 and OGA was increased and that the overall O-GlcNAc level was decreased following GATA-1 activation; inhibition of OGA with TMG diminished these changes. Acute TMG treatment impaired expression of neutrophil genes during ATRA-mediated differentiation of myeloid cells. Interestingly, sustained TMG treatment in G1E-ER4 reduced ␤ major -globin gene expression but increased expression of Fndc5 and Laptm5. Together, these data suggest that O-GlcNAc plays a role in regulating a subset of GATA-1-controlled erythroid genes and that disruption of O-GlcNAc homeostasis impairs erythropoietic cell differentiation, potentially maintaining the cells at a more progenitor-like stage.

O-GlcNAc homeostasis regulates hematopoiesis
We previously demonstrated that GATA-1 interacts with OGT and OGA in mouse erythroleukemia BirA cells (20), and our current results validate these findings. O-GlcNAcylation of GATA-1 was not detected, but GATA-1interacting proteins are likely modified by O-GlcNAc because following differentiation, the O-GlcNAc antibody pulled down GATA-1. The interaction of OGT and OGA with GATA-1-ER and O-GlcNAcylation of GATA-1-ERinteracting proteins during erythropoiesis could change protein-protein interactions within GATA-1-containing transcriptional complexes, thereby contributing to GATA-1-mediated gene activation or repression. Our RNA-Seq analysis suggests that inhibition of OGA by TMG alters gene transcription during erythropoiesis and that these changes may reflect either reprogramming of the cell to a more myeloid-leukocyte-like lineage or delayed progenitor stage differentiation. Thus, O-GlcNAc may play a role in determining hematopoietic cell fate at an earlier developmental stage. Inhibition of OGA by TMG in hematopoietic stem cells prior to lineage commitment may change the O-GlcNAcylation pattern and skew hematopoietic stem cell differentiation toward a particular lineage. Previous reports show that elevated O-GlcNAc in mouse embryonic stem cells delays the onset of neuronal differentiation (35), whereas increased O-GlcNAcylation skews human pluripotent stem cell differentiation toward adipose mesoderm markers and away from ectoderm markers (29). Because myeloid leukemia cells also show reduced O-GlcNAcylation during differentiation, OGA inhibition might reduce hematopoietic differentiation in general or shift commitment toward another hematopoietic lineage.
Of the 1,173 differentially expressed genes (E 2 ϩ acute TMG compared with E 2 only), we determined that 47 of these genes are erythroid-specific genes, including Gfi1b (61), Alas2 (62), and Gata2 (3), which are critical for normal erythroid differentiation. The Gfi1b gene encodes for a transcription factor important in promoting erythrocyte differentiation, Gata2 expression inhibits erythrocyte differentiation, and the Alas2 gene is critical in porphyrin biosynthesis. We also found that 433 of the differentially expressed genes were under control of GATA-1 with 243 genes up-regulated and 190 down-regulated. Notably, 84% of activated GATA-1 target genes became more activated in the presence of TMG, whereas 79% of repressed GATA-1 target genes became more repressed, indicating that inhibition of OGA amplifies GATA-1-mediated transcriptional effects. One of the GATA-1 target genes is Laptm5, which is expressed in hematopoietic cells (58). The Laptm5 gene encodes a transmembrane receptor that is associated with lysosomes and has been shown to positively modulate inflammatory signaling pathways and cytokine secretion in macrophages (59). Transcription of Laptm5 was activated following GATA-1 restoration in G1E-ER4 cells, and the activation was enhanced in the presence of TMG. Strikingly, the occupancy of the Laptm5 GATA-binding site by GATA-1-ER, OGT, and OGA after inhibition of OGA decreased as did the O-GlcNAc level. Even though GATA-1-ER was decreased following E 2 ϩ acute TMG treatment, TMG amplified the increase of Laptm5 transcription compared with E 2 treatment only. Thus, TMG treatment and subsequent retention of O-GlcNAc-modified proteins may stabilize the interaction of GATA-1 with other transcription factors, thereby reducing the recruitment or turnover of these proteins. Hence, even with reduced occupancy by GATA-1 at the GATA-binding site, the remaining GATA-1 complex may be more stable or more active in promoting transcription. These data further argue for the importance of the rate of O-GlcNAc cycling at GATA-binding sites. Reduction of O-GlcNAc cycling at the Laptm5 gene with E 2 ϩ acute TMG treatment may lead to the recruitment of other transcriptional coactivators or stabilize the GATA-1-containing activator complex at the Laptm5 GATA-binding site, resulting in increased transcription.
Based on the Laptm5 data, we propose the following general model for OGT, OGA, and GATA-1 regulation of transcription during erythropoiesis. Following GATA-1 restoration by E 2 , GATA-1 transcriptional activator complexes are recruited to GATA-binding sites located in noncoding regions of genes, such as promoters and introns (7, 57), resulting in increased gene expression (Fig. 9a). OGT interacts with and modifies a component or components of these complexes, thereby promoting the recruitment of other transcriptional coactivators. OGA then removes O-GlcNAc from one or more of the proteins in the complexes, leading to complete recruitment of the promoter (Fig. 9b). When OGA is inhibited by TMG, O-GlcNAc cycling at the activator complexes is disrupted, leading to decreased cycling of OGT, OGA, and GATA-1 at the GATA-binding sites; however, a stable population of activator or repressor complex stays bound to the promoter, amplifying the regulation on the gene (Fig. 9c). In addition, it is possible that the components of the activator complex might also change in the absence of OGA function. Regardless, we hypothesize that this alternate activator/repressor complex is likely to be more stable and active due to the changes in protein-protein interactions within the complexes or alterations of the nearby chromatin structure, resulting in enhanced transcription or repression of the target genes.

O-GlcNAc homeostasis regulates hematopoiesis
We hypothesize that O-GlcNAcylation is required for proper hematopoietic cell differentiation. Our data are supported by other published studies of differentiation (20,(35)(36)(37)(38) that also show a decrease of global O-GlcNAcylation level, suggesting that in most cells differentiation requires this reduction. The relationship of reduced O-GlcNAc levels to differentiation deserves further study. Many cells enter the cell cycle before differentiation, suggesting that a decrease in O-GlcNAc might enhance cell cycle progression (63). Alternately, cells might redistribute metabolites away from the hexosamine biosyn-thetic pathway that makes UDP-GlcNAc, the metabolic substrate for OGT, to other metabolic pathways needed for growth (64).

Cell culture
G1E and G1E-ER4 cells were kind gifts from Dr. Soumen Paul at the University of Kansas Medical Center. These lines were cultured as described previously (6,41). Cells were incubated at 37°C in 5% CO 2 in a 95% humidified incubator. Cells were treated with 10 M TMG (S.D. Specialty Chemicals) and/or 1 M E 2 (E2758, Sigma) (41). After 30 h, cells were harvested for subsequent analysis. For sustained TMG treatment, TMG was added to cells every day (10 M) for 2 weeks prior to E 2 induction. NB4 and HL60 cells were grown as described previously (66). ATRA (R2625, Sigma) was used at 1 M in DMSO.

Western blotting and IP
Cells were lysed on ice as described previously (20). Wholecell lysates were used for electrophoresis and subjected to Western blotting as described previously (20). All Western blotting results were repeated in three independent experiments, and the representative images are shown. IP was performed as described previously (20). We used 2 mg of cell lysates and 2 g of antibody for each reaction. OGA and OGT relative protein levels were measured by analyzing the band densities using ImageJ 1.49 (67) and then normalized to the density of GAPDH.

Benzidine staining and May-Grünwald/Giemsa stain
Benzidine staining and May-Grünwald/Giemsa stain were performed as described previously (6) to assess the hemoglobin content and morphology changes, respectively. Cell stain images were visualized and photographed using a Nikon Eclipse 80i digital microscope (Nikon Instruments Inc., Melville, NY).

Total RNA isolation and RT-PCR
Total RNA was isolated from 5 ϫ 10 6 cells using TRI Reagent (T9424, Sigma) according to the manufacturer's instructions. For the RT reactions, 0.5 g of total RNA was reverse transcribed to cDNA using iScript Reverse Transcription Supermix (170-8841, Bio-Rad) following the manufacturer's instructions. Reactions were incubated in a thermal cycler (Model 2720, Applied Biosystems) using the following protocol: priming for 5 min at 25°C, RT for 20 min at 46°C, RT inactivation for 1 min at 95°C, and hold at 4°C.

RNA sequencing
RNA was extracted using TRI Reagent according to manufacturer's instructions, and RNA-Seq libraries were prepared using the Illumina TruSeq Stranded mRNA Library Preparation kit Set A (RS-122-2101, Illumina) according to the manufacturer's instructions. RNA libraries were sequenced using an Illumina HiSeq 2500 Sequencing System.

RNA-Seq data analysis
FastQC (0.11.2) and RSEM (1.2.22) software were used to assess the quality of the RNA-Seq results, align the reads to the mouse genome reference GRCm38/mm10, and calculate the gene expression values. R (3.2.2) and EdgeR (3.14.0) software were used to normalize the expression values using the TMM method (weighted trimmed mean of M-values) followed by differential expression analysis. To reduce the burden of multiple testing in differential gene expression analyses, a filter was initially applied to reduce the number of genes. Genes were removed if they did not present a meaningful gene expression across all samples; only genes with counts per million of Ͼ10 in at least two samples were considered in differential expression analyses. The Benjamini and Hochberg procedure was used to control the false discovery rate (FDR). The following R packages were utilized for calculations and visualizations: EdgeR and Gplot.

O-GlcNAc homeostasis regulates hematopoiesis
Thermo Fisher Scientific) were added to the mixture according to the manufacturer's instructions followed by rotation at 4°C for 2 h. Dynabeads were washed with 1 ml of cold IP buffer six times and then mixed with 10% (w/v) Chelex 100 slurry (1421253, Bio-Rad); for input samples, 10% of sheared chromatin was used and processed at the same time as the ChIP samples. Samples were boiled at 95°C for 10 min followed by treatment with RNase A and proteinase K. Samples were boiled again for 10 min to inactivate proteinase K and centrifuged at 12,000 ϫ g for 1 min at 4°C. The supernatants were transferred to fresh tubes and stored at Ϫ20°C for subsequent qPCR analysis.

RT-qPCR
RT-qPCR was performed as described previously (20). The primer sequences for measuring target gene transcription levels are listed in Table S5. The reactions were run in a CFX96 Touch Real-Time PCR Detection System (185-5195, Bio-Rad).

Flow cytometry
Cell surface markers were characterized using flow cytometry analysis. Cells were first counted with a hemocytometer to ensure that 5 ϫ 10 5 cells were used per surface marker. Cells were centrifuged, the supernatant was discarded, and the pellets were washed twice in ice-cold PBS and 0.1% BSA (w/v). Pellets were then resuspended in 100 l of ice-cold PBS and 0.1% BSA (w/v) per surface marker test and then blocked with 1% (v/v) normal rat serum for 30 min. After washing twice in 1 ml of ice-cold PBS and 0.1% BSA (w/v), cells were again resuspended in 100 l of ice-cold PBS and 0.1% BSA (w/v) per surface marker test and then incubated with 0.5 g of antibody in the dark for 40 min. Cells were then washed twice in 1 ml of ice-cold PBS and 0.1% BSA (w/v) and resuspended in 350 l of PBS and 0.1% BSA (w/v) for flow cytometry analysis. 10,000 cells were analyzed per surface marker using an LSR II BD Biosciences flow cytometer (FACSDiva v.8 software). Dead cells were gated from the analysis.

qPCR data analysis
Quantification cycle (Cq) values were calculated by CFX Manager TM software (Bio-Rad). For cDNA RT-qPCR data, the dynamic ranges of RT and amplification efficiency were evaluated before applying the ⌬⌬Cq method to calculate relative gene expression change. The transcription level of the target gene was normalized to the internal control as -fold change. For ChIP DNA qPCR data, the Cq value was normalized as a percentage of input. Data generated in at least three independent experiments are presented as mean Ϯ S.E.; the two-tailed Student t test statistic was applied with p Ͻ 0.05 considered to be a significant difference.

GSEA
GSEA was performed according to the GSEA user guide (53,54). For gene ontology analysis, biological process gene sets from the Molecular Signature Database (MsigDB) were used. Pathways with nominal p value Ͻ5% and FDR q value Ͻ25% were included in the analysis.