Hostname: page-component-7c8c6479df-r7xzm Total loading time: 0 Render date: 2024-03-19T08:18:48.103Z Has data issue: false hasContentIssue false

Genetic and environmental predictors of serum 25-hydroxyvitamin D concentrations among middle-aged and elderly Chinese in Singapore

Published online by Cambridge University Press:  15 May 2012

Kim Robien*
Affiliation:
Division of Epidemiology and Community Health, School of Public Health, University of Minnesota, 1300 South 2nd Street, Suite 300, Minneapolis, MN 55455, USA Masonic Cancer Center, University of Minnesota, MMC 807 Mayo, 8807A, 420 Delaware Street, SE, Minneapolis, MN 55455, USA
Lesley M. Butler
Affiliation:
Environmental and Radiological Health Sciences, Colorado State University, 148B EH Building, 1681 Campus Delivery, Fort Collins, CO 80523-1681, USA
Renwei Wang
Affiliation:
University of Pittsburgh Cancer Institute, UPMC Cancer Pavillion, Suite 4C-464, 5150 Center Avenue, Pittsburgh, PA 15232, USA
Kenneth B. Beckman
Affiliation:
Biomedical Genomics Center, University of Minnesota, MMC 501 Mayo, 8501, 420 Delaware Street, SE, Minneapolis, MN 55455, USA
Dinesha Walek
Affiliation:
Biomedical Genomics Center, University of Minnesota, MMC 501 Mayo, 8501, 420 Delaware Street, SE, Minneapolis, MN 55455, USA
Woon-Puay Koh
Affiliation:
Saw Swee Hock School of Public Health, National University of Singapore, MD3, 16 Medical Drive, Singapore117597, Singapore
Jian-Min Yuan
Affiliation:
University of Pittsburgh Cancer Institute, UPMC Cancer Pavillion, Suite 4C-464, 5150 Center Avenue, Pittsburgh, PA 15232, USA Department of Epidemiology, Graduate School of Public Health, University of Pittsburgh, Pittsburgh, PA, USA
*
*Corresponding author: Dr K. Robien, fax +1 612 624 0315, E-mail: robie004@umn.edu
Rights & Permissions [Opens in a new window]

Abstract

Vitamin D is known for maintaining Ca homeostasis and bone structure, and may also decrease susceptibility to chronic and infectious diseases. However, data on vitamin D status and its predictors among Southeast Asian populations are limited. We evaluated the distribution and determinants (genetic and environmental) of serum 25-hydroxyvitamin D (25(OH)D) concentrations among 504 middle-aged and elderly participants (aged 45–74 years) in the Singapore Chinese Health Study. Data on dietary and other lifestyle factors were collected by trained interviewers. Serum 25(OH)D concentrations and genetic polymorphisms in vitamin D metabolism pathway enzymes (cytochrome P450 (CYP) 2R1, 3A4, 27B1, 24A1; vitamin D binding protein (also known as group-specific component, GC); and vitamin D receptor) were measured using stored biospecimens. Mean 25(OH)D concentration was 68·8 nmol/l. Serum 25(OH)D concentrations were positively associated with dietary vitamin D intake, and inversely associated with hours spent sitting at work. BMI was not associated with 25(OH)D concentrations. CYP2R1 rs10741657, rs12794714, rs1993116; CYP3A4 rs2242480; and GC rs4588, rs7041, rs16847015, rs2298849 were statistically significantly associated with 25(OH)D concentrations. Individuals with the Gc2-2 haplotype (rs4588AA/rs7041TT) had statistically significantly lower 25(OH)D concentrations compared to all other Gc haplotypes (P-trend < 0·001). The majority of participants (86 %) had 25(OH)D concentrations ≥ 50 nmol/l, which is consistent with the 2011 Institute of Medicine (US) recommendation for bone health, and 32 % had concentrations of ≥ 75 nmol/l that are thought to be required for broader health effects. Dietary vitamin D intake, hours spent indoors at work and genetic variation in CYP2R1, CYP3A4 and GC are significant predictors of 25(OH)D concentrations among Singapore Chinese.

Type
Full Papers
Copyright
Copyright © The Authors 2012

Vitamin D (as the 1,25-dihydroxyvitamin D metabolite) is a steroid hormone that is well known for its role in maintaining Ca homeostasis and normal bone structure. Recent evidence suggests that in addition to Ca homeostasis, the vitamin may also play a role in a variety of other physiological processes such as modulation of inflammatory pathways(Reference Cohen-Lahav, Douvdevani and Chaimovitz1) and susceptibility to diabetes(Reference Boucher2), cancer(3), and infectious(Reference Baeke, Takiishi and Korf4) and cardiovascular(Reference Reddy Vanga, Good and Howard5) diseases. Thus, the nutrient could play a significant role in public health.

In the USA, the Institute of Medicine recently proposed ≥ 50 nmol/l (20 ng/ml) as the definition of vitamin D adequacy based solely on requirements to optimise bone health, due to a lack of data to support recommendations for the prevention of other disease endpoints(Reference Ross, Taylor, Yaktine and Del Valle6). However, many leading vitamin D researchers continue to recommend serum 25-hydroxyvitamin D (25(OH)D) concentrations of ≥ 75 nmol/l (30 ng/ml) to achieve the broader health benefits(Reference Vieth, Bischoff-Ferrari and Boucher7Reference Schwalfenberg and Whiting10).

Season, UVB exposure, skin pigmentation, age, race, sex, obesity and dietary/supplemental vitamin D intake have all been previously reported to influence serum 25(OH)D concentrations(Reference Holick11). However, the effect of genetic variation in the vitamin D synthesis and metabolism pathway on circulating concentrations is less well understood. Vitamin D enters the circulation through the activation of vitamin D precursors by UV radiation on the skin to produce cholecalciferol, or via absorption of dietary or supplemental ergo- or cholecalciferol from the intestinal tract. It is then converted to 25(OH)D via 25-hydroxylases (cytochrome P450 (CYP) 2R1, 27A1 and 3A4) in the liver(Reference Holick11). Further hydroxylation of 25(OH)D via 1α-hydroxylase (CYP27B1) in the kidney or at the local tissue level produces 1,25-dihydroxycholecalciferol(Reference Holick11). Catabolism of vitamin D metabolites occurs via 24-hydroxylase (CYP24A1)(Reference Holick11). Vitamin D binding protein (also known as group-specific component, GC) is the transport protein for vitamin D metabolites in circulation(Reference Speeckaert, Huang and Delanghe12). Genetic variation in any of these steps has the potential to alter serum 25(OH)D concentrations.

