Severe sickle cell anemia is associated with increased plasma levels of TNF‐R1 and VCAM‐1
Authors: Daniel A. Dworkis, Elizabeth S. Klings, Nadia Solovieff, Guihua Li, Jacqueline N. Milton, Stephen W. Hartley, Efthymia Melista, Jason Parente, Paola Sebastiani, Martin H. Steinberg, Clinton T. Baldwin
Journal: American Journal of Hematology (2010)

Sickle cell anemia (SCA, HBB glu6val) is characterized by multiple complications and a high degree of phenotypic variability: some subjects have only sporadic pain crises and few acute hospitalizations, while others experience multiple serious complications, high levels of morbidity, and accelerated mortality [1]. The tumor necrosis factor-α (TNF-α) signaling pathway plays important roles in inflammation and the immune response; variation in this pathway might be expected to modify the overall severity of SCA through the pathway's effects on the vascular endothelium [2, 3]. We examined plasma biomarkers of TNF-α activity and endothelial cell activation for associations with SCA severity in 24 adults (12 mild, 12 severe). Two biomarkers, tumor necrosis factor-α receptor-1 (TNF-R1) and vascular cell adhesion molecule-1 (VCAM-1) were significantly higher in subjects with severe SCA. Along with these biomarker differences, we also examined data from a genome-wide association study (GWAS) using SCA severity as a disease phenotype, and found evidence of genetic association between disease severity and a single nucleotide polymorphism (SNP) in VCAM1, which codes for VCAM-1, and several SNPs in ARFGEF2, a gene involved in TNF-R1 release [4]. Exposure to TNF-α activates endothelial cells, leading to up-regulation and enhanced cell-surface expression of multiple proteins including VCAM-1 and intracellular adhesion molecule-1 (ICAM-1) [5]. During acute infections, these proteins facilitate the attachment of leukocytes to the vessel wall, ultimately resulting in transmigration into infected tissue [6]. In SCA, variable TNF-α signaling might alter disease severity by activating systemic and pulmonary endothelium, and by modulating rates of sickle vasoocclusion secondary to altered leukocyte-endothelial cell adhesion [7, 8]. Associations between some individual clinical complications of SCA and plasma levels of VCAM-1 and ICAM-1, but not TNF-R1, have been described previously [9, 10]. Here, we report associations between increased plasma levels of VCAM-1 and TNF-R1 and severe SCA, as measured by the SCA severity score, a validated metric which incorporates information about steady-state laboratory data and medical history into a Bayesian network in order to estimate a subject's probability of dying within the next five years due to complications of SCA [11, 12]. In a population of 52 adult subjects with SCA recruited for a study of pulmonary hypertension, we used the severity score to identify the 12 subjects with the most mild clinical presentations (probability of death within five years = 7.5 ± 2.6%, mean ± SD) and the 12 subjects with the most severe presentations (probability = 80.3 ± 18%). No significant differences were observed between the two groups in terms of age or hydroxyurea treatment (P-values 0.35 and 0.68, respectively), while non-significant trends toward being female and having lower glomerular filtration rates (GFR) were observed in the severe group (P-values 0.1, 0.1). Subjects in the severe group had lower levels of hemoglobin than subjects in the mild group (median 8.3 vs. 9.45 g/dL, P-value = 0.046), though the difference was not significant after adjusting for multiple testing. To look for associations between SCA severity and TNF-α axis activity and endothelial cell activation, steady-state plasma samples from these two groups were tested for levels of E-selectin, ICAM-1, nitric oxide (NO), P-selectin, TNF-α, TNF-R1, and VCAM-1. Only TNF-R1 and VCAM-1 showed significant differences between the two groups after correcting for multiple testing; both biomarkers were significantly higher in severe cases of SCA (P-value = 0.0015 and P-value = 0.0019, respectively, Table I). Plasma levels of TNF-R1 and VCAM-1 showed a significant correlation after log transformation (Figure 1, P-value = 0.036). Comparison of plasma VCAM-1 and TNF-R1. Scatter-plot of log-transformed plasma levels of VCAM-1 and TNF-R1 among 24 subjects with SCA. Full circles represent subjects in the severe group (n = 12), and empty boxes represent subjects in the mild group (n = 12). Dotted lines mark the median log-transformed values for all 24 subjects. Renal disease, a common occurrence in SCA, can lead to increased plasma levels of TNF-R1 secondary to decreased renal clearance [13, 14]. To evaluate this and other possible confounding factors, we tested for associations between TNF-R1 or VCAM-1 levels and GFR, gender, age, and hydroxyurea treatment status. We also looked for potential association with plasma nitric oxide (NO), since endothelial cell-surface expression of leukocyte adhesion proteins is altered by NO bioavailibility [15-17]. TNF-R1 levels were found to be associated with gender (P-value = 0.048), but were not significantly associated with GFR, hydroxyurea use, NO levels, or with age (P-value >0.05 in all cases). VCAM-1 levels were not significantly associated with any of these factors (P-values >0.05). The relationship between TNF-R1 and SCA severity remained significant in a linear regression model involving both TNF-R1 and gender (P-value = 4.6 × 10−03), indicating that TNF-R1 as a plasma biomarker is significantly related to increased SCA severity beyond its association with gender. In contrast to our results, multiple groups have shown that hydroxyurea therapy was associated with decreased levels of multiple inflammatory biomarkers, including VCAM-1; this may be due to our small sample size, or to differences in efficacy of hydroxyurea treatment [9, 18, 19]. We identified 155 SNPs in 16 candidate genes which either code for one of the six biomarkers, or are involved in their release into the plasma. SNP data was obtained from a previous GWAS, which examined 1265 subjects with SCA from the Cooperative Study of Sickle Cell Disease (CSSCD) using the severity score as a phenotype. While none of the 155 SNPs met stringent criteria for genome-wide significance, several of these SNPs did show significant association with SCA severity at moderate levels (BF >3; Table II) [4]. Of the 47 SNPs examined in the six candidate genes coding for the plasma biomarkers (ICAM1, SELE, SELP, TNF, TNFRSF1A, and VCAM1), only one (rs1041163), located in the promoter of VCAM1, was significantly associated with SCA severity (BF = 3.5). This same SNP was previously reported to be associated with small vessel stroke in the CSSCD [20]. We also examined 108 SNPs in 10 genes involved in the release of these biomarkers or in the generation of NO (ADAM17, ARFGEF2, ARF1, ARF3, ERAP-1, IL1B, NOS1, NOS2A, NOS3, and NUCB2) [21-25]. Of these SNPs, the top association was with rs2273102, located in intron 12 of ARFGEF2, a gene known to be involved in endothelial release of TNF-R1 into circulation in an ADAM17-independent manner (BF = 70 under a dominant model) [24]. Three of the other 15 SNPs in ARFGEF2 (rs6019548, rs1115535, and rs6019566) also showed significant associations with SCA severity (BF >3). Using SNP set enrichment (SSE) analysis, which summarizes genetic association across multiple SNPs within a single gene, ARFGEF2 as a whole was found to be significantly related to SCA severity (P-value = 5.5 × 10−05) [4]. Significant association with SCA severity was also observed with rs1143634, located in the gene IL1B (BF >30 under a recessive model). IL1B codes for the cytokine IL-1B, which induces ADAM17-dependent release of TNF-R1 into circulation, among other roles [21]. None of the other 15 genes were significantly related to SCA severity by SSE analysis. Our data showing severity-related variability in plasma levels of biomarkers of TNF-α signaling, together with genetic association between SCA severity and genes involved in biomarker generation and release, supports the hypothesis that genetic variability along the TNF-α signaling pathway may modulate the severity of SCA. Further investigation, including the confirmation of both the genetic and biomarker associations with SCA severity in a single population, is necessary to develop a more complete understanding of the ways in which genetic variation in TNF-α signaling might affect the clinical manifestations of this disease. With such confirmation, the associations shown here between SCA severity and variability along the TNF-α signaling pathway might be used to further prospectively stratify subjects with SCA into different levels of risk, and, ultimately, may lead to the ability to more rationally designed therapeutic plans for subjects with different clinical presentations of SCA. Subjects included in the biomarker analysis were recruited as part of a study conducted at Boston Medical Center (BMC) to identify genetic modifiers associated with pulmonary hypertension in sickle cell disease [10]. Peripheral blood samples were collected when the subjects were at least two weeks removed from a vasoocclusive crisis, and at least four weeks removed from an episode of acute chest syndrome. Platelet-poor plasma samples were prepared as described previously, and stored at −70°C until the time of the study [10]. Of the 52 subjects in the BMC group available at the time of this study, the 24 with the most extreme severity scores were selected for complete biomarker analysis. For these 24 subjects, GFR was estimated using the Modification of Diet in Renal Disease formula adjusted for use in African-Americans [26]. In one case, GFR estimation was not possible, as a creatinine level was not recorded at the time of sample collection. Subjects included in the genetic association study were described fully in a previous report [4]. Briefly, the genetic association population set was comprised of 1265 subjects with SCA (1088 mild and 177 severe), originally recruited as part of the Cooperative Study of Sickle Cell Disease (CSSCD) [27]. These studies were approved by the Institutional Review Board at Boston Medical Center. Severity score was calculated using a Bayesian network described previously, and available for download at http://www.bu.edu/sicklecell/downloads/Projects [11]. For subjects in the biomarker analysis, severity score was calculated at study enrollment, coincident with phlebotomy for plasma samples. For subjects in the CSSCD, severity scores were calculated using the median values of data acquired over multiple longitudinal visits, with thresholds for defining severe cases of SCA described previously [11, 14]. TNF-α, TNF-R1, and P-selectin were assayed by enzyme-linked immunoassay (ELISA) using commercially available kits from R & D Biosystems (Minneapolis, MN). VCAM-1, ICAM-1, and E-selectin levels were assayed by ELISA using kits from Invitrogen (Carlsbad, CA). Nitric oxide metabolites were assayed as nitrate/nitrite levels using a colorimetric assay (Cayman Chemical, Ann Arbor, MI). Univariate associations between individual biomarker levels and SCA severity, gender, or hydroxyurea therapy status were examined using Mann-Whitney U tests. Correlations between biomarker levels and age or estimated GFR were examined using Spearman correlation tests or Pearson correlation for log-transformed values. An overall significance level of 0.05 was set as significant for all univariate tests, and a Bonferonni correction was used to adjust for multiple testing (Bonferonni corrected P-value = 0.0029). For multivariate analysis, biomarker levels were natural-log-transformed and fit to linear models involving SCA severity with or without other variables, as described in the text. All univariate tests and linear models were performed in R (available at www.r-project.org). DNA isolated from peripheral blood samples from the 1265 subjects in the genetic association population was analyzed for >600k SNPs using the Illumina Human 610-Quad chip, according to the manufacturer's protocol. Potential associations between SCA severity and individual SNPs were tested using Bayesian methods described in detail elsewhere [28, 29]. SNPs with Bayes Factor (BF) > 3 under any model were considered to show evidence of association with SCA severity with a false positive rate (FPR) of ∼5%; SNPs with BF >10 have a FPR of ∼1%, and SNPs with BF >1,000 were considered to show association at a level of genome-wide significance [4, 30]. In addition, associations between SCA severity and individual genes were tested using the SSE score, which uses a hypergeometric distribution to assess the probability that multiple SNPs within a single gene are all related to a phenotype [4]. Work funded by NIH grant HL R01068970 and Hematology training grant ST32 HL 007501. Daniel A. Dworkis* †, Elizabeth S. Klings*, Nadia Solovieff‡, Guihua Li*, Jacqueline N. Milton‡, Stephen W. Hartley‡, Efthymia Melista* †, Jason Parente†, Paola Sebastiani‡, Martin H. Steinberg*, Clinton T. Baldwin* †, * Department of Medicine, Boston University School of Medicine, Boston, Massachusetts, † Center for Human Genetics, Boston University School of Medicine, Boston, Massachusetts, ‡ Department of Biostatistics, Boston University School of Public Health, Boston, Massachusetts.