Identification of gene networks associated with erythroid differentiation

Comparison of gene expression profiles between human erythroid cells derived from fetal liver and adult peripheral blood

Background

A key event in human development is the establishment of erythropoietic progenitors in the bone marrow, which is accompanied by a fetal-to-adult switch in hemoglobin expression. Understanding of this event could lead to medical application, notably treatment of sickle cell disease and β-thalassemia. The changes in gene expression of erythropoietic progenitor cells as they migrate from the fetal liver and colonize the bone marrow are still rather poorly understood, as primary fetal liver (FL) tissues are difficult to obtain.

Methods

We obtained human FL tissue and adult peripheral blood (AB) samples from Thai subjects. Primary CD34⁺ cells were cultured in vitro in a fetal bovine serum-based culture medium. After 8 days of culture, erythroid cell populations were isolated by flow cytometry. Gene expression in the FL- and AB-derived cells was studied by Affymetrix microarray and reverse-transcription quantitative PCR. The microarray data were combined with that from a previous study of human FL and AB erythroid development, and meta-analysis was performed on the combined dataset.

Results

FL erythroid cells showed enhanced proliferation and elevated fetal hemoglobin relative to AB cells. A total of 1,391 fetal up-regulated and 329 adult up-regulated genes were identified from microarray data generated in this study. Five hundred ninety-nine fetal up-regulated and 284 adult up-regulated genes with reproducible patterns between this and a previous study were identified by meta-analysis of the combined dataset, which constitute a core set of genes differentially expressed between FL and AB erythroid cells. In addition to these core genes, 826 and 48 novel genes were identified only from data generated in this study to be FL up- and AB up-regulated, respectively. The in vivo relevance for some of these novel genes was demonstrated by pathway analysis, which showed novel genes functioning in pathways known to be important in proliferation and erythropoiesis, including the mitogen-activated protein kinase (MAPK) and the phosphatidyl inositol 3 kinase (PI3K)-Akt pathways.

Discussion

The genes with upregulated expression in FL cells, which include many novel genes identified from data generated in this study, suggest that cellular proliferation pathways are more active in the fetal stage. Erythroid progenitor cells may thus undergo a reprogramming during ontogenesis in which proliferation is modulated by changes in expression of key regulators, primarily MYC, and others including insulin-like growth factor 2 mRNA-binding protein 3 (IGF2BP3), neuropilin and tolloid-like 2 (NETO2), branched chain amino acid transaminase 1 (BCAT1), tenascin XB (TNXB) and proto-oncogene, AP-1 transcription factor subunit (JUND). This reprogramming may thus be necessary for acquisition of the adult identity and switching of hemoglobin expression.

Global transcriptome analysis for identification of interactions between coding and noncoding RNAs during human erythroid differentiation

Studies on coding genes, miRNAs, and lncRNAs during erythroid development have been performed in recent years. However, analysis focusing on the integration of the three RNA types has yet to be done. In the present study, we compared the dynamics of coding genes, miRNA, and lncRNA expression profiles. To explore dynamic changes in erythropoiesis and potential mechanisms that control these changes in the transcriptome level, we took advantage of high throughput sequencing technologies to obtain transcriptome data from cord blood hematopoietic stem cells and the following four erythroid differentiation stages, as well as from mature red blood cells. Results indicated that lncRNAs were promising cell marker candidates for erythroid differentiation. Clustering analysis classified the differentially expressed genes into four subtypes that corresponded to dynamic changes during stemness maintenance, mid-differentiation, and maturation. Integrated analysis revealed that noncoding RNAs potentially participated in controlling blood cell maturation, and especially associated with heme metabolism and responses to oxygen species and DNA damage. These regulatory interactions were displayed in a comprehensive network, thereby inferring correlations between RNAs and their associated functions. These data provided a substantial resource for the study of normal erythropoiesis, which will permit further investigation and understanding of erythroid development and acquired erythroid disorders.

Novel methods for studying normal and disordered erythropoiesis

Erythropoiesis is a process during which multipotential hematopoietic stem cells proliferate, differentiate and eventually form mature erythrocytes. Interestingly, unlike most cell types, an important feature of erythropoiesis is that following each mitosis the daughter cells are morphologically and functionally different from the parent cell from which they are derived, demonstrating the need to study erythropoiesis in a stage-specific manner. This has been impossible until recently due to lack of methods for isolating erythroid cells at each distinct developmental stage. This review summarizes recent advances in the development of methods for isolating both murine and human erythroid cells and their applications. These methods provide powerful means for studying normal and impaired erythropoiesis associated with hematological disorders.

