Serum creatinine levels are lower in diabetic patients compared with their nondiabetic counterparts. Therefore, estimated glomerular filtration rate (eGFR) is higher in the former than in the latter group. Factors associated with overestimation of renal function in diabetic patients were examined, and new formulae reflecting precise eGFR were created.
Eighty subjects (age 56.5 ± 15.4 years; 35 males [43.8%]; 40 patients with diabetes and 40 nondiabetic subjects) were enrolled. GFR was evaluated by inulin clearance (Cin). eGFR values were calculated based on serum creatinine and/or serum cystatin C levels. The factors related to the dissociation between eGFR and Cin in diabetic patients and the agreement among each of three eGFR and Cin were compared.
Although Cin was not significantly different between the diabetic and nondiabetic subjects (P = 0.2866), each of three eGFR measures from the diabetic patients was significantly higher than that of the nondiabetic subjects (P < 0.01). There were significant and positive correlations between the ratio of each eGFR/Cin, hemoglobin A1c, and glycated albumin. The intraclass correlation coefficients in diabetic patients were weaker than those in the nondiabetic subjects, and the intercepts of the regression lines between each eGFR measure and Cin in the diabetic patients were significantly higher than those of the nondiabetic subjects. New formulae for the calculation of eGFR corrected by the glycemic control indices were better than the original eGFR, particularly in diabetic patients.
eGFR overestimates Cin as glycemic controls worsen. eGFR corrected by hemoglobin A1c is considered to be clinically useful and feasible.
Introduction
With the continuing increase in the number of patients with end-stage renal disease, diabetic nephropathy, particularly resulting from type 2 diabetes mellitus (DM), has become the single most common cause of end-stage renal disease in Japan (1), as well as in the U.S. and Europe (2,3). Therefore, early diagnosis and strict glycemic control has emerged as a key issue to prevent the development of diabetic renal disease or at least to slow down the disease process (4,5). However, some diabetic patients often need to undergo renal replacement therapies at a relatively higher glomerular filtration rate (GFR) level, even before meeting the criteria for dialysis (6).
In addition to microalbuminuria, evaluation of the GFR by a creatinine-based formula is widely believed to provide a clinically useful index for the assessment of the progression of renal dysfunction and cardiovascular risk in diabetic nephropathy (7,8). Although serum creatinine is widely used to estimate GFR, serum creatinine is affected by muscle mass, sex, and age. Serum creatinine levels are significantly lower in the diabetic patients (9,10). Low serum creatinine levels were associated with type 2 DM in a recent study of nonobese middle-aged Japanese men (11), leading to the hypothesis that low creatinine might reflect low muscle mass volume (12,13). Low serum creatinine may also be related to significantly higher estimated GFR (eGFR) when calculated based on serum creatinine. Low serum creatinine levels in patients with diabetic nephropathy who require dialysis treatment are likely due to their low muscle volume (13). We previously demonstrated a significant decrease of serum creatinine in anuric diabetic hemodialysis patients compared with their nondiabetic counterparts (9), which was due to a reduction in the amount of creatine in muscle mass in the diabetic hemodialysis patients (9). Alternatively, it is increasingly recognized that the fraction of creatinine excreted from the renal tubuli that is not attributable to glomerular creatinine filtration might lead to overestimation of renal function when assessed by serum creatinine, particularly when renal function has deteriorated (10).
It is known that cystatin C is filtered by the glomeruli and metabolized in proximal tubular cells by Megalin in a Ca2+-dependent manner (14). Cystatin C is not affected by muscle mass, sex, or age. Accordingly, serum cystatin C levels and equations that estimate GFR based on serum cystatin C have recently been proposed as better markers of GFR (15–18).
The gold standard in the determining GFR is the measurement of inulin clearance (Cin). To date, no data exist regarding the relationship between glycemic control and Cin or eGFR in type 2 diabetic patients. Thus, we evaluated the correlation among three eGFR measures (i.e., eGFR based on serum creatinine [eGFRcr], serum cystatin C [eGFRcys], and both serum creatinine and cystatin C [eGFRcr-cys]) and Cin in both diabetic patients and nondiabetic subjects. We developed new formulae for the calculation of eGFR corrected by the glycemic control indices to more precisely assess renal function in diabetic nephropathy, after elucidating the degree to which renal function is overestimated in diabetic patients.
Research Design and Methods
Subjects and Assays
The study protocol was approved by the Ethics Committee of Osaka City University Graduate School of Medicine (#1444), Osaka, Japan, and the study was performed between January 2009 and March 2013. Written informed consent was obtained from each patient. The design was a single-center, randomized study that was conducted at Osaka City University Hospital. Eighty subjects (age 56.6 ± 15.4 years; 35 males and 45 females; 40 diabetic patients and 40 nondiabetic subjects) were enrolled. The diagnosis of DM was based on a history of diabetes or criteria according to the Report of the Expert Committee on the Diagnosis and Classification of Diabetes Mellitus (19). The patients enrolled were restricted to those with chronic kidney disease (CKD) stage 1–5 and were not on dialysis. Patients who were taking drugs that might affect tubular creatinine secretion and who had thyroid dysfunction that might affect cystatin C metabolism were excluded from the study.
Glycated albumin (GA) was measured by an enzymatic method using the Lucica GA-L kit (Asahi Kasei Pharma Corp., Tokyo, Japan) (20). GA was calculated as the percentage of GA relative to total albumin, which was measured in the same serum sample using a new bromocresol purple method (20). Serum cystatin C was measured by a colloidal gold immunoassay (Alfresa Pharma) that was traceable to ERM-DA471/IFCC or calibrated to the previously reported colloidal gold immunoassay (15).
Assessment of Renal Function by eGFR Based on Serum Creatinine, Serum Cystatin C, and the Combination of Serum Creatinine–Cystatin C
To develop an accurate eGFR estimation equation for the Japanese population, the Japanese Society of Nephrology launched a specific formula of eGFR based on serum creatinine, with the simultaneous measurement of Cin in 763 Japanese subjects, which was obtained by nationwide cooperation of nephrologists from December 2006 to July 2007, as described below (21):
Recently, the Japanese Society of Nephrology developed GFR-estimating equations based on serum cystatin C (CKD-EPIcys, eGFRcys) (22) and serum creatinine plus serum cystatin C (CKD-EPIcr-cys, eGFRcr-cys) (15).
In the current study, cystatin C was measured in 70 patients, and eGFRcys and eGFRcr-cys were calculated.
Measurements of Cin
Cin was determined by the constant input clearance technique with inulin. According to the method by Horio et al. (23), continuous intravenous infusion of inulin from a forearm antecubital vein was performed in the morning under a fasting state (i.e., a simple method of determining Cin by a single urine collection). In brief, the patients received 500 mL of water orally 15 min before the infusion. A 1% inulin solution (weight/volume) in saline was infused at 300 mL/h for the first 30 min and at 100 mL/min for the following 90 min. Patients completely voided their bladders at 45 min. Blood samples for the measurement of plasma inulin were collected at the same time. To maintain hydration, 180 mL of water was provided orally at the time of voiding the bladder. Blood and urine samples were taken at the end of the clearance period to measure plasma and urine inulin, respectively. A urine collection period of 90 min was set in order to increase the accuracy of the clearance study. Cin was calculated by the UinV/Pin method, where Uin represents the urinary inulin concentration, V the urinary volume, and Pin the plasma inulin concentration. Pin was the mean value of the plasma inulin concentration at the beginning and end of the clearance period. Pin was determined colorimetric ally using the N-1 naphthylethylenediamine and the anthrone method with a Corning 258 spectrophotometer (Corning) (24). In 59 patients, consisting of 23 diabetic patients and 36 nondiabetic subjects, creatinine clearance (Ccr) was also measured simultaneously during Cin measurement.
Statistical Methods
The data are expressed as the mean ± SD. Correlations between two variables were examined using simple regression analysis. Multiple regression analyses were performed to evaluate the relationships between eGFR/Cin and clinical parameters. The differences of the intercepts of the regression lines between Cin and eGFR were evaluated by ANCOVA. All analyses were performed using StatView 5 (SAS Institute Inc., Cary, NC) on a Windows computer (Microsoft, Redmond, WA). The level of statistical significance was set at P < 0.05. The agreement between Cin and each eGFR measure with and without correction by hemoglobin A1c was evaluated using intraclass correlation coefficient (ICC).
Results
Clinical Characteristics and Renal Function
The baseline characteristics of the 80 subjects enrolled in the current study are shown in Table 1. The average age of all subjects was 56.6 ± 15.4 years; 35 (43.8%) were male. Mean serum creatinine, serum cystatin C, eGFRcr, eGFRcys, eGFRcr-cys, and Cin were 1.0 ± 0.5 mg/dL, 1.2 ± 0.7 mg/dL, 64.1 ± 23.7 mL/min/1.73 m2, 68.9 ± 26.8 mL/min/1.73 m2, 67.7 ± 27.3 mL/min/1.73 m2, and 65.3 ± 23.8 mL/min/1.73 m2, respectively. The diabetic patients were significantly older than the nondiabetic subjects. There was no significant difference in Cin between the diabetic patients and nondiabetic subjects (P = 0.2866). Serum creatinine levels in the diabetic patients were significantly lower than the nondiabetic subjects (P < 0.0001). Serum cystatin C in the diabetic patients was significantly lower than the nondiabetic subjects (P < 0.0001). eGFRcr, eGFRcys, and eGFRcr-cys in the diabetic patients were all significantly higher than those in the nondiabetic subjects. These data indicated that the assessment of renal function on the basis of serum creatinine and cystatin C might lead to overestimation of renal function in diabetic patients.
Relationship Between Each eGFR Measure (eGFRcr, eGFRcys, and eGFRcr-cys) and Cin in Diabetic Patients and Nondiabetic Subjects
As shown in Fig. 1, there was a significant and positive correlation between eGFRcr and Cin in both diabetic patients (r = 0.683; P < 0.0001) and nondiabetic subjects (r = 0.930; P < 0.0001). There was a significant and positive correlation between eGFRcys and Cin in both diabetic patients (r = 0.584; P < 0.0001) and nondiabetic subjects (r = 0.845; P < 0.0001). There was also a significant and positive correlation between eGFRcr-cys and Cin in both diabetic patients (r = 0.712; P < 0.0001) and nondiabetic subjects (r = 0.930; P < 0.0001).
In each eGFR measure of eGFRcr, eGFRcys, and eGFRcr-cys, the correlation coefficients for the nondiabetic subjects were greater than those of the diabetic patients, and the intercepts of the regression line for diabetic patients were significantly higher than those of the nondiabetic subjects (eGFRcr, P < 0.0001; eGFRcys, P = 0.0027; and eGFRcre-cys, P < 0.0001, respectively). These findings indicated significant overestimation of renal function in diabetic patients on the basis of each of the three eGFRs measurements.
Correlation of the eGFRcr/Cin, eGFRcys/Cin, and eGFRcr-cys/Cin With Clinical Parameters
We next evaluated the factors associated with the overestimation of renal function in diabetic patients. Since the dissociation of the regression line between each eGFR measure and Cin could be represented by the eGFR/Cin ratio, we first examined a simple regression analysis of each eGFR/Cin ratio with the various clinical parameters. Although there were no significant correlations between each eGFR/Cin ratio and age, sex, BMI, or blood pressure, there were significant and positive correlations between each eGFR/Cin ratio and the glycemic control indices of hemoglobin A1c (eGFRcr/Cin: r = 0.605, P < 0.0001; eGFRcys/Cin: r = 0.340, P = 0.0042; and eGFRcr-cys/Cin: r = 0.603, P < 0.0001) and GA (eGFRcr/Cin: r = 0.565, P < 0.0001; eGFRcys/Cin: r = 0.372, P = 0.0092; and eGFRcr-cys/Cin: r = 0.548, P < 0.0001) (Fig. 2).
Multiple regression analyses were performed to examine whether the glycemic control indices of hemoglobin A1c and GA were associated with eGFR/Cin ratio of eGFRcr/Cin, eGFRcys/Cin, and eGFRcr-cys/Cin after adjustment for age, sex, BMI, and blood pressure. Glycemic control indices of hemoglobin A1c and GA were significantly and positively associated with each eGFR/Cin ratio of eGFRcr/Cin (n = 80), eGFRcys/Cin (n = 70), and eGFRcr-cys/Cin (n = 70) after the adjustment.
Relationship Between the Secretory Component of Ccr and Glycemic Control Indices
We next evaluated the mechanism of overestimation of the renal function as evaluated by eGFRcr in 59 subjects (23 diabetic patients and 36 nondiabetic subjects), in whom Ccr was measured concomitant with the Cin measurements. To quantify the secretory component of Ccr (SFcr), SFcr was calculated as follows (25):
Relationships between the SFcr and glycemic control indices of hemoglobin A1c and GA were examined. There were significant positive correlations between SFcr and hemoglobin A1c (r = 0.359; P = 0.0053) and between SFcr and GA (r = 0.536; P = 0.0009).
The Formulae for eGFR Corrected by the Glycemic Control Index
From the regression line of Fig. 2, the formulae for each eGFR measure corrected by the glycemic control indices of hemoglobin A1c and GA were considered. A linear function, y = ax + b, was applied to elucidate eGFR corrected by each of the glycemic control indices of HbA1c and GA, where y = each eGFR/Cin, x = HbA1c or GA, a = slope, and b = intercept.
Based on the results of Fig. 2, when using HbA1c, eGFR/Cin=slope × HbA1c + intercept, then, Cin = eGFR/(slope × HbA1c + intercept); the calculated Cin was then considered to be the eGFR corrected by HbA1c.
Based on the results of Fig. 2; when using GA, eGFR/Cin = slope × GA + intercept; then, Cin = eGFR/(slope × HbA1c + intercept); the calculated Cin was then considered to be the eGFR corrected by GA.
Using the above calculation, the following formulae were considered for each eGFR measure corrected by the glycemic control indices of HbA1c and GA.
Cin, eGFRcr, eGFRcys, eGFRcr-cys, and Each eGFR Measure Corrected by Hemoglobin A1c and GA in Diabetic Patients and Nondiabetic Subjects
As shown in Table 1, although there were no significant differences in the actual Cin between the diabetic patients and nondiabetic subjects, there were significant differences in the original eGFRcr (P < 0.0001). However, as expected, eGFRcr corrected by the glycemic control indices of hemoglobin A1c and GA did not differ significantly between the diabetic patients and the nondiabetic subjects (P = 0.1841 and P = 0.2493, respectively). Furthermore, although the original eGFRcys was significantly different between the diabetic patients and the nondiabetic subjects (P = 0.0017), eGFRcys corrected by the glycemic control indices of hemoglobin A1c and GA did not differ between groups (P = 0.1393 and P = 0.2571, respectively). Finally, although the original eGFRcr-cys was significantly different between the diabetic patients and the nondiabetic subjects (P < 0.0001), eGFRcr-cys corrected by the glycemic control indices of hemoglobin A1c and GA did not differ between groups (P = 0.0562 and P = 0.1149, respectively).
Relationship and Agreement Between Each eGFR Measure (eGFRcr, eGFRcys, and eGFRcr-cys) and Each eGFR Measure Corrected by the Glycemic Control Indices and Cin
Correlations and ICCs between Cin, each of the uncorrected eGFR measures (eGFRcr, eGFRcys, and eGFRcr-cys) and each eGFR measure corrected by hemoglobin A1c were examined in order to validate the suitability of the correction. As shown in Table 2, the correlation coefficients and ICCs between Cin and each eGFR measure corrected by hemoglobin A1c were stronger than those between Cin and the original eGFR measures in all subjects. The correlation coefficients and ICCs between Cin and each eGFR measure corrected by hemoglobin A1c became stronger than those between Cin and the original eGFR measures in the diabetic patients. Similarly, the correlation coefficients and ICCs between Cin and each eGFR measure corrected by GA became stronger than those between Cin and the original each eGFR measure in all subjects, particularly in diabetic patients (data not shown).
CONCLUSIONS
In the current study, we demonstrated that the estimation of renal function by eGFRcr, eGFRcys, and eGFRcr-cys led to overestimation of renal function in diabetic patients. Our study also showed that the correlation coefficients and ICCs of each of three eGFR measures and Cin in diabetic patients were lower than those of nondiabetic subjects and that the intercepts of the regression lines in diabetic patients were significantly higher than those in nondiabetic subjects, indicating that each of three eGFR measurements in diabetic patients are more inaccurate than in the nondiabetic subjects. Further, the values of each of the three eGFR/Cin ratios (eGFRcr/Cin, eGFRcys/Cin, and eGFRcr-cys/Cin) correlated significantly with the glycemic control indices of hemoglobin A1c and GA, suggesting that the three values of eGFR estimation (i.e., eGFRcr, eGFRcys, and eGFRcr-cys) overestimate renal function as glycemic controls worsened. Thus, the apparent increase in each of three eGFR values relative to Cin might be associated with poor glycemic condition in diabetic patients.
In general, Ccr overestimates Cin (10,25–27). This phenomenon has been reported to be caused by creatinine secretion from the renal tubuli (25,28), in addition to glomerular filtration of creatinine. To overcome the weakness of Ccr compared with Cin, the Modification of Diet in Renal Disease (MDRD) formula of eGFR was developed in 2006, in consideration of serum creatinine, age, sex, and race (29,30). The eGFR equation based on serum creatinine, age, and sex was also constructed by the Japanese Society of Nephrology by directly measuring Cin in 763 Japanese subjects (21). Furthermore, the Chronic Kidney Disease Epidemiology (CKD-EPI) formula was recently developed (15). However, since serum creatinine is affected by muscle mass, sex, and age and since the serum creatinine levels are significantly lower in diabetic patients (9,10), cystatin C has recently been proposed as a better marker of renal function. Cystatin C is produced by the nucleated cells of the body and acts as a circulating cysteine proteinase inhibitor (31). It is known that cystatin C is filtered by the glomeruli and is not affected by muscle mass, sex, and age, as is observed with serum creatinine levels (18). Serum cystatin C levels and an eGFR equation based on serum cystatin C have been recently proposed as more relevant markers of renal function (16–18). Horio et al. (15) recently reported eGFR equations by using serum cystatin C and eGFR by using the combination of serum creatinine and cystatin C. Recently, Inker et al. (32) reported that the equation with the combination of serum creatinine and cystatin C performed better than that with either of these markers alone. Consistent with the findings of Inker et al. (32), in our study, the ICC of eGFRcr-cys was stronger than that of eGFRcr and eGFRcys in the nondiabetic subjects, when compared with Cin.
However, in the current study, eGFRcr, as well as eGFRcys and eGFRcr-cys, were higher in diabetic patients compared with the nondiabetic subjects. This result suggests that sustained elevations of plasma glucose might lead to increased eGFR in diabetic patients. Accordingly, in the current study, eGFR/Cin ratios were examined in relation to glycemic control indices of hemoglobin A1c and GA. We found that there were significant and positive correlations between each of three eGFR/Cin ratios and hemoglobin A1c and between each of three eGFR/Cin ratios and GA, as shown in Fig. 2. Serum creatinine is excreted into urine via glomerular filtration and tubular secretion (25). It is considered that Ccr represents the net effect of glomerular creatinine clearance and tubular creatinine secretion (25). Cin is considered to be equivalent to glomerular creatinine clearance, since inulin is not secreted from the renal tubuli. Accordingly, the formula of the secretory component of Ccr (SFcr: [Ccr − Cin]/Ccr) is considered to be a ratio of tubular creatinine secretion (25). The secretory components of Ccr correlated significantly and positively with glycemic control indices of hemoglobin A1c and GA, suggesting that increased secretion of creatinine was related to the hyperglycemic condition. Organic cation transporters, which have been reported to be associated with creatinine secretion in renal tubuli (33), are modulated by high glucose in rats via the increased oxidative stress of advanced glycation end products (34). The results of the current study may reflect these experimental findings. Thus, from the findings of the current study, we consider that the overestimation of renal function in eGFR based on serum creatinine in diabetic patients may be caused by increased tubular secretion of creatinine.
It should be noted that the apparent overestimation of renal function on the basis of higher eGFR, not only eGFRcr but also eGFRcys and eGFRcr-cys, relative to Cin might confound the precise assessment of renal function and may thus delay the initiation of treatment when renal function declines in diabetic patients, particularly in those with poorer glycemic control. In this study, we developed new formulae that were corrected by the glycemic control indices. Although eGFRcr, eGFRcys, and eGFRcr-cys in the diabetic patients were significantly higher than those in nondiabetic subjects, the eGFR measures after correction by hemoglobin A1c and GA were not significantly different between the diabetic patients and the nondiabetic subjects, which was consistent with the result of Cin, which was not significantly different between diabetic patients and nondiabetic subjects (Table 1). With correction of each eGFR measure by hemoglobin A1c, the correlation coefficients and the ICCs between Cin and each eGFR measure corrected by hemoglobin A1c became stronger than those between Cin and the uncorrected eGFR measure in all subjects. The correlation coefficients and the ICCs became even stronger in the diabetic patients (Table 2). These results indicate that, in order to accurately evaluate renal function, eGFR should be corrected by the glycemic control indices in diabetic patients, particularly in those with poorer glycemic control. Since serum creatinine and hemoglobin A1c are among the laboratory data that we generally measure in clinical practice, compared with serum cystatin C and GA, we consider that eGFRcr corrected by hemoglobin A1c is the most feasible and useful in the evaluation of renal function, particularly in diabetic patients. From the current study, the formula, eGFRcr/(0.428 + 0.085 × HbA1c), could be used for better and more accurate evaluation of renal function of diabetic patients in order to accurately evaluate renal function.
This study has some limitations. Firstly, the study was performed in a relatively small number of patients, and a further large-scale study with a greater number of patients is needed to confirm the clinical validity of eGFR correction by hemoglobin A1c. Rognant et al. (35) evaluated the three equations for GFR (i.e., the Cockcroft and Gault [CG] equation, MDRD equations, and CKD-EPI collaboration equation) in 246 diabetic patients in six institutions, and they found that the simplified MDRD formula for estimating GFR performed better than the CKD-EPI and CG equations. In the current study, we evaluated renal function in 40 diabetic patients and 40 nondiabetic subjects by Cin in a single institution, which could prevent the variation of the values among several institution. We compared three eGFR equations, based on serum creatinine and serum cystatin C (eGFRcr, eGFRcys, and eGFRcr-cys), and also explored the factors that affected dissociation of eGFR and Cin. We found that glycemic control indices were significantly associated with the dissociations between Cin and three eGFR values. Furthermore, we constructed new equations, which used hemoglobin A1c or GA, and showed that these equations using hemoglobin A1c or GA were better than the original eGFR equations, particularly in diabetic patients. Although we created additional formulae for correction by hemoglobin A1c and GA that were only applied to diabetic patients (data not shown), these measures did not differ substantially from those in all patients and subjects. Secondly, although the higher eGFRcr in diabetic patients was suggested to be caused by tubular creatinine secretion, we could not determine the mechanism of overestimation of the renal function in the diabetic patients when it was estimated by eGFRcys, since we were not able to measure urine cystatin C. However, it has been reported that urinary cystatin C appears even before the increase in the classical biomarkers of diabetic nephropathy, such as albuminuria and urinary protein (36). It has also been reported that cystatin C is filtered by glomeruli and metabolized in proximal tubular cells by binding to megalin in a Ca2+-dependent manner (14). Ogasawara et al. (37) reported that urinary megalin was increased in correlation with the severity of type 2 diabetic nephropathy, likely leading to reduced metabolism of cystatin C in tubuli. These reports suggest that increases in urinary cystatin C excretion in diabetic patients may be one of the mechanisms involved in the dissociation of eGFRcys and Cin. Thirdly, this was a cross-sectional study. Further studies may be needed to explore whether the dissociation of eGFR and Cin could be reduced by improving glycemic control. However, this is the first study in which poor glycemic control has been considered to cause the dissociation of eGFR and Cin in diabetic patients. Finally, our study, which developed new formulae corrected by hemoglobin A1c or GA, examined Japanese subjects. Further studies are needed to construct formulae that adjust for glycemic control indices, in other races, and further large-scale studies with a greater number of patients are needed to confirm the clinical validity of eGFR correction by hemoglobin A1c.
In conclusion, in CKD patients, each of three eGFR measures (eGFRcr, eGFRcys, and eGFRcr-cys) overestimated renal function when glycemic control was poor, suggesting that eGFR is needed to be corrected by the glycemic control indices. The new formulae for the correction of eGFR by hemoglobin A1c are better than the original eGFR. We consider that eGFRcr corrected by hemoglobin A1c is clinically feasible and useful.
Article Information
Duality of Interest. No potential conflicts of interest relevant to this article were reported.
Author Contributions. A.T. and E.I. researched data, contributed to the discussion, and wrote and reviewed the manuscript. Y.O., M.Ic., S.N., Y.M., K.M., J.U., S.F., M.E., T.N., and M.In. contributed to the discussion and reviewed the manuscript. A.T. and E.I. are the guarantors of this work and, as such, had full access to all the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis.