Introduction

Distal sensorimotor polyneuropathy (DSPN) is often considered a late complication of diabetes with a prevalence of up to 50% in people with longer disease duration [1]. However, increasing evidence suggests that DSPN is also more prevalent in people with prediabetes compared to those with normal glucose tolerance [2, 3] and in individuals with obesity compared to normal-weight individuals [4, 5]. Given the high burden of DSPN [6] and the lack of effective therapeutic options, other modifiable risk factors should also be investigated.

Among these modifiable risk factors, vitamin D is determined by exposure to sunlight, dietary sources or supplementation [7]. Increasing evidence shows associations between low vitamin D levels, as assessed by serum 25-hydroxyvitamin D (25(OH)D), and diseases of the central nervous system [8]. In vitro and in vivo studies indicated that vitamin D is also an essential player in the peripheral nervous system [9] that regulates neural development and function by promoting neuronal cell differentiation [10], enhancing neurotrophin expression [11, 12] and improving myelinisation after injury [13]. Based on these preclinical findings, and given the link between low vitamin D levels and diabetes, particularly in older individuals [14], low vitamin D has been singled out as a potential risk factor in the pathogenesis of DSPN.

Epidemiological studies in individuals with diabetes showed an association of low vitamin D levels with both signs and symptoms of DSPN [15,16,17,18,19,20,21], impaired nerve conduction velocity [18, 19, 21,22,23,24,25] and neuropathic pain [19, 26]. A meta-analysis including 1,048 Caucasian individuals with diabetes reported an association between vitamin D deficiency and prevalence of DSPN [27]. These findings were confirmed in Asian individuals with diabetes [28]. Moreover, two recent meta-analyses also reported a link between poor vitamin D status and various lower extremity complications such as diabetic foot ulcers [29] and diabetic foot infection [30].

As the association between vitamin D and DSPN has mostly been investigated in cross-sectional studies and only in individuals with diabetes, it is yet unknown whether there is a temporal relationship between vitamin D and DSPN and whether such an association also exists in individuals with prediabetes. Also, most previous studies were small and did not adjust for the season of blood sampling, an important determinant of vitamin D levels. Moreover, despite the established causal relationship between obesity and low vitamin D status [31] and evidence indicating a positive association between obesity and DSPN [4], no previous study examined effect modification by obesity. Given these limitations, evidence from large-scale, well-designed prospective studies taking into account obesity status in older individuals with different glucose tolerance categories is required.

Therefore, this study aimed to examine cross-sectional and prospective associations between circulating levels of 25(OH)D and DSPN and explore a possible modifying effect of obesity and glucose tolerance categories on these associations in a large population-based cohort of individuals aged 62–81 years from Germany.

Material and methods

Study design and study population

This study was based on data from the Cooperative Health Research in the Region of Augsburg (KORA) F4 study (2006–2008) and KORA FF4 study (2013–2014), both follow-up examinations of the population-based KORA S4 study (1999–2001) conducted in Augsburg and two adjacent counties in Southern Germany. The design of the KORA study has been described elsewhere [32]. The study was carried out in accordance with the Declaration of Helsinki, including obtaining written informed consent from all participants and approved by the ethics board of the Bavarian Chamber of Physicians (Munich, Germany).

We included 1,161 participants of the KORA F4 study aged 62–81 years (i.e., the age group for whom DSPN was assessed). As described in Fig. 1, exclusions for various reasons resulted in a sample size of 1065 participants for the complete case analysis of prevalent DSPN. For the prospective analysis, we further excluded 443 participants who did not participate at FF4 for various reasons (i.e., illness, unavailability, non-interest or change of the geographic location), 19 participants with missing data for the Michigan Neuropathy Screening Instrument (MNSI) score at FF4, and 181 participants with prevalent DSPN at F4. These exclusions resulted in a sample size of 422 participants for the complete case analysis of incident DSPN.

Fig. 1: Flow chart describing the selection of the study population.
figure 1

25(OH)D, 25-hydroxyvitamin D; DSPN, distal sensorimotor polyneuropathy; KORA, Cooperative Health Research in the Augsburg Region; MNSI, Michigan Neuropathy Screening Instrument.

Assessment of glucose tolerance categories