Characterization of transcription factor networks involved in umbilical cord blood CD34+ stem cells-derived erythropoiesis

Fetal stem cells isolated from umbilical cord blood (UCB) possess a great capacity for proliferation and differentiation and serve as a valuable model system to study gene regulation. Expanded knowledge of the molecular control of hemoglobin synthesis will provide a basis for rational design of therapies for β-hemoglobinopathies. Transcriptome data are available for erythroid progenitors derived from adult stem cells, however studies to define molecular mechanisms controlling globin gene regulation during fetal erythropoiesis are limited. Here, we utilize UCB-CD34⁺ stem cells induced to undergo erythroid differentiation to characterize the transcriptome and transcription factor networks (TFNs) associated with the γ/β-globin switch during fetal erythropoiesis. UCB-CD34⁺ stem cells grown in the one-phase liquid culture system displayed a higher proliferative capacity than adult CD34⁺ stem cells. The γ/β-globin switch was observed after day 42 during fetal erythropoiesis in contrast to adult progenitors where the switch occurred around day 21. To gain insights into transcription factors involved in globin gene regulation, microarray analysis was performed on RNA isolated from UCB-CD34⁺ cell-derived erythroid progenitors harvested on day 21, 42, 49 and 56 using the HumanHT-12 Expression BeadChip. After data normalization, Gene Set Enrichment Analysis identified transcription factors (TFs) with significant changes in expression during the γ/β-globin switch. Forty-five TFs were silenced by day 56 (Profile-1) and 30 TFs were activated by day 56 (Profile-2). Both GSEA datasets were analyzed using the MIMI Cytoscape platform, which discovered TFNs centered on KLF4 and GATA2 (Profile-1) and KLF1 and GATA1 for Profile-2 genes. Subsequent shRNA studies in KU812 leukemia cells and human erythroid progenitors generated from UCB-CD34⁺ cells supported a negative role of MAFB in γ-globin regulation. The characteristics of erythroblasts derived from UCB-CD34⁺ stem cells including prolonged γ-globin expression combined with unique TFNs support novel mechanisms controlling the γ/β-globin switch during UCB-derived erythropoiesis.

Developmental transcriptome analysis of human erythropoiesis

To globally survey the changes in transcriptional landscape during terminal erythroid differentiation, we performed RNA sequencing (RNA-seq) on primary human CD34(+) cells after ex vivo differentiation from the earliest into the most mature erythroid cell stages. This analysis identified thousands of novel intergenic and intronic transcripts as well as novel alternative transcript isoforms. After rigorous data filtering, 51 (presumptive) novel protein-coding transcripts, 5326 long and 679 small non-coding RNA candidates remained. The analysis also revealed two clear transcriptional trends during terminal erythroid differentiation: first, the complexity of transcript diversity was predominantly achieved by alternative splicing, and second, splicing junctional diversity diminished during erythroid differentiation. Finally, 404 genes that were not known previously to be differentially expressed in erythroid cells were annotated. Analysis of the most extremely differentially expressed transcripts revealed that these gene products were all closely associated with hematopoietic lineage differentiation. Taken together, this study will serve as a comprehensive platform for future in-depth investigation of human erythroid development that, in turn, may reveal new insights into multiple layers of the transcriptional regulatory hierarchy that controls erythropoiesis.

Global transcriptome analyses of human and murine terminal erythroid differentiation

We recently developed fluorescence-activated cell sorting (FACS)-based methods to purify morphologically and functionally discrete populations of cells, each representing specific stages of terminal erythroid differentiation. We used these techniques to obtain pure populations of both human and murine erythroblasts at distinct developmental stages. RNA was prepared from these cells and subjected to RNA sequencing analyses, creating unbiased, stage-specific transcriptomes. Tight clustering of transcriptomes from differing stages, even between biologically different replicates, validated the utility of the FACS-based assays. Bioinformatic analyses revealed that there were marked differences between differentiation stages, with both shared and dissimilar gene expression profiles defining each stage within transcriptional space. There were vast temporal changes in gene expression across the differentiation stages, with each stage exhibiting unique transcriptomes. Clustering and network analyses revealed that varying stage-specific patterns of expression observed across differentiation were enriched for genes of differing function. Numerous differences were present between human and murine transcriptomes, with significant variation in the global patterns of gene expression. These data provide a significant resource for studies of normal and perturbed erythropoiesis, allowing a deeper understanding of mechanisms of erythroid development in various inherited and acquired erythroid disorders.

