Abstract
Background:
25-Hydroxyvitamin D [25(OH)D] deficiency has been implicated as a possible risk factor for the onset and progression of diabetes kidney disease (DKD). The aim of this study was to evaluate the interaction between levels of 25(OH)D and DKD in type 2 diabetes mellitus (DM) patients.
Methods:
Cross-sectional design, outpatient type 2 DM. Glomerular filtration rate (GFR) was measured by 51Cr-EDTA and estimated by Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI), urinary albumin excretion (UAE) by immunoturbidimetry, and 25(OH)D by chemiluminescence. Receiver operating characteristic (ROC) curve analysis and generalized linear model (Poisson robust regression estimator) were used to assess the interaction between 25(OH)D levels and renal function.
Results:
A total of 114 type 2 DM patients aged 60 ± 10 years, 49 males (43%), DM duration 22 ± 10 years, with GFR > 60 ml/min/1.73 m2 were evaluated. Patients with GFRs 60–90 (n = 50) had significantly lower 25(OH)D levels than individuals with GFRs > 90 ml/min/1.73 m2 (n = 64), respectively 40 ± 20 versus 48 ± 20 nmol/l, p = 0.027. This difference was more pronounced for older individuals (39 ± 20 versus 54 ± 23 nmol/l, respectively), and Poisson robust regression disclosed that lower 25(OH)D [Poisson regression (PR) = 0.989, confidence interval (CI): 0.978–0.999, p = 0.034], and advanced age (PR = 1.050, CI: 1.007–1.096, p = 0.023) were significantly associated with the lower GFR category, adjusted for seasons. ROC curve analysis showed that the cutoff point of 25(OH)D of 41 nmol/l was associated with lower GFR [area under the curve (AUC) = 0.694, p = 0.009]. CKD-EPI estimated GFR (eGFR) was not associated with 25(OH)D in any analysis. There was no difference in 25(OH)D levels between patients with elevated UAE as compared with normoalbuminuric ones (44 ± 21 versus 46 ± 19 nmol/l, p = 0.587).
Conclusion:
Lower levels of 25(OH)D are associated with decreased GFR in patients with type 2 DM, especially in older patients, with no evidence of interaction with UAE levels.
Keywords
Introduction
A recent meta-analysis of data on diabetes mellitus (DM) patients revealed an annual incidence of elevated albuminuria of 2–8%; the incidence of estimated glomerular filtration rate (eGFR) < 60 ml/min/1.73 m2 was as high as 4% per year, and that of end-stage renal disease (ESRD) up to 2% per year. 1 Currently, diabetes kidney disease (DKD) is responsible for about 40% of chronic kidney disease (CKD) cases worldwide.2–4 These findings indicate that DKD is a noteworthy health problem all over the world, even though data from developing countries are notably scarce.
Current data suggest that cardiovascular and all-cause rates of mortality are declining in DM, so an increased prevalence of DKD is an upcoming trend.5,6 Sharp increases of 40–700% in the incidence of ESRD in diabetes have been reported between 2002 and 2015 in Russia, Mexico, Bosnia, Australia, and several Asian countries, also reflecting the increasing prevalence of diabetes itself.2,6 Whereas these increases have been reported for the overall population, diabetic patients have tended to show a decrease in rates of micro- and macrovascular complications over the last 20 years. Notwithstanding, the smallest reductions have been described for ESRD, in contrast with the marked reductions in acute myocardial infarction, stroke, and amputations. 3
Unfavorable DKD evolution stems from the fact that current treatment is not able to completely control the progression of the disease, since the pathogenesis of kidney damage has not yet been fully elucidated. 4 Vitamin D deficiency has raised special interest because some cross-sectional studies in DM patients have demonstrated an association between lower levels of 25(OH)D and elevated urinary albumin excretion (UAE).7–9 25(OH)D deficiency is highly prevalent among patients with type 2 diabetes, especially in the presence of DKD, when up to 50–70% of patients are affected.7,8,10–15 The direct interaction between decreased 1,25(OH)2D levels and renal damage is expected and well established, but it is rather intriguing when it comes to 25(OH)D. Either 25(OH)D deficiency would be causing kidney injury, or occurring as a consequence of renal involvement. Alternatively, confounding factors such as obesity, high blood pressure, and hyperglycemia – usually coexisting with lower 25(OH)D levels – would be the real promoters of kidney disease.4,16 However, not all studies confirm the association between 25(OH)D levels and albuminuria,17,18 and, as a matter of fact, no previous study has properly investigated measured GFR in this scenario. Therefore, the aim of this study was to investigate the interaction between 25(OH)D levels and the degree of renal function impairment, analyzing urinary albumin excretion (UAE) and GFR measured with a reference method.
Subjects and methods
Subjects
This cross-sectional study included type 2 DM outpatients with preserved renal function (GFR > 60 ml/min/1.73 m2) from the Endocrine Division of the South University Hospital. Informed written consent was obtained from all patients, as approved by Grupo de Pesquisa e Pós-Graduação (GPPG), project number 140501, and ethical approval was obtained from the Research Ethics Committee of Hospital de Clínicas de Porto Alegre (Comitê de Ética e Pesquisa: CEP). It complied with the Declaration of Helsinki. Type 2 DM was defined by onset after 35 years of age, no insulin therapy in the first 5 years since diagnosis, no previous episodes of ketoacidosis, and compatible C-peptide levels. 19 Skin color was self-defined by the patient. No patient was taking vitamin D supplements.
The patients underwent a complete physical examination. Body mass index (BMI) was determined as the weight in kilograms divided by height in meters squared. Waist circumference was measured along a line midway between the rib cage and iliac crest. Blood pressure was taken twice, after 5 min of sitting at rest, with a validated digital sphygmomanometer (OMRON, Model HEM-705CP); the mean of two measurements was used.
Methods
Serum samples were collected after 12-h fasting; glucose was determined by the glucose oxidase method, cholesterol by an enzymatic method, and HbA1c was analyzed with certified HPLC (Merck-Hitachi L-9100 Analyzer, Merck, Darmstadt, Germany).
GFR was measured by the 51Cr-EDTA single-injection method, after a single intravenous dose of 150 µCi, and calculated as: volume of distribution × 0.693 × 0.87 × 1000/t1/2, where 0.87 is a correction factor developed by Chantler. 20 Blood samples were collected at 0, 2, 3, and 4 h after the injection, and the results were expressed in ml/min/1.73 m2. The coefficient of variation established in our institution is 12%. 21 GFR was also estimated by the CKD-EPI (Chronic Kidney Disease Epidemiology Collaboration) equation and expressed as eGFR (ml/min/1.73 m²), using a traceable Jaffe creatinine assay. 22 Patients were categorized according to 51Cr-EDTA GFR (51Cr GFR) as ⩾90 or 60–90 ml/min/1.73 m2, according to Kidney Disease Improving Global Outcomes (KDIGO). 4
UAE was measured by immunoturbidimetry in random urine samples, and values above 14 mg/l were considered elevated. 23 DKD was defined by elevated UAE and/or decreased GFR.
25(OH)D was measured by chemiluminescence, which detects levels of vitamin D between 10 and 375 nmol/l. The inter-assay coefficients of variation (CV) were 7% and 6.3% for values of 45 and 92.5 nmol/l, respectively, with a 5% intra-assay CV. Vitamin D deficiency was defined by levels of 25(OH)D < 50 nmol/l. 24
Statistical analysis
The results were expressed as mean ± standard deviation (SD), percentages or median [interquartile range (IQR)]. Continuous parametric variables were assessed using Student’s t test or ANOVA; variables without normal distribution were analyzed using non-parametric tests or Mann–Whitney test; qualitative variables were analyzed using the chi-square test. Receiver operating characteristic (ROC) curves were used to assess the association between GFR levels and 25(OH)D cutoff points; Youden’s test was used to check the equilibrium point. Poisson analysis was performed to assess the possible association of categorized binary GFR with 25(OH)D, adjusted for season, gender, and HbA1c.
To evaluate the relation between elevated albuminuria and 25(OH)D levels, the sample size was calculated on the website www.openepi.com; considering an alpha error of 5% and a power of 80%, 50 individuals were required in each group to find a putative difference. 25 As for GFR and 25(OH)D interaction, since no studies were available in the literature, the power of our findings, and also the Cohen effect size, were calculated to analyze the impact of 25(OH)D on renal function. 26
A p value < 0.05 was considered statistically significant. Statistical analyzes were performed using SPSS software, version 25 (SPSS, Chicago, IL, USA).
Results
This study included 114 type 2 DM patients, with a mean age of 60 ± 10 years, 49 males (43%), mean DM duration of 22 ± 10 years, and BMI of 30 ± 4 kg/m2.
Table 1 presents the clinical and laboratory characteristics of patients classified according to GFR levels: >90 versus 60–90 ml/min/1.73 m2. Patients with GFR 60–90 were older and displayed significantly lower 25(OH)D values compared with those with GFR > 90 ml/min/1.73 m2 (Figure 1A). In contrast to 51Cr-EDTA measured GFR, estimated CKD-EPI GFR was not able to recognize any association with lower 25(OH)D values, since, when the patients were divided by estimated CKD-EPI GFR values of >90 versus 60–90 ml/min/1.73 m2, vitamin D levels were similar between the groups thus classified (45 ± 18 versus 45 ± 22 nmol/l, p = 0.911). Even when analyzed according to each stage of estimated GFR in older patients, it was also not possible to see any association, as opposed to measured GFR. Furthermore, there was a positive significant correlation between 51Cr-EDTA GFR and 25(OH)D (r = 0.22, p = 0.021), not observed for estimated GFR (r = 0.09, p = 0.356).
Clinical and laboratory characteristics of patients with type 2 DM according to GFR (51Cr-EDTA GFR) categorized as >90 versus 60–90 ml/min/1.73 m2.
Data expressed as mean ± SD, percent, median (IQR).
χ2 test was used for qualitative variables and Student’s t test and Mann–Whitney test were used for quantitative variables. Statistical significance was assumed at p < 0.05.
ACEi, angiotensin converting enzyme inhibitor; ARB, angiotensin receptor blocker; BMI, body mass index; CKD-EPI, Chronic Kidney Disease Epidemiology Collaboration; 25(OH)D, hydroxyvitamin D; DM, diabetes mellitus; eGFR, estimated glomerular filtration rate; GFR, glomerular filtration rate; HbA1c, glycated hemoglobin; IQR, interquartile range; SD, standard deviation; UAE, urinary albumin excretion.

