Mass spectrometry-based proteomics now enables the absolute quantification of thousands of proteins in individual cell types. We used this technology to analyze the dynamic proteome changes occurring during human erythropoiesis. We quantified the absolute expression of 6,130 proteins during erythroid differentiation from late burst-forming units-erythroid (BFU-Es) to orthochromatic erythroblasts. A modest correlation between mRNA and protein expression was observed. We identified several proteins with unexpected expression patterns in erythroid cells, highlighting a breakpoint in the erythroid differentiation process at the basophilic stage. We also quantified the distribution of proteins between reticulocytes and pyrenocytes after enucleation. These analyses identified proteins that are actively sorted either with the reticulocyte or the pyrenocyte. Our study provides the absolute quantification of protein expression during a complex cellular differentiation process in humans, and it establishes a framework for future studies of disordered erythropoiesis.

Myelinated axons are organized into specialized domains critical to their function in saltatory conduction, i.e., nodes, paranodes, juxtaparanodes, and internodes. Here, we describe the distribution and role of the 4.1B protein in this organization. 4.1B is expressed by neurons, and at lower levels by Schwann cells, which also robustly express 4.1G. Immunofluorescence and immuno-EM demonstrates 4.1B is expressed subjacent to the axon membrane in all domains except the nodes. Mice deficient in 4.1B have preserved paranodes, based on marker staining and EM in contrast to the juxtaparanodes, which are substantially affected in both the PNS and CNS. The juxtaparanodal defect is evident in developing and adult nerves and is neuron-autonomous based on myelinating cocultures in which wt Schwann cells were grown with 4.1B-deficient neurons. Despite the juxtaparanodal defect, nerve conduction velocity is unaffected. Preservation of paranodal markers in 4.1B deficient mice is associated with, but not dependent on an increase of 4.1R at the axonal paranodes. Loss of 4.1B in the axon is also associated with reduced levels of the internodal proteins, Necl-1 and Necl-2, and of alpha-2 spectrin. Mutant nerves are modestly hypermyelinated and have increased numbers of Schmidt-Lanterman incisures, increased expression of 4.1G, and express a residual, truncated isoform of 4.1B. These results demonstrate that 4.1B is a key cytoskeletal scaffold for axonal adhesion molecules expressed in the juxtaparanodal and internodal domains that unexpectedly regulates myelin sheath thickness.

Recent studies have shown that transcription and alternative splicing can be mechanistically coupled. In the EPB41 (protein 4.1R) and EPB41L3 (protein 4.1B) genes, we showed previously that promoter/alternative first exon choice is coupled to downstream splicing events in exon 2. Here we demonstrate that this coupling is conserved among several vertebrate classes from fish to mammals. The EPB41 and EPB41L3 genes from fish, bird, amphibian, and mammal genomes exhibit shared features including alternative first exons and differential splice acceptors in exon 2. In all cases, the 5'-most exon (exon 1A) splices exclusively to a weaker internal acceptor site in exon 2, skipping a fragment designated as exon 2'. Conversely, alternative first exons 1B and 1C always splice to the stronger first acceptor site, retaining exon 2'. These correlations are independent of cell type or species of origin. Since exon 2' contains a translation initiation site, splice variants generate protein isoforms with distinct N-termini. We propose that these genes represent a physiologically relevant model system for mechanistic analysis of transcription-coupled alternative splicing.

Stabilization of the lipid bilayer membrane in red blood cells by its association with an underlying membrane-associated cytoskeleton has long been recognized as critical for proper red blood cell function. One of the principal connections between skeleton and bilayer is via linkages between band 3, the integral membrane protein that transports anions across the cell surface, and membrane skeletal elements including ankyrin, adducin, spectrin, and the junctional complex of the skeleton. Here, we use membrane tether formation coupled with fluorescent labeling of membrane components to examine the importance of band 3 in stabilizing the bilayer-skeletal association. In membranes from a patient deficient in band 3, the energy associated with the bilayer skeleton is approximately zero, whereas when band 3 is immobilized by ligation with the monoclonal antibody R10, the energy of association approximately doubles. Fluorescence images of tethers reveal that approximately 40% of the band 3 on the normal cell surface can be pulled into the tether, confirming a lateral segregation of membrane components during tether formation. These results validate a critical role for band 3 in stabilizing the bilayer-skeletal association in red cells.

