Hypoxic response contributes to altered gene expression and precapillary pulmonary hypertension in patients with sickle cell disease.

BACKGROUND: We postulated that the hypoxic response in sickle cell disease (SCD) contributes to altered gene expression and pulmonary hypertension, a complication associated with early mortality. METHODS AND RESULTS: To identify genes regulated by the hypoxic response and not other effects of chronic anemia, we compared expression variation in peripheral blood mononuclear cells from 13 subjects with SCD with hemoglobin SS genotype and 15 subjects with Chuvash polycythemia (VHL(R200W) homozygotes with constitutive upregulation of hypoxia-inducible factors in the absence of anemia or hypoxia). At a 5% false discovery rate, 1040 genes exhibited >1.15-fold change in both conditions; 297 were upregulated and 743 downregulated including MAPK8 encoding a mitogen-activated protein kinase important for apoptosis, T-cell differentiation, and inflammatory responses. Association mapping with a focus on local regulatory polymorphisms in 61 patients with SCD identified expression quantitative trait loci for 103 of these hypoxia response genes. In a University of Illinois SCD cohort, the A allele of a MAPK8 expression quantitative trait locus, rs10857560, was associated with precapillary pulmonary hypertension defined as mean pulmonary artery pressure ≥25 mm Hg and pulmonary capillary wedge pressure ≤15 mm Hg at right heart catheterization (allele frequency, 0.66; odds ratio, 13.8; n=238). This association was confirmed in an independent Walk-Treatment of Pulmonary Hypertension and Sickle Cell Disease With Sildenafil Therapy cohort (allele frequency, 0.65; odds ratio, 11.3; n=519). The homozygous AA genotype of rs10857560 was associated with decreased MAPK8 expression and present in all 14 of the identified precapillary pulmonary hypertension cases among the combined 757 patients. CONCLUSIONS: Our study demonstrates a prominent hypoxic transcription component in SCD and a MAPK8 expression quantitative trait locus associated with precapillary pulmonary hypertension.

Iron deficiency modifies gene expression variation induced by augmented hypoxia sensing.

In congenital Chuvash polycythemia (CP), VHL(R200W) homozygosity leads to elevated hypoxia inducible factor (HIF) levels at normoxia. CP is often treated by phlebotomy resulting in iron deficiency, permitting us to examine the separate and synergistic effects of iron deficiency and HIF signaling on gene expression. We compared peripheral blood mononuclear cell gene expression profiles of eight VHL(R200W) homozygotes with 17 wildtype individuals with normal iron status and found 812 up-regulated and 2120 down-regulated genes at false discovery rate of 0.05. Among differential genes we identified three major gene regulation modules involving induction of innate immune responses, alteration of carbohydrate and lipid metabolism, and down-regulation of cell proliferation, stress-induced apoptosis and T-cell activation. These observations suggest molecular mechanisms for previous observations in CP of lower blood sugar without increased insulin and low oncogenic potential. Studies including 16 additional VHL(R200W) homozygotes with low ferritin indicated that iron deficiency enhanced the induction effect of VHL(R200W) for 50 genes including hemoglobin synthesis loci but suppressed the effect for 107 genes enriched for HIF-2 targets. This pattern is consistent with potentiation of HIF-1α protein stability by iron deficiency but a trend for down-regulation of HIF-2α translation by iron deficiency overriding an increase in HIF-2α protein stability.

Gene expression changes in sickle cell reticulocytes and their clinical associations.

Transcriptional changes in compensatory erythropoiesis in sickle cell anemia (SCA) and their disease modulation are unclear. We detected 1226 differentially expressed genes in hemoglobin SS reticulocytes compared to non-anemic hemoglobin AA controls. Assessing developmental expression changes in hemoglobin AA erythroblasts for these genes suggests heightened terminal differentiation in early erythroblasts in SCA that diminishes toward the polychromatic to orthochromatic stage transition. Comparison of reticulocyte gene expression changes in SCA with that in Chuvash erythrocytosis, a non-anemic disorder of increased erythropoiesis due to constitutive activation of hypoxia inducible factors, identified 453 SCA-specific changes attributable to compensatory erythropoiesis. Peripheral blood mononuclear cells (PBMCs) in SCA contain elevated proportions of erythroid progenitors due to heightened erythropoiesis. Deconvolution analysis in PBMCs from 131 SCA patients detected 54 genes whose erythroid expression correlated with erythropoiesis efficiency, which were enriched with SCA-specific changes (OR = 2.9, P = 0.00063) and annotation keyword "ubiquitin-dependent protein catabolic process", "protein ubiquitination", and "protein polyubiquitination" (OR = 4.2, P = 7.5 × 10-5). An erythroid expression quantitative trait locus of one of these genes, LNX2 encoding an E3 ubiquitin ligase, associated with severe pain episodes in 774 SCA patients (OR = 1.7, P = 3.9 × 10-5). Thus, erythroid gene transcription responds to unique conditions within SCA erythroblasts and these changes potentially correspond to vaso-occlusive manifestations.