Previous studies have identified SNP in the vitamin D receptor (VDR)(Reference Ma, Stampfer and Gann13Reference Smolders, Damoiseaux and Menheere20), CYP2R1 (Reference Ramos-Lopez, Bruck and Jansen21Reference Bu, Armas and Lappe26), CYP27B1 (Reference Orton, Morris and Herrera22, Reference Ramos-Lopez, Kahles and Weber23, Reference Hypponen, Berry and Wjst27) and GC (Reference Engelman, Fingerlin and Langefeld19, Reference Wang, Zhang and Richards24Reference Bu, Armas and Lappe26, Reference Kurylowicz, Ramos-Lopez and Bednarczuk28Reference Janssens, Bouillon and Claes35) genes to be associated with alterations in serum 25(OH)D concentrations. These studies were primarily conducted among Caucasian populations living at higher latitudes with significant seasonal variation in UV exposure, and only a small number of studies have considered environmental vitamin D exposures or personal characteristics as potential effect modifiers.

Data on the distribution and determinants of serum 25(OH)D concentrations in Southeast Asian populations are limited. In this cross-sectional observational study, we evaluated the distribution of serum 25(OH)D concentrations, and identified genetic, dietary and lifestyle predictors of serum 25(OH)D concentrations among middle-aged and older Chinese men and women in Singapore. As Singapore is 1° north of the equator, this study population provides a unique opportunity to evaluate the factors associated with vitamin D status in the absence of seasonal variation in UV exposure, which confounds studies conducted at higher latitudes.

Methods

The Singapore Chinese Health Study (SCHS) is a population-based prospective cohort study of 63 257 Singapore Chinese men and women (aged 45–74 years) that was assembled between 1993 and 1998 to elucidate the role of diet and genetic factors in the causation of human cancer. Participants in the study were recruited from among permanent residents or citizens of Singapore who resided in government-built housing estates, and were one of the two major dialect groups of Singapore Chinese, Hokkien or Cantonese. At recruitment, each study subject was interviewed in person by a trained interviewer using a structured questionnaire that included questions on lifestyle, health behaviours and sociodemographic factors, as well as a 165-item FFQ. The FFQ was specifically designed to assess the dietary habits of Chinese in Singapore, and was subsequently validated using multiple 24 h dietary recalls(Reference Hankin, Stram and Arakawa36). However, because vitamin D is technically challenging to measure in food(Reference Byrdwell, Devries and Exler37), data on vitamin D contained in the nutrient databases used to analyse FFQ for epidemiological studies (including the US Department of Agriculture database used for the SCHS vitamin D analysis) are known to be incomplete(Reference Holden, Lemar and Exler38). Thus, our estimates probably underestimate the participants' actual dietary vitamin D intake. The only form of supplemental vitamin D intake to be assessed was cod liver oil; data on frequency of use over the year before the interview were collected and considered for use in this analysis. Participants were not asked specifically about time spent indoors v. outdoors; however, participants were asked about the average number of hours spent ‘sitting at work’ each day, and hours spent doing ‘vigorous work such as moving heavy furniture, loading or unloading trucks, shovelling or equivalent manual labour’. Responses to these questions were included in our analyses as surrogates for time spent indoors (‘sitting at work’) and outdoors (‘vigorous work’).

Beginning in April 1994, a random 3 % sample of cohort participants was also asked to provide blood or buccal cells, and spot urine samples as a pilot study to determine the feasibility of a larger biospecimen collection effort within the cohort. Details of the biospecimen collection, processing and storage procedures have been described(Reference Koh, Yuan and Sun39). The first 504 SCHS participants who provided biospecimens were included in this substudy. This study was conducted according to the guidelines laid down in the Declaration of Helsinki and all procedures involving human subjects were approved by the Institutional Review Boards of the National University of Singapore and the University of Minnesota. Written informed consent was obtained from all participants.

Serum 25-hydroxyvitamin D assay

Quantitative determination of serum 25(OH)D concentrations was performed using a direct, competitive chemiluminescence immunoassay on the DiaSorin LIAISON platform(Reference Ersfeld, Rao and Body40) by Heartland Assays, a laboratory participating in the Vitamin D External Quality Assessment Scheme(Reference Carter, Carter and Gunter41). This assay measures total 25(OH)D (both the ergocalciferol (D2) and cholecalciferol (D3) derived 25(OH)D metabolites; the assay does not discriminate between the two forms). Two types of blinded controls were randomly interspersed among the study samples for quality control purposes: (1) pooled blood samples (n 10), and (2) vitamin D standard (standard reference material 972, Vitamin D in Human Serum) from the US National Institute of Standards and Technology (NIST)(Reference Tai, Bedner and Phinney42) (n 16). Mean inter- and intra-batch CV for the 25(OH)D concentrations in the pooled blood samples was 6·1 and 5·4 %, respectively, which is consistent with those of previously published reports(Reference Giovannucci, Liu and Rimm43Reference Cauley, Lacroix and Wu47). For the NIST standard level 1 (58·0 nmol/l 25(OH)D3 in unaltered human serum), the mean inter-batch CV was 6·2 %. Similarly, for the NIST standard level 2 (4·18 nmol/l 25(OH)D2 and 30·0 nmol/l 25(OH)D3 in diluted human serum), the mean inter-batch CV was 5·2 %. The CV for the NIST standards were well below those of previously published reports(Reference Gallicchio, Helzlsouer and Chow48). A number of studies have demonstrated that 25(OH)D is extremely stable with long-term storage (as long as 40 years) and under a variety of storage conditions(Reference Zerwekh49Reference Bodnar, Catov and Wisner52). Thus, our samples, which were stored for an average of 13 (sd 1·2) years, should have experienced little, if any, 25(OH)D degradation, and any degradation would have been consistent across all samples.

SNP selection