After an overnight fast, all participants without known diabetes in KORA F4/FF4 received a standard 75-g oral glucose tolerance test (OGTT). Four glucose tolerance categories were defined using fasting and 2-h glucose levels according to the American Diabetes Association criteria: normal glucose tolerance (NGT), impaired fasting glucose (IFG), impaired glucose tolerance (IGT), and type 2 diabetes (T2D). IFG and/or IGT were considered as “prediabetes”. Known T2D in KORA F4 was defined as a self-report that was validated by the responsible physician or as the current self-reported use of glucose-lowering drugs.

Measurement of serum 25(OH)D

At F4, serum 25-hydroxyvitamin D (25(OH)D) concentrations were measured with the LIAISON® 25OH Vitamin D TOTAL Assay (DiaSorin Inc., Stillwater, MN, USA). The minimum detectable limit was 10 nmol/L, the interassay coefficient of variation (CV) was 8.7% and 9.1% for target values of 36.75 nmol/L and 118.75 nmol/L, respectively.

Assessment of DSPN

DSPN was assessed using the foot examination part of the MNSI, which includes four items: the appearance of feet (0 = normal, 1 = abnormal), foot ulcers (0 = normal, 1 = abnormal), ankle reflexes (0 = present, 0.5 = present with reinforcement, 1 = absent) and the vibration perception at the great toes (0 = present, 0.5 = reduced, 1 = absent) [33]. Age-dependent limits were used to classify the normal vibration perception threshold [34]. MNSI was extended by bilateral examination of sensory perception using a 10-g monofilament (Neuropen) (0 = normal, 1 = abnormal), resulting in a total MNSI score ranging from 0 (all aspects are normal) to a maximum of 10 points. An MNSI score >2 points was used to define prevalent and incident DSPN as previously suggested [32, 35]. This definition satisfies the diagnostic criteria for “possible DSPN” according to the Toronto Diabetic Neuropathy Expert Group [36].

Assessment of covariates

At the F4 visit, body weight and height were measured using standardized equipment and procedures. Body mass index (BMI) was calculated as weight (in kilograms) divided by height (in meters squared). Waist circumference (WC) was measured midway between the lower rib margin and the iliac crest. Obesity (BMI ≥ 30 kg/m2) and abdominal obesity (WC ≥ 102 cm for men and ≥88 cm for women) were defined following the NIH cutoff points [37].

For each participant, morning blood was collected and serum samples were aliquoted and stored until analysis. The date of blood draw was used to determine the season for 25(OH)D measurement: fall (September–November), winter (December–February), spring (March–May) and summer (June–August). HbA1c was measured using cation-exchange high-performance liquid chromatographic photometric assays on an Adams HA-8160 hemoglobin analysis system (Menarini Diagnostics, Florence, Italy) [38]. Total cholesterol and triacylglycerol were measured using a Dimension RxL instrument (Dade Behring) [38]. The estimated glomerular filtration rate (eGFR, mL/min/1.73m2) was calculated according to the Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) equation (2012) based on both serum creatinine and serum cystatin C [39]. Beta-nerve growth factor (beta-NGF) was measured with the inflammation panel from OLINK proteomics (Uppsala, Sweden) as described before [40].

Information on age, lifestyle and medical history was obtained during the F4 visit by trained interviewers using standardized questionnaires. Alcohol consumption was graded as none (0 g/day), moderate (<20 g/day for women, <40 g/day for men), high (20–<60 g/day for women, 40–<80 g/day for men) and heavy (≥60 g/day for women, ≥80 g/day for men). Smoking status was categorized into current smokers (including regular and irregular smokers), former smokers and never smokers. Physical activity was assessed according to duration of leisure time activities ([1] >2 h/week, [2] 1–2 h/week, [3] <1 h/week, [4] none) separately in winter and in summer. Participants who had a total sum score for physical activity in winter and in summer ≥5 were classifed as “physically inactive”, otherwise “physically active”. Hypertension was defined as systolic blood pressure ≥140 mmHg and/or diastolic blood pressure ≥90 mmHg and/or use of anti-hypertensive medication given that participants were aware of having hypertension. Self-reported history of stroke and myocardial infarction (yes/no) were also recorded.

