Research article | Open | Open Peer Review | Published:
Association of plasma F2-isoprostanes and isofurans concentrations with erythropoiesis-stimulating agent resistance in maintenance hemodialysis patients
BMC Nephrologyvolume 16, Article number: 79 (2015)
In patients undergoing maintenance hemodialysis (HD), hyporesponsiveness to erythropoiesis stimulating agents (ESAs) is associated with adverse clinical outcomes. Systemic inflammation is highly prevalent in HD patients and is associated with ESA hyporesponsiveness. Oxidative stress is also highly prevalent in HD patients, but no previous study has determined its association with ESA response. This study assessed the association of plasma markers of oxidative stress and inflammation with ESA resistance in patients undergoing maintenance HD.
We analyzed data from 165 patients enrolled in the Provision of Antioxidant Therapy in Hemodialysis study, a randomized controlled trial evaluating antioxidant therapy in prevalent HD patients. Linear and mixed-effects regression were used to assess the association of baseline and time-averaged high sensitivity F2-isoprostanes, isofurans, C-reactive protein (hsCRP), and interleukin-6 (IL-6) with ESA resistance index (ERI), defined as the weekly weight-adjusted ESA dose divided by blood hemoglobin level. Unadjusted models as well as models adjusted for potential confounders were examined. Predicted changes in ERI per month over study follow-up among baseline biomarker quartiles were also assessed.
Patients with time-averaged isofurans in the highest quartile had higher adjusted mean ERI compared with patients in the lowest quartile (β = 14.9 ng/ml; 95 % CI 7.70, 22.2; reference group <0.26 ng/ml). The highest quartiles of hsCRP and IL-6 were also associated with higher adjusted mean ERI (β = 10.8 mg/l; 95 % CI 3.52, 18.1 for hsCRP; β = 10.2 pg/ml; 95 % CI 2.98, 17.5 for IL-6). No significant association of F2-isoprostanes concentrations with ERI was observed. Analyses restricted to baseline exposures and ERI showed similar results. Baseline hsCRP, IL-6, and isofurans concentrations in the highest quartiles were associated with greater predicted change in ERI over study follow-up compared to the lowest quartiles (P = 0.008, P = 0.004, and P = 0.04, respectively). There was no association between baseline F2-isoprostanes quartile and change in ERI.
In conclusion, higher concentrations of isofurans, hsCRP and IL-6, but not F2-isoprostanes, were associated with greater resistance to ESAs in prevalent HD patients. Further research is needed to test whether interventions that successfully decrease oxidative stress and inflammation in patients undergoing maintenance HD improve ESA responsiveness.
Anemia leading to treatment with erythropoiesis-stimulating agents (ESAs) is highly prevalent among patients with end stage renal disease (ESRD) and both the severity of anemia and variability in hemoglobin have been associated with poor clinical outcomes in patients on maintenance dialysis [1–4]. Hyporesponsiveness to ESAs – manifested as an increase in the ESA dose requirement to maintain a constant blood hemoglobin level, a drop in blood hemoglobin on a constant ESA dose, or a failure to achieve a target hemoglobin level on very high doses of ESAs – is also associated with increased risk of adverse outcomes [5–7]. Furthermore, ESA therapy is costly, which is further amplified in patients with ESA hyporesponsiveness .
The etiology of ESA resistance in patients undergoing maintenance hemodialysis (HD) is complex and incompletely understood. One likely contributor is chronic systemic inflammation, which has been found to be present in up to 80 % of HD patients [9, 10]. Several prior studies have identified an association between ESA resistance and blood levels of systemic inflammatory markers such as C-reactive protein and interleukin-6 (IL-6) [11–13]. Furthermore, small-scale investigations of anti-inflammatory interventions have suggested a potential benefit on improving ESA responsiveness [14, 15].
Growing evidence suggests that an exaggerated oxidative stress burden may also play a role in inducing ESA resistance in patients with ESRD . Plasma F2-isoprostanes, an established measure of oxidative stress, has been found to be up to four times higher in patients undergoing maintenance HD compared with healthy subjects . Isofurans, a more recently described biomarker of oxidative stress, has been linked to increased risk for cardiovascular disease and kidney dysfunction [18, 19]. Oxidative stress may promote ESA resistance by increasing red blood cell fragility, decreasing red cell life span, and/or reducing endogenous erythropoietin synthesis [20, 21]. No prior study has examined the association of plasma F2-isoprostanes or isofurans with ESA resistance in HD patients. Given the poor clinical outcomes associated with ESA hyporesponsiveness, there is a compelling need to investigate oxidative stress as a novel pathway leading to resistance to ESAs.
Prior studies have shown that oxidative stress and inflammation are synergistically linked. During an inflammatory response, stimulated leukocytes release superoxide and myeloperoxidase, catalyzing oxidative injury and in turn leading to further release of pro-inflammatory cytokines and acute phase proteins through activation of the transcription nuclear factor κB . However, whether oxidative stress and inflammation may act independently to promote ESA hyporesponsiveness in patients with ESRD remains unknown. The primary objective of this study was to test the hypothesis that higher concentrations of inflammatory and oxidative stress biomarkers are each independent associate with greater resistance to ESAs in patients undergoing maintenance HD.
Patients were identified from the Provision of Antioxidant Therapy in Hemodialysis (PATH) study, a randomized clinical trial comparing the effect of combination antioxidant therapy with mixed tocopherols plus α-Lipoic acid (ALA) versus matched placebo on biomarkers of inflammation and oxidative stress in patients undergoing maintenance dialysis . Inclusion criteria for the PATH study included ESRD receiving thrice-weekly maintenance HD for at least 120 days, age >18 years, and life expectancy >1 year. Exclusion criteria included AIDS; active malignancy; recent or anticipated kidney transplant; hospitalization over the past 30 days; vitamin E, vitamin C, or anti-inflammatory medications over the past 30 days; and history of poor adherence to hemodialysis. Of the 353 participants randomized in the PATH trial, we identified the 165 individuals with completed baseline measurements of F2-isoprostanes, isofurans, high-sensitivity C-reactive protein (hsCRP), interleukin-6 (IL-6), and who were treated with an ESA at all study visits, for inclusion in the present study. All identified patients completed the study. Baseline patient characteristics for participants included in the present study compared to those excluded are presented in Additional file 1: Table S1. The PATH trial was approved by the Western and Vanderbilt University Medical Center institutional review boards, and all patients provided written informed consent before study enrollment.
Patients underwent study visits at baseline and then monthly for 6 months. Baseline assessment included collection of demographic characteristics, body mass index (BMI), comorbid conditions, etiology of ESRD, type of vascular access, and blood for routine chemistries. In addition, blood was obtained at baseline and at month 6 for measurement of biomarkers of inflammation and oxidative stress, which represented the exposures of interest. Oxidative stress was quantified by simultaneous measurement of plasma F2-isoprostane and isofurans by gas chromatography–mass spectrometry analysis . Inflammation was quantified by measurement of hsCRP using ELISA (Diagnostic Systems Laboratories, Webster, TX) and IL-6 cytokine concentrations also using ELISA (BioSource International, Camarillo, CA).
Patients underwent blood collection for measurement of blood hemoglobin at baseline, and at each monthly follow-up visit. In addition, at each study visit, each patient’s last ESA dose and average hemodialysis frequency over the prior month was ascertained. The primary outcome was the ESA resistance index (ERI), calculated as the weekly weight-adjusted ESA dose (units per kilogram of body weight per week) divided by the hemoglobin level (grams per deciliter). For each patient, the ERI was calculated at baseline, and at each monthly follow-up visit.
Baseline data were summarized, stratified by quartile of ERI. Data are presented as mean ± standard deviation, medians with interquartile ranges, and proportions as appropriate. Data were complete for all baseline covariates. Scatter plots were created and Pearson’s correlation coefficients calculated to summarize the relationship among the primary exposure inflammatory and oxidative stress biomarkers. For both the oxidative stress markers (F2-isoprostanes and isofurans) and the inflammatory markers (hsCRP and IL-6), baseline measurements of exposure variables were grouped into quartiles. Time-averaged exposures were generated by taking the mean of the baseline and 6-month follow-up measurements, and were also grouped into quartiles for analysis. Two main analyses were performed. First, multivariable linear regression was performed to determine the relationship of quartiles of baseline F2-isoprostanes, isofurans, hsCRP, and IL-6 with baseline ERI. Referent categories for comparisons were the lowest quartile of each baseline exposure variable: F2-isoprostanes of <0.049 ng/ml, isofurans of <0.25 ng/ml, hsCRP of <19.8 mg/l, and IL-6 of <8.3 pg/ml. Second, the relationship between quartiles of time-averaged oxidative stress and inflammatory markers with mean ERI over 6-month follow-up was assessed using a multilevel mixed-effects linear regression model with unstructured covariance accounting for correlation among repeated measures within a patient. Referent categories for comparisons were the lowest quartile of each time-averaged exposure variable: F2-isoprostanes of <0.049 ng/ml, isofurans of <0.26 ng/ml, hsCRP of <19.8 mg/l, and IL-6 of <8.3 pg/ml.
Linear and mixed effects regressions using biomarker quartiles as the exposure were selected for the main analysis after model diagnostics were applied to regressions using exposures as continuous variables and substantial violations of key model assumptions were identified. Specifically, plots of residuals versus fitted values were examined and revealed both curvature and unequal variances for all four exposures. Additionally, plots of sorted residuals versus theoretical quantiles of the normal distribution revealed evidence of non-normality for all four exposures. Outliers were also identified, but deletion diagnostics did not reveal substantially influential observations. To address these concerns, log transformations of exposures and ERI were performed, and regressions repeated using robust standard errors. These changes did not satisfactorily address the violation of linearity. Thus, for the primary analyses, regression using exposure quartiles was performed to account for non-linearity.
For each primary analysis, three models were examined: (1) unadjusted; (2) adjusted for potential demographic and laboratory confounders including for age, sex, race, BMI, etiology of ESRD, diabetes status, transferrin saturation (serum iron/total iron binding capacity, single-pool Kt/V, and serum parathyroid hormone); and (3) adjusted for model 2 covariates plus log-transformed hsCRP (for analyses using oxidative stress markers as the exposures), or log-transformed isofurans (for analyses using inflammatory markers as the exposures). Additionally, for analyses using mixed-effects linear models, study visit number and randomly assigned treatment group were added as covariates. Adjusted restricted cubic splines with four degrees of freedom were constructed to illustrate the relationships between baseline measures of oxidative stress and inflammation with baseline ERI. Covariates included in the adjusted spline models were the same as for regression model 3. For illustrative purposes, splines and 95 % confidence intervals were superimposed on scatter plots of baseline biomarker values (x-axis) and ERI (y-axis), where each dot represents an individual patient. Three patients with outlying values for IL-6 > 100 pg/ml were excluded only for the purposes of illustrative spline construction. For spline figures, a logarithmic scale was used on the x-axis given the right-skewed distribution of oxidative stress and inflammatory markers.
A secondary analysis was performed to determine the association of quartiles of baseline markers of inflammation and oxidative stress with predicted change in ERI over the 6-month study follow-up. This analysis was conducted using a generalized estimating equation model including interaction terms for exposure quartile and study visit, and was conducted on cohort of 253 participants with non-missing baseline oxidative stress and inflammatory markers, including subjects not treated with ESA during at study visit. Adjustment covariates were identical to primary analysis regression model 3. Statistical comparisons were made between quartiles 4 and 1 for each biomarker exposure. Additionally, statistical significance of the linear trend was over all quartiles was determined.
All analyses were performed using Stata, version 11 (StataCorp LP, College Station, TX; www.stata.com).
Baseline characteristics of the study cohort, stratified by quartile of baseline ERI are summarized in Table 1. Compared to patients with baseline ERI <5.49, patients with baseline ERI >19.2 were more likely to be African American, less likely to dialyze using an arteriovenous fistula compared with a graft, and had higher prevalence of cardiovascular disease and greater maintenance HD duration. Additionally, patients with higher baseline ERI had lower serum iron and lower transferrin saturation.
Baseline F2-isoprostane concentrations were positively correlated with hsCRP (r = 0.30) (Fig. 1). In contrast, isofurans concentrations were only very weakly correlated with hsCRP (r = 0.09), and were not correlated with IL-6 levels (r = 0.01). Plots illustrating the relationship of F2-isoprostanes with isofurans and hsCRP with IL-6 are shown in Additional file 1: Figure S1.
Association of markers of oxidative stress with ERI
Table 2 shows the associations of baseline F2-isoprostanes and isofurans with baseline ERI. In unadjusted analyses, patients in the highest quartile of isofurans concentrations (>0.71 ng/ml) had significantly higher ERI compared to patients with concentrations <0.25 ng/ml (β = 7.49; 95 % CI 1.14, 13.8; P = 0.02). This association persisted following adjustment for demographic variables and transferrin saturation (Model 2), as well as following further adjustment for hsCRP (Model 3). Additionally, patients in the second quartile (0.25–< 0.45 ng/ml) but not the third quartile (0.45–< 0.73 ng/ml) also had significantly higher ERI relative to those in the lowest quartile. This association remained significant across adjustment models, though the significance was attenuated with greater adjustment. In contrast, no significant association between quartile of F2-isoprostanes and ERI was observed, regardless of level of adjustment. Figure 2 illustrates the continuous adjusted relationships between baseline concentrations of F2-isoprostanes and isofurans with ERI.
Table 3 shows the associations of time-averaged markers of oxidative stress with mean ESA resistance over follow-up, using linear mixed-effects models. The results of the analyses using time-averaged isofurans as the exposure variable were not substantially different than those using only the baseline measurements, with the second and highest quartiles of isofurans concentration associated with higher mean ERI (reference group: patients with time-averaged isofurans concentration <0.26 ng/ml), regardless of the specific adjustment model. In contrast to the analyses using baseline measurements, patients with time-averaged F2-isoprostanes in the highest quartile had significantly higher mean ERI compared to patients with concentrations <0.049 ng/ml (β = 9.82; 95 % CI 1.78, 17.9; P = 0.02) with adjustment for baseline demographic and laboratory variables. However, following full adjustment for potential confounders including for hsCRP (Model 3), this association was no longer significant.
Comparing predicted change in ERI per month over six-month study follow-up, patients with baseline isofurans concentrations >0.73 ng/ml had a significantly higher predicted monthly ERI increase compared to patients with isofurans <0.25 n/gml (6.26 ng/ml/month; 95 % CI 0.19, 12.33; P = 0.04) (Table 4). In contrast, there was no significant difference among quartiles of baseline F2-isoprostanes concentrations with respect to monthly ERI change. Neither oxidative stress marker was associated with a linear trend in ERI change across quartiles.
Association of markers of inflammation with ERI
In unadjusted analyses, similar to the observed associations between isofurans and ERI, patients with baseline hsCRP ≥101 mg/l had higher mean baseline ERI compared with patients with hsCRP <19.8 mg/dl (β = 8.45; 95 % CI 1.63, 15.3; P = 0.01) (Table 5). Conversely, patients with baseline hsCRP falling into the middle two quartiles did not have significantly higher baseline ERI compared with patients with hsCRP in the lowest quartile. Patients with baseline IL-6 concentrations in the highest quartile of ≥20.9 pg/ml, but not those with measurements in the middle two quartiles, had significantly higher ERI compared to patients with IL-6 < 8.3 pg/ml (β = 9.15; 95 % CI 2.78, 15.5; P = 0.005). Upon adjustment for demographic factors and baseline transferrin saturation, and further adjustment for oxidative stress as measured by isofurans concentration, these observed associations remained statistically significant for IL-6, but were only of borderline significance for hsCRP (Table 5).
Table 6 presents the results of analyses assessing the association of time-averaged markers of inflammation with mean ERI over follow-up, using linear mixed-effects models. For time-averaged IL-6, the results were substantially similar to the analyses using only baseline measurements, with higher mean ERI observed in patients with IL-6 in the highest quartile, regardless of the level of adjustment. In contrast, for time-averaged hsCRP, patients with marker concentrations in the highest quartile had significantly higher mean ERI not only in the unadjusted model, but also in the fully adjusted model (β = 10.8 mg/l; 95 % CI 3.52, 18.1; P = 0.004).
Figure 3 illustrates the continuous unadjusted relationship between hsCRP and IL-6 with ERI. Examination of the restricted cubic splines further supported the use of quantile regression rather than continuous linear regression, as the relationship between inflammatory markers and ERI was curvilinear. Similar to the results of the regression analyses, an association between greater inflammation and higher ERI was evident only at higher concentrations of hsCRP and IL-6.
Comparing predicted change in ERI per month over six-month study follow-up, patients with baseline hsCRP concentrations >101 mg/l had a significantly higher predicted monthly ERI increase compared to patients with hsCRP <19.8 mg/l (8.27 mg/l/month; 95 % CI 2.17, 14.38; P = 0.008) (Table 4). Similarly, patients with baseline IL-6 in the highest quartile had significantly greater monthly change in ERI over follow-up compared to patients in the lowest quartile (P = 0.004). Neither inflammatory marker was associated with a linear trend in ERI change across quartiles.
Two sets of sensitivity analyses were performed. In the first set, the study cohort was expanded to include participants of PATH who had one or more study visits during which they were untreated with ESAs. Inclusion of these additional participants increased the total number of subjects from 165 to 253. Linear regressions assessing the association between baseline biomarker concentrations and baseline ERI, and linear mixed-effects regressions assessing the association between time-averaged exposures and mean ERI were repeated with the expanded cohort. The results of these analyses were similar to those obtained in the primary analysis (Additional file 1: Tables S2–S5). On exception was that patients with baseline IL-6 in the highest quartile had significantly higher ERI compared to those in the lowest quartile only in an unadusted analysis; following full adjustment, this association was no longer presents (Additional file 1: Table S4).
In the second set of sensitivity analyses, only patients who were treated with ESA at each study visit, and who were assigned to the control group in the PATH trial were included in the study cohort. This restricted cohort to 87 participants was then analyzed using linear mixed effects regression to determine the association of time-averaged inflammatory and oxidative stress exposures with ERI. This approach was employed to eliminate the possible impact of antioxidant treatment on the time-averaged exposures. The results of this set of analyses are reported in Additional file 1: Tables S6 and S7. Patients with time-averaged hsCRP in the top two quartiles had significantly greater mean ERI compared to those in the lowest quartile (β = 11.9 mg/l; 95 % CI 2.91, 20.8; P = 0.01). There were no significant associations between quartile of IL-6, F2-isoprostanes or isofurans and ERI.
In this study examining the association of markers of inflammation and oxidative stress on ESA resistance in patients undergoing maintenance dialysis, we found that higher concentrations of hsCRP, IL-6 and isofurans, but not F2-isoprostanes, were associated with increased resistance to ESAs, as measured by ERI. These associations were evident in unadjusted models, and also following adjustment for potential demographic and laboratory confounders. Additionally, these associations were not linear, with higher ERI only observed for patients with inflammatory biomarkers in the top quartile, and for patients with isofurans in the second and fourth but not the third quartile. Finally, the association between increased concentrations of isofurans and ERI remained significant following adjustment for hsCRP, and that for the markers of inflammation and ERI remained significant with adjustment for isofurans. This observation suggests possible distinct biological mechanisms leading from inflammation versus oxidative stress to ESA resistance.
Multiple mechanisms may explain how systemic inflammation leads to inhibition of erythropoiesis and thus to resistance to ESAs in patients undergoing maintenance HD. Patient with ESRD have many contributory factors to systemic inflammation, including uremic toxin buildup, dialysate microcontamination, exposure to indwelling catheters and synthetic grafts, and high prevalence of comorbid diabetes and atherosclerotic disease . Inflammation has in turn been shown to lead to activation of macrophages with removal of circulating erythrocytes, down-regulation of erythropoietin receptors, and dysregulated cell signaling following activation of the erythropoietin receptor [26, 27]. Additionally, inflammatory cytokines have been demonstrated to decrease erythropoietin production, suppress proliferation of red cell precursors in the bone marrow, and decrease release of mature erythrocytes into the circulation [28–30].
The data on the contribution of oxidative stress to anemia in patients with ESRD undergoing maintenance dialysis are far more limited. Oxidative stress has been shown to be associated with anemia through downregulation of erythropoietin synthesis, decreased red cell production, and increased erythrocyte fragility. In particular, human erythrocyte nitric oxide synthesis and production of peroxynitrite are activated by oxidative stress leading to upregulation of superoxide dismutase, increased erythrocyte membrane rigidity, and shortening of red cell lifespan [21, 31]. In clinical studies, lipid peroxidation of the erythrocyte cell membrane has been positively associated with markers of red cell destruction, and with resistance to administered ESAs in patients undergoing maintenance HD . Additionally, elevated circulating 8-hydroxy-2’-deoxyguanosine, a product of oxidized DNA, has been associated with greater ESA resistance in HD patients . Furthermore, oxidative stress activates cationic erythrocyte membrane channels permeable to ionic calcium, anionic channels permeable to chloride ions, and apoptotic caspases, all of which may lead to upregulation of programmed erythrocyte cell death [33, 34]. In our study, we found that isofurans were poorly correlated with both hsCRP and IL-6, and there remained a predictive relationship between isofurans and ERI even after adjustment for hsCRP as a biomarker of inflammation. These observations provide evidence that the pathway leading from oxidative stress to ESA resistance is likely distinct from the more well-characterized pathway linking inflammation with ESA resistance. Our findings are in line with those of other studies that have identified pathways leading to ESA resistance that are independent of inflammation. Kalim et al. have recently demonstrated the association of non-enzymatic carbamylation of serum albumin with erythropoietin resistance in ESRD patients, even after adjustment for serum markers of systemic inflammation .
In our study, we found an unexpected non-linear relationship between quartiles of isofurans and ESA resistance that was particular notable in the time-averaged analyses. Specifically, we observed higher ERI in patients with isofurans falling not only in the highest quartile, but also in the middle-low quartile. No significant difference was found among patients with isofurans in the third quartile. Certainly one possible explanation for this finding is a statistical one, namely absence of a sufficient sample size to provide adequate power to detect a difference of 5.5 ng/ml in isofurans concentration between the first and third quartiles. Indeed, the β coefficient estimates for the first and third quartiles were 7.3 versus 5.5 ng/ml, numerically quite similar. Alternatively, a true non-linear relationship may exist. Studies of DNA mutagens that act via oxidative stress have shown non-linear dose-effect relationships, with multiple threshold “plateaus” that correspond to the action of cellular defense mechanisms such as reactive oxygen species scavengers and antioxidant enzymes . Similarly, a non-linear relationship between cigarette smoke exposure and cardiovascular risk has also been demonstrated . Thus, in our study, stepping from the second to the third quartile of oxidative stress could theoretically trigger an additional protective mechanism that may reinstate the balance beteween oxidants and antioxidants. In the case of ESA resistance in patients undergoing maintenance HD, such a trigger could partially restore ESA responsiveness. Stepping to the highest level of oxidative stress may overwhelm cellular defense mechanisms leading to ESA resistance. Given that such an explanation remains speculative, further basic and translational research is needed to elucidate the connection between oxidative stress, anemia and ESA resistance.
We observed that high concentrations of isofurans but not F2-isoprostanes were associated with higher ERI. This difference in the predictive role of these two products of free radical-induced peroxidation of arachadonic acid on the degree of ESA resistance was an unexpected finding. Prior studies have shown that relative formation of lipid peroxidation products, including F2-isoprostanes and isofurans, is profoundly influenced by ambient oxygen concentration . In particular, as oxygen concentration increases, F2-isoprostanes formation becomes disfavored, and isofurans formation becomes more favored. In highly oxygenated tissues, such as the central nervous system and the kidneys, isofurans concentrations can exceed F2-isoprostanes levels by greater than two-fold . At the cellular level, elevated levels of reactive oxygen species, as have been demonstrated in ESRD, have been shown to be associated with impaired mitochondrial respiration [39, 40]. Impaired mitochondrial function with blockade of the electron transport chain may then lead to impaired oxygen utilization and buildup of intracellular oxygen molecules with concomitant increased intracellular oxygen concentrations. Independent of blood hemoglobin levels, such increased oxygen concentrations favor isofurans production over F2-isoprostanes production. Furthermore, impaired signaling along the erythropoietin axis is associated with decreased mitochondrial biogenesis, and greater mitochondrial dysfunction, also leading to impaired oxygen utilization, buildup of oxygen tension, and subsequent predominance of isofurans [41–43]. These observations together are one possible explanation for why we observed an association between ERI and isofurans, but not F2-isoprostanes. However, this explanation is theoretical and further research is needed to determine the cellular pathways leading to differential oxidative stress biomarker production.
Our study has multiple strengths. The PATH trial is the largest published trial of antioxidant therapy conducted to date in patients undergoing maintenance HD. All analyzed patients thus had excellent follow-up and were well characterized. Additionally, the lipid peroxidation biomarkers (F2-isoprostanes and isofurans) were analyzed at a reference quality laboratory using gas chromatography/mass spectrometry, thus providing gold standard measures of oxidative stress. Furthermore, repeated measurements of erythropoietin dose response over a 6-month period allowed use of a linear mixed-effects model to incorporate longitudinal data on ESA resistance, increasing the reliability of the observed associations.
Our study also has some limitations. First, though we identified a number of important associations with ESA resistance, it is unclear whether these are causal. Second, the use of data from a clinical trial population may decrease the generalizability of the observed findings. Third, though the results from the PATH trial showed no significant impact of antioxidant treatment status on ESA resistance over study follow-up, it is possible that treatment status may impact the association between inflammatory or oxidative stress markers and ERI. We have attempted to account for this by controlling for trial group assignment in our multivariable regression models, by the possibility of residual confounding remains. Fourth, patients enrolled in the PATH trial were slightly younger and had a higher mean BMI than the general United States ESRD population, which may limit generalizeabilty. Recent work, however, has shown that age does not predict ESA response in either pediatric or adult patients undergoing dialysis . Additionally, though BMI is associated with ESA dose in HD patients, the correlation is weak, and calculation of ERI incorporates body weight [45, 46]. Additionally, we have adjusted for both age and BMI in our final adjustment models. A final limitation is that given the observational nature of our study, there remains the possibility of residual confounding by unmeasured variables on the association between the measured biomarkers inflammation and oxidative stress with ESA resistance.
In conclusion, we found that in patients undergoing maintenance HD, higher concentrations of isofurans, hsCRP, and IL-6 but not F2-isoprostanes, were associated with increased resistance to erythropoiesis-stimulating agents. The association between isofurans and ERI persisted even after adjustment for levels of inflammation. These findings support that oxidative stress is an important contributor to ESA resistance in the ESRD population, and suggest a pathway independent of systemic inflammation. Future studies should test whether interventions that successfully decrease oxidative stress and inflammation in patients undergoing maintenance dialysis lead to reductions in ESA dose requirements.
Erythropoiesis-stimulating agent resistance index
End stage renal disease
High-sensitivity C-reactive protein
Body mass index
Enzyme-linked immunosorbent assay
U.S. Renal Data System. USRDS 2013 Annual Data Report: Atlas of Chronic Kidney Disease and End Stage Renal Disease in the United States. Bethesda, MD: National Institutes of Health, National Institute of Diabetes and Digestive and Kidney Diseases; 2013.
Brookhart MA, Schneeweiss S, Avorn J, Bradbury BD, Liu J, Winkelmayer WC. Comparative mortality risk of anemia management practices in incident hemodialysis patients. JAMA. 2010;303:857–64.
Molnar MZ, Mehrotra R, Duong U, Kovesdy CP, Kalantar-Zadeh K. Association of hemoglobin and survival in peritoneal dialysis patients. Clin J Am Soc Nephrol. 2011;6:1973–81.
Regidor DL, Kopple JD, Kovesdy CP, Kilpatrick RD, McAllister CJ, Aronovitz J, et al. Associations between changes in hemoglobin and administered erythropoiesis-stimulating agent and survival in hemodialysis patients. J Am Soc Nephrol. 2006;17:1181–91.
Duong U, Kalantar-Zadeh K, Molnar MZ, Zaritsky JJ, Teitelbaum I, Kovesdy CP, et al. Mortality associated with dose response of erythropoiesis-stimulating agents in hemodialysis versus peritoneal dialysis patients. Am J Nephrol. 2012;35:198–208.
Panichi V, Rosati A, Bigazzi R, Paoletti S, Mantuano E, Beati S, et al. Anaemia and resistance to erythropoiesis-stimulating agents as prognostic factors in haemodialysis patients: results from the RISCAVID study. Nephrol Dial Transplant. 2011;26:2641–8.
Solomon SD, Uno H, Lewis EF, Eckardt K-U, Lin J, Burdmann EA, et al. Erythropoietic response and outcomes in kidney disease and type 2 diabetes. N Engl J Med. 2010;363:1146–55.
Charytan C. Bundled-rate legislation for medicare reimbursement for dialysis services: implications for anemia management with ESAs. Clin J Am Soc Nephrol. 2010;5:2355–62.
Stenvinkel P, Heimbürger O, Paultre F, Diczfalusy U, Wang T, Berglund L, et al. Strong association between malnutrition, inflammation, and atherosclerosis in chronic renal failure. Kidney Int. 1999;55:1899–911.
Kaysen GA, Dubin JA, Müller H-G, Mitch WE, Rosales LM, Levin NW. Relationships among inflammation nutrition and physiologic mechanisms establishing albumin levels in hemodialysis patients. Kidney Int. 2002;61:2240–9.
Rattanasompattikul M, Molnar MZ, Zaritsky JJ, Hatamizadeh P, Jing J, Norris KC, et al. Association of malnutrition-inflammation complex and responsiveness to erythropoiesis-stimulating agents in long-term hemodialysis patients. Nephrol Dial Transplant. 2013;28:1936–45.
Schneider A, Schneider MP, Scharnagl H, Jardine AG, Wanner C, Drechsler C. Predicting erythropoietin resistance in hemodialysis patients with type 2 diabetes. BMC Nephrol. 2013;14:67.
Gunnell J, Yeun JY, Depner TA, Kaysen GA. Acute-phase response predicts erythropoietin resistance in hemodialysis and peritoneal dialysis patients. Am J Kidney Dis. 1999;33:63–72.
Sirken G, Kung S-C, Raja R. Decreased erythropoietin requirements in maintenance hemodialysis patients with statin therapy. ASAIO J Am Soc Artif Intern Organs. 2003;49:422–5.
Cooper A, Mikhail A, Lethbridge MW, Kemeny DM, Macdougall IC. Pentoxifylline improves hemoglobin levels in patients with erythropoietin-resistant anemia in renal failure. J Am Soc Nephrol. 2004;15:1877–82.
Himmelfarb J, Stenvinkel P, Ikizler TA, Hakim RM. The elephant in uremia: oxidant stress as a unifying concept of cardiovascular disease in uremia. Kidney Int. 2002;62:1524–38.
Ikizler TA, Morrow JD, Roberts LJ, Evanson JA, Becker B, Hakim RM, et al. Plasma F2-isoprostane levels are elevated in chronic hemodialysis patients. Clin Nephrol. 2002;58:190–7.
Alkazemi D, Egeland GM, Roberts LJ, Kubow S. Isoprostanes and isofurans as non-traditional risk factors for cardiovascular disease among Canadian Inuit. Free Radic Res. 2012;46:1258–66.
Ware LB, Fessel JP, May AK, Roberts LJ. Plasma biomarkers of oxidant stress and development of organ failure in severe sepsis. Shock Augusta Ga. 2011;36:12–7.
Gallucci MT, Lubrano R, Meloni C, Morosetti M, Manca di Villahermosa S, Scoppi P, et al. Red blood cell membrane lipid peroxidation and resistance to erythropoietin therapy in hemodialysis patients. Clin Nephrol. 1999;52:239–45.
Suliman HB, Ali M, Piantadosi CA. Superoxide dismutase-3 promotes full expression of the EPO response to hypoxia. Blood. 2004;104:43–50.
Himmelfarb J. Linking oxidative stress and inflammation in kidney disease: which is the chicken and which is the egg? Semin Dial. 2004;17:449–54.
Himmelfarb J, Ikizler TA, Ellis C, Wu P, Shintani A, Dalal S, et al. Provision of Antioxidant Therapy in Hemodialysis (PATH): A Randomized Clinical Trial. J Am Soc Nephrol. 2013;25:623–33.
Milne GL, Gao B, Terry ES, Zackert WE, Sanchez SC. Measurement of F2- isoprostanes and isofurans using gas chromatography–mass spectrometry. Free Radic Biol Med. 2013;59:36–44 [Methods in Lipid Oxidation].
Bamgbola OF. Pattern of resistance to erythropoietin-stimulating agents in chronic kidney disease. Kidney Int. 2011;80:464–74.
Macdougall IC, Cooper AC. Erythropoietin resistance: the role of inflammation and pro‐inflammatory cytokines. Nephrol Dial Transplant. 2002;17 suppl 11:39–43.
Nemeth E, Rivera S, Gabayan V, Keller C, Taudorf S, Pedersen BK, et al. IL-6 mediates hypoferremia of inflammation by inducing the synthesis of the iron regulatory hormone hepcidin. J Clin Invest. 2004;113:1271–6.
Jelkmann W. Proinflammatory cytokines lowering erythropoietin production. J Interferon Cytokine Res. 1998;18:555–9.
Allen DA, Breen C, Yaqoob MM, Macdougall IC. Inhibition of CFU-E colony formation in uremic patients with inflammatory disease: role of IFN-gamma and TNF-alpha. J Investig Med. 1999;47:204–11.
Dallalio G, Law E, Robert T, Means J. Hepcidin inhibits in vitro erythroid colony formation at reduced erythropoietin concentrations. Blood. 2006;107:2702–4.
Deliconstantinos G, Villiotou V, Stavrides JC, Salemes N, Gogas J. Nitric oxide and peroxynitrite production by human erythrocytes: a causative factor of toxic anemia in breast cancer patients. Anticancer Res. 1995;15:1435–46.
Kato A, Odamaki M, Hishida A. Blood 8‐hydroxy‐2′‐deoxyguanosine is associated with erythropoietin resistance in haemodialysis patients. Nephrol Dial Transplant. 2003;18:931–6.
Lang F, Qadri SM. Mechanisms and significance of eryptosis, the suicidal death of erythrocytes. Blood Purif. 2012;33:125–30.
Abed M, Artunc F, Alzoubi K, Honisch S, Baumann D, Föller M, et al. Suicidal erythrocyte death in end-stage renal disease. J Mol Med. 2014;92:871–9.
Kalim S, Tamez H, Wenger J, Ankers E, Trottier CA, Deferio JJ, et al. Carbamylation of Serum Albumin and Erythropoietin Resistance in End Stage Kidney Disease. Clin J Am Soc Nephrol. 2013;8:1927–34.
Platel A, Nesslany F, Gervais V, Marzin D. Study of oxidative DNA damage in TK6 human lymphoblastoid cells by use of the in vitro micronucleus test: determination of no-observed-effect levels. Mutat Res. 2009;678:30–7.
Ambrose JA, Barua RS. The pathophysiology of cigarette smoking and cardiovascular disease: an update. J Am Coll Cardiol. 2004;43:1731–7.
Fessel JP, Jackson Roberts L. Isofurans: novel products of lipid peroxidation that define the occurrence of oxidant injury in settings of elevated oxygen tension. Antioxid Redox Signal. 2004;7:202–9.
Rowley S, Liang L-P, Fulton R, Shimizu T, Day B, Patel M. Mitochondrial respiration deficits driven by reactive oxygen species in experimental temporal lobe epilepsy. Neurobiol Dis. 2015;75:151–8.
Victor VM, Rovira-Llopis S, Saiz-Alarcon V, Sangüesa MC, Rojo-Bofill L, Bañuls C, et al. Altered mitochondrial function and oxidative stress in leukocytes of anorexia nervosa patients. PLoS ONE. 2014;9, e106463.
Mille-Hamard L, Billat VL, Henry E, Bonnamy B, Joly F, Benech P, et al. Skeletal muscle alterations and exercise performance decrease in erythropoietin-deficient mice: a comparative study. BMC Med Genomics. 2012;5:29.
Stoyanoff TR, Todaro JS, Aguirre MV, Zimmermann MC, Brandan NC. Amelioration of lipopolysaccharide-induced acute kidney injury by erythropoietin: involvement of mitochondria-regulated apoptosis. Toxicology. 2014;318:13–21.
Wu C-L, Chen S-D, Yin J-H, Hwang C-S, Yang D-I. Erythropoietin and sonic hedgehog mediate the neuroprotective effects of brain-derived neurotrophic factor against mitochondrial inhibition. Neurobiol Dis. 2010;40:146–54.
Bamgbola OF, Kaskel FJ, Coco M. Analyses of age, gender and other risk factors of erythropoietin resistance in pediatric and adult dialysis cohorts. Pediatr Nephrol. 2009;24:571–9.
Locatelli F, Andrulli S, Memoli B, Maffei C, Vecchio LD, Aterini S, et al. Nutritional-inflammation status and resistance to erythropoietin therapy in haemodialysis patients. Nephrol Dial Transplant. 2006;21:991–8.
Vega A, Ruiz C, Abad S, Quiroga B, Velázquez K, Yuste C, et al. Body composition affects the response to erythropoiesis-stimulating agents in patients with chronic kidney disease in dialysis. Ren Fail. 2014;36:1073–7.
This work was supported by a Clinical Translational Science Award 1UL-1RR024975 from the National Center for Research Resources, K24 DK 62849 and T32DK007467 from the National Institute of Diabetes and Digestive and Kidney Diseases, R01 HL070938 from the National Heart, Lung, and Blood Institute and Center in Molecular Toxicology, P30 ES000267 from National Institute of Environmental Health Sciences, and an unrestricted gift from the Northwest Kidney Centers Foundation.
The authors declare that they have no competing interests.
MR designed the analytic plan, performed the statistical analysis, and drafted the manuscript. TI and JH conceived of the study, participated in its design and coordination, and helped draft the manuscript. CE participated in the design and coordination of the study, and provided critical review of the manuscript. RM contributed to design of the analytic plan, and helped draft the manuscript. All authors read and approved the final manuscript.
Baseline Characteristics of Participants Included Compared to Those Excluded. Table S2. Association of quartiles of baseline markers of oxidative stress with baseline ESA resistance index (ERI) in maintenance hemodialysis patients, including subjects not treated with ESA during one or more study visits during follow-up (n = 253). Table S3. Association of quartiles of time-averaged markers oxidative stress with mean ESA resistance index (ERI) over 6-month follow-up in maintenance hemodialysis patients, including subjects not treated with ESA during one or more study visits during follow-up (n = 253). Table S4. Association of quartiles of baseline markers of inflammation with baseline ESA resistance index (ERI) in maintenance hemodialysis patients, including subjects not treated with ESA during one or more study visits during follow-up (n = 253). Table S5. Association of quartiles of time-averaged markers of inflammation with mean ESA resistance index (ERI) over 6-month follow-up in maintenance hemodialysis patients, including subjects not treated with ESA during one or more study visits during follow-up (n = 253). Table S6. Association of quartiles of time-averaged markers oxidative stress with mean ESA resistance index (ERI) over 6-month follow-up in maintenance hemodialysis patients, excluding subjects in PATH antioxidant treatment arm (n = 87). Table S7. Association of quartiles of time-averaged markers of inflammation with mean ESA resistance index (ERI) over 6-month follow-up in maintenance hemodialysis patients, excluding subjects in PATH antioxidant treatment arm (n = 87). Figure S1. Correlation of baseline inflammatory and oxidative stress markers in maintenance hemodialysis patients.