Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Open Access

Peer-reviewed

Research Article

Vitamin D status is inversely associated with markers of risk for type 2 diabetes: A population based study in Victoria, Australia

  • Poonam K. Pannu,

    Affiliation Food, Nutrition & Health, School of Public Health, Curtin Health Innovation Research Institute, Curtin University, Perth, Western Australia, Australia

    ⨯
  • Leonard S. Piers,

    Affiliation Health Intelligence Unit, System Intelligence and Analytics Branch, Health Strategy, Productivity and Analytics Division, Department of Health, Melbourne, Victoria, Australia

    ⨯
  • Mario J. Soares ,

    m.soares@curtin.edu.au

    Affiliation Food, Nutrition & Health, School of Public Health, Curtin Health Innovation Research Institute, Curtin University, Perth, Western Australia, Australia

    ⨯
  • Yun Zhao,

    Affiliation Occupation and the Environment, School of Public Health, Faculty of Health Sciences, Curtin University, Perth, Western Australia, Australia

    ⨯
  • Zahid Ansari

    Affiliation Health Intelligence Unit, System Intelligence and Analytics Branch, Health Strategy, Productivity and Analytics Division, Department of Health, Melbourne, Victoria, Australia

    ⨯

Abstract

A growing body of evidence suggests a protective role of vitamin D on the risk of type 2 diabetes mellitus (T2DM). We investigated this relationship in a population sample from one Australian state. The data of 3,393 Australian adults aged 18–75 years who participated in the 2009–2010 Victorian Health Monitor survey was analyzed. Socio-demographic information, biomedical variables, and dietary intakes were collected and fasting blood samples were analyzed for 25, hydroxycholecalciferol (25OHD), HbA1c, fasting plasma glucose (FPG), and lipid profiles. Logistic regression analyses were used to evaluate the association between tertiles of serum 25OHD and categories of FPG (<5.6 mmol/L vs. 5.6–6.9 mmol/L), and HbA1c (<5.7% vs. 5.7–6.4%). After adjusting for social, dietary, biomedical and metabolic syndrome (MetS) components (waist circumference, HDL cholesterol, triglycerides, and blood pressure), every 10 nmol/L increment in serum 25OHD significantly reduced the adjusted odds ratio (AOR) of a higher FPG [AOR 0.91, (0.86, 0.97); p = 0.002] and a higher HbA1c [AOR 0.94, (0.90, 0.98); p = 0.009]. Analysis by tertiles of 25OHD indicated that after adjustment for socio-demographic and dietary variables, those with high 25OHD (65–204 nmol/L) had reduced odds of a higher FPG [AOR 0.60, (0.43, 0.83); p = 0.008] as well as higher HbA1c [AOR 0.67, (0.53, 0.85); p = 0.005] compared to the lowest 25OHD (10–44 nmol/L) tertile. On final adjustment for other components of MetS, those in the highest tertile of 25OHD had significantly reduced odds of higher FPG [AOR 0.61, (0.44, 0.84); p = 0.011] and of higher HbA1c [AOR 0.74, (0.58, 0.93); p = 0.041] vs. low 25OHD tertile. Overall, the data support a direct, protective effect of higher 25OHD on FPG and HbA1c; two criteria for assessment of risk of T2DM.

Citation: Pannu PK, Piers LS, Soares MJ, Zhao Y, Ansari Z (2017) Vitamin D status is inversely associated with markers of risk for type 2 diabetes: A population based study in Victoria, Australia. PLoS ONE 12(6): e0178825. https://doi.org/10.1371/journal.pone.0178825

Editor: Leng Huat Foo, University Sains Malaysia, MALAYSIA

Received: September 26, 2016; Accepted: May 21, 2017; Published: June 2, 2017

Copyright: © 2017 Pannu et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: The data belong to the Victorian Government, Australia and we have signed a legal contract to analyze the de-identified data. Any interested researcher can access the data by contacting the Department of Health and Human Services at health.intelligence@dhhs.vic.gov.au or by phoning: The Manager, Health Intelligence Unit, Department of Health & Human Services, Victoria, Australia, on +61 3 9096 5242. Additional information may be found at https://www2.health.vic.gov.au/public-health/population-health-systems/health-status-of-victorians/survey-data-and-reports.

Funding: PKP is an Australian postgraduate research scholar.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Vitamin D status, as judged from circulating concentrations of 25, hydroxycholecalciferol (25OHD), is a worldwide concern. In many countries across all continents, approximately 50% of those populations have an inadequate 25OHD status (<50 nmol/L) [1]. Countries like India and China have some of the highest rates of vitamin D deficiency (<25 nmol/L) [1]. The 25OHD status of Australians is also surprisingly low for a country blessed with abundant sunshine. Current estimates indicate that ~31% of the population have inadequate 25OHD levels [2], with a higher prevalence in older Australians [3]. The prevalence of type 2 diabetes mellitus (T2DM) has also risen tremendously in the last 10 years, with projections that countries like India and China will have the highest numbers by 2030 at 79.4 and 42.3 million respectively [4].

There is an ongoing interest in the extra-skeletal effects of vitamin D including its potential to blunt the risk of developing T2DM. Positive outcomes would present a tangible public health solution, if causality is accepted. Such a relationship was first suggested in 1967 by Milner and Hales, who found that insulin secretion in rabbits was dependent on calcium and magnesium, which are tightly regulated by vitamin D [5]. Accumulating evidence has indicated that higher 25OHD status may have several anti-diabetic effects, including improvement in insulin sensitivity, stabilizing HbA1c levels [6], and improving beta cell function [7], whereas low 25OHD status may increase risk of T2DM [8]. Thus, in the current environment of increasing rates of T2DM [9], their close parallelism with insufficient levels of 25OHD deserves investigation in population based studies. There are several lifestyle factors that modulate the risk of T2DM, including dietary components and patterns [10, 11], physical activity, and smoking [12]. The risk of developing T2DM over 20 years appears to increase with the accumulation of metabolic syndrome (MetS) components. The risk of T2DM increased by: 11.9% in those with zero Mets components, 31.2% in those with three MetS components and 40.8% in those with four or five MetS components [13]. Though the presence of MetS components increases the risk of T2DM, glucose is the most strongly correlated factor in predicting the development of diabetes in the future [14]. In a study of more than 58,000 adults, as the number of components increased, so did the incidence of diabetes [14]. However, some gaps may exist with one study in Hispanic Americans finding that impaired glucose tolerance had a greater predictive power than the individual MetS components [15]. Thus, the presence of MetS is another major risk predictor of increased T2DM [16].

