Application of DSM-BIA in dry weight assessment in continuous ambulatory peritoneal dialysis

doi:10.21203/rs.3.rs-770259/v1

Objectives

Dry weight assessment (DWA) is an important part of dialysis and fluid management in patients receiving renal replacement therapy. With the development of bioimpedance analysis (BIA), the development of the direct segmental multi-frequency BIA (DSM-BIA) has provided a more convenient measure for DWA of dialysis patients, but its accuracy remains unclear. This study was designed to evaluate the application of DSM-BIA in DWA of continuous ambulatory peritoneal dialysis (CAPD) patients.

Design

This is a cross-sectional study. Using the conventional BIA as a reference, we examined the accuracy of the DSM-BIA technique for assessing dry weight in CAPD patients and analyzed the potential factors influencing their fluid volume status.

Setting and Participants:

A total of 31 patients with end-stage renal disease receiving CAPD and 310 healthy volunteers were recruited for this study.

Methods

The intraclass correlation coefficients (ICC) and Bland-Altman plots were used to assess the consistency between DSM-BIA and the conventional BIA for DWA. Univariate and multivariate linear regression analyses were used to explore the influencing factors associated with the edema index.

Results

DSM-BIA and the conventional BIA technology were consistent in DWA in CAPD patients (ICC female 0.972, ICC male 0.882, ICC total 0.960). Similarly, Bland-Altman plots showed good agreements between the two methods in DWA for both genders. Univariate and multivariate linear regression analysis showed both eGFR level (P = 0.04) and serum NT-pro BNP concentration (P = 0.007) were positively correlated with the ratio of extracellular water to total body water (ECW/TCW).

Conclusions

DSM-BIA in DWA has good accuracy in clinical applications and has potential application value for fluid volume management in CAPD patients.

Biological Chemistry

Health Policy

Urology & Nephrology

Dry weight assessment

direct segmental multi-frequency BIA

continuous ambulatory peritoneal dialysis

fluid management

Adequate control of extracellular fluid volume in patients receiving renal replacement therapy (RRT) is one of the main objectives of dialysis therapy, so clinicians need to make an accurate assessment of the patient's fluid status. Dry weight assessment (DWA), defined as a patient's post-dialysis weight while maintaining normal blood pressure during dialysis without the use of hypertension medications[1], is an important method for fluid volume management in dialysis patients. Inaccurate DWA will lead to corresponding hypervolemia (headache, hypertension, dyspnea, and orthopnea) or hypovolemia (cramps, fatigue, and orthostatic hypotension). However, low sensitivity to these symptoms and high variability among patients tend to lead to more severe volume depletion or volume overload, which are associated with a higher risk of cardiovascular events and death[2].

Traditionally, DWA has been achieved through some clinical examinations, including the usual blood pressure measurement and target dry weight titration until normal blood volume, which is entirely dependent on the expertise and skills of the clinician[5].The clinical assessment is far more complex and includes vascular stiffness, cardiac insufficiency, hypoproteinemia, and morbidity, as well as the desire to better preservation of residual kidney function, so more objective tools are needed to assess dry weight for patients receiving RRT.

Several objective methods have been proposed to support the correct estimation of dry weight in dialysis patients, including inferior vena cava ultrasound, radionuclide dilution techniques, and echocardiography. However, these methods have many disadvantages, such as time consumption, difficulty in daily practice, and inability to quantify excess or deficiency of body fluids. So far, most dialysis centers still rely on subjective clinical criteria to assess dry weight. Studies have concluded that bioelectrical impedance analysis (BIA) and bioimpedance spectroscopy (BIS) may be useful for the clinical assessment of fluid overload[5, 6] in patients receiving RRT due to their objectivity and portability. Both BIA and BIS are mainly based on the electrical conductivity of body water composition, calculating the water volume of body by measuring the impedance value of water by turning on a small amount of alternating current. The difference between BIA and BIS is that BIA is mainly based on built-in regression equations, while BIS is based on Cole modeling and mixture theories[5]. The early BIA is a single-frequency whole-body impedance measurement method, which must be compensated by empirical estimates. When the empirical estimation formula is applied to other types of groups, the accuracy is greatly reduced and can not be widely used. The whole ankle-wrist method, as the convitonal BIA technique, mainly measures the electrical impedance from the ipsilateral ankle to the wrist, and multiplies it to get the whole-body electrical impedance, which provides the possibility for the application of BIA technology in clinic. The defect of this convetional BIA is that the human body is regarded as a simple cylinder and only two electrodes are used to measure the electrical impedance from the wrist to the ankle of the unilateral limb, ignoring the physiological complexity of the human body structure[6]. Although the traditional single-frequency BIA method is not considered to be superior than BIS,with the improvement and development of direct segmental multi-frequency BIA (DSM-BIA) technology, DSM-BIA has certain advantages in the evaluation of human body composition. DSM-BIA divides the human body into five segments (left and right upper limbs, left and right lower limbs, trunk) to ensure the accuracy of impedance values of each segment, multi-frequency measurement (1-1000kHz) fully ensures the accuracy of extracellular and intracellular water volume, and the eight-point contact measurement method (left and right thumb contact point, left and right palm contact point, left and right forefoot contact point, left and right rear heel contact point) directly measures the trunk impedance without using the empirical estimation, which ensures the high repeatability and accuracy of the results[7]. Recently, DSM-BIA has shown superior results in the estimation of body composition[8] and fluid volume management of hemodialysis patients^[9]༌but the results in DWA of continuous ambulatory peritoneal dialysis (CAPD) patients are not very clear.. Therefore, we conducted this - study to investigate the accuracy of DSM-BIA in DWA of continuous ambulatory peritoneal dialysis (CAPD) patients and to explore the related factors affecting fluid status in patients with CAPD.

Study design and population

The study design was a prospective cross-sectional study of patients with ESRD who had received CAPD for more than 3 months. Participants of CAPD were mainly from patients regularly followed up in the dialysis center of our hospital. Inclusion criteria were as follows: 1) The patient received CAPD due to renal function decline to stage 5 of chronic kidney disease, and the duration of CAPD was longer than 3 months; 2) All patients were older than 18 years and younger than 75 years; 3) Serum albumin ≥ 30g/L; 4) Body mass index (BMI) ≤ 40kg/m². Exclusion criteria were as follows: 1) Patients with heart failure with cardiac function ≥NYHA III grade, and patients with liver failure or advanced tumor; 2) Patients with skin infection of limbs; 3) Patients with metal prosthesis or pacemaker installed in their bodies; 4) Vascular diseases of lower extremities (such as arteriosclerosis obliterans, diabetic foot and varicose veins); 5) A history of central venous and inferior vena cava stenosis. The research program was approved by the Institutional Review Board(the Ethics Committee of Dongzhimen Hospital Affiliated with Beijing University of Chinese Medicine). Each patient signed a written informed consent to participate in the study before being assessed. The study was performed in accordance with the ethical standards of the Declaration of Helsinki (1964) and its subsequent amendments.

Measurement of dry weight assessment by body composition analyzer

Measurement of dry weight assessment was performed using the Inbody S10 Multifrequency Bioelectrical Impedance analyzer (Biospace, Seoul, Korea). All participants were required to empty the gastrointestinal tract, avoid strenuous exercise or showers, and stop the infusion of nutrient solution for more than 3 hours before measurement. The indoor temperature should be kept at 20-25 ℃. All measurements were taken under the guidance of trained research nurses.

Direct segmental multi-frequency bioelectrical impedance analysis (DSM-BIA)

The DSM-BIA technique assumes that the human body consists of five interconnected cylinders that measure the impedance of the subject's trunk, arms, and legs at six different frequencies (1 kHz, 5 kHz, 50 kHz, 250 kHz, 500 kHz, and 1000 kHz). The electrical frequency spectrum was used to predict the intracellular water (ICW) and extracellular water (ECW) of total body water (TBW) in different body parts, and immediate and extensive quantitative values of various body composition parameters were also measured.

Conventional bioelectric impedance analysis

Conventional bioelectric impedance analysis relies on formulas to estimate dry weight. In this study, we used the formula proposed by Chamney to directly calculate dry weight using the results of InbodyS10 Multifrequency Bioelectrical Impedance analyzer[8], which is based on the assumption that patients will gradually approach the volumetric state of a healthy person as the body is dehydrated with dialysis, and the dry weight is calculated when the ECW% in body weight intersects with the slope of the normal value. The Chamney formula is as follows, where ECW% (ECW/TBW%) needs to be obtained from the corresponding normal population in the same BMI area group, Wgt_m is the bodyweight of patients before receiving CAPD, and ECW_m is the weight of ECW in patients measured by bioelectrical impedance analyzer.

