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

The associations between VDR BsmI polymorphisms and risk of vitamin D deficiency, obesity and insulin resistance in adolescents residing in a tropical country

  • Rayinda Rahmadhani ,

    Contributed equally to this work with: Rayinda Rahmadhani, Nur Lisa Zaharan

    Affiliation Department of Pharmacology, Faculty of Medicine, University of Malaya, Kuala Lumpur, Malaysia

    ⨯
  • Nur Lisa Zaharan ,

    Contributed equally to this work with: Rayinda Rahmadhani, Nur Lisa Zaharan

    nurlisazaharan@um.edu.my

    Affiliation Department of Pharmacology, Faculty of Medicine, University of Malaya, Kuala Lumpur, Malaysia

    ⨯
  • Zahurin Mohamed,

    Affiliation Department of Pharmacology, Faculty of Medicine, University of Malaya, Kuala Lumpur, Malaysia

    ⨯
  • Foong Ming Moy,

    Affiliation Julius Centre University of Malaya, Department of Social and Preventive Medicine, Faculty of Medicine, University of Malaya, Kuala Lumpur, Malaysia

    ⨯
  • Muhammad Yazid Jalaludin

    Affiliation Department of Paediatrics, Faculty of Medicine, University of Malaya, Kuala Lumpur, Malaysia

    ⨯

Abstract

Background

The vitamin D receptor (VDR) gene is expressed abundantly in different tissues; including adipocytes and pancreatic beta cells. The rs1544410 or BsmI single nucleotide polymorphism (SNP) in the intronic region of the VDR gene has been previously associated with vitamin D levels, obesity and insulin resistance.

Aims

This study was aimed to examine the association between BsmI polymorphism and risk of vitamin D deficiency, obesity and insulin resistance in adolescents living in a tropical country.

Methods

Thirteen-year-old adolescents were recruited via multistage sampling from twenty-three randomly selected schools across the city of Kuala Lumpur, Malaysia (n = 941). Anthropometric measurements were obtained. Obesity was defined as body mass index higher than the 95th percentile of the WHO chart. Levels of fasting serum vitamin D (25-hydroxyvitamin D (25(OH)D)), glucose and insulin were measured. HOMA-IR was calculated as an indicator for insulin resistance. Genotyping was performed using the Sequenom MassARRAY platform (n = 807). The associations between BsmI and vitamin D, anthropometric parameters and HOMA-IR were examined using analysis of covariance and logistic regression.

Result

