Hydroxyurea therapy in children with SCA decreases plasma markers of endothelial activation, notably angiopoietin-2.
Angiopoietin-2 levels are associated with TCD velocities, a surrogate marker for primary stroke risk in SCA.
Visual Abstract
Hydroxyurea reduces morbidity and mortality in children with sickle cell anemia (SCA). The endothelium is central to SCA-related complications including stroke. However, hydroxyurea’s impact on the endothelium is not well described. To address this gap, we measured plasma levels of endothelial activation markers (angiopoietin-2, P-selectin, soluble endothelial selectin [sE-selectin], soluble intercellular cellular adhesion molecule 1, and soluble vascular endothelial cellular adhesion molecule) by enzyme-linked immunosorbent assay after initiation of hydroxyurea therapy. Samples were collected from Ugandan children with SCA enrolled in a clinical trial evaluating hydroxyurea vs placebo (NOHARM trial). Samples were collected at enrollment; and then after 2, 4, and 12 months of follow-up. Longitudinal changes in biomarker levels were evaluated using linear mixed effects models. Transcranial Doppler (TCD) velocities were measured at 10 to 12 months follow-up to assess cerebral blood flow and primary stroke risk. Mediation analysis was used to explore causal pathways of hydroxyurea-mediated effects on TCD velocities. In total, 798 plasma samples were tested from 205 children (mean enrollment age, 2.2 years). At enrollment, higher levels of angiopoietin-2 were associated with a previous medical history of dactylitis, vaso-occlusive crises, acute chest syndrome, and transfusion (P < .05 for all). Hydroxyurea therapy at a fixed dose of 20 mg/kg per day decreased plasma angiopoietin-2, P-selectin, and sE-selectin levels over the study period (P < .05 for all). Angiopoietin-2 and sE-selectin were associated with higher TCD velocities. Mediation analysis suggests that hydroxyurea decreases TCD velocities through an increase in fetal and total hemoglobin. Increased fetal and total hemoglobin, and decreased white blood cell count may decrease TCD velocity, in part, through an angiopoiten-2–mediated pathway. This trial was registered at www.ClinicalTrials.gov as #NCT01976416.
Introduction
Sickle cell anemia (SCA) is a hemoglobinopathy caused by the homozygous recessive βS mutation in the β-globin gene resulting in hemoglobin S (HbS) expression. Most births affected by SCA occur in sub-Saharan Africa, where mortality remains high.1-3 In SCA, the polymerization and precipitation of deoxygenated HbS results in dense sickled red blood cells, increased red blood cell rigidity, damage to the cell membrane, and hemolysis. SCA leads to a chronic inflammatory state associated with vascular dysfunction. This interplay of red blood cell sickling, hemolysis, and chronic inflammation at the level of the endothelium results in progressive vaso-occlusive end-organ damage.4,5 Significant sequelae include cerebrovascular complications involving large cerebral arterial vessels in patients with SCA. SCA is a leading cause of stroke in children, with a life-time risk of stroke between 25% and 30% without disease-modifying therapy.6
Endothelial activation can be assessed by measuring circulating levels of soluble cellular adhesion molecules (CAMs).7 CAMs regulate leukocyte trafficking and undergo proteolytic cleavage to generate soluble CAMs, including soluble intercellular CAM 1 (sICAM-1), soluble vascular endothelial CAM (sVCAM-1), and soluble endothelial selectin (sE-selectin).4,8,9 Endothelial activation also results in the rapid mobilization and release of Weibel-Palade bodies containing P-selectin and angiopoietin-2 (Angpt-2).10 Angpt-2 is a vascular growth factor that is released in the context of endothelial activation and can lead to increased vascular responsiveness to inflammatory cytokines through the upregulation of CAMs, and can lead to vascular leak and endothelial cell apoptosis when produced in excess.11
Hydroxyurea is a disease-modifying therapy that functions by increasing fetal Hb (HbF), reducing red blood cell polymerization and hemolysis, and reducing white blood cell (WBC) and platelet counts, at least partially, through myelosuppression.12 Hydroxyurea therapy is associated with fewer sickle-related complications including strokes, and reduces all-cause mortality in African children.13,14 We hypothesized that hydroxyurea may also improve clinical outcomes through improved microvascular function. We tested this hypothesis in a randomized, double-blind, placebo-controlled clinical trial with Ugandan children with SCA receiving fixed-dose hydroxyurea or placebo for 1 year (NOHARM).13 An a priori defined secondary end point was to evaluate changes in markers of endothelial activation and inflammation, and transcranial Doppler (TCD) velocities associated with hydroxyurea therapy. Using serial samples collected over the first year of the study, we evaluated whether markers of endothelial activation and inflammation changed with hydroxyurea treatment, and used mediation analysis to assess the extent to which hydroxyurea improved TCD velocities as a surrogate measure of stroke risk through changes in endothelial activation as an intermediate, or mediating, pathway.15
Methods
Study participants
NOHARM was a randomized, double-blinded phase 3 placebo-controlled clinical trial of oral hydroxyurea or placebo at a fixed dose of 20 ± 2.5 mg/kg of body weight per day, as previously described (ClinicalTrials.gov identifier: NCT01976416).13,16 Inclusion criteria were documented SCA (HbSS by Hb electrophoresis); age of between 1.00 and 3.99 years at enrollment; weight of at least 5.0 kg at enrollment; residence within the area of Kampala, Uganda, where malaria is endemic; attending the Mulago Hospital Sickle Cell Clinic; and a willingness to comply with all study-related treatments, evaluations, and follow-up. Exclusion criteria included active use of hydroxyurea; known history of HIV, malignancy, or active clinical tuberculosis; severe malnutrition, defined as a weight-for-length z score of less than −3 using the World Health Organization growth standards; preexisting hematologic toxicity; a transfusion within 30 days; or alanine aminotransferase or creatinine level of >2× the upper limit of normal for age.
All children enrolled in the study received routine care for SCA according to the Mulago Hospital Sickle Cell Clinic guidelines including daily folic acid, penicillin prophylaxis daily for children aged <5 years, sulfadoxine-pyrimethamine prophylaxis for malaria prevention, and antihelminth medication according to the Ministry of Health guidelines. All children received an insecticide-treated bed net to prevent malaria.
Study procedures
On enrollment, children received a standard physical examination and had a detailed history recorded to evaluate sickle-related complications within the previous year. A blood sample was collected for confirmation of sickle cell status by Hb electrophoresis, a complete blood count and reticulocyte count, and biochemistries to assess liver and kidney function. EDTA–anticoagulated whole blood was centrifuged at 1200g for 20 minutes, and then the plasma layer was carefully collected using a sterile transfer pipette, and aliquots were stored in barcoded microtubes at −80°C until testing. Plasma samples were stored for biomarker evaluation at enrollment, and at 2-, 4- and 12-months follow-up. Between 10 and 12 months of follow-up, children had cerebral arterial blood flow velocities assessed by a single, certified examiner, with results confirmed off-site by central review, as previously described.17,18
Study approval
The study was approved in Uganda by the Makerere School of Medicine research ethics committee, the Mulago Hospital research ethics committee, the Uganda National Drug Authority, and the Uganda National Council of Science and Technology, and in the United States by the institutional review boards at Indiana University and Cincinnati Children's Hospital Medical Center. Caregivers of the children provided written informed consent for participation in the trial with the use of forms in English or the local language.
Immunoassays
A secondary outcome of the clinical trial included a change in markers of endothelial activation with hydroxyurea treatment.16 Enzyme-linked immunosorbent assays were conducted in batch on plasma samples using DuoSet enzyme-linked immunosorbent assays by R&D Systems (Minneapolis, MN) using methods previously described.19 We measured levels of Angpt-2 (1:5 dilution; assay range, 0.24-60 ng/mL; coefficient of variation (CV), 4.97%), P-selectin (1:50 dilution; assay range, 3.13-800 ng/mL; CV, 5.17%), C-reactive protein (1:10 000 to 1:20 000 dilution; assay range, 0.16-80 mg/dL; CV, 9.87%), soluble VCAM-1 (sVCAM-1, 1:5000 dilution; assay range, 76.8-20 000 ng/mL; CV, 6.38%), sICAM-1 (1:1000 dilution; assay range, 29.8-4000 ng/mL; CV, 7.47%), and sE-selectin (1:50 dilution; assay range, 2.35-600 ng/mL; CV, 3.87%). All sample testing was conducted in duplicate by technicians blinded to participant details.
Statistical analysis
Analysis was conducted using STATA version 17 (StataCorp LP, College Station, TX) and GraphPad Prism version 7.03 (GraphPad Software Inc, San Diego, CA). Clinical and laboratory data were described using the mean and standard deviation, or number (percentage). To evaluate the relationship between log10-transformed biomarker levels and a history of sickle-related complications, logistic regression was used. Differences in longitudinal biomarker levels based on study treatment were evaluated using linear mixed effects models in which observations within participants were correlated using a patient-specific intercept, and time points were treated as categorical variables and adjusted for participant sex and age at enrollment. A banded diagonal covariance matrix was used to model the within-subject variance–covariance errors, and the mixed model was fit by restricted maximum likelihood, and the Kenward-Roger approximations was used to estimate the denominator degrees of freedom. Spearman correlation was used to evaluate the relationship between biomarker levels and complete blood count measures. The Benjamini-Hochberg false discovery rate was applied at a threshold of 0.05 to adjust for multiple comparisons, as indicated.
We leveraged the serial collection of data to explore multiple physiologically possible pathways through which endothelial activation may affect TCD velocities through the inclusion of intermediate, or mediating, variables. Mediation models aim to elucidate the underlying mechanisms or pathways through which the independent variable affects the outcome. The direct effect describes the relationship between the independent variable and study outcome without considering the mediating variables, whereas the indirect effect captures the influence of the independent variable on the study outcome via its effect on the mediating variables, which then influence the study outcome. Together, the direct and indirect effects encompass the total effect of the independent variable on the outcome of interest. Models were constructed using measures of HbF at 2 months, total Hb and WBC count at 3 months, angiopoietin-2 at 4 months, and TCD velocity at 10 to 12 months. We started by examining the temporal relationships between (1) hydroxyurea and HbF at 2 months; (2) Hb or WBC count at 3 months, and angiopoietin-2 at 4 months; and (3) angiopoietin-2 at 4 months, and TCD velocities at 10 to 12 months using linear regression models. Following the method of VanderWeele and Vansteelandt for estimating mediation effects for multiple mediators,20 we estimated the indirect effects of hydroxyurea on TCD velocities through Hb and through WBC count. We used the Stata Structural Equation Modeling Builder for mediation analysis, and all analyses were adjusted for age and sex.21
Results
Plasma samples were available for 205 children in the NOHARM trial at enrollment and comprised 798 samples, 400 from children treated with hydroxyurea and 398 from children randomized to receive placebo. The mean (± standard deviation) age of children at enrollment was 2.2 years (±0.87), and 46.3% of children were female (Table 1). A history of sickle-cell related complications was common, with 57.1% of children hospitalized in the year before enrollment; 54.2% with a history of blood transfusion, 78.1% reporting a history of dactylitis, 81.5% reporting vaso-occlusive crises (VOC), and 15.6% reporting a history of acute chest syndrome. None of the children had a history of clinical stroke on enrollment.
Description of children included in the endothelial analysis
. | Placebo (n = 102) . | Hydroxyurea (n = 103) . | P value . |
---|---|---|---|
Demographics | |||
Age at enrollment, mean (±SD), y | 2.3 (0.85) | 2.2 (0.88) | .383 |
Female, n (%) | 46 (44.7) | 49 (47.6) | .722 |
Height, cm | 86.2 (8.2) | 85.6 (8.9) | .616 |
Weight, kg | 11.5 (2.1) | 11.3 (2.1) | .453 |
Weight-for-height z score | −0.41 (1.1) | −0.40 (1.1) | .959 |
Past medical history at enrollment, n (%) | |||
Dactylitis | 81 (79.4) | 79 (76.7) | .639 |
VOC | 80 (78.4) | 87 (84.5) | .266 |
Stroke | 0 (0.0) | 0 (0.0) | -- |
Splenomegaly | 5 (4.9) | 6 (5.8) | .769 |
Acute chest syndrome | 12 (11.8) | 20 (19.4) | .131 |
Transfusion | 56 (54.9) | 55 (53.4) | .829 |
Hospitalization in previous 1 y | 53 (52.0) | 64 (62.1) | .141 |
Laboratory measures | |||
Enrollment | |||
Hb, g/dL | 7.5 (0.96) | 7.5 (1.1) | .562 |
MCV, fL | 79.6 (9.4) | 79.4 (9.1) | .885 |
HbF/(HbF + HbS) | 13.3 (6.0) | 14.5 (7.2) | .331 |
ARC, × 109/μL | 0.38 (0.11) | 0.38 (0.12) | .687 |
WBC, × 109/μL | 18.7 (5.4) | 19.0 (7.2) | .906 |
ANC, × 109/μL | 6.2 (2.7) | 6.5 (3.1) | .375 |
Platelets, × 109/μL | 419 (137) | 357 (172) | .001∗ |
ALT, U/L | 18.8 (9.2) | 18.1 (8.8) | .451 |
Creatinine, μmol/L | 24.6 (6.0) | 25.2 (8.4) | .619 |
12-month follow-up | |||
Hb, g/dL | 7.4 (0.92) | 8.7 (1.3) | <.001∗ |
MCV, fL | 80.5 (8.1) | 88.1 (8.7) | <.001∗ |
HbF/(HbF + HbS) | 10.4 (4.9) | 22.4 (9.2) | <.001∗ |
ARC, × 109/μL | 0.39 (0.12) | 0.24 (0.11) | <.001∗ |
WBCs, × 109/μL | 18.0 (5.2) | 13.7 (5.2) | <.001∗ |
ANC, × 109/μL | 6.6 (2.6) | 5.1 (2.5) | <.001∗ |
Platelets, × 109/μL | 446 (144) | 371 (167) | .001∗ |
ALT, U/L | 17.3 (6.7) | 18.6 (8.3) | .448 |
Creatinine, μmol/L | 0.33 (0.09) | 0.31 (0.10) | .189 |
. | Placebo (n = 102) . | Hydroxyurea (n = 103) . | P value . |
---|---|---|---|
Demographics | |||
Age at enrollment, mean (±SD), y | 2.3 (0.85) | 2.2 (0.88) | .383 |
Female, n (%) | 46 (44.7) | 49 (47.6) | .722 |
Height, cm | 86.2 (8.2) | 85.6 (8.9) | .616 |
Weight, kg | 11.5 (2.1) | 11.3 (2.1) | .453 |
Weight-for-height z score | −0.41 (1.1) | −0.40 (1.1) | .959 |
Past medical history at enrollment, n (%) | |||
Dactylitis | 81 (79.4) | 79 (76.7) | .639 |
VOC | 80 (78.4) | 87 (84.5) | .266 |
Stroke | 0 (0.0) | 0 (0.0) | -- |
Splenomegaly | 5 (4.9) | 6 (5.8) | .769 |
Acute chest syndrome | 12 (11.8) | 20 (19.4) | .131 |
Transfusion | 56 (54.9) | 55 (53.4) | .829 |
Hospitalization in previous 1 y | 53 (52.0) | 64 (62.1) | .141 |
Laboratory measures | |||
Enrollment | |||
Hb, g/dL | 7.5 (0.96) | 7.5 (1.1) | .562 |
MCV, fL | 79.6 (9.4) | 79.4 (9.1) | .885 |
HbF/(HbF + HbS) | 13.3 (6.0) | 14.5 (7.2) | .331 |
ARC, × 109/μL | 0.38 (0.11) | 0.38 (0.12) | .687 |
WBC, × 109/μL | 18.7 (5.4) | 19.0 (7.2) | .906 |
ANC, × 109/μL | 6.2 (2.7) | 6.5 (3.1) | .375 |
Platelets, × 109/μL | 419 (137) | 357 (172) | .001∗ |
ALT, U/L | 18.8 (9.2) | 18.1 (8.8) | .451 |
Creatinine, μmol/L | 24.6 (6.0) | 25.2 (8.4) | .619 |
12-month follow-up | |||
Hb, g/dL | 7.4 (0.92) | 8.7 (1.3) | <.001∗ |
MCV, fL | 80.5 (8.1) | 88.1 (8.7) | <.001∗ |
HbF/(HbF + HbS) | 10.4 (4.9) | 22.4 (9.2) | <.001∗ |
ARC, × 109/μL | 0.39 (0.12) | 0.24 (0.11) | <.001∗ |
WBCs, × 109/μL | 18.0 (5.2) | 13.7 (5.2) | <.001∗ |
ANC, × 109/μL | 6.6 (2.6) | 5.1 (2.5) | <.001∗ |
Platelets, × 109/μL | 446 (144) | 371 (167) | .001∗ |
ALT, U/L | 17.3 (6.7) | 18.6 (8.3) | .448 |
Creatinine, μmol/L | 0.33 (0.09) | 0.31 (0.10) | .189 |
Data presented as mean (±SD) unless otherwise indicated. P values obtained from Student t test for Gaussian continuous variables, Wilcoxon rank-sum test for non-Gaussian continuous variables, and χ2 test for categorical variables.
ANC, automated neutrophil count; ALT, alanine aminotransferase; ARC, automated reticulocyte count; MCV, mean corpuscular volume; SD, standard deviation.
Significant after adjustment for 30 comparisons, as described in “Methods.”
Plasma levels of endothelial activation markers were assessed in enrollment samples, and we evaluated whether a history of sickle-related complications was associated with differences in endothelial activation using regression. A history of dactylitis, VOC, acute chest syndrome, and transfusion was associated with increased Angpt-2 levels at enrollment after adjustment for age and sex (all P < .05 before correction). After correction for multiplicity, only a history of dactylitis and transfusion remained significantly associated with Angpt-2 levels at enrollment (Table 2). A history of dactylitis was also associated with increased sE-selectin on admission although the relationship was not significant after adjusting for multiple comparisons.
Plasma markers of endothelial activation and inflammation at enrollment and a history of SCA-related complications
. | n/N (%) . | Adjusted β (95% CI) . | P value . |
---|---|---|---|
Angpt-2 | |||
Dactylitis | 160/205 (78.1) | 0.30 (0.14-0.47) | <.001∗ |
VOC | 167/205 (81.5) | 0.22 (0.03-0.40) | .022 |
ACS | 32/205 (15.6) | 0.21 (0.02-0.41) | .034 |
Transfusion | 111/205 (54.2) | 0.23 (0.09-0.37) | .002∗ |
Hb | |||
Dactylitis | 160/205 (78.1) | −0.05 (−0.10 to −0.01) | .023 |
VOC | 167/205 (81.5) | 0.01 (−0.04 to 0.05) | .789 |
ACS | 32/205 (15.6) | −0.01 (−0.06 to 0.05) | .823 |
Transfusion | 111/205 (54.2) | −0.07 (−0.10 to −0.03) | <.001∗ |
HbF | |||
Dactylitis | 160/205 (78.1) | −0.12 (−0.29 to 0.06) | .189 |
VOC | 167/205 (81.5) | −0.16 (−0.35 to 0.02) | .083 |
ACS | 32/205 (15.6) | 0.06 (−0.13 to 0.26) | .524 |
Transfusion | 111/205 (54.2) | −0.24 (−0.38 to −0.10) | .001∗ |
P-selectin | |||
Dactylitis | 160/205 (78.1) | 0.01 (−0.11 to 0.14) | .822 |
VOC | 167/205 (81.5) | 0.002 (−0.14 to 0.14) | .981 |
ACS | 32/205 (15.6) | 0.06 (−0.09 to 0.20) | .456 |
Transfusion | 111/205 (54.2) | −0.004 (−0.11 to 0.10) | .947 |
CRP | |||
Dactylitis | 160/205 (78.1) | 0.001 (−0.46 to 0.46) | .998 |
VOC | 167/205 (81.5) | −0.11 (−0.61 to 0.39) | .665 |
ACS | 32/205 (15.6) | −0.47 (−1.00 to 0.05) | .076 |
Transfusion | 111/205 (54.2) | 0.17 (−0.21 to 0.56) | .378 |
sVCAM-1 | |||
Dactylitis | 160/205 (78.1) | 0.02 (−0.15 to 0.18) | .845 |
VOC | 167/205 (81.5) | 0.02 (−0.16 to 0.20) | .853 |
ACS | 32/205 (15.6) | 0.02 (−0.17 to 0.21) | .819 |
Transfusion | 111/205 (54.2) | −0.04 (−0.17 to 0.10) | .607 |
sICAM-1 | |||
Dactylitis | 160/205 (78.1) | −0.001 (−0.13 to 0.13) | .990 |
VOC | 167/205 (81.5) | −0.10 (−0.24 to 0.04) | .167 |
ACS | 32/205 (15.6) | −0.004 (−0.15 to 0.14) | .955 |
Transfusion | 111/205 (54.2) | −0.05 (−0.15 to 0.06) | .397 |
sE-selectin | |||
Dactylitis | 160/205 (78.1) | 0.18 (0.04-0.32) | .014 |
VOC | 167/205 (81.5) | −0.01 (−0.17 to 0.14) | .859 |
ACS | 32/205 (15.6) | −0.07 (−0.24 to 0.10) | .407 |
Transfusion | 111/205 (54.2) | 0.14 (0.03-0.26) | .018 |
. | n/N (%) . | Adjusted β (95% CI) . | P value . |
---|---|---|---|
Angpt-2 | |||
Dactylitis | 160/205 (78.1) | 0.30 (0.14-0.47) | <.001∗ |
VOC | 167/205 (81.5) | 0.22 (0.03-0.40) | .022 |
ACS | 32/205 (15.6) | 0.21 (0.02-0.41) | .034 |
Transfusion | 111/205 (54.2) | 0.23 (0.09-0.37) | .002∗ |
Hb | |||
Dactylitis | 160/205 (78.1) | −0.05 (−0.10 to −0.01) | .023 |
VOC | 167/205 (81.5) | 0.01 (−0.04 to 0.05) | .789 |
ACS | 32/205 (15.6) | −0.01 (−0.06 to 0.05) | .823 |
Transfusion | 111/205 (54.2) | −0.07 (−0.10 to −0.03) | <.001∗ |
HbF | |||
Dactylitis | 160/205 (78.1) | −0.12 (−0.29 to 0.06) | .189 |
VOC | 167/205 (81.5) | −0.16 (−0.35 to 0.02) | .083 |
ACS | 32/205 (15.6) | 0.06 (−0.13 to 0.26) | .524 |
Transfusion | 111/205 (54.2) | −0.24 (−0.38 to −0.10) | .001∗ |
P-selectin | |||
Dactylitis | 160/205 (78.1) | 0.01 (−0.11 to 0.14) | .822 |
VOC | 167/205 (81.5) | 0.002 (−0.14 to 0.14) | .981 |
ACS | 32/205 (15.6) | 0.06 (−0.09 to 0.20) | .456 |
Transfusion | 111/205 (54.2) | −0.004 (−0.11 to 0.10) | .947 |
CRP | |||
Dactylitis | 160/205 (78.1) | 0.001 (−0.46 to 0.46) | .998 |
VOC | 167/205 (81.5) | −0.11 (−0.61 to 0.39) | .665 |
ACS | 32/205 (15.6) | −0.47 (−1.00 to 0.05) | .076 |
Transfusion | 111/205 (54.2) | 0.17 (−0.21 to 0.56) | .378 |
sVCAM-1 | |||
Dactylitis | 160/205 (78.1) | 0.02 (−0.15 to 0.18) | .845 |
VOC | 167/205 (81.5) | 0.02 (−0.16 to 0.20) | .853 |
ACS | 32/205 (15.6) | 0.02 (−0.17 to 0.21) | .819 |
Transfusion | 111/205 (54.2) | −0.04 (−0.17 to 0.10) | .607 |
sICAM-1 | |||
Dactylitis | 160/205 (78.1) | −0.001 (−0.13 to 0.13) | .990 |
VOC | 167/205 (81.5) | −0.10 (−0.24 to 0.04) | .167 |
ACS | 32/205 (15.6) | −0.004 (−0.15 to 0.14) | .955 |
Transfusion | 111/205 (54.2) | −0.05 (−0.15 to 0.06) | .397 |
sE-selectin | |||
Dactylitis | 160/205 (78.1) | 0.18 (0.04-0.32) | .014 |
VOC | 167/205 (81.5) | −0.01 (−0.17 to 0.14) | .859 |
ACS | 32/205 (15.6) | −0.07 (−0.24 to 0.10) | .407 |
Transfusion | 111/205 (54.2) | 0.14 (0.03-0.26) | .018 |
Linear regression models predicting natural log–transformed biomarker levels with sickle-related complications as predictor variables adjusting for child age and sex.
ACS, acute chest syndrome; CRP, C-reactive protein.
Significant after adjustment for 32 comparisons (models) as described in “Methods.”
Impact of hydroxyurea on endothelial activation and inflammation
To evaluate whether hydroxyurea was associated with changes in endothelial activation over the first year, we conducted longitudinal analyses to evaluate changes in plasma biomarker levels from treatment initiation (Figure 1). At enrollment there were no differences in levels of endothelial activation markers or inflammation by treatment arm (P > .05 for all). However, children randomized to hydroxyurea had lower levels of sE-selectin, P-selectin, and Angpt-2 than children treated with placebo, independent of child age and sex, with differences in markers (sE-selectin or P-selectin and Angpt-2) evident after 2 to 4 months of treatment, respectively, and sustained over 1-year follow-up (Figure 1).
Longitudinal changes in markers of endothelial activation over 1-year follow-up by trial arm. Scatterplot depicting levels of endothelial activation markers at baseline (0), and at 2-, 4-, and 12-months follow-up by treatment arm with hydroxyurea (red) and placebo (blue). Biomarkers measured included Angpt-2 (A), P-selectin (B), CRP (C), sVCAM-1 (D), sICAM-1 (E), and sE-selectin (F). Linear regression was used to plot the changes in markers over time by treatment arm, and the effect of hydroxyurea was evaluated using linear mixed effects models evaluating changes in longitudinal levels of the marker of endothelial activation (dependent variable) including a participant-specific random intercept, a categorical time variable (0, 2, 4, and 12 months), and adjusting for age at enrollment and sex. Angpt-2, P-selectin, and sE-selectin remain significant after adjustment for comparisons, as described in “Methods.”
Longitudinal changes in markers of endothelial activation over 1-year follow-up by trial arm. Scatterplot depicting levels of endothelial activation markers at baseline (0), and at 2-, 4-, and 12-months follow-up by treatment arm with hydroxyurea (red) and placebo (blue). Biomarkers measured included Angpt-2 (A), P-selectin (B), CRP (C), sVCAM-1 (D), sICAM-1 (E), and sE-selectin (F). Linear regression was used to plot the changes in markers over time by treatment arm, and the effect of hydroxyurea was evaluated using linear mixed effects models evaluating changes in longitudinal levels of the marker of endothelial activation (dependent variable) including a participant-specific random intercept, a categorical time variable (0, 2, 4, and 12 months), and adjusting for age at enrollment and sex. Angpt-2, P-selectin, and sE-selectin remain significant after adjustment for comparisons, as described in “Methods.”
Endothelial activation is associated with increased TCD velocities
Because the endothelium is an important regulator of vascular reactivity and function, we evaluated the relationships between markers of endothelial activation and TCD velocities after 10 to 12 months of study treatment. A log10 increase in enrollment values of both Angpt-2 (adjusted β coefficient, 31.5 cm/s; 95% confidence interval [CI], 13.4-49.7) and sE-selectin (adjusted β coefficient, 38.7 cm/s; 95% CI, 17.5-59.9) was independently associated with a higher TCD velocity (adjusted P < .05).
Of 161 children with TCD measurements between 10 and 12 months, 34 children (21.1%) had a conditional or abnormal TCD velocity (≥170 cm/s), which confers an increased risk of primary stroke.22 After adjustment for child age, sex, and hydroxyurea therapy, a log10 increase in Angpt-2 or C-reactive protein on enrollment was associated with increased odds of a conditional or abnormal TCD at 10 to 12 months, with an adjusted odds ratio of 6.42 (95% CI, 1.03-39.99) and 2.00 (95% CI, 1.01-3.93), respectively (P < .05 for both before correction). However, after correcting for multiple comparisons, these adjusted odds ratios lost significance.
Relationship between endothelial activation and inflammation and hematologic parameters
We evaluated the relationship between our markers of endothelial activation and inflammation and hematologic parameters at enrollment and 12-months follow-up with results stratified by treatment arm (Figure 2). sE-selectin and Angpt-2 were negatively correlated with Hb and HbF and positively correlated with WBC and neutrophil counts, whereas P-selectin was positively correlated with WBC, neutrophil, and platelet counts.
Longitudinal correlations between biomarkers of inflammation and endothelial activation and laboratory variables by treatment arm. Heat map constructed using the Spearman ρ of paired correlations across visits of markers of inflammation and endothelial activation (C-reactive protein [CRP], sICAM-1, P-selectin [P-sel], sVCAM-1, sE-selectin [sE-Sel], and Angpt-2) and hematologic variables including Hb, automated reticulocyte count (ARC), mean corpuscular volume (MCV), HbF %, WBC count, automated neutrophil count (ANC), and platelet count (Plt). Treatment group is indicated along the x-axis as hydroxyurea (H) or placebo (P). Scheduled follow-up visits are indicated along the y-axis as the number of months since enrollment (0, 2, 4, 12). Biomarkers that were significantly changed are indicated with a bracket, and an asterisk represents the P value in the linear mixed effects models in which ∗∗P < .01 and ∗∗∗P < .001. The cellular sources of the biomarkers measured are indicated below the graph with WBC, Plt, and endothelium (Endo) indicated, with a plus sign indicating that it is produced by the cell and a minus sign indicating no cellular production.
Longitudinal correlations between biomarkers of inflammation and endothelial activation and laboratory variables by treatment arm. Heat map constructed using the Spearman ρ of paired correlations across visits of markers of inflammation and endothelial activation (C-reactive protein [CRP], sICAM-1, P-selectin [P-sel], sVCAM-1, sE-selectin [sE-Sel], and Angpt-2) and hematologic variables including Hb, automated reticulocyte count (ARC), mean corpuscular volume (MCV), HbF %, WBC count, automated neutrophil count (ANC), and platelet count (Plt). Treatment group is indicated along the x-axis as hydroxyurea (H) or placebo (P). Scheduled follow-up visits are indicated along the y-axis as the number of months since enrollment (0, 2, 4, 12). Biomarkers that were significantly changed are indicated with a bracket, and an asterisk represents the P value in the linear mixed effects models in which ∗∗P < .01 and ∗∗∗P < .001. The cellular sources of the biomarkers measured are indicated below the graph with WBC, Plt, and endothelium (Endo) indicated, with a plus sign indicating that it is produced by the cell and a minus sign indicating no cellular production.
By leveraging the longitudinal collection of data over the course of the study, we were able to evaluate causal pathways associated with hydroxyurea effect on TCD velocity as a surrogate measure of stroke risk (Figure 3). We ultimately found 2 statistically significant pathways through which hydroxyurea may indirectly affect TCD velocities. The first pathway showed that hydroxyurea was associated with an increase in HbF at 2 months (β = 6.91; P < .001), which, in turn, was associated with an increased total Hb at 3 months (β = 0.09; P < .001), which, in turn, was associated with a decrease in transcranial velocities at 10 to 12 months (β = −7.33; P = .001). The second showed that hydroxyurea was associated with an increase in total Hb at 3 months (β = 0.43; P = .008) independent of HbF, which, in turn, was associated with a decrease in TCD velocities at 10 to 12 months (β = −7.33; P = .001). Apart from these 2 indirect pathways, hydroxyurea was not shown to have a direct effect on TCD velocities (β = 1.73; P = .698). Of note, multiple pathways involving angiopoiten-2 at 4 months approached statistical significance, including increased total Hb at 3 months associated with a decrease in angiopoiten-2 at 4 months (β = −0.17; P = .063) and decreased WBC count at 3 months associated with a decrease in angiopoietin-2 at 4 months (β = 0.03; P = .057).
Causal diagram showing hydroxyurea-mediated effect on endothelial activation and TCD velocities. Hydroxyurea (HU) fixed-dose treatment at 20 mg/kg was initiated at baseline (0 months); HbF was assessed at 2 months after HU treatment; total Hb and WBC count were measured at 3 months after the start of HU therapy; Angpt-2 was measured at 4 months after the start of study treatment; and TCD velocities measured between 10 and 12 months from the start of HU treatment. Associations between variables were derived from a structural equation modeling mediation model, presented as arrows. Significant relationships between variables are depicted using solid lines, whereas nonsignificant relationships are presented using dotted lines. Black dotted lines indicate that the relationship approached significance (P < .10). There was an indirect effect of hydroxyurea on TCD velocity through its sequential effect on HbF and Hb (indirect effect, −4.34; 95% CI, −7.35 to −1.33; P = .005) and through Hb alone (indirect effect, −3.14; 95% CI, −6.06 to −0.22; P = .035). Through these indirect pathways, hydroxyurea treatment had a total negative effect on TCD velocity (total effect, −8.55; 95% CI, –16.72 to −0.38; P = .040). †P < .10.
Causal diagram showing hydroxyurea-mediated effect on endothelial activation and TCD velocities. Hydroxyurea (HU) fixed-dose treatment at 20 mg/kg was initiated at baseline (0 months); HbF was assessed at 2 months after HU treatment; total Hb and WBC count were measured at 3 months after the start of HU therapy; Angpt-2 was measured at 4 months after the start of study treatment; and TCD velocities measured between 10 and 12 months from the start of HU treatment. Associations between variables were derived from a structural equation modeling mediation model, presented as arrows. Significant relationships between variables are depicted using solid lines, whereas nonsignificant relationships are presented using dotted lines. Black dotted lines indicate that the relationship approached significance (P < .10). There was an indirect effect of hydroxyurea on TCD velocity through its sequential effect on HbF and Hb (indirect effect, −4.34; 95% CI, −7.35 to −1.33; P = .005) and through Hb alone (indirect effect, −3.14; 95% CI, −6.06 to −0.22; P = .035). Through these indirect pathways, hydroxyurea treatment had a total negative effect on TCD velocity (total effect, −8.55; 95% CI, –16.72 to −0.38; P = .040). †P < .10.
Discussion
Hydroxyurea increases HbF production and decreases inflammation, but its impact on endothelial function is not well characterized. In this study, hydroxyurea decreased plasma levels of sE-selectin, P-selectin, and Angpt-2 in children with SCA; lower levels of Angpt-2 and sE-selectin at enrollment and 12-month follow-up were associated with decreased TCD velocities, a strong predictor of stroke. Mediation analysis suggested that hydroxyurea therapy decreased TCD velocities, in part, by increasing HbF, which then increases Hb, and that there may be some Angpt-2–mediated pathways involving Hb and WBC count. Together, these findings provide new insight into how hydroxyurea may decrease SCA-related complications, including stroke, and suggest that Angpt-2 may be both a surrogate marker for hydroxyurea therapeutic response and potentially a target for SCA therapies.
The endothelium is a highly dynamic cell layer that serves as a critical regulator of health and disease.23 Endothelial activation in SCA occurs through multiple pathways. Upregulated CAMs recruit immune cells to tissues leading to inflammation and adhesion of sickled erythrocytes to the endothelium, resulting in obstructed capillary flow, ischemia-reperfusion injury, and hemolysis.24,25 Angpt-2 is a proangiogenic and proinflammatory protein that destabilizes the endothelium. Systemic increases in Angpt-2 are associated with sepsis,26,27 retinopathy,28,29 diabetes,30-32 severe malaria,19,33 neurocognitive injury and decline,34-36 and ischemic stroke.37 In the present study, Angpt-2 levels at enrollment were associated with a history of prior transfusion, dactylitis, acute chest syndrome, and VOC. The relationship between Angpt-2 and SCA-related complications on enrollment is consistent with previous studies reporting elevated Angpt-2 levels in children with SCA in steady state and during VOC,38 in adults with SCA and a history of leg ulcers,39 and in malaria.40 Increases in Angpt-2 levels may reflect an underlying proangiogenic state triggered by ischemic insults, spanning from small-vessel complications (eg, dactylitis) to large-vessel complications (eg, stroke).
Several studies have evaluated hydroxyurea’s impact on endothelial function using in vitro or animal models of SCA with supporting evidence from plasma of people with SCA.4,41 In the HUSTLE study (ClinicalTrials.gov identifier: NCT00305175), hydroxyurea decreased levels of sE-selectin but not sICAM-1 or sVCAM-1.42 This effect was associated with reduced pneumococcal burden and improved survival in an experimental SCA mouse mice.42 Hydroxyurea was associated with reduced levels of sE-selectin and sVCAM-1 in a US pediatric population of children with SCA, and reduced Angpt-1 and vascular endothelial growth factor in adults with SCA.43,44 In this study, hydroxyurea decreased sE-selectin 2 months into treatment, and Angpt-2 and P-selectin 4 months into treatment, with these decreases persisting through 1 year follow-up. Consistent with the HUSTLE study, we did not observe differences in sVCAM-1 or sICAM-1. These data demonstrate that hydroxyurea selectively affects the endothelium and may provide insights into the mechanisms of protection.
Hydroxyurea may indirectly affect Angpt-2 levels by increasing Hb and HbF, resulting in less erythrocyte sickling and reducing WBC counts. However, hydroxyurea has also been shown to have direct effects on the endothelium through decreased neutrophil adhesion and increased nitric oxide bioavailability.45 E-selectin expression is specific to the endothelium and is transcriptionally induced by inflammatory cytokines (tumor necrosis factor α and interleukin-1β) and P-selectin. Murine data confirm a selective decrease in sE-selectin associated with hydroxyurea treatment, independent of changes in markers of inflammation, and associated with reduced neutrophil recruitment and tissue injury in pneumococcal infection.42 In contrast to E-selectin, P-selectin and Angpt-2 are stored in Weibel-Palade bodies and rapidly released in response to endothelial activation in response to exogenous stimuli (ie, tumor necrosis factor α and thrombin) and environmental cues (eg, hypoxia),11,46 and are inhibited in response to nitric oxide production.47
There was a strong association between Angpt-2 levels at enrollment and TCD velocities measured at 1 year follow-up; to our knowledge, the first such report in a pediatric SCA cohort. TCD measures blood flow velocity in the middle cerebral and internal carotid arteries, with increases in blood flow velocity associated with increased stroke risk in children, which may be attributable to hyperemia related to impaired vasodilatation.48,49 Hydroxyurea is associated with reduced TCD velocity and reduced stroke risk.17,50,51 In the NOHARM cohort, hydroxyurea was associated with a median TCD velocity reduction of 12 cm/s after 1 year of treatment.18 Reduced TCD velocities are thought to be driven by increased HbF and total Hb expression and myelosuppression (reduced leukocyte and platelet counts).45 Although our findings require validation, Angpt-2 levels represent a potential biomarker to risk stratify children in settings in which TCD testing is not available.
Mediation analysis suggests that hydroxyurea reduces TCD velocities, in part, through increased HbF and total Hb. Mediation analysis further suggested possible Angpt-2–mediated pathways involving both Hb and WBC count, which approached statistical significance. Mechanistically, we hypothesize that both the hydroxyurea-mediated myelosuppression and increased Hb/HbF results in improved oxygenation of red blood cells and less endothelial activation. It is important to note that hydroxyurea may also reduce TCD velocities by directly affecting the endothelium.
Study strengths include the randomized design and serial measures of endothelial activation. We identified possible causal pathways, suggesting that hydroxyurea’s impact on Angpt-2 may be mediated through increased Hb. Although this mediation analysis has limitations on implied causality, we posit a biologically plausible mechanism through which hydroxyurea-mediated increases in Hb, decrease Angpt-2 levels and stroke risk. More study is needed to validate these potential pathways and explore other possible mechanisms.
In conclusion, in this study we found that hydroxyurea treatment in children with SCA leads to persistently decreased markers of endothelial activation. These data suggest that serial measures of Angpt-2 could be used to measure clinical response to hydroxyurea therapy, may serve as a risk stratification marker for stroke, and could be an additional target for SCA therapies.
Acknowledgments
The authors thank the study participants and their families who volunteered to participate in this study, the Mulago Hospital Sickle Cell Clinic, Global Health Uganda, and the Data Coordinating Center at Cincinnati Children’s Hospital Medical Center.
The Doris Duke Charitable Foundation provided funds for the study (ICRA 2013139), and Addmedica Inc (Paris, France) donated hydroxyurea (Siklos) for the NOHARM trial. Biomarker evaluation was supported by a Doris Duke Charitable Foundation fellowship (A.C.) and a Fogarty International Center Global Health fellow grant (D43 TW009345 to F.S.).
Authorship
Contribution: R.O.O., R.E.W., and C.C.J. designed the study, supervised the trial, and analyzed the results; A.L.C. designed and conducted the biomarker testing with A.C. and T.F.S.; T.L. coordinated many critical aspects of the trial to ensure its safe and successful operational execution and participated in the analysis; M.N. conducted the TCD analysis; A.L., J.M.S., K.M., and A.L.C. performed statistical analyses; T.F.S. and A.L.C. wrote the first draft of the manuscript; and all authors participated in the editing of the manuscript and approved the final version.
Conflict-of-interest disclosure: R.E.W. is a consultant for Nova Laboratories; receives research support (drug donation) from Bristol Myers Squibb and Addmedica; and serves on 2 data and safety monitoring boards. A.L.C. is an inventor on a patent related to life-threatening responses to infection that includes angiopoietin-2. None of these disclosures are relevant to the results and conclusions of the NOHARM trial. The remaining authors declare no competing financial interests.
Correspondence: Andrea L. Conroy, Indiana University School of Medicine, Ryan White Center for Pediatric Infectious Diseases and Global Health, 1044 W. Walnut St, Indianapolis, IN 46202; email: [email protected].
References
Author notes
Original data may be obtained by sending a request including proposed analysis and rationale for analysis to the corresponding author, Andrea L. Conroy ([email protected]); individual participant data will not be shared.
The online version of this article contains a data supplement.