High-level embryonic globin production with efficient erythroid differentiation from a K562 erythroleukemia cell line

Linking regulatory variants to target genes by integrating single-cell multiome methods and genomic distance

Methods that analyze single-cell paired RNA-seq and ATAC-seq multiome data have shown great promise in linking regulatory elements to genes. However, existing methods differ in their modeling assumptions and approaches to account for biological and technical noise-leading to low concordance in their linking scores-and do not capture the effects of genomic distance. We propose pgBoost, an integrative modeling framework that trains a non-linear combination of existing linking strategies (including genomic distance) on fine-mapped eQTL data to assign a probabilistic score to each candidate SNP-gene link. We applied pgBoost to single-cell multiome data from 85k cells representing 6 major immune/blood cell types. pgBoost attained higher enrichment for fine-mapped eSNP-eGene pairs (e.g. 21x at distance >10kb) than existing methods (1.2-10x; p-value for difference = 5e-13 vs. distance-based method and < 4e-35 for each other method), with larger improvements at larger distances (e.g. 35x vs. 0.89-6.6x at distance >100kb; p-value for difference < 0.002 vs. each other method). pgBoost also outperformed existing methods in enrichment for CRISPR-validated links (e.g. 4.8x vs. 1.6-4.1x at distance >10kb; p-value for difference = 0.25 vs. distance-based method and < 2e-5 for each other method), with larger improvements at larger distances (e.g. 15x vs. 1.6-2.5x at distance >100kb; p-value for difference < 0.009 for each other method). Similar improvements in enrichment were observed for links derived from Activity-By-Contact (ABC) scores and GWAS data. We further determined that restricting pgBoost to features from a focal cell type improved the identification of SNP-gene links relevant to that cell type. We highlight several examples where pgBoost linked fine-mapped GWAS variants to experimentally validated or biologically plausible target genes that were not implicated by other methods. In conclusion, a non-linear combination of linking strategies, including genomic distance, improves power to identify target genes underlying GWAS associations.

Variable characteristics overlooked in human K-562 leukemia cell lines with a common signature

K-562 is a well-known in vitro cellular model that represents human leukemia cell lines. Although the K-562 cells have been extensively characterized, there are inconsistencies in the data across publications, showing the presence of multiple K-562 cell lines. This suggests that analyzing a single K-562 cell line is insufficient to provide reliable reference data. In this study, we compared three K-562 cell lines with different IDs (RCB0027, RCB1635, and RCB1897) to investigate the fundamental characteristics of K-562 cells. Amplifications of the BCR-ABL1 fusion gene and at 13q31 were detected in all three cell lines, whereas each genome exhibited distinctive features of sequence variants and loss of heterozygosity. This implies that each K-562 cell line can be characterized by common and unique features through a comparison of multiple K-562 cell lines. Variations in transcriptome profiles and hemoglobin synthesis were also observed among the three cell lines, indicating that they should be considered sublines that have diverged from the common ancestral K-562 despite no changes from the original cell name. This leads to unintentional differences in genotypes and/or phenotypes among cell lines that share the same name. These data show that characterizing a single K-562 cell line does not necessarily provide data that are applicable to other K-562 cells. In this context, it is essential to modify cell names in accordance with changes in characteristics during cell culture. Furthermore, our data could serve as a reference for evaluating other K-562 sublines, facilitating the discovery of new K-562 sublines with distinct characteristics. This approach results in the accumulation of K-562 sublines with diverged characteristics and expands the options available, which may help in selecting the most suitable K-562 subline for each experiment.

Cellular and animal models for the investigation of β-thalassemia