Those with AA genotype of BsmI had significantly lower levels of 25(OH)D (p = 0.001) compared to other genotypes. No significant differences was found across genotypes for obesity parameters. The AA genotype was associated with higher risk of vitamin D deficiency (p = 0.03) and insulin resistance (p = 0.03) compared to GG. The A allele was significantly associated with increased risk of vitamin D deficiency compared to G allele (adjusted odds ratio (OR) = 1.63 (95% Confidence Interval (CI) 1.03–2.59, p = 0.04). In those with concurrent vitamin D deficiency, having an A allele significantly increased their risk of having insulin resistance compared to G allele (adjusted OR = 2.66 (95% CI 1.36–5.19, p = 0.004).

Conclusion

VDR BsmI polymorphism was significantly associated with vitamin D deficiency and insulin resistance, but not with obesity in this population.

Citation: Rahmadhani R, Zaharan NL, Mohamed Z, Moy FM, Jalaludin MY (2017) The associations between VDR BsmI polymorphisms and risk of vitamin D deficiency, obesity and insulin resistance in adolescents residing in a tropical country. PLoS ONE 12(6): e0178695. https://doi.org/10.1371/journal.pone.0178695

Editor: Zhiming Zhu, Third Military Medical University Daping Hospital and Research Institute of Surgery, CHINA

Received: August 3, 2016; Accepted: May 17, 2017; Published: June 15, 2017

Copyright: © 2017 Rahmadhani 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 availability for this study is subjected to ethical approval from the Ethics Committee of the University of Malaya Medical Centre. Data will be available upon request to all interested researchers. Request of data can be made to the corresponding author, Dr Nur Lisa Zaharan and the; Ethics Committee of the University of Malaya Medical Centre, University Malaya Medical Centre, Lembah Pantai, 59100 Kuala Lumpur email: ummc@ummc.edu.my.

Funding: This study is funded under the High Impact Research, Ministry of Higher Education Grant, Malaysia (HIR H00021-00-E000082) and the University Malaya Research Grant (RG310-11HTM). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

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

Introduction

The vitamin D endocrine system takes part in various biological processes including musculoskeletal development, erythropoiesis and blood pressure regulation [1]. Yet, it was estimated that one billion people had deficient or insufficient levels of vitamin D based on population studies [2]. The prevalence of vitamin D deficiency varies between different regions of the world whilst showing seasonal variations [3–7]. Factors that have been associated with increased risk of vitamin D deficiency include low skin exposure to sunlight, low dietary intake of vitamin D, high body mass index (BMI), and genetic predispositions [8–12].

Interestingly, approximately two thirds of adults in Malaysia, a country in South East Asia, had vitamin D deficiency despite receiving sunlight all-year round (it is located at 3 degree north of the equator) [13]. More worryingly, a study on thirteen-year old Malaysian adolescents revealed that nearly eighty percent of them fulfilled the criteria for vitamin D deficiency [14]. In both Malaysian adults and adolescents, adiposity was associated with vitamin D deficiency [13, 14]. Furthermore, the prevalence of obesity among children and adolescents in this country was also rather disturbing, being the highest in the South East Asian region with boys at 22% and girls at 19% [15]. Studies that examined factors that may contribute to the relationship between vitamin D levels and obesity are therefore required to further understand the possible mechanisms that may link them both, especially in our young population.

In addition to environmental contributions, genetic factors may account for 23% to 80% of variability in serum vitamin D levels, as observed in twin studies [16, 17]. There is, however, currently a dearth in studies that examined possible genetic contributors to vitamin D deficiency in the Malaysian population. Amongst single nucleotide polymorphisms (SNP) that may be associated with vitamin D levels are those located on the vitamin D receptor (VDR) gene (18). The VDR gene plays a crucial role in the modulation of vitamin D pathways and regulation of hormone responsive genes [18]. Interestingly, the VDR gene is also expressed in adipocytes and pancreatic beta cells and thus may influence body composition by directly regulating the differentiation and metabolism of adipocytes; or indirectly by insulin modulation [19, 20]. One of the SNPs in the VDR gene associated with vitamin D levels is the rs1544410 SNP, located in the intronic region (intron 8 near the 3´end). This SNP is a restriction fragment length polymorphism of the restriction endonucleases BsmI [21] and is also commonly known as the BsmI polymorphism. It is thought to affect VDR translational activity due to its strong linkage disequilibrium with a polyadenosine (poly (A)) microsatellite repeat in the 3´untranslated region [22]. In addition to vitamin D levels, the BsmI polymorphism has also been shown to be associated with obesity, insulin resistance and type 2 diabetes in some population [23–25]. Unfortunately, very few studies examined the relative influence of this polymorphism on vitamin D levels, adiposity and insulin resistance in the paediatric population [26–28]. Thus, this cross-sectional study was aimed to examine the association between BsmI (rs1544410) VDR gene polymorphism with vitamin D deficiency, adiposity and insulin resistance in our adolescents.

Methodology

Study participants

Ethical approval for this study was obtained from the Ethics Committee of the University of Malaya Medical Centre (MEC 896.123). The permission to carry out this study was obtained from the Ministry of Education and the principals of the respective participating schools. Thirteen-year-old participants were recruited via multistage sampling from 23 randomly selected government-funded secondary schools across the city of Kuala Lumpur, Malaysia from January 2012 to July 2012. Written informed consent was obtained by the researchers from either parents or guardians of each of the participants prior to the study.

Anthropometric and pubertal assessments

Anthropometric measurements were taken by trained researchers following a standard protocol. Body weight of the participants was measured with participants in light clothing and shoes removed using a digital calibrated floor scale (SECA 813; Seca GmbH&Co., Hamburg, Germany) to the nearest 0.1kg. Height with shoes removed was measured using a portable stadiometer (SECA 813; Seca GmbH&Co., Hamburg, Germany) to the nearest 0.1cm. BMI was calculated as baseline body weight divided by height squared (kg/m2). Waist circumference (WC) was measured at the midway between iliac crest and the 10th rib, according to WHO STEPS protocol [29]. Hip circumference (HC) was measured at the widest extension of the buttocks. Both WC and HC were measured using inelastic measurement tape (SECA 203, Hamburg, Germany). Waist-Hip Ratio (WHR) was calculated as WC divided by HC (WC/HC). Body fat percentage (BF%) was measured using a validated portable body composition analyzer (Inbody 230; Biospace Co.Ltd, Seoul, Korea) according to the instruction manual provided. Pubertal stage was self-assessed using a coloured Tanner stages illustration [30].

Vitamin D level measurement

Fasting venous blood (after at least 8 hours of fasting) was drawn from the participants by trained physicians in a sitting position between 8 am to 9 am in the morning. Vitamin D levels (25-hydroxyvitamin D (25(OH)D)) were determined by an accredited clinical diagnostic laboratory (CDL) at the University of Malaya Medical Centre (UMMC), Kuala Lumpur, Malaysia adhering to standard clinical laboratory protocols. The kits used were Elecsys® Vitamin D total assay (Cobas, Roche Diagnostics Limited, Switzerland) which implements an electro-chemiluminescence binding assay (ECLIA) to determine the total 25(OH)D in-vitro. Measurements were conducted according to manufacturer’s protocol.

Measures of insulin resistance

Fasting blood glucose (FBG) and fasting insulin (FI) were measured by the clinical diagnostic laboratory as mentioned above. Insulin resistance was calculated using a homeostatic model assessment of insulin resistance (HOMA-IR); HOMA-IR = FI (μU/mL) x FBG (mmol/L)/ 22.5.

Genotyping

Genomic DNA from blood samples was extracted using QiAmp DNA Mini Kit (Qiagen, Hilden, Germany). The quality of the extracted DNA was such that the absorbance ratios of at least 1.8 was attained for both 260/280 and 260/230 readings. All DNA samples and duplicates were diluted to 10ng/μL and 20ng/μL, respectively, before being transferred to the respective wells. Quality controls included a blank and five duplicates. The VDR SNP of rs1544410 (BsmI) was genotyped using a Sequenom MassARRAY platform with iPLEX GOLD chemistry (Sequenom, San Diego, California) following the manufacturer’s protocols. The MassARRAY system was created based on the technology of MALDI-TOF (matrix-assisted laser desorption/ionization time-of-flight) mass spectrometry. The genotyping call rate was >95% for the variant.

Statistical analyses

There were 941 participants, all of which had complete anthropometric data. Vitamin D levels were available in 678 participants, fasting insulin for 795 participants and genotyping for 807 participants. Obesity was defined as BMI higher than the 95th percentile according to the WHO 2007 Growth Reference Data for 5–19 years [31]. Vitamin D status was categorized into three groups according to Misra et al [10] whereby vitamin D deficiency was defined as those with a 25(OH)-D level of ≤15 ng/mL (≤37.5 nmol/L), vitamin D insufficiency as those with a 25(OH)-D level between 15 ng/mL and 20 ng/mL (37.5–50 nmol/L) and vitamin D sufficiency as those with a 25(OH)-D level of ≥20 ng/mL (≥50 nmol/L). Insulin resistance was categorised based on HOMA-IR values as proposed by Yin et al [32], in which different cut-offs were used based on pubertal stages. A HOMA-IR >2.6 was used to define insulin resistance in pre-pubertal adolescents, while for pubertal adolescents the cut-off was a HOMA-IR >3.2.

Continuous variables were presented as means with standard deviations. All variables were tested for normal distribution using the Kolmogorov-Smirnov test. The Mann-Whitney U test was used to compare demographic and anthropometric measurements, vitamin D levels and insulin resistance between boys and girls. Categorical variables were presented as frequency and percentages and comparisons between genders were performed using chi-square testing. Correlation analyses were performed using Pearson’s correlation to examine the relationship between vitamin D levels and BMI as well as fasting insulin. In addition, binary logistic regression was used to examine the risk of obesity and insulin resistance in those with vitamin D deficiency compared to those with sufficient levels of vitamin D. Gender, ethnicities, maternal education and pubertal stage (and BMI for insulin resistance) were included as covariates. Maternal education was used as a proxy for socio-economic status [33].

Genotype distribution was assessed for Hardy-Weinberg equilibrium (HWE) by using the χ2 test. A p-value of more than 0.05 signifies agreement with HWE. Comparisons of adiposity parameters, 25(OH)D levels and metabolic status between genotypes (GG, GA and AA) were made using an analysis of covariance (ANCOVA) general linear model. The associations between VDR BsmI polymorphism (both genotypes GG, GA and AA as well as alleles G and A) and risk of vitamin D deficiency, obesity and insulin resistance were examined using binary logistic regression; adjusting for covariates mentioned above. In addition, the risks of obesity and insulin resistance as stratified according to vitamin D status comparing A allele (mutant) to the G allele (wild) were examined using binary logistic regression. A p-value of less than 0.05 was considered as statistically significant. Statistical analysis was performed using SPSS 24.0 software (IBM SPSS Statistics).

Results

Participants’ demographics, anthropometric data, vitamin D status and measures of insulin resistance

Overall, there were 941 (13 year old) participants, of which the majority were girls (72%) and of Malay ethnicity (75%) as presented in Table 1. The prevalence of obesity in our adolescent population was 21%, with 19% of girls and 29% of boys were categorised as obese while the prevalence of morbid obesity was 11%. The proportion of those considered as having insulin resistance was similar to obesity at 21%, with no significant differences found between boys and girls (20% vs 21% respectively). Of the 648 participants with 25(OH)D levels, 45% of them fulfilled the criteria for vitamin D deficiency (12% boys and 56% girls).

thumbnail
Download:
Table 1. Demographic, anthropometric and clinical profiles of thirteen-year old Malaysian adolescents in comparing boys and girls.

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

The majority of girls had already reached puberty compared to boys (p = 0.03). Comparing girls and boys in terms of adiposity parameters, boys have significantly higher waist circumferences and waist hip ratios while girls have significantly higher percentages of body fat (p<0.0001). In terms of biochemistry profiles, girls were found to have significantly lower 25(OH)D levels (p<0.0001) with higher levels of fasting insulin (p<0.0001) as compared to boys (Table 1).

The 25(OH)D levels were negatively correlated with both BMI (r = -0.130, p = 0.002) and fasting insulin (r = -0.127, p = 0.02), albeit weakly. Those with vitamin D deficiency was associated with increased risk of being obese (adjusted OR = 3.12, 95% CI 1.37–7.11, p = 0.007) but not of insulin resistance (adjusted OR = 1.78, 95% CI 0.80–3.97, p = 0.16) compared to those with sufficient levels after adjusting for covariates such as gender, ethnicity, maternal education, puberty status and BMI (for insulin resistance).

The association between BsmI polymorphism and anthropometric measurements, vitamin D levels and measure of insulin resistance

The frequencies of GG, GA and AA genotypes of VDR BsmI SNP in our adolescent population met the criteria for the Hardy Weinberg Equilibrium. Those with a GG genotype had significantly higher levels of 25(OH)D, while those with the AA genotype had higher levels of fasting insulin and HOMA-IR values as compared to other genotypes (Table 2). However, after adjustment of covariates, there was no significant differences between genotypes in terms of fasting insulin and HOMA-IR. In addition, no significant differences in adiposity parameters were found across the three genotypes.

thumbnail
Download:
Table 2. Anthropometric measurements, vitamin D Levels and measures of insulin resistance according to rs1544410 (BsmI) genotypes.

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

Those with an AA genotype of VDR BsmI were significantly associated with increased risk of vitamin D deficiency compared to the GG genotype (OR = 8.37 (95% CI 1.07, 65.71)) as presented in Table 3. No significant association was observed between BsmI polymorphism and risk of obesity in our adolescents (Table 4). The AA genotype was also found to confer a higher risk of insulin resistance compared to GG (OR = 2.75 (95% CI 1.13, 6.67)) as shown in Table 5.

thumbnail
Download:
Table 3. Association between VDR BsmI and risk of vitamin D deficiency in Malaysian adolescents presented as OR (unadjusted and adjusted) with 95% CI.

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

thumbnail
Download:
Table 4. Association between VDR BsmI and risk of obesity in Malaysian adolescents presented as OR (unadjusted and adjusted) with 95% CI.

https://doi.org/10.1371/journal.pone.0178695.t004

thumbnail
Download:
Table 5. Association between VDR BsmI and risk of insulin resistance in Malaysian adolescents presented as OR (unadjusted and adjusted) with 95% CI.

https://doi.org/10.1371/journal.pone.0178695.t005

Having the A allele of the VDR BsmI was significantly associated with increased risk of vitamin D deficiency in the Malaysian adolescent compared to the G allele (adjusted OR = 1.63 (95% CI 1.03, 2.59, p = 0.04) as presented in Table 6. No significant increased risk of obesity and insulin resistance was found with the A allele compared to the G allele in this population. However, when stratified according to vitamin D status, it was shown that in those with vitamin D deficiency, having an A allele significantly increased their risk of having insulin resistance compared to G allele (adjusted OR = 2.66 (95% CI 1.36, 5.19, p = 0.004) (Table 7).

thumbnail
Download:
Table 6. Risk of vitamin D deficiency, obesity and insulin resistance with A allele of VDR BsmI compared to G allele in Malaysian Adolescent Population presented as OR (unadjusted and adjusted) with 95% CI.

https://doi.org/10.1371/journal.pone.0178695.t006

thumbnail
Download:
Table 7. Risk of obesity and insulin resistance with A allele of VDR BsmI polymorphism compared to G allele in Malaysian adolescents when stratified according to Vitamin D status presented as adjusted OR with 95% CI.

https://doi.org/10.1371/journal.pone.0178695.t007

Discussion

Results of this study indicated that more than a fifth of Malaysian adolescents in the capital city of Kuala Lumpur, were obese. More worryingly, a tenth were considered morbidly obese while a fifth were demonstrated to have insulin resistance. The percentage of those with vitamin D deficiency was forty-five percent with girls found to have significantly lower 25(OH)D levels compared to boys. We found that adolescents carrying the AA genotype of VDR BsmI were associated with increased risk of both vitamin D deficiency and insulin resistance compared to the GG genotype. The A allele of the BsmI was associated with significantly higher risk of vitamin D deficiency compared to the G allele. Moreover, in those with concurrent vitamin D deficiency, having an A allele was significantly associated with increased risk of insulin resistance in this population. Although vitamin D deficiency was associated with increased risk of obesity, we were not able to demonstrate any significant association between BsmI polymorphisms and adiposity parameters or risk of obesity in our adolescent population.

The prevalence of obesity in our study was comparable to the recent multi-country study by Marie Ng et al [15], although higher than the 9% prevalence documented from the Malaysian Health and Adolescents Longitudinal Research Team study (MyHeARTs) involving Malaysian adolescents of similar age [34]. The participants for our study were recruited from the urban city of Kuala Lumpur whilst those of MyHeARTs were recruited from three different states in Malaysia, including both urban and rural areas [34]. This study may be the first to examine insulin resistance in an adolescent population in Malaysia. Although the previously reported prevalence of metabolic syndrome was rather low at 3% in our adolescents [35], our findings on insulin resistance are of public health concern. More studies on insulin resistance in adolescents of this country are required to examine possible contributors, both environmental and genetic factors, as it can lead to increasing number of younger patients with type 2 diabetes mellitus. The proportion of adolescents with vitamin D deficiency, were much lower than that found in the MyHeARTs (78%) [14]. The discrepancy may be due to geographical differences as mentioned above. In addition, different assays were used to measure 25(OH)D levels and thus may contribute to the variability. Higher occurrence of vitamin D deficiency in our adolescent girls as compared to boys was similar to findings in primary-school children (7–12 years old) [36] and those of MyHeARTs [14]. Nevertheless, vitamin D deficiency in this population was still comparatively higher than other populations of children and adolescents from different regions in the world; for example, Denmark at 8% [4], Mexico at 18% [6] and Turkey at 25% during summer [37]. Possible contributors to vitamin D deficiency in the Malaysian population have been discussed in great depth by Sadat et al [14].

We demonstrated that those with AA (mutant) genotype of the BsmI had the lowest levels of vitamin D compared to other genotypes and had higher risk of vitamin D deficiency compared to the GG genotype (wild). The A allele itself was significantly associated with vitamin D deficiency in our population. Our finding is consistent with studies on Brazilian children where BsmI polymorphism was found to be associated with serum vitamin D level [26, 28] although it was the wild variant (GG) that was associated with lower vitamin D levels in Brazilian girls [26]. In contrast, the association between BsmI polymorphism and vitamin D levels was not supported in other populations such as the Arabs (Saudi) [38], Uygur and Kazakhs [39] as well as the European adolescents [40], suggesting that this SNP may influence vitamin D levels rather differently in different population.

The allele of the BsmI polymorphism is located at the intronic region of the VDR gene (intron 8) and near the 3’ end of the VDR gene and is thought to possibly alter the expression of other genes’ by affecting the stability of VDR mRNA and its gene transcription [41]. Interestingly, Ogunkolade et al demonstrated that although VDR polymorphisms were significant determinant of VDR mRNa and VDR protein levels in peripheral blood, there was no correlation found between VDR polymorphisms and vitamin D levels. Possible contribution of other VDR genes as well as interaction with environmental factors such as dietary supplements and physical activities needs to be further investigated to fully understand the role of this SNP in vitamin D deficiency in our population. It is worth noting that besides vitamin D levels, this SNP was also associated with calcium absorption especially in younger women [42]. However, a recent meta-analysis did not demonstrate the risk of osteoporosis with BsmI in postmenopausal women [43]. Hence, the role of this SNP in modulating vitamin D levels and its consequences on bone health need to be further elucidated in our population with a longitudinal study design.

As there were no significant differences found when comparing fasting insulin and HOMA-IR across the genotypes after adjustment of covariates, it is possible that other factors such as gender (as shown in the results), ethnicity [44], pubertal stages [32], and socioeconomic status [45] to also influence levels of insulin in this population. However, when insulin resistance was examined as a categorical variable with the cut-off for HOMA-IR chosen based on pubertal stage [32], it was found that the AA genotype of BsmI conferred higher risk of insulin resistance as compared to the GG genotype. The cohort of Brazilian children also reported to demonstrate association between this SNP and insulin resistance [28]. Thus, our study highlights the possible role of VDR BsmI polymorphism in mediating insulin resistance in early adolescence in non-Caucasian populations. In adults, similar findings on the relationship between BsmI and insulin resistance were reported in Caucasians [46] and Bangladeshi Asians [24] while studies on other population such as Brazilians [47], Egyptians and Polish women [48, 49] did not demonstrate significant association. The pancreatic beta cell’s expression of VDR [19] supports the idea that polymorphisms in the VDR region such as the BsmI may exert genomic actions possibly influencing insulin secretion [19, 50]. Ogunkalade et al demonstrated that VDR expression did determine insulin secretory capacity in a functional study [51]. While direct action on the pancreas has been postulated, this SNP has not been demonstrated to be associated with an increased susceptibility to type 1 diabetes mellitus on its own [52]. Although studies on BsmI suggested an association with type 2 diabetes mellitus, a consequence of insulin resistance, this was not reflected by meta-analyses as many of these studies were of small sample sizes [53, 54]. As insulin resistance and type 2 diabetes often arises as a result of a complex interplay between genetic and lifestyle influences, larger population based studies need to be carefully designed to further evaluate the association between this SNP and early onset insulin resistance.

We found a significant association between the A allele of BsmI and risk of insulin resistance in those with concurrent vitamin D deficiency. This is interesting as vitamin D deficiency on its own was not found to be associated with increased risk of insulin resistance. The VDR gene may modify the relationship between vitamin D and insulin resistance [51]. Jain et al previously reported that in New Zealand, South Asian women responded differently in terms of insulin sensitivity when given vitamin D supplementation according to VDR gene polymorphism (FokI) [55]. It was suggested that genotyping of the VDR gene may predict response to vitamin D intervention to improve insulin sensitivity [55]. Thus, future clinical study incorporating the VDR polymorphisms and vitamin D intervention may be desirable in our population of young adolescents with insulin resistance.

We did not find any significant differences between the genotypes of BsmI in terms of adiposity parameters. Furthermore, no significant association was observed between this SNP and risk of obesity, even when stratified according to vitamin D status. This was supported from findings in other populations such as Polish females (post-menopausal) [49], non-Hispanic white US adults [56] and US females [57]. On the contrary, BsmI was significantly associated with obesity in other populations such as the Arabs [38], French Caucasians [58], Polish males [24] and Swedish females [59]. VDR was shown to play a critical role in mediating the inhibitory actions of vitamin D on adipogenesis [60]. The expression of VDR mRNA also changes during adipocyte differentiation [61]. It was hypothesized by Ochs-Balcom et al that linkage disequilibrium with other functional SNPs in the 3’ VDR region may explain the association of BsmI with obesity [57].

This study provides additional data on vitamin D deficiency, obesity and insulin resistance among an adolescent population in Southeast Asia, a cohort on which data is still currently lacking. We implemented direct measurements of anthropometric profiles by trained researchers, therefore reducing self-reporting bias. Even though the schools were selected randomly across different zones in Kuala Lumpur, Malaysia, girls may have been overrepresented possibly due to volunteer bias. Although Malaysia is made up of three major ethnicities; the Malays, Chinese and Indians, only a third of participants were Chinese or Indian, thus inter-ethnic genetic comparisons were not examined. We measured total 25(OH)D which may not truly reflect vitamin D status. It was demonstrated that although total 25(OH)D was different between Americans of African and European decent, the free serum 25(OH)D was comparable [11, 62]. Thus, vitamin D deficiency diagnosed solely on total 25(OH)D measurements may be misleading in certain populations. Hence, more studies are required in our population to establish our vitamin D status.

This study focused solely on the BsmI polymorphism. There are other well studied VDR gene polymorphisms such as rs731236 (TaqI), rs7975232 (ApaI) and rs10735810 (FokI) [18] which may influence vitamin D levels, adiposity and insulin resistance in this population. The FokI polymorphism has been associated with increased risk of type 2 diabetes in the Asian population [54] and hence the role of this SNP needs to be explored further. In addition, the interactions between the genotypes of BsmI and vitamin D deficiency with risk of obesity and insulin resistance were not able to be examined as the sample size was relatively small for this purpose to allow meaningful analysis.

In conclusion, the VDR BsmI polymorphism was significantly associated with vitamin D deficiency and insulin resistance, but not with obesity. Functional studies are needed in future to further characterize the contribution of this polymorphism to risk of vitamin D deficiency and metabolic dysregulation in the younger population.

Acknowledgments

The study was supported by the High Impact Research, Ministry of Higher Education Grant, Malaysia (HIR H00021-00-E000082) and the University Malaya Research Grant (RG310-11HTM). The authors acknowledge the contributions of other co-investigators, postgraduate students and research assistants involved in this community based study including Associate Professor Dr Mohd Nahar Azmi Mohamed, Associate Professor Dr Ivy Chung, Mrs Fatin Fauzi, Mr Amir Hakim, Miss Debbie Ann Loh, Dr Jerri Ling, Dr Sandeep Singh, Miss Siti Aliwiyah Wan Jusoh and Mrs Nor Fazida Rasidi, The authors are grateful to all the students and their parents or guardians for their participation in making this study a success. We would also like to extend our gratitude to all the school teachers and colleagues for their support and assistance in this study. This manuscript has been proofread by the Research Management Centre, Faculty of Medicine, University of Malaya.

Author Contributions

  1. Conceptualization: RR NLZ ZM FMM MYJ.
  2. Data curation: RR NLZ.
  3. Formal analysis: RR NLZ FMM MYJ.
  4. Funding acquisition: NLZ ZM.
  5. Investigation: RR NLZ.
  6. Methodology: RR NLZ ZM FMM MYJ.
  7. Project administration: RR NLZ ZM MYJ.
  8. Resources: RR NLZ ZM MYJ.
  9. Supervision: NLZ ZM MYJ.
  10. Validation: RR NLZ.
  11. Visualization: RR NLZ.
  12. Writing – original draft: RR NLZ FMM MYJ.
  13. Writing – review & editing: RR NLZ FMM ZM MYJ.

References

  1. 1. Santoro D, Caccamo D, Lucisano S, Buemi M, Sebekova K, Teta D, et al. Interplay of Vitamin D, Erythropoiesis, and the Renin-Angiotensin System. BioMed Res Int. 2015;2015: 145828. pmid:26000281
  2. 2. Holick MF. Vitamin D Deficiency. N Engl J Med. 2007;357(3):266–81. pmid:17634462
  3. 3. Erol M, Yigit O, Kucuk SH, Gayret OB. Vitamin D Deficiency in Children and Adolescents in Bagcilar, Istanbul. J Clin Res Pediatr Endocrinol. 2015;7(2):134–9. pmid:26316436
  4. 4. Madsen KH, Rasmussen LB, Mejborn H, Andersen EW, Molgaard C, Nissen J, et al. Vitamin D status and its determinants in children and adults among families in late summer in Denmark. Br J Nutr. 2014;112(5):776–84. pmid:24932732
  5. 5. Djennane M, Lebbah S, Roux C, Djoudi H, Cavalier E, Souberbielle JC. Vitamin D status of schoolchildren in Northern Algeria, seasonal variations and determinants of vitamin D deficiency. Osteoporos Int. 2014;25(5):1493–502. pmid:24566583
  6. 6. Flores M, Macias N, Lozada A, Sánchez LM, Díaz E, Barquera S. Serum 25-hydroxyvitamin D levels among Mexican children ages 2 y to 12 y: A national survey. Nutrition. 2013;29(5):802–4. pmid:23422537
  7. 7. Vierucci F, Del Pistoia M, Fanos M, Gori M, Carlone G, Erba P, et al. Vitamin D status and predictors of hypovitaminosis D in Italian children and adolescents: a cross-sectional study. Eur J Pediatr. 2013;172(12):1607–17. pmid:23959324
  8. 8. Mangin M, Sinha R, Fincher K. Inflammation and vitamin D: the infection connection. Inflamm Res. 2014;63(10):803–19. pmid:25048990
  9. 9. Wise LA, Ruiz-Narváez EA, Haddad SA, Rosenberg L, Palmer JR. Polymorphisms in vitamin D–related genes and risk of uterine leiomyomata. Fertil Steril. 2014;102(2):503–10.e1. pmid:24890271
  10. 10. Misra M PD, Petryk A, Collett-Solberg PF, Kappy M. Vitamin D Defficiency in Children and Its Management: Review of Current Knowledge and Recommendations. Pediatrics. 2008;122:398–417. pmid:18676559
  11. 11. Powe CE, Evans MK, Wenger J, Zonderman AB, Berg AH, Nalls M, et al. Vitamin D–Binding Protein and Vitamin D Status of Black Americans and White Americans. N Engl J Med. 2013;369(21):1991–2000. pmid:24256378
  12. 12. Kuhn T, Kaaks R, Teucher B, Hirche F, Dierkes J, Weikert C, et al. Dietary, lifestyle, and genetic determinants of vitamin D status: a cross-sectional analysis from the European Prospective Investigation into Cancer and Nutrition (EPIC)-Germany study. Eur J Nutr. 2014;53(3):731–41. pmid:24005870
  13. 13. Shafinaz I, Moy F. Vitamin D level and its association with adiposity among multi-ethnic adults in Kuala Lumpur, Malaysia: a cross sectional study. BMC Public Health. 2016 16::232. pmid:26951992
  14. 14. Al-Sadat N, Majid H, Sim P, Su T, Dahlui M, Abu Bakar M, et al. Vitamin D deficiency in Malaysian adolescents aged 13 years: findings from the Malaysian Health and Adolescents Longitudinal Research Team study (MyHeARTs). BMJ Open. 2016 6(8):e010689. pmid:27540095
  15. 15. Ng M, Fleming T, Robinson M, Thomson B, Graetz N, Margono C, et al. Global, regional, and national prevalence of overweight and obesity in children and adults during 1980–2013: a systematic analysis for the Global Burden of Disease Study 2013. Lancet. 2014;384(9945):766–81. pmid:24880830
  16. 16. El-Hajj Fuleihan G, Bouillon R, Clarke B, Chakhtoura M, Cooper C, McClung M, et al. Serum 25-Hydroxyvitamin D Levels: Variability, Knowledge Gaps, and the Concept of a Desirable Range. J Bone Miner Res. 2015;30(7):1119–33. pmid:25952470
  17. 17. Karohl C, Su SY, Kumari M, Tangpricha V, Veledar E, Vaccarino V, et al. Heritability and seasonal variability of vitamin D concentrations in male twins. Am J Clin Nutr. 2010;92(6):1393–8. pmid:20943799
  18. 18. Pike JW, Meyer MB, Benkusky NA, Lee SM, St John H, Carlson A, et al. Genomic Determinants of Vitamin D-Regulated Gene Expression. Vitam Horm. 2016;100:21–44. pmid:26827947
  19. 19. Lee S, Clark SA, Gill RK, Christakos S. 1,25-Dihydroxyvitamin D3 and pancreatic beta-cell function: vitamin D receptors, gene expression, and insulin secretion. Endocrinology. 1994;134(4):1602–10. pmid:8137721
  20. 20. Mutt SJ, Hypponen E, Saarnio J, Jarvelin MR, Herzig KH. Vitamin D and adipose tissue-more than storage. Front Physiol. 2014;5:228. pmid:25009502
  21. 21. Morrison N, Yeoman R, Kelly P, Eisman J. Contribution of trans-acting factor alleles to normal physiological variability: vitamin D receptor gene polymorphism and circulating osteocalcin. Proc Natl Acad Sci USA. 1992;89(15):6665–9 pmid:1353882
  22. 22. Ingles S, Haile R, Henderson B, Kolonel L, Nakaichi G, Shi C, et al. Strength of linkage disequilibrium between two vitamin D receptor markers in five ethnic groups: implications for association studies. Cancer Epidemiol Biomarkers Prev. 1997 6(2):93–8. pmid:9037559
  23. 23. Al-Daghri NM, Al-Attas OS, Alkharfy KM, Khan N, Mohammed AK, Vinodson B, et al. Association of VDR-gene variants with factors related to the metabolic syndrome, type 2 diabetes and vitamin D deficiency. Gene. 2014;542(2):129–33. pmid:24680778
  24. 24. Filus A, Trzmiel A, Kuliczkowska-Plaksej J, Tworowska U, Jedrzejuk D, Milewicz A, et al. Relationship between vitamin D receptor BsmI and FokI polymorphisms and anthropometric and biochemical parameters describing metabolic syndrome. Aging Male. 2008;11(3):134–9. pmid:18821289
  25. 25. Zhao Y, Liao S, He J, Jin Y, Fu H, Chen X, et al. Association of vitamin D receptor gene polymorphisms with metabolic syndrome: a case-control design of population-based cross-sectional study in North China. Lipids Health Dis. 2014;13:129. pmid:25106919
  26. 26. Santos BR, Mascarenhas LPG, Satler F, Boguszewski MCS, Spritzer PM. Vitamin D deficiency in girls from South Brazil: a cross-sectional study on prevalence and association with vitamin D receptor gene variants. BMC Pediatr. 2012;12.
  27. 27. Ferrarezi DAF, Bellili-Munoz N, Nicolau C, Cheurfa N, Guazzelli IC, Frazzatto E, et al. Allelic variations in the vitamin D receptor gene, insulin secretion and parents' heights are independently associated with height in obese children and adolescents. Metabolism. 2012;61(10):1413–21. pmid:22551951
  28. 28. Cobayashi F, Lourenço B, Cardoso M. 25-Hydroxyvitamin D3 Levels, BsmI Polymorphism and Insulin Resistance in Brazilian Amazonian Children. Int J Mol Sci. 2015;16(6):12531. pmid:26047339
  29. 29. World Health Organization. Waist Circumference and Waist–Hip Ratio: Report of a WHO Expert Consultation. Geneva: World Health Organization; 2008
    • 30. Sun Y, Tao F, Su P, China Puberty Research Collaboration. Self-assessment of pubertal Tanner stage by realistic colour images in representative Chinese obese and non-obese children and adolescents. Acta Paediatr. 2012 101(4):e163–6. pmid:22176343
    • 31. World Health Organization. Growth reference for 5–19 years 2007 [updated 2015; cited 2015. Available from: http://www.who.int/growthref/en/.
      • 32. Yin J, Li M, Xu L, Wang Y, Cheng H, Zhao X, et al. Insulin resistance determined by Homeostasis Model Assessment (HOMA) and associations with metabolic syndrome among Chinese children and teenagers. Diabetol Metab Syndr. 2013;5(1):71. pmid:24228769
      • 33. Sanders T, Feng X, Fahey PP, Lonsdale C, Astell-Burt T. Greener neighbourhoods, slimmer children? Evidence from 4423 participants aged 6 to 13 years in the Longitudinal Study of Australian children. Int J Obes. 2015;39(8):1224–9.
      • 34. Hazreen MA, Su TT, Jalaludin MY, Dahlui M, Chinna K, Ismail M, et al. An exploratory study on risk factors for chronic non-communicable diseases among adolescents in Malaysia: overview of the Malaysian Health and Adolescents Longitudinal Research Team study (The MyHeART study). BMC Public Health. 2014;14 Suppl 3:S6.
      • 35. Fadzlina AA, Harun F, Nurul Haniza MY, Al Sadat N, Murray L, Cantwell MM, et al. Metabolic syndrome among 13 year old adolescents: prevalence and risk factors. BMC Public Health. 2014;14 Suppl 3:S7.
      • 36. Khor GL, Chee WS, Shariff ZM, Poh BK, Arumugam M, Rahman JA, et al. High prevalence of vitamin D insufficiency and its association with BMI-for-age among primary school children in Kuala Lumpur, Malaysia. BMC Public Health. 2011;11:95. pmid:21310082
      • 37. Andýran N, Çelik N, Akça H, Doðan G. Vitamin D Deficiency in Children and Adolescents. J Clin Res Pediatr Endocrinol. 2012;4(1):25–9. pmid:22394709
      • 38. Al-Daghri NM, Al-Attas OS, Alkharfy KM, Khan N, Mohammed AK, Vinodson B, et al. Association of VDR-gene variants with factors related to the metabolic syndrome, type 2 diabetes and vitamin D deficiency. Gene. 2014;542(2):129–33. pmid:24680778
      • 39. Xu XJ, Mao JF, Zhang MC, Liu HM, Li HX, Lei H, et al. Vitamin D Deficiency in Uygurs and Kazaks Is Associated with Polymorphisms in CYP2R1 and DHCR7/NADSYN1 Genes. Med Sci Monit. 2015;21:1960–8. pmid:26149120
      • 40. Valtuena J, Gonzalez-Gross M, Huybrechts I, Breidenassel C, Ferrari M, Mouratidou T, et al. Factors associated with vitamin D deficiency in European adolescents: the HELENA study. J Nutr Sci Vitaminol (Tokyo). 2013;59(3):161–71.
      • 41. Speer G, Cseh K, Winkler G, Vargha P, Braun E, Takacs I, et al. Vitamin D and estrogen receptor gene polymorphisms in type 2 diabetes mellitus and in android type obesity. Eur J Endocrinol. 2001;144(4):385–9. pmid:11275948
      • 42. Kiel DP, Myers RH, Cupples LA, Kong XF, Zhu XH, Ordovas J, et al. The BsmI vitamin D receptor restriction fragment length polymorphism (bb) influences the effect of calcium intake on bone mineral density. J Bone Miner Res. 1997;12(7):1049–57. pmid:9200004
      • 43. Jia F, Sun RF, Li QH, Wang DX, Zhao F, Li JM, et al. Vitamin D receptor BsmI polymorphism and osteoporosis risk: a meta-analysis from 26 studies. Genet Test Mol Biom arkers. 2013;17(1):30–4.
      • 44. Kodama K, Tojjar D, Yamada S, Toda K, Patel CJ, Butte AJ. Ethnic differences in the relationship between insulin sensitivity and insulin response: a systematic review and meta-analysis. Diabetes Care. 2013;36(6):1789–96. pmid:23704681
      • 45. Pulkki L, Keltikangas-Jarvinen L, Ravaja N, Viikari J. Child-rearing attitudes and cardiovascular risk among children: moderating influence of parental socioeconomic status. Prev Med. 2003;36(1):55–63. pmid:12473425
      • 46. Oh JY, Barrett-Connor E. Association between vitamin D receptor polymorphism and type 2 diabetes or metabolic syndrome in community-dwelling older adults: the Rancho Bernardo Study. Metabolism. 2002;51(3):356–9. pmid:11887173
      • 47. Schuch NJ, Garcia VC, Vivolo SR, Martini LA. Relationship between Vitamin D Receptor gene polymorphisms and the components of metabolic syndrome. Nutr J. 2013;12:96. pmid:23855914
      • 48. Mackawy AM, Badawi ME. Association of vitamin D and vitamin D receptor gene polymorphisms with chronic inflammation, insulin resistance and metabolic syndrome components in type 2 diabetic Egyptian patients. Meta gene. 2014;2:540–56. pmid:25606437
      • 49. Tworowska-Bardzinska U, Lwow F, Kubicka E, Laczmanski L, Jedzrzejuk D, Dunajska K, et al. The vitamin D receptor gene BsmI polymorphism is not associated with anthropometric and biochemical parameters describing metabolic syndrome in postmenopausal women. Gynecol Endocrinol. 2008;24(9):514–8. pmid:18958772
      • 50. Christakos S, Raval-Pandya M, Wernyj RP, Yang W. Genomic mechanisms involved in the pleiotropic actions of 1,25-dihydroxyvitamin D3. Biochemical J. 1996;316 (Pt 2):361–71.
      • 51. Ogunkolade BW, Boucher BJ, Prahl JM, Bustin SA, Burrin JM, Noonan K, et al. Vitamin D receptor (VDR) mRNA and VDR protein levels in relation to vitamin D status, insulin secretory capacity, and VDR genotype in Bangladeshi Asians. Diabetes. 2002;51(7):2294–300. pmid:12086963
      • 52. Tizaoui K, Kaabachi W, Hamzaoui A, Hamzaoui K. Contribution of VDR polymorphisms to type 1 diabetes susceptibility: Systematic review of case-control studies and meta-analysis. J Steroid Biochem Mol Biol. 2014;143:240–9. pmid:24742873
      • 53. Yu F, Cui LL, Li X, Wang CJ, Ba Y, Wang L, et al. The genetic polymorphisms in vitamin D receptor and the risk of type 2 diabetes mellitus: an updated meta-analysis. Asia Pac J Clin Nutr. 2016;25(3):614–24. pmid:27440697
      • 54. 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–41. pmid:23506721
      • 55. Jain R, von Hurst PR, Stonehouse W, Love DR, Higgins CM, Coad J. Association of vitamin D receptor gene polymorphisms with insulin resistance and response to vitamin D. Metabolism. 2012;61(3):293–301. pmid:21871642
      • 56. Beydoun MA, Tanaka T, Beydoun HA, Ding EL, Ferrucci L, Zonderman AB. Vitamin D receptor and megalin gene polymorphisms are associated with central adiposity status and changes among US adults. J Nutr Sci. 2013;2:e33. pmid:25191583
      • 57. Ochs-Balcom HM, Chennamaneni R, Millen AE, Shields PG, Marian C, Trevisan M, et al. Vitamin D receptor gene polymorphisms are associated with adiposity phenotypes. Am J Clin Nutr. 2011;93(1):5–10. pmid:21048058
      • 58. Ye WZ, Reis AF, Dubois-Laforgue D, Bellanne-Chantelot C, Timsit J, Velho G. Vitamin D receptor gene polymorphisms are associated with obesity in type 2 diabetic subjects with early age of onset. Eur J Endocrinol. 2001;145(2):181–6. pmid:11454514
      • 59. Grundberg E, Brandstrom H, Ribom EL, Ljunggren O, Mallmin H, Kindmark A. Genetic variation in the human vitamin D receptor is associated with muscle strength, fat mass and body weight in Swedish women. Eur J Endocrinol. 2004;150(3):323–8. pmid:15012617
      • 60. Kong J, Li YC. Molecular mechanism of 1,25-dihydroxyvitamin D3 inhibition of adipogenesis in 3T3-L1 cells. Am J Physiol Endocrinol Metab. 2006;290(5):E916–24. pmid:16368784
      • 61. Burton GR, Guan Y, Nagarajan R, McGehee RE Jr. Microarray analysis of gene expression during early adipocyte differentiation. Gene. 2002;293(1–2):21–31. pmid:12137940
      • 62. Shieh A, Aloia JF. Assessing Vitamin D Status in African Americans and the Influence of Vitamin D on Skeletal Health Parameters. Endocrinol Metab Clin North Am. 2017;46(1):135–52. pmid:28131129