Transcriptomic and phospho-proteomic analyzes of erythroblasts expanded in vitro from normal donors and from patients with polycythemia vera

Erythropoiesis is a tightly regulated process which becomes decoupled from its normal differentiation program in patients with polycythemia vera (PV). Somatic mutations in JAK2 are commonly associated with this myeloid proliferative disorder. To gain insight into the molecular events that are required for abnormally developing erythroid cells to escape dependence on normal growth signals, we performed in vitro expansion of mature erythroblasts (ERY) from seven normal healthy donors and from seven polycythemic patients in the presence of IL3, EPO, SCF for 10, 11, or 13 days. Normal ERYs required exposure to the glucocorticoid dexamethasone (Dex) for expansion, while PV-derived ERYs expanded in the absence of dexamethasone. RNA expression profiling revealed enrichment of two known oncogenes, GPR56 and RAB4a, in PV-derived ERYs along with reduced expression levels of transcription factor TAL1 (ANOVA FDR < 0.05). While both normal and polycythemic-derived ERYs integrated signaling cascades for growth, they did so via different signaling pathways which are represented by their differential phospho-profiles. Our results show that normal ERYs displayed greater levels of phosphorylation of EGFR, PDGFRβ, TGFβ, and cKit, while PV-derived ERYs were characterized by increased phosphorylation of cytoplasmic kinases in the JAK/STAT, PI3K, and GATA1 pathways. Together these data suggest that PV erythroblast expansion and maturation may be maintained and enriched in the absence of dexamethasone through reduced TAL1 expression and by accessing additional signaling cascades. Members of this acquired repertoire may provide important insight into the pathogenesis of aberrant erythropoiesis in myeloproliferative neoplasms such as polycythemia vera.

Characterization of the transcriptome profiles related to globin gene switching during in vitro erythroid maturation

Background

The fetal and adult globin genes in the human β-globin cluster on chromosome 11 are sequentially expressed to achieve normal hemoglobin switching during human development. The pharmacological induction of fetal γ-globin (HBG) to replace abnormal adult sickle β^S-globin is a successful strategy to treat sickle cell disease; however the molecular mechanism of γ-gene silencing after birth is not fully understood. Therefore, we performed global gene expression profiling using primary erythroid progenitors grown from human peripheral blood mononuclear cells to characterize gene expression patterns during the γ-globin to β-globin (γ/β) switch observed throughout in vitro erythroid differentiation.

Results

We confirmed erythroid maturation in our culture system using cell morphologic features defined by Giemsa staining and the γ/β-globin switch by reverse transcription-quantitative PCR (RT-qPCR) analysis. We observed maximal γ-globin expression at day 7 with a switch to a predominance of β-globin expression by day 28 and the γ/β-globin switch occurred around day 21. Expression patterns for transcription factors including GATA1, GATA2, KLF1 and NFE2 confirmed our system produced the expected pattern of expression based on the known function of these factors in globin gene regulation. Subsequent gene expression profiling was performed with RNA isolated from progenitors harvested at day 7, 14, 21, and 28 in culture. Three major gene profiles were generated by Principal Component Analysis (PCA). For profile-1 genes, where expression decreased from day 7 to day 28, we identified 2,102 genes down-regulated > 1.5-fold. Ingenuity pathway analysis (IPA) for profile-1 genes demonstrated involvement of the Cdc42, phospholipase C, NF-Kβ, Interleukin-4, and p38 mitogen activated protein kinase (MAPK) signaling pathways. Transcription factors known to be involved in γ-and β-globin regulation were identified.

The same approach was used to generate profile-2 genes where expression was up-regulated over 28 days in culture. IPA for the 2,437 genes with > 1.5-fold induction identified the mitotic roles of polo-like kinase, aryl hydrocarbon receptor, cell cycle control, and ATM (Ataxia Telangiectasia Mutated Protein) signaling pathways; transcription factors identified included KLF1, GATA1 and NFE2 among others. Finally, profile-3 was generated from 1,579 genes with maximal expression at day 21, around the time of the γ/β-globin switch. IPA identified associations with cell cycle control, ATM, and aryl hydrocarbon receptor signaling pathways.

Conclusions

The transcriptome analysis completed with erythroid progenitors grown in vitro identified groups of genes with distinct expression profiles, which function in metabolic pathways associated with cell survival, hematopoiesis, blood cells activation, and inflammatory responses. This study represents the first report of a transcriptome analysis in human primary erythroid progenitors to identify transcription factors involved in hemoglobin switching. Our results also demonstrate that the in vitro liquid culture system is an excellent model to define mechanisms of global gene expression and the DNA-binding protein and signaling pathways involved in globin gene regulation.

