Next Article in Journal
The Prevalence of Dietary Supplement Use Among Elementary, Junior High, and High School Students: A Nationwide Survey in Japan
Next Article in Special Issue
High Fat Mixed Meal Tolerance Test Leads to Suppression of Osteocalcin Decrease in Obese Insulin Resistant Subjects Compared to Healthy Adults
Previous Article in Journal
Nontoxic Glucomoringin-Isothiocyanate (GMG-ITC) Rich Soluble Extract Induces Apoptosis and Inhibits Proliferation of Human Prostate Adenocarcinoma Cells (PC-3)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Nutrients
Volume 10
Issue 9
10.3390/nu10091175
nutrients-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Association between Vitamin D Deficiency and Single Nucleotide Polymorphisms in the Vitamin D Receptor and GC Genes and Analysis of Their Distribution in Mexican Postmenopausal Women

by
Berenice Rivera-Paredez
1,
Nayeli Macías
2,
Mayeli M. Martínez-Aguilar
3,
Alberto Hidalgo-Bravo
4,
Mario Flores
2,
Amado D. Quezada-Sánchez
5,
Edgar Denova-Gutiérrez
2,
Miguel Cid
6,
Angelica Martínez-Hernández
6,
Lorena Orozco
6,
Manuel Quiterio
7,
Yvonne N. Flores
8,9,
Jorge Salmerón
1,7 and
Rafael Velázquez-Cruz
3,*
1
Academic Unit in Epidemiological Research, Research Center in Policies, Population and Health, School of Medicine, National Autonomous University of Mexico, Mexico City 04510, Mexico
2
Nutrition and Health Research Center, National Institute of Public Health (INSP), Cuernavaca, Morelos 62100, Mexico
3
Genomics of Bone Metabolism Laboratory, National Institute of Genomic Medicine (INMEGEN), Mexico City 14610, Mexico
4
Department of Genetics, National Institute of Rehabilitation (INR), Mexico City 14389, Mexico
5
Center for Evaluation and Surveys Research, National Institute of Public Health (INSP), Cuernavaca, Morelos 62100, Mexico
6
Inmunogenomics and Metabolic Diseases Laboratory, National Institute of Genomic Medicine (INMEGEN), Mexico City 14610, Mexico
7
Center for Population Health Research, National Institute of Public Health (INSP), Cuernavaca, Morelos 62100, Mexico
8
Epidemiology and Health Services Research Unit, Mexican Institute of Social Security, Cuernavaca, Morelos 62000, Mexico
9
UCLA Department of Health Policy and Management, UCLA Kaiser Permanente Center for Health Equity, Fielding School of Public Health and Jonsson Comprehensive Cancer Center, Los Angeles, CA 90001, USA
*
Author to whom correspondence should be addressed.
Nutrients 2018, 10(9), 1175; https://doi.org/10.3390/nu10091175
Submission received: 24 July 2018 / Revised: 18 August 2018 / Accepted: 23 August 2018 / Published: 27 August 2018
(This article belongs to the Special Issue Nutritional Status and Bone Health)
Browse Figures
Versions Notes

Abstract

:
Genome-wide association studies in people with European ancestry suggest that polymorphisms in genes involved in vitamin D (VD) metabolism have an effect on serum concentrations of 25-hydroxyvitamin D. However, nothing is known about these polymorphisms in populations with Amerindian ancestry. Our aim was to evaluate the association between genetic variants on the vitamin D receptor (VDR) and the vitamin D binding protein (GC) genes, involved in the VD pathway, and VD deficiency in 689 unrelated Mexican postmenopausal women. We also described the frequencies of these variants in 355 postmenopausal women from different ethnic groups. Based on our preliminary results of 400 unrelated Mexican postmenopausal women, three single nucleotide polymorphisms (SNPs) were selected for genotyping. The SNPs rs4516035 in VDR and rs2282679 in GC were associated with VD deficiency. Additionally, women who carried three risk alleles had a 3.67 times higher risk of suffering VD deficiency, compared to women with no risk alleles (p = 0.002). The rs4516035-C allele frequency in the Amerindian population was enriched in the South East region of Mexico. In contrast, the highest frequency of the rs2298850-C allele, a proxy for the tag SNP rs2282679, was observed in the South region. Our results indicate that genetic variants in VDR and GC genes are associated with VD deficiency in Mexican postmenopausal women. Moreover, an association was observed for the variants rs3794060 and rs4944957 of the DHCR7/NADSYN1 gene with osteopenia/osteoporosis.

Graphical Abstract

1. Introduction

Mexico is undergoing demographic changes characterized by a growing number of people over the age of 50 years. The 2015 population census reported an estimated 23 million people within this age group. According to The National Population Council of Mexico (CONAPO), the proportion of individuals over 50 years of age will increase more than 55% by 2020, and up to 200% by 2050 [1]. The aging of the population, accompanied by an extended life expectancy, may result in a higher prevalence of osteopenia and osteoporosis. In Mexico, the estimated prevalence of osteopenia and osteoporosis is 56% and 16%, respectively [2]. Furthermore, osteoporosis increases fragility fractures, resulting in higher morbidity, mortality and economic burden, which is a growing public health problem [2]. Osteoporosis also causes approximately 30,000 femur and/or hip fractures per year and accounts for 19.5% of vertebral fractures in women [3].
Vitamin D (VD) regulates calcium (Ca) and the intestinal absorption of phosphorous. As part of the aging process, serum parathyroid hormone increases and serum 25-hydroxyvitamin D (25(OH) D) diminishes [4,5]. Low serum concentrations of 25(OH)D are associated with low bone mineral density (BMD) [4,5]. A reduction of BMD can lead to osteopenia and osteoporosis; therefore, detection of 25(OH)D deficiency should be considered a public health priority. Globally, an estimated 1 billion people have some form of VD deficiency [6]. In Mexico, the prevalence of VD deficiency ranges from 8% to 24% in all age groups. However, among post-menopausal women, the prevalence has been reported to be as high as 46.7% [7,8,9]. The main source of 25(OH)D is 7-dehydrocholesterol which is activated by sunlight ultraviolet B (UVB) radiation. In order to become its active molecule, 7-dehydrocholesterol goes through hydroxylation, first in the liver and then in the kidneys [10]. The main food sources of VD are milk, oily fish and fortified foods. Some factors that can affect VD concentrations are age, body adiposity, race, and genetic variants [11]. Genetic factors have been found to contribute to VD concentration by as much as 53% [12]. The genetic variants derived from genome-wide association studies (GWAS) conducted on European populations that are associated with VD concentrations are engaged in cholesterol synthesis (NADSYN1/DHCR7 gene), hydroxylation (CYP2R1 gene) and VD transport (GC gene) [13,14]. To the best of our knowledge, there are no studies that have evaluated the relationship between VD metabolism and its deficiency with the genetic variants of the above-mentioned genes in Mexican population. The aim of the present study was to evaluate the association between genetic variants in the genes involved in the VD pathway and VD deficiency in Mexican postmenopausal women. In addition, the frequencies of these variants are described in Mexican Amerindian ethnic groups.

2. Materials and Methods

2.1. Participants—Second Wave of the Health Workers Cohort Study

The present study is a cross-sectional analysis. The data was collected during a second assessment of the “The Health Workers Cohort Study” (HWCS) from 2010–2012. The HWCS is a prospective study conducted in workers from the Mexican Institute of Social Security (IMSS) in Cuernavaca, Morelos (located in the central region of Mexico). The HWCS focuses on the association between different health endpoints with genetic background and lifestyle. A total of 1855 HWCS participants aged 18 to 85 years were considered for the second evaluation. The variables and procedures of the HWCS have been described in detail somewhere else [15]. We analyzed the SNPs and VD concentrations of 689 out of 754 postmenopausal women evaluated during 2010–2012 (Figure S1). We excluded a total of 65 women; 35 related women, 27 without SNP quality control and three women with missing 25(OH)D concentrations. In order to know the frequencies of the genetic variants of interest in Mexican Amerindian populations, we evaluated an independent sample of 355 postmenopausal women from 37 different ethnic groups taken from the Metabolic Analysis in an Indigenous Sample (MAIS) cohort [16]. Women from the MAIS cohort were included according to the following criteria: (1) women who identified themselves as indigenous, (2) her parents and grandparents spoke the same language and (3) they were born in the same area as their parents and grandparents. The Mexican population is highly heterogenic and according to previous data, there is a gradient of European ancestry from North to South [16,17,18]. Based on this, the Amerindian populations were classified into five geographic regions: North (N), Central East (CE), Central West (CW), South (S) and Southeast (SE) [16]. This study was planned and performed according to the guidelines of the Declaration of Helsinki. The Research Ethics Committee of the IMSS (No. 12CEI 09 006 14) approved the study protocol and informed consent forms. Written informed consent was obtained from all participants.