Dry weight (kg) =（Wgt_m - ECW_m）/(1-ECW%)

Echocardiography and laboratory analysis

All patients were given 2.5% peritoneal dialysis solution at night for 8 hours. The peritoneal dialysis solution was released after the morning of the second day, and another 2.5% peritoneal dialysis solution was injected. Blood samples were collected 2 hours later, and the peritoneal dialysis solution was released 4 hours later, and the solution specimens were retained. The blood samples were used to test the blood routine and serum biochemical parameters, and the peritoneal dialysate samples were used to test the corresponding biochemical parameters. Echocardiography was performed in all patients after CAPD.

Statistical analyses

Normally distributed continuous variables are expressed as mean ± standard deviation (SD), skewed continuous variables are expressed as median and interquartile ranges, and categorical variables are expressed as frequency and percentage. The unpaired Student’s t-test was used to compare the continuous variables with normal distribution, and Mann-Whitney U test was used to test the continuous variables with non-normal distribution. The Bland-Altman plots were used to examine the consistency of the two methods in DWA. Pearson correlation was used to quantify the deviations seen on the Bland-Altman curve. To improve the clinical utility, univariate and multivariate linear regression analysis were used to assess the relevant factors influencing patients' fluid volume status. A two-sided P＜0.05 was considered statistically significant. GraphPad Prism 8.4.3 (GaphPadSoftware Inc, LaJolla, CA), and SPSS 16.0 (IBM Corp., Armonk, NY) were used for all statistical analyses.

Clinical characteristics of Study Participants

A total of 31 CAPD patients were enrolled in our study, and their clinical characteristics according to gender are summarized in Table 1. As expected, female subjects had significantly lower height, weight, and body surface area than male subjects, but no significant gender differences in body mass index (P=0.07). Moreover, the smoking and drinking histories of females were significantly lower than those of males. It should be noted that in our study, there was a significant imbalance in the solute transport capacity in the peritoneum of males and females, and the detailed distribution of peritoneal transport type is shown in Supplement table 1.

Agreements between DSM-BIA and Conventional BIA in DWA

We collected data from 310 healthy adults as a reference control group, whose clinical characteristics and body composition analysis results are presented in Supplement Table 2.

The DWA calculated by the DSM-BIA and Conventional BIA according to the different genders was shown in Table 2. By comparison, there are good agreements between two methods in DWA for total CAPD patients (intraclass correlation coefficients (ICC) 0.960, 95%CI 0.918 - 0.981; Table 2) and both genders. Furthermore, it was noted that the agreement of the two dry weight measures in women (ICC 0.970, 95% CI 0.880 - 0.993; Table 2) was slightly stronger than in men (ICC 0.882, 95%CI 0.709 - 0.952; Table 2) and in total CAPD patients.

The Bland-Altman plots were used to visually evaluated the agreement between the DSM-BIA and Conventional BIA in DWA, as shown in Figure 1. For all CAPD patients, 93.5% (29/31) of the scatter plots are within the estimated 95% CI range Figure 1A, suggesting good agreement between the two assessment methods in the measurements of DWA. The mean difference of DWA between DSM-BIA and the Conventional BIA is 3.2kg (SD of 4.0kg), which indicates that DSM-BIA overestimates the dry weight compared to Conventional BIA. There is no association was found between the difference (y-axis) across the range of -11.9kg to 5.2kg and the average (x-axis) across the range of 46.6kg to 82.7kg between the two methods for DWA (R² = 0.006, P=0.67; Figure 1A). The Bland-Altman plots were also used to examine the agreement of the two methods in the measurements of DWA in male and female CAPD patients, respectively. The results showed that all the data points (100%) of the two assessment methods are limited in the estimated 95% CI range (Figure 1B and Figure 1C), revealing the excellent agreement of the two assessment methods in DWA for both genders.

Factors associated with edema index in CAPD patients