Human erythropoiesis serves as a paradigm of physiologic cellular differentiation. This process is also of considerable interest for better understanding anemias and identifying new therapies. Here, we apply deep transcriptomic and accessible chromatin profiling to characterize a faithful ex vivo human erythroid differentiation system from hematopoietic stem and progenitor cells. We reveal stage-specific transcriptional states and chromatin accessibility during various stages of erythropoiesis, including 14,260 differentially expressed genes and 63,659 variably accessible chromatin peaks. Our analysis suggests differentiation stage-predominant roles for specific master regulators, including GATA1 and KLF1. We integrate chromatin profiles with common and rare genetic variants associated with erythroid cell traits and diseases, finding that variants regulating different erythroid phenotypes likely act at variable points during differentiation. In addition, we identify a regulator of terminal erythropoiesis, TMCC2, more broadly illustrating the value of this comprehensive analysis to improve our understanding of erythropoiesis in health and disease.

Kindlin-3 (K3) is critical for the activation of integrin adhesion receptors in hematopoietic cells. In humans and mice, K3 deficiency is associated with impaired immunity and bone development, bleeding, and aberrant erythrocyte shape. To delineate how K3 deficiency (K3KO) contributes to anemia and misshaped erythrocytes, mice deficient in erythroid (K3KO∖EpoR-cre) or myeloid cell K3 (K3KO∖Lyz2cre), knockin mice expressing mutant K3 (Q597W598 to AA) with reduced integrin-activation function (K3KI), and control wild-type (WT) K3 mice were studied. Both K3-deficient strains and K3KI mice showed anemia at baseline, reduced response to erythropoietin stimulation, and compromised recovery after phenylhydrazine (PHZ)-induced hemolytic anemia as compared with K3WT. Erythroid K3KO and K3 (Q597W598 to AA) showed arrested erythroid differentiation at proerythroblast stage, whereas macrophage K3KO showed decreased erythroblast numbers at all developmental stages of terminal erythroid differentiation because of reduced erythroblastic island (EBI) formation attributable to decreased expression and activation of erythroblast integrin α4β1 and macrophage αVβ3. Peripheral blood smears of K3KO∖EpoR-cre mice, but not of the other mouse strains, showed numerous aberrant tear drop-shaped erythrocytes. K3 deficiency in these erythrocytes led to disorganized actin cytoskeleton, reduced deformability, and increased osmotic fragility. Mechanistically, K3 directly interacted with F-actin through an actin-binding site K3-LK48. Taken together, these findings document that erythroid and macrophage K3 are critical contributors to erythropoiesis in an integrin-dependent manner, whereas F-actin binding to K3 maintains the membrane cytoskeletal integrity and erythrocyte biconcave shape. The dual function of K3 in erythrocytes and in EBIs establish an important functional role for K3 in normal erythroid function.

