- 1Division of Cardiology, Department of Internal Medicine, National Taiwan University Hospital Hsinchu Branch, Hsinchu, Taiwan
- 2College of Medicine, National Taiwan University, Taipei, Taiwan
- 3Department of Biological Science and Technology, National Yang Ming Chiao Tung University, Hsinchu, Taiwan
- 4ANSN Clinic, Hsinchu, Taiwan
- 5Hemodialysis Center, National Taiwan University Hospital Hsinchu Branch, Hsinchu, Taiwan
- 6School of Medicine, National Yang Ming Chiao Tung University, Hsinchu, Taiwan
- 7School of Medicine, National Yang-Ming University, Taipei, Taiwan
- 8Division of Nephrology, Department of Medicine, Taipei Veterans General Hospital, Taipei, Taiwan
- 9Institute of Population Health Science, National Health Research Institutes, Miaoli County, Taiwan
- 10Center of Quality Management, National Taiwan University Hospital Hsinchu Branch, Hsinchu, Taiwan
- 11Institute of Biomedical Engineering, National Tsing Hua University, Hsinchu, Taiwan
- 12Institute of Cellular and System Medicine, National Health Research Institutes, Miaoli County, Taiwan
Background: Blood pressure variability (BPV) is an important risk factor for cardiovascular events in hemodialysis patients. We sought to determine the impact of BPV on hemodialysis access thrombosis.
Methods: We enrolled 1,011 prevalent hemodialysis patients from 12 hemodialysis centers since January 2018 and followed them until December 2020. Predialysis blood pressure (BP) was assessed at 12-week intervals. The coefficient of variation derived from 36 consecutive BP measurements was used as the metric for variability. The primary outcome was incident hemodialysis access thrombosis. Linear regression models were used to assess factors associated with BPV at baseline. Kaplan-Meier curves of the time until vascular access events were drawn and log-rank tests were calculated. Cox proportional hazards models were performed to assess the association of BPV with incident vascular access events.
Results: The average coefficient of variance for systolic BPV was 10.9%. BPV was associated with age, body mass index, mean BP, diabetes, coronary and peripheral artery disease, history of access dysfunction, graft access, intradialytic hypotension, and use of antihypertensive medications. There were 194 access thrombosis events and 451 access stenosis events during a median follow-up period of 30 months. After adjustment of potential confounding factors, BPV was associated with increased risk of access thrombosis [hazard ratio = 1.27, 95% confidence interval (CI), 1.18–1.44, per 1 standard deviation increase in BPV]. The patients in the highest BPV quartile had 2.45 times the risk of thrombosis (CI, 1.62–3.70). The association was independent of average BP, intradialytic hypotension, and comorbidities. Similar trends of association were found in the subgroups analyzed. Comparative analysis using a time-varying variable model and different metrics of BPV showed consistent results.
Conclusion: Our findings underscored the impact of BP fluctuation on vascular access thrombosis.
Introduction
Hemodialysis access dysfunction continues to be a major source of morbidity and mortality in patients with end-stage renal disease (1). Thrombosis is the most common cause of access failure and also causes substantial stress on hemodialysis patients and their caregivers. Although stenosis is the most common cause of access thrombosis, 20–40% of thrombotic events occur in the absence of stenosis (2, 3). Moreover, various factors may precipitate the occurrence of thrombosis in an access with underlying stenosis (4). Factors other than stenosis are less well studied, but identifying these factors is critical for the prevention of dialysis access thrombosis.
Hypertension is a well-known risk factor for vascular diseases in the general population. Nonetheless, most studies found the association of vascular access thrombosis with hypotension (5–7). Blood pressure (BP) is the primary driving force of blood flow in the vascular access, and a low-flow state is one of the major determinants for thrombosis (8). Although this causal relationship makes sense, studies have rarely examined the relationship between BP and access thrombosis in a comprehensive manner. Most previous studies used average BP to investigate the relationship between BP and thrombosis. However, flow-related thrombosis should be susceptible to BP fluctuation as well, not only the usual BP. Hemodialysis patients are particularly prone to wide BP fluctuations, either due to comorbidities or unique dialysis factors. Although average BP is traditionally a risk factor for systemic vascular events, growing evidence suggests that BP instability and variability may contribute to the development of vascular events as well (9–11).
Both short-term (intradialytic or interdialytic period) or long-term (day-by-day or visit-by-visit) BP variability (BPV) are informative measures to predict vascular events (12). The effect of short-term BPV on the prognosis of hemodialysis patients has been reported. Intradialytic hypotension and intradialytic BP variability were associated with vascular access outcomes (6, 13, 14). Long-term BPV is reproducible over time but is only modestly related to short-term variability (15, 16). Growing evidence shows that long-term BPV predicts mortality and cardiovascular events in hemodialysis patients (9–11, 17–21). Nonetheless, the impact of long-term BPV on dialysis access remained unknown. Long-term BPV measurements are more feasible for hemodialysis patients than non-dialysis patients, which could be derived from BPs recorded during dialysis sessions.
We hypothesized that higher visit-to-visit BPV would increase the risk of vascular access thrombosis. Therefore, we conducted a prospective study to investigate the relation between BPV and dialysis access thrombosis.
Materials and methods
Study participants and design
This is a prospective multicenter study (Hsinchu V.A. study, ClinicalTrials.gov identifier: NCT04692636) designed to investigate clinical factors related to cardiovascular events of maintenance hemodialysis patients. The study participants were recruited from January 1, 2018 to December 31, 2018, from 12 hemodialysis centers in the Hsinchu district of Taiwan. The Hsinchu district is composed of 1 city and 13 townships, with a total population of about 1,000,000. Four of the 12 hemodialysis centers are hospital-based centers and eight of them are dialysis clinics. The inclusion criteria were as follows: (1) age between 18 and 90 years old, (2) maintenance dialysis of more than 6 months, and (3) no hospitalization during the previous 3 months. For the analysis of hemodialysis access outcomes in this study, the following exclusion criteria were used during the exposure period: (1) use of central venous catheters in dialysis, (2) death, facility transfer, or modality change, (3) missing dialysis session, and (4) missing BP data > 10% (Figure 1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/881454/fcvm-09-881454-HTML/image_m/fcvm-09-881454-g001.jpg)
Figure 1. Study design (A) and flow diagram (B). Flow diagram of study participants and schematic presentation of outcomes models for primary analysis (baseline BPV model) and comparative analysis (time-varying BPV model). BP, blood pressure; BPV, blood pressure variability; CVC, central vein catheter; HD, hemodialysis; and lab, laboratory data.
We obtained information on patients’ demographics, comorbidities, dialysis-related information, medications, and laboratory data from the medical records of dialysis centers. These data were collected at baseline and updated every 3 months by trained study coordinators. BP data were extracted from a computerized dialysis information system (INDAS, Integrative Data Acquisition System, TW) in which vital signs and dialysis parameters from each dialysis session are stored digitally. Dialysis accesses were followed at all respective hemodialysis centers under the same protocol, and were referred to our angiographic unit for evaluation or management under the same criteria. Events of dialysis access were collected by reviewing the dialysis records, and angiography or surgery reports at a 3-month interval. Subjects were followed until December 31, 2020, and censored at the time of death, kidney transplant, or transfer to a non-study center. The study was approved by the Institutional Review Board and all patients signed an informed consent to participate in the study.
Exposure variables
Blood pressure measurements were conducted before (pre-dialysis), during, and after (post-dialysis) each hemodialysis session, in a seated position, using automated oscillometric devices as per the dialysis unit routine. BP values were extracted from a computerized database every 90 days. We calculated the baseline mean BP and BPV using all the pre-dialysis BP measurements from week 1 to the end of week 12. During the follow-up period, we assessed the mean BP and BPV every 12 weeks. Allowable ranges were 50–300 mm Hg for systolic BP and 10 to 150 mm Hg for diastolic BP. When a value was outside the range, the data point was coded as missing. We only analyzed patients with BP data of more than 30 occasions to maintain excellent reproducibility of BPV data (Supplementary Figure 1). The BPV was calculated over a 12-week interval as the exposure for the next 12-week period at risk (Figure 1). The raw BP data were transformed into four BPV metrics: standard deviation (SD), coefficient of variance (CV), absolute real variation (ARV), and variance independent of the mean (VIM). ARV is calculated as the average absolute difference between consecutive measurements. VIM is a transformation of SD and is by definition uncorrelated to mean BP. VIM is calculated by fitting a curve through a plot of SD systolic BP (y-axis) against mean systolic BP (x-axis) with the parameter d estimated from the curve [VIM = (SD/meanx)] (22). We used the CV of consecutive 36 predialysis systolic BPs in the 12-week exposure period as the primary measure of BPV for analysis. We also used SD, VIM, and ARV of the predialysis BPs in the exposure period of the additional analysis.
Covariates
The following domains of covariates were assessed over the first 12 weeks as the baseline level and updated every 12 weeks: demographic and anthropometric, clinical, dialysis (dialysis vintage, vascular access age, and vascular access type), laboratory, and medication. Demographic and anthropometric domains include age, sex, and body mass index. Clinical domains included diabetes, coronary artery disease, cerebrovascular disease, heart failure, and history of vascular access dysfunction. To account for other comorbidities, a comorbidity index previously validated in dialysis patients was calculated (23). Dialysis domains included dialysis vintage, shunt age, vascular access type (fistula vs. graft), urea clearance, relative fluid removed per dialysis session [(pre-dialysis weight-post-dialysis weight)/post-dialysis weight], and dialysis institutions (dialysis clinics or hospital-based dialysis center). Laboratory domains included serum albumin, cholesterol, hemoglobin, calcium, and phosphate product were collected every 12 weeks. Medication domains included antiplatelet, anticoagulant, statin, antihypertensive medicines. We also included mean predialysis SBP and intradialytic hypotension as covariates. Intradialytic hypotension was defined as a nadir SBP less than 90 mm Hg in at least 30% of dialysis sessions (24).
Outcome measurements
The surveillance protocol of dialysis access included physical examination, dynamic venous pressure monitoring at each dialysis session and dialysis dose, transonic examination of access flow (if available) at monthly interval. Patients were referred to our angiographic unit for evaluation on the basis of the following criteria: (1) clinical manifestations suggesting vascular access dysfunction (decreased thrill, abnormal bruit, increased pulsatility, and prolonged bleeding from the puncture site); (2) a reduction in flow rate of >25% from baseline; (3) total access blood flow <500 ml/min by ultrasound dilution method; and (4) increased venous pressure during dialysis (dynamic venous pressure exceeding threshold level measured three consecutive times); and (5) unexplained decrease (>0.2 unit) in the delivered dialysis dose. Our primary outcome of interest was the time to the first vascular access thrombosis, defined as the access that had clotted without blood flow. Information of vascular access thrombosis was obtained from documentation in the dialysis records, angiography reports, or surgery reports. The secondary outcome was the time to the first vascular access stenosis, defined as the anatomical stenosis of more than 50% by angiography with corresponding clinical or hemodynamic abnormalities of the referring criteria. Because the patient may change vascular access throughout the follow-up period, the only outcome of first vascular access was used for analysis. In the primary baseline model analysis, vascular access outcomes were defined as the first event that occurred from week 13 to the end of the study or censoring. In the companion timing-varying analyses, vascular outcomes were defined as the first event within 12 weeks after the end of the 12-week exposure period (Figure 1).
Statistical analysis
We reported patients’ baseline characteristics across quantiles of systolic BPV. We assessed factors associated with BPV at baseline using linear regression models. Kaplan-Meier curves of the time until vascular access events were drawn and log-rank tests were calculated. Cox proportional hazards models were performed to assess the association of BPV with incident vascular access events. Three sets of analyses were examined to access the relationship between BPV and vascular access events. The primary analysis describes the relationship between baseline BPV (derived from a 12-week exposure period) and vascular access outcomes. The first sensitivity analysis used time-varying BPV (calculated over a 12-week interval) as well as other time-varying covariates for analysis. The second sensitivity analysis excluded thrombosis or stenosis events during hospitalization in order to eliminate the hemodynamic effect due to acute illness. Univariate analyses were done before proceeding to multivariate models. Unadjusted and adjusted models with hazard ratios (HRs) and 95% CI were presented per 1-SD increase of BPV. Age, sex, and mean SBP were considered as potential a priori confounders. Significance thresholds of 0.1 in the univariate linear regression model were selected as potential confounding factors in multivariate models. Proportional hazard assumption was checked graphically using the log-log plot and was found to be acceptable for the factors of interest. Subgroup analyses were performed on patients defined by age, sex, systolic BP, diabetes, coronary artery disease, intradialytic hypotension, vascular access, and institution. These interaction of these factors on the associations between BPV and thrombosis events were assessed by adding interaction terms in the model. We considered two-tailed P values < 0.05 as statistically significant. All analyses were performed using the SPSS software (IBM SPSS Statistics for Windows, Version 22.0. Armonk, NY, United States: IBM Corp).
Results
Baseline characteristics
A total of 1,137 patients were enrolled. After the exclusion of 77 patients dialyzed via central vein catheters and 49 patients for other reasons, the final cohort included for analysis consisted of 1,011 patients (Figure 1). Baseline characteristics of the study population are presented in Table 1. The mean age of the study participants was 66 years (SD of 14 years) and 521 of them were male (52%). Their median hemodialysis duration was 59 months, median shunt age was 56 months, and 192 patients (19%) using arteriovenous graft (AVG) as their regular vascular access at enrollment. Baseline characteristics of the final cohort were similar to the hemodialysis patients of a nationwide survey in 2017 (Supplementary Table 1). These patients were followed up for a median of 30 months and 243 patients died during the follow-up period. Five patients received a kidney transplant and 107 patients were lost to follow-up or transferred to non-study centers. A total of 451 patients had vascular access stenosis events and 194 patients had vascular access thrombosis events.
Blood pressure and blood pressure variability metrics
The mean pre-dialysis systolic BP was 141 mm Hg and diastolic BP was 71 mm Hg. The average CV for systolic BP was 10.9% (SD, 3.0%). Other metrics of BPV, such as SD, ARV, and VIM, were provided in Supplementary Table 2. Table 2 shows that baseline BPV was correlated with age, body mass index, systolic BP, diabetes, coronary artery disease, peripheral artery disease, vascular access dysfunction, comorbidity index, AVG, intradialytic hypotension, albumin level, cholesterol level, and use of antihypertensive medications.
Effect of systolic blood pressure and blood pressure variability on outcomes of vascular access
Baseline mean BP was not associated with stenosis or thrombosis of dialysis access. When the mean BP was stratified into quartiles, the lowest BP group (quartile 1) had a higher risk of access thrombosis (Table 3). The BP quartiles were not associated with vascular access stenosis. Baseline BPV was associated with vascular access stenosis (Figure 2). In univariate Cox regression analysis, per 1 SD increase of BPV was associated with a 17% increased risk of stenosis [HR, 1.17; confidence interval (CI), 1.08–1.28, P < 0.001]. When BPV was stratified into quartiles, the highest BPV quartile had 1.60 times (CI, 1.22–2.08, P < 0.001) risk of stenosis than the lowest quartile (Table 3). After multivariate adjustment, BPV did not associate with risk of stenosis (Table 4). Baseline BPV was associated with thrombosis of vascular access (Figure 2). In univariate Cox regression analysis, every 1 SD increase in BPV was associated with a 35.7% increased risk of thrombosis (HR, 1.36; CI, 1.21–1.52, P < 0.001). When BPV was stratified into quartiles, the highest BPV quintile had 2.45 times the risk of thrombosis (1.62–3.70, P < 0.001) than the lowest quartile, with a graded effect among the quartiles (Table 3). After multivariate adjustment, BPV remained significantly associated with the risk of thrombosis (Table 4). The relationship between BP and BPV quartiles with vascular access outcomes is displayed in Figure 3.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/881454/fcvm-09-881454-HTML/image_m/fcvm-09-881454-t003.jpg)
Table 3. Univariate Cox regression analysis of predictors for incident vascular access thrombosis or stenosis during follow-up period.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/881454/fcvm-09-881454-HTML/image_m/fcvm-09-881454-g002.jpg)
Figure 2. Kaplan-Meier plots of vascular access outcomes. Kaplan-Meier plots of stenosis-free patency (left panel) and thrombosis-free patency (right panel) rates by quartiles of predialysis systolic blood pressure variability (SBPV). Q1 (blue line), first (lowest) SBPV quartile; Q2 (red line), second quartile; Q3 (green line), third quartile; Q4 (yellow line), fourth (highest) SBPV quartile. No. at risk, the number of patients at risk.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/881454/fcvm-09-881454-HTML/image_m/fcvm-09-881454-t004.jpg)
Table 4. Multivariate Cox regression analysis of the association between BPV and outcomes of vascular access, displayed by the whole cohort and stratified by access types.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/881454/fcvm-09-881454-HTML/image_m/fcvm-09-881454-g003.jpg)
Figure 3. Risk of vascular access outcomes by systolic blood pressure and variability. Unadjusted and adjusted hazard ratio of thrombosis (A,C) and stenosis (B,D) events by quantiles of mean systolic blood pressure in panels (A,B) and systolic blood pressure variability in panels (C,D). The risk was adjusted for age, sex, and all the baseline factors with a p-value < 0.1 in the univariate Cox regression analysis. The squares indicate the hazard ratio and the error bars indicate 95% confidence intervals. CI, confidence interval; HR, hazard ratio; SBP, systolic blood pressure; and SBPV, systolic blood pressure variability.
Subgroup analyses
Because the risk of events differed significantly between arteriovenous fistula (AVF) and AVG, we analyzed the association of BPV and thrombosis by access types (Table 4). After adjustment, BPV remains associated with thrombosis both in patients with AVF and AVG. The results in the univariable model were stratified by potential confounders, shown in Figure 4. Subgroup analysis showed that the association between BPV and thrombosis remained significant in most subgroups, except for patients with history of vascular access dysfunction.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/881454/fcvm-09-881454-HTML/image_m/fcvm-09-881454-g004.jpg)
Figure 4. Subgroup analyses. Unadjusted hazard ratio of thrombosis events per one standard deviation increase of predialysis systolic blood pressure variability (BPV) in different subgroups of patients. The interaction of factors in each subgroup on the association of BPV with thrombosis were presented by the “p for interaction.”
Sensitivity analysis
To account for the change of BPV over time, we analyzed the association of BPV and dialysis access outcomes by treating BPV as a time-dependent variable, rolling over 12 weeks. The effect of BPV on vascular access outcomes, regardless of stenosis or thrombosis, was similar to that observed in the baseline BPV model (Supplementary Table 3). To account for the effect of mean BP, different metrics of BPV, such as VIM and ARV, were used in the analysis. Analysis of these BP-independent indices with vascular access outcomes showed similar results, both for stenosis events and thrombosis events (Supplementary Table 4). After excluding 5 stenosis events and 4 thrombosis events during hospitalization, the effect of systolic BPV on vascular outcomes was similar to that in the primary analysis (Supplementary Table 5).
Discussion
Visit-to-visit BPV is emerging as a relevant risk factor for cardiovascular events in dialysis and non-dialysis patients (9, 11, 25, 26). We showed in this cohort that the average predialysis systolic BPV was 10.9%, high than 6.1% in the general populations and similar to 9.9–10.7% reported in hemodialysis patients (11, 27, 28). Our analysis shows for the first time that visit-to-visit BPV is significantly associated with dialysis access thrombosis, either fistulas or grafts. This association is independent of baseline BP levels or intradialytic hypotension. After the adjustment of several confounding factors, BPV remained a significant predictor of thrombosis. Each SD increase in BPV was associated with a 26.8% increased risk of thrombosis, and the highest BPV quartile (CV > 14.9%) had a 2.45 times higher risk of thrombosis. Thus, visit-to-visit BPV seems to be a potential predictor and a probable modifiable risk factor for thrombosis of vascular access.
A low-flow state is a critical determinant for access thrombosis. Arterial BP is the primary driving force of access blood flow. Low BP has been shown to decrease access flow and increase thrombosis events. In the absence of structural abnormalities, low BP accounted for at least 20–40% access thrombosis (5). In previous studies, mean BP was used to evaluate the effect of BP on the risk of thrombosis, with an increase from 3% to 11% per SD decrease of BP. Nonetheless, greater fluctuations and lower nadir of BP theoretically increase the exposure to low-flow state. We showed that the highest BPV quartile had a 2.45 times higher risk of thrombosis than the lowest quartile. In the subgroup analysis, the effect was more prominent among individuals with systolic BP < 140 mm Hg, suggesting a pathogenic link to low BP burden. Although BP and BPV were closely related, the effect of BPV remained even after meticulous adjustment by different methods. Taking the above evidence together, BP fluctuation is a relevant risk factor for vascular access thrombosis.
Previous studies of the relationship between BP and vascular access focused on short-term variability or instability during dialysis. For example, a secondary analysis of the Hemodialysis (HEMO) Study revealed frequent episodes of intradialytic hypotension associated with AV fistula thrombosis. Two recent retrospective studies also found that the variability of intradialytic BP can predict access thrombosis (13, 14). Nonetheless, in most circumstances, access thrombosis occurred in the inter-dialytic period. The multivariate and subgroup analysis in our analyses suggest that the predictive value of long-term BPV was independent of intradialytic hypotension. In addition to arterial stiffness, long-term variability was more affected by environmental and behavioral factors. In contrast, short-term variability during dialysis was susceptible to cardiac function, arterial stiffness, and dialysis factors (6, 29). Previous studies demonstrated that a greater interdialytic weight gain or ultrafiltration volume was associated with short-term variability during dialysis. In our study, long-term BPV was not affected by fluid removal or urea clearance, as was also demonstrated in previous studies (9, 11). Therefore, the predictive value of long-term BPV was independent and complementary to intradialytic variability (27, 30).
There were other possible reasons for the association between BPV and thrombosis. Baseline BPV was associated with risk factors of arterial stiffness, such as aging, obesity, diabetes, and vascular diseases (31). Diabetes and vascular diseases also increased the risk of access thrombosis. BPV may only be an intermediate factor in the thrombosis caused by diabetes and vascular diseases. Nonetheless, BPV predicted thrombosis after adjusting for diabetes and vascular diseases. The subgroup analysis also demonstrated that the association was independent of diabetes or vascular diseases. Therefore, BPV may not only be a marker of underlying vascular diseases but also a maker of thrombosis events. Dialysis access types are associated with BPV, as demonstrated in our results and previous studies (32). We analyzed the effect of BPV by access types because the anatomical configuration and hence hemodynamics of fistulas and grafts differ greatly. Fistulas require only one surgical anastomosis and are almost entirely endothelialized that do not often thrombose. Grafts are composed of unendothelialized foreign material that are inherently thrombogenic and had five-times risk of thrombosis in this study. Nonetheless, even after stratification by access types, high BPV remains a significant risk factor of thrombosis, either for fistulas or grafts.
Blood pressure variability was associated with stenosis of vascular access in univariate analysis. Intimal hyperplasia at the outflow vein is the most common cause of stenosis for either fistulas or grafts (33). BPV was reported to be associated with inflammation and endothelial dysfunction (31, 34), that may be the links to the development of intimal hyperplasia (35, 36). Arterial stiffness is one of the mechanisms for visit-to-visit BPV (37, 38). We indeed found that BPV was associated with risk factors of arterial stiffness, such as age, diabetes, and vascular diseases. Therefore, the association between BPV and stenosis may be mediated via these common pathogenic pathways as well (35, 36, 39). After multivariate adjustment, the risk estimates and significance of the association between BPV and stenosis decreased. Therefore, BPV is more likely to be a risk marker of stenosis rather than a risk maker.
Currently, no standardized method is available for assessing long-term BPV, either in dialysis or non-dialysis patients. The density of BP measurements varied widely among previous studies. In our study, high-density measurement (36 times in 12 weeks) was used to estimate BPV based on data from a digitized dialysis information system. Although such a high-density measurement is not easy for non-dialysis patients, it can be obtained from regular dialysis sessions without additional loadings on the staff. Home BP is addressed in recent guidelines for BPV assessment. Nonetheless, device standardization, measurement, timing, and data transformation are all potential obstacles to the application of home BP in clinical practice. In contrast, using clinic BP in dialysis centers is more feasible and the prognostic value was well-documented in previous studies (9–11). We used CV, which also accounted for mean BP as the primary BPV metric for analysis. Other BPV metrics, such as VIM and ARV, are available to reflect the fluctuation independent of average BP. Nonetheless, these metrics require complex calculations, which limits their clinical applicability.
Our study has several strengths. It was conducted in a typical clinical setting using a representative community-based cohort, with characteristics comparable to that of a nationwide registry. The detailed time-varying patient-level data on key variables were available, including fluid removal and anti-hypertensive medications, detailed ascertainment of comorbidities using medical records, and prospective ascertainment of vascular access events using dialysis records or angiographic reports. It is a multi-center study with a longitudinal study design and a long-term follow-up period, allowing us to collect repeated BP measurements, relevant confounders, and reliable outcomes. Unlike most studies that focused on a certain BPV indicator, we used several BPV indices to strengthen our findings. We used high-density measurements, which rendered a stable and reproducible BPV.
Some limitations should be addressed. Our sample was relatively small compared to nationwide registries or claims databases. From the research point of view, VIM is an ideal approach to exclude the influence of usual BP. However, in terms of clinical decision making, a CV may be a more practical option. Although the use of medications was adjusted, we did not have a direct assessment of medication adherence. BPV measurements were performed as part of routine clinical practice and we did not standardize the device for BP measurement. Different devices used for BP measurement could be a source of extra variability. The validity for lower-density BP measurements was not clear. A high-density BPV measurement is laborious and the assistance of information technology is important in clinical practice. We did not examine variability of diastolic BP or mean BP because studies on diastolic BPV are relatively few and have conflicting results (12, 28, 40). Finally, the relative few vascular thrombosis events may have contributed to the lack of a significant findings among the subgroups.
Given the high cost and prevalence of vascular access thrombosis, finding modifiable risk factors and implementable strategies may improve patient outcomes and healthcare burden. The current practice already dictates the association of BP variability and instability with mortality and cardiovascular events. Our results expanded the role of BPV in vascular access thrombosis events, either in fistulas or grafts. The impact of BPV is more pronounced in dialysis patients with normal BP. Uncertainty regarding appropriate BP targets in hemodialysis patients continues, and whether higher BP targets are warranted in patients at risk of access thrombosis deserved further investigation. Evidence derived from post hoc analyses of previous studies suggested that certain classes of antihypertensive medicines were more effective on BPV than others (41). Definite answers can only be provided by randomized clinical trials. Our analysis underscores the importance of including vascular access patency as a relevant outcome in future studies on BP management of hemodialysis patients. The importance of BP fluctuation on vascular access thrombosis should be emphasized, rather than focusing on average BP alone. Further studies are needed to assess whether maintaining the stability of BP can be a promising target in the prevention of vascular access thrombosis.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving human participants were reviewed and approved by National University Hospital, Hsinchu Branch, IRB. The patients/participants provided their written informed consent to participate in this study.
Author contributions
M-YH and C-CW: research idea and study design. M-YH, S-YC, Chi-HC, and Chiu-HC: data acquisition. M-TL, S-YC, and C-CW: data analysis/interpretation. M-TL and S-YC: statistical analysis. C-CW: supervision or mentorship. All author contributed important intellectual content during manuscript drafting or revision and agrees to be personally accountable for the individual’s own contributions and to ensure that questions pertaining to the accuracy or integrity of any portion of the work, even one in which the author was not directly involved, are appropriately investigated and resolved, including with documentation in the literature if appropriate.
Funding
This study was supported by grants from National Taiwan University Hospital Hsinchu Branch (107-HCH025, 108-HCH046, 109-HCH001, and 109-HCH057) and the Ministry of Science and Technology (MOST-106-2314-B-002-173-MY3 and MOST-109-2314-B-002-241-MY3).
Acknowledgments
We thank the staff of the Biotechnology R&D Center of National Taiwan University Hospital Hsinchu Branch for their assistance in the statistical analysis. We thank the staff of ANSN clinics, Chiahua clinic, Catholic Mercy Hospital, Nanmen General Hospital, and National Taiwan University Hospital Hsinchu and Chutong branches for their assistance in enrolling and following the patients. We would like to thank Ann Charis Tan for her assistance in reviewing and revising the manuscript for grammar and syntax.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcvm.2022.881454/full#supplementary-material
References
1. Saran R, Robinson B, Abbott KC, Bragg-Gresham J, Chen X, Gipson D, et al. US renal data system 2019 annual data report: epidemiology of kidney disease in the United States. Am J Kidney Dis. (2020) 75(1 Suppl. 1):A6–7. doi: 10.1053/j.ajkd.2019.09.003
2. Moist LM, Churchill DN, House AA, Millward SF, Elliott JE, Kribs SW, et al. Regular monitoring of access flow compared with monitoring of venous pressure fails to improve graft survival. J Am Soc Nephrol. (2003) 14:2645–53. doi: 10.1097/01.ASN.0000089562.98338.60
3. Schwab SJ, Oliver MJ, Suhocki P, McCann R. Hemodialysis arteriovenous access: detection of stenosis and response to treatment by vascular access blood flow. Kidney Int. (2001) 59:358–62. doi: 10.1046/j.1523-1755.2001.00498.x
4. Wu CC, Hsieh MY, Hung SC, Kuo KL, Tsai TH, Lai CL, et al. Serum indoxyl sulfate associates with postangioplasty thrombosis of dialysis grafts. J Am Soc Nephrol. (2016) 27:1254–64. doi: 10.1681/ASN.2015010068
5. May RE, Himmelfarb J, Yenicesu M, Knights S, Ikizler TA, Schulman G, et al. Predictive measures of vascular access thrombosis: a prospective study. Kidney Int. (1997) 52:1656–62. doi: 10.1038/ki.1997.499
6. Chang TI, Paik J, Greene T, Desai M, Bech F, Cheung AK, et al. Intradialytic hypotension and vascular access thrombosis. J Am Soc Nephrol. (2011) 22:1526–33. doi: 10.1681/ASN.2010101119
7. Goldwasser P, Avram MM, Collier JT, Michel MA, Gusik SA, Mittman N. Correlates of vascular access occlusion in hemodialysis. Am J Kidney Dis. (1994) 24:785–94. doi: 10.1016/S0272-6386(12)80672-X
9. Shafi T, Sozio SM, Bandeen-Roche KJ, Ephraim PL, Luly JR, St Peter WL, et al. Predialysis systolic BP variability and outcomes in hemodialysis patients. J Am Soc Nephrol. (2014) 25:799–809. doi: 10.1681/ASN.2013060667
10. Brunelli SM, Thadhani RI, Lynch KE, Ankers ED, Joffe MM, Boston R, et al. Association between long-term blood pressure variability and mortality among incident hemodialysis patients. Am J Kidney Dis. (2008) 52:716–26. doi: 10.1053/j.ajkd.2008.04.032
11. Chang TI, Flythe JE, Brunelli SM, Muntner P, Greene T, Cheung AK, et al. Visit-to-visit systolic blood pressure variability and outcomes in hemodialysis. J Hum Hypertens. (2014) 28:18–24. doi: 10.1038/jhh.2013.49
12. Rothwell PM, Howard SC, Dolan E, O’Brien E, Dobson JE, Dahlöf B, et al. Prognostic significance of visit-to-visit variability, maximum systolic blood pressure, and episodic hypertension. Lancet. (2010) 375:895–905. doi: 10.1016/S0140-6736(10)60308-X
13. Zheng CM, Chiu YP, Hou YC, Liu YM, Wu MS, Lin YF, et al. Influence of intradialytic systolic blood pressure changes on arteriovenous access thrombosis in maintenance hemodialysis patients. Int J Clin Pract. (2020) 75:e13799. doi: 10.1111/ijcp.13799
14. Kim JY, Seo HM, Kim M, Kim H. A relationship of intradialytic blood pressure variability with vascular access outcomes in patients on hemodialysis. Hemodial Int. (2019) 23:158–66. doi: 10.1111/hdi.12720
15. Brueren MM, Limpt PV, Schouten HJA, de Leeuw PW, van Ree JW. Is a series of blood pressure measurements by the general practitioner or the patient a reliable alternative to ambulatory blood pressure measurement? A study in general practice with reference to short-term and long-term between-visit variability. Am J Hypertens. (1997) 10:879–85. doi: 10.1016/S0895-7061(97)00125-8
16. Howard SC, Rothwell PM. Reproducibility of measures of visit-to-visit variability in blood pressure after transient ischaemic attack or minor stroke. Cerebrovasc Dis. (2009) 28:331–40. doi: 10.1159/000229551
17. Liao R, Li J, Xiong Y, Lin L, Wang L, Sun S, et al. Association of peridialysis blood pressure and its variability with cardiovascular events in hemodialysis patients. Kidney Blood Press Res. (2018) 43:1352–62. doi: 10.1159/000492595
18. Di Iorio B, Di Micco L, Torraca S, Sirico ML, Guastaferro P, Chiuchiolo L, et al. Variability of blood pressure in dialysis patients: a new marker of cardiovascular risk. J Nephrol. (2013) 26:173–82. doi: 10.5301/jn.5000108
19. Amari Y, Morimoto S, Iida T, Yurugi T, Oyama Y, Aoyama N, et al. Characteristics of visit-to-visit blood pressure variability in hemodialysis patients. Hypertens Res. (2019) 42:1036–48. doi: 10.1038/s41440-019-0231-9
20. Kim HY, Kang YU, Kim CS, Choi JS, Bae EH, Ma SK, et al. Association of age and BP variability with long-term mortality in hemodialysis patients. Kidney Blood Press Res. (2013) 38:172–80. doi: 10.1159/000355765
21. Rossignol P, Cridlig J, Lehert P, Kessler M, Zannad F. Visit-to-visit blood pressure variability is a strong predictor of cardiovascular events in hemodialysis: insights from FOSIDIAL. Hypertension. (2012) 60:339–46. doi: 10.1161/HYPERTENSIONAHA.111.190397
22. Dolan E, O’Brien E. Blood pressure variability: clarity for clinical practice. Hypertension. (2010) 56:179–81. doi: 10.1161/HYPERTENSIONAHA.110.154708
23. Liu J, Huang Z, Gilbertson DT, Foley RN, Collins AJ. An improved comorbidity index for outcome analyses among dialysis patients. Kidney Int. (2010) 77:141–51. doi: 10.1038/ki.2009.413
24. Flythe JE, Xue H, Lynch KE, Curhan GC, Brunelli SM. Association of mortality risk with various definitions of intradialytic hypotension. J Am Soc Nephrol. (2015) 26:724–34. doi: 10.1681/ASN.2014020222
25. Mehlum MH, Liestøl K, Kjeldsen SE, Julius S, Hua TA, Rothwell PM, et al. Blood pressure variability and risk of cardiovascular events and death in patients with hypertension and different baseline risks. Eur Heart J. (2018) 39:2243–51. doi: 10.1093/eurheartj/ehx760
26. Gosmanova EO, Mikkelsen MK, Molnar MZ, Lu JL, Yessayan LT, Kalantar-Zadeh K, et al. Association of systolic blood pressure variability with mortality, coronary heart disease, stroke, and renal disease. J Am Coll Cardiol. (2016) 68:1375–86. doi: 10.1016/j.jacc.2016.06.054
27. Wang Y, Qin Y, Fan X, Cai J, Ye W, Xia J, et al. Variability in predialysis systolic blood pressure and long-term outcomes in hemodialysis patients. Kidney Blood Press Res. (2018) 43:115–24. doi: 10.1159/000487111
28. Muntner P, Shimbo D, Tonelli M, Reynolds K, Arnett DK, Oparil S. The relationship between visit-to-visit variability in systolic blood pressure and all-cause mortality in the general population: findings from NHANES III, 1988 to 1994. Hypertension. (2011) 57:160–6. doi: 10.1161/HYPERTENSIONAHA.110.162255
29. Flythe JE, Inrig JK, Shafi T, Chang TI, Cape K, Dinesh K, et al. Association of intradialytic blood pressure variability with increased all-cause and cardiovascular mortality in patients treated with long-term hemodialysis. Am J Kidney Dis. (2013) 61:966–74. doi: 10.1053/j.ajkd.2012.12.023
30. Parati G, Ochoa JE, Lombardi C, Bilo G. Blood pressure variability: assessment, predictive value, and potential as a therapeutic target. Curr Hypertens Rep. (2015) 17:537. doi: 10.1007/s11906-015-0537-1
31. Parati G, Liu X, Ochoa JE, Bilo G. Prognostic relevance of blood pressure variability. Am J Kidney Dis. (2013) 62:682–4. doi: 10.1161/HYPERTENSIONAHA.113.01801
32. Liao R, Wang L, Li J, Lin L, Sun S, Xiong Y, et al. Hemodialysis access type is associated with blood pressure variability and echocardiographic changes in end-stage renal disease patients. J Nephrol. (2019) 32:627–34. doi: 10.1007/s40620-018-00574-y
33. Roy-Chaudhury P, Sukhatme VP, Cheung AK. Hemodialysis vascular access dysfunction: a cellular and molecular viewpoint. J Am Soc Nephrol. (2006) 17:1112–27. doi: 10.1681/ASN.2005050615
34. Cecchi E, Giglioli C, Valente S, Lazzeri C, Gensini GF, Abbate R, et al. Role of hemodynamic shear stress in cardiovascular disease. Atherosclerosis. (2011) 214:249–56. doi: 10.1016/j.atherosclerosis.2010.09.008
35. Wu CC, Wen SC, Yang CW, Pu SY, Tsai KC, Chen JW. Plasma ADMA predicts restenosis of arteriovenous fistula. J Am Soc Nephrol. (2009) 20:213–22. doi: 10.1681/ASN.2008050476
36. Wu CC, Chen TY, Hsieh MY, Lin L, Yang CW, Chuang SY, et al. Monocyte chemoattractant protein-1 levels and postangioplasty restenosis of arteriovenous fistulas. Clin J Am Soc Nephrol. (2017) 12:113–21. doi: 10.2215/CJN.04030416
37. Nagai M, Dote K, Kato M, Sasaki S, Oda N, Kagawa E, et al. Visit-to-visit blood pressure variability, average BP level and carotid arterial stiffness in the elderly: a prospective study. J Hum Hypertens. (2017) 31:292–8. doi: 10.1038/jhh.2016.77
38. Mancia G, Facchetti R, Parati G, Zanchetti A. Visit-to-visit blood pressure variability, carotid atherosclerosis, and cardiovascular events in the european lacidipine study on atherosclerosis. Circulation. (2012) 126:569–78. doi: 10.1161/CIRCULATIONAHA.112.107565
39. Diaz KM, Veerabhadrappa P, Kashem MA, Thakkar SR, Feairheller DL, Sturgeon KM, et al. Visit-to-visit and 24-h blood pressure variability: association with endothelial and smooth muscle function in African Americans. J Hum Hypertens. (2013) 27:671–7. doi: 10.1038/jhh.2013.33
40. Li H, Xue J, Dai W, Liao X, Zhu P, Zhou Q, et al. Blood pressure variability and outcomes in end-stage renal disease patients on dialysis: a systematic review and meta-analysis. Kidney Blood Press Res. (2020) 45:631–44. doi: 10.1159/000508979
Keywords: blood pressure variability, hemodialysis, vascular access, thrombosis, stenosis, blood pressure, follow-up studies, renal dialysis
Citation: Hsieh M-Y, Cheng C-H, Chen C-H, Liao M-T, Lin C-C, Yang T-F, Chuang S-Y and Wu C-C (2022) The association of long-term blood pressure variability with hemodialysis access thrombosis. Front. Cardiovasc. Med. 9:881454. doi: 10.3389/fcvm.2022.881454
Received: 22 February 2022; Accepted: 18 July 2022;
Published: 04 August 2022.
Edited by:
Jeroen Peter Kooman, Maastricht University Medical Centre, NetherlandsReviewed by:
Ana Luisa Palacios Acedo, INSERM U1263 Centre de Recherche en Cardiovasculaire et Nutrition, FranceDana Ciobanu, Iuliu Haţieganu University of Medicine and Pharmacy, Romania
Copyright © 2022 Hsieh, Cheng, Chen, Liao, Lin, Yang, Chuang and Wu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Chih-Cheng Wu, chihchengwumd@gmail.com
†These authors have contributed equally to this work