2.2. Data Collection and Blood Sample Collection

Participants answered a self-administered questionnaire focused on characteristics such as birth date, education, marital status, medical family history, past medical history, current medication use, lifestyle information (e.g., diet, physical activity, smoking status, alcohol consumption). Anthropometric body composition measurements and blood samples were collected using standard procedures. Physical activity information was obtained with a questionnaire validated in Spanish in a sample with similar characteristics [19,20]. Participants reported recreational physical activity such as walking, running, etc. during a typical week of the previous year.

2.3. Anthropometric Measurements

Body weight was measured with a calibrated electronic scale (model BC-533, Tanita, IL, USA) with participants wearing underwear and bare feet. Height was taken using a conventional stadiometer (Seca HH, Hamburg, Germany), participants were standing barefoot with relaxed shoulders. Height was taken while touching the top of the head at the moment of maximum inspiration. Waist circumference was measured at the highest point of the iliac crest [21]. All measurements had concordance coefficients between 0.83 and 0.90. All measurements were performed by trained personnel using standardized techniques [21].

2.4. Body Composition Measurements

Bone mineral density (BMD), bone mineral content, lean body mass, appendicular lean body mass, and fat mass, were obtained with a Lunar DPX NT instrument (Lunar Radiation Corp., Madison, WI, USA) with standardized procedures and trained personnel [22]. Briefly, daily quality assurance scans were conducted using the phantom provided by the manufacturer; the daily variation coefficient was within normal operational standards and the in vivo variation coefficient was lower than 1.0%–1.5% [15]. Women attended the body composition evaluation after 12 h of fasting, neither drinking alcohol nor exercising. A dual-energy X-ray absorptiometry (DEXA) measurement was performed at the lumbar spine (L2–L4), hip, femoral neck, Ward’s angle and total body, to obtain density values expressed as g/cm2. We used World Health Organization criteria, T-score, to classify women with osteopenia or osteoporosis. A T-score > −1 standard deviation (SD) was considered as normal, from −1.0 to −2.5 SD as osteopenia, and under −2.5 SD as osteoporosis [23].

2.5. Dietary Measurements

We used a semi-quantitative food frequency questionnaire, previously validated in a Mexican population [24]. The questionnaire collects data about the frequency of consumption of 116 food items during the previous year. The instrument specifies commonly used size portions. We calculated the energy and nutrient intake by multiplying the frequency of consumption of each food by the nutrient content. We took the information from a comprehensive database of food contents [25].

2.6. Biological Measurements

Blood samples, approximately 20 mL, were taken after 8 h of fasting. We evaluated aminotransferase, triglycerides, glucose, total cholesterol, high-density lipoproteins, low-density lipoproteins, uric acid, creatinine, etc. All biomedical assays were performed using a Selectra XL instrument (Randox Laboratories Ltd., Antrim, UK), in concordance with the proceedings of the International Federation of Clinical Chemistry and Laboratory Medicine [26]. Genomic DNA was isolated from peripheral blood leukocytes using the Gentra Puregene Blood Kit (QIAGEN, Minneapolis, MN, USA), according to the manufacturer’s instructions.

2.7. Measurement of Plasma 25(OH)D

Serum 25-hydroxyvitamin D was measured with an Abbott Architect® Chemiluminescent Microparticle Immuno Assay (CMIA, Lake Bluff, IL, USA). The intra- and inter-assay variation coefficients of the method have been reported as <10%. This method has shown high reactivity (~100%) with 25-hydroxyvitamin D3 and an acceptable performance compared to LC/MS/MS (r = 0.73) [27]. Chemiluminescence has been reported as an accurate and reproducible method [28]. VD deficiency was defined as serum 25(OH)D levels <20 ng/mL, as reported in a previous study in an elderly Mexican population [9].

2.8. Single Nucleotide Polymorphisms Selection and Genotyping

Single nucleotide polymorphisms (SNPs) identified through genome-wide association studies of 25(OH)D serum concentrations and previous association studies in European populations were considered for inclusion. Twenty-nine SNPs involved in the VD metabolic pathway were selected (Table S1) [11,29,30]. In the first stage, a GoldenGate design (Illumina, San Diego, CA, USA) with a total of 384 SNPs (including the 29 SNPs within the VD metabolic pathway) was performed on 400 postmenopausal women from the HWCS. The information obtained was used to test associations between the SNPs and BMD. Information of the parental frequencies and genotype data from a panel of 96 Ancestry Informative Markers (AIMs) was included in the analysis. This information was included in order to control for the effect of false associations due to population stratification. Genotype data of the 29 SNPs was extracted from the initial study for further analysis. The primary analysis of the genotyping data was carried out using the Illumina GenomeStudio software v.2011.1 (Illumina, Inc., San Diego, CA, USA), as described previously [31]. Based on this analysis, the SNPs rs10783219 and rs4516035 on the VDR gene and rs2282679 on the GC gene were selected, considering the statistical significance (p < 0.05) and that they were not in linkage disequilibrium, for further genotyping in the remaining 289 postmenopausal women. The SNPs were genotyped using commercial predesigned TaqMan Probes (Applied Biosystems, Foster City, CA, USA) and called automatically using the SDS 2.2.1 software (Applied Biosystems, Foster City, CA, USA). In both stages, quality control exclusions were implemented because of low SNP call rate (<97%), deviation from the Hardy–Weinberg equilibrium (p ≤ 0.05), related women and gender concordance. After quality control, 689 postmenopausal women with both genetic data and 25(OH)D measurements were available for analysis. The MAIS cohort samples, were previously genotyped with the Genome-Wide Human SNP 6.0 Microarray (Affymetrix, Santa Clara, CA, USA) [16].

2.9. Construction of the Genetic Risk Score (GRS)

To evaluate the combined effect of the SNPs that were significantly associated with vitamin D deficiency, we constructed a genetic risk score (GRS) [32] for each individual, which included two SNPs (rs2282679-GC and rs4516035-VDR). The GRS (ranging from 0 to 4) was constructed by summing the number of risk alleles from these two SNPs for each individual. Genotypes for each SNP were scored using an additive model (0 for homozygous for the non-risk allele, 1 for heterozygous, and 2 for homozygous for the risk allele).

2.10. Statistical Analysis

Samples were grouped according to VD serum concentrations. A t-test and chi-square test were used for analyzing sociodemographic and clinical variables between subjects with deficient and normal concentrations of VD. Genotype and allele frequencies were calculated for every SNP. The Hardy-Weinberg equilibrium was tested for each SNP using a chi-square test. A multiple logistic regression model was fit using vitamin D deficiency as the dependent variable; all the independent variables were simultaneously included. The association between SNPs or GRS and VD deficiency was evaluated through logistic regression models adjusted by co-variables of interest in order to obtain odds ratio (ORs) and 95% confidence intervals (95% CI). The first model was adjusted by age (<60, 60−74, >74 years), body mass index (normal <25 kg/m2, overweight 25 to <30 kg/m2, obesity ≥30 kg/m2), VD intake (tertiles), season of the year when the blood sample was collected (spring, summer, winter, autumn) and physical activity (active ≥30 min/day). A second model using body fat percentage instead of body mass index (BMI) was generated. However, because the estimators did not change, only the second model, adjusted by body fat percent, is presented. A plot of the prevalence of vitamin D deficiency or OR (CI 95%) and GRS was generated. In addition, a sensitivity analysis was performed using a linear regression model to evaluate the association between serum 25-(OH) D levels and the SNPs of the GC and VDR genes; adjusting for co-variables. To approximate normality in vitamin D levels, a square root transformation was performed, as described in Xu et al. [33]. The association between SNPs and BMD was evaluated through logistic regression adjusted by age and BMI. In addition, VD concentrations (non-deficient and deficient), blood collection season (winter, spring, summer and autumn) and physical activity (active ≥30 min/day) were included in a second logistic regression model. A value of p < 0.05 was considered to be statistically significant. All statistical tests were performed with STATA 13.0 (Stata Corp, College Station, TX, US) for Windows [34]. We used the HaploView 4.1 software (Broad Institute, Cambridge, MA, USA) to evaluate the linkage disequilibrium [35]. Allelic frequency distribution in Mexican Amerindian (MA) populations was represented on maps taken from the National Commission of Knowledge and Use of Biodiversity (CONABIO). Maps were analyzed using the software R and QGIS 2.14 (Quantum Geographical Information Systems) [16]. To allow for multiple testing (3 SNPs), we report the results of a Bonferroni-corrected p value of 0.02.

3. Results

3.1. Participant Characteristics