To determine the possible factors associated with the fluid volume status of patients with CAPD, univariate and multivariate linear regression analyses were performed in the entire CAPD population in the study. The edema index (ECW/TBW) measured by the body composition analyzer was taken as the dependent variable, and the patient's blood test results, dialysate test results, and the parameters of the cardiac function evaluated by echocardiography were taken as the independent variables. On univariate analysis. as shown in Table 3 and Supplement Table 3, a statistically significant positive linear correlation was found on ECW/TBW with age (P=0.004), estimated glomerular filtration rate (eGFR) (P=0.04), C-reactive protein (P=0.04), NT-pro BNP (P=0.02), left ventricular end-systolic diameter (LVESD) (P=0.03) and left atrial anteroposterior diameter (LAAD) (P=0.005), whereas a statistically significant negative linear correlation was found on ECW/TBW with Platelets (P=0.04), left ventricular ejection fraction (LVEF) (P=0.02), and Left ventricular fractional shortening (LVFS) (P=0.02). No statistically significant correlation was found for other factors in CAPD patients (Supplement Table 3). In multivariable linear stepwise analysis, as can be seen in Table 3, NT-proBNP (P=0.007) and eGFR (P=0.04) were independently and positively correlated with ECW/TCW among the various factors examined. Additionally, the R² value is 0.409, which means that NT-proBNP and eGFR can explain a 40.9% change of ECW/TCW.

This study demonstrated that the DSM-BIA assessment of dry weight is in good agreement with the conventional BIA assessment method in CAPD patients regardless of the gender, and that volume overload control represented as edema index (ECW/TBW ratio) under the guidance of DSM-BIA is closely related to the eGFR level and serum NT-pro BNP concentration in CAPD patients. These results confirmed the value of the application of DSM-BIA in DWA and its usefulness in fluid volume management of CAPD patients.

It is well known that in patients with end-stage renal disease, the increase in systemic sodium and fluid volume is an inevitable consequence of the decrease in renal function, which is followed by the elevation of blood pressure and damage in blood vessels. These factors are closely related to the risk of serious cardiovascular complications[1]. Hypertension control in end-stage renal disease is one of the important objectives of treatment management. Studies have shown that reduced sodium excretion as the decline in renal function leads to hypertension in about 90% of dialysis patients, and that good control of sodium and fluid balance in patients can control hypertension without the use of antihypertensive drugs[9]. Given that positive sodium balance and the resulting hypervolemia are two major causes of hypertension and increased cardiovascular mortality in patients with end-stage renal disease, it is important to accurately assess the amount of fluid discharged from the body during dialysis. Therefore, achieving and maintaining dry weight is an effective strategy to control blood pressure and reduce the risk of cardiovascular events in patients receiving dialysis. Our study confirms that the DSM-BIA method has good accuracy in dry weight assessment in CAPD patients, which will provide clinicians with a more convenient and effective way to assess dry weight in clinical practice. A previous randomized controlled trial (RCT) involving hemodialysis patients showed that the guidance under BIA for dialysis significantly reduced all-cause mortality over 2.5 years follow-up period, and a recent three-years-extended exploratory study of a prospective RCT concluded that DSM-BIA guided fluid management improved hard outcomes in terms of patient survival and non-cardiovascular survival[6]. These findings provide more evidence and confidence for the clinical application of DSM-BIA and merit further study in longer follow-up periods and more carefully designed clinical studies.

BIA is an ideal tool for measuring the impact of interventions aimed at changing fluid status, so it is a good way to routinely monitor fluid status and nutrition in dialysis patients, helping them control fluid overload more effectively and more quickly. Clinicians need to be very cautious and combine clinical judgment in setting target bodyweight because there is increasing evidence from bioimpedance data that excessive hydration predicts poorer survival[5, 10]. A cross-sectional study has shown that DSM-BIA provides a good assessment of total body water to obtain an accurate assessment of hydration status[11]. DSM-BIA can accurately measure water in and out of cells and produce an "edema index", the ratio of ECW to TBW. Several cohort studies have shown that ECW/TBW ratio is strongly associated with poorer survival in PD patients, even after adjusting for known confounders such as comorbidities, inflammation, and hypoalbuminemia[12]. In our study, we observed that a lower ECW/TBW ratio was associated with a lower level of eGFR in CAPD patients. Studies have shown that excessive water itself leads to a decrease in residual renal function[13, 6], so it is necessary to better control fluid overload clinically to better preserve the residual renal function. Importantly, we also found that higher ECW/TBW ratios were strongly associated with higher NT-pro BNP levels. Elevated NT-pro BNP is known to be strongly associated with cardiac abnormalities[14], such as left ventricular hypertrophy, atrial enlargement, or estimated right ventricular systolic pressure. However, the correlation between cardiac biomarkers and hydration status remained highly heterogeneous, with some studies finding no correlation with hydration status, while others observed a correlation between change in ECW and NT-pro BNP[15], which also was consistent with our results.

