- Mehdi Nouraie1,
- Janet S. Lee2,3,
- Yingze Zhang2,
- Tamir Kanias2,3,
- Xuejun Zhao4,
- Zeyu Xiong2,
- Timothy B. Oriss3,
- Qilu Zeng2,
- Gregory J. Kato4,
- J. Simon R. Gibbs5,
- Mariana E. Hildesheim3,
- Vandana Sachdev4,
- Robyn J. Barst6,
- Roberto F. Machado7,
- Kathryn L. Hassell8,
- Jane A. Little9,
- Dean E. Schraufnagel7,
- Lakshmanan Krishnamurti10,
- Enrico Novelli2,
- Reda E. Girgis11,
- Claudia R. Morris12,
- Erika Berman Rosenzweig6,
- David B. Badesch8,
- Sophie Lanzkron11,
- Oswaldo L. Castro1,
- Jonathan C. Goldsmith13,
- Victor R. Gordeuk7,* and
- Mark T. Gladwin2,3,*⇓ Walk-PHASST Investigators and Patients
- 1Howard University, Washington, USA
- 2Vascular Medicine Institute, University of Pittsburgh, USA
- 3Division of Pulmonary, Allergy and Critical Care Medicine, University of Pittsburgh, USA
- 4Cardiovascular and Pulmonary Medicine Branch, NHLBI, Bethesda, MD, USA
- 5National Heart & Lung Institute, Imperial College London, UK
- 6Columbia University, New York, USA
- 7University of Illinois, Chicago, IL, USA
- 8University of Colorado HSC, Denver, CO, USA
- 9Case Western Reserve University, Cleveland, OH, USA
- 10Children’s Hospital of Pittsburgh, Pittsburgh, PA, USA
- 11Johns Hopkins University, Baltimore, MD, USA
- 12Children’s Hospital & Research Center Oakland, Oakland, CA, USA
- 13National Heart Lung and Blood Institute/NIH, Bethesda, MD, USA
- Correspondence: Mark T. Gladwin. E-mail:
The intensity of hemolytic anemia has been proposed as an independent risk factor for the development of certain clinical complications of sickle cell disease, such as pulmonary hypertension, hypoxemia and cutaneous leg ulceration. A composite variable derived from several individual markers of hemolysis could facilitate studies of the underlying mechanisms of hemolysis. In this study, we assessed the association of hemolysis with outcomes in sickle cell anemia. A hemolytic component was calculated by principal component analysis from reticulocyte count, serum lactate dehydrogenase, aspartate aminotransferase and total bilirubin concentrations in 415 hemoglobin SS patients. Association of this component with direct markers of hemolysis and clinical outcomes was assessed. As primary validation, both plasma red blood cell microparticles and cell-free hemoglobin concentration were higher in the highest hemolytic component quartile compared to the lowest quartile (P≤0.0001 for both analyses). The hemolytic component was lower with hydroxyurea therapy, higher hemoglobin F, and alpha-thalassemia (P≤0.0005); it was higher with higher systemic pulse pressure, lower oxygen saturation, and greater values for tricuspid regurgitation velocity, left ventricular diastolic dimension and left ventricular mass (all P<0.0001). Two-year follow-up analysis showed that a high hemolytic component was associated with an increased risk of death (hazard ratio, HR 3.44; 95% confidence interval, CI: 1.2–9.5; P=0.02). The hemolytic component reflects direct markers of intravascular hemolysis in patients with sickle cell disease and allows for adjusted analysis of associations between hemolytic severity and clinical outcomes. These results confirm associations between hemolytic rate and pulse pressure, oxygen saturation, increases in Doppler-estimated pulmonary systolic pressures and mortality (Clinicaltrials.gov identifier: NCT00492531).
Sickle cell anemia is a hemoglobinopathy that is characterized by hemolysis. After adjustment for known modulators of disease severity such as hemoglobin F and alpha thalassemia, the baseline rate of hemolysis within an individual with sickle cell anemia is stable over time but heterogeneous among individuals.1,2 We have used principal component analysis to derive a hemolytic component from reticulocyte count, lactate dehydrogenase, aspartate transaminase and total bilirubin in prior epidemiological studies of sickle cell disease.3 This standard statistical data reduction approach uses conventional clinical measurements to explain the maximum shared variance among these indirect measures of hemolysis.4 Principal component analysis is a long-established statistical method for studying underlying mechanisms reflected in individual biological or physical measurements4,5 and, therefore, useful when the objective of the study is to investigate the role of an underlying mechanism.5 Some authors have questioned the validity of the hemolytic component as a variable to investigate underlying mechanisms of hemolysis in sickle cell anemia based on a lack of validation with direct markers of hemolysis; they have also challenged the relationship between the intensity of hemolytic anemia and clinical manifestations of sickle cell disease, such as pulmonary and systemic hypertension, low oxygen saturation, and history of cutaneous leg ulceration.6,7 Direct effects of intravascular hemolysis include increased cell-free plasma hemoglobin and release of sub-micron red blood cell (RBC) microparticles which contain a substantial amount of hemoglobin.8 The walk-PHaSST study (treatment of pulmonary hypertension and sickle cell disease with sildenafil therapy) was designed to develop a large, multi-center, observational cohort of patients with sickle cell disease and to conduct a smaller interventional trial of treatment with sildenafil in the subset of patients with elevated tri-cuspid regurgitation velocity.9 In this report, we first sought to confirm whether the hemolytic component reflects direct indicators of hemolysis in the Walk-PHaSST cohort and then to determine the clinical manifestations of a heightened level of hemolysis as reflected in the component in this new and large clinical cohort.
Design and Methods
Subject selection and clinical evaluation
This report focuses on patients with hemoglobin SS aged 12 years and over who were studied at steady state at 9 United States Centers and one United Kingdom Center as previously described.10 Local institutional review boards or ethics committees approved the protocol and written informed consent was obtained. Hemoglobin oxygen saturation was measured by pulse oximetry at room temperature and atmospheric oxygen pressure. Echocardiography was performed at the participating institutions and read centrally in the NHLBI echocardiography core laboratory. Cardiac measurements were performed according to the American Society of Echocardiography guidelines.11 Percentages of hemoglobin S, F and A were measured by high-performance liquid chromatography (HPLC) (Ultra Resolution System, Trinity Biotech). The -α3.7 and -α4.2 thalassemia alleles were detected by molecular methodology based on polymerase chain reaction at the University of Pittsburgh as previously described.12 Serum N-terminal pro-brain natriuretic peptide (NT-pro BNP) concentration was measured by a sandwich immunoassay using polyclonal antibodies that recognize epitopes located in the N-terminal segment (1– 6) of pro-BNP (1–108) (Elecsys analyser; Roche Diagnostics, Mannheim, Germany).
Derivation of a hemolytic component
Principal component analysis5 was used to derive a hemolytic component from lactate dehydrogenase (site-adjusted values from linear regression), aspartate aminotransferase (site-adjusted values from linear regression), total bilirubin and reticulocyte percent in 415 hemoglobin SS patients using the Stata 10.1 software package. For primary validation, it was also derived in a subset of 235 of these 415 hemoglobin SS patients without recent blood transfusion as documented by the absence of hemoglobin A measured by HPLC analysis. Note that both of these populations represent new validation cohorts, designed to independently test prior observations from a large pediatric sickle cell patient population.3
Plasma cell-free hemoglobin and red blood cell microparticles
A standardized protocol was instituted for obtaining and processing heparinized venous blood samples from each of the coordinating centers within one hour of the blood draw. Venous blood samples were centrifuged at ~1400× g for 10 min at either 4°C or room temperature, depending upon the availability of a refrigerated centrifuge at each of the centers. Plasma samples were collected, stored at −80°C, and subsequently shipped on dry ice to a central biorepository at the University of Pittsburgh. There, the frozen plasma samples were thawed carefully on ice, individual aliquots prepared for plasma hemoglobin ELISA and RBC microparticle analysis, and the aliquots were immediately refrozen at −80°C until the day of experimentation. Preliminary studies performed on fresh and thawed frozen plasma samples from the same healthy volunteer showed a similar relationship in microparticle counts across varying conditions (Online Supplementary Figure S1).
From the subgroup of patients who did not have hemoglobin A in their HPLC analysis of hemoglobin, we selected those in the highest (n=57) and lowest (n=57) quartiles of the hemolytic component distribution for measurement of plasma cell-free hemoglobin concentration by an anti-hemoglobin ELISA13 and RBC microparticles by flow cytometry.8 We have recently detailed a method of RBC microparticle enumeration for clinical plasma samples utilizing a flow cytometric approach.14 Briefly, a master mixture of Flow Cytometry Absolute Counting Standard microbeads (7.6 micron size) (Bangs Laboratories Inc., Fishers, IN, USA) and MegamixTM fluorescent beads of three diameters (0.5, 0.9 and 3 microns) (BioCytex, Marseille, France) at a known volumetric ratio were spiked into each sample and vigorously vortexed immediately prior to flow cytometric analysis on a FACSAria flow cytometer (BD Biosciences, San Jose, CA, USA). We defined RBC microparticles based upon the presence of a cell-specific antigen (glycophorin A) and surface phosphatidylserine using Annexin V binding. Although it is recognized that glycophorin A positive (A+) but Annexin V negative (A–) events of less than 1 micron in size occur and likely represent RBC microparticles, the definition we used is in agreement with that set forth by the International Society on Thrombosis and Hemostasis Scientific and Standardization Subcommittees for platelet microparticles.15
Glycophorin A+ Annexin V+ events, as defined by glycophorin A-PE and Annexin V-FITC immunostaining above isotype fluorescence threshold, were further validated as RBC microparticles by their light scatter distribution relative to the 0.5 and 0.9 micron fluorescent beads spiked into each sample. RBC microparticle (mp) counts were calculated using the following equation: [number of 7.6 micron counting beads added to each sample X (number of mp events/number of 7.6 micron bead events)]. As 10 μL of each sample was assayed, the calculated RBC mp count was divided by 10 in order to express units as “RBC mp counts per μL”.
Follow up for mortality
Study subjects from the Walk-PHASST observational study from 10 participating centers were followed for a median of 2.4 years (IQR 2.1–2.8 years). Over this period, life status was determined for a total of 395 out of 415 hemoglobin SS patients aged 12 years and over, and 16 deaths were observed. The association of the hemolytic component at study entry with mortality was examined using Kaplan-Meier survival curves and Cox’s proportional hazards models comparing the upper tertile with the lower two tertiles of hemolytic component (66.7th percentile,1.28 relative units). The potentially confounding effects of age, gender, hydroxyurea use, hemoglobin, fetal hemoglobin, and white blood cell count on the association between hemolytic component and mortality were examined using Cox’s proportional hazards regression.
Linear regression analysis was used to adjust lactate dehydrogenase and aspartate aminotransferase for site. Continuous variables with a skewed distribution were converted to a normal distribution using natural log or square root. Student’s t-test was used to compare plasma cell-free hemoglobin and plasma RBC microparticles between the highest and lowest quartiles of hemolytic component in patients with no hemoglobin A. Step-wise multiple linear regression was used to assess the independent effect of hemolytic component on clinical and echocardiographic parameters. Pair-wise correlation between predictors and Variance Inflation Factor was used to assess co-linearity in each model. A potential interaction between site and predictors was evaluated in each model. Missing values were the result of lack of laboratory or echocardiographic observations in random subgroups of patients. Analyses were performed using Stata 10.1 software (StataCorp, College Station, TX, USA).
The hemolytic component derived from lactate dehydrogenase, aspartate aminotransferase, reticulocyte percentage and total bilirubin in 415 hemoglobin SS patients had a mean of 0 (SD=1.50) and predicted 55% of the variation among all four variables (Eigenvalue=2.20). Figure 1 demonstrates the relationship of the hemolytic component with the variables from which it was derived and the relationships among these variables.
The hemolytic component has strong relationships with direct and indirect markers of intravascular hemolysis
To determine whether the hemolytic component reflects direct indicators of intravascular hemolysis, we examined plasma concentrations of both RBC microparticles and cell free hemoglobin in the lowest and highest quartiles of hemolytic component. RBC microparticle counts (Figure 2; P=0.0001) and cell-free hemoglobin levels (Figure 3; P<0.0001) were significantly higher in the highest quartile of hemolytic component values compared to the lowest quartile. Out of 415 patients, 199 (48%) were receiving hydroxyurea treatment. The hemolytic component was lower with hydroxyurea therapy (Figure 4A; P=0.0002), α-thalassemia genotype (Figure 4B; P for trend=0.0005), and increasing percentage of hemoglobin F (Figure 4C; P<0.0001), supporting the hemolytic component’s relationship with the intensity of total hemolysis. The hemolytic component also correlated significantly and inversely with hemoglobin concentration (Figure 4D; P<0.0001). In multiple linear regression analysis, α-thalassemia single deletion (standardized beta = −0.15, P=0.002), double deletion (standardized beta = −0.23, P<0.0001), hemoglobin F (standardized beta = −0.40, P<0.0001), hemoglobin A (representing recent blood transfusion; standardized beta = −0.31, P<0.0001), female gender (standardized beta = −0.11, P=0.024), and hydroxyurea therapy (standardized beta = −0.11, P=0.038), were each independently associated with lower hemolytic component (Table 1).
Association of the hemolytic component with clinical findings among all sickle cell anemia patients
The hemolytic component was associated directly with systemic pulse pressure and inversely with hemoglobin oxygen saturation determined by pulse oximetry after accounting for multiple comparisons. The associations with oxygen saturation persisted after additionally adjusting for hemoglobin concentration and hydroxyurea therapy. The hemolytic component did not correlate with a history of leg ulcers among all patients, but it did so in the subgroup of patients not on hydroxyurea (P=0.001) (Tables 2 and 3).
Association of hemolytic component with BNP and echocardiographic findings among all sickle cell anemia patients
The frequency of elevated tricuspid regurgitation velocity (>3.0 m/s) was 34 (9.2%) and the frequency of extremely elevated velocity (>3.5 m/s) was 11 (3.0%). The hemolytic component was associated positively with the NT-proBNP concentration, tricuspid regurgitation velocity, left ventricular mass index, cardiac output, left ventricular end-diastolic volume, left atrial volume index and left ventricular diastolic dimension after accounting for multiple comparisons. Correlation with tricuspid regurgitation velocity, left ventricular mass index and left ventricular diastolic dimension persisted after adjusting for hemoglobin concentration and hydroxyurea therapy (Table 4)
Association of the hemolytic component with mortality
Using Kaplan-Meier survival analysis, a hemolytic component value in the upper tertile was associated with an increase in the risk of death over a median of 2.4 years of follow up (Figure 4; P=0.011). The risk of mortality associated with a hemolytic component value in the upper tertile was more than three times higher relative to a value in the lower two tertiles as determined using Cox’s proportional-hazards regression (HR 3.44; 95% CI: 1.2–9.5). Adjustment for age, gender, hydroxyurea use, hemoglobin, and white cell count had little effect on the magnitude or significance of the association between mortality and hemolytic component, while adjustment for fetal hemoglobin resulted in an elevated but not significant risk of death (HR 2.21; 95% CI: 0.7–6.7).
The results of this study support the concept that a statistical approach to identify shared variability among indirect hemolytic markers can produce a variable that reflects intravascular hemolytic rate. This provides a validated value for hemolytic rate that facilitates studies of the underlying mechanisms of hemolysis and that allows for adjustment of confounding variables, including hemoglobin level, and analysis of associated clinical phenotypes. While this study focuses on sickle cell anemia, it is likely that the hemolytic component could be applied to the clinical study of other hemolytic disorders, including extravascular hemolytic disorders.
There is considerable data to suggest that the rate of hemolysis and the degree of anemia have both overlapping, and divergent pathophysiological and clinical consequences in sickle cell anemia. A decreased hemoglobin concentration leads to decreased delivery of oxygen to the tissues. Intravascular hemolysis leads to the release of hemoglobin and arginase-1 into the plasma, the scavenging of plasma NO by cell-free hemoglobin, and the depletion by arginase-1 of plasma arginine, the obligate substrate for the NO synthases.16,17 These processes may contribute to reduced NO bioavailability and vascular dysfunction. Recent studies suggest that hemolysis also increases plasma levels of hemin, which may activate the innate immune system and drive inflammatory respons-es.18,19 Hemolysis also drives platelet and hemostatic activation and is proposed to generate reactive oxygen species and activate vascular oxidases.20–23 A progressive vasculopathy related to hemolytic severity develops in some patients with sickle cell disease24,25 and is characterized by systemic hypertension, Doppler-echocardiography estimated elevation in systolic pulmonary artery pressure, endothelial dysfunction, intimal and smooth muscle proliferative changes in conduit blood vessels, and increased risk of death.26,27 While this hypothesis for phenotypic manifestations determined in part by the process of intravascular hemolysis has been challenged,6,7,28 our study provides substantial support for this model. The calculation of a hemolytic component resolves the problem of dealing with correlated predictors in multivariate analyses, allows for adjustment of important potential confounders, and allows for adjustment in multi-center studies for site variability in laboratory assay protocols and standards.
The current validated analysis of hemolytic rates in this large cohort of patients with homozygous hemoglobin-S disease suggests that hemolytic anemia in general, and intravascular hemolysis specifically, may constitute an independent risk factor for the development or manifestation of certain clinical complications or endophenotypes (subphenotypes).29,30 This study also confirms prior reported associations using single indirect biomarkers of hemolysis such as lactate dehydrogenase, including high pulse pressure, Doppler-estimated increase in systolic pulmonary artery pressure, low hemoglobin-oxygen saturation and prospective risk of death.27,31 The mechanism for the association between high degree of hemolysis and lower oxygen saturation is not known, but it is hypothesized to represent dysregulation of blood flow to the lung, causing imbalance in ventilation and perfusion matching.30,32 We were unable to confirm prior observations on the association between self-reported history of priapism and hemolytic anemia. In addition, there were not sufficient stroke events in our study to evaluate this complication. These studies also confirm a lack of association with self-reported history of rates of vaso-occlusive pain crisis and acute chest syndrome, which have been prospectively shown to correlate with higher steady state hemoglobin values in prior registry studies of sickle cell disease patients.33,34 Interestingly, this paradox was evident in the Co-operative Study of Sickle Cell Anemia, in which a high steady state hemoglobin value was associated with increased risk of vaso-occlusive crisis and the acute chest syndrome, but low hemoglobin was associated with increased risk of stroke and death.33–36
The use of a validated hemolytic component may be useful in future clinical trials to stratify patients and direct personalized treatment approaches to disease prognostication and therapeutic intervention. For example, a drug that inhibits hemolysis directly, such as the Gardos channel inhibitor senicapoc, may not be effective in reducing vaso-occlusive events and, in fact, increased vaso-occlusive events significantly in patients not taking hydroxyurea.37 This drug may have been effective in patients with more severe hemolytic anemia, as characterized by a high quartile hemolytic component and, if taken with hydroxyurea to limit vaso-occlusive complications, could conceivably improve hemoglobin oxygen saturations and reduce the risk of pulmonary and systemic hypertension, leg ulcers, and death.38,39 We included patients on hydroxyurea in order to have a representative sample of sickle cell anemia patients and to ensure that our findings have practical application. Furthermore, a validated hemolytic component may be useful to explore as a potential biomarker to assess efficacy in future clinical therapeutic trials. RBC microparticle numbers are increased in both steady-state and sickle cell crisis,40 inversely correlate with Hb levels in crisis and steady-state,40 and positively correlate with plasma free Hb levels in steady state.41 We have recently validated our methods for enumerating RBC microparticles in clinical plasma specimens.25 As detailed in the Design and Methods section, a standardized protocol was instituted for obtaining and processing heparinized venous blood samples from each of the co-ordinating centers within one hour of the blood draw. All samples were processed and handled in the same manner, and batched for flow cyto-metric analysis in a blinded fashion thus making relative quantitation of RBC microparticles in stored plasma specimens feasible. We do not have the data to investigate any role of auto-splenectomy on hemolytic variables and number of circulating red cell microparticles in this study, but this should be evaluated in future studies with defined splenic status or before and after splenectomy. It has been reported that splenectomy in thalassemia intermedia patients is associated with higher plasma hemoglobin and red cell microparticle levels.42
There are a number of limitations to this study. Episodes of painful vaso-occlusive crisis, acute chest syndrome and priapism were self-reported and, therefore, subject to recall bias. The laboratory measurements used were not centrally analyzed, hence the need for ‘site-adjusted’ values. The total bilirubin concentration rather than the indirect bilirubin concentration was used to calculate the hemolytic component because direct bilirubin concentrations were not measured. Patients receiving treatment with hydroxyurea were included in this study, and hydroxyurea is likely to have other effects on sickle cell complications than just reducing hemolysis. The pulse ox test was performed at all participating centers without standardization and, therefore, may have limited accuracy in monitoring oxygen saturation43 in certain circumstances such as the presence of methemoglobinemia. However, a number of studies have indicated that pulse oximetry is reliable for diagnosing hypoxia during steady state conditions and with acute complications in sickle cell disease patients.44-46 Details of deaths were not available from Walk-PHASST because of lack of funds for follow up beyond alive or dead status. There is a possibility that some of the fatal events were not related to sickle cell disease and this could not be accounted for by the survival analysis.
In summary, we found that the hemolytic component correlated with direct markers of intravascular hemolysis in patients with sickle cell anemia, namely plasma red blood cell microparticles and cell-free hemoglobin concentration. It also correlated with measures known to influence total hemolytic rate, namely hemoglobin F percent and the presence of alpha-thalassemia. Thus, this variable produced by a data reduction statistical methodology can potentially help in studies of the underlying mechanisms of hemolysis. Even after adjustment for hemoglobin concentration and hydroxyurea therapy, higher hemolytic component was associated with higher systemic pulse pressure, lower oxygen saturation, greater tricuspid regurgitation velocity and left ventricular diastolic dimension and mass, and increased risk of death. Thus, the hemolytic component reflects direct markers of intravascular hemolysis in patients with sickle cell disease, correlates with important clinical outcomes, and allows for adjusted analysis of associations between hemolytic severity and clinical outcomes.
↵* Dr. Mark T. Gladwin and Victor R. Gordeuk share senior authorship.
Authorship and Disclosures
Information on authorship, contributions, and financial & other disclosures was provided by the authors and is available with the online version of this article at www.haematologica.org.
This work was supported by funds from the National Heart, Lung, and Blood Institute, National Institutes of Health, Department of Health and Human Services, under contract HHSN268200617182C. This work was also supported in part by NIH CTSA grant UL1 RR024131 (to CRM), grant ns. 2 R25 HL003679–08 and 1 R01 HL079912–02 from NHLBI (to VRG), by Howard University GCRC grant n. 2MOI RR10284– 10 from NCRR, NIH, Bethesda, MD, by R01HL086884 and R01HL086884 03S1 (to JSL), by NIH grants R01HL098032 and RO1HL096973 (to MTG), by the Institute for Transfusion Medicine and the Hemophilia Center of Western Pennsylvania and by the intramural research program of the National Institutes of Health.
- Received May 2, 2012.
- Accepted August 14, 2012.
- Copyright© Ferrata Storti Foundation