The fetal-to-adult hemoglobin switch is regulated in a developmental stage-specific manner and reactivation of fetal hemoglobin (HbF) has therapeutic implications for treatment of β-thalassemia and sickle cell anemia, two major global health problems. Although significant progress has been made in our understanding of the molecular mechanism of the fetal-to-adult hemoglobin switch, the mechanism of epigenetic regulation of HbF silencing remains to be fully defined. Here, we performed whole-genome bisulfite sequencing and RNA sequencing analysis of the bone marrow-derived GYPA+ erythroid cells from β-thalassemia-affected individuals with widely varying levels of HbF groups (HbF ≥ 95th percentile or HbF ≤ 5th percentile) to screen epigenetic modulators of HbF and phenotypic diversity of β-thalassemia. We identified an ETS2 repressor factor encoded by ERF, whose promoter hypermethylation and mRNA downregulation are associated with high HbF levels in β-thalassemia. We further observed that hypermethylation of the ERF promoter mediated by enrichment of DNMT3A leads to demethylation of γ-globin genes and attenuation of binding of ERF on the HBG promoter and eventually re-activation of HbF in β-thalassemia. We demonstrated that ERF depletion markedly increased HbF production in human CD34+ erythroid progenitor cells, HUDEP-2 cell lines, and transplanted NCG-Kit-V831M mice. ERF represses γ-globin expression by directly binding to two consensus motifs regulating γ-globin gene expression. Importantly, ERF depletion did not affect maturation of erythroid cells. Identification of alterations in DNA methylation of ERF as a modulator of HbF synthesis opens up therapeutic targets for β-hemoglobinopathies.

Racial differences in the pathophysiology of atherothrombosis are poorly understood. We explored the function and transcriptome of platelets in healthy black (n = 70) and white (n = 84) subjects. Platelet aggregation and calcium mobilization induced by the PAR4 thrombin receptor were significantly greater in black subjects. Numerous differentially expressed RNAs were associated with both race and PAR4 reactivity, including PCTP (encoding phosphatidylcholine transfer protein), and platelets from black subjects expressed higher levels of PC-TP protein. PC-TP inhibition or depletion blocked PAR4- but not PAR1-mediated activation of platelets and megakaryocytic cell lines. miR-376c levels were differentially expressed by race and PAR4 reactivity and were inversely correlated with PCTP mRNA levels, PC-TP protein levels and PAR4 reactivity. miR-376c regulated the expression of PC-TP in human megakaryocytes. A disproportionately high number of microRNAs that were differentially expressed by race and PAR4 reactivity, including miR-376c, are encoded in the DLK1-DIO3 locus and were expressed at lower levels in platelets from black subjects. These results suggest that PC-TP contributes to the racial difference in PAR4-mediated platelet activation, indicate a genomic contribution to platelet function that differs by race and emphasize a need to consider the effects of race when developing anti-thrombotic drugs.

Hereditary spherocytosis (HS) is the most common, nonimmune, hereditary, chronic hemolytic anemia after hemoglobinopathies. The genetic defects in membrane function causing HS lead to perturbation of the RBC metabolome, with altered glycolysis. In mice genetically lacking protein 4.2 (4.2-/-; Epb42), a murine model of HS, we showed increased expression of pyruvate kinase (PK) isoforms in whole and fractioned RBCs in conjunction with abnormalities in the glycolytic pathway and in the glutathione (GSH) system. Mitapivat, a PK activator, metabolically reprogrammed 4.2-/- mouse RBCs with amelioration of glycolysis and the GSH cycle. This resulted in improved osmotic fragility, reduced phosphatidylserine positivity, amelioration of RBC cation content, reduction of Na/K/Cl cotransport and Na/H-exchange overactivation, and decrease in erythroid vesicles release in vitro. Mitapivat treatment significantly decreased erythrophagocytosis and beneficially affected iron homeostasis. In mild-to-moderate HS, the beneficial effect of splenectomy is still controversial. Here, we showed that splenectomy improves anemia in 4.2-/- mice and that mitapivat is noninferior to splenectomy. An additional benefit of mitapivat treatment was lower expression of markers of inflammatory vasculopathy in 4.2-/- mice with or without splenectomy, indicating a multisystemic action of mitapivat. These findings support the notion that mitapivat treatment should be considered for symptomatic HS.