Taken together, our findings suggest that DSM-BIA assessment of dry weight has good accuracy in clinical applications and has potential application value for fluid volume management in CAPD patients. It is necessary to conduct prospective large-scale randomized controlled trials to determine the impact of DSM-BIA guidance on long-term outcomes in CAPD patients.

Acknowledgements: We thank all of the CAPD patients and healthy volunteers who participated in this study, as well as all of the healthcare professionals who provided assistance for this research.

Conflicts of Interest: There are no conflicts of interest.

Funding sources: This research did not receive any funding from agencies in the public, commercial, or not-for-profit sectors.

Gunal, A. I. (2013) How to determine 'dry weight'? Kidney Int Suppl (2011) 3 (4):377–379. doi:10.1038/kisup.2013.81
Mehrotra, R., Devuyst, O., Davies, S. J. & Johnson, D. W. The Current State of Peritoneal Dialysis. J Am Soc Nephrol, 27 (11), 3238–3252 https://doi.org/10.1681/asn.2016010112 (2016).
Davies, S. J. Peritoneal dialysis–current status and future challenges. Nat Rev Nephrol, 9 (7), 399–408 https://doi.org/10.1038/nrneph.2013.100 (2013).
Jain, A. K., Blake, P., Cordy, P. & Garg, A. X. Global trends in rates of peritoneal dialysis. J Am Soc Nephrol, 23 (3), 533–544 https://doi.org/10.1681/asn.2011060607 (2012).
Davies, S. J. & Davenport, A. The role of bioimpedance and biomarkers in helping to aid clinical decision-making of volume assessments in dialysis patients. Kidney Int, 86 (3), 489–496 https://doi.org/10.1038/ki.2014.207 (2014).
Tian, N. et al. Bioimpedance Guided Fluid Management in Peritoneal Dialysis: A Randomized Controlled Trial. Clin J Am Soc Nephrol, 15 (5), 685–694 https://doi.org/10.2215/cjn.06480619 (2020).
Ling, C. H. et al. Accuracy of direct segmental multi-frequency bioimpedance analysis in the assessment of total body and segmental body composition in middle-aged adult population. Clin Nutr, 30 (5), 610–615 https://doi.org/10.1016/j.clnu.2011.04.001 (2011).
Chamney, P. W. et al. A whole-body model to distinguish excess fluid from the hydration of major body tissues. Am J Clin Nutr, 85 (1), 80–89 https://doi.org/10.1093/ajcn/85.1.80 (2007).
Li, Z. et al. The epidemiology of systolic blood pressure and death risk in hemodialysis patients. Am J Kidney Dis, 48 (4), 606–615 https://doi.org/10.1053/j.ajkd.2006.07.005 (2006).
Udayaraj, U. P. et al. Blood pressure and mortality risk on peritoneal dialysis. Am J Kidney Dis, 53 (1), 70–78 https://doi.org/10.1053/j.ajkd.2008.08.030 (2009).
Bedogni, G. et al. Accuracy of an eight-point tactile-electrode impedance method in the assessment of total body water. Eur J Clin Nutr, 56 (11), 1143–1148 https://doi.org/10.1038/sj.ejcn.1601466 (2002).
Koh, K. H., Wong, H. S., Go, K. W. & Morad, Z. Normalized bioimpedance indices are better predictors of outcome in peritoneal dialysis patients. Perit Dial Int, 31 (5), 574–582 https://doi.org/10.3747/pdi.2009.00140 (2011).
Tian, N. et al. The Impact of Fluid Overload and Variation on Residual Renal Function in Peritoneal Dialysis Patient. PLoS One, 11 (4), e0153115 https://doi.org/10.1371/journal.pone.0153115 (2016).
Wang, A. Y. et al. Diagnostic potential of serum biomarkers for left ventricular abnormalities in chronic peritoneal dialysis patients. Nephrol Dial Transplant, 24 (6), 1962–1969 https://doi.org/10.1093/ndt/gfp067 (2009).
Davenport, A. Changes in N-terminal pro-brain natriuretic peptide correlate with fluid volume changes assessed by bioimpedance in peritoneal dialysis patients. Am J Nephrol, 36 (4), 371–376 https://doi.org/10.1159/000343286 (2012).