(A) 25(OH)D values according to GFR stages, and (B) UAE values. Student’s t test was used in (A) and Mann–Whitney test was used in (B).
Table 2 shows the results of patients divided by UAE levels, categorized as normo (<14 mg/l) versus micro/macroalbuminuria (⩾14 mg/l). The values of 25(OH)D were similar between patients with normal versus elevated UAE (Figure 1B). Likewise, the proportion of vitamin D deficiency was similar between the groups.
Clinical and laboratory characteristics of patients with type 2 DM according to UAE categorized as normo (<14 mg/l) versus micro/macroalbuminuria (⩾14 mg/l).
Data expressed as mean ± SD, percent, median (IQR).
χ2 test was used for qualitative variables and Student’s t test and Mann–Whitney test were used for quantitative variables. Statistical significance was assumed at p < 0.05.
ACEi, angiotensin converting enzyme inhibitor; ARB, angiotensin receptor blocker; BMI, body mass index; 25(OH)D, hydroxyvitamin D; GFR, glomerular filtration rate; HbA1c, glycated hemoglobin; IQR, interquartile range; SD, standard deviation; UAE, urinary albumin excretion.
Since age seemed to influence the association of GFR with 25(OH)D levels, we categorized the patients according to median age, which was 60 years. Figure 2 shows the levels of 25(OH)D according to GFR in patients classified by age ⩾ 60 (n = 62) versus <60 years (n = 52). Interestingly, 25(OH)D levels were lower in patients with lower GFRs, but only in patients ⩾60 years (39 ± 20; n = 34 versus 54 ± 23 nmol/l; n = 28), in contrast to patients <60 years (43 ± 17; n = 16 versus 44 ± 17 nmol/l; n = 36, in each GFR stage). The calculated power for this finding was 77%. In the same way, 25(OH)D was positively and moderately correlated with GFR only in older patients (r = 0.36, p = 0.004), as opposed to younger patients (r = 0.12, p = 0.396). Furthermore, Poisson robust analyses taking binary GFR category as the outcome, demonstrated that, in older subjects, lower 25(OH)D [ Poisson regression (PR) = 0.989, confidence interval (CI): 0.978–0.999, p = 0.034], and advanced age (PR = 1.050, CI: 1.007–1.096, p = 0.023) were significantly associated with the lower GFR category, adjusted for seasons. This association was not observed in younger patients [PR for 25(OH)D was 0.994, CI: 0.969–1.020, p = 0.628], age: PR 1.016, CI: 0.956–1.080, and season of the year, p = 0.610. These results remained the same even after removing an outlier patient that had elevated levels of 25(OH)D.