The aim of this study was to investigate the association of 25OHD, and the risk of T2DM. We [17] as well as others [18–21], have shown that higher 25OHD status significantly reduced the risk of MetS and its components. Hence it was essential to adjust for several lifestyle factors and components of MetS, other than glucose, in order to correctly identify any independent association between 25OHD and risk of T2DM. One other study [8] has also investigated the association between 25OHD and T2DM and adjusted for three out of the four MetS components. Thus our study appears to be one of the first to adjust for all MetS components. Impaired fasting plasma glucose (FPG) and HbA1c levels are now recommended as key determinants of early risk of T2DM [22]. While high FPG is an immediate indicator of poor glucose homeostasis, HbA1c is a better indicator of longer term control of blood glucose, and recommended cut-offs for both these biomarkers are used to diagnose T2DM [22]. The underlying hypothesis of the present investigation was that increases in 25OHD would reduce the odds ratio of a high FPG and a high HbA1c after adjustment for socio-demographic, dietary and biomedical confounders.

Materials and methods

Sample

The Victorian Health Monitor (VHM) survey was a state-wide cross-sectional population based study [23] conducted in Victoria, Australia. Victoria lies in the south-east of Australia, and has a latitude of 37°47'S and longitude of 144°58'E. Data was collected between May 2009 and April 2010 including: physical information, dietary behaviour information and biomedical information. The physical and biomedical information of participants were collected by trained staff at four training sites. The VHM employed a stratified cluster sample selection method of Census Collection Districts within eight Department of Health regions in Victoria. Data were collected on 3,653 adults aged 18–75 years. From this sample, we excluded participants with: missing HbA1c and FPG data (n = 31), those with HbA1c ≥6.5% (n = 39) and FPG >7 mmol/L (n = 16) as they were classified as having T2DM as per the American Diabetes Association (ADA) cut-offs [22], those with T2DM (n = 140), those with type 1 diabetes (n = 9), and those on diabetic medications (n = 25). A total of 3,393 subjects were included in this analysis. Further details on physical, dietary, and biomedical data collection and analysis have been previously described [17, 23, 24].

Biomedical measurements

Participants attended a testing site after an overnight fast of at least 10 hours. Blood samples were collected by venepuncture, and were subsequently transported to a central laboratory in Melbourne, Australia. Bloods were analysed for: FPG, HbA1c, 25OHD, high density lipoprotein cholesterol (HDL-C), and triglycerides (TG). The components in the blood were measured as follows: FPG using the hexokinase method, HbA1c using immunosassay (Roche Integra chemistry analyser), 25OHD concentration were measured based on the DiaSorin Corporation Liaison® 25OHD total assay HDL-C using elimination/catalase method; and TG using GPO Trinder reagent set with serum blank. The blood pressure (BP) of participants was measured by survey staff using an automated BP monitor, which was calibrated regularly [23].

Physical measurements

The anthropometric measurements were made at the testing sites by trained staff, and included height, weight and waist circumference. Height was measured without shoes using a stadiometer. Weight was measured without shoes and light clothing, using a digital weighting scale. Waist circumference was measured using a steel measuring tape. Body mass index (BMI) was calculated from the weight and height measurements [23].

Dietary and physical activity measurements

Dietary information was collected by multiple-pass 24 h diet recall using computer assisted telephone interviews. Dietary recall interviews were conducted by dietitians from the Department of Nutrition and Dietetics at Monash University in Melbourne, Australia. The FoodWorks® nutrition software (FoodWorks® Interview) were used for conducting the dietary recalls. Based on the three dietary recalls, the mean intake for each nutrient was calculated and used in the analysis [23]. Further information on the assessment of dietary intake data can be found in our previous publications [17, 23]. Physical activity information was collected via interviews with the participant. The time spent in physical activity was calculated based on the sum of the time spent walking or performing moderate activity plus double the time spent in vigorous activity (to indicate its greater intensity) [23, 25].

FPG and HbA1c

Fasting plasma glucose and HbA1c were the two dependent variables. The ADA cut-offs were used to identify those who were at a risk of T2DM based on FPG and HbA1c levels[22]. A binary variable was used to categorize subjects as being at low or high risk for T2DM: FPG <5.6 mmo/L (low risk, normal), vs. 5.6–6.9 mmol/L (high risk), and HbA1c <5.7% (low risk, normal) vs. 5.7–6.4% (high risk).

25OHD concentration

25OHD concentration was the primary independent variable. 25OHD concentration were categorized as tertiles: low 25OHD (median 33 nmol/L; range 10–44 nmol/l), medium 25OHD (median 54 nmol/L; range 45–65 nmol/L) and high 25OHD (median 77 nmol/L; range 65–204 nmol/L).

Statistical analysis

Socio-demographic factors.

In our analysis we considered a number of confounders, based on our [26] and others experience in the area [12, 27]. We adjusted for the following socio-demographic factors: age, gender, county of birth, Index of Relative Socio-economic Disadvantage (IRSED), physical activity, smoking status, and season. Age were entered as continuous variables. Country of birth was categorized as those born in Australia or overseas. The socio-economic indicator used was the IRSED, which is an index based on the social and economic conditions of individuals within an area [28]. Subjects were categorized into IRSED quartiles: quartile 1 (most disadvantaged), quartile 2 (disadvantaged), quartile 3 (less disadvantaged), and quartile 4 (least disadvantaged). Physical activity levels were classified into three categories: sufficient activity (≥150 minutes/week), insufficient activity (1–149 minutes/week), and inactive (0 minutes/week). There is a known seasonal variation to FPG and HbA1c, so we had to adjust for season [29, 30]. Season of biomedical assessment refers to the season of the year that the participant attended the testing site, and had their bloods collected for assessment of 25OHD status. Season of biomedical assessment were grouped as summer, autumn, winter, and spring. Smoking status were categorized into three categories: current smoker, ex-smoker, and non-smoker.