Information on the regular use of lipid-lowering drugs and non-steroidal anti-inflammatory drugs (NSAIDs) during the last 7 days was collected through a database-supported computer software as described in detail elsewhere [32]. Briefly, participants were asked to bring all medications to the study center for verification of use. Intake of vitamin D supplements or cod liver oil during the last 6 weeks was assessed by a separate interview question.

Statistical analysis

Summary results are presented as medians (25th percentile, 75th percentile) or percentages. Differences in characteristics between participants with and without DSPN were calculated with the Kruskal–Wallis non-parametric rank test and the chi-square test. Correlations of serum 25(OH)D with continuous variables were estimated with Spearman’s correlation coefficients.

The associations of baseline serum 25(OH)D with prevalent and incident DSPN (as main dependent variables) were assessed with Poisson regression with robust error variance. Four models with increasing complexity were built to assess the impact of the covariates. Model 1 (minimally adjusted) included sex, age and season of blood draw as covariates. Model 2 was additionally adjusted for BMI. Model 3 was further adjusted for HbA1c, total cholesterol, log-transformed triacylglycerol, physical activity, alcohol consumption, smoking, eGFR, hypertension, stroke, myocardial infarction, use of lipid-lowering drugs and use of non-steroidal anti-inflammatory drugs. Lastly, we adjusted for vitamin D supplementation in model 4 (fully adjusted main model). All covariates used in the four regression models were measured at baseline (F4). Results are presented as risk ratios (RRs) with their 95% confidence intervals (CIs). Baseline serum 25(OH)D levels were modeled as a continuous variable (per 10 nmol/L decrease) and as a categorical variable categorized into tertiles. In the analysis using tertiles, tertile three was used as the reference category, and the linear trend test was performed by treating serum 25(OH)D as a continuous variable.

To explore the potential effect modification of obesity and glucose tolerance categories (T2D vs. prediabetes/NGT or NGT vs. prediabetes/T2D) on the associations between baseline serum 25(OH)D and prevalent and incident DSPN, we performed stratified analyses in the fully adjusted model (model 4). p-values for interaction were obtained by adding a cross-product term between 25(OH)D and the potential effect modifier.

We used restricted cubic splines to assess the shape of the association between 25(OH)D and the MNSI score (continuous outcome) [41]. Knots were placed at the 5th, 50th, and 95th percentiles of baseline serum 25(OH)D.

All statistical tests were 2-sided with an alpha level of 0.05, and all analyses were carried out using SAS statistical software (version 9.4; SAS Institute, Inc., Cary, NC, USA).

Results

Characteristics of the study population

Baseline characteristics of the 1065 participants included in the cross-sectional analysis are presented in Table 1. The median age was 70 years, and 51% were men. Obesity was present in 33% of the study sample, abdominal obesity in 60%, T2D in 22% and serum 25(OH)D levels <25 nmol/L in 17%. The median baseline serum 25(OH)D level was 44.3 nmol/L. DSPN was present in 35% of the study sample. Compared to individuals without prevalent DSPN, those with prevalent DSPN were older, taller and more likely to have obesity and abdominal obesity and higher HbA1c values. However, they had lower cholesterol levels and lower eGFR values. In addition, they reported higher alcohol consumption, less physical activity and had a higher prevalence of T2D and myocardial infarction.

Table 1 Baseline characteristics of the study population used in the analysis of prevalent DSPN.

Baseline characteristics of the 422 participants included in the prospective analysis are presented in Table 2. Participants were followed for a mean duration of 6.5 (standard deviation 0.5) years, during which 189 (45%) developed DSPN. Individuals with incident DSPN were older, taller and had higher BMI and waist circumference, more often hypertension but lower eGFR values. They were also more likely to be physically inactive.

Table 2 Baseline characteristics of the study population used in the analysis of incident DSPN.

Spearman correlation showed that older age, higher BMI, poor glycaemic control, higher cholesterol, and triacylglycerol levels were all negatively correlated with serum 25(OH)D (all p < 0.05). Higher eGFR values were positively correlated with serum 25(OH)D (p = 0.005). There was no correlation between 25(OH)D and beta-NGF (p = 0.96).

Association of serum 25(OH)D with prevalent DSPN