Using a candidate pathway approach, we included common genetic variants (minor allele frequency ≥ 20 %) of genes involved in 25(OH)D synthesis (CYP2R1, CYP3A4, CYP27B1), transport (GC), gene transcription (VDR) and catabolism (CYP24A1). Data from the International HapMap Project (Tagger Pairwise method, HapMap Data Release 27 Phase II+III, February 2009, on NCBI B36 assembly, dbSNP b126 for the Han Chinese in Beijing, China (CHB) population) were used to identify haplotype tagging SNP for the Han Chinese. No HapMap tag SNP were identified for CYP27A1 in the CHB population. Also, any SNP that had been documented in the literature as having functional and/or phenotypic relevance were included.

A total of fifty-five SNP and eight proxies met the inclusion criteria (by rs (reference SNP) number, proxies in parentheses): VDR: 731236 also known as TaqI, 1540339, 1544410 also known as BsmI, 2107301 (12717991), 2189480, 2228570 also known as FokI (and previously reported as rs10735810), 2238136, 2239180, 2239182, 2239184, 2239186, 2254210, 2283342, 2525044, 2853564, 2853559, 3782905, 3847987, 4760658, 7305032, 739837, 757343, 7965281, 7975232 also known as ApaI, 10783215, 10875695 (7136534), 11168268 (11168266), 11168275, 11168287, 11574143, 12721364. CYP2R1: 1993116, 12794714 (10500804), 10741657. CYP3A4: 2242480, 2246709, 3735451. CYP27B1: 4646536. CYP24A1: 912505, 2181874, 2209314, 2296241 (2585428), 2762941, 3787555, 4809958 (6068816), 6022999. GC: 4588 (2298850), 7041, 10488854, 12512631, 16847015, 222016, 2298849 (3733359), 705117, 705120.

DNA extraction and genotype determinations

DNA extraction and genotype determinations were performed by the University of Minnesota's BioMedical Genomics Center. DNA was extracted from buffy coats using a Qiagen Kit (Qiagen, Inc.). Genotype determinations were performed on the commercially available high-throughput genotyping SequenomMassARRAY platform (Sequenom, Inc.). Uniquely located negative controls were routinely included in each plate. These wells were used as controls for genotyping assays, and their unique locations serve as a fingerprint to identify the plate and its orientation. For quality control purposes, only SNP with >95 % call rates were included in the analyses. All SNP were found to be in Hardy–Weinberg equilibrium, except GC rs4588. Genotyping for rs4588 was repeated using a TaqMan assay (Invitrogen). The data from the TaqMan assay were found to be in Hardy–Weinberg equilibrium, and were used for all analyses.

Data analysis

The distribution of serum 25(OH)D in the study population was markedly skewed towards low values. Thus, the statistical analyses were performed on logarithmically transformed values, and geometric means are presented. ANOVA methods were used to compare mean serum concentrations of 25(OH)D by potential lifestyle, sociodemographic, dietary and genetic predictors of 25(OH)D concentration. Dialect group, education level, menopausal status (women), BMI, height, weight, body surface area, physical activity, smoking status, hours spent sitting at work, season of blood draw, use of cod liver oil supplements and dietary intake of vitamin D, Ca, fish, dairy products and alcohol were considered as potential predictors. Age, sex and time interval between last meal and blood draw were included as covariates in all analyses. Dietary vitamin D and hours spent sitting at work were also included as covariates in ANOVA for the assessment of genetic predictors. To test for linear trend, the potential predictor was included as an ordinal variable in general linear models. Each potential predictor was examined individually by assessing its effect on the overall model fit (R 2, F-test). The final multivariate model for non-genetic predictors of serum 25(OH)D concentrations included age, sex, dietary vitamin D intake, hours spent sitting at work, and time interval between last meal and blood draw.

Associations between single SNP and serum 25(OH)D concentrations were evaluated individually in the multivariate-adjusted models. Trends in serum 25(OH)D concentrations across each SNP genotype were tested for statistical significance by including the SNP as a three-level variable (homozygous wildtype, heterozygous and homozygous variant).

All P values quoted are two-tailed, and significance was defined as P< 0·05. For the SNP analyses, statistical significance was defined as P≤ 0·01 to minimise the likelihood of reporting false-positive findings due to multiple comparisons. Calculations were performed using the SAS statistical software system (SAS Institute).

Results

The study population was 56 % women, and the mean age of study participants was 55·7 years (Table 1). Most women were postmenopausal at baseline. Compared to men, women were less educated, less likely to be a smoker, spent less time sitting at work or in vigorous work, and consumed less dietary vitamin D. Mean serum 25(OH)D concentration was 68·6 nmol/l overall, and lower in women (mean: 64·2 nmol/l) than in men (74·3 nmol/l, P< 0·001), and a greater percentage of women (18 %) had 25(OH)D concentrations < 50 nmol/l compared to men (9 %).

Table 1 Characteristics of subjects within the Singapore Chinese Health Study, overall and by sex (Mean values, standard deviations, ranges, number of subjects and percentages)

25(OH)D, 25-hydroxyvitamin D.

* χ2 or t test P-values for differences by sex for categorical and continuous variables, respectively.

Serum 25(OH)D concentrations were statistically significantly associated with dietary vitamin D, Ca and dairy product intake among women, but not men (Table 2). Serum 25(OH)D levels decreased with increasing number of hours spent sitting at work for both men and women, although the linear relationship was not statistically significant for women. There were no associations with 25(OH)D concentrations for cod liver oil supplement use, fish intake or time between last meal and blood draw, regardless of sex (data not shown). Women engaging in vigorous work for at least half an hour per week showed significantly higher serum 25(OH)D level compared to their counterparts with no vigorous work, whereas there was no difference in serum 25(OH)D level between men with and without vigorous work. Age, BMI, alcohol intake and working status were not associated with serum 25(OH)D concentration in this population. When age, sex, dietary vitamin D intake, hours spent sitting at work, and time interval between last meal and blood draw were considered simultaneously in the final multivariate model, 10·2 % of the variation in serum 25(OH)D concentrations was explained (P< 0·01).

Table 2 Geometric means of 25-hydroxyvitamin D (25(OH)D) by potential predictors overall and by sex (Mean values and standard deviations)

* Mean values were adjusted for age (years), sex (among overall) and time interval from last meal to blood draw.