25(OH)D values according to GFR values >90 versus 60–90 ml/min/1.73 m2 in patients aged <60 versus >60 years of age. ANOVA test followed by Tukey pairwise post hoc was used. *p = 0.004.
To allow a more meaningful interpretation of clinical relevance of 25(OH)D levels, we also calculated the effect size of Cohen. 26 The effect size of 25(OH)D levels on GFR was 0.4 for the whole group (small effect), and 0.69 for older group (intermediate effect), according to the interpretation table for different effect sizes. 26
To define the cutoff of 25(OH)D capable of identifying a link between the level of GFR and 25(OH)D levels, a ROC curve was constructed. By means of the Youden test, the value of 41 nmol/l of vitamin D (Figure 3, AUC = 0.694, p = 0.009) was determined as the cutoff point that identified the presence of a decline in measured GFR, with a calculated power of 78%.

The value of 41 nmol/l of 25(OH)D was the identified Youden cutoff point that was associated with the presence of a decline in measured GFR in patients >60 years; AUC = 0.694, p = 0.009, power 78%.
Discussion
In the present study, lower levels of 25(OH)D were associated with lower GFR in type 2 DM patients with preserved renal function; no evidence of interaction between 25(OH)D and increased UAE was found.
Regarding GFR, we were able to demonstrate that 25(OH)D levels were lower in patients with GFR values of 60–90 when compared with >90 ml/min/1.73 m2, primarily in older patients. As far as we know, our study was the first to use a reference GFR measurement method to demonstrate this association. Previous studies investigated only equation-based eGFR. When analyzing eGFR, some studies described a direct interaction between 25(OH)D and GFR,15,16 whereas other authors did not confirm this finding.7,8,12,15,18,27 As we and others have previously demonstrated, the accuracy of GFR estimated by equations is poor, especially in patients with diabetes, markedly underestimating true GFR.21,28–30 This finding is probably due to either analytical factors, such as interference of glucose levels on creatinine measurement, or to the lack of sensitivity of creatinine to identify glomerular hyperfiltration, a peculiar scenario in diabetes. 21 Thus, the studies that determined eGFR could have failed to demonstrate such an association because of eGFR limitations. Accordingly, when GFR was evaluated by CKD-EPI equation in our study, no association between 25(OH)D and renal function was found, in contrast to the direct correlation obtained with measured 51Cr-EDTA GFR. Nevertheless, this significant correlation was limited herein to older patients, probably because age reduces renal functioning mass, making the kidney more susceptible to damage. Several factors, from genetic predisposition to environmental factors, create a condition where the renal status of elderly individuals shows different susceptibility to injury. Renal aging is a complex process expressed by anatomical and functional changes, such as impaired function of glomerular, tubular, and endocrine systems, such as 25(OH)D activation. 31 Therefore, this reduced activation, on top of low levels of cholecalciferol, might explain the increased vulnerability of elderly kidneys. Through a presumed increased state of chronic inflammation, 25(OH)D deficiency would represent an additional source of aggression to the aging kidney. In a cohort of elderly patients with various degrees of renal impairment, reduced 25(OH)D levels were independent predictors of progression to dialysis and death. 32 A recent investigation showed that the kidney-secreted hormone Klotho would be a central player at the ageing-inflammation interface. 33
Regarding albuminuria, some studies have demonstrated that patients with elevated UAE have lower 25(OH)D levels compared with normoalbuminuric individuals.34,35 However, herein, no association between UAE and 25(OH)D was found. This divergence could be attributed to the fact that about 70% of our patients were using angiotensin converting enzyme inhibitor (ACEi) or angiotensin receptor blocker (ARB), thereby masking increased UAE levels. Another explanation could be that urinary albumin levels in our cohort were not as high as those of previous studies; an inverse correlation has been previously demonstrated between levels of albuminuria and 25(OH)D.7,12,35 The importance of this putative association is still a matter of debate; whereas cross-sectional and longitudinal studies suggest an unfavorable effect of vitamin D deficiency on albuminuria, studies of vitamin D supplementation show controversial effects on UAE.34,36,37 Further clarification could come from ongoing larger randomized studies. 38
Typical characteristics of type 2 DM, such as obesity, high blood pressure and hyperglycemia, usually occur concomitantly with lower 25(OH)D levels.16,18 These parameters could be the real factors promoting kidney damage. Indeed, HbA1c has been found to be inversely related to serum 25(OH)D levels. 39 However, whereas two recent meta-analyses showed that vitamin D supplementation improved glycemic control in vitamin D-deficient type 2 diabetes patients, two other studies showed conflicting results.40,41 In our study, there was a weak inverse correlation between HbA1c and 25(OH)D; nevertheless, this was not evident in multivariable analysis. The same dispute applies to blood pressure levels. A number of studies have demonstrated that both systolic and diastolic pressures are increased in patients with vitamin D deficiency,16,42 and a meta-analysis described that supplementation was able to reduce blood pressure in patients with diabetes. 43 This question is still open to debate, and the independent roles of vitamin D deficiency and elevated blood pressure on kidney damage are still to be elucidated.
An important issue to be considered in our study is the season of the year in which the samples for vitamin D analyses were collected. The seasonal variation of ultra-violet B (UVB) solar radiation causes a variation in circulating levels of 25(OH)D. In a cohort of 4758 individuals inhabitants of northern Finland, individuals whose blood samples were collected during periods of low sunlight exposure (winter months) showed an increased risk of showing low levels of 25(OH)D. 44 Likewise, vitamin D deficiency was observed in winter and spring as compared with summer and fall in a group of 11,150 Italian individuals, 45 and the same pattern has been demonstrated in a German population. 46 A retrospective analysis of 25(OH)D including 39,004 individuals in Brazil, aged 2–95 years, revealed that vitamin D deficiency was more evident in winter when compared with summer. 47 All these studies confirm the importance of seasonality as a potential bias in determining 25(OH)D levels. In uni and multivariable analyses, the season of collection was not related to renal function in our study.
The present study has strengths and limitations. First, the patients were examined with a GFR reference method, in contrast to all previous studies, which evaluated only eGFR, which is particularly inaccurate in diabetic patients.23,28,29 Indeed, in contrast to measured GFR, CKD-EPI eGFR was not able to recognize any association with 25(OH)D values in our study. Furthermore, the use of a reference GFR method allows a more precise analysis of patients with relatively preserved renal function, to reinforce that lower 25(OH)D levels were present even in early stages of kidney damage. 48 As for limitations, the cross-sectional nature of data does not allow the analysis of causality, and only speculations can be made in this sense. Another point is the relatively restricted size of the sample. Nonetheless, we carefully described our results showing a calculated power of about 80%, also describing details of Cohen effect size to improve data interpretation.
In conclusion, our study was the first to demonstrate an interrelationship between lower levels of 25(OH)D and lower GFR in older patients with type 2 DM, in a relatively preserved range of GFR function. No evidence of interaction was found with UAE, although the broad use of ACEi/ARB could have interfered with this analysis.
Footnotes
Author contribution(s)
Funding
The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by grants from Fundação de Amparo à Pesquisa do Estado do Rio Grande do Sul (FAPERGS), Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Programa de Apoio a Núcleos de Excelência (PRONEX process number 16/2551-0000-476-5), Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), Fundo de Incentivo à Pesquisa e Eventos (FIPE) Hospital de Clínicas de Porto Alegre.
Conflict of interest statement
The authors declare that there is no conflict of interest.