Defects in genes encoding ribosomal proteins cause Diamond Blackfan Anemia (DBA), a red cell aplasia often associated with physical abnormalities. Other bone marrow failure syndromes have been attributed to defects in ribosomal components but the link between erythropoiesis and the ribosome remains to be fully defined. Several lines of evidence suggest that defects in ribosome synthesis lead to "ribosomal stress" with p53 activation and either cell cycle arrest or induction of apoptosis. Pathways independent of p53 have also been proposed to play a role in DBA pathogenesis. We took an unbiased approach to identify p53-independent pathways activated by defects in ribosome synthesis by analyzing global gene expression in various cellular models of DBA. Ranking-Principal Component Analysis (Ranking-PCA) was applied to the identified datasets to determine whether there are common sets of genes whose expression is altered in these different cellular models. We observed consistent changes in the expression of genes involved in cellular amino acid metabolic process, negative regulation of cell proliferation and cell redox homeostasis. These data indicate that cells respond to defects in ribosome synthesis by changing the level of expression of a limited subset of genes involved in critical cellular processes. Moreover, our data support a role for p53-independent pathways in the pathophysiology of DBA.

Alternative pre-messenger RNA splicing remodels the human transcriptome in a spatiotemporal manner during normal development and differentiation. Here we explored the landscape of transcript diversity in the erythroid lineage by RNA-seq analysis of five highly purified populations of morphologically distinct human erythroblasts, representing the last four cell divisions before enucleation. In this unique differentiation system, we found evidence of an extensive and dynamic alternative splicing program encompassing genes with many diverse functions. Alternative splicing was particularly enriched in genes controlling cell cycle, organelle organization, chromatin function and RNA processing. Many alternative exons exhibited differentiation-associated switches in splicing efficiency, mostly in late-stage polychromatophilic and orthochromatophilic erythroblasts, in concert with extensive cellular remodeling that precedes enucleation. A subset of alternative splicing switches introduces premature translation termination codons into selected transcripts in a differentiation stage-specific manner, supporting the hypothesis that alternative splicing-coupled nonsense-mediated decay contributes to regulation of erythroid-expressed genes as a novel part of the overall differentiation program. We conclude that a highly dynamic alternative splicing program in terminally differentiating erythroblasts plays a major role in regulating gene expression to ensure synthesis of appropriate proteome at each stage as the cells remodel in preparation for production of mature red cells.

Protein 4.1R, a member of the 4.1 family, functions as a bridge between cytoskeletal and plasma membrane proteins. It is expressed in T cells, where it binds to a linker for activation of T cell (LAT) family member 1 and inhibits its phosphorylation and downstream signaling events after T cell receptor triggering. The role of the 4.1R protein in cell activation through other immunoreceptors is not known. In this study, we used 4.1R-deficient (4.1R-KO) and 4.1R wild-type (WT) mice and explored the role of the 4.1R protein in the high-affinity IgE receptor (FcεRI) signaling in mast cells. We found that bone marrow mast cells (BMMCs) derived from 4.1R-KO mice showed normal growth in vitro and expressed FcεRI and c-KIT at levels comparable to WT cells. However, 4.1R-KO cells exhibited reduced antigen-induced degranulation, calcium response, and secretion of tumor necrosis factor-α. Chemotaxis toward antigen and stem cell factor (SCF) and spreading on fibronectin were also reduced in 4.1R-KO BMMCs, whereas prostaglandin E2-mediated chemotaxis was not affected. Antibody-induced aggregation of tetraspanin CD9 inhibited chemotaxis toward antigen in WT but not 4.1R-KO BMMCs, implying a CD9-4.1R protein cross-talk. Further studies documented that in the absence of 4.1R, antigen-mediated phosphorylation of FcεRI β and γ subunits was not affected, but phosphorylation of SYK and subsequent signaling events such as phosphorylation of LAT1, phospholipase Cγ1, phosphatases (SHP1 and SHIP), MAP family kinases (p38, ERK, JNK), STAT5, CBL, and mTOR were reduced. Immunoprecipitation studies showed the presence of both LAT1 and LAT2 (LAT, family member 2) in 4.1R immunocomplexes. The positive regulatory role of 4.1R protein in FcεRI-triggered activation was supported by in vivo experiments in which 4.1R-KO mice showed the normal presence of mast cells in the ears and peritoneum, but exhibited impaired passive cutaneous anaphylaxis. The combined data indicate that the 4.1R protein functions as a positive regulator in the early activation events after FcεRI triggering in mast cells.

Ubiquitination is an enzymatic post-translational modification that affects protein fate. The ubiquitin-proteasome system (UPS) was first discovered in reticulocytes where it plays important roles in reticulocyte maturation. Recent studies have revealed that ubiquitination is a dynamic and reversible process and that deubiquitylases are capable of removing ubiquitin from their protein substrates. Given the fact that the UPS is highly active in reticulocytes, it is speculated that deubiquitylases may play important roles in erythropoiesis. Yet, the role of deubiquitylases in erythropoiesis remains largely unexplored. In the present study, we found that the expression of deubiquitylase USP7 is significantly increased during human terminal erythroid differentiation. We further showed that interfering with USP7 function, either by short hairpin RNA-mediated knockdown or USP7-specific inhibitors, impaired human terminal erythroid differentiation due to decreased GATA1 level and that restoration of GATA1 levels rescued the differentiation defect. Mechanistically, USP7 deficiency led to a decreased GATA1 protein level that could be reversed by proteasome inhibitors. Furthermore, USP7 interacts directly with GATA1 and catalyzes the removal of K48-linked poly ubiquitylation chains conjugated onto GATA1, thereby stabilizing GATA1 protein. Collectively, our findings have identified an important role of a deubiquitylase in human terminal erythroid differentiation by stabilizing GATA1, the master regulator of erythropoiesis.

TSPO2 (translocator protein 2) is a transmembrane protein specifically expressed in late erythroblasts and has been postulated to mediate intracellular redistribution of cholesterol. We identified TSPO2 as the causative gene for the HK (high-K+) trait with immature red cell phenotypes in dogs and investigated the effects of the TSPO2 defects on erythropoiesis in HK dogs with the TSPO2 mutation and Tspo2 knockout (Tspo2-/-) mouse models. Bone marrow-derived erythroblasts from HK dogs showed increased binucleated and apoptotic cells at various stages of maturation and shed large nuclei with incomplete condensation when cultured in the presence of erythropoietin, indicating impaired maturation and cytokinesis. The canine TSPO2 induces cholesterol accumulation in the endoplasmic reticulum and could thereby regulate cholesterol availability by changing intracellular cholesterol distribution in erythroblasts. Tspo2-/- mice consistently showed impaired cytokinesis with increased binucleated erythroblasts, resulting in compensated anemia, and their red cell membranes had increased Na,K-ATPase, resembling the HK phenotype in dogs. Tspo2-deficient mouse embryonic stem cell-derived erythroid progenitor (MEDEP) cells exhibited similar morphological defects associated with a cell-cycle arrest at the G2/M phase, resulting in decreased cell proliferation and had a depletion in intracellular unesterified and esterified cholesterol. When the terminal maturation was induced, Tspo2-/- MEDEP cells showed delays in hemoglobinization; maturation-associated phenotypic changes in CD44, CD71, and TER119 expression; and cell-cycle progression. Taken together, these findings imply that TSPO2 is essential for coordination of maturation and proliferation of erythroblasts during normal erythropoiesis.