Asian-specific BMI cut-points were used(Reference Choo54).

Median values of dietary vitamin D intake for tertiles (μg/4·184 kJ) were: 0·8, 1·4 and 2·5, respectively, for all subjects; 0·9, 1·4 and 2·4, respectively, among men; and 0·8, 1·6 and 2·8, respectively, among women.

§ Median values of dietary Ca intake for tertiles (mg/4·184 kJ) were: 177, 214 and 362, respectively, for all subjects; 168, 217 and 309, respectively, among men; and 185, 260 and 395, respectively, among women.

Median values of dairy product intake (g/4·184 kJ) were: 1·3, 14·4 and 118, respectively, for all subjects; 1·1, 10·8 and 77·0, respectively, among men; and 1·5, 17·3 and 134, respectively, among women.

Among women, only geometric means for none or any alcohol and beer intake is presented, as only nine women reported drinking any alcohol.

Of the fifty-five SNP assessed, eight SNP in CYP2R1, CYP3A4 and GC were associated with 25(OH)D concentrations (Table 3). For five of the SNP, the major allele was associated with lower 25(OH)D concentrations (e.g. CYP2R1 rs10741657 and rs1993116; and GC rs7041, rs2298849 and rs1687015), while for the remaining three SNP, the major allele was associated with higher 25(OH)D concentrations (e.g. CYP2R1 rs12794714; CYP3A4 rs2242480; and GC rs4588). The strongest association was with the GC SNP rs4588, where the decrease in copies of the major allele (e.g. from 2 to 0) was associated with a 11·5 nmol/l decrease in serum 25(OH)D concentration (P< 0·001). Including genotype information for individual SNP into the multivariable model explained an additional 0·8–3·7 % of the variation in serum 25(OH)D concentrations in this cohort (Table 3).

Table 3 Geometric means of serum 25-hydroxyvitamin D (25(OH)D) by genotype (Geometric mean values and 95 % confidence intervals)

CYP2R1, cytochrome P450 2R1; rs, refSNP or reference SNP; CYP3A4, cytochrome P450 3A4, GC, group-specific component (vitamin D binding protein).

* 25(OH)D concentrations are multivariate-adjusted for: age, sex, dietary vitamin D, hours spent sitting at work, and time interval between last meal and blood draw.

Variance explained is the model R 2× 100. The partial R 2 for each of the covariates in the individual multivariate models was: age ( < 0·001), sex (0·070–0·076), dietary vitamin D (0·005–0·008), hours spent sitting at work (0·020–0·026), and time interval between last meal and blood draw ( < 0·001).

The variance explained by genotype is the partial R 2× 100 for the individual genotype in a linear regression model with variables for age, sex, dietary vitamin D, hours spent sitting at work, and time interval between last meal and blood draw.

In addition to the single-SNP associations with serum 25(OH)D concentrations, we also evaluated the combined effect of two well-described GC SNP, rs4588 and rs7041. These SNP have been previously shown to jointly determine three well-described protein transcripts: Gc1s, Gc1f and Gc2(Reference Speeckaert, Huang and Delanghe12). As shown in Table 4, the mean 25(OH)D concentration was the highest for individuals with two copies of the Gc1s allele (Gc1s-1s: rs4588CC, rs7041GG), the lowest for individuals with two copies of the Gc2 allele (Gc2-2: rs4588AA, rs7041TT), and intermediate for those with any other Gc haplotype (P-trend < 0·001). When stratified by median dietary vitamin D intake, the trends with 25(OH)D by haplotype are similar to the overall pattern currently presented in Table 4 for both high and low dietary vitamin D intake (data not shown). However, when stratified by hours spent sitting at work, the trend of decreasing 25(OH)D concentration from Gc1s-1s to Gc2-2 haplotype was only evident among those who reported no hours spent sitting at work (n 268, P-trend < 0·001, Table 5). The trend was less evident and was not statistically significant among those who reported any sitting hours at work (data not shown), although the interaction between Gc haplotype and hours spent sitting at work was not statistically significant (P-interaction = 0·24).

Table 4 Serum 25-hydroxyvitamin D (nmol/l) by group-specific component (GC) haplotype (n 467) (Adjusted geometric means, 95 % confidence intervals and number of subjects)

* Adjusted geometric mean value was significantly different from that of the Gc2-2 haplotype (P< 0·01).

From ANCOVA with adjustments for age, sex, dietary vitamin D intake, hours spent sitting at work and time interval between the last meal and blood draw.

Table 5 Serum 25-hydroxyvitamin D (nmol/l) by group-specific component (GC) haplotype for subjects who reported no hours sitting at work (n 268) (Adjusted geometric means, 95 % confidence intervals and number of subjects)

From ANCOVA with adjustments for age, sex, dietary vitamin D intake, hours spent sitting at work and time interval between the last meal and blood draw.

Discussion

On average, participants in this study had serum 25(OH)D concentrations that would be considered sufficient for optimal bone health according to the Institute of Medicine recommendations(Reference Ross, Taylor, Yaktine and Del Valle6), but somewhat below the recommendations from leading vitamin D researchers of 75 nmol/l to address a broader range of health concerns(Reference Vieth, Bischoff-Ferrari and Boucher7). However, the serum 25(OH)D concentrations for participants in this study were higher, on average, than has been reported for comparably aged men (58·0 nmol/l) and women (54·8 nmol/l) in the USA(Reference Looker, Johnson and Lacher53), which might be expected of individuals living near the equator.

BMI was not associated with serum 25(OH)D concentrations in this population, mostly probably due to the fact that the mean BMI for the study participants was within the normal range for Asian populations ( < 23 kg/m2)(Reference Choo54), and the range was narrow. Mean dietary vitamin D intake for these middle- and older-aged Singaporean men (2·8 (sd 1·8) μg/d) and women (2·4 (sd 1·7) μg/d) was lower than comparably aged men (5·1 (sd 0·3) μg/d) and women (3·9 (sd 0·4) μg/d) in the USA according to the 2005–6 National Health and Nutrition Examination Survey data(Reference Bailey, Dodd and Goldman55), and somewhat lower than adult men (3·1 (sd 0·1) μg/d) and women (2·7 (sd 0·1) μg/d) in the UK according to the 2008–9 National Diet and Nutrition Survey data(Reference Bates, Lennox and Swan56).

