Abstract:Abstract Reticulocytosis in sickle cell disease (SCD) is driven by tissue hypoxia from hemolytic anemia and vascular occlusion. Gene expression changes caused by hypoxia and other factors during reticulocytosis may impact SCD outcomes. We detected 1226 differentially expressed genes in SCD reticulocyte transcriptome compared to normal Black controls. To assess the role of hypoxia-mediating HIFs from other regulation of changes of the SCD reticulocyte transcriptome, we compared differential expression in SCD to that in Chuvash erythrocytosis (CE), a disorder characterized by constitutive upregulation of HIFs in normoxia. Of the SCD differentially expressed genes, 28% were shared between CE and SCD and thus classified as HIF-mediated. The HIF-mediated changes were generally in genes promoting erythroid maturation. We found that genes encoding the response to endoplasmic reticulum stress generally lacked HIF mediation. We then investigated the clinical correlation of erythroid gene expression for the 1226 differentially expressed genes detected in SCD reticulocytes, using clinical measures and gene expression data previously profiled in peripheral blood mononuclear cells (PBMCs) of 157 SCD patients at the University of Illinois at Chicago (UIC). Normal PBMCs contain only a small number of erythroid progenitors, but in SCD or CE PBMCs the erythroid transcriptome is enriched due to elevated circulating erythroid progenitors from heightened erythropoiesis (PMID: 32399971). We applied deconvolution analysis to assess the clinical correlation of erythroid gene expression, using a 16-gene expression signature of erythroid progenitors previously identified in SCD PBMCs. Deconvolution analysis uses the proportion of cell/tissue or specific marker genes (here the erythroid specific 16-gene signature) to dissect gene expression variation in biological samples with cell/tissue type heterogeneity. We correlated, in the 157 UIC patients, erythroid gene expression with i) degree of anemia as indicated by hemoglobin concentration, ii) vaso-occlusive severe pain episodes per year, and iii) degree of hemolysis measured by a hemolysis index. The analysis identified 231 genes associated with at least one of the complications. Increased expression of 40 erythroid specific genes, including 15 HIF-mediated genes, was associated with all three complications. These 40 genes are all upregulated in SCD reticulocytes and correlated with low hemoglobin concentration, frequent severe pain episodes, and high hemolysis index, suggesting that these manifestations may share a relationship to stress erythropoiesis-driven transcriptional activity. Expression quantitative trait loci (eQTL) contain genetic polymorphisms that associate with gene expression level, which can be viewed as a natural experiment to investigate the causal relations between gene expression change and phenotypic outcomes. To assess the causal effect of erythroid gene expression, we tested association between erythroid eQTL and the clinical manifestations in 906 SCD patients from the Walk-PHaSST and PUSH cohorts. We first mapped erythroid eQTL in the 157 UIC patients, who were previously genotyped by array, applying deconvolution algorithm on the same PBMC data for the 1226 differential genes in SCD reticulocytes, and detected 54 distinct eQTL for 30 genes at 5% false discovery rate. After adjusting for multiple comparisons, we found that the C allele of rs16911905, located in the β-globin cluster and associated with increased erythroid expression of HBD (encodes δ-globin of hemoglobin A 2), significantly correlated with lower hemoglobin concentration (β=-0.064, 95% CI -0.092 - -0.036, P=6.7×10 -6). The C allele was also associated with higher hemolytic rate (P=0.031), less frequent pain episodes (P=0.045), and increased erythroid expression of HBB here encoding sickle β-globin (P=5.1x10 -5). The association of the C allele with lower hemoglobin concentration was then validated in 242 patients from the UIC cohort (β=-0.071, 95% CI -0.13 - -0.011, P=0.023), as was the trend of association with higher hemolytic rate (P=0.0031) and less pain episodes (P=0.034). Our findings reveal HIF- and non-HIF-mediated genes in SCD stress erythropoiesis, and identify novel clinical associations for a HBD eQTL. Our study highlights the correlation of altered erythroid gene expression with SCD hemolytic and vaso-occlusive manifestations. Disclosures Saraf: Global Blood Therapeutics: Membership on an entity's Board of Directors or advisory committees, Research Funding; Novartis: Membership on an entity's Board of Directors or advisory committees, Research Funding; Pfizer: Research Funding. Gordeuk: Modus Therapeutics: Consultancy; Novartis: Research Funding; Incyte: Research Funding; Emmaus: Consultancy, Research Funding; Global Blood Therapeutics: Consultancy, Research Funding; CSL Behring: Consultancy.

Reduced GATA1 levels are associated with ineffective erythropoiesis in sickle cell anemia

Senp1-Mediated Gata1 Desumoylation is Critical for Definitive Erythropoiesis

Gene expression changes in sickle cell reticulocytes and their clinical associations

Low expression of stem cell antigen-1 on mouse haematopoietic precursors is associated with erythroid differentiation

Transcriptomic profiling of an ex vivo model of erythropoiesis uncovers potential genes and pathways involved into the pathogenesis of sickle cell disease

HIF-Mediated and Non-HIF-Mediated Differential Gene Expressions in Sickle Cell Reticulocyte and Their Impact on Clinical Manifestations

Single-cell analyses demonstrate that a heme-GATA1 feedback loop regulates red cell differentiation

Aberrant GPA expression and regulatory function of red blood cells in sickle cell disease.

Impaired DNA damage response and inflammatory signalling underpins hematopoietic stem cell defects in Gata2 haploinsufficiency

Genetic Association of Clinical Complications for Genes Differentially Expressed in Reticulocytes of Sickle Cell Anemia

Complete absence of GLUT1 does not impair human terminal erythroid differentiation

Evaluation of GBT021601 As a Therapeutic Agent to Restore Bone Marrow Health and Effective Erythropoiesis in a Sickle Mouse Model

Gata-1, -2 and -3 Genes Expression in Bone Marrow Microenviroment with Chronic Aplastic Anemia

Pathological Role of Halted Proteasome Machinery in Sickle Cell Disease

An activin receptor IIA ligand trap corrects ineffective erythropoiesis in β-thalassemia

Elevated CDKN1A (P21) mediates β-thalassemia erythroid apoptosis, but its loss does not improve β-thalassemic erythropoiesis

Ineffective erythropoiesis and its treatment

Structural and Functional Insight of Sphingosine 1-Phosphate-mediated Pathogenic Metabolic Reprogramming in Sickle Cell Disease

Deubiquitinase Usp7 Regulates Erythroid Development Via Deubiquinating Gata1

The Novel GATA1-interacting Protein HES6 is an Essential Transcriptional Cofactor for Human Erythropoiesis.