Dietary factors.

Based on previous research [31–33] dietary factors included in the analyses were: dietary fiber, magnesium, alcohol, calcium, zinc, carbohydrate intake, energy intake and under/over reporting of energy intake. Self-reported energy intake may result in the under or over reporting of true energy intake (EI), and this may confound the estimation of any diet and disease related outcomes [34]. We predicted basal metabolic rate (BMR) using the Henry/Oxford equations based on a range of ages (18–30, 30–60, 60–70, and 70 and above years), gender and body weight [35]. Rather than use the Goldberg cut-offs to identify under-reporters and over-reporters [36], we calculated the ratio of energy intake to BMR (EI:BMR) and treated it as a confounder. Dietary fiber, magnesium, alcohol, calcium, zinc, carbohydrate intake, energy intake, and under/over reporting of energy intake were all entered as continuous variables.

Biomedical factors.

Biomedical factors included in the analyses were: MetS components including waist circumference, high density lipoprotein (HDL) cholesterol, triglycerides (TG), BP, as well as body mass index (BMI), and haemoglobin levels (HbA1c model only). Haemoglobin levels were adjusted for in the HbA1c model only, as recent findings have indicated that haemoglobin levels may increase HbA1c levels [37, 38]. BMI and haemoglobin levels were entered as a continuous variable. Those with MetS tend to be at higher risk of developing T2DM, thus we adjusted for MetS components in the HbA1c and FPG model [39]. MetS components were each classified as binary variables, as defined by the joint interim statement [40]. Waist circumference were: normal waist circumference (<94cm for males or if Aboriginal or Torres Strait Islander (ATSI), Asian or South American <90cm; <80cm for females), or elevated waist circumference (≥ 94cm for males or if ATSI, Asian or South American ≥90cm; ≥80cm for females). HDL cholesterol were: normal HDL (≥1.0 mmol/L for males; ≥1.3 mmol/L for females), or low HDL (<1,0 mmol/l for males; <1.3 mmol/L for females). TG were: normal TG (<1.7 mmol/l), or hypertriglyceridaemia (≥1.7 mmol/l). BP were: normal BP (<130/85 mmHg and no anti-hypertensive medications), or high BP (≥130/85 mmHG or on anti-hypertensive medications) [40].

Statistical analysis.

The statistical analysis was conducted in the following two stages:

Step 1: Descriptive statistics for HbA1c and FPG were obtained and normality was assessed for variables of interest (natural logarithm transformation was applied if variable was skewed). Differences between groups were then examined by Independent samples t test and χ2 test.

Step 2: Multiple logistic regression analyses were employed to obtain adjusted odds ratio (AOR) and 95% CI for the associations between serum 25OHD and having higher HbA1c and FPG, respectively. Three categories of variables were used in the regression models including:socio-demographic variables (age, sex, country of birth, IRSED, physical activity, smoking status, and season), dietary factors (dietary fiber, magnesium, alcohol, calcium, zinc, carbohydrate intake, energy intake, and under/over reporting of energy intake) and biomedical variables (MetS components: waist circumference (ethnic specific cut-offs), HDL cholesterol (normal/low), TG (normal/high), BP (normal/high), BMI and haemoglobin levels for the HbA1c model only). All these variables have a known association with FPG and HbA1c and have been previously reported in the literature.

The primary research variable of interest, serum 25OHD, was entered into the multiple regression model as 1) a continuous variable (10 nmol/L increments) and also as 2) a categorical variable based on 25OHD tertiles. For both continuous and categorical 25OHD, the multiple regression model was initially adjusted for the socio-demographic variables in Model 1; secondly both socio-demographic variables, dietary factors and haemoglobin levels (HbA1c model only) in Model 2; and finally all socio-demographic, dietary factors and biomedical variables altogether in Model 3. The IBM SPSS Statistics for Windows, Version 21.0, was used for the statistical analyses. The VHM survey employed the use of the multistage stratified cluster-sampling procedure. In order to adjust for the unequal selection probability due to this sampling method, complex samples analysis was used. Complex samples approach takes into account the complex survey sampling and selection probability used in the VHM survey. Variables (strata variable, weighting variable and clustering variable) describing the survey design in terms of stratification, clustering and multistage sampling were entered into the SPSS complex samples approach for generating sampling weights in estimation and standard errors. A two-tailed p value of less than 0.05 was accepted as statistical significance.

Results

An overview of socio-demographic, dietary and clinical characteristics of subjects based on FPG and HbA1c levels are shown in Table 1 (please refer to S1 Table for the full socio-demographic, dietary and clinical characteristics of subjects). The prevalence of those with normal FPG (<5.6 mmol/L) was 84% and for high FPG (5.6–6.9 mmol/L) was 16%. 61% of the population had normal HbA1c (<5.6%), while 39% had high HbA1c (5.7–6.4%).

thumbnail
Download:
Table 1. Socio-demographic and clinical characteristics of participants by FPG and HbA1c.

https://doi.org/10.1371/journal.pone.0178825.t001

Association between serum 25OHD and FPG