In the total sample, there was no significant association between serum 25(OH)D and prevalent DSPN in all models using 25(OH)D as a continuous variable or tertiles (Table 3). However, in the fully adjusted model, the analysis stratified by obesity (Table 4) identified a significant association between a lower serum 25(OH)D levels and prevalent DSPN in participants with obesity (RR (95% CI) per 10 nmol/L; 1.08 (1.01, 1.16)), whereas there was no evidence for an association in participants with normal weight (pinteraction = 0.002). When 25(OH)D was modeled as tertiles, the association of the lowest vs. the highest tertile was also stronger in participants with obesity (tertile1: RR (95% CI) 1.17 (0.98, 1.39); tertile 2: RR (95% CI) 1.38 (0.97, 1.95) than in normal-weight participants (pinteraction = 0.028).

Table 3 Associations of baseline serum 25(OH)D levels with prevalent and incident DSPN.
Table 4 Associations of baseline serum 25(OH)D levels with prevalent and incident DSPN according to obesity.

To assess the shape of the association between 25(OH)D and DSPN, we incorporated restricted cubic splines into a multivariable regression model with MNSI score at F4 as an outcome in individuals with obesity and normal-weight individuals separately (Fig. 2). In participants with obesity, the association between serum 25(OH)D and MNSI score was linear (p for nonlinearity = 0.790; p for overall association = 0.022), whereas in the normal-weight participants, the association was also linear (p for nonlinearity = 0.373) but not significant (p for overall association = 0.663).

Fig. 2: Dose–response relationship plots of the association between serum 25-hydroxyvitamin D levels and MNSI score at F4 in individuals with obesity (left panel) and individuals without obesity (right panel).
figure 2

Black lines, estimate; dashed lines, confidence limits; dots, knots set at the 5th, 50th, and 95th percentile of serum 25-hydroxyvitamin D. The reference was the 50th percentile. Associations were adjusted for sex, age (continuous), season of blood draw, HbA1c (continuous), total cholesterol (continuous), triacylglycerol (log-transformed, continuous), physical activity (active/inactive), alcohol consumption (none/moderate/high), smoking (current/former/never), eGFR (continuous), hypertension (yes/no), stroke (yes/no), myocardial infarction (yes/no), use of lipid-lowering drugs (yes/no), use of non-steroidal anti-inflammatory drugs (yes/no), and vitamin D supplementation (yes/no).

There was no interaction between serum 25(OH)D and glucose tolerance categories (pinteraction = 0.939 for NGT vs. prediabetes/T2D, pinteraction = 0.279 for NGT/prediabetes vs. T2D).

Association between 25(OH)D and incident DSPN

Table 3 shows the associations between serum 25(OH)D and incident DSPN in the total sample (n = 422). No significant association was found between 25(OH)D as a continuous variable and incident DSPN. However, in model 1 adjusted for age, sex and season, serum 25(OH)D levels in the first and second tertiles were associated with a higher risk of incident DSPN (RR (95% CI) 1.20 (1.03, 1.40) and RR (95% CI) 1.46 (1.08, 1.97), respectively) compared to tertile 3. These associations remained robust after additional adjustment for BMI but were attenuated after full adjustment for covariates.

Stratified analysis by obesity did not show any difference between participants with obesity and those with normal weight (pinteraction = 0.269 for 10 nmol/L and pinteraction = 0.209 for tertiles).

There was also no interaction between serum 25(OH)D and glucose tolerance categories (pinteraction = 0.094 for NGT vs. prediabetes/T2D, pinteraction = 0.722 for NGT/prediabetes vs. T2D).

Discussion

This study showed that lower serum 25(OH)D was associated with prevalent DSPN in individuals with obesity but not in normal-weight individuals in the older population. There were no differences in the associations according to glucose tolerance categories. Prospective results showed an association between lower serum 25(OH)D and higher risk of DSPN in the total sample, which was abolished after extensive adjustment for confounders.

Previous epidemiological studies reported an association between low serum 25(OH)D and DSPN [15,16,17,18,19,20,21,22,23,24,25,26]. However, all these studies were predominantly conducted in individuals with type 2 diabetes, and none examined the association between serum 25(OH)D and DSPN in individuals with and without obesity. Our study is the first to use a population-based design and to assess effect modification by obesity. Studies examining the association between 25(OH)D and DSPN only in individuals with diabetes might have limitations. First, they do not consider recent research indicating the increasing prevalence of DSPN in prediabetes [2, 42]. Second, observational studies [14], Mendelian randomization studies [43] and several RCTs [44] suggest a role of low 25(OH)D levels in increasing insulin resistance, decreasing β-cell function and promoting the development of T2D. Thus, even after adjusting for diabetes duration and glycaemic control, it is difficult to know whether the observed associations between low 25(OH)D levels and DSPN are genuine associations or rather an epiphenomenon as previously suggested [18].

Obesity has long been associated with low vitamin D status irrespective of age and ethnicity [45]. Using a genetic approach that limits confounding, a Mendelian randomization study including 21 cohorts and up to 42,000 participants confirmed a causal association between higher BMI and lower circulating 25(OH)D levels [31]. Several mechanisms were suggested to explain lower vitamin D in individuals with obesity. Formerly, vitamin D was thought to sequestrate in the excessive adipose tissue in individuals with obesity, decreasing its release in the circulation [45]. Recent studies suggested that a volumetric dilution rather than sequestration best explains the low vitamin D status in individuals with obesity [46] and that obesity decreases the activity of the hepatic 25-hydroxylase enzyme responsible for converting vitamin D into 25(OH)D [47].

Of note, obesity increases the risk of DSPN both in older individuals from the general population [4] and in those with type 2 diabetes [5]. Obesity is associated with chronic systemic inflammation where enlarged adipocytes release a variety of pro-inflammatory adipokines [48]. Using data from the KORA study and a mediation analysis approach, inflammation has been proposed as a biologically plausible link between obesity and DSPN [4]. In particular, three biomarkers of subclinical inflammation, including two chemokines (C-C motif chemokine ligand 7 [CCL7], C-X-C motif chemokine ligand 10 [CXCL10]) and the soluble form of the transmembrane receptor (Delta/Notch-like epidermal growth factor-related receptor [DNER]) partially mediate the association between obesity and DSPN [4].

Our data did not show a significant link between serum 25(OH)D and DSPN in the total sample in the cross-sectional and prospective analyses. Moreover, we believe that the magnitude of the associations in individuals with obesity is small. Nevertheless, given previous epidemiologic evidence linking vitamin D to DSPN, we cannot rule out a possible role of vitamin D in DSPN.

Several mechanisms may explain the role of vitamin D in the pathogenesis of DSPN. First, vitamin D has neuroprotective actions, which might be reduced in low 25(OH)D status. In vitro studies using murine [12] and rat cells [49] showed that treatment of Schwann cells with 1,25(OH)2D3, the active form of vitamin D, stimulates NGF and insulin-like growth factor (IGF-1) expression. NGF is a neurotrophin that plays an essential role in repairing nerve damage [50]. IGF-1 is a growth factor that stimulates neurogenesis and accelerates myelination and axonal regrowth [11]. We observed higher levels of beta-NGF in individuals with prevalent and incident DSPN, which might be explained by an upregulation to promote nerve regeneration and repair. However, no correlation between 25(OH)D and beta-NGF was found in our study, in contrast to a study in individuals with type 1 diabetes, which reported a positive correlation [17]. A study in rodent sciatic nerves showed that treatment with vitamin D protects Schwann cells against apoptosis induced by advanced glycation end-products [51]. In vitro studies using human cells demonstrated that treating neuroblastoma cells with vitamin D and other growth factors induced neuronal cell differentiation and generated sustainable and well-differentiated neurons [10]. Furthermore, studies in vitamin D receptor-deficient mice underlined the role of vitamin D in axonal homogeneity [49].

Second, vitamin D has anti-inflammatory and immunomodulatory actions, which might be reduced in low 25(OH)D status contributing to the pathogenesis of DSPN. Notably, a severe vitamin D deficiency has been found associated with increased expression of inflammatory cytokines [52], particularly interleukin-6 and tumor necrosis factor-α, both previously reported to increase the risk of DSPN [53].

In individuals with obesity and low vitamin D, higher circulating levels of biomarkers of inflammation due to obesity along with reduced anti-inflammatory and neuroprotective actions of vitamin D could represent mechanisms that might explain the association between low 25(OH)D levels and prevalent DSPN in our study.

Studies comparing serum 25(OH)D in individuals with diabetes showed lower serum or plasma 25(OH)D levels in individuals with DSPN than those without DSPN [15, 16, 21, 22, 24, 25] in contrast to our results. However, we still did find an association between low 25(OH)D levels and incident DSPN in the initial models. Differences in findings might be explained by the study population (individuals with diabetes vs. individuals from the general population), by between-study heterogeneity in the assessment of DSPN (i.e., self-report, neuropathy symptom score, neuropathy disability score, vibration perception threshold, Semmes-Weinstein monofilament testing and nerve conduction studies) and by the more extensive adjustment for potential confounders in our study. In particular, it is important to note that our study is the only study that took into account the confounding effect of the season on the associations between 25(OH)D and DSPN. Despite the known seasonal fluctuation of vitamin D levels, none of the previous studies had adjusted for the season of blood draw for 25(OH)D measurements except one study, which adjusted for sunlight exposure and daily activity. However, it included only individuals for whom blood was collected during the summer [26].

Our study shows that only 6% of the participants took vitamin D supplements or cod liver oil. This proportion was below our expectations given the widespread use of vitamin D supplements among older individuals. An explanation could be the neglect of vitamin D supplementation in older individuals, often having a high proportion of polypharmacy to manage ageing-related diseases. Intervention studies testing the effect of vitamin D supplementation on neuropathic pain indicated that supplementation with high, single or weekly dosages of vitamin D was associated with a significant decrease in the signs and symptoms of DSPN [54, 55]. However, these trials were small-scale [56], not placebo-controlled, restricted to individuals with T2D and painful diabetic neuropathy, and supplemented with vitamin D doses based on a “one-size-fits-all” approach without considering individual differences. However, it remains unknown whether vitamin D supplementation could be an effective therapeutic regimen that can change the pathological features of DSPN and its progression. It is worth noting that vitamin D supplementation in other neurological diseases did not seem effective in modifying the disease phenotype for reasons probably related to patient selection. Therefore, future trials assessing the efficacy of vitamin D supplementation as a disease-modifying treatment for DSPN should carefully consider individual differences, particularly differences in obesity status, to increase their ability to detect the supplementation’s beneficial effect.

Strengths and limitations

This study includes a large sample size, and its design allows for both cross-sectional and prospective analyses with extended follow-up and comprehensive adjustment for confounders. Also, KORA is a population-based study, which is optimal to test the hypothesis of effect modifications by obesity or glucose tolerance categories on the associations between 25(OH)D and DSPN. Some limitations of our study should be noted. First, the inclusion of white and older individuals from one country limits the generalizability of our findings. As vitamin D is mainly synthesized in the skin and strongly depends on skin color and age, future studies should aim to replicate our findings in populations of other ethnicities and among young and middle-aged individuals. Another limitation of this study is that 25(OH)D was not measured with gold standard methods, i.e., high-performance liquid chromatography–electrospray ionization–mass spectroscopy. DSPN was not assessed with nerve conduction studies, thus does not satisfy the definition of “confirmed DSPN”. A further limitation was the smaller sample size of the prospective compared to the cross-sectional analysis. This limitation indicates that our prospective analysis may not have been adequately powered, which may explain the lack of association between serum 25(OH)D levels and incident DSPN in the fully adjusted models. Finally, the observational study design cannot rule out potential residual confounding. Future studies based on Mendelian Randomization (MR) analyses that use instrumental variables (IVs) are needed to strengthen evidence of a causal association between 25(OH)D and DSPN because they are not subject to biases with confounding and reverse causation.

Conclusion

Our findings from a large, well-designed observational cohort found no consistent associations between serum 25(OH)D and DSPN in cross-sectional and prospective analyses. However, our data suggests an association between lower serum 25(OH)D levels and higher DSPN prevalence in individuals with obesity. Further large-scale prospective population-based studies are required to confirm and extend our findings. In a clinical setting, intervention studies are needed to clarify whether body weight loss (i.e., by bariatric surgery or lifestyle interventions) in individuals with obesity can effectively increase vitamin D levels and subsequently prevent DSPN. Alternatively, well-designed randomized clinical trials are needed to determine the benefits of vitamin D supplementation, taking into account body weight, in preventing the development of DSPN.