Sun exposure is known to be a major determinant of vitamin D status. Singapore receives 12 h sunlight/d throughout the year, with a midday solar zenith angle that ranges from a minima of 0–3° (March, September) to a peak of 22–25° (June, December)(57). Despite a small amount of variation in solar zenith angle, we did not observe significant seasonal variation in serum 25(OH)D concentrations in this cohort. The UV index ranges from 10 (December) to 13 (March/April), indicating very high ambient UV radiation levels(58). Given the average daily high temperatures of 31°C (88°F)(59), many Singaporeans avoid the heat of the midday sun. In our study, as the reported number of hours spent sitting at work increased, serum 25(OH)D concentrations decreased (P-trend < 0·001).

Three other reports have described serum 25(OH)D concentrations among healthy adults living within 10° of the equator. Rahman et al. (Reference Rahman, Chee and Yassin60) evaluated 276 post-menopausal women living near Kuala Lumpur, Malaysia (2°N). They found that ethnic Chinese women had significantly higher mean 25(OH)D concentrations compared to Malay women (68·8 (sd 15·7) v. 44·4 (10·6) nmol/l, P< 0·05). Dietary vitamin D intake did not differ between the two ethnic groups. The Chinese women had a significantly lower mean BMI, and reported more regular physical activity than the Malay women. The Malay also tend to have more skin pigmentation than the Chinese, and many Malay women follow Muslim dress codes that further limit UV exposure. Serum 25(OH)D concentrations were significantly correlated with BMI, fat mass, parathyroid hormone concentrations and physical activity scores. Green et al. (Reference Green, Skeaff and Rockell61) evaluated 378 younger women living in Kuala Lumpur (mean age: 25·2 years) and 126 women living in Jakarta, Indonesia (6°S, mean age: 30·0 years). Among the women in Malaysia, they also found higher serum 25(OH)D concentrations among the ethnic Chinese (mean: 58·0 nmol/l, 95 % CI 55·0, 61·0) compared to the Malay (mean: 43·0 nmol/l, 95 % CI 40·0, 46·0) or Indian (mean: 45·0 nmol/l, 95 % CI 43·0, 48·0) women (P< 0·01). The Indonesian women had serum 25(OH)D concentrations that were comparable to the Malay and Indian women (mean: 46·0 nmol/l, 95 % CI 43·0, 48·0). Moy & Bulgiba(Reference Moy and Bulgiba62) recently reported on the vitamin D status of 380 Malay study participants (158 men, 222 women) in a voluntary health screening programme in Kuala Lumpur. The women had significantly lower serum 25(OH)D concentrations (mean: 36·2 nmol/l, 95 % CI 34·5, 38·0) compared to the men (mean: 56·2 nmol/l, 95 % CI 53·2, 59·2, P< 0·001), which could be partially explained by differences in religious dress codes. Age, sex, BMI and abdominal obesity were found to be statistically significantly associated with vitamin D insufficiency in this study cohort. The women in our study had serum 25(OH)D concentrations that were similar to the ethnic Chinese participants in both the Rahman et al.(Reference Rahman, Chee and Yassin60) and Green et al. (Reference Green, Skeaff and Rockell61) studies.

Genetic variants in CYP2R1, CYP3A4 and GC were significantly associated with serum 25(OH)D concentrations in our study. The CYP2R1 rs10741657, rs12794714, and rs1993116 and GC rs4588 and rs7041 findings are consistent with several recent reports(Reference Engelman, Fingerlin and Langefeld19, Reference Ramos-Lopez, Bruck and Jansen21, Reference Bu, Armas and Lappe26, Reference Kurylowicz, Ramos-Lopez and Bednarczuk28Reference Fu, Yun and Oczak32, Reference Fang, van Meurs and Arp34, Reference Janssens, Bouillon and Claes35), including two large genome-wide association studies(Reference Wang, Zhang and Richards24, Reference Ahn, Yu and Stolzenberg-Solomon25). Several CYP enzymes have been shown to have 25-hydroxylase activity, and CYP2R1 has emerged as the predominant 25-hydroxylase, with the highest binding affinity and specificity for vitamin D(Reference Schuster63). CYP2R1 rs10741657 lies in the promoter region, rs12794714 is a synonymous SNP in exon 1, and rs1993116 is in intron 1. To our knowledge, no previous studies have evaluated the association between genetic variation in CYP3A4 and serum 25(OH)D concentrations in human subjects. CYP3A4 rs2242480 is in intron 10 near the exon/intron boundary. Although the functional relevance of this SNP is unclear, a recent pharmacokinetic study suggests that individuals with the homozygous variant rs2242480TT genotype have significantly lower CYP3A4 activity(Reference Zhang, Chang and Kan64).

The GC SNP, rs4588 and rs7041, are both in exon 11. Consistent with the findings of several previous studies(Reference Engelman, Fingerlin and Langefeld19, Reference Kurylowicz, Ramos-Lopez and Bednarczuk28, Reference Daiger, Miller and Chakraborty65Reference Constans, Arlet and Viau68), we also observed that individuals with two copies of the Gc2 allele (Gc2-2) have significantly lower 25(OH)D concentrations compared to other GC genotypes. In vitro data have shown that the Gc2 protein has a significantly lower affinity constant for 25(OH)D3 compared to the Gc1s or Gc1f proteins(Reference Speeckaert, Huang and Delanghe12). GC rs16847015 and rs2298849 both lie in intron 1, and their functional relevance is unclear.

While several previous reports identified genetic variants in the VDR (Reference Ma, Stampfer and Gann13Reference Smolders, Damoiseaux and Menheere20) and CYP27B1 (Reference Orton, Morris and Herrera22, Reference Ramos-Lopez, Kahles and Weber23, Reference Hypponen, Berry and Wjst27) genes as being associated with serum 25(OH)D concentrations, we did not observe any statistically significant associations for those genes in our study population. The studies reporting significant findings for VDR and CYP27B1 variants tended to be among smaller study populations, conducted at higher latitudes, and many failed to adjust for factors known to alter serum 25(OH)D concentrations such as season of blood draw, BMI and dietary/supplemental vitamin D intake.