Enucleation is a key event in mammalian erythropoiesis responsible for the generation of enucleated reticulocytes. Although progress is being made in developing mechanistic understanding of enucleation, our understanding of mechanisms for enucleation is still incomplete. The MAPK pathway plays diverse roles in biological processes, but its role in erythropoiesis has yet to be fully defined. Analysis of RNA-sequencing data revealed that the MAPK pathway is significantly upregulated during human terminal erythroid differentiation. The MAPK pathway consists of 3 major signaling cassettes: MEK/ERK, p38, and JNK. In the present study, we show that among these 3 cassettes, only ERK was significantly upregulated in late-stage human erythroblasts. The increased expression of ERK along with its increased phosphorylation suggests a potential role for ERK activation in enucleation. To explore this hypothesis, we treated sorted populations of human orthochromatic erythroblasts with the MEK/ERK inhibitor U0126 and found that U0126 inhibited enucleation. In contrast, inhibitors of either p38 or JNK had no effect on enucleation. Mechanistically, U0126 selectively inhibited formation/accumulation of cytoplasmic vesicles and endocytosis of the transferrin receptor without affecting chromatin condensation, nuclear polarization, or enucleosome formation. Treatment with vacuolin-1 that induces vacuole formation partially rescued the blockage of enucleation by U0126. Moreover, phosphoproteomic analysis revealed that inactivation of the ERK pathway led to downregulation of the endocytic recycling pathway. Collectively, our findings uncovered a novel role of ERK activation in human erythroblast enucleation by modulating vesicle formation and have implications for understanding anemia associated with defective enucleation.

Spectrin tetramers of the membranes of enucleated mammalian erythrocytes play a critical role in red blood cell survival in circulation. One of the spectrins, αI, emerged in mammals with enucleated red cells after duplication of the ancestral α-spectrin gene common to all animals. The neofunctionalized αI-spectrin has moderate affinity for βI-spectrin, whereas αII-spectrin, expressed in nonerythroid cells, retains ancestral characteristics and has a 10-fold higher affinity for βI-spectrin. It has been hypothesized that this adaptation allows for rapid make and break of tetramers to accommodate membrane deformation. We have tested this hypothesis by generating mice with high-affinity spectrin tetramers formed by exchanging the site of tetramer formation in αI-spectrin (segments R0 and R1) for that of αII-spectrin. Erythrocytes with αIIβI presented normal hematologic parameters yet showed increased thermostability, and their membranes were significantly less deformable; under low shear forces, they displayed tumbling behavior rather than tank treading. The membrane skeleton is more stable with αIIβI and shows significantly less remodeling under deformation than red cell membranes of wild-type mice. These data demonstrate that spectrin tetramers undergo remodeling in intact erythrocytes and that this is required for the normal deformability of the erythrocyte membrane. We conclude that αI-spectrin represents evolutionary optimization of tetramer formation: neither higher-affinity tetramers (as shown here) nor lower affinity (as seen in hemolytic disease) can support the membrane properties required for effective tissue oxygenation in circulation.

Whole-exome sequencing is increasingly used for diagnosis and identification of appropriate therapies in patients. Here, we present the case of a 3-yr-old male with a lifelong and severe transfusion-dependent anemia of unclear etiology, despite an extensive clinical workup. Given the difficulty of making the diagnosis and the potential side effects from performing interventions in patients with a congenital anemia of unknown etiology, we opted to perform whole-exome sequencing on the patient and his parents. This resulted in the identification of homozygous loss-of-function mutations in the EPB41 gene, encoding erythrocyte protein band 4.1, which therefore causes a rare and severe form of hereditary elliptocytosis in the patient. Based on prior clinical experience in similar patients, a surgical splenectomy was performed that resulted in subsequent transfusion independence in the patient. This case illustrates how whole-exome sequencing can lead to accurate diagnoses (and exclusion of diagnoses where interventions, such as splenectomy, would be contraindicated), thereby resulting in appropriate and successful therapeutic intervention-a major goal of precision medicine.

