Skip to main content

The relationship between vitamin D status and islet function in patients with type 2 diabetes mellitus



The aim of the study was to explore the relationship between vitamin D status and islet function in patients with type 2 diabetes mellitus.


The participants were recruited from Hebei General Hospital. Basic characteristics and blood indicators were collected after fasting overnight. The data were analyzed statistically using SPSS 22.0. Analysis of variance, a nonparametric test, or a trend Chi-square test was used for the comparisons. The association between 25-hydroxy vitamin D and modified homeostasis model assessment-β was assessed using multivariate ordinal logistic regression.


One hundred seventy-four patients aged 26 to 79 years with type 2 diabetes mellitus were included in this study. Patients with vitamin D deficiency had a lower modified homeostasis model assessment-β level compared with those without vitamin D deficiency. There were differences in body mass index, diabetes course, glycosylated hemoglobin, fasting blood glucose, fasting blood C-peptide, triglyceride, and 25-hydroxy vitamin D among different modified homeostasis model assessment-β groups based upon the tertiles. 25-hydroxy vitamin D, as continuous or categorical variables, was positively related to modified homeostasis model assessment-β whether or not cofounding factors were adjusted.


There is an association between increased 25-hydroxy vitamin D levels and improvement in modified homeostasis model assessment-β function in patients with type 2 diabetes mellitus.

Trial registration

Cross-sectional trails ChiCTR2000029391, Registration Date: 29/01/2020.

Peer Review reports


Widespread type 2 diabetes mellitus (T2DM) has caused a burden worldwide. Evaluation of islet function is important for patients with T2DM because it is usually used to assess the progress of T2DM and to guide drug usage. Traditionally, homeostasis model assessment (HOMA)-β has been used to evaluate islet function using fasting blood glucose (FBG) and insulin levels [1], but this is only used in patients who have not received insulin therapy; therefore, its application is limited for patients who have received injected insulin. Li [2] found that a modified HOMA-β that replaces insulin with fasting C-peptide (FCP) could be used in all patients whether or not they have taken insulin.

The relationship between vitamin D and numerous diseases including metabolic abnormalities such as insulin resistance (IR) and T2DM has been identified over time [3]. It is generally recognized that 25-hydroxy vitamin D (25OHD) represents the patient’s vitamin D status. Vitamin D deficiency has been defined as 25OHD < 20 ng/mL and vitamin D insufficiency has been defined as 25OHD < 30 ng/mL [4]. Vitamin D insufficiency or deficiency was strongly related with accelerated development of IR [5]. A low 25OHD level was related to an increased risk of T2DM and poor glucose control. A meta-analysis of randomized controlled trials (RCTs) showed that vitamin D supplementation decreased glycosylated hemoglobin (HbA1c) and FBG levels in T2DM patients with vitamin D deficiency or without obesity However, for patients with vitamin D insufficiency or sufficiency, vitamin D supplementation did not change the status of glycemic control [6]. Furthermore, low 25OHD level increased the risk or the severity of complications, such as diabetic peripheral neuropathy [7] and diabetic retinopathy [8].

Previous studies explored the relationship between vitamin D status and islet function using traditional HOMA-β. The aim of the present study was to explore the association between 25OHD and the modified HOMA-β in patients with T2DM.

Materials and methods

Study design

This was a cross-sectional study. This study was registered at Clinical Trial (Registration Number: ChiCTR2000029391, Registration Date: 29/01/2020) and approved by the Hebei General Hospital Ethics Committee. It was performed in accordance with the Declaration of Helsinki.

Study population

The population was from the Department of Endocrinology of Hebei General Hospital, China. All participants were hospitalized between June 2018 and December 2019. We included adult patients (≥18 years) with T2DM. A diagnosis of diabetes mellitus (DM) was made on the basis of the World Health Organization definition, as follows: FPG level ≥ 7.0 mmol/L or 2-h post-load value ≥11.1 mmol/L. Patients who had incomplete data or were diagnosed with severe diseases (e.g., myocardial infarction, cerebral hemorrhage, or malignancy) or infectious diseases (e.g., pneumonia or urinary infection) within 3 months were excluded. Two hundred and forty-eight patients were enrolled and 174 patients were included in this cross-sectional study. On the basis of vitamin D deficiency, the participants were grouped into control group (n = 69) and vitamin D deficiency group (n = 105). According to tertiles, modified HOMA-β was grouped three groups and 58 participants were included in each group.


