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.

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

Gene expression changes in sickle cell reticulocytes and their clinical associations

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

Gene-Gene Interactions And The Pathophysiology Of Sickle Cell Disease: Modeling The Effects Of Snps On Sickle Cell-Associated Vasoocclusive Events Using Classification And Regression Trees And Stochastic Gradient Boosting

Identification of genetic variants associated with clinical features of sickle cell disease

Polymorphisms (snps) in Multiple Genes of the Tgf-ß/Bmp Pathway Are Associated with a Global Measure of Sickle Cell Disease Severity.

Genetic contributions to pain modulation in sickle cell: A focus on single nucleotide polymorphisms

Precision Medicine For Sickle Cell Disease Through Whole Genome Sequencing

Leg Ulcers in Sickle Cell Anemia Are Associated with Laboratory Markers of Hemolysis and SNPs in KL and Genes of the TGF-β/BMP Pathway.

A polygenic score for acute vaso-occlusive pain in pediatric sickle cell disease

Differences in the Circulating Proteome in Individuals with versus without Sickle Cell Trait

Multiple-ancestry genome-wide association study identifies 27 loci associated with measures of hemolysis following blood storage

New Insights into Clinical Management for Sickle Cell Disease: Uncovering the Significant Pathways Affected by the Involvement of Sickle Cell Disease

Nontargeted Plasma Proteomic Analysis of Renal Disease and Pulmonary Hypertension in Patients with Sickle Cell Disease

Hypoxic Response Contributes to Altered Gene Expression and Precapillary Pulmonary Hypertension in Patients with Sickle Cell Disease.

Hemolytic Biomarkers Predict Adhesiveness of Sickle Blood Cells in a Clinical Adhesion Bioassay

TNFSF/TNFRSF cytokine gene expression in sickle cell anemia: Up-regulated TNF-like cytokine 1A (TL1A) and its decoy receptor (DcR3) in peripheral blood mononuclear cells and plasma

Association of Circulating Transcriptomic Profiles with Mortality in Sickle Cell Disease

Biomarkers of clinical severity in treated and untreated sickle cell disease: a comparison by genotypes of a single center cohort and African Americans in the NHANES study.

In Vivoevaluation of the Effect of Sickle Cell Hemoglobin S, C and Therapeutic Transfusion on Erythrocyte Metabolism and Cardiorenal Dysfunction

In vivo evaluation of the effect of sickle cell hemoglobin S, C and therapeutic transfusion on erythrocyte metabolism and cardiorenal dysfunction