When the serum 25OHD was entered as a continuous covariate to the multiple regression model, for every increment in serum 25OHD of 10 nmol/L, the odds of having higher FPG reduced by 9% (AOR 0.91, (0.86, 0.97); p = 0.002) after adjusting for socio-demographic variables, dietary factors, and biomedical variables in Model 3. For all models when serum 25OHD was entered as a categorical factor, compared with those people in the low 25OHD tertile, those with the high 25OHD tertile had a significantly reduced risk of higher FPG. More specifically, after adjustment for both socio-demographic, dietary factors and BMI in Model 2 the odds of having higher FPG reduced by 40% for those in high 25OHD tertile (AOR 0.60, (0.43, 0.83); p = 0.008) vs. low 25OHD tertile. After further adjustment for MetS components in Model 3, the AOR appeared relatively stable with a 39% reduced odds of higher FPG in the high 25OHD vs. low 25OHD tertile (AOR 0.61, (0.44, 0.84); p = 0.011) (Table 2).

thumbnail
Download:
Table 2. The association of serum 25OHD and FPG: Crude and adjusted odds ratio and their 95% CI based on logistic regression.

https://doi.org/10.1371/journal.pone.0178825.t002

Association between serum 25OHD and HbA1c

In Model 3 after adjustment for socio-demographic, dietary, and biomedical variables every 10 nmol/L increment in 25OHD also significantly reduced the odds of higher HbA1c by 6% (AOR 0.94, (0.90, 0.98); p = 0.009). In Model 2, after adjustment for socio-demographic variables, dietary factors, BMI and haemoglobin levels, there was a significantly reduced odds of having higher HbA1c by 33% in those with high vs. low 25OHD tertile. After further adjustment for MetS components in Model 3, a 26% reduced odds of higher HbA1c (AOR 0.74, (0.58, 0.93); p = 0.041) were found in the high 25OHD tertile group, compared to the low 25OHD tertile group (Table 3).

thumbnail
Download:
Table 3. The association of serum 25OHD and HbA1c: Crude and adjusted odds ratio and their 95% CI based on logistic regression.

https://doi.org/10.1371/journal.pone.0178825.t003

Discussion

The results from this population based study of adults from Victoria, Australia indicate that higher 25OHD levels were significantly related to a reduced risk of higher FPG and HbA1c levels. This significant inverse association persisted after the adjustment for a number of socio-demographic, dietary, and biomedical variables and MetS components. We found that those in the high 25OHD tertile had a 39% reduced risk of higher FPG and 26% reduced risk of higher HbA1c, when compared to the low 25OHD tertile. This was independent of MetS components, and haemoglobin levels in the HbA1c model, which have been found to confound the association between 25OHD and T2DM. Two recent meta-analyses [41–42] found no significant effect of vitamin D supplementation on IR [41], FPG [41, 42] or HbA1c [42]. Though, another found a beneficial effect of vitamin D supplementation on FPG and HbA1c, despite no effect on IR [43]. Our recent review summarised that although there were inverse associations between 25OHD and IR, the systematic reviews and meta-analysis in that review did not favour a casual role [44]. The present findings are in line with previous observational studies that found an inverse association between 25OHD status and high FPG levels [20, 45] and high HbA1c levels [46–51]. However, these studies adjusted for fewer variables, with smaller sample sizes than our study and were not population based [45, 49–51]. The differences in findings between both interventional and observational based studies indicates a need for good quality randomized controlled trials.

A number of inter-related factors contribute to the pathogenesis of T2DM, including the presence of MetS components. One prospective study found that those with IR and high FPG (5.6–6.9 mmol/L) had double the risk of worsened cardio-metabolic profile after nine years [52]. Obesity, dyslipidemia and hypertension are MetS components which are often found in those with pre-diabetes [53]. We adjusted for these variables in our analysis and found that the association between 25OHD and risk of T2DM still existed, irrespective of MetS. Other studies investigating 25OHD and risk of T2DM have adjusted for either none [47–49], or one MetS components [46, 54, 55], with one study accounting for all of the MetS components [8].

Potential mechanisms

The beneficial effect of 25OHD in reducing risk of T2DM is likely due to its effect on insulin action. The expression of the vitamin D receptor (VDR) in pancreatic beta cells indicates the importance of vitamin D in beta cell function [56, 57] and insulin secretion [5]. During times of vitamin D deficiency, beta cell function is blunted and insulin secretion is diminished [58]. Vitamin D may also indirectly influence insulin action via a calcium mediated effect. Vitamin D tightly regulates calcium homeostasis, whereby intracellular calcium levels are required to ensure effective action of insulin within different tissues [56]. Vitamin D is also involved in the regulation of the renin angiotensin system [59], endothelial vasodilation and lipid levels [60] which are mechanisms relating to the MetS components. T2DM and MetS are inter-twined, wherein IR appears to be a key player in the development of both conditions [61, 62]. Though there are commonalties in the mechanisms underlying both conditions, our study found that even on adjustment for MetS components, the association between the high 25OHD tertile and lower odds of higher FPG and HbA1c levels persisted. This may potentially suggest that vitamin D has a beneficial role in T2DM, independent of MetS. However, the cross-sectional nature of this study does not provide further insight into this observation.

Low grade chronic inflammation is a hallmark of many chronic diseases and may precede T2DM [63] possibly via initiation of IR [64]. Vitamin D has a role in immunity, and cellular studies show consistent reduction of inflammatory markers following cholecalciferol supplementation [65]. Adequate circulating 25OHD levels are required to obtain optimal anti-inflammatory responses in the body [66], especially in those tissues like immune cells, where the enzymes for conversion of 25OHD to 1,25OH2D are present. However, the optimum 25OHD levels for modulating inflammation responses are yet to be determined [67]. There is plausible but not confirmatory evidence to suggest the value is around 75–80 nmol/L [68, 69]; a point where maximum suppression of parathyroid hormone is also expected [70]. In support, a recent review [71], and emerging randomized controlled trials [72, 73] found beneficial effects of vitamin D supplementation on IR and fat loss in those individuals who reached this value over the period of the trial [72, 73]. In the present study those in the highest 25OHD tertile had a median of 77 nmol/L, where reduced risk of higher FPG and HbA1c was observed.

Limitations