The present study has several strengths including being the largest study of vitamin D status among Southeast Asians to date, lack of seasonal UV variation, and consideration of dietary vitamin D exposures, lifestyle and sociodemographic factors as well as genetic variation as potential factors contributing to serum 25(OH)D concentrations. We also took a comprehensive approach to assessing the effect of genetic variation in the entire vitamin D metabolism pathway, as opposed to the evaluation of single candidate genes or SNP. Limitations of the present study include incomplete assessment of time spent outdoors during daylight hours and degree of skin pigmentation, factors which may contribute to variation in serum 25(OH)D concentrations. Due to technical challenges in accurately measuring the vitamin D content of foods(Reference Byrdwell, Devries and Exler37, Reference Holden, Lemar and Exler38), our assessment of the participants' dietary vitamin D intake is probably underestimated. This underestimation of dietary vitamin D intake should occur to all study subjects non-differentially, which would result in underestimating the association between dietary vitamin D intake and serum 25(OH)D concentrations. Supplemental vitamin D intake, other than cod liver oil, was not specifically assessed. However, only 8 % of our study cohort reported taking any vitamins or minerals at least once a week.

Our findings confirm the growing body of literature documenting an association between GC and CYP2R1 genetic variation and serum 25(OH)D concentrations. Future studies of both predictors of 25(OH)D concentrations and disease outcomes thought to be associated with vitamin D status should include an assessment of GC and CYP2R1 genotype. Further research is needed to confirm our findings related to CYP3A4 rs2242480.

Acknowledgements

The authors' contributions to the present study were as follows: K. R., L. M. B. and J.-M. Y. were responsible for the study concept and design and obtainment of funding; K. R. was responsible for study supervision; K. R., K. B. B., D. W., W.-P. K. and J.-M. Y. for the acquisition of data; L. M. B., R. W., K. R. and J.-M. Y. for the data analysis and interpretation; and K. R. and L. M. B. for drafting of the manuscript. All authors contributed to critical revision of the manuscript. None of the authors reported a conflict of interest. The authors thank Siew-Hong Low of the National University of Singapore for supervising the fieldwork of the SCHS, and Kazuko Arakawa for the development and maintenance of the cohort study database. Further, the authors acknowledge the founding, long-standing Principal Investigator of the SCHS – Mimi C. Yu. This study was funded by the University of Minnesota Masonic Cancer Center and the National Cancer Institute (R01 CA144034).

References