The study included a total of 689 non-related postmenopausal women between 45–92 years. The average age was 62 years (SD ± 8.8), more than 45% were overweight and around 30% were obese. Average serum concentration of VD was 21.0 ng/mL and mean VD intake was 199 UI/day. The prevalence of VD deficiency was near 44%. The most frequent chronic diseases in the study population were hypertension (45%), osteopenia/osteoporosis (38%) and hyperuricemia (33%) (Table 1).
There was a positive association between VD status and age; VD-deficient women were older, had a lower mean BMD of the femoral neck and spine and had a higher prevalence of diabetes in comparison to women without VD deficiency (p < 0.05). Moreover, VD deficient-women had a lower average of VD intake in comparison to women with no deficiency (p < 0.05) (Table 2).

3.2. Factors Asociated with Vitamin D Deficiency

Logistic regression analysis was adjusted for education, smoking, VD intake, season of blood collection, leisure time physical activity and body mass index. Results revealed that women between 60 and 74 years and women older than 74 years old had an OR of 1.95 (95% CI 1.36−2.79) and 3.07 (95% CI 1.70−5.56), respectively of having VD deficiency compared to woman between 45–59 years old. Additionally, it was observed that women with the highest intake of VD had a 39% (OR = 0.61, 95% CI: 0.39−0.95) lower probability of presenting with VD deficiency in comparison to women in the lowest tertile of VD intake. Results were adjusted by the same variables mentioned above (Table S2).

3.3. Association between SNPs and Vitamin D Deficiency

In the first stage, we identified five SNPs lying on the GC gene, which is involved in VD transport. We also identify three SNPs on the vitamin D receptor gene (VDR) that were positively associated to VD deficiency (Table S3). This result was conserved after excluding 29 women with a VD serum concentration beyond 29 ng/mL; for example, the SNP rs4516035 had an OR of 3.50 (95% CI: 1.37−8.89) for VD deficiency (data not shown). The results of crude and adjusted analysis for the VDR variants were statistically significant. We observed an increased risk for VD deficiency with the variants rs10783219, rs4516035 and rs7139166 under co-dominant and additive genetic models (Table S3).
The association for the GC gene variants (rs17467825, rs2282679, rs3755967, rs2298850 and rs1155563), under an additive model, was statistically significant in both the crude and adjusted models. For the SNPs rs17467825, rs2282679 and rs3755967, the associations with VD deficiency remained statistically significant for the risk genotype under a co-dominant model. However, for the rs2298850 and rs1155563 variants, the association did not remain statistically significant under the co-dominant model. The SNPs rs7041 and rs12512631 had an inverse association with VD deficiency (p = 0.03), only with the unadjusted model. The polymorphisms associated with VD deficiency had strong linkage disequilibrium. The SNPs rs17467825, rs2282679, rs3755967, and rs2298850 of the GC gene had an r2 value of >0.93 and for the VDR SNPs, rs4516035 and rs7139166 had an r2 value of 0.95 (Figure S2). Two SNPs were chosen as the tagging SNPs (rs4516035 and rs2282679).
In the second stage, the SNPs rs10783219, rs4516035 and rs2282679 were selected to be evaluated in the entire postmenopausal women cohort. For SNPs rs4516035 and rs2282679, under the additive model, it was observed that the OR increased with each additional copy of the minor allele (OR = 1.40, 95% CI: 1.08−4.43 and OR = 1.53, 95% CI: 1.15−2.04, respectively). However, for the SNP rs10783219, the association was no longer statistically significant (p = 0.377) (Table 3). Additionally, we evaluated the serum 25-(OH)D levels and the SNPs of the GC and VDR genes. However, only the SNP rs2282679 (β = −0.12 (95% CI: −0.20, 0.04), p = 0.004) of the GC gene was found to be associated with serum 25(OH)D levels under the additive model (data not shown).

3.4. Association between GRS and Vitamin D Deficiency

The results of the genetic risk score generated with the combination of two variants showed that the greater the number of risk alleles, the higher the VD deficiency prevalence. For individuals with no risk alleles, the prevalence of VD deficiency was 36.8 compared to 66.1 for individuals with three or more risk alleles (p = 0.002) (Figure 1A). In addition, the genetic risk score (combining the two variants) was associated with the risk of VD deficiency, with an adjusted OR in the top category of 3.34 (95% CI: 1.46–7.68, p = 0.004) (Figure 1B).

3.5. Association between SNPs and Bone Mineral Density

We also explored whether the variants have some relationship with bone mineral density (BMD). A statistically significant association was observed for the variants rs3794060 and rs4944957 of the DHCR7/NADSYN1 gene with osteopenia/osteoporosis. After adjusting for age and body mass index, the association persisted for the two variants (p < 0.05) (Table S4). Additionally, the model was adjusted for VD concentrations, blood collection season and physical activity, and the association was not affected (p < 0.05) (data not shown).

3.6. Minor Allele Frequency

Table S1 shows the minor allele frequency (MAF) distribution between Utah residents with Northern and Eastern European ancestry (CEU), Los Angeles residents with Mexican ancestry (MXL), and the Mexican cohort HWCS. MAFs of the tagging SNPs rs4516035 (on the VDR gene) and rs2282679 (on the GC gene) were compared between CEU, MXL and HWCS populations. The “C” allele of the rs4516035 SNP was more frequent in European populations (0.37), while the frequency in the HWCS cohort was very similar to the frequency observed in the population with Mexican ancestry from Los Angeles (0.267 versus 0.27). With respect to the frequency of the “C” allele of rs2298850, which is a proxy for rs2282679, it was more frequent in residents from Los Angeles with Mexican ancestry (0.27); the frequency in the HWCS cohort was similar to the frequency observed in European populations (0.23 in both) (Table S1).

3.7. Minor Allele Frequency from Different Ethnic Groups

To obtain a global view of the frequency of the rs4516035 and rs2298850 SNPs in the Mexican population, we were able to evaluate the allele and genotype frequencies in 37 different Mexican Amerindian groups (Table 4). The genotype distribution for the two SNPs evaluated was in Hardy–Weinberg equilibrium in Mexican Amerindian populations. Ancestry analysis showed an average Amerindian ancestry of 95 ± 5.7% in the MAIS sample [16]. The 37 Mexican Amerindian (MA) groups were sorted into five major geographic regions (North (N), Central East (CE), Central West (CW), South (S) and South East (SE)). This analysis showed a high heterogeneity of the allele and genotype frequencies for both polymorphisms between regions. Frequencies followed a geographic gradient across Mexico, similar to previous reports [36]. The frequency of the rs4516035-C allele was enriched in the SE region (29%), followed by the S (23%) and CE regions (20%). The lowest frequency was found in the N region (8%) (Figure 2A). The rs4516035-C allele showed the lowest frequency in the N (0–20%) and in a few groups from the other geographic regions, such as Tlahuicas (CE) and Popolucas from the Sierra (CE). Of note, the CC genotype was absent in groups from the N and CW regions, and was only present in a few groups such as Nahuatls from Morelos (CE); Chontals from Oaxaca, Huaves and Mixtecos from the Costa (S); Chujs, Kanjobals, Kaqchikel, Mayas, Mochos and Tseltals (SE) (Table 4 and Figure 2B). Notably, the highest frequencies of the rs4516035-C allele were observed among Nahuatls from Morelos (56%) in the CE region; Zapotecos (50%) and Chontals from Oaxaca (33%) in the S region; Tseltals (50%), Kanjobals (45%), Jakaltecos (30%) and Kaqchikels (28%), inhabiting the SE (Table 4 and Figure 2B). The geographic distribution of the rs2298850-C allele frequency, the proxy for rs2282679 (the tag SNP), followed a gradient opposite to that observed for rs4516035 (Figure 3). In general, frequencies were very similar through geographic areas: N (18%), CW (11%), CE (10%) and SE (16%), with the exception of the S region that showed the highest frequency (21%) (Table 4 and Figure 3A). The CC genotype was only present in a few groups, such as Guarijio (N), Matlaltzinca (CE); Chontal and Mazateco from Oaxaca, and Mixteco from the Costa in the S region and Mam (SE) (Table 4). The highest frequencies of the rs2298850-C allele were found in the Guarijío (56%), Mixteco-Costa (30%), Chontal from Oaxaca (33%) and Mam (33%) (Figure 3B).

4. Discussion