The cross-sectional design does not afford causality of association, though we have controlled for several known confounders. The VHM did not collect information on supplement use, so we cannot separate the potential effects of vitamin D supplement and increased sun exposure on the higher 25OHD levels. Family history of T2DM may increase the risk of development of the disease [74], but unfortunately such information was not collected as part of the survey thus is a potential confounder. Approximately >75% of our sample were born in Australia, so potential effects of country of birth/ethnicity could not be determined.

Studies have indicated that 25OHD may vary due to the genetic variation of three polymorphisms in the vitamin D genes, including the vitamin D binding protein, VDR and the 25OHD activating enzyme [75, 76]. The concentration of the vitamin D binding protein may vary between individuals [77, 78] as well as fluctuations in its binding affinity [79]. There may also be a genetic association between VDR polymorphisms and T2DM. Recent evidence has found that the polymorphism of certain VDRs may increase susceptibility to T2DM [80] and certain MetS components [81]. Lastly, the conversion of 25OHD to its active form requires CYP27B1, an enzyme that may be affected by genetic variation [82]. These genetic variants were not studied as part of the VHM survey but certainly provide avenues for novel insights into the associations described here.

Strengths

This study used a representative sample of adult Australians from one state of the country. We adjusted for a wide range of socio-demographic, biomedical and dietary factors. Our study appears to be one of the few [8] which has accounted for all MetS components when investigating associations between 25OHD and T2DM. We also adjusted for the possible misreporting of energy intakes [35]. That being said 24 h recalls in this study were obtained from a five-pass method which is considered the gold standard for dietary information. We also adjusted the HbA1c analysis for haemoglobin concentrations, since concomitant anaemia may be associated with inaccurate HbA1c levels [37, 38]. Other studies in the area have not adjusted for this factor [8, 46–49, 54, 55].

Conclusions

Higher 25OHD status was associated with lower prevalent FPG as well as lower HbA1c concentrations after accounting for socio-demographic, lifestyle variables and MetS components. Such outcomes could suggest a direct role for the vitamin in the prevention of T2DM.

Supporting information

S1 Table. Socio-demographic and clinical characteristics of participants by FPG and HbA1c.

Data are presented as mean estimate (weighted) (%) for categorical variables, and mean estimate (weighted) and (SE) for continuous variables. Difference in the continuous and categorical variables between groups were assessed by independent samples t-test and Chi-square test, respectively. Legend: d, day; SE, standard error; min, minutes; wk, week.

https://doi.org/10.1371/journal.pone.0178825.s001

(DOCX)

Acknowledgments

MJS acknowledges the School of Public Health, Curtin University for infrastructure support, and the Victorian Department of Health and Human Services for use of the Victorian Health Monitor survey dataset. PKP is an Australian postgraduate research scholar. The opinions and analysis in this manuscript are those of the authors and not those of: the Department of Health and Human Services Victoria; the Victorian Government; the Secretary to the Department of Health Victoria or the Victorian Minister for Health.

Author Contributions

  1. Conceptualization: MJS.
  2. Formal analysis: PKP YZ.
  3. Methodology: PKP LSP MJS YZ ZA.
  4. Project administration: MJS.
  5. Software: PKP YZ.
  6. Supervision: MJS.
  7. Validation: YZ.
  8. Visualization: PKP LSP YZ ZA.
  9. Writing – original draft: PKP MJS.
  10. Writing – review & editing: LSP YZ ZA.