1Cohen-Lahav, M, Douvdevani, A, Chaimovitz, C, et al. (2007) The anti-inflammatory activity of 1,25-dihydroxyvitamin D3 in macrophages. J Steroid Biochem Mol Biol 103, 558562.Google Scholar
2Boucher, BJ (2011) Vitamin D insufficiency and diabetes risks. Curr Drug Targets 12, 6187.CrossRefGoogle ScholarPubMed
3World Health Organization, International Agency for Research on Cancer (IARC) (2008) Vitamin D and Cancer. Working Group Report 5. Lyon: WHO.Google Scholar
4Baeke, F, Takiishi, T, Korf, H, et al. (2010) Vitamin D: modulator of the immune system. Curr Opin Pharmacol 10, 482496.CrossRefGoogle ScholarPubMed
5Reddy Vanga, S, Good, M, Howard, PA, et al. (2010) Role of vitamin D in cardiovascular health. Am J Cardiol 106, 798805.Google Scholar
6Institute of Medicine (2011) Dietary Reference Intakes for Vitamin D and Calcium [Ross, AC, Taylor, CL, Yaktine, AL and Del Valle, HB, editors]. Washington, DC: Institute of Medicine.Google Scholar
7Vieth, R, Bischoff-Ferrari, H, Boucher, BJ, et al. (2007) The urgent need to recommend an intake of vitamin D that is effective. Am J Clin Nutr 85, 649650.Google Scholar
8Heaney, RP & Holick, MF (2011) Why the IOM recommendations for vitamin D are deficient. J Bone Miner Res 26, 455457.Google Scholar
9Norman, AW (2011) Vitamin D nutrition is at a crossroads. Public Health Nutr 14, 744745.CrossRefGoogle Scholar
10Schwalfenberg, GK & Whiting, SJ (2011) A Canadian response to the 2010 Institute of Medicine vitamin D and calcium guidelines. Public Health Nutr 14, 746748.Google Scholar
11Holick, MF (2007) Vitamin D deficiency. N Engl J Med 357, 266281.Google Scholar
12Speeckaert, M, Huang, G, Delanghe, JR, et al. (2006) Biological and clinical aspects of the vitamin D binding protein (Gc-globulin) and its polymorphism. Clin Chim Acta 372, 3342.Google Scholar
13Ma, J, Stampfer, MJ, Gann, PH, et al. (1998) Vitamin D receptor polymorphisms, circulating vitamin D metabolites, and risk of prostate cancer in United States physicians. Cancer Epidemiol Biomarkers Prev 7, 385390.Google Scholar
14Ogunkolade, BW, Boucher, BJ, Prahl, JM, et al. (2002) 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 51, 22942300.Google Scholar
15d'Alesio, A, Garabedian, M, Sabatier, JP, et al. (2005) Two single-nucleotide polymorphisms in the human vitamin D receptor promoter change protein–DNA complex formation and are associated with height and vitamin D status in adolescent girls. Hum Mol Genet 14, 35393548.Google Scholar
16Zajickova, K, Hill, M, Vankova, M, et al. (2006) Low-density lipoprotein receptor-related protein 5 and vitamin D receptor gene polymorphisms in relation to vitamin D levels in menopause. Clin Chem Lab Med 44, 10661069.Google Scholar
17Vupputuri, MR, Goswami, R, Gupta, N, et al. (2006) Prevalence and functional significance of 25-hydroxyvitamin D deficiency and vitamin D receptor gene polymorphisms in Asian Indians. Am J Clin Nutr 83, 14111419.Google Scholar
18Baroncelli, GI, Bereket, A, El Kholy, M, et al. (2008) Rickets in the Middle East: role of environment and genetic predisposition. J Clin Endocrinol Metab 93, 17431750.Google Scholar
19Engelman, CD, Fingerlin, TE, Langefeld, CD, et al. (2008) Genetic and environmental determinants of 25-hydroxyvitamin D and 1,25-dihydroxyvitamin D levels in Hispanic and African Americans. J Clin Endocrinol Metab 93, 33813388.Google Scholar
20Smolders, J, Damoiseaux, J, Menheere, P, et al. (2009) Fok-I vitamin D receptor gene polymorphism (rs10735810) and vitamin D metabolism in multiple sclerosis. J Neuroimmunol 207, 117121.Google Scholar
21Ramos-Lopez, E, Bruck, P, Jansen, T, et al. (2007) CYP2R1 (vitamin D 25-hydroxylase) gene is associated with susceptibility to type 1 diabetes and vitamin D levels in Germans. Diabetes Metab Res Rev 23, 631636.Google Scholar
22Orton, SM, Morris, AP, Herrera, BM, et al. (2008) Evidence for genetic regulation of vitamin D status in twins with multiple sclerosis. Am J Clin Nutr 88, 441447.CrossRefGoogle ScholarPubMed
23Ramos-Lopez, E, Kahles, H, Weber, S, et al. (2008) Gestational diabetes mellitus and vitamin D deficiency: genetic contribution of CYP27B1 and CYP2R1 polymorphisms. Diabetes Obes Metab 10, 683685.Google Scholar
24Wang, TJ, Zhang, F, Richards, JB, et al. (2010) Common genetic determinants of vitamin D insufficiency: a genome-wide association study. Lancet 376, 180188.Google Scholar
25Ahn, J, Yu, K, Stolzenberg-Solomon, R, et al. (2010) Genome-wide association study of circulating vitamin D levels. Hum Mol Genet 19, 27392745.CrossRefGoogle ScholarPubMed
26Bu, FX, Armas, L, Lappe, J, et al. (2010) Comprehensive association analysis of nine candidate genes with serum 25-hydroxy vitamin D levels among healthy Caucasian subjects. Hum Genet 128, 549556.Google Scholar
27Hypponen, E, Berry, DJ, Wjst, M, et al. (2009) Serum 25-hydroxyvitamin D and IgE – a significant but nonlinear relationship. Allergy 64, 613620.Google Scholar
28Kurylowicz, A, Ramos-Lopez, E, Bednarczuk, T, et al. (2006) Vitamin D-binding protein (DBP) gene polymorphism is associated with Graves’ disease and the vitamin D status in a Polish population study. Exp Clin Endocrinol Diabetes 114, 329335.Google Scholar
29Lauridsen, AL, Vestergaard, P, Hermann, AP, et al. (2005) 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 77, 1522.Google Scholar
30Abbas, S, Linseisen, J, Slanger, T, et al. (2008) The Gc2 allele of the vitamin D binding protein is associated with a decreased postmenopausal breast cancer risk, independent of the vitamin D status. Cancer Epidemiol Biomarkers Prev 17, 13391343.Google Scholar
31Sinotte, M, Diorio, C, Berube, S, et al. (2009) Genetic polymorphisms of the vitamin D binding protein and plasma concentrations of 25-hydroxyvitamin D in premenopausal women. Am J Clin Nutr 89, 634640.Google Scholar
32Fu, L, Yun, F, Oczak, M, et al. (2009) Common genetic variants of the vitamin D binding protein (DBP) predict differences in response of serum 25-hydroxyvitamin D [25(OH)D] to vitamin D supplementation. Clin Biochem 42, 11741177.Google Scholar
33Ahn, J, Albanes, D, Berndt, SI, et al. (2009) Vitamin D-related genes, serum vitamin D concentrations and prostate cancer risk. Carcinogenesis 30, 769776.Google Scholar
34Fang, Y, van Meurs, JB, Arp, P, et al. (2009) Vitamin D binding protein genotype and osteoporosis. Calcif Tissue Int 85, 8593.Google Scholar
35Janssens, W, Bouillon, R, Claes, B, et al. (2010) Vitamin D deficiency is highly prevalent in COPD and correlates with variants in the vitamin D-binding gene. Thorax 65, 215220.Google Scholar
36Hankin, JH, Stram, DO, Arakawa, K, et al. (2001) Singapore Chinese Health Study: development, validation, and calibration of the quantitative food frequency questionnaire. Nutr Cancer 39, 187195.Google Scholar
37Byrdwell, WC, Devries, J, Exler, J, et al. (2008) Analyzing vitamin D in foods and supplements: methodologic challenges. Am J Clin Nutr 88, 554S557S.Google Scholar
38Holden, JM, Lemar, LE & Exler, J (2008) Vitamin D in foods: development of the US Department of Agriculture database. Am J Clin Nutr 87, 1092S1096S.CrossRefGoogle ScholarPubMed
39Koh, WP, Yuan, JM, Sun, CL, et al. (2003) Angiotensin I-converting enzyme (ACE) gene polymorphism and breast cancer risk among Chinese women in Singapore. Cancer Res 63, 573578.Google Scholar
40Ersfeld, DL, Rao, DS, Body, JJ, et al. (2004) Analytical and clinical validation of the 25 OH vitamin D assay for the LIAISON automated analyzer. Clin Biochem 37, 867874.Google Scholar
41Carter, GD, Carter, CR, Gunter, E, et al. (2004) Measurement of Vitamin D metabolites: an international perspective on methodology and clinical interpretation. J Steroid Biochem Mol Biol 89–90, 467471.Google Scholar
42Tai, SS, Bedner, M & Phinney, KW (2010) Development of a candidate reference measurement procedure for the determination of 25-hydroxyvitamin D3 and 25-hydroxyvitamin D2 in human serum using isotope-dilution liquid chromatography-tandem mass spectrometry. Anal Chem 82, 19421948.Google Scholar
43Giovannucci, E, Liu, Y, Rimm, EB, et al. (2006) Prospective study of predictors of vitamin D status and cancer incidence and mortality in men. J Natl Cancer Inst 98, 451459.Google Scholar
44Bertone-Johnson, ER, Chen, WY, Holick, MF, et al. (2005) Plasma 25-hydroxyvitamin D and 1,25-dihydroxyvitamin D and risk of breast cancer. Cancer Epidemiol Biomarkers Prev 14, 19911997.Google Scholar
45Freedman, DM, Chang, SC, Falk, RT, et al. (2008) Serum levels of vitamin D metabolites and breast cancer risk in the prostate, lung, colorectal, and ovarian cancer screening trial. Cancer Epidemiol Biomarkers Prev 17, 889894.Google Scholar
46Abbas, S, Linseisen, J, Slanger, T, et al. (2008) Serum 25-hydroxyvitamin D and risk of post-menopausal breast cancer-results of a large case–control study. Carcinogenesis 29, 9399.Google Scholar
47Cauley, JA, Lacroix, AZ, Wu, L, et al. (2008) Serum 25-hydroxyvitamin D concentrations and risk for hip fractures. Ann Intern Med 149, 242250.Google Scholar
48Gallicchio, L, Helzlsouer, KJ, Chow, WH, et al. (2010) Circulating 25-hydroxyvitamin D and the risk of rarer cancers: design and methods of the Cohort Consortium Vitamin D Pooling Project of Rarer Cancers. Am J Epidemiol 172, 1020.Google Scholar
49Zerwekh, JE (2004) The measurement of vitamin D: analytical aspects. Ann Clin Biochem 41, 272281.CrossRefGoogle ScholarPubMed
50Hollis, BW (2008) Measuring 25-hydroxyvitamin D in a clinical environment: challenges and needs. Am J Clin Nutr 88, 507S510S.CrossRefGoogle Scholar
51Antoniucci, DM, Black, DM & Sellmeyer, DE (2005) Serum 25-hydroxyvitamin D is unaffected by multiple freeze-thaw cycles. Clin Chem 51, 258261.Google Scholar
52Bodnar, LM, Catov, JM, Wisner, KL, et al. (2009) Racial and seasonal differences in 25-hydroxyvitamin D detected in maternal sera frozen for over 40 years. Br J Nutr 101, 278284.Google Scholar
53Looker, AC, Johnson, CL, Lacher, DA, et al. (2011) Vitamin D Status: United States, 2001–2006. NCHS Data Brief, No. 59. Hyattsville, MD: National Center for Health Statistics.Google Scholar
54Choo, V (2002) WHO reassesses appropriate body-mass index for Asian populations. Lancet 360, 235.Google Scholar
55Bailey, RL, Dodd, KW, Goldman, JA, et al. (2010) Estimation of total usual calcium and vitamin D intakes in the United States. J Nutr 140, 817822.Google Scholar
56Bates, B, Lennox, A & Swan, G (editors) (2010) National Diet Nutrition Survey (2008/2009). London: Food Standards Agency, Department of Health.Google Scholar
57 Cornwall C, Horiuchi A & Lehman C (2011) NOAA Solar Position Calculator. National Oceanic and Atmospheric Administration Earth System Research Lab. http://www.srrb.noaa.gov/highlights/sunrise/azel.html (accessed 23 April 2011).Google Scholar
58 World Health Organization (2011) Ultraviolet Radiation and the INTERSUN Programme: UV Index. http://www.who.int/uv/intersunprogramme/activities/uv_index/en/index3.html (accessed 23 April 2011).Google Scholar
59Department of Statistics, Ministry of Trade and Industry, Republic of Singapore (2011) Monthly Digest of Statistics Singapore, March 2011. Republic of Singapore: Ministry of Trade and Industry.Google Scholar
60Rahman, SA, Chee, WS, Yassin, Z, et al. (2004) Vitamin D status among postmenopausal Malaysian women. Asia Pac J Clin Nutr 13, 255260.Google Scholar
61Green, TJ, Skeaff, CM, Rockell, JE, et al. (2008) Vitamin D status and its association with parathyroid hormone concentrations in women of child-bearing age living in Jakarta and Kuala Lumpur. Eur J Clin Nutr 62, 373378.Google Scholar
62Moy, FM & Bulgiba, A (2011) High prevalence of vitamin D insufficiency and its association with obesity and metabolic syndrome among Malay adults in Kuala Lumpur, Malaysia. BMC Public Health 11, 735.Google Scholar
63Schuster, I (2011) Cytochromes P450 are essential players in the vitamin D signaling system. Biochim Biophys Acta 1814, 186199.CrossRefGoogle ScholarPubMed
64Zhang, W, Chang, YZ, Kan, QC, et al. (2010) CYP3A4*1G genetic polymorphism influences CYP3A activity and response to fentanyl in Chinese gynecologic patients. Eur J Clin Pharmacol 66, 6166.CrossRefGoogle ScholarPubMed
65Daiger, SP, Miller, M & Chakraborty, R (1984) Heritability of quantitative variation at the group-specific component (Gc) locus. Am J Hum Genet 36, 663676.Google Scholar
66Lauridsen, AL, Vestergaard, P & Nexo, E (2001) Mean serum concentration of vitamin D-binding protein (Gc globulin) is related to the Gc phenotype in women. Clin Chem 47, 753756.Google Scholar
67Brown, IR, Carter, ND & Sood, A (1979) Vitamin D binding globulin phenotypes in liver disease. Clin Chim Acta 95, 7582.Google Scholar
68Constans, J, Arlet, P, Viau, M, et al. (1983) Unusual sialilation of the serum DBP associated with the Gc 1 allele in alcoholic cirrhosis of the liver. Clin Chim Acta 130, 219230.Google Scholar
Figure 0

Table 1 Characteristics of subjects within the Singapore Chinese Health Study, overall and by sex (Mean values, standard deviations, ranges, number of subjects and percentages)

Figure 1

Table 2 Geometric means of 25-hydroxyvitamin D (25(OH)D) by potential predictors overall and by sex (Mean values and standard deviations)

Figure 2

Table 3 Geometric means of serum 25-hydroxyvitamin D (25(OH)D) by genotype (Geometric mean values and 95 % confidence intervals)

Figure 3

Table 4 Serum 25-hydroxyvitamin D (nmol/l) by group-specific component (GC) haplotype (n 467)† (Adjusted geometric means, 95 % confidence intervals and number of subjects)

Figure 4

Table 5 Serum 25-hydroxyvitamin D (nmol/l) by group-specific component (GC) haplotype for subjects who reported no hours sitting at work† (n 268) (Adjusted geometric means, 95 % confidence intervals and number of subjects)