Basic data were collected by questionnaire, which included sex, age, DM course, diabetes history, smoking, alcohol drinking, hypertension, and insulin use. Height in cm and weight in kg were measured without shoes and with light clothes, respectively, and the body mass index (BMI) was calculated using the body mass index (BMI; kg/m2), as follows: weight (kg) / height2 (m2). Height and weight were rounded up to the nearest 0.1 and the mean value was obtained by measuring twice and taking the average of the two measurements.

Blood samples from all participants were collected from the antecubital vein after 8 h of fasting. The same indicators for all samples were tested using one machine in one laboratory. These indicators included FBG, FCP, HbA1c, total protein, albumin, the lipid profile (total cholesterol [TC], triglyceride [TG], high-density lipoprotein cholesterol [HDL-C], low-density lipoprotein cholesterol [LDL-C]), and bone turnover biomarkers (25OHD, β-C-terminal cross-linked telopeptide of type I collagen [β-CTX], procollagen type 1 N-terminal propeptide [P1NP], and parathyroid hormone [PTH]). The modified HOMA-β was calculated as follows: 270 × FCP / (FBG − 3.5) × 0.333 (FCP in ng/mL and FBG in mmol/L).

All data were inputted by one researcher and checked by a second researcher. If errors were observed, data were checked and revised under two researchers’ supervision.

Statistical analysis

The data were analyzed statistically using SPSS 22.0 (IBM, Inc., Armonk, NY, USA). The modified HOMA-β results were grouped into tertiles. Differences in the measurement data were compared using an analysis of variance (ANOVA) if the data had a normal distribution, and these data were expressed as the mean ± standard deviation (SD). A nonparametric test was used if the data had a skewed distribution, and these data were expressed as the median (interquartile range). A trend Chi-square test was used for counting data, and these data were expressed as the number (percentage). The association between 25OHD and modified HOMA-β was assessed using linear correlation analysis. After comparisons among different groups, cofounding factors were screened based on a P-value < 0.1. We used multivariate ordinal logistic regression analysis to explore the relationship of vitamin D and modified HOMA-β. Results were shown as the odds ratio (OR) and 95% confidence interval (95%CI). In Model 1, no cofounding factors were adjusted. In Model 2, the DM course, age, and BMI were adjusted. In Model 3, the DM course, age, BMI, HbA1C, TG, β-CTX, and PINP were adjusted. A P-value < 0.05 was considered to be statistically significant.


Characteristics of the included population

A flow diagram for patient recruitment is shown in Fig. 1. The study population consisted 174 patients (122 men and 52 women) aged 26 to 79 with T2DM. The average HbA1c was 9.1% (7.6, 10.9). 25OHD levels ranged from 4.37 to 38.6 ng/mL. The participants’ baseline characteristics and blood indicators are shown in Table 1.

Fig. 1

Flow diagram of participant recruitment into the study

Table 1 Participants’ basic characteristics

Modified HOMA-β levels based on vitamin D deficiency

Vitamin D deficiency was defined as 25OHD < 20 ng/mL. On the basis of the presence or absence of vitamin D deficiency, the modified HOMA-β level decreased significantly in the vitamin D deficiency group compared with the control group (no vitamin D) (394.7 [259.42, 564.76] vs. 536.44 [318.65, 762.12]; p = 0.027; Fig. 2).

Fig. 2

Comparison of modified HOMA-β between the control group and the vitamin D deficiency group. Note: HOMA, homeostasis model assessment

Characteristics of the participants based on modified HOMA-β tertiles

Table 2 shows the comparisons among different modified HOMA-β groups on the basis of the tertiles. BMI was higher in T2 and T3 groups than that in T1 group (P-value = 0.041). There were differences in the DM course among different modified HOMA-β groups (P-value = 0.007). Compared with T1 group, the HbA1c level was decreased in the other two groups (P-value < 0.001). The FBG level decreased from the T1 group to the T3 group (P-value < 0.001), and FCP showed the opposite trend (P-value < 0.001). TG and 25OHD levels were higher in the T3 group than in the T1 and T2 groups (P-value = 0.003 and 0.005, respectively). We screened BMI, DM course, age, HbA1c, TG, 25OHD, β-CTX, and P1NP as confounding factors on the basis of P < 0.1. Because FBG and FCP were needed to calculate the modified HOMA-β, they were not seen as confounding factors.

Table 2 Comparisons of clinical characteristics in type 2 diabetic patients on the basis of the modified HOMA-β grouping

Linear correlation analysis of 25OHD and modified HOMA-β

Linear correlation analysis showed that as the 25OHD level increased, the modified HOMA-β level also increased (r = 0.199, P = 0.009; Fig. 3).

Fig. 3

Linear correlation between 25OHD and modified HOMA-β. Note: 25OHD, 25-hydroxy vitamin D; HOMA, homeostasis model assessment.

Correlation of 25OHD and modified HOMA-β

The first ordinal logistic regression analysis was conducted where 25OHD was seen as a continuous variable (shown in Table 3). In Model 1, no confounding factors were adjusted, and this analysis showed that 25OHD increased by 1 ng/mL, and the modified HOMA-β increased by 5.7%. In Model 2, the DM course, age, and BMI were adjusted, and the modified HOMA-β increased by 6.9%. In Model 3, HbA1C, TG, β-CTX, and PINP were adjusted based on Model 2, and the age increased by 1 ng/mL and modified HOMA-β improved by 5.9%.

Table 3 Multivariate ordinal logistic regression analysis of relationship between 25OHD and modified HOMA-β (25OHD as continuous variable)

A second ordinal logistic regression analysis was conducted where 25OHD was considered to be a categorical variable (no vitamin D deficiency and vitamin D deficiency) (shown in Table 4). In Model 1, no confounding factors were adjusted, and the modified HOMA-β in patients with no vitamin D deficiency increased by 2.184-fold compared with patients with vitamin D deficiency. In Model 2, when DM course, age, and BMI were adjusted, the OR was 2.520. In Model 3, HbA1c, TG, CTX, and PINP were adjusted based on Model 2, and the modified HOMA-β in patients with no vitamin D deficiency increased 2.234-fold compared with compared with those with vitamin D deficiency.

Table 4 Multivariate ordinal logistic regression analysis of relationship between 25OHD and modified HOMA-β (25OHD as categorical variable)


In this study, we investigated the relationship between vitamin D status and modified HOMA-β. To eliminate potential confounding factors, the P-value was limited to 0.1. Although there was a weak relationship between 25OHD levels and modified HOMA-β, the relationship may also be affected by other factors and sample size. Previous studies have shown a correlation between these two variables [2]. We further explored the relationship between 25OHD levels and modified HOMA-β through logistic regression analysis, while considering potential confounding factors. In the analysis process, we considered 25OHD to be a continuous and categorical variable, and we conducted a multivariate ordinal logistic regression analysis twice to obtain better and more plausible results. Whether or not confounding factors were adjusted, the results showed that 25OHD was positively related to the modified HOMA-β, meaning that, with an increase in 25OHD, the modified HOMA-β function improved. In addition, the relationship between 25OHD and fat mass is controversial. Some studies have shown that 25OHD levels and fat mass are negatively related [9], whereas other studies have not shown this relationship. Therefore, fat mass may be a confounder for the final results [10]. However, we did not measure the patients’ fat mass in our study.

Beyond the effect of vitamin D on calcium and bone metabolism, the role of vitamin D on IR, glucose metabolism, and T2DM has received widespread attention. Vitamin D can regulate insulin secretion and mediate pancreatic β cell survival, which was shown in previous studies [11, 12].

The mechanisms of vitamin D that regulate islet function may be involved in several aspects. First, the vitamin D receptor (VDR) was expressed in pancreatic β cells, and vitamin D can exert its effect by directly binding to the VDR [13]. Mice lacking VDR showed impairment in insulin secretion after glucose loading [14]. Vitamin D can also promote insulin secretion because it can identify the insulin gene promotor’s vitamin D-responsive elements (VDRE). The VDRE regulates insulin gene transcription and is involved in intracellular junctions and cellular growth in pancreatic β cells [15]. Second, 25(OH)D-1α-hydroxylase (CYP27B1) is also expressed in pancreatic β cells, which is an enzyme that mediates 25OHD to transform it into the active 1,25(OH)2D3. Third, intracellular calcium (Ca2+) levels are essential for insulin action. Vitamin D is involved in regulating Ca2+ flux in pancreatic cells [16]. Fourth, vitamin D can act on pancreatic β cells and on insulin-sensitive tissues, such as white adipose tissue, muscle, and liver to regulate insulin sensitivity [11]. For example, 1,25(OH)2D3 improves glucose metabolism via the SIRT1/IRS1/GLUT4 signaling cascade in C2C12 myotubes [17]. In addition, vitamin D is indirectly related to oxidative stress [18], sub-inflammation [19], and epigenetics [20].

The hyperinsulinemic–euglycemic clamp is the gold standard that is currently used to quantitatively evaluate insulin secretion and IR. Its accuracy and stability are widely recognized. However, because of the invasiveness of surgery and technical difficulties, its use has been limited. In clinical studies, HOMA-IR and HOMA-β are an evaluation method that uses fasting glucose and insulin levels. Patients with T2DM, especially patients with poor islet function, usually inject insulin subcutaneously to have good glycemic control. Therefore, the modified HOMA-β is need for these patients.

The current study has some limitations. First, for cross-sectional studies, our sample size was small. The results are more convincing when there is a large sample size in a cross-sectional observational study. Second, exposure to sunshine, physical exercise, vitamin D, calcium intake, or drugs that treat osteoporosis may be confounding factors, but we did not take these factors into account. Third, because of this study type, we could not obtain a causal relationship between vitamin D and islet function. Fourth, patients with or without insulin therapy were all included in this study. It is more appropriate to recruit patients without insulin therapy to compare the traditional HOMA-β and modified HOMA-β.


In conclusion, taking all the potential confounding factors into account, the vitamin D status and modified HOMA-β showed a positive relationship. Thus, an increase in or supplementation of vitamin D may be related to improved modified HOMA-β function and delayed pancreatic β cell functional decline.

Availability of data and materials

The original data can be obtained by email request (



Analysis of variance


Apolipoprotein A1


Apolipoprotein B


Body mass index


Confidence interval




Diabetes mellitus


Fasting blood glucose


Fasting C-peptide


Glycated hemoglobin


High-density lipoprotein cholesterol


Homeostasis model assessment


Low-density lipoprotein cholesterol


25-hydroxyvitamin D


Standard deviation


Total cholesterol






Odds ratio


β-C-terminal cross-linked telopeptide of type I collagen


Procollagen type 1 N-terminal propeptide


Parathyroid hormone


Vitamin D receptor


Vitamin D-responsive elements


  1. 1.

    Fasipe OJ, Ayoade OG, Enikuomehin AC. Severity grade assessment classifications for both insulin resistance syndrome and status of pancreatic beta cell function in clinical practice using homeostasis model assessment method indices. Can J Diabetes. 2020;44(7):663–9.

    Article  PubMed  Google Scholar 

  2. 2.

    Li X, Zhou ZG, Qi HY, Chen XY, Huang G. Replacement of insulin by fasting C-peptide in modified homeostasis model assessment to evaluate insulin resistance and islet β cell function. Zhong Nan Da Xue Xue Bao Yi Xue Ban. 2004;29(4):419–23.

    CAS  PubMed  Google Scholar 

  3. 3.

    Hanmin W, Weiwen C, Dongqing L, Xiaoe Y, Xiaode Z, Nancy O, et al. Vitamin D and chronic diseases. Aging Dis. 2017;8(3):346–53.

    Article  Google Scholar 

  4. 4.

    Holick MF, Binkley NC, Bischoff-Ferrari HA, Gordon CM, Hanley DA, Heaney RP, et al. Evaluation, treatment, and prevention of vitamin D deficiency: an endocrine society clinical practice guideline. J Clin Endocrinol Metab. 2011;96(7):1911–30.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Szymczak-Pajor I, Śliwińska A. Analysis of association between vitamin D deficiency and insulin resistance. Nutrients. 2019;11(4).

  6. 6.

    Wu C, Qiu S, Zhu X, Li L. Vitamin D supplementation and glycemic control in type 2 diabetes patients: a systematic review and meta-analysis. Metabolism. 2017;73:67–76.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    He R, Hui YH, Zhao ZJ, Zhao J. Vitamin D deficiency increases the risk of peripheral neuropathy in Chinese patients with type 2 diabetes. Diabetes Metab Res Rev. 2017;33.

  8. 8.

    Ashinne B, Rajalakshmi R, Anjana RM, Narayan VKM, Jayashri R, Mohan V, et al. Association of serum vitamin D levels and diabetic retinopathy in Asian Indians with type 2 diabetes. Diabetes Res Clin Pract. 2018;139:308–13.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Oliai Araghi S, van Dijk SC, Ham AC, Brouwer-Brolsma EM, Enneman AW, Sohl E, et al. BMI and body fat mass is inversely associated with vitamin D levels in older individuals. J Nutr Health Aging. 2015;19(10):980–5.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    Forney LA, Earnest CP, Henagan TM, Johnson LE, Castleberry TJ, Stewart LK. Vitamin D status, body composition, and fitness measures in college-aged students. J Strength Cond Res. 2014;28(3):814–24.

    Article  PubMed  Google Scholar 

  11. 11.

    Pramono A, Jocken JWE, Blaak EE. Vitamin D deficiency in the aetiology of obesity-related insulin resistance. Diabetes Metab Res Rev. 2019;35(5):e3146.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Jayanarayanan S, Anju TR, Smijin S, Cheramadathikudiyil S, Paulose. Vitamin D3 supplementation increases insulin level by regulating altered IP3 and AMPA receptor expression in the pancreatic islets of streptozotocin-induced diabetic rat. J Nutr Biochem. 2015;26:1041–9.

    CAS  Article  Google Scholar 

  13. 13.

    Morró M, Vilà L, Franckhauser S, Mallol C, Casellas A. Vitamin D receptor overexpression in β-cells ameliorates diabetes in mice. Diabetes. 2020;69(5):927–39.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  14. 14.

    Zeitz U, Weber K, Soegiarto DW, Wolf E, Erben RG. Impaired insulin secretory capacity in mice lacking a functional vitamin D receptor. FASEB J. 2003;17(3):509–11.

    CAS  Article  PubMed  Google Scholar 

  15. 15.

    Wolden-Kirk H, Overbergh L, Gysemans C, Brusgaard K, Naamane N, Van Lommel L, et al. Unraveling the effects of 1,25OH2D3 on global gene expression in pancreatic islets. J Steroid Biochem Mol Biol. 2013;136:68–79.

    CAS  Article  PubMed  Google Scholar 

  16. 16.

    Rey D, Miranda Sulis P, Alves Fernandes T, Goncalves R, Silva Frederico MJ, Costa GM, et al. Astragalin augments basal calcium influx and insulin secretion in rat pancreatic islets. Cell Calcium. 2019;80:56–62.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Manna P, Achari AE, Jain SK. 1,25(OH)2-vitamin D3 upregulates glucose uptake mediated by SIRT1/IRS1/GLUT4 signaling cascade in C2C12 myotubes. Mol Cell Biochem. 2018;444(1-2):102–8.

    CAS  Article  Google Scholar 

  18. 18.

    Lee WC, Mokhtar SS, et al. Vitamin D status and oxidative stress in diabetes mellitus. Cell Mol Biol. 2018;64:60–9.

    Article  Google Scholar 

  19. 19.

    Filgueiras MS, Rocha NP, Novaes JF, Bressan J. Vitamin D status, oxidative stress, and inflammation in children and adolescents: a systematic review. Crit Rev Food Sci Nutr. 2018;60(4):660–9.

    CAS  Article  PubMed  Google Scholar 

  20. 20.

    Bahrami A, Sadeghnia HR, Tabatabaeizadeh SA, Bahrami-Taghanaki H, Behboodi N, Esmaeili H, et al. Genetic and epigenetic factors influencing vitamin D status. J Cell Physiol. 2017;233(5):4033–43.

    CAS  Article  PubMed  Google Scholar 

Download references


We thank Jodi Smith, PhD, from Liwen Bianji, Edanz Editing China (, for editing the English text of a draft of this manuscript.


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

Author information




Conceptualization: HZ and SCC. Data collection, Methodology and Analysis: HZ, CZ and MMZ. Writing and preparation of the original draft: HZ. Supervision: SCC. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Shuchun Chen.

Ethics declarations

Ethics approval and consent to participate

This study was approved by Hebei General Hospital Ethics Committee and participants signed consent.

Consent for publication

Not Applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.


Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Zhao, H., Zheng, C., Zhang, M. et al. The relationship between vitamin D status and islet function in patients with type 2 diabetes mellitus. BMC Endocr Disord 21, 203 (2021).

Download citation


  • Vitamin D
  • Islet function
  • Homeostasis model assessment-β
  • Type 2 diabetes