Table 1. Comparisons of clinical Characteristics of Study Participants according to gender.

Characteristic	Male (N=21)	Female (N=10)	Total (N=31)	P value
Age, years	56.3±10.6	54.9±14.1	55.9±11.9	0.76
Dialysis vintage, months	31.5±18.2	39.2±23.1	34.0±19.6	0.32
History of smoking	14(67.0%)	1(10.0%)	15(48.3%)	0.006
History of drinking	12(57.1%)	1(10.0%)	13(41.9%)	0.02
Height, cm	172.4±3.8	154.5±4.7	166.0±9.4	<0.001
Weight, kg	71.3±87.6	53.6±6.3	64.8±11.0	<0.001
BMI, kg/m²	24.0±2.6	22.2±2.1	23.3±2.3	0.07
BSA, m²	1.8±0.1	1.5±0.1	1.7±0.2	<0.001
eGFR, ml/min per 1.73 m²	5.3±1.7	5.6±1.8	5.4±1.7	0.67
Ccr, ml/min	7.6±2.2	6.5±1.3	5.5±1.8	0.14
PET	0.7±0.1	0.6±0.1	0.6±0.1	0.01
UFV, ml	205.1±85.3	209.1±166.2	215.4±100.1	0.42
Systolic BP, mm Hg	144.6±16.8	142.0±20.7	146.7±17.7	0.72
Diastolic BP, mm Hg	84.5±10.6	84.6±11.0	83.8±10.5	0.98

Abbreviation: BMI, body mass index; BSA, body surface area; eGFR, estimated glomerular filtration rate; Ccr, creatinine clearance rate; PET, peritoneal equilibration test; UFV, Ultrafiltration volume; BP, blood pressure.

Table2. Dry weight assessment calculated by different methods according to gender.

Items	Male (N=21)	Female (N = 10)	Total (N = 31)
DSM-BIA	70.9±6.2	54.2±5.4	65.5±9.8
Conventional BIA	67.4±7.5	51.9±6.2	62.4±10.1
ICC	0.88(0.71 - 0.95)	0.97(0.88 - 0.99)	0.96(0.92 - 0.98）

Abbreviation: ICC, intraclass correlation coefficient.

Table 3. Linear regression analysis of parameters possibly associated with ECW/TCW.

Parameter	ANOVA（F）	Univariate Linear Regression （Std. β）	P value	Multivariate Linear Regression（Std. β）	P value
Age, year	9.6	136.3	0.004
Gender	0.4	6.1	0.53
eGFR, ml/min per 1.73 m²	4.5	21.2	0.04	0.003	0.04
BMI, kg/m²	0.2	34.4	0.65
MAP, mmHg	6.9	126.4	0.01
UFA, ml	0.5	1313.0	0.50
Dialysis vintage, month	0.1	258.1	0.76
PET	1.0	1.2	0.32
Clinical edema	5.6	8.1	0.03
Blood test results
Platelets, ×10⁹/L	4.5	-2070.0	0.04
C-reactive protein, mg/L	4.5	232.4	0.04
NT-pro BNP, pg/ml	5.9	16049.0	0.02	˂0.001	0.007
M-mode and Two-Dimensional Echocardiography
LVESD, mm	5.4	183.9	0.03
LVEF, %	6.4	-345.4	0.02
LVFS, %	7.0	-201.6	0.01
LAAD, mm	9.5	199.9	0.005

Abbreviation: eGFR, estimated glomerular filtration rate; BMI, body mass index; MAP, mean arterial pressure; UFV, Ultrafiltration volume; PET, peritoneal equilibration test; LVESD, Left ventricular end systolic diameter; LVEF, Left ventricular ejection fraction; LVFS, Left ventricular fractional shortening; LAAD, Left atrial anteroposterior diameter;

No competing interests reported.

Supplementarymaterial.docx

Application of DSM-BIA in dry weight assessment in continuous ambulatory peritoneal dialysis

Status:

Version 1

Abstract

Figures

Introduction

Methods

Results

Discussion

Declarations

References

Tables

Additional Declarations

Supplementary Files

Status:

Version 1