The aim of this study was to evaluate the association between SNPs lying on genes involved in VD metabolism with vitamin D deficiency in a sample of Mexican postmenopausal women. We identify two genetic variants of genes participating in the VD metabolic pathway associated with VD deficiency. Under an additive model, we observed a 40% increased risk of suffering VD deficiency for each risk allele of the rs4516035 SNP. Similarly, the presence of two risk alleles of GC, the rs2282679 SNP, was highly associated with VD deficiency. Furthermore, as the number of alleles increased, there was a higher risk of having VD deficiency.
VD deficiency was highly prevalent in our study at almost forty percent. This prevalence is similar to that reported by the Nutrition and Health National Survey in Mexico performed in 2012 [37]. Other studies have shown similar findings; in a sample of Mexican elderly women, there was a prevalence of VD deficiency of 46.8%, compared with the prevalence of 63.0% in Hispanic women and 29.5% in women with osteoporosis [9,38,39].
It has been documented that age, BMI and VD intake are the main variables associated with VD deficiency [9,40,41]. We observed a higher frequency of obesity in VD deficient women compared to the group with no deficiency, however, this difference was not statistically significant. A previous study in the Mexican elderly population observed that for each unit of BMI increment, the OR for suffering VD deficiency increased 1.03 times (95% CI: 1.01–1.06) [9]. Additionally, a study in older Mexican adults showed that overweight/obese individuals had 1.78 times risk of having VD concentrations within the first tertile (≤20.4 ng/mL) and 1.94 times within second tertile (20.5–26.6 ng/mL), compared to individuals with normal weight [41]. Although a relationship between BMI and VD status has been documented, some studies have not found such association. For example, Pathak et al., did not find an effect of VD supplementation on adiposity in their meta-analysis. Thus, the association between adiposity and VD remains controversial [42]. A possible explanation for not finding an association between BMI and VD could be the genetic heterogeneity of the Mexican population [16,17,18,37].
VD concentrations depend on many factors, including environmental and genetics. Likewise, Forrest et al. [38], showed that Hispanic women were 3.2 more likely to have VD deficiency than Caucasian women. Furthermore, having no college education, being obese, aging and no daily milk consumption are some of the factors also related to VD deficiency in addition to ethnicity. In our study, we found that women between 45 to 59 years and 60 to 74 years had a frequency of VD deficiency of 30% and 50%, respectively. Another factor associated with vitamin D deficiency is its intake. In our study, women with VD deficiency had a lower intake. However, these findings do not agree with the results from Contreras-Manzano et al. [37]. They did not observe differences in average VD consumption between the group with and without VD deficiency. Possible explanations for their findings could be the age of the women evaluated, and the food-frequency questionnaire used. Another variable frequently associated with VD status is the season of the year when evaluations are performed [40]. However, in this sample of Mexican postmenopausal women, there was no relationship between these two variables. This may be explained by the geographic situation of Mexico. The tropic of Cancer crosses above the central region of the country; therefore, changes in sunlight hours between summer and winter are minimal. This effect is especially important in the central and south regions.
An important finding was that women with VD deficiency had higher levels of fasting serum glucose and a higher prevalence of type 2 diabetes. VD has been previously related to insulin resistance and type 2 diabetes in observational studies [9,38]. Carrillo-Vega, et al. [9] observed that increased levels of glycated hemoglobin were associated with an OR of 1.16 (95% CI: 1.07–1.25) for VD deficiency in older Mexican adults. Although in this study, we only have data on glucose levels and self-reports of type 2 diabetes, in other studies these parameters are frequently used for the diagnosis of type 2 diabetes [43]. Our findings corroborate the association between type 2 diabetes and VD deficiency (OR 1.94 95% CI 1.27–2.94, data not shown).
To the best of our knowledge, this is the first study identifying the association of GC and VDR variants with VD deficiency in Mexican postmenopausal women. Vitamin D-binding protein (VDBP) is encoded by the GC gene. It binds to VD sterol metabolites for transporting them through the circulatory system towards target organs. Previous studies have demonstrated that the concentration of VD metabolites was strong and positively associated with serum VDBP levels [11,44,45,46]. In this study, the tag SNP rs2282679 on GC was associated with VD deficiency. Furthermore, we found almost a two fold risk increment for having VD deficiency in women who were homozygous for the G allele in the GC polymorphism rs2282679. Our data resembles previous associations of these SNPs through genome-wide association studies [11,31,47,48] in other ethnic populations, such as the Chinese population, where rs2298849 in the GC gene was significantly associated with serum 25(OH)D3 levels in postmenopausal women [33]. In addition, a positive association between plasma VDBP (Vitamin D Binding Protein) and 1,25(OH)2D concentration has been reported in a cohort from Denmark [49]. A recent report in non-Hispanic populations (White and Black), showed evidence that polymorphisms on the GC gene may be related to 25(OH)D status during pregnancy. The minor allele of the rs7041 SNP was related to increased 25(OH)D and the rs4588 SNP was associated with decreased 25(OH)D among pregnant women [50]. Although in our study, these polymorphisms were not observed to be associated with vitamin D levels, we cannot rule out the influence on VDBP levels in our population.
The VDR polymorphisms have been previously associated to VD status. Vitamin D receptor binds to 1-25-(OH)2D3, allowing transcription of genes that modulate VD metabolism [51]. A study in Hispanic individuals observed an inverse association between a VDR SNP (rs10783219) and 25(OH)D concentrations adjusted for age and gender (p = 0.004) [48]. Our results differ from Engelman et al. [30], as we did not observe an association between rs10783219 and VD deficiency in Mexican postmenopausal women. In contrast, an association was observed for the SNP rs4516035 (p = 0.012). This variant is located at a GATA binding site. The C allele is known to eliminate the GATA binding site, resulting in lower VDR promoter activity [52]. This observation leads us to postulate that the effect of some variants among the Mexican Mestizo population may be more significant than in Europeans. Additional studies in other populations are need to explore the extent of the association between these variants and VD deficiency.
We did not find significant association between variants in CYP2R1, CYP27B1, CYP24A1 and NADSYN1/DHCR7 genes and VD serum concentrations, while several studies have claimed to find association [52,53,54]. The negative results might be mainly attributed to ethnic differences between the Mexican population and other groups around the world. Additionally, this might be explained by the lacking functionality of these SNPs in Mexican populations.
On the other hand, studies regarding the relationship between VDR polymorphisms and BMD are inconclusive. A meta-analysis considering different populations showed an association between some VDR polymorphisms and low BMD [55]. The VDR ApaI polymorphism decreased the risk of osteoporosis in Caucasian postmenopausal women, while in Asian populations, VDR BsmI and VDR FokI were associated with an increased risk of osteoporosis. In this study, VDR gene polymorphisms (ApaI, BsmI, Cdx2, FokI and TaqI) were not analyzed. However, an association of the NADSYN1/DHCR7 gene with osteopenia/osteoporosis was identified. This gene encodes the 7-dehydrocholesterol reductase, which catalyzes the production of cholesterol from 7-dehydrocholesterol, using NADPH, in the de novo synthesis pathway of 25(OH)D in the skin. More studies in Mexican population with other polymorphisms on more genes involved in VD metabolism are needed to confirm or rule out the association with bone mineral density observed in other populations.
The highest frequency (29%) of the rs4516035-C allele observed among ethnic groups from the South East (SE) of Mexico agrees with previous studies, which have reported that Amerindian ancestry shows an increasing gradient from North to South. These findings contrast markedly with the distribution of VD deficiency along the country reported in the National Health and Nutrition Survey-2012 (ENSANUT-2012). The ENSANUT observed that the frequency of VD deficiency was higher in indigenous women from the South region of the country when compared to those from the North [8,38]. We hypothesized that the highest VD deficiency frequencies observed in the South of the country may be related to a genetic predisposition, low dietary intake of VD and the more intense pigmentation of the skin.
This study has some limitations that should be considered. First, it was based on a cross-sectional design and the analyzed SNPs were selected from previous reports conducted in Hispanic and European populations [11,30,31]. Second, it does not include measurements of 25(OH)D in Amerindian groups or sun exposure in Mexican postmenopausal women. However, we use the blood collection season as a proxy for sun exposure. Third, the women in the present study are part of a cohort study that included health workers of the Mexican Institute of Social Security, who are likely more educated and healthier than the general population of Mexico. Therefore, our findings related to sociodemographic and dietary variables cannot predict the association at the national level. However, we consider that the HWCS population could be representative of adults living in urban areas in Central Mexico. Fourth, after adjusting our analyses for multiple testing, no variant remained associated with VD deficiency, probably due to our limited sample size.

5. Conclusions

In conclusion, this is the first study of genetic variants related to VD deficiency in the Mexican population. We show that polymorphisms on the VDR and GC genes are associated with VD deficiency in postmenopausal women. These results suggest an important role for genetic variants in the high prevalence of VD deficiency observed in a sunny country with a strong indigenous ethnic background, such as Mexico.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-6643/10/9/1175/s1, Table S1: Allele frequencies of SNPs lying on selected genes involved in the metabolism of vitamin D, Table S2: Multiple logistic regressions for Vitamin D deficiency as the dependent variable, Table S3: Odds ratio for the association of vitamin D metabolism genetic variants and vitamin D deficiency in a sample of 400 postmenopausal women. Table S4: Multivariate ORs and 95% CIs for the association between DHCR7/NADSYN1 tag SNP and osteopenia/osteoporosis (n = 689), Figure S1: Flowchart of study population, Figure S2: Linkage disequilibrium map of single nucleotide polymorphisms in the VDR and GC genes.