β-Thalassemia is a genetic form of anemia due to mutations in the β-globin gene, that leads to ineffective and extramedullary erythropoiesis, abnormal red blood cells and secondary iron-overload. The severity of the disease ranges from mild to lethal anemia based on the residual levels of globins production. Despite being a monogenic disorder, the pathophysiology of β-thalassemia is multifactorial, with different players contributing to the severity of anemia and secondary complications. As a result, the identification of effective therapeutic strategies is complex, and the treatment of patients is still suboptimal. For these reasons, several models have been developed in the last decades to provide experimental tools for the study of the disease, including erythroid cell lines, cultures of primary erythroid cells and transgenic animals. Years of research enabled the optimization of these models and led to decipher the mechanisms responsible for globins deregulation and ineffective erythropoiesis in thalassemia, to unravel the role of iron homeostasis in the disease and to identify and validate novel therapeutic targets and agents. Examples of successful outcomes of these analyses include iron restricting agents, currently tested in the clinics, several gene therapy vectors, one of which was recently approved for the treatment of most severe patients, and a promising gene editing strategy, that has been shown to be effective in a clinical trial. This review provides an overview of the available models, discusses pros and cons, and the key findings obtained from their study.

Regulation of erythroid differentiation in K562 cells by the EPAS1-IRS2 axis under hypoxic conditions

Red blood cells (RBCs) produced in vitro have the potential to alleviate the worldwide demand for blood transfusion. Hematopoietic cell differentiation and proliferation are triggered by numerous cellular physiological processes, including low oxygen concentration (<5%). In addition, hypoxia inducible factor 2α (HIF-2α) and insulin receptor substrate 2 (IRS2) were found to be involved in the progression of erythroid differentiation. However, the function of the HIF-2α-IRS2 axis in the progression of erythropoiesis is not yet fully understood. Therefore, we used an in vitro model of erythropoiesis generated from K562 cells transduced with shEPAS1 at 5% O₂ in the presence or absence of the IRS2 inhibitor NT157. We observed that erythroid differentiation was accelerated in K562 cells by hypoxia. Conversely, knockdown of EPAS1 expression reduced IRS2 expression and erythroid differentiation. Intriguingly, inhibition of IRS2 could impair the progression of hypoxia-induced erythropoiesis without affecting EPAS1 expression. These findings indicated that the EPAS1-IRS2 axis may be a crucial pathway that regulates erythropoiesis and that drugs targeting this pathway may become promising agents for promoting erythroid differentiation.

Metabolic sensor O-GlcNAcylation regulates erythroid differentiation and globin production via BCL11A

Background

Human erythropoiesis is a tightly regulated, multistep process encompassing the differentiation of hematopoietic stem cells (HSCs) toward mature erythrocytes. Cellular metabolism is an important regulator of cell fate determination during the differentiation of HSCs. However, how O-GlcNAcylation, a posttranslational modification of proteins that is an ideal metabolic sensor, contributes to the commitment of HSCs to the erythroid lineage and to the terminal erythroid differentiation has not been addressed.

Methods

Cellular O-GlcNAcylation was manipulated using small molecule inhibition or CRISPR/Cas9 manipulation of catalyzing enzyme O-GlcNAc transferase (OGT) and removing enzyme O-GlcNAcase (OGA) in two cell models of erythroid differentiation, starting from: (i) human umbilical cord blood-derived CD34⁺ hematopoietic stem/progenitor cells (HSPCs) to investigate the erythroid lineage specification and differentiation; and (ii) human-derived erythroblastic leukemia K562 cells to investigate the terminal differentiation. The functional and regulatory roles of O-GlcNAcylation in erythroid differentiation, maturation, and globin production were investigated, and downstream signaling was delineated.

Results

First, we observed that two-step inhibition of OGT and OGA, which were established from the observed dynamics of O-GlcNAc level along the course of differentiation, promotes HSPCs toward erythroid differentiation and enucleation, in agreement with an upregulation of a multitude of erythroid-associated genes. Further studies in the efficient K562 model of erythroid differentiation confirmed that OGA inhibition and subsequent hyper-O-GlcNAcylation enhance terminal erythroid differentiation and affect globin production. Mechanistically, we found that BCL11A is a key mediator of O-GlcNAc-driven erythroid differentiation and β- and α-globin production herein. Additionally, analysis of biochemical contents using synchrotron-based Fourier transform infrared (FTIR) spectroscopy showed unique metabolic fingerprints upon OGA inhibition during erythroid differentiation, supporting that metabolic reprogramming plays a part in this process.

Conclusions

The evidence presented here demonstrated the novel regulatory role of O-GlcNAc/BCL11A axis in erythroid differentiation, maturation, and globin production that could be important in understanding erythropoiesis and hematologic disorders whose etiology is related to impaired erythroid differentiation and hemoglobinopathies. Our findings may lay the groundwork for future clinical applications toward an ex vivo production of functional human reticulocytes for transfusion from renewable cell sources, i.e., HSPCs and pluripotent stem cells.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13287-022-02954-5.

Cabozantinib promotes erythroid differentiation in K562 erythroleukemia cells through global changes in gene expression and JNK activation

Cabozantinib is a potent tyrosine kinase inhibitor with multiple targets including MET, VEGFR2, RET, KIT, and FLT3. Cabozantinib is widely used for the treatment of medullary thyroid cancer and renal cell carcinoma. We recently suggested cabozantinib as a potential therapeutic alternative for acute myeloid leukemia (AML) patients with FLT3-internal tandem duplication (FLT3-ITD). Here, we report that cabozantinib can promote differentiation in erythroid leukemia cells. We found that K562 erythroid leukemia cells treated with 1 μM cabozantinib for 72 h underwent erythroid lineage differentiation. Transcriptomic analysis revealed that various pathways associated with heme biosynthesis, hemoglobin production, and GATA1 targets were upregulated, whereas cell survival pathways were downregulated. Further examination revealed that cabozantinib-induced erythroid differentiation is at least in part regulated by JNK activation and phosphorylation. Levels of phosphorylated BCR-ABL, AKT, STAT5, ERK, and p38 also decreased following cabozantinib treatment. Therefore, we indicate that cabozantinib has dual functions. First, it induces K562 cell differentiation toward the erythroid lineage by upregulating heme biosynthesis, globin synthesis, and erythroid-associated reactions. Second, cabozantinib inhibits K562 cell proliferation by inhibiting the phosphorylation of BCR-ABL and the downstream MAPK, PI3K-AKT, and JAK-STAT signaling pathways.

The opposing roles of the mTOR signaling pathway in different phases of human umbilical cord blood-derived CD34+ cell erythropoiesis

As an indispensable, even lifesaving practice, red blood cell (RBC) transfusion is challenging due to several issues, including supply shortage, immune incompatibility, and blood-borne infections since donated blood is the only source of RBCs. Although large-scale in vitro production of functional RBCs from human stem cells is a promising alternative, so far, no such system has been reported to produce clinically transfusable RBCs due to the poor understanding of mechanisms of human erythropoiesis, which is essential for the optimization of in vitro erythrocyte generation system. We previously reported that inhibition of mammalian target of rapamycin (mTOR) signaling significantly decreased the percentage of erythroid progenitor cells in the bone marrow of wild-type mice. In contrast, rapamycin treatment remarkably improved terminal maturation of erythroblasts and anemia in a mouse model of β-thalassemia. In the present study, we investigated the effect of mTOR inhibition with rapamycin from different time points on human umbilical cord blood-derived CD34⁺ cell erythropoiesis in vitro and the underlying mechanisms. Our data showed that rapamycin treatment significantly suppressed erythroid colony formation in the commitment/proliferation phase of erythropoiesis through inhibition of cell-cycle progression and proliferation. In contrast, during the maturation phase of erythropoiesis, mTOR inhibition dramatically promoted enucleation and mitochondrial clearance by enhancing autophagy. Collectively, our results suggest contrasting roles for mTOR in regulating different phases of human erythropoiesis.

βT87Q-Globin Gene Therapy Reduces Sickle Hemoglobin Production, Allowing for Ex Vivo Anti-sickling Activity in Human Erythroid Cells

