Abstract
Background/Objective: Calcium loading has been associated with cardiovascular risk in hemodialysis (HD) patients. However, it remains to be elucidated whether alterations of intradialytic calcium buffering add to the increased cardiovascular disease burden in this high-risk population. Methods: Intradialytic calcium kinetics was evaluated in a cross-sectional observational study by measuring dialysate-sided ionized calcium mass balance (iCaMB), calcium buffer capacity, and change in serum calcium levels in 40 chronic HD patients during a routine HD session. A dialysate calcium of 3.5 mEq/L was used to adequately challenge calcium buffer mechanisms. Aortic pulse wave velocity and serum osteocalcin levels were measured prior to the HD session. Presence of cardiovascular disease and diabetes was assessed. Results: The mean dialysate-sided iCaMB, extracellular fluid ionized calcium mass gain, and buffered ionized calcium mass were 469 (±154), 111 (±49), and 358 (±145) mg/HD, respectively. The mean ionized serum calcium increase (∆iCa) was 0.42 (±0.14) mEq/L per HD. The mean intradialytic calcium buffer capacity was 73 (±18)%. Multivariate regression analysis revealed significant independent association of (1) iCaMB with the dialysate-to-blood calcium gradient at HD start and (2) intradialytic calcium buffer capacity with undercarboxylated osteocalcin. The presence of coronary heart disease was associated with higher ∆iCa but not iCaMB in the multivariate model. Conclusions: In line with our proof-of-concept study, we provide clinical evidence for a rapidly accessible and exchangeable calcium pool involved in intradialytic calcium regulation and for the role of osteocalcin as a potential biomarker. Our findings argue for evaluating the prognostic potential of intradialytic calcium kinetics in prospective clinical trials.
Introduction
Current KDIGO guidelines [1] recommend a dialysate calcium (dCa) concentration of 2.5–3 mEq/L in order to avoid excessive intradialytic calcium (Ca) loading as this has been associated with an elevated cardiovascular (CV) risk [2]. However, we and others have previously demonstrated that even under these precautions calcium loading during dialysis occurs in a relevant number of chronic hemodialysis (HD) patients. For example, Gotch et al. [3] calculated that up to 40% of chronic HD patients are still in a positive overall Ca mass balance (CaMB) with a dCa concentration of 2.5 mEq/L when considering both inter- and intradialytic balance. However, even without taking interdialytic calcium uptake into account, a substantial interindividual variability of intradialytic CaMB ranging from −500 up to +400 mg/HD session has been observed with the use of 2.5 mEq/L dCa [4-6]. As the main determinant for calcium transfer during HD is the predialysis serum calcium concentration, these differences are most likely due to the varying use of active vitamin D and/or calcimimetic drugs but different CaMB assessment approaches have also to be taken into consideration [7, 8]. By utilizing a dialysate-sided ionized CaMB (iCaMB) assessment approach, we previously demonstrated that the use of both 2.5 and 3.5 mEq/L dCa regularly results in a highly positive intradialytic iCaMB, which is primarily predicted by the predialysis dialysate-to-blood ionized Ca gradient [4]. The CaMB was not reflected by changes in serum Ca concentration corroborating previous findings that the latter is not a reliable marker of Ca burden [9]. This is because a rapidly accessible and exchangeable Ca pool counteracts fluctuations of serum Ca deviations on a minute-to-minute basis in the event of high intradialytic Ca loading independent of calcitropic hormones, intestinal and renal Ca handling as well as cellular bone remodeling. This concept has been previously postulated [10] and corroborated for the first time in vivo in our clinical proof-of-concept study by reporting individual Ca buffer capacity in 28 chronic HD patients, ranging from 60 to 100% of whole individual intradialytic Ca load [4]. The lack of sensitivity of changes in serum Ca concentration and differences in the individual ability to remove excess Ca from the extracellular fluid (ECF) compartment implicate that dialysate-sided iCaMB assessment is less error-prone compared with blood-sampling based CaMB assessment strategies but rather requires partial or complete dialysate fluid collection.
If indeed an altered intradialytic Ca buffering has clinical consequences, such as a higher associated CV risk, is not completely clear. For example, Yamada et al. [11] previously demonstrated that intradialytic change in serum total Ca (∆tCa) is associated with aortic calcification in HD patients and Tagawa et al. [12] recently showed that higher dCa is associated with myocardial infarction among diabetic HD patients, which was statistically explained through an effect on ∆tCa. However, in the absence of CaMB assessment in both studies, it remains to be elucidated whether high intradialytic Ca loading and/or impaired Ca buffering per se increase CV risk. We, thus, conducted a cross-sectional, observational clinical study assessing intradialytic Ca kinetics (i.e., dialysate-sided iCaMB, intradialytic Ca buffer capacity, and change in ionized and total serum Ca levels) and its association with the prevalence of CV disease in 40 chronic HD patients.
Materials and Methods
Study Population and Procedures
Fourty chronic end-stage renal disease patients who underwent bicarbonate HD for more than 1 month were recruited to undergo noninvasive assessment of intradialytic Ca kinetics (i.e., dialysate-sided iCaMB, intradialytic Ca buffer capacity, and change in serum Ca levels; see definitions below), serum osteocalcin measurement as well as assessment of aortic pulse wave velocity (PWV) and blood- and pulse-pressure during a regular HD session. Exclusion criteria were inability or unwillingness to provide written informed consent, age below 18 years, and predialysis serum ionized Ca (iCa) concentration above 2.7 mEq/L; however, no patient was excluded based on these a priori established criteria. Study patients were not part of our previous proof-of-concept study [4]. HD treatment followed patients’ routine prescriptions with a dCa set to 3.5 mEq/L for a single HD session in order to adequately challenge the intradialytic Ca buffer mechanisms. Aiming at a constant ultrafiltration during the session, only mid- and end of the week HD sessions were used for study procedures. The vascular access was evaluated by measuring recirculation using the dialysis device (Fresenius 5008). Only subjects with recirculation values below 10% were qualified for study participation.
Blood- and Dialysate-Based Measurements
To assess intradialytic Ca kinetics, blood was drawn before the dialyzer (i.e., “precapillary”) prior to the beginning and then every 30 min until the end of each HD session to determine the respective precapillary iCa concentrations. In parallel, iCa concentrations were determined after 15, 30, 45, 60 min, and then every 30 min until the end of the dialysis session both in the fresh and spent dialysate. Ionized dCa in fresh and spent dialysate as well as precapillary serum iCa concentrations were analyzed by using an on-site blood gas analysis device (ABL 800 Flex, Drott Medizintechnik GmbH, Wiener Neudorf, Austria) with the appropriate settings for dialysate and serum specimen according to the manufacturer’s instructions. In addition, serum samples were collected before the HD session in each study participant and stored at −80°C for ELISA-based osteocalcin analysis. Serum carboxylated osteocalcin (cOC) and undercarboxylated osteocalcin (ucOC) were measured using respective high-sensitive enzyme immunoassay kits (Takara Bio Europe/Clontech, Saint-Germain-en-Laye, France) according to the manufacturer’s instructions. Serum total Ca (tCa) was measured using a colorimetric NM-BAPTA chromophore-based assay (Roche Diagnostics GmBH, Mannheim, Germany) by the Central Institute for Medical and Chemical Laboratory Diagnostics (ZIMCL) at the Medical University Innsbruck/University Hospital Innsbruck.
Assessment of Dialysate-Sided iCaMB
iCaMB was assessed according to a previously published approach [4, 5, 13], which estimates net Ca mass differences between the inlet and outlet dialysate tubing for each 30-min HD session time period and includes both the diffusive and convective component of iCaMB:
where iCaMB (mEq/HD) is the ionized calcium mass balance; dCa mean in (mEq/L) is the mean precapillary dCa over the dialysis session; QD (L/min) is the dialysate flow rate; t (min) is the duration of dialysis session; k is the postcapillary dCa (mEq/L) measured at 0 min (k = 1), 30 min (k = 2), 60 min (k = 3), 90 min (k = 4), 120 min (k = 5), 150 min (k = 6), 180 min (k = 7), 210 min (k = 8), and 240 min (k = 9); UF (L/min) is the ultrafiltration rate; and CaMB (mg/HD) = CaMB in mEq/HD × 80.16.
Calculation of Intradialytic Calcium Buffer Capacity
Intradialytic Ca buffer capacity was determined by setting the intradialytic change in ECF ionized Ca mass (∆iCaECF) in relation to iCaMB:
Intradialytic calcium buffer capacity (% of whole Ca load/HD session) = (1 – ∆iCAECF/iCAMB) × 100
The intradialytic increase in serum ionized Ca (iCa) concentration was used for calculation of ∆iCaECF (mg/HD). Similar iCa kinetics in ECF and serum was assumed according to Andersen et al. who showed only a small difference (approx. 5%) between serum and interstitial fluid iCa concentrations [14]:
∆iCaECF (mg/HD) = post-HD iCa (mg/L) × post-HD ECF (L) – pre-HD iCa (mg/L) × pre-HD ECF (L)
ECF volume was calculated at the beginning and end of the HD session based on patient weight before and after the HD session as 20% of patient weight before and after the HD session, respectively, to account for ultrafiltration-based ECF volume changes. Of note, the cellular compartment does not add to the distribution volume of iCa as an extremely low intracellular Ca concentration is mandatory for electrophysiologic reason.
Calculation of Intradialytic Change in Serum Ca Levels
The intradialytic change in serum iCa (∆iCa) and tCa (∆tCa) levels was calculated as follows:
∆iCa (mEq/L) = iCa post-HD (mEq/L) – iCa pre-HD (mEq/L)
∆tCa (mEq/L) = tCa post-HD (mEq/L) – tCa pre-HD (mEq/L)
Data Collection for Cross-Sectional Analysis
Clinical parameters of vascular calcification, namely, aortic PWV and peripheral blood and pulse pressure, were measured noninvasively before the midweek dialysis session in a calm environment with the patient in supine position after 5 min of supine rest. Brachial systolic blood pressure (SBP) and diastolic blood pressure (DBP) were calculated as the mean values of a minimum of 3 measurements taken at least 3 min apart on the nonfistula arm using a validated sphygmomanometer (HEM-780-D, Omron, Kyoto, Japan). Aortic PWV was measured by an oscillometric method (Mobil-O-Graph 24 h PWA device, IEM GmbH, Stolberg, Germany) utilizing an ARCsolver algorithm-based pulse wave analysis [15].
Ultrafiltration rate, dialysate flow rate, effective blood flow, and duration of dialysis session were obtained from the HD device (Fresenius 5008). Medical treatment (type and quantity of phosphate binders used, type and quantity of vitamin D substitution, and quantity of calcimimetic drugs) over the last 6 months or since dialysis onset in case of more recent dialysis start, comorbidities (diabetes status, smoking status, coronary heart disease, congestive heart failure, and peripheral arterial occlusive disease), demographic parameters (e.g., age and gender), and routinely collected laboratory parameters (e.g., hemoglobin, hematocrit, pH, serum bicarbonate, serum albumin, serum protein, serum phosphate, serum Ca, and serum PTH levels) were obtained from our local computerized clinical documentation systems and documented using a standardized case report form.
Definition of Comorbidities and Clinical Events
Coronary heart disease (CHD)/coronary artery disease = history of myocardial infarction and/or coronary artery bypass graft and/or coronary artery angioplasty.
Peripheral arterial occlusive disease (PAOD) = history of documented Fontaine Stage II, III, or IV.
Congestive heart failure (CHF) = history of hospital admission due to nonpneumopathological dyspnea and/or acute cardiogenic pulmonary edema documented on chest X-ray and left ventricular ejection fraction (LVEF) < 40% as documented by echocardiography.
Diabetes mellitus (DM) = history of diabetes type I or II and use of a hypoglycemic agent or insulin.
Statistical Analysis
A priori sample size calculation was based on the results of our previous proof-of-concept study in 28 HD patients [4]: with a reported mean intradialytic calcium buffer capacity of 78 (±7)%, 40 patients would allow for a 95% CI of ±2.2%. The results of intradialytic Ca kinetics assessment and baseline characteristics are presented as mean values (±SD) or absolute frequencies. Associations among intradialytic Ca kinetics parameters and osteocalcin levels, CV disease, comorbidities, laboratory and demographic parameters, and medical treatment as well as clinical parameters of vascular calcification (PWV and peripheral blood and pulse pressure) were assessed by means of correlation and regression analysis. Parametric distribution of the data was assessed by the use of Kolmogorov-Smirnov test. Pearson’s correlation coefficient (r) or Kendall’s Tau-b coefficient were computed for parametric and nonparametric correlation of metric variables, respectively. The nonparametric Mann-Whitney U test was used to compare the distribution of metric variables among dichotomized nominal parameters. To identify independent predictors of intradialytic Ca kinetics and CV disease, only parameters that showed a significant bivariate correlation (metric variables) or a significantly different distribution (nominal variables) were entered into a multiple linear (metric variables) or logistic (nominal variables) regression model, respectively. The results are presented as standardized beta coefficients (β) or regression coefficient B, respectively. Statistical analysis was performed with SPSS version 24.0. The level of significance (p value) was set to 0.05. Adjustment for multiple testing was not performed.
Results
Intradialytic Ca kinetics was assessed in 40 chronic HD patients during a routine dialysis session (see Table 1 for patient cohort characteristics). Using 3.5 mEq/L dCa allowed to adequately challenge intradialytic Ca buffer mechanisms and led to a positive iCaMB in all patients with a Ca load of 469 (±154) mg/HD on average. The mean intradialytic ∆iCa and ∆tCa was 0.42 (±0.14) and 0.9 (±0.32) mEq/L, respectively (Table 2). With a calculated mean ∆iCaECF of 111 (±49) mg/HD, the averaged amount of buffered iCa mass was 358 (±145) mg/HD. Thus, the mean intradialytic iCa buffer capacity (i.e., relation between ∆iCaECF and iCaMB) was 73 (±18)%. In order to identify predictors of intradialytic Ca kinetics parameters, univariate correlation analysis was conducted: iCaMB correlated positively with the dialysate-to-blood iCa gradient at HD start (p < 0.05) and negatively with serum bicarbonate at HD start and mean hemoglobin (Hb) levels (p < 0.01 each). For intradialytic Ca buffer capacity, we found a positive correlation with serum ucOC (p < 0.05) and a negative correlation with dry weight (p < 0.01). ∆iCa was positively associated with the dialysate-to-blood iCa gradient at HD start (p < 0.01) and negatively associated with serum bicarbonate at HD end (p < 0.05) (Table 3).
These predictors of intradialytic calcium kinetics were entered in a multiple linear regression model: Here, we could demonstrate an independent association of iCaMB with measured dialysate-to-blood iCa gradient at HD start (β = 0.341 and p = 0.01), that is, a higher dialysate-to-blood iCa gradient relates to higher iCaMB. In addition, iCaMB was independently associated with serum bicarbonate at HD start (β = −0.361 and p = 0.007) and mean Hb levels (β = 0.439 and p = 0.001). Intradialytic Ca buffer capacity was found to be independently associated with ucOC (β = 0.461 and p = 0.004) and dry weight (β = −0.398 and p = 0.011). ∆iCa was independently associated with the dialysate-to-blood iCa gradient at dialysis start only (β = 0.772 and p < 0.001) (Table 4).
We next aimed at identifying potential associations between calcium kinetics parameters and CV endpoints (CHD, PAOD, CHF, and CVD) as well as DM status by applying univariate correlation analysis: The presence of CHD was associated with higher iCaMB and higher ∆iCa (p < 0.05 each). The presence of DM was associated with lower intradialytic Ca buffer capacity (p < 0.01) (Table 5). Aortic PWV, an established clinical parameter of vascular calcification, correlated negatively with serum cOC (p < 0.05) (Table 3). Multiple logistic regression analysis showed an independent positive association of ∆iCa with the presence of CHD (B = 12.8 and p = 0.04); however, this did not hold true for iCaMB (Table 6). No further significant associations were found between calcium kinetics parameters, CV endpoints, and any other reported patients characteristics (Table 1).
Discussion
The present cross-sectional clinical study assessed intradialytic Ca kinetics (i.e., dialysate-sided iCaMB, intradialytic Ca buffer capacity, ∆iCa, and ∆tCa) and its association with CV disease in 40 chronic HD patients. We demonstrate that iCaMB is independently determined by the dialysate-to-blood gradient at HD start and provide robust evidence that a rapidly accessible and exchangeable Ca pool exists in vivo that acutely counteracts serum Ca deviations in HD patients in the event of high intradialytic CaMB: using 3.5 mEq/L dCa, the mean intradialytic Ca buffer capacity was found to be 73 (±18)%, which is strikingly similar to a mathematical modeling approach by Gotch et al. demonstrating that approximately 76% of positive CaMB is acutely buffered in chronic HD patients [3, 16] and analog to the observed 78% in our previous proof-of-concept study [4]. In line with these findings, the use of a novel single pool variable-volume calcium kinetics model proposed by di Filippo et al. [8, 17] suggested the existence of intradialytic calcium fluxes from the ECF compartment to a non-modeled pool. By demonstrating an independent positive association of ∆iCa with the presence of CHD in chronic HD patients, we present novel evidence for the clinical relevance of intradialytic Ca kinetics assessment in vivo. Importantly, a significant negative correlation was found between ∆iCa and intradialytic Ca buffer capacity (p < 0.01) in our study. So far, evidence for a potential role of ∆Ca as CV risk factor is rather limited and solely based on ∆tCa data, which has been previously demonstrated to be associated with aortic calcification [11] and statistically related to the increased risk for MI with the use of higher dCa among diabetic HD patients [12]. Although it was previously suggested that serum Ca kinetics assessment should be rather based on iCa than on tCa levels [18], only a single CaMB study previously reported both total and ionized serum calcium data when evaluating the use of 2.7 mEq/L dCa [19]. In the present study, we measured both serum ionized and total ∆Ca and found that the presence of CHD is independently associated with ∆iCa but not ∆tCa. In contrast to the study by Tagawa et al. [12], our results are not limited to diabetic HD patients and based on an established dialysate-sided iCaMB assessment approach [4, 20]. Unlike previous CaMB studies, we measured both pre- and postfilter ionized dCa in order to account for known differences between prescribed label value prefilter dCa versus measured prefilter dCa [4]. Indeed, the mean ionized prefilter dCa was found to be 2.92 (±0.12) mEq/L, that is, 83% of 3.5 mEq/L dCa label value (Table 1). As the dialysate-to-blood iCa gradient is the main driving force of iCaMB, the use of label value prefilter dCa for iCaMB calculation would result in a substantial overestimation of the latter. We, thus, accounted for the observed discrepancy between label value prefilter dCa and measured prefilter dCa that most likely results from Ca complexation with anions, for example, bicarbonate, in the dialysate fluid. Interestingly, no relevant differences between ionized and total dCa-based CaMB assessment have been observed in previous studies [5, 21]. Without assessing CaMB, however, previous studies focussing on ∆Ca could not answer the question whether high intradialytic Ca load or altered extracellular Ca buffering per se affects CV risk. Importantly, iCaMB did not correlate with intradialytic Ca buffer capacity, ∆iCa, ∆tCa, or ∆iCaECF in our study, suggesting that intradialytic Ca burden – an established CV risk factor [22] – and Ca buffering are distinct processes that should be considered as independent CV risk factors. Based on a 2.2-fold increase in serum tCa compared with serum iCa, we approximatively calculated that 30% of buffered iCaMB is actually protein-bound; however, the remaining 70% are buffered outside the ECF compartment. In this regard, the involvement of osteocalcin, the most abundant noncollagenous bone protein [23, 24], in the rapidly accessible and exchangeable Ca pool has been suggested previously [25] and first experimentally described in vivo by our workgroup [4]. In this regard, osteocalcin has been supposed to stabilize bone surface brushite layer, which might act as rapidly accessible and exchangeable Ca pool [10]. By confirming the positive independent association of ucOC, that is, a circulating form of OC that does not undergo carboxylation or is actively decarboxylated [26, 27], with intradialytic Ca buffer capacity in the present study, we further support the hypothesis of bone being – at least partly – involved in acute extracellular Ca regulation [10]. It is tempting to speculate that extraosseous Ca storage, for example, in the vasculature and soft tissues, acts synergistically in terms of intradialytic Ca buffering, but – in contrast to rapid osseous Ca storage – rather increases CV risk. In this context, we found a negative correlation of cOC levels with PWV, an established marker of vascular calcification and CV disease in HD patients [28], thereby indicating a potential shift from bone to vascular Ca storage: cOC is synthesized by mature osteoblasts, stimulated by active vitamin D, and constitutively present in bone matrix; however, after carboxylation, a proportion can effuse in the circulation [24]. In this regard, we found a positive correlation of active vitamin D dose with ucOC (p < 0.01) and cOC (ns) in our study cohort. Assuming its role as a brushite stabilizer, reduced cOC levels might indicate reduced bone Ca buffer capacity and higher CV risk due to compensatory extraosseous Ca deposition, which is reflected by increased PWV (but not necessarily by reduced overall Ca buffer capacity). This is in line with several in vivo studies demonstrating an inverse relationship between circulating osteocalcin levels and CV risk both in the general [29-33] and HD population, where osteocalcin was suggested to have a vasculoprotective role [34-37]. While a few studies found contradictory results [38-40], recent evidence even suggested a U-shaped relationship [41, 42]. Animal studies also revealed conflicting results, with some demonstrating an inhibitory effect of osteocalcin on bone formation [43] and others identifying osteocalcin as potential promotor of osteochondrogenic differentiation of vascular smooth muscle cells during vitamin D-induced vascular calcification [44]. Growing evidence now suggests a dual role of osteocalcin, with different circulating osteocalcin isoforms that regulate vascular calcification and act as endocrine mediator, for example, by regulating insulin secretion and sensitivity, as well as tissue-expressed osteocalcin directly affecting the calcification process [27, 37]. Likewise, our study provides in vivo evidence for a predictive potential of different circulating osteocalcin isoforms in HD patients, that is, ucOC for intradialytic Ca buffering and cOC for a shift from bone to vascular calcium storage. Based on the strong independent correlation of aortic PWV with age, a much larger HD cohort would be statistically required to elucidate an independent association of cOC with elevated PWV. Our present study validates previous study findings from smaller HD cohorts [4] in an independent larger patient population and presents novel evidence for the clinical relevance of assessing intradialytic calcium kinetics by showing an independent association of ∆iCa with CHD. Our data also support the potential role of OC as a marker of altered intradialytic Ca buffer capacity and vascular calcification thereby reflecting the complex interplay between bone- and vascular Ca storage in HD patients. Overall, our results strengthen the concept of a rapidly accessible and exchangeable Ca pool involved in intradialytic Ca regulation in vivo. Cross-sectional study design and the small patient number are main study limitations. Thus, prospective studies in larger HD patient cohorts are needed to evaluate the prognostic potential of intradialytic Ca kinetics parameters.
Acknowledgement
We thank the patients for study participation, the dialysis staff for supporting the study, and Dr. Herbert Schramek for his valuable contribution during manuscript preparation.
Statement of Ethics
This study was conducted in accordance with the World Medical Association Declaration of Helsinki. Written informed consent was obtained from each study participant prior to study inclusion using a standardized patient information and consensus form according to Good Clinical Practice (GCP) guidelines, and the study protocol was approved by the Innsbruck Medical University ethics committee prior to study initiation (protocol number AN2014-0313 343/4.7 391/5.3). All patient-associated information was managed entirely coded. Data collection was conducted using a standardized evaluation form (case report form) according to GCP recommendations. All patient-associated samples and clinical data are subject to privacy protection according to the current European General Data Protection Regulation. This study complies with the STROBE standards for reporting of observational studies.
Conflict of Interest Statement
The authors have no conflicts of interest to declare.
Funding Sources
The Austrian National Bank supported this work by providing an educational research grant (OENB Jubiläumsfonds Project No. 15366) to M.P. The funders did not have any role in study design; collection, analysis, and interpretation of data; writing the report; and the decision to submit the report for publication.
Author Contributions
M.P.: research idea and study design; L.F. and R.H.: data acquisition; M.P., L.F., R.H., and T.R.: data analysis/interpretation; M.P. and T.R.: statistical analysis; M.P. and G.M.: manuscript preparation and drafting; M.P., L.F., R.H., T.R., and G.M: approval of the final version.