- Open Access
Risk factors associated with hemodialysis central venous catheter malfunction; a retrospective analysis of a randomized controlled trial
Canadian Journal of Kidney Health and Diseasevolume 1, Article number: 15 (2014)
We previously reported a reduction in central venous catheter (CVC) malfunction when using once-weekly recombinant tissue-plasminogen activator (rt-PA) as a locking solution, compared with thrice-weekly heparin.
To identify risk factors for CVC malfunction to inform a targeted strategy for rt-PA use.
Canadian hemodialysis (HD) units.
Adults with newly placed tunnelled upper venous system CVCs randomized to a locking solution of rt-PA(1 mg/mL) mid-week and heparin (5000 u/ml) on the other HD sessions, or thrice-weekly heparin (5000 u/ml).
CVC malfunction (the primary outcome) was defined as: peak blood flow less than 200 mL/min for thirty minutes during a HD session; mean blood flow less than 250 mL/min for two consecutive HD sessions; inability to initiate HD.
Cox regression was used to determine the association between patient demographics, HD session CVC-related variables and the outcome of CVC malfunction.
Patient age (62.4 vs 65.4 yr), proportion female sex (35.6 vs 48.4%), and proportion with first catheter ever (60.7 vs 61.3%) were similar between patients with and without CVC malfunction. After multivariate analysis, risk factors for CVC malfunction were mean blood processed < 65 L when compared with ≥ 85 L in the prior 6 HD sessions (HR 4.36; 95% CI, 1.59 to 11.95), and mean blood flow < 300 mL/min, or 300 – 324 mL/min in the prior 6 HD sessions (HR 7.65; 95% CI, 2.78 to 21.01, and HR 5.52; 95% CI, 2.00 to 15.23, respectively) when compared to ≥ 350 mL/min.
This pre-specified post-hoc analysis used a definition of CVC malfunction that included blood flow, which may result in an overestimate of the effect size. Generalizability of results to HD units where trisodium citrate locking solution is used may also be limited.
HD session characteristics including mean blood processed and mean blood flow were associated with CVC malfunction, while patient characteristics were not. Whether targeting these patients at greater risk of CVC malfunction with rt-PA as a locking solution improves CVC longevity remains to be determined.
Nous avons précédemment fait rapport d’une réduction du dysfonctionnement du cathéter veineux central (CVC) lors de l’utilisation hebdomadaire de l’activateur tissulaire du plasminogène obtenu par génie génétique (rt-PA) comme solution verrou, plutôt que l’administration d’héparine trihebdomadaire.
Déterminer les facteurs de risques de dysfonctionnement du CVC afin d’indiquer une stratégie visée quant à l’utilisation du rt-PA.
Les services canadiens d’hémodialyse (HD).
Adultes à qui on a nouvellement installé un CVC tunnellisé intravasculaire dans le système veineux supérieur, et qui ont reçu soit une solution verrou de rt-PA (1 mg/ml) en milieu de semaine et de l’héparine (5 000 u/ml) lors des autres séances d’HD, soit de l’héparine trois fois par semaine (5 000 u/ml).
On a défini le dysfonctionnement du CVC (résultat primaire) comme étant : un débit sanguin de pointe inférieur à 200 mL/min durant trente minutes, au cours d’une séance d’HD; un débit sanguin moyen inférieur à 250 mL/min lors de deux séances d’HD consécutives; l’impossibilité d’entamer l’HD.
On a eu recours au modèle de régression de Cox pour déterminer l’association entre les données démographiques des participants, les variables relatives au CVC lors des séances d’HD et le résultat d’un dysfonctionnement du CVC.
L’âge des participants (62,4 c. 65,4 ans), la proportion des participants de sexe féminin (35,6% c. 48,4%), et la proportion de ceux à qui on a installé un cathéter pour la première fois (60,7% c. 61,3%) étaient similaires entre les patients qui ont subi un dysfonctionnement du CVC et ceux qui n’en ont pas subi. L’analyse multifactorielle révèle que les facteurs de risque liés au dysfonctionnement du CVC sont un traitement moyen du sang < 65 L, comparativement à ≥ 85 L au cours des 6 séances d’HD précédentes (HR 4,36; 95% CI, 1,59 à 11,95), et un débit sanguin moyen < 300 mL/min, ou 300 – 324 mL/min lors des 6 séances d’HD précédentes (HR 7,65; 95% CI, 2,78 à 21,01, et HR 5,52; 95% CI, 2,00 à 15,23, respectivement), comparativement à ≥ 350 mL/min.
Limites de l’étude
Cette analyse prédéterminée et post-hoc reposait sur une définition du dysfonctionnement du CVC qui comprenait le débit sanguin, ce qui pourrait résulter en une surestimation de l’ampleur de l’effet réel. La validité externe des résultats pour les services d’HD qui utilisent le citrate trisodique comme solution verrou pourrait aussi être limitée.
Les caractéristiques des séances d’HD comprenant le traitement moyen du sang et le débit sanguin moyen ont été associées au dysfonctionnement du CVC, alors que les caractéristiques des participants ne l’ont pas été. Il reste à déterminer si le fait de cibler les patients courant un risque accru de dysfonctionnement du CVC avec rt-PA comme solution verrou améliore la longévité du CVC.
What was known before
Central venous catheters (CVC) are commonly used to establish access in adult North American hemodialysis patients. CVC malfunction is a common clinically encountered complication of the use of tunnelled CVC. Previous authors have identified patient and CVC characteristics associated with malfunction.
What this adds
In this pre-specified, post-hoc analysis of data from a randomized controlled study we found that patient characteristics were not associated with risk of CVC malfunction after adjusting for hemodialysis session characteristics. Further research is required to predict CVC malfunction and strategies for prevention.
Current guidelines recommend hemodialysis (HD) vascular access via an arterio-venous fistula. However, central venous catheters (CVC) are used by approximately 70% of incident HD patients in North America [1–5]. Complications associated with CVC use include thrombosis and decreased dialysis adequacy due to catheter malfunction [6–8]. These complications are associated with increased healthcare resource utilization; up to 50% of all tunnelled CVC’s fail within one year [9–12].
CVC locking solutions are used in the interdialytic period to decrease the risk of complications; the optimal strategy to minimize CVC complications remains to be determined [13–16]. We previously reported that using once-weekly rt-PA as a locking solution leads to a two-fold reduction in CVC malfunction, compared with thrice-weekly heparin . However, given the potential costs associated with this strategy, its use and uptake may be limited. Identifying patients at greatest risk of CVC malfunction would permit targeted use of this strategy. Using data from the prior randomized trial  we sought to determine the association between patient characteristics and CVC-related variables during a dialysis session and the outcome of CVC malfunction.
The subjects and design of the original trial have been described previously . In summary, adult HD patients with newly inserted tunnelled, upper venous system HD catheters were eligible. Subjects were excluded if they were on systemic anticoagulation, undergoing current antibiotic use for catheter related bacteremia, the CVC was inserted via guidewire exchange, or they had risk factors for increased bleeding. If patients with catheter-related bacteremia completed antibiotic treatment, and remained bacteremia free for three hemodialysis sessions, they were eligible. CVC reversal was permitted for up to three consecutive hemodialysis sessions after which the study protocol mandated the CVC run in the non-reversed orientation and the patient was then eligible for the primary outcome.
The study was conducted in adherence with the Declaration of Helsinky, approved by the research ethics board of each institution, and patients provided written, informed consent. Patients were randomized to either thrice-weekly catheter locking solution with heparin (5000 u/mL, full luminal volume), or substitution of the mid-week heparin with recombinant tissue plasminogen activator (rt-PA) (1 mg in 1 mL with saline to full luminal volume) and standard heparin catheter locking solution following the other dialysis sessions. All other HD care was per site-specific protocols.
CVC malfunction (the primary outcome) was defined as: peak blood flow less than 200 mL/min for thirty minutes during a HD session; mean blood flow less than 250 mL/min for two consecutive HD sessions; inability to initiate HD . Each patient was followed-up for six months. Patients were censored at the first of: study end, removal of the CVC, initiation of alternative renal replacement therapy, or transfer to a centre not involved in the trial.
Candidate risk factors for CVC malfunction
Potential risk factors included patient characteristics (age, gender, race, duration of dialysis, cause of ESRD), baseline blood work (hemoglobin, platelets, albumin), comorbidity (prior myocardial infarction, congestive heart failure, diabetes, pulmonary embolism, gastrointestinal bleed, deep venous thrombosis, hypertension, cancer, current smoker), type of treatment (heparin-only or rt-PA), and CVC characteristics (any previous CVC, number of prior CVCs, insertion site, indication). We also considered HD-specific variables (mean blood processed, mean blood flow, mean arterial pressure, mean venous pressure) in the immediate HD session, as well as the 6 HD sessions, prior to CVC malfunction (for those with a primary outcome) or the same time period prior to censoring (for those without a primary outcome). We adjusted for CVC reversal given its potential role as a confounder, and assessed its role as a risk factor in a univariate analysis.
Cox proportional hazards regression analyses for time to catheter malfunction (for those with an event) or censoring (for those without an event) were performed, initially with all candidate risk factors for CVC malfunction assessed in a univariate fashion, followed by multivariate analyses with adjustment for type of treatment, age, gender, first dialysis catheter ever, diabetes, history of rt-PA use with a prior CVC and CVC reversal. All statistical tests were two-sided, and reported P values were considered significant if less than 0.05. All analyses were performed using STATA 11.2 (Stata Corporation, College Station, Texas).
We included all 225 patients from the original randomized trial in this study. Baseline characteristics were similar between patients with CVC malfunction and patients with no CVC malfunction (Table 1). One hundred sixty-three patients (72.4%) had no CVC malfunction and sixty-two (27.6%) patients had malfunction. The median duration of follow-up was 162 days for patients with no malfunction and 48.5 days for patients with CVC malfunction.
The study period included 24,085 CVC days; 3,048 in the patients with CVC malfunction and 21,037 in the patients with no CVC malfunction. Of the 62 patients with CVC malfunction, 40 (64.5%) were treated with heparin only. Of the 163 patients with no CVC malfunction, 75 (46.0%) were treated with heparin only. Univariate analysis found a significant, increased risk of CVC malfunction when HD lines were reversed at the session immediately prior to malfunction (hazard ratio [HR] 11.2, 95% CI 6.35-19.7, p < 0.001), as well as at least once in the six sessions prior to malfunction (HR 9.50, 95% CI 4.49 – 20.10, p < 0.001). Within the 6 antecedent HD sessions to the primary outcome, risk increased with an increasing number of sessions (to a maximum of three as per protocol) with CVC reversal (HR 1.72; 95% CI 1.51 – 1.96, p < 0.001), and in patients with a history of rt-PA use for catheter malfunction (HR 1.89 95% CI 1.04-3.46 p = 0.038). The risk of CVC malfunction decreased as blood processed or mean blood flow increased (p < 0.001). A trend to increased risk was found for female patients (HR 1.57 95% CI 0.95-2.59 p = 0.075). Age, diabetes, dialysis duration, first CVC status, hemoglobin, and number of CVC in the prior year were not significantly associated with malfunction (Table 2).
In a multivariable analysis the only variables significantly associated with CVC malfunction were HD session characteristics. After adjusting for baseline characteristics (age, gender, diabetes, first dialysis catheter ever), treatment allocation, history of rt-PA use, and CVC line reversal, patients with mean blood processed < 65 L had significantly higher risk of CVC malfunction when compared with patients with mean blood processed ≥ 85 L in the prior 6 HD sessions (HR 4.36; 95% CI, 1.59 to 11.95). When compared with patients with mean blood flow (averaged over the antecedent six sessions) ≥ 350 mL/min, those with mean blood flow < 300 mL/min and those with mean blood flow 300 to 324 mL/min had significantly higher risk of CVC malfunction (HR 7.65; 95% CI, 2.78 to 21.01 for the former comparison, and HR 5.52; 95% CI, 2.00 to 15.23 for the latter comparison) (Table 3).
In this pre-specified post-hoc analysis of data from a randomized controlled study we quantified the risk for CVC malfunction conferred by routinely observed and monitored patient and HD session characteristics. We found a significant increased risk when a reduction in litres of blood processed and mean blood flow over a two-week period occurred, however there was no association between baseline patient characteristics and risk of CVC malfunction.
CVC malfunction is the most common cause for CVC removal [20–22]. However, few studies have identified and quantified risk factors or predictors of CVC malfunction. Previous studies of patients with tunnelled CVCs reported that diabetic status and prior CVC exposure were associated with malfunction [12, 23]. However, data lacking on malfunction definition, locking solution, and the inclusion of infections and non-infectious criteria for CVC removal make extrapolation of these results difficult. While CVC location has been identified as a risk factor for malfunction, this study examined non-tunnelled catheters only . An observational cohort study of 3,364 incident and prevalent HD patients with a CVC in the United States reported the risk of CVC malfunction (blood flow <300 mL/min) was lower for males and black race (compared with white) and higher if CVC was not the first access modality . Male patients in our study were also at decreased risk but, possibly due to small numbers, we saw no significant effect of first CVC ever or ethnicity.
We analysed six sessions (a two week period) as we felt it would encompass a reasonable timeframe for clinicians to review dialysis characteristics in real-world practice; potentially guiding decision-making to pre-empt malfunction. In univariate analysis we found the strongest risk factor for CVC malfunction was CVC reversal in the prior HD session, or any of the six sessions prior to malfunction. However, our study protocol limited the number of consecutive sessions that could be performed in the reversed position to three, therefore this variable was included as a potential confounder in the adjusted analysis rather than an independent predictor given its potential relationship with the outcome. KDIGO Guidelines  suggest that CVC reversal be a temporary solution to managing malfunction, and authors have found that one third of patients may routinely undergo dialysis with their catheters in the reversed position . This prevalence may be influenced by studies questioning the impact of CVC line reversal or blood flow on dialysis adequacy. Carson et al. noted CVC reversal improved clearance, despite increased recirculation, if the blood flow through the CVC concurrently increased to an average 315 mL/min  and Moist et al. reported decreased blood flow through reversed CVCs did not significantly affect online urea clearance if blood flow was greater than 250 mL/min . Our results corroborate these findings as adjusted subgroup analysis by blood flow or blood processed quartile found a statistically significant increased risk when blood flows were less than 325 mL/min or total blood processed was less than 65 L per session in the preceding six sessions to a malfunction event.
We sought to pragmatically identify HD session characteristics that may signal an increased risk of CVC malfunction, and warrant consideration for early intervention. Decreased blood flow below 325 mL/min or total blood processed less than 65 L per session averaged over six sessions may therefore signal early or mild CVC malfunction. Mechanistically, it is plausible that a developing CVC thrombus or fibrin sheath, partially obstructing the CVC, and will decrease the average blood flow and blood processed. Interventions such as thrombolytics, catheter stripping, or catheter exchange may be required to improve CVC function, although these interventions are associated with costs, and long term catheter patency rates are suboptimal irrespective of intervention chosen [29, 30]. CVC reversal, decreased blood flow, or decreased blood processed may herald the future onset of severe malfunction (as defined by blood flow criteria), although whether interventions to identify and treat these patients is feasible, improves outcomes, or is cost effective requires further investigation.
Our study has limitations. First, this is a post-hoc analysis, and the definition of CVC malfunction included blood flow as a criteria (a significant predictor of the malfunction outcome), which may lead to an overestimate of the true effect of blood flow as a predictor. Secondly, our protocol limited the number of sessions that could be performed with a CVC in the reversed position. At the time of our study, CVC reversal was believed to result in a reduction in dialysis adequacy. CVC reversal was the strongest risk factor in our univariate analysis, therefore future studies could prospectively evaluate this association over longer time periods. Finally, the generalizability of these results to other units where trisodium citrate is used as a locking solution may be limited.
In conclusion, we found that HD session characteristics related to blood flow and blood processed were associated with an increased risk of CVC malfunction, while patient characteristics were not. Further study is needed to determine if the routine use of rt-PA as a locking solution in patients at increased risk of CVC malfunction will result in clinical and economic benefit.
Canadian Organ Replacement Register Annual Report: Treatment of End-Stage Organ Failure in Canada, 2001 to 2010. In Book Canadian Organ Replacement Register Annual Report: Treatment of End-Stage Organ Failure in Canada, 2001 to 2010. Ottawa, Ont: CIHI; 2011.
System USRD: USRDS 2012 Annual Data Report: Atlas of Chronic Kidney Disease and End-Stage Renal Disease in the United States. Book USRDS 2012 Annual Data Report: Atlas of Chronic Kidney Disease and End-Stage Renal Disease in the United States 2012.
III. NKF-K/DOQI: Clinical Practice Guidelines for Vascular Access: update 2000. Am J Kidney Dis 2001, 37: S137-S181.
Jindal K, Chan CT, Deziel C, Hirsch D, Soroka SD, Tonelli M, Culleton BF: Hemodialysis clinical practice guidelines for the Canadian Society of Nephrology. J Am Soc Nephrol 2006, 17: S1-S27. 10.1681/ASN.2005121372
Mendelssohn DC, Ethier J, Elder SJ, Saran R, Port FK, Pisoni RL: Haemodialysis vascular access problems in Canada: results from the Dialysis Outcomes and Practice Patterns Study (DOPPS II). Nephrol Dial Transplant 2006, 21: 721–728. 10.1093/ndt/gfi281
Vanholder V, Hoenich N, Ringoir S: Morbidity and mortality of central venous catheter hemodialysis: a review of 10 years' experience. Nephron 1987, 47: 274–279. 10.1159/000184523
Develter W, De Cubber A, Van Biesen W, Vanholder R, Lameire N: Survival and complications of indwelling venous catheters for permanent use in hemodialysis patients. Artif Organs 2005, 29: 399–405. 10.1111/j.1525-1594.2005.29067.x
Schon D, Whittman D: Managing the complications of long-term tunneled dialysis catheters. Semin Dial 2003, 16: 314–322. 10.1046/j.1525-139X.2003.16060.x
Polkinghorne KR, McDonald SP, Atkins RC, Kerr PG: Vascular access and all-cause mortality: a propensity score analysis. J Am Soc Nephrol 2004, 15: 477–486. 10.1097/01.ASN.0000109668.05157.05
Dhingra RK, Young EW, Hulbert-Shearon TE, Leavey SF, Port FK: Type of vascular access and mortality in U.S. hemodialysis patients. Kidney Int 2001, 60: 1443–1451. 10.1046/j.1523-1755.2001.00947.x
Lee H, Manns B, Taub K, Ghali WA, Dean S, Johnson D, Donaldson C: Cost analysis of ongoing care of patients with end-stage renal disease: the impact of dialysis modality and dialysis access. Am J Kidney Dis 2002, 40: 611–622. 10.1053/ajkd.2002.34924
Little MA, O'Riordan A, Lucey B, Farrell M, Lee M, Conlon PJ, Walshe JJ: A prospective study of complications associated with cuffed, tunnelled haemodialysis catheters. Nephrol Dial Transplant 2001, 16: 2194–2200. 10.1093/ndt/16.11.2194
Ivan DM, Smith T, Allon M: Does the heparin lock concentration affect hemodialysis catheter patency? Clin J Am Soc Nephrol 2010, 5: 1458–1462. 10.2215/CJN.01230210
Thomas CM, Zhang J, Lim TH, Scott-Douglas N, Hons RB, Hemmelgarn BR: Concentration of heparin-locking solution and risk of central venous hemodialysis catheter malfunction. Asaio J 2007, 53: 485–488. 10.1097/MAT.0b013e3180619519
Weijmer MC, Van Den Dorpel MA, Van De Ven PJ, Ter Wee PM, Van Geelen JA, Groeneveld JO, Van Jaarsveld BC, Koopmans MG, Le Poole CY, Schrander Van Der Meer AM, Siegert CE, Stas KJ: Randomized, clinical trial comparison of trisodium citrate 30% and heparin as catheter-locking solution in hemodialysis patients. J Am Soc Nephrol 2005, 16: 2769–2777. 10.1681/ASN.2004100870
Yevzlin AS, Sanchez RJ, Hiatt JG, Washington MH, Wakeen M, Hofmann RM, Becker YT: Concentrated heparin lock is associated with major bleeding complications after tunneled hemodialysis catheter placement. Semin Dial 2007, 20: 351–354. 10.1111/j.1525-139X.2007.00294.x
Hemmelgarn BR, Moist LM, Lok CE, Tonelli M, Manns BJ, Holden RM, LeBlanc M, Faris P, Barre P, Zhang J, Scott-Douglas N: Prevention of dialysis catheter malfunction with recombinant tissue plasminogen activator. N Engl J Med 2011, 364: 303–312. 10.1056/NEJMoa1011376
Hemmelgarn BR, Moist L, Pilkey RM, Lok C, Dorval M, Tam PY, Berall MJ, LeBlanc M, Toffelmire EB, Manns BJ, Scott-Douglas N: Prevention of catheter lumen occlusion with rT-PA versus heparin (Pre-CLOT): study protocol of a randomized trial [ISRCTN35253449]. BMC Nephrol 2006, 7: 8. 10.1186/1471-2369-7-8
National Kidney Foundation: KDOQI Clinical Practice Guidelines and Clinical Practice Recommendations for 2006 Updates: Hemodialysis Adequacy, Peritoneal Dialysis Adequacy and Vascular Access. Am J Kidney Dis 2006,48(Suppl 1):S1-S322.
Swartz RD, Messana JM, Boyer CJ, Lunde NM, Weitzel WF, Hartman TL: Successful use of cuffed central venous hemodialysis catheters inserted percutaneously. J Am Soc Nephrol 1994, 4: 1719–1725.
Trerotola SO, Johnson MS, Harris VJ, Shah H, Ambrosius WT, McKusky MA, Kraus MA: Outcome of tunneled hemodialysis catheters placed via the right internal jugular vein by interventional radiologists. Radiology 1997, 203: 489–495. 10.1148/radiology.203.2.9114110
Moureau N, Poole S, Murdock MA, Gray SM, Semba CP: Central venous catheters in home infusion care: outcomes analysis in 50,470 patients. J Vasc Interv Radiol 2002, 13: 1009–1016. 10.1016/S1051-0443(07)61865-X
Cetinkaya R, Odabas AR, Unlu Y, Selcuk Y, Ates A, Ceviz M: Using cuffed and tunnelled central venous catheters as permanent vascular access for hemodialysis: a prospective study. Ren Fail 2003, 25: 431–438. 10.1081/JDI-120021155
Hryszko T, Brzosko S, Mazerska M, Malyszko J, Mysliwiec M: Risk factors of nontunneled noncuffed hemodialysis catheter malfunction. A prospective study. Nephron Clin Pract 2004, 96: c43-c47. 10.1159/000076398
Griffiths RI, Newsome BB, Block GA, Herbert RJ, Danese MD: Patterns of Hemodialysis Catheter Dysfunction Defined According to National Kidney Foundation Guidelines As Blood Flow <300 mL/min. Int J Nephrol 2011, 2011: 891259.
Pannu N, Jhangri GS, Tonelli M: Optimizing dialysis delivery in tunneled dialysis catheters. Asaio J 2006, 52: 157–162. 10.1097/01.mat.0000202081.13974.39
Carson RC, Kiaii M, MacRae JM: Urea clearance in dysfunctional catheters is improved by reversing the line position despite increased access recirculation. Am J Kidney Dis 2005, 45: 883–890. 10.1053/j.ajkd.2005.01.029
Moist LM, Hemmelgarn BR, Lok CE: Relationship between blood flow in central venous catheters and hemodialysis adequacy. Clin J Am Soc Nephrol 2006, 1: 965–971. 10.2215/CJN.01090306
Janne D'Othee B, Tham JC, Sheiman RG: Restoration of patency in failing tunneled hemodialysis catheters: a comparison of catheter exchange, exchange and balloon disruption of the fibrin sheath, and femoral stripping. J Vasc Interv Radiol 2006, 17: 1011–1015. 10.1097/01.RVI.0000223692.99800.2E
Macrae JM, Loh G, Djurdjev O, Shalansky S, Werb R, Levin A, Kiaii M: Short and long alteplase dwells in dysfunctional hemodialysis catheters. Hemodial Int 2005, 9: 189–195. 10.1111/j.1492-7535.2005.01131.x
We would like to thank Dr Adeera Levin, MD, Department of Medicine, University of British Columbia, Vancouver, British Columbia Canada for her input and advice regarding the project design.
The original PreCLOT trial was funded by Hoffman La-Roche Ltd.
DW assisted in the interpretation of data and drafting and critical revision of the manuscript. LM assisted in the conception and design of the study and participated in critical revision of the manuscript. JM participated in the conception and design of the study and participated in critical revision of the manuscript. NSD participated in the conception and design of the study and participated in critical revision of the article. JZ participated in the data collection and analysis and participated in critical revision of the manuscript. MT participated in the conception and design of the study and participated in critical revision of the manuscript. CL participated in conception and design of the study and participated in critical revision of the manuscript. SS participated in conception and design of the study and participated in critical revision of the manuscript. BH assisted in the conception and design of the study, data interpretation, and participated in drafting and critical revision of the manscript. All authors read and approved the final manuscript.