Ex-vivo expansion of red blood cells: how real for transfusion in humans?

Blood transfusion is indispensable for modern medicine. In developed countries, the blood supply is adequate and safe but blood for alloimmunized patients is often unavailable. Concerns are increasing that donations may become inadequate in the future as the population ages prompting a search for alternative transfusion products. Improvements in culture conditions and proof-of-principle studies in animal models have suggested that ex-vivo expanded red cells may represent such a product. Compared to other cell therapies transfusion poses the unique challenge of requiring great cell doses (2.5×10(12) cells vs 10(7) cells). Although production of such cell numbers is theoretically possible, current technologies generate red cells in numbers sufficient only for safety studies. It is conceived that by the time these studies will be completed, technical barriers to mass cell production will have been eliminated making transfusion with ex-vivo generated red cells a reality.

Global gene expression analysis of human erythroid progenitors

Understanding the pattern of gene expression during erythropoiesis is crucial for a synthesis of erythroid developmental biology. Here, we isolated 4 distinct populations at successive erythropoietin-dependent stages of erythropoiesis, including the terminal, pyknotic stage. The transcriptome was determined using Affymetrix arrays. First, we demonstrated the importance of using defined cell populations to identify lineage and temporally specific patterns of gene expression. Cells sorted by surface expression profile not only express significantly fewer genes than unsorted cells but also demonstrate significantly greater differences in the expression levels of particular genes between stages than unsorted cells. Second, using standard software, we identified more than 1000 transcripts not previously observed to be differentially expressed during erythroid maturation, 13 of which are highly significantly terminally regulated, including RFXAP and SMARCA4. Third, using matched filtering, we identified 12 transcripts not previously reported to be continuously up-regulated in maturing human primary erythroblasts. Finally, using transcription factor binding site analysis, we identified potential transcription factors that may regulate gene expression during terminal erythropoiesis. Our stringent lists of differentially regulated and continuously expressed transcripts containing many genes with undiscovered functions in erythroblasts are a resource for future functional studies of erythropoiesis. Our Human Erythroid Maturation database is available at https://cellline.molbiol.ox.ac.uk/eryth/index.html. [corrected].

An antimicrobial peptide regulates tumor-associated macrophage trafficking via the chemokine receptor CCR2, a model for tumorigenesis

BACKGROUND

Tumor-associated macrophages (TAMs) constitute a significant part of infiltrating inflammatory cells that are frequently correlated with progression and poor prognosis of a variety of cancers. Tumor cell-produced human beta-defensin-3 (hBD-3) has been associated with TAM trafficking in oral cancer; however, its involvement in tumor-related inflammatory processes remains largely unknown.

METHODOLOGY

The relationship between hBD-3, monocyte chemoattractant protein-1 (MCP-1), TAMs, and CCR2 was examined using immunofluorescence microscopy in normal and oral carcinoma in situ biopsy specimens. The ability of hBD-3 to chemoattract host macrophages in vivo using a nude mouse model and analysis of hBD-3 on monocytic cell migration in vitro, applying a cross-desensitization strategy of CCR2 and its pharmacological inhibitor (RS102895), respectively, was also carried out.

CONCLUSIONS/FINDINGS

MCP-1, the most frequently expressed tumor cell-associated chemokine, was not produced by tumor cells nor correlated with the recruitment of macrophages in oral carcinoma in situ lesions. However, hBD-3 was associated with macrophage recruitment in these lesions and hBD-3-expressing tumorigenic cells induced massive tumor infiltration of host macrophages in nude mice. HBD-3 stimulated the expression of tumor-promoting cytokines, including interleukin-1alpha (IL-1alpha), IL-6, IL-8, CCL18, and tumor necrosis factor-alpha (TNF-alpha) in macrophages derived from human peripheral blood monocytes. Monocytic cell migration in response to hBD-3 was inhibited by cross-desensitization with MCP-1 and the specific CCR2 inhibitor, RS102895, suggesting that CCR2 mediates monocyte/macrophage migration in response to hBD-3. Collectively, these results indicate that hBD-3 utilizes CCR2 to regulate monocyte/macrophage trafficking and may act as a tumor cell-produced chemoattractant to recruit TAMs. This novel mechanism is the first evidence of an hBD molecule orchestrating an in vivo outcome and demonstrates the importance of the innate immune system in the development of tumors.