Lentiviral addition of βT87Q-globin, a modified β-globin with an anti-sickling mutation, is currently being used in gene therapy trials for sickle cell disease (SCD) and β-thalassemia patients. βT87Q-globin interferes with sickle hemoglobin (HbS) polymerization. Here, we generated the SCD mutation in an immortalized human erythroid cell line (HUDEP-2) to investigate the anti-sickling activity of βT87Q-globin. Sickle HUDEP-2 (sHUDEP-2) cells produced robust HbS after differentiation and sickled under deoxygenated conditions, comparable with SCD CD34⁺ progeny. Lentiviral transduction provided 9.5–26.8 pg/cell βT87Q-globin (R² = 0.83) in a vector copy number (VCN)-dependent manner, resulting in a significant reduction of sickling ratios (R² = 0.92). Interestingly, βT87Q-globin transduction markedly reduced endogenous β^S-globin (R² = 0.84) to an undetectable level (0.4–16.8 pg/cell) in sHUDEP-2 cells, as well as endogenous β-globin in human CD34⁺ cell-derived erythroid cells. RNA sequencing (RNA-seq) analysis with βT87Q-transduced sHUDEP-2 and human CD34⁺-derived cells revealed activation of inflammation- and proliferation-related programs, suggesting minimal changes in background gene expression except for βT87Q-globin expression and endogenous β/β^S-globin suppression. In summary, using sHUDEP-2 and CD34⁺-derived cells, we demonstrated that lentiviral addition of βT87Q-globin strongly reduced endogenous β-/β^S-globin expression, resulting in an anti-sickling effect. Our findings should be helpful to understand the anti-sickling effects of therapeutic genes in SCD gene therapy.

Biallelic correction of sickle cell disease-derived induced pluripotent stem cells (iPSCs) confirmed at the protein level through serum-free iPS-sac/erythroid differentiation

New technologies of induced pluripotent stem cells (iPSCs) and genome editing have emerged, allowing for the development of autologous transfusion therapies. We previously demonstrated definitive β‐globin production from human embryonic stem cell (hESC)‐derived erythroid cell generation via hemangioblast‐like ES‐sacs. In this study, we demonstrated normal β‐globin protein production from biallelic corrected sickle cell disease (SCD) iPSCs. We optimized our ES/iPS‐sac method for feeder cell‐free hESC maintenance followed by serum‐free ES‐sac generation, which is preferred for electroporation‐based genome editing. Surprisingly, the optimized protocol improved yields of ES‐sacs (25.9‐fold), hematopoietic‐like spherical cells (14.8‐fold), and erythroid cells (5.8‐fold), compared with our standard ES‐sac generation. We performed viral vector‐free gene correction in SCD iPSCs, resulting in one clone with monoallelic and one clone with biallelic correction, and using this serum‐free iPS‐sac culture, corrected iPSC‐generated erythroid cells with normal β‐globin, confirmed at DNA and protein levels. Our serum‐free ES/iPS‐sac protocol with gene correction will be useful to develop regenerative transfusion therapies for SCD.

Robust generation of erythroid and multilineage hematopoietic progenitors from human iPSCs using a scalable monolayer culture system

One of the most promising objectives of clinical hematology is to derive engraftable autologous hematopoietic stem cells (HSCs) from human induced pluripotent stem cells (iPSCs). Progress in translating iPSC technologies to the clinic relies on the availability of scalable differentiation methodologies. In this study, human iPSCs were differentiated for 21 days using STEMdiff™, a monolayer-based approach for hematopoietic differentiation of human iPSCs that requires no replating, co-culture or embryoid body formation. Both hematopoietic and non-hematopoietic cells were functionally characterized throughout differentiation. In the hematopoietic fraction, an early transient population of primitive CD235a⁺ erythroid progenitor cells first emerged, followed by hematopoietic progenitors with multilineage differentiation activity in vitro but no long-term engraftment potential in vivo. In later stages of differentiation, a nearly exclusive production of definitive erythroid progenitors was observed. In the non-hematopoietic fraction, we identified a prevalent population of mesenchymal stromal cells and limited arterial vascular endothelium (VE), suggesting that the cellular constitution of the monolayer may be inadequate to support the generation of HSCs with durable repopulating potential. Quantitative modulation of WNT/β-catenin and activin/nodal/TGFβ signaling pathways with CHIR/SB molecules during differentiation enhanced formation of arterial VE, definitive multilineage and erythroid progenitors, but was insufficient to orchestrate the generation of engrafting HSCs. Overall, STEMdiff™ provides a clinically-relevant and readily adaptable platform for the generation of erythroid and multilineage hematopoietic progenitors from human pluripotent stem cells.