The metabolic state of quiescent hematopoietic stem cells (HSCs) is an important regulator of self-renewal, but it is unclear whether or how metabolic parameters contribute to HSC lineage specification and commitment. Here, we show that the commitment of human and murine HSCs to the erythroid lineage is dependent upon glutamine metabolism. HSCs require the ASCT2 glutamine transporter and active glutamine metabolism for erythroid specification. Blocking this pathway diverts EPO-stimulated HSCs to differentiate into myelomonocytic fates, altering in vivo HSC responses and erythroid commitment under stress conditions such as hemolytic anemia. Mechanistically, erythroid specification of HSCs requires glutamine-dependent de novo nucleotide biosynthesis. Exogenous nucleosides rescue erythroid commitment of human HSCs under conditions of limited glutamine catabolism, and glucose-stimulated nucleotide biosynthesis further enhances erythroid specification. Thus, the availability of glutamine and glucose to provide fuel for nucleotide biosynthesis regulates HSC lineage commitment under conditions of metabolic stress.

Artemisinins are the cornerstone of anti-malarial drugs. Emergence and spread of resistance to them raises risk of wiping out recent gains achieved in reducing worldwide malaria burden and threatens future malaria control and elimination on a global level. Genome-wide association studies (GWAS) have revealed parasite genetic loci associated with artemisinin resistance. However, there is no consensus on biochemical targets of artemisinin. Whether and how these targets interact with genes identified by GWAS, remains unknown. Here we provide biochemical and cellular evidence that artemisinins are potent inhibitors of Plasmodium falciparum phosphatidylinositol-3-kinase (PfPI3K), revealing an unexpected mechanism of action. In resistant clinical strains, increased PfPI3K was associated with the C580Y mutation in P. falciparum Kelch13 (PfKelch13), a primary marker of artemisinin resistance. Polyubiquitination of PfPI3K and its binding to PfKelch13 were reduced by the PfKelch13 mutation, which limited proteolysis of PfPI3K and thus increased levels of the kinase, as well as its lipid product phosphatidylinositol-3-phosphate (PI3P). We find PI3P levels to be predictive of artemisinin resistance in both clinical and engineered laboratory parasites as well as across non-isogenic strains. Elevated PI3P induced artemisinin resistance in absence of PfKelch13 mutations, but remained responsive to regulation by PfKelch13. Evidence is presented for PI3P-dependent signalling in which transgenic expression of an additional kinase confers resistance. Together these data present PI3P as the key mediator of artemisinin resistance and the sole PfPI3K as an important target for malaria elimination.

Burst-forming unit-erythroid (BFU-E) and colony-forming unit-erythroid (CFU-E) cells are erythroid progenitors traditionally defined by colony assays. We developed a flow cytometry-based strategy for isolating human BFU-E and CFU-E cells based on the changes in expression of cell surface markers during in vitro erythroid cell culture. BFU-E and CFU-E are characterized by CD45(+)GPA(-)IL-3R(-)CD34(+)CD36(-)CD71(low) and CD45(+)GPA(-)IL-3R(-)CD34(-)CD36(+)CD71(high) phenotypes, respectively. Colony assays validated phenotypic assignment giving rise to BFU-E and CFU-E colonies, both at a purity of ∼90%. The BFU-E colony forming ability of CD45(+)GPA(-)IL-3R(-)CD34(+)CD36(-)CD71(low) cells required stem cell factor and erythropoietin, while the CFU-E colony forming ability of CD45(+)GPA(-)IL-3R(-)CD34(-)CD36(+)CD71(high) cells required only erythropoietin. Bioinformatic analysis of the RNA-sequencing data revealed unique transcriptomes at each differentiation stage. The sorting strategy was validated in uncultured primary cells isolated from bone marrow, cord blood, and peripheral blood, indicating that marker expression is not an artifact of in vitro cell culture, but represents an in vivo characteristic of erythroid progenitor populations. The ability to isolate highly pure human BFU-E and CFU-E progenitors will enable detailed cellular and molecular characterization of these distinct progenitor populations and define their contribution to disordered erythropoiesis in inherited and acquired hematologic disease. Our data provides an important resource for future studies of human erythropoiesis.