References

  1. 1. Lips P. Worldwide status of vitamin D nutrition. J Steroid Biochem Mol Biol. 2010;121(1–2):297–300. pmid:20197091
  2. 2. Daly RM, Gagnon C, Lu ZX, Magliano DJ, Dunstan DW, Sikaris KA, et al. Prevalence of vitamin D deficiency and its determinants in Australian adults aged 25 years and older: a national, population-based study. Clin Endocrinol (Oxf). 2012;77(1):26–35.
  3. 3. Nowson CA, McGrath JJ, Ebeling PR, Haikerwal A, Daly RM, Sanders KM, et al. Vitamin D and health in adults in Australia and New Zealand: a position statement. Med J Aust. 2012;196(11):686–687. pmid:22708765
  4. 4. Wild S, Roglic G, Green A, Sicree R, King H. Golbal prevalence of diabetes: estimates for the year 2000 and projections for 2030. Diabetes Care. 2004;27(5):1047–1053. pmid:15111519
  5. 5. Milner RDG, Hales CN. The role of calcium and magnesium in insulin secretion from rabbit pancreas studies in vitro. Diabetologia. 1967;3:47–49. pmid:4903553
  6. 6. Nagpal J, Pande JN, Bhartia A. A double-blind, randomized, placebo-controlled trial of the short-term effect of vitamin D3 supplementation on insulin sensitivity in apparently healthy, middle-aged, centrally obese men. Diabet Med. 2009;26(1):19–27. pmid:19125756
  7. 7. Mitri J, Dawson-Hughes B, Hu FB, Pittas AG. Effects of vitamin D and calcium supplementation on pancreatic beta cell function, insulin sensitivity, and glycemia in adults at high risk of diabetes: the Calcium and Vitamin D for Diabetes Mellitus (CaDDM) randomized controlled trial. Am J Clin Nutr. 2011;94(2):486–494. pmid:21715514
  8. 8. Gagnon C, Lu ZX, Magliano DJ, Dunstan DW, Shaw JE, Zimmet PZ, et al. Serum 25-hydroxyvitamin D, calcium intake, and risk of type 2 diabetes after 5 years: results from a national, population-based prospective study (the Australian Diabetes, Obesity and Lifestyle study). Diabetes Care. 2011;34(5):1133–1138. pmid:21430082
  9. 9. International Diabetes Federation. IDF diabetes atlas, 6th edn. Brussels, Belgium: International Diabetes Federation; 2013.
    • 10. Cespedes EM, Hu FB, Tinker L, Rosner B, Redline S, Garcia L, et al. Multiple healthful dietary patterns and type 2 diabetes in the Women's Health Initiative. Am J Epidemiol. 2016;183(7):622–633. pmid:26940115
    • 11. Ley SH, Hamdy O, Mohan V, Hu FB. Prevention and management of type 2 diabetes: dietary components and nutritional strategies. The Lancet. 2014;383(9933):1999–2007.
    • 12. Shamshirgaran SM, Jorm L, Bambrick H, Hennessy A. Independent roles of country of birth and socioeconomic status in the occurrence of type 2 diabetes. BMC Public Health. 2013.
    • 13. Wannamethee SG, Shaper AG, Lennon L, Morris RW. Metabolic syndrome vs Framingham risk score for prediction of coronary heart disease, stroke, and type 2 diabetes mellitus. Arch Intern Med. 2005;165:2644–2650. pmid:16344423
    • 14. Nichols GA, Moler EJ. Diabetes incidence for all possible combinations of metabolic syndrome components. Diabetes Res Clin Pract. 2010;90(1):115–121. pmid:20630609
    • 15. Lorenzo C, Okoloise M, Williams K, Stern MP, Haffner SM. The metabolic syndrome as a predictor of type 2 diabetes. Diabetes Care. 2003;26(11):3153–3159. pmid:14578254
    • 16. Bonora E, Targher G, Formentini G, Calcaterra F, Lombardi S, Marini F, et al. The metabolic syndrome is an independent predictor of cardiovascular disease in type 2 diabetic subjects. Prospective data from the Verona Diabetes Complications Study. Diabet Med. 2004;21(1):52–58. pmid:14706054
    • 17. Pannu PK, Zhao Y, Soares MJ, Piers LS, Ansari Z. The associations of vitamin D status and dietary calcium with the metabolic syndrome: an analysis of the Victorian Health Monitor survey. Public Health Nutr. 2016:1–12.
    • 18. Chon SJ, Yun BH, Jung YS, Cho SH, Choi YS, Kim SY, et al. Association between vitamin D status and risk of metabolic syndrome among Korean postmenopausal women. PLoS One. 2014;9(2):e89721. pmid:24586986
    • 19. Ford ES, Ajani UA, McGuire LC, Liu S. Concentrations of serum vitamin D and the metabolic syndrome among US adults. Diabetes Care. 2005;28(5):1128–1230.
    • 20. Gagnon C, Lu ZX, Magliano DJ, Dunstan DW, Shaw JE, Zimmet PZ, et al. Low serum 25-hydroxyvitamin D is associated with increased risk of the development of the metabolic syndrome at five years: results from a national, population-based prospective study (The Australian Diabetes, Obesity and Lifestyle Study: AusDiab). J Clin Endocrinol Metab. 2012;97(6):1953–1961. pmid:22442263
    • 21. Vitezova A, Zillikens MC, van Herpt TT, Sijbrands EJ, Hofman A, Uitterlinden AG, et al. Vitamin D status and metabolic syndrome in the elderly: the Rotterdam Study. Eur J Endocrinol. 2015;172(3):327–335. pmid:25468955
    • 22. American Diabetes Association. Diagnosis and classification of diabetes mellitus. Diabetes Care. 2010;33 (S1):S62–69.
    • 23. Department of Health. The Victorian health monitor. Melbourne, Victoria: State Government of Victoria; 2012.
      • 24. Department of Health. The Victorian health monitor food and nutrition report. Melbourne, Victoria: State Government of Victoria; 2012.
        • 25. Armstrong T, Bauman A, Davies J. Physical Activity Patterns of Australian Adults: Results of the 1999 National Physical Activity Survey. Canberra, ACT: Australian Institute of Health and Welfare; 2000.
          • 26. Markwick A, Ansari Z, Sullivan M, McNeil J. Social determinants and psychological distress among Aboriginal and Torres Strait islander adults in the Australian state of Victoria: a cross-sectional population based study. Soc Sci Med. 2015;128:178–187. pmid:25618607
          • 27. Krishnan S, Rosenberg L, Palmer JR. Physical activity and television watching in relation to risk of type 2 diabetes: the Black Women's Health Study. Am J Epidemiol. 2009;169(4):428–434. pmid:19056835
          • 28. Australian Bureau of Statistics. Socio-economic indexes for areas (SEIFA). Canberra, ACT: Australian Bureau of Statistics; 2011.
            • 29. Gikas A, Sotiropoulos A, Pastromas V, Papazafiropoulou A, Apostolou O, Pappas S. Seasonal variation in fasting glucose and HbA1c in patients with type 2 diabetes. Prim Care Diabetes. 2009;3(2):111–114. pmid:19535310
            • 30. Tseng CL, Brimacombe M, Xie M, Rajan M, Wang H, Kolassa J, et al. Seasonal patterns in monthly hemoglobin A1c values. Am J Epidemiol. 2005;161(6):565–574. pmid:15746473
            • 31. Dong JY, Xun P, He K, Qin LQ. Magnesium intake and risk of type 2 diabetes: meta-analysis of prospective cohort studies. Diabetes Care. 2011;34(9):2116–2122. pmid:21868780
            • 32. Yang L, Shu L, Jiang J, Qiu H, Zhao G, Zhou Y, et al. Long-term effect of dietary fibre intake on glycosylated haemoglobin A1c level and glycaemic control status among Chinese patients with type 2 diabetes mellitus. Public Health Nutr. 2014;17(8):1858–1864. pmid:23883570
            • 33. Mitri J, Muraru MD, Pittas AG. Vitamin D and type 2 diabetes: a systematic review. Eur J Clin Nutr. 2011;65(9):1005–1015. pmid:21731035
            • 34. Goldberg GR, Black AE. Assessment of the validity of reported enery intakes—review and recent developments. Scandinavian Journal of Nutrition. 1998;42:6–9.
            • 35. Henry CJK. Basal metabolic rate studies in humans: measurement and development of new equations. Public Health Nutr. 2007;8(7a):1133–1152.
            • 36. Goldberg GR, Black AE, Jebb SA, Coles TJ, Murgatroyd PR, Coward WA, et al. Critical evaluation of energy intake data using fundamental principles of energy physioology: 1. derivation of cut-off limits to identidy under-recording. Eur J Clin Nutr. 1991;45:569–581. pmid:1810719
            • 37. Attard SM, Herring AH, Wang H, Howard AG, Thompson AL, Adair LS, et al. Implications of iron deficiency/anemia on the classification of diabetes using HbA1c. Nutr Diabetes. 2015;5:e166. pmid:26098445
            • 38. Kim C, Bullard KM, Herman WH, Beckles GL. Association between iron deficiency and hba1c levels among adults without diabetes in the National Health and Nutrition Examination Survey, 1999–2006. Diabetes Care. 2010;33(4):780–785. pmid:20067959
            • 39. Grundy SM. Pre-diabetes, metabolic syndrome, and cardiovascular risk. J Am Coll Cardiol. 2012;59(7):635–643. pmid:22322078
            • 40. Alberti KG, Eckel RH, Grundy SM, Zimmet PZ, Cleeman JI, Donato KA, et al. Harmonizing the metabolic syndrome: a joint interim statement of the International Diabetes Federation Task Force on Epidemiology and Prevention; National Heart, Lung, and Blood Institute; American Heart Association; World Heart Federation; International Atherosclerosis Society; and International Association for the Study of Obesity. Circulation. 2009;120(16):1640–1645. pmid:19805654
            • 41. Jamka M, Woźniewicz M, Jeszka J, Mardas M, Bogdański P, Stelmach-Mardas M. The effect of vitamin D supplementation on insulin and glucose metabolism in overweight and obese individuals: systematic review with meta-analysis. Sci Rep. 2015;176(1):R1–R14.
            • 42. Krul-Poel YH, Ter Wee MM, Lips P, Simsek S. The effect of vitamin D supplementation on glycaemic control in patients with type 2 diabetes mellitus: a systematic review and meta-analysis. Eur J Endocrinol. 2017;176(1):R1–R14. pmid:27484453
            • 43. Poolsup N, Suksomboon N, Plordplong N. Effect of vitamin D supplementation on insulin resistance and glycaemic control in prediabetes: a systematic review and meta-analysis. Diabet Med. 2016;33(3):290–299. pmid:26308752
            • 44. Soares MJ, Pannu PK, Calton EK, Reid CM, Hills AP. Vitamin D status and calcium intake in systemic inflammation, insulin resistance and the metabolic syndrome: an update on current evidence. Trends in Food Sci Technol. 2017;62:79–90.
            • 45. Forouhi NG, Luan J, Cooper A, Boucher BJ, Wareham NJ. Baseline serum 25-hydroxy vitamin D is predictive of future glycemic status and insulin resistance: the Medical Research Council Ely Prospective Study 1990–2000. Diabetes. 2008;57(10):2619–2625. pmid:18591391
            • 46. Hutchinson MS, Figenschau Y, Njolstad I, Schirmer H, Jorde R. Serum 25-hydroxyvitamin D levels are inversely associated with glycated haemoglobin (HbA(1c). The Tromso Study. Scand J Clin Lab Invest. 2011;71(5):399–406. pmid:21517716
            • 47. Hypponen E, Power C. Vitamin D status and glucose homeostasis in the 1958 British birth cohort: the role of obesity. Diabetes Care. 2006;29(10):2244–2246. pmid:17003300
            • 48. Kositsawat J, Freeman VL, Gerber BS, Geraci S. Association of A1C levels with vitamin D status in U.S. adults: data from the National Health and Nutrition Examination Survey. Diabetes Care. 2010;33(6):1236–1238. pmid:20215453
            • 49. Manickam B, Neagu V, Kukreja SC, Barengolts E. Relationship between glycated hemoglobin and circulating 25-hydroxyvitamin D concentration in African American and Caucasian American men. Endocr Pract. 2013;19(1):73–80. pmid:23186960
            • 50. Need AG, O'Loughlin PD, Horowitz M, Nordin BE. Relationship between fasting serum glucose, age, body mass index and serum 25 hydroxyvitamin D in postmenopausal women. Clin Endocrinol (Oxf). 2005;62(6):738–741.
            • 51. Zoppini G, Galletti A, Targher G, Brangani C, Pichiri I, Negri C, et al. Glycated haemoglobin is inversely related to serum vitamin D levels in type 2 diabetic patients. PLoS One. 2013;8(12):e82733. pmid:24358222
            • 52. Salazar MR, Carbajal HA, Espeche WG, Aizpurua M, Leiva Sisnieguez CE, Leiva Sisnieguez BC, et al. Insulin resistance: the linchpin between prediabetes and cardiovascular disease. Diab Vasc Dis Res. 2016;13(2):157–163. pmid:26802220
            • 53. Ferrannini E, Gastaldelli A, Iozzo P. Pathophysiology of prediabetes. Med Clin North Am. 2011;95(2):327–339, vii-viii. pmid:21281836
            • 54. Mattila C, Knekt P, Mannisto S, Rissanen H, Laaksonen MA, Montonen J, et al. Serum 25-hydroxyvitamin D concentration and subsequent risk of type 2 diabetes. Diabetes Care. 2007;30(10):2569–2570. pmid:17626891
            • 55. Pittas AG, Dawson-Hughes B, Li T, Van Dam RM, Willett WC, Manson JE, et al. Vitamin D and calcium intake in relation to type 2 diabetes in women. Diabetes Care. 2006;29(3):650–656. pmid:16505521
            • 56. Cavalier E, Delanaye P, Souberbielle JC, Radermecker RP. Vitamin D and type 2 diabetes mellitus: where do we stand? Diabetes Metab. 2011;37(4):265–272. pmid:21345709
            • 57. Norman AW. Intestinal calcium absorption: a vitamin D-hormone-mediated adaptive response. Am J Clin Nutr. 1990;51:290–300. pmid:2407100
            • 58. Boucher BJ. Inadequate vitamin D status: does it contribute to the disorders comprising syndrome ‘X’? Br J Nutr. 1998;79:315–327. pmid:9624222
            • 59. Al Mheid I, Patel RS, Tangpricha V, Quyyumi AA. Vitamin D and cardiovascular disease: is the evidence solid? Eur Heart J. 2013;34(48):3691–3698. pmid:23751422
            • 60. Dalan R, Liew H, Tan WKA, Chew DEK, Leow MK-S. Vitamin D and the endothelium: basic, translational and clinical research updates. IJC Metabolic & Endocrine. 2014;4:4–17.
            • 61. Ferrannini E, Haffner SM, Mitchell BD, Stern MP. Hyperinsulinaemia: the key feature of a cardiovascular and metabolic syndrome. Diabetologia. 1991;34:416–422. pmid:1884900
            • 62. Reaven GM. Role of insulin resistance in human disease. Diabetes. 1988;37:1595–1607. pmid:3056758
            • 63. Pickup JC. Inflammation and activated innate immunity in the pathgenesis of type 2 diabetes. Diabetes Care. 2004;27:813–823. pmid:14988310
            • 64. Keane KN, Cruzat VF, Carlessi R, de Bittencourt PI Jr., Newsholme P. Molecular events linking oxidative stress and inflammation to insulin resistance and beta-cell dysfunction. Oxid Med Cell Longev. 2015; 2015.
            • 65. Calton EK, Keane KN, Newsholme P, Soares MJ. The impact of vitamin D levels on inflammatory status: a systematic review of immune cell studies. PLoS One. 2015;10(11):e0141770. pmid:26528817
            • 66. Zhang Y, Leung DY, Richers BN, Liu Y, Remigio LK, Riches DW, et al. Vitamin D inhibits monocyte/macrophage proinflammatory cytokine production by targeting MAPK phosphatase-1. J Immunol. 2012;188(5):2127–2135. pmid:22301548
            • 67. Bikle DD. Vitamin D: newly discovered actions require reconsideration of physiologic requirements. Trends Endocrinol Metab. 2010;21(6):375–384. pmid:20149679
            • 68. Heaney RP. Vitamin D in health and disease. Clin J Am Soc Nephrol. 2008;3(5):1535–1541. pmid:18525006
            • 69. Vieth R. Why the minimum desirable serum 25-hydroxyvitamin D level should be 75 nmol/L (30 ng/ml). Best Pract Res Clin Endocrinol Metab. 2011;25(4):681–691. pmid:21872808
            • 70. Durazo-Arvizu RA, Dawson-Hughes B, Sempos CT, Yetley EA, Looker AC, Cao G, et al. Three-phase model harmonizes estimates of the maximal suppression of parathyroid hormone by 25-hydroxyvitamin D in persons 65 years of age and older. J Nutr. 2010;140:595–599. pmid:20089790
            • 71. Pannu PK, Calton EK, Soares MJ. Calcium and vitamin D in obesity and related chronic disease. In: Henry J, editor. Advances in Food and Nutrition Research, volume 77. London: Academic Press; 2016. pp. 57–88. https://doi.org/10.1016/bs.afnr.2015.11.001 pmid:26944102
              • 72. Rosenblum JL, Castro VM, Moore CE, Kaplan LM. Calcium and vitamin D supplementation is associated with decreased abdominal visceral adipose tissue in overweight and obese adults. Am J Clin Nutr. 2012;95(1):101–108. pmid:22170363
              • 73. von Hurst PR, Stonehouse W, Coad J. Vitamin D supplementation reduces insulin resistance in South Asian women living in New Zealand who are insulin resistant and vitamin D deficient—a randomised, placebo-controlled trial. Br J Nutr. 2010;103(4):549–555. pmid:19781131
              • 74. InterAct C, Scott RA, Langenberg C, Sharp SJ, Franks PW, Rolandsson O, et al. The link between family history and risk of type 2 diabetes is not explained by anthropometric, lifestyle or genetic risk factors: the EPIC-InterAct study. Diabetologia. 2013;56(1):60–69. pmid:23052052
              • 75. Wang TJ, Zhang F, Richards JB, Kestenbaum B, van Meurs JB, Berry D, et al. Common genetic determinants of vitamin D insufficiency: a genome-wide association study. The Lancet. 2010;376(9736):180–188.
              • 76. Boucher BJ. Is vitamin D status relevant to metabolic syndrome? Dermatoendocrinol. 2012;4(2):212–224. pmid:22928079
              • 77. Fu L, Yun F, Oczak M, Wong BY, Vieth R, Cole DE. Common genetic variants of the vitamin D binding protein (DBP) predict differences in response of serum 25-hydroxyvitamin D [25(OH)D] to vitamin D supplementation. Clin Biochem. 2009;42(10–11):1174–1177. pmid:19302999
              • 78. Lauridsen AL, Vestergaard P, Hermann AP, Brot C, Heickendorff L, Mosekilde L, et al. Plasma concentrations of 25-hydroxy-vitamin D and 1,25-dihydroxy-vitamin D are related to the phenotype of Gc (vitamin D-binding protein): a cross-sectional study on 595 early postmenopausal women. Calcif Tissue Int. 2005;77(1):15–22. pmid:15868280
              • 79. Arnaud J, Constans J. Affinity differences for vitamin D metabolites associated with the genetic isoforms of the human serum carrier protein (DBP). Hum Genet. 1993;92:183–188. pmid:8370586
              • 80. Li L, Wu B, Liu JY, Yang LB. Vitamin D receptor gene polymorphisms and type 2 diabetes: a meta-analysis. Arch Med Res. 2013;44(3):235–241. pmid:23506721
              • 81. Martinez-Hernandez A, Cordova EJ, Rosillo-Salazar O, Garcia-Ortiz H, Contreras-Cubas C, Islas-Andrade S, et al. Association of HMOX1 and NQO1 polymorphisms with metabolic syndrome components. PLoS One. 2015;10(5):e0123313. pmid:25933176
              • 82. Hypponen E, Berry DJ, Wjst M, Power C. Serum 25-hydroxyvitamin D and IgE—a significant but nonlinear relationship. Allergy. 2009;64(4):613–620. pmid:19154546