Author Contributions

B.R.-P., J.S. and R.V.-C. were involved in the study design; B.R.-P., M.M.M.-A., A.H.-B., E.D.-G., M.C., M.Q., J.S. and R.V-C. conducted the research; M.Q., J.S., M.C., A.M-H. and L.O. were responsible of the data collection of the populations included in the study; B.R.-P. supervised the data collection and performed the statistical analyses; N.M., Y.N.F., A.H.-B., M.F. and A.D.Q.-S. contributed to the drafts of the manuscript; B.R.-P., J.S. and R.V.-C. responsibility for the final content of the manuscript; and all authors: substantial contributions to the manuscript, writing of the manuscript, and reading and approval of the final manuscript.

Funding

This work was supported by National Council of Science and Technology (grants INFR-2016-01-270405, SALUD-2010-C01-139795 and partially supported by INMEGEN (64-18/2011/I).

Acknowledgments

Berenice Rivera Paredez is a student from the Programa de Doctorado en Ciencias en Epidemiología, Instituto Nacional de Salud Pública (INSP), and received fellowship from the Mexican Council of Science and Technology (CONACYT). The authors acknowledge the staff of the Unidad de Investigación Epidemiológica y en Servicios de Salud-IMSS, Cuernavaca, Morelos. We also acknowledge the technical assistance provided for Jeny Flores Morales (National Institute of Genomic Medicine, INMEGEN).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ojeda-Lavín, A. Consejo Nacional de Población (CONAPO). Available online: http://www.conapo.gob.mx/es/CONAPO/Proyecciones_Datos (accessed on 2 April 2018).
  2. Clark, P.; Tamayo, J.A.; Cisneros, F.; Rivera, F.C.; Valdés, M. Epidemiology of osteoporosis in Mexico. Present and future directions. Rev. Investig. Clin. 2013, 65, 183–191. [Google Scholar]
  3. Clark, P.; Lavielle, P.; Franco-Marina, F.; Ramírez, E.; Salmerón, J.; Kanis, J.A.; Cummings, S.R. Incidence rates and life-time risk of hip fractures in Mexicans over 50 years of age: A population-based study. Osteoporos. Int. 2005, 16, 2025–2030. [Google Scholar] [CrossRef] [PubMed]
  4. Kuchuck, N.O.; Pluijm, S.M.; Van Schoor, N.M.; Smit, J.H.; Lips, P. Relationships of serum 25-hydroxyvitamin D to bone mineral density and serum parathyroid hormone and markers of bone turnover in older persons. J. Clin. Endocrinol. Metab. 2009, 94, 1244–1250. [Google Scholar] [CrossRef] [PubMed]
  5. Hwang, S.; Choi, H.S.; Kim, K.M.; Rhee, Y.; Lim, S.K. Association between serum 25-hydroxyvitramin D and bone mineral density and proximal geometry in Korean: The Korean National Health and Nutrition Examinations Surveys (KNHANES) 2008–2009. Osteoporos. Int. 2014, 26, 163–171. [Google Scholar] [CrossRef] [PubMed]
  6. Holick, M.F.; Chen, T.C. Vitamin D deficiency: A worldwide problem with health consequences. Am. J. Clin. Nutr. 2008, 87, S1080–S1086. [Google Scholar] [CrossRef] [PubMed]
  7. Elizondo-Montemayor, E.; Ugalde-Asas, P.; Serrano-González, M.; Cuello-García, C.; Borbolla-Escoboza, J. Serum 25-hyroxyvitamin D concentration, life factors and obesity in Mexican Children. Obesity 2010, 18, 1805–1811. [Google Scholar] [CrossRef] [PubMed]
  8. Flores, M.; Barquera, S.; Sánchez, L.; Macías, N.; Lozada, A.; Díaz, E.; Barquera, S. Concentraciones de Vitaminas D en niños, Adolescentes y Adultos Mexicanos Resultados de la ENSANUT-2006. Available online: http://www.cdi.salud.gob.mx:8080/BasesCDI/Archivos/EstudioseInvestigaciones/ReporteVitaminaD.pdf (accessed on 10 January 2018).
  9. Carrillo-Vega, M.F.; García-Peña, C.; Gutiérrez-Robledo, L.M.; Pérez-Zepeda, M.U. Vitamin D deficiency in older adults and its associated factors: A cross-sectional analysis of the Mexican Health and Aging Study. Arch. Osteoporos. 2017, 12, 8. [Google Scholar] [CrossRef] [PubMed]
  10. McMillan, A.; Hicks, J.; Isabella, C.; Higa, G.M. A critical analysis of the (near) legendary status of vitamin D. Expert Rev. Endocrinol. 2012, 7, 103–119. [Google Scholar] [CrossRef]
  11. Nair, R.; Maseeh, A. Vitamin D: The “sunshine” vitamin. J. Pharmacol. Pharmacother. 2012, 3, 118–126. [Google Scholar] [CrossRef] [PubMed]
  12. Hunter, D.; de Lange, M.; Snieder, H.; MacGregor, A.J.; Swaminathan, R.; Thakker, R.V.; Spector, T.D. Genetic contribution to bone metabolism, calcium excretion, and vitamin D and parathyroid hormone regulation. J. Bone Miner. Res. 2001, 16, 371–378. [Google Scholar] [CrossRef] [PubMed]
  13. Wang, T.J.; Zhang, F.; Richards, J.B.; Kestenbaum, B.; Meurs, J.B.; Berry, D.; Kiel, D.; Streeten, E.A.; Ohlsson, C.; Koller, D.L.; et al. Common genetic determinant of vitamin d insufficiency: A genome-wide association study. Lancet 2010, 376, 180–188. [Google Scholar] [CrossRef]
  14. Wang, W.; Ingles, S.A.; Torres-Mejía, G.; Stern, M.C.; Stanczyc, F.; Schwartz, G.G.; Nelson, D.O.; Fejerman, L.; Wolff, R.K.; Slattery, M.L.; et al. Genetic variants and non-genetic factors predict circulating vitamin D levels in Hispanic and non-Hispanic white women: The breast cancer disparities study. In. J. Mol. Epidemiol. Genet. 2014, 5, 31–46. [Google Scholar]
  15. Denova-Gutiérrez, E.; Flores, Y.N.; Gallegos-Carrillo, K.; Ramírez-Palacios, P.; Rivera-Paredez, B.; Muñoz-Aguirre, P.; Velázquez-Cruz, R.; Torres-Ibarra, L.; Meneses-León, J.; Méndez-Hernández, P.; et al. Health workers cohort study: Methods and study design. Salud Publica Mex. 2016, 58, 708–716. [Google Scholar] [CrossRef] [PubMed]
  16. Contreras-Cubas, C.; Sánchez-Hernández, B.E.; García-Ortiz, H.; Martínez-Hernández, A.; Barajas-Olmos, F.; Cid, M.; Mendoza-Caamal, E.C.; Centeno-Cruz, F.; Ortiz-Cruz, G.; Jiménez-López, J.C. Heterogenous distribution of mthfr gene variants among mestizos and diverse amerindian groups from mexico. PLoS ONE 2016, 11, e0163248. [Google Scholar] [CrossRef] [PubMed]
  17. Ojeda-Granados, C.; Panduro, A.; Gonzalez-Aldaco, K.; Sepulveda-Villegas, M.; Rivera-Iñiguez, I.; Roman, S. Tailoring nutritional advice for mexicans based on prevalence profiles of diet-related adaptive gene polymorphisms. J. Pers. Med. 2017, 7, 16. [Google Scholar] [CrossRef] [PubMed]
  18. Rubi-Castellanos, R.; Martinez-Cortes, G.; Munoz-Valle, J.F.; Gonzalez-Martin, A.; Cerda-Flores, R.M.; Anaya-Palafox, M.; Rangel-Villalobos, H. Pre-Hispanic Mesoamerican demography approximates the present-day ancestry of Mestizos throughout the territory of Mexico. Am. J. Phys. Anthropol. 2009, 139, 284–294. [Google Scholar] [CrossRef] [PubMed]
  19. Wolf, A.M.; Hunter, D.J.; Colditz, G.A.; Manson, J.E.; Stampfer, M.J.; Corsano, K.A.; Rosner, B.; Kriska, A.; Willett, W.C. Reproducibility and validity of a self-administered physical activity questionnaire. Int. J. Epidemiol. 1994, 23, 991–999. [Google Scholar] [CrossRef] [PubMed]
  20. Martínez, G.M.A.; López, F.C.; Varo, J.J.; Sánchez, V.A.; Martinez, J.A. Validation of the Spanish version of the physical activity questionnaire used in the Nurses “Health Study and the Health Professionals” Follow-up Study. Public Health Nutr. 2005, 8, 920–927. [Google Scholar] [CrossRef]
  21. Lohman, T.; Roche, A.; Martorell, R. Anthropometric Standardization Reference Manual; Human Kinetics: Champaign, IL, USA, 1998. [Google Scholar]
  22. National Center for Health Statistics. National Health and Nutrition Examination Survey. Available online: https://www.cdc.gov/nchs/data/nhanes/nhanes_07_08/manual_dexa.pdf consulted (accessed on 15 January 2018).[Green Version]
  23. WHO Scientific Group on the Assessment of Osteoporosis at Primary Health Care Level. Available online: http://www.who.int/chp/topics/Osteoporosis.pdf (accessed on 15 February 2018).
  24. Hernández-Avila, M.; Romieu, I.; Parra, S.; Hernández-Avila, J.; Madrigal, H.; Willett, W. Validity and reproducibility of a food frequency questionnaire to assess dietary intake of women living in Mexico City. Salud Publica Mex. 1998, 40, 133–140. [Google Scholar] [CrossRef] [PubMed]
  25. Hernández-Avila, M.; Resoles, M.; Parra, S.; Romieu, I. Software: Sistema de Evaluación de Hábitos Nutricionales y Consumo de Nutrimentos. Available online: http://www.insp.mx/snut2003/ (accessed on 18 August 2017).
  26. Tate, J.R.; Berg, K.; Couderc, R.; Dati, F.; Kostner, G.M.; Marcovina, S.M.; Rifai, N.; Sakurabayashi, I.; Steinmetz, A. International Federation of Clinical Chemistry and Laboratory Medicine (IFCC) Standardization Project for the Measurement of Lipoprotein(a). Phase 2: Selection and properties of a proposed secondary reference material for lipoprotein(a). Clin. Chem. Lab. Med. 1999, 37, 949–958. [Google Scholar] [CrossRef] [PubMed]
  27. Freeman, J.; Wilson, K.; Spears, R.; Shalhoub, V.; Sibley, P. Performance evaluation of four 25-hydroxyvitamin D assays to measure 25-hydroxyvitamin D2. Clin. Biochem. 2015, 48, 1097–1104. [Google Scholar] [CrossRef] [PubMed]
  28. Flores, A.; Flores, M.; Macias, N.; Hernández-Barrera, L.; Rivera, M.; Contreras, A.; Villalpando, S. Vitamin D deficiency is common and is associated with overweight in Mexican children aged 1–11 years. Public Health Nutr. 2017, 20, 1807–1815. [Google Scholar] [CrossRef] [PubMed]
  29. Ahn, J.; Yu, K.; Stolzenberg-Solomon, R.; Simon, K.C.; McCullough, M.L.; Gallicchio, L.; Jacobs, E.J.; Ascherio, A.; Helzlsouer, K.; Jacobs, K.B.; et al. Genome-wide association study of circulating vitamin D levels. Hum. Mol. Genet. 2010, 19, 2739–2745. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Engelman, C.D.; Fingerlin, T.E.; Langefeld, C.D.; Hicks, P.J.; Rich, S.S.; Wagenknecht, L.E.; Bowden, D.W.; Norris, J.M. Genetic and environmental determinants of 25-hydroxyvitamin D and 1,25-dihydroxyvitamin D levels in Hispanic and African Americans. J. Clin. Endocrinol. Metab. 2008, 93, 3381–3388. [Google Scholar] [CrossRef] [PubMed]
  31. Velázquez-Cruz, R.; García-Ortiz, H.; Castillejos-López, M.; Quiterio, M.; Valdés-Flores, M.; Orozco, L.; Villarreal-Molina, T.; Salmerón, J. WNT3A gene polymorphisms are associated with bone mineral density variation in postmenopausal mestizo women of an urban Mexican population: Findings of a pathway-based high-density single nucleotide screening. Age 2014, 36, 9635. [Google Scholar] [CrossRef] [PubMed]
  32. Dudbridge, F. Power and predictive accuracy of polygenic risk scores. PLoS Genet. 2013, 9, e1003348. [Google Scholar] [CrossRef]
  33. Xu, W.; Sun, J.; Wang, W.; Wang, X.; Jiang, Y.; Huang, W.; Zheng, X.; Wang, Q.; Ning, Z.; Pei, Y.; et al. Association of genetic variants of vit D binding protein (DBP/GC) and of the enzyme catalyzing its 25-hydroxylation (DCYP2R1) and serum vit D in postmenopausal women. Hormones 2014, 13, 345–352. [Google Scholar] [CrossRef] [PubMed]
  34. StataCorp. Stata Statistical Software: Release 13. 2013. Available online: https://www.stata.com/ (accessed on 1 January 2017).
  35. Barrett, J.C.; Fry, B.; Maller, J.; Daly, M.J. Haploview: Analysis and visualization of LD and haplotype maps. Bioinformatics 2005, 21, 263–265. [Google Scholar] [CrossRef] [PubMed]
  36. Moreno-Estrada, A.; Gignoux, C.R.; Fernández-López, J.C.; Zakharia, F.; Sikora, M.; Contreras, A.V.; Acuña-Alonzo, V.; Sandoval, K.; Eng, C.; Romero-Hidalgo, S.; et al. Human genetics. The genetics of Mexico recapitulates Native American substructure and affects biomedical traits. Science 2014, 344, 1280–1285. [Google Scholar] [CrossRef] [PubMed]
  37. Contreras-Manzano, A.; Villalpando, S.; Robledo-Pérez, R. Vitamin D status by sociodemographic factors and body mass index in Mexican women at reproductive age. Salud Publica Mex. 2017, 59, 518–525. [Google Scholar] [CrossRef] [PubMed]
  38. Forrest, K.Y.; Stuhldreher, W.L. Prevalence and correlates of vitamin D deficiency in US adults. Nutr. Res. 2011, 31, 48–54. [Google Scholar] [CrossRef] [PubMed]
  39. Lips, P.; Hosking, D.; Lippuner, K.; Norquist, J.M.; Wehren, L.; Maalouf, G.; Ragi-Eis, S.; Chandler, J. The prevalence of vitamin D inadequacy amongst women with osteoporosis: An international epidemiological investigation. J. Intern. Med. 2006, 260, 245–254. [Google Scholar] [CrossRef] [PubMed]
  40. Yu, H.J.; Kwon, M.J.; Woo, H.Y.; Park, H. Analysis of 25-hydroxyvitamin D status according to age, gender, and seasonal variation. J. Clin. Lab. Anal. 2016, 30, 905–911. [Google Scholar] [CrossRef] [PubMed]
  41. Rontoyanni, V.G.; Avila, J.C.; Kaul, S.; Wong, R.; Veeranki, S.P. Association between Obesity and Serum 25(OH)D Concentrations in Older Mexican Adults. Nutrients 2017, 9, 97. [Google Scholar] [CrossRef] [PubMed]
  42. Pathak, K.; Soares, M.J.; Calton, E.K.; Zhao, Y.; Hallett, J. Vitamin D supplementation and body weight status: A systematic review and meta-analysis of randomized controlled trials. Obes. Rev. 2014, 15, 528–537. [Google Scholar] [CrossRef] [PubMed]
  43. Schneider, A.L.; Pankow, J.S.; Heiss, G.; Selvin, E. Validity and reliability of self-reported diabetes in the atherosclerosis risk in communities study. Am. J. Epidemiol. 2012, 176, 738–743. [Google Scholar] [CrossRef] [PubMed]
  44. Tint, G.S.; Irons, M.; Elias, E.R.; Batta, A.K.; Frieden, R.; Chen, T.S.; Salen, G. Defective cholesterol biosynthesis associated with the Smith-Lemli-Opitz syndrome. N. Engl. J. Med. 1994, 330, 107–113. [Google Scholar] [CrossRef] [PubMed]
  45. Heaney, R.P. Vitamin D and calcium interactions: Functional outcomes. Am. J. Clin. Nutr. 2008, 88, 541S–544S. [Google Scholar] [CrossRef] [PubMed]
  46. Lasky-Su, J.; Lyon, H.N.; Emilsson, V.; Heid, I.M.; Molony, C.; Raby, B.A.; Lazarus, R.; Klanderman, B.; Soto-Quiros, M.E.; Avila, L.; et al. On the replication of genetic associations: Timing can be everything! Am. J. Hum. Genet. 2008, 82, 849–858. [Google Scholar] [CrossRef] [PubMed]
  47. Jiang, X.; O′Reilly, P.F.; Aschard, H.; Hsu, Y.H.; Richards, J.B.; Dupuis, J.; Ingelsson, E.; Karasik, D.; Pilz, S.; Berry, D.; et al. Genome-wide association study in 79,366 European-ancestry individuals informs the genetic architecture of 25-hydroxyvitamin D levels. Nat. Commun. 2018, 9, 260. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Hong, J.; Hatchell, K.E.; Bradfield, J.P.; Andrew, B.; Alessandra, C.; Chao-Qiang, L.; Langefeld, C.D.; Lu, L.; Lu, Y.; Lutsey, P.L.; et al. Trans-ethnic evaluation identifies novel low frequency loci associated with 25-hydroxyvitamin D. J. Clin. Endocrinol. Metab. 2018, 103, 1380–1392. [Google Scholar] [CrossRef] [PubMed]
  49. Lauridsen, A.L.; Vestergaard, P.; Hermann, A.P.; Brot, C.; Heickendorff, L.; Mosekilde, L.; Nexo, E. 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, 15–22. [Google Scholar] [CrossRef] [PubMed]
  50. Baca, K.M.; Govil, M.; Zmuda, J.M.; Simhan, H.N.; Marazita, M.L.; Bodnar, L.M. Vitamin D metabolic loci and vitamin D status in Black and White pregnant women. Eur. J. Obstet. Gynecol. Reprod. Biol. 2018, 220, 61–68. [Google Scholar] [CrossRef] [PubMed]
  51. Barry, E.L.; Rees, J.R.; Peacock, J.L.; Mott, L.A.; Amos, C.I.; Bostick, R.M.; Figueiredo, J.C.; Ahnen, D.J.; Bresalier, R.S.; Burke, C.A.; et al. Genetic variants in CYP2R1, CYP24A1, and VDR modify the efficacy of vitamin D3 supplementation for increasing serum 25-hydroxyvitamin D levels in a randomized controlled trial. J. Clin. Endocrinol. Metab. 2014, 99, E2133–E2137. [Google Scholar] [CrossRef] [PubMed]
  52. Fang, Y.; van Meurs, J.B.J.; d’Alesio, A. Promoter and 3′-untranslated-region haplotypes in the vitamin D receptor gene predispose to osteoporotic fracture: The rotterdam study. Am. J. Hum. Genet. 2005, 77, 807–823. [Google Scholar] [CrossRef] [PubMed]
  53. Moon, R.J.; Harvey, N.C.; Cooper, C.; D′Angelo, S.; Curtis, E.M.; Crozier, S.R.; Barton, S.J.; Robinson, S.M.; Godfrey, K.M.; Graham, N.J.; et al. Response to antenatal cholecalciferol supplementation is associated with common vitamin d-related genetic variants. J. Clin. Endocr. Metab. 2017, 102, 2941–2949. [Google Scholar] [CrossRef] [PubMed]
  54. Sollid, S.T.; Hutchinson, M.Y.; Fuskevag, O.M.; Joakimsen, R.M.; Jorde, R. Large individual differences in serum 25-hydroxyvitamin D response to vitamin D supplementation: Effects of genetic factors, body mass index, and baseline concentration. Results from a randomized controlled trial. Horm. Metab. Res. 2016, 48, 27–34. [Google Scholar] [CrossRef] [PubMed]
  55. Zhang, L.; Yin, X.; Wang, J.; Xu, D.; Wang, Y.; Yang, J.; Tao, Y.; Zhang, S.; Feng, X.; Yan, C. Associations between VDR gene polymorphisms and osteoporosis risk and bone mineral density in postmenopausal women: A systematic review and meta-analysis. Sci. Rep. 2018, 8, 981. [Google Scholar] [CrossRef] [PubMed]
Nutrients 10 01175 g001 550
Figure 1. Dose-response relationship between the number of risk alleles in rs4516035, rs2282679 and vitamin D deficiency. * Lines indicate 95% confidence intervals. (A) Prevalence of vitamin D deficiency status by GRS. (B) Odds Ratio of vitamin D deficiency status by GRS.
Figure 1. Dose-response relationship between the number of risk alleles in rs4516035, rs2282679 and vitamin D deficiency. * Lines indicate 95% confidence intervals. (A) Prevalence of vitamin D deficiency status by GRS. (B) Odds Ratio of vitamin D deficiency status by GRS.
Nutrients 10 01175 g001
Nutrients 10 01175 g002 550
Figure 2. Geographic distribution of the allelic frequencies (T/C) for the rs4516035 polymorphism in Mexican Amerindian groups. MEX, Mexico State; MOR, Morelos; OAX, Oaxaca; PUE, Puebla; SLP, San Luis Potosi, VER, Veracruz. (A) Distribution of the allelic frequencies for the rs4516035 polymorphism within each geographic region. (B) Allelic frequencies for the rs4516035 polymorphism in the Amerindian groups.
Figure 2. Geographic distribution of the allelic frequencies (T/C) for the rs4516035 polymorphism in Mexican Amerindian groups. MEX, Mexico State; MOR, Morelos; OAX, Oaxaca; PUE, Puebla; SLP, San Luis Potosi, VER, Veracruz. (A) Distribution of the allelic frequencies for the rs4516035 polymorphism within each geographic region. (B) Allelic frequencies for the rs4516035 polymorphism in the Amerindian groups.
Nutrients 10 01175 g002
Nutrients 10 01175 g003 550
Figure 3. Geographic distributions of the allelic frequencies (G/C) for the rs2298850 polymorphism, (proxy for rs2282679) in the Mexican Amerindian groups. MEX, State of Mexico; MOR, Morelos; OAX, Oaxaca; PUE, Puebla; SLP, San Luis Potosi, VER, Veracruz. (A) Distribution of the allelic frequencies for the rs2298850 polymorphism within each geographic region. (B) Allelic frequencies for the rs2298850 polymorphism in the Amerindian groups.
Figure 3. Geographic distributions of the allelic frequencies (G/C) for the rs2298850 polymorphism, (proxy for rs2282679) in the Mexican Amerindian groups. MEX, State of Mexico; MOR, Morelos; OAX, Oaxaca; PUE, Puebla; SLP, San Luis Potosi, VER, Veracruz. (A) Distribution of the allelic frequencies for the rs2298850 polymorphism within each geographic region. (B) Allelic frequencies for the rs2298850 polymorphism in the Amerindian groups.
Nutrients 10 01175 g003
Table
Table 1. Characteristics of postmenopausal women of the Health Workers Cohort Study (2010).
Table 1. Characteristics of postmenopausal women of the Health Workers Cohort Study (2010).
Characteristicsn = 689
Age (years) *62.2 (8.8)
45−59, n (%)279 (40.9)
60−74, n (%)340 (49.9)
>74, n (%)63 (9.2)
BMI (kg/m2) *28.1 (4.7)
Nutritional Status, n (%)
Overweight, λ306 (44.9)
Obesity, λ204 (29.9)
Waist circumference (cm) *95.9 (11.3)
Body fat proportion *45.6 (6.1)
Leisure time physical activity (min/day) *20.9 (27.3)
Active (≥30 min/day), %186 (27.3)
Serum 25(OH)D (ng/mL) *21.0 (5.7)
25(OH)D Deficiency (<20ng/mL), n (%)298 (43.8)
Bone Femoral neck Density (g/cm2) *0.87 (0.13)
Bone Femoral neck status, n (%)
Osteopenia224 (32.8)
Osteoporosis32 (4.7)
Bone Spine Density (g/cm2) *1.00 (0.16)
Bone Spine status, n (%)
Osteopenia324 (47.5)
Osteoporosis165 (24.2)
Season of blood collection, n (%)
Winter112 (16.4)
Spring412 (60.4)
Summer113 (16.6)
Autumn45 (6.6)
Diet
Total energy intake (kcal/day) *1805 (854)
Alcohol (g/day) *1.6 (4.2)
Vitamin D intake (IU/day) *199 (157)
Calcium intake (mg/day) *918 (545)
Calcium supplements, n (%)194 (28.5)
Medication use, n (%)
Hormone replacement therapy, n (%)44 (6.5)
* Mean (standard deviation). λ Overweight (25 to <30 kg/m2) and Obesity (≥30 kg/m2).
Table
Table 2. Characteristics of postmenopausal women of the Health Workers Cohort Study by vitamin D status (2010).
Table 2. Characteristics of postmenopausal women of the Health Workers Cohort Study by vitamin D status (2010).
Non-Deficient: 25-OH-D ≥20 ng/mL
n = 386 (56.2%)		Deficient: 25-OH-D <20 ng/mL
n = 301 (43.8%)		p Value
Age(years) *61.2 (8.4)63.4 (9.2)0.0008
Age categories, n (%)
45–59180 (46.6)101 (33.6)0.006
60–74179 (46.4)163 (54.2)0.04
>7427 (7.0)37 (12.3)0.02
BMI (kg/m2) *28.1 (4.7)28.2 (4.8)0.76
Nutritional Status, n (%)
Overweight, λ175 (45.3)133 (44.2)0.71
Obesity, λ114 (29.5)91 (30.2)0.84
Waist circumference (cm) *95.5 (11.1)96.4 (11.6)0.30
Body fat proportion *45.7 (6.2)45.4 (6.1)0.53
Leisure time physical activity (min/day) *20.5 (27.7)19.4 (26.0)0.60
Active (≥30 min/day), %28.525.80.43
Femoral neck (g/cm2) *0.89 (0.13)0.86 (0.13)0.004
Femoral status, n (%)
Osteopenia121 (31.4)107 (35.7)0.23
Osteoporosis17 (4.4)16 (5.3)0.57
Spine BMD (g/cm2) *1.01 (0.16)0.99 (0.16)0.004
Spine status, n (%)
Osteopenia187 (48.6)139 (46.3)0.55
Osteoporosis83 (21.6)84 (28.0)0.05
Season of blood collection, n (%)
Winter54 (14.8)55 (18.3)0.22
Spring249 (64.5)163 (54.1)0.006
Summer58 (15.0)61 (20.3)0.07
Autumn22 (5.7)22 (7.3)0.40
Diet
Total energy intake (kcal/day) *1807 (826)1814 (895)0.92
Alcohol intake (g/day) *1.5 (3.5)1.8 (4.9)0.28
Vitamin D intake (IU/day) *209 (160)188 (151)0.02
Calcium intake (mg/day) *941 (537)890 (554)0.22
Calcium supplements, n (%)116 (30.1)79 (26.3)0.29
Medication use
Hormone replacement therapy, n (%)23 (6.0)21 (7.0)0.60
Calcium, n (%)9 (2.3)5 (1.7)0.55
* Mean (standard deviation). p values from t-test (continuous variables) or chi2 test (categorical variables). λ Overweight (25 to <30 kg/m2) and Obesity (≥ 30 kg/m2).
Table
Table 3. Odds ratios for the association of vitamin D metabolism genetic variants and vitamin D deficiency in a sample of 689 postmenopausal women.
Table 3. Odds ratios for the association of vitamin D metabolism genetic variants and vitamin D deficiency in a sample of 689 postmenopausal women.
GeneSNPGenotypeNot Deficient *
n (%)		Deficient *
n (%)		Unadjuste
d OR (CI 95%)		p ValueAdjusted Model
OR (CI 95%)		p Value
VDRrs10783219TT135 (35.3)94 (31.5)1.0 1.0
TA170 (44.5)145 (48.7)1.22 (0.87–1.73)0.2481.20 (0.84–1.72)0.315
AA77 (20.2)59 (19.8)1.10 (0.72–1.69)0.6621.13 (0.72–1.75)0.597
Additive model1.07 (0.87–1.32)0.5371.10 (0.89–1.38)0.377
VDRrs4516035TT219 (57.5)146 (49.0)1.0 1.0
TC141 (37.0)130 (43.6)1.38 (1.01–1.90)0.0451.41 (1.01–1.95)0.041
CC21 (5.5)22 (7.4)1.57 (0.83–2.96)0.1621.75 (0.91–3.35)0.091
Additive model1.31 (1.02–1.68)0.0321.40 (1.08–4.43)0.012
GCrs2282679TT247 (64.6)164 (55.0)1.0 1.0
TG127 (33.3)119 (39.9)1.41 (1.03–1.94)0.0341.48 (1.07–2.06)0.019
GG8 (2.1)15 (2.0)2.82 (1.17–6.81)0.0213.00 (1.20–7.49)0.019
Additive model1.50 (1.14–1.97)0.0031.53 (1.15–2.04)0.003
* Non-deficient: serum 25(OH)D ≥20 ng/mL; deficient: serum 25(OH)D <20 ng/mL. Adjusted model: adjusted for age (<60, 60–74, >74 years), body fat proportion (tertiles), vitamin D intake IU/day (tertiles), season of blood collection (winter, spring, summer, autumn) and leisure time physical activity (>30 min/day).
Table
Table 4. Genotype frequencies of GC and VDR variants in 37 Mexican Amerindian groups.
Table 4. Genotype frequencies of GC and VDR variants in 37 Mexican Amerindian groups.
PopulationsGenotypic Frequencies
Geographical RegionEthnic Grouprs2298850 (Proxy for rs2282679)rs4516035
GGGCCCTTCTCC
CountFreqCountFreqCountFreqCountFreqCountFreqCountFreq
North (N)Guarijio20.2240.4430.3380.8910.1100
Mayo40.8010.200030.6020.4000
Seri71.0000081.00000
Tarahumara190.56150.4400330.9710.0300
Yaqui41.0000030.7510.2500
Total360.72190.2230.067550.8450.1600
Central West (CW)Purepecha70.7820.220070.7820.2200
Total70.7820.220070.7820.2200
Central East (CE)Huasteco70.7030.300070.7030.3000
Matlaltzinca80.8010.1010.1060.6040.4000
Mazateco PUE40.4450.560040.4450.5600
Nahuatl MOR80.8910.11000080.8910.11
Nahuatl PUE50.7120.290040.5730.4300
Nahuatl SLP61.0000040.6720.3300
Otomi280.7880.2200220.61110.3130.08
Pame80.8910.110070.7030.3000
Popoluca de la Sierra70.8810.120070.8810.1200
Tlahuica51.0000051.00000
Totonaco PUE60.5550.450060.5550.4500
Totonaco VER70.7030.300060.6040.4000
Total990.79300.20610.009780.601490.38140.018
South (S)Chinanteco50.5050.500070.7030.3000
Chontal OAX30.5020.3310.1730.5020.3310.17
Huave60.8610.140050.7210.1410.14
Mazateco OAX40.5720.2910.1450.7120.2900
Mixe80.8910.10040.4450.5600
Mixteco Alto80.7330.270060.5550.4500
Mixteco Costa50.5040.4010.1070.7030.3000
Zapoteco60.7520.250030.3820.2430.38
Total450.66200.2930.05400.587230.32850.086
South East (SE)Chuj50.8310.170020.3320.3320.33
Jakalteco30.6020.400020.4030.6000
Kanjobal90.9010.100020.2070.7010.10
Kaqchikel60.6730.330050.5630.3310.11
Mam40.4440.4410.1270.7820.2200
Maya120.8030.200090.6040.2720.13
Mocho50.8310.170050.830010.17
Tojolabal50.6330.370050.6330.3700
Tseltal30.6020.400010.2030.6010.20
Tsotsil50.6330.370060.7520.2500
Zoque30.6720.330040.8010.2000
Total600.68250.3110.01480.552300.35380.095
Mexico State; MOR, Morelos; OAX, Oaxaca; PUE, Puebla; SLP, San Luis Potosi; VER, Veracruz.

Share and Cite

MDPI and ACS Style

Rivera-Paredez, B.; Macías, N.; Martínez-Aguilar, M.M.; Hidalgo-Bravo, A.; Flores, M.; Quezada-Sánchez, A.D.; Denova-Gutiérrez, E.; Cid, M.; Martínez-Hernández, A.; Orozco, L.; et al. Association between Vitamin D Deficiency and Single Nucleotide Polymorphisms in the Vitamin D Receptor and GC Genes and Analysis of Their Distribution in Mexican Postmenopausal Women. Nutrients 2018, 10, 1175. https://doi.org/10.3390/nu10091175

AMA Style

Rivera-Paredez B, Macías N, Martínez-Aguilar MM, Hidalgo-Bravo A, Flores M, Quezada-Sánchez AD, Denova-Gutiérrez E, Cid M, Martínez-Hernández A, Orozco L, et al. Association between Vitamin D Deficiency and Single Nucleotide Polymorphisms in the Vitamin D Receptor and GC Genes and Analysis of Their Distribution in Mexican Postmenopausal Women. Nutrients. 2018; 10(9):1175. https://doi.org/10.3390/nu10091175

Chicago/Turabian Style

Rivera-Paredez, Berenice, Nayeli Macías, Mayeli M. Martínez-Aguilar, Alberto Hidalgo-Bravo, Mario Flores, Amado D. Quezada-Sánchez, Edgar Denova-Gutiérrez, Miguel Cid, Angelica Martínez-Hernández, Lorena Orozco, and et al. 2018. "Association between Vitamin D Deficiency and Single Nucleotide Polymorphisms in the Vitamin D Receptor and GC Genes and Analysis of Their Distribution in Mexican Postmenopausal Women" Nutrients 10, no. 9: 1175. https://doi.org/10.3390/nu10091175

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop