Next Article in Journal
Dietary Habits and Nutrition in Rheumatoid Arthritis: Can Diet Influence Disease Development and Clinical Manifestations?
Next Article in Special Issue
Interaction between Metabolic Genetic Risk Score and Dietary Fatty Acid Intake on Central Obesity in a Ghanaian Population
Previous Article in Journal
The Effects on Immune Function and Digestive Health of Consuming the Skin and Flesh of Zespri® SunGold Kiwifruit (Actinidia Chinensis var. Chinensis ‘Zesy002’) in Healthy and IBS-Constipated Individuals
Previous Article in Special Issue
A Nutrigenetic Approach to Investigate the Relationship between Metabolic Traits and Vitamin D Status in an Asian Indian Population
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Nutrients
Volume 12
Issue 5
10.3390/nu12051455
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

Independent and Interactive Influences of Environmental UVR, Vitamin D Levels, and Folate Variant MTHFD1-rs2236225 on Homocysteine Levels

by
Patrice Jones
1,2,*,
Mark Lucock
1,
Charlotte Martin
1,
Rohith Thota
3,4,
Manohar Garg
3,
Zoe Yates
5,
Christopher J. Scarlett
1,
Martin Veysey
6 and
Emma Beckett
1,2
1
School of Environmental & Life Sciences, University of Newcastle, Ourimbah, NSW 2258, Australia
2
Hunter Medical Research Institute, New Lambton Heights, NSW 2305, Australia
3
Nutraceuticals Research Group, University of Newcastle, Callaghan, NSW 2308, Australia
4
Riddet Institute, Massey University, Palmerston North 4442, New Zealand
5
Biomedical Sciences & Pharmacy, University of Newcastle, Callaghan, NSW 2308, Australia
6
Hull-York Medical School, University of Hull, Hull YO10 5DD, UK
*
Author to whom correspondence should be addressed.
Nutrients 2020, 12(5), 1455; https://0-doi-org.brum.beds.ac.uk/10.3390/nu12051455
Submission received: 15 April 2020 / Revised: 11 May 2020 / Accepted: 14 May 2020 / Published: 18 May 2020
(This article belongs to the Special Issue Nutrient Gene Interactions)
Browse Figure
Versions Notes

Abstract

:
Elevated homocysteine (Hcy) levels are a risk factor for vascular diseases. Recently, increases in ultraviolet radiation (UVR) have been linked to decreased Hcy levels. This relationship may be mediated by the status of UVR-responsive vitamins, vitamin D and folate, and/or genetic variants influencing their levels; however, this has yet to be examined. Therefore, the independent and interactive influences of environmental UVR, vitamin D and folate levels and related genetic variants on Hcy levels were examined in an elderly Australian cohort (n = 619). Red blood cell folate, 25-hydroxyvitamin D (25(OH)D), and plasma Hcy levels were determined, and genotyping for 21 folate and vitamin D-related variants was performed. Erythemal dose rate accumulated over six-weeks (6W-EDR) and four-months (4M-EDR) prior to clinics were calculated as a measure of environmental UVR. Multivariate analyses found interactions between 6W-EDR and 25(OH)D levels (pinteraction = 0.002), and 4M-EDR and MTHFD1-rs2236225 (pinteraction = 0.006) in predicting Hcy levels. The association between 6W-EDR and Hcy levels was found only in subjects within lower 25(OH)D quartiles (<33.26 ng/mL), with the association between 4M-EDR and Hcy occurring only in subjects carrying the MTHFD1-rs2236225 variant. 4M-EDR, 6W-EDR, and MTHFD1-rs2236225 were also independent predictors of Hcy. Findings highlight nutrient–environment and gene–environment interactions that could influence the risk of Hcy-related outcomes.

1. Introduction

Homocysteine (Hcy) is a non-essential amino acid produced as an intermediate product in the synthesis of methionine and cysteine. Elevated levels of Hcy (i.e., hyperhomocysteinemia) can lead to the onset of multiple cardiovascular and neurovascular outcomes, such as atherosclerosis, stroke, and Alzheimer’s disease [1,2]. Several genetic, dietary, and other lifestyle factors are established determinants of Hcy status. The effects of these determinants are commonly examined in isolation; however, interactions may exist between determinants that lead to variation in Hcy levels and related risk of disease, and this is understudied. Investigation into such environment–gene–nutrient interactions in Hcy metabolism is limited but may offer new avenues in modifying Hcy and its related risk of pathology.
An inverse relationship between folate and Hcy is well established due to the central role of folate in converting Hcy to methionine [3,4]. Several environmental and genetic factors that regulate systemic folate status, such as dietary intake of folate, or the MTHFR-rs1801133 (commonly MTHFR-C677T) variant are recognised as key modulators of Hcy status. An emerging determinant of blood folate status—ultraviolet radiation (UVR)—is also a potential modulator of Hcy. UVR exposure has been shown to cause the breakdown of folates in vitro [5,6], and increases in UVR exposure have been associated with a decreased folate status in several human cohort studies [6,7,8,9,10]. More recently, we have reported that increases in environmental UVR levels are also associated with decreases in Hcy levels [11]. This association is counter-intuitive given the well-known inverse relationship between folate and Hcy levels, but it may be modulated in part by folate-related genetic factors, which require further investigation. Previous investigation into the negative association between UVR levels and red blood cell folate levels have found that this association was dependent on the presence of the MTHFR-rs1801133 variant. Furthermore, the distribution of multiple folate-related polymorphisms are linked to latitude [12] and skin phototype [13], indicating multiple folate polymorphisms that may interact with UVR exposure and potentially modulate the relationship between UVR and Hcy levels.
The relationship between Hcy and UVR may in part be dictated by changes in vitamin D levels. The relationship between UVR environment and vitamin D levels is well known, with vitamin D levels also having been reported to be inversely associated with Hcy. An inverse association between blood vitamin D and Hcy was first shown in data from the National Health and Nutrition Examination Survey (2001–2006) [14] and has since been replicated in several observational and intervention studies examining both healthy subjects [15,16] and subjects with Hcy-related conditions [17,18]. The mechanism by which vitamin D is linked to Hcy has not been elucidated, but likely relates to the influences of vitamin D on cystathionine B-synthase activity, a key enzyme in the folate-independent transsulfuration pathway that interconverts cysteine and Hcy [19]. However this relationship has only been demonstrated in vitro [19]. Hcy may be influenced by changes in vitamin D levels, as well as variants in related genes that influence overall vitamin D function. Multiple genes control vitamin D synthesis, metabolism and activity, with the frequency of multiple common variants in these genes known to vary by UVR environment and latitude [20,21,22,23,24,25].
Therefore, the aim of the present study was to further examine the negative association between environmental UVR and Hcy status reported in an elderly Australian cohort [26], examining whether this association is modulated by independent or interactive influences of vitamin D and/or folate levels, and related genetic factors.

2. Materials and Methods

2.1. Subjects

This study was a secondary analysis using samples and pre-existing data from the Retirement Health and Lifestyle Study (RHLS), a cross-sectional study examining the health and lifestyle of older Australians (>65 years) living in the Central Coast region of NSW, Australia (n = 650) [11,26,27,28]. This Australian elderly cohort was an appropriate focus for this investigation and previous studies due to Australians, in general, being exposed to high levels of environmental UVR [29], with elderly Australians a sub-population particularly at risk of adverse UVR-related effects as biological photoprotective factors reduce with age [30]. Briefly, subjects were eligible to be included in the initial RHLS if they were ≥65 years and lived independently within the community or resided in retirement villages on the Central Coast, NSW, Australia for at least the last 12 months. Subjects were ineligible if they did not live independently, if another household member was already taking part in the study, if they had language/communicative difficulties and/or were unable to provide informed consent. Subjects were eligible to be included in this secondary analysis if data was available for plasma Hcy levels (n = 619). Subjects provided written informed consent, and ethics approval for the study was obtained from the University of Newcastle Human Research Ethics Committee (reference no. H-2008-0431). Further details on the original study population and design have been reported previously [11,26,27,28].

2.2. Blood Biochemistry; Hcy and Vitamin Levels

Fasting blood samples of subjects were collected in EDTA-lined tubes. Blood samples were stored at −20 °C as whole blood, or stored at −80 °C as plasma or serum components in lithium heparin tubes or tubes containing clot activator respectively. Red blood cell (RBC) folate, serum vitamin B12 and creatinine were assessed by the Hunter Area Pathology Service via standardised assays. Total plasma Hcy was measured by a selective fluorescence assay (JD Biotech Corp, Taipei, Taiwan; linear range 1–60 umol/L, CV < 8%) [31,32]. All samples were tested in triplicate, and control samples were tested between batches as per manufacturer’s protocol (JD Biotech Corp, Taipei, Taiwan). Vitamin D levels were assessed via HPLC methods as the total serum concentration of 25-hydroxyvitamin D (25(OH)D) [33]. The detection range for this method was 9–193 ng/mL with a CV% < 10%. Method precision and accuracy were examined for quality control via repetitive assessments of triplicates of quality control samples prepared with different concentrations of vitamin D (20 nmol/L and 40 nmol/L 25(OH)D).

2.3. Dietary Intake, Smoking Status, and Body Mass Index

Dietary intake of subjects was estimated using a validated self-administered food frequency questionnaire (FFQ) covering 225 food items and all food groups, with questionnaires analysed using FoodworksTM software (V.2.10.146; Xyris Software, Brisbane, QLD, Australia). These food frequency questionnaires were as published in previous studies [28,31] and were adapted from the Commonwealth Scientific and Industrial Research Organisation (CSIRO) FFQ [34]. Additional intake of nutrients via supplement use was included in calculated nutrient intakes as total dietary equivalents when applicable. Subject dietary data was not used if FFQs were invalid, with FFQs considered invalid if they were incomplete, or subjects reported excess (>30,000 kj/day) or deficient (<3500 kj/day) energy consumption, or excessive consumption of a single food group (≥11 serves/day) [28,31].
The smoking status of subjects was evaluated via an interviewer-administered questionnaire and subjects were categorised as non-smokers, ex-smokers or current smokers. Only 3% of subjects were current smokers, and therefore current smokers and those with a history of smoking were combined in the analyses. Subject’s body mass index (BMI (kg/m2)) was calculated using anthropometric measures taken during clinic appointments via standard procedures [27], and subjects were classified as underweight (BMI; <18.5), normal weight (BMI; 18.5–24.9), overweight (BMI; 25–29.9), or obese (BMI > 30). Only 2% of subjects were classified as underweight and were considered with those of normal BMIs in analyses.

2.4. Genotyping of Vitamin D and Folate Genetic Variants

Genomic DNA was isolated from whole blood samples using Qiagen QIAamp mini-kits following manufacturer’s protocols for blood samples (Qiagen, Hilden, Germany). Genotypes were assessed for 10 folate-related genetic variants (MTRR-rs1801394, MTR-rs1805087, MTHFR-rs1801133, MTHFR-rs1801131, SHMT-rs1979277, MTHFD1-rs2236225, RFC1-rs1051266, TYMS-rs11280056, TYMS-rs45445694, and DHFR-rs70991108) and 11 vitamin D-related variants (GC-rs4588, CYP2R1-rs10741657, DHCR7/NADSYN1-rs12785878, CYP24A1-rs17216707, VDR-rs4516035, VDR-rs757343, VDR-rs2228570, VDR-rs731236, VDR-rs7975232, VDR-rs11568820, and VDR-rs1544410). These variants were selected for assessment as they are common and well-characterized variants related to changes in vitamin status and activity (Figure 1).
Genotyping was undertaken via mixed PCR methods. The majority of variants were genotyped following previously outlined RFLP-PCR methods [35,36,37,38,39,40,41,42,43,44,45,46,47]; MTRR-rs1801394, MTR-rs1805087, MTHFR-rs1801133, MTHFR-rs1801131, SHMT-rs1979277 MTHFD1-rs2236225, RFC1-rs1051266, TYMS-rs11280056, GC-rs4588, CYP2R1-rs10741657, VDR-rs4516035, VDR-rs757343, VDR-rs2228570, VDR-rs731236, VDR-rs7975232, and VDR-rs1544410. Variants DHFR-rs70991108 and VDR-rs11568820 were assessed following allele-specific PCR [48,49], with TYMS-rs45445694 evaluated via PCR and gel electrophoresis [50]. Genotyping for DHCR7/NADSYN1-rs12785878 and CYP24A1-rs17216707 were performed through qPCR methods using Taqman genotyping assays (Assay IDs; C_32063037_10 and C_33659702_10).

2.5. Estimation of Environmental UVR Levels: Accumulated Area Erythemal Dose Rate

Coordinates of subjects’ reported residential location were used to gather data on erythemal dose rate (EDR) in subject location prior to clinic appointments, gathered from NASA’s Total Ozone Mapping Spectrometer (Accessed via https://giovanni.gsfc.nasa.gov/giovanni/). EDR is a measure of the potential for biological damage to be caused by UVR [51]. Information was gathered for the total amount of EDR accumulated over two time periods of interest—over the six weeks (6W-EDR) and four months (4M-EDR) prior to the subject’s clinic appointments. 4M-EDR was assessed in a previous investigation examining relationships between environmental UVR and Hcy and/or RBC folate levels [11], and was chosen due to four months being the approximate lifespan of an RBC (i.e., relevant to RBC folate turnover), respectively. 6W-EDR was also assessed in the current study, with six weeks being indicative of the time-lag between UVR exposure and serum 25(OH)D changes [52,53]. EDR values ranged considerably between RHLS subjects and captured varying seasons due to original subject recruitments and clinic appointments spanning > 18 months (>80 clinic dates). Further details on the method have been published previously [11,54,55].

2.6. Statistical Analyses

Statistical analyses were performed using JMP software (V.14.2.0; SAS Institute Inc., Cary, NC, USA). Descriptive statistics (means, 95% confidence intervals, ranges, and frequencies) were calculated and presented as appropriate. Multifactorial modelling (standard least squares regression) was used to examine associations between variables of interest, with multiplicative interaction terms included in models where appropriate. Adjustment variables included age, sex, and key determinants of Hcy levels—dietary intake of vitamin B6 and alcohol, RBC folate, vitamin B12 and creatinine levels, BMI category, and smoking status [56,57,58]. Mixed direction stepwise regression was undertaken using significance levels of p ≤ 0.250 or p > 0.250 to enter or remove variables from models where appropriate. Mallow’s Cp criterion was used for selecting the model where Cp first approaches p variables. Multiple comparisons of least-squares means were made using Tukey’s HSD post hoc tests. Adjusted R2 values and p-values are reported for final models, with standardised parameter estimates (β) and p-values reported for individual variables. Significance level (p < 0.05) was adjusted using the Bonferroni method [59] to account for multiple testing where appropriate.

3. Results

3.1. Subject Characteristics

The mean age of subjects was 77 years (Table 1), with 56% being female. Mean RBC folate, 25(OH)D and Hcy levels were within the reference ranges, with reference ranges defined as 317–1422 nmol/L for RBC folate [60] and 5–15 μmol/L for Hcy levels [61]. The mean 25(OH)D level was higher than the range considered adequate (20–28 ng/mL) in the Australian population [62]. However, 14% of the cohort had 25(OH)D levels indicative of vitamin D deficiency (<20 ng/mL) [62]. Furthermore, 16% of the cohort had mild elevations in Hcy levels (15–35 μmol/L); however Hcy levels beyond the reference range for Hcy (5–15 μmol/L) are expected in elderly populations [61]. There were no incidences of folate deficiency (RBC folate < 317 nmol/L).
Reported intake of alcohol, creatinine levels, vitamin B12 levels, and smoking status differed by sex. Mean reported alcohol intake and creatinine levels were higher in males compared to females (10.4 (95% CI; 9.0–11.4) vs. 4.4 (3.7–5.1); p < 0.001, and 10.6 (10.0–11.2) vs. 8.4 (7.8–8.9); p < 0.001). Vitamin B12 levels were lower in males 224.6 (207.4–241.8) vs. 251.0 (235.6–266.2); p = 0.03). More males had a history of smoking (64% vs. 32%; p < 0.001).
The frequencies of examined folate and vitamin D genetic variants in this cohort are outlined in Table S1. Genotype frequencies did not deviate from Hardy–Weinberg expectations, except for the GC-rs4588 variant (Table S1). Genotypes were further classified as binary variables, with analyses comparing presence vs. absence of the polymorphic allele. DHCR7/NADSYN1-rs12785878 and VDR-rs7975232 variants were excluded from further analysis due to uneven groups, with >95% of subjects carrying polymorphic alleles for these variants.

3.2. Independent and Interactive Influences of EDR and Levels of 25(OH)D and RBC Folate on Hcy Levels

Increases in EDR accumulated over six weeks prior to clinic appointment (6W-EDR) was associated with decreased Hcy levels (unadjusted model; β = −0.25, p < 0.001; Table 2). A similar negative association was also shown between 4M-EDR and Hcy levels as previously reported [11] (unadjusted model; β = −0.29, p < 0.001; Table 2). Both these associations remained following adjustments with 25(OH)D and RBC folate levels (β = −0.24, p < 0.001 and β = −0.29, p < 0.001; Model 1, Table 2), and following further adjustments with potential confounders (β = −0.24, p < 0.001 and β = −0.30, p < 0.001; Model 2, Table 2). Notably, neither 25(OH)D or RBC folate levels were significant independent predictors of Hcy levels in these models (Table 2; Models 1–2) or when the binary associations between Hcy and either 25(OH)D or RBC folate were considered (Figure S1).
Potential environment–nutrient interactions between EDR measures and levels of 25(OH)D and/or folate in determining Hcy levels were examined to consider if the effects of EDR measures were dependent on levels of UVR-sensitive vitamins. A significant interaction between 6M-EDR and 25(OH)D in predicting Hcy was observed (unadjusted model; pinteraction = 0.01; Table 3). This interaction remained following adjustments with potential Hcy determinants and when considering multiple testing corrections using the Bonferroni method (Adjusted model; pinteraction = 0.002; Table 2). The direct effect of 6W-EDR was also an independent predictor of Hcy levels in unadjusted and adjusted models (β = −0.24, p < 0.001 in both cases; Table 3). When stratifying the association between 6W-EDR and Hcy by quartiles of 25(OH)D levels, this association was only significant for subjects in the two lower quartiles of 25(OH)D levels (Q1; β = −0.35, p = 0.001, Q2; β = −0.35, p = 0.001; Table 4). However, the overall interaction between 6W-EDR and 25(OH)D quartiles (when assessed as a categorical variable) did not remain when applying strict adjustments to the p-value threshold via the Bonferroni method (pinteraction = 0.004; Table 4).

3.3. Independent and Interactive Influences of Vitamin D and/or Folate Genetic Variants on Hcy Levels

Stepwise regression for all examined vitamin D and folate genetic variants resulted in two folate-related variants being entered into the predictive model for Hcy (DHFR-rs70991108 and MTHFD1-rs2236225). Neither variant was found to be a significant independent predictor of Hcy levels in unadjusted and adjusted models (Table S2).
Potential interactions between these variants and RBC folate and/or EDR measures were also assessed to consider if influences of EDR or folate were modified by DHFR-rs70991108 and MTHFD1-rs2236225 variants. No significant interactions were observed between variants and RBC folate in predicting Hcy. However, MTHFD1-rs2236225 was a significant independent predictor of Hcy levels following consideration of interactions with RBC folate (β = −0.11, p = 0.02; adjusted model; Table S3). Presence of MTHFD1-rs2236225 polymorphic allele was associated with higher Hcy levels relative to ancestral homozygotes (10.7 μmol/L (95% CI; 8.9–12.6 μmol/L) vs. 9.3 μmol/L (7.4–11.3 μmol/L).
A significant interaction was found between MTHFD1-rs2236225 and 4M-EDR in predicting Hcy levels (adjusted models; pinteraction = 0.006; Table 5). MTHFD1-rs2236225, 4M-EDR, and 6M-EDR were also direct independent predictors of Hcy levels in these models (adjusted models; Table 5). When the association between 4M-EDR and Hcy levels was stratified by carriage of the MTHFD1-rs2236225 polymorphic allele, a significant association was only observed in subjects carrying the polymorphic allele (β = −0.28, p < 0.001).
Significant interactions identified between EDR measures, 25(OH)D levels, and/or MTHFD1-rs2236225 in previous analyses were considered together in multivariate models. When considering interactions between 6W-EDR and both 25(OH)D levels and MTHFD1-rs2236225 in a single multivariate model, only the interaction between 6W-EDR and 25(OH)D levels were significant (pinteraction ≤ 0.001; adjusted model; Table 6). The direct effect of 6W-EDR was also a significant predictor of Hcy levels in this model (β = −0.19, p = 0.001), with the overall model predicting 9% of variance in Hcy levels. When considering interactions between 4M-EDR and both 25(OH)D levels and MTHFD1-rs2236225, only the interaction between 4M-EDR and MTHFD1-rs2236225 remained significant (pinteraction = 0.006; adjusted model; Table 6). The direct effects of MTHFD1-rs2236225 and 4M-EDR were also significantly associated with Hcy levels in this model (β = −0.09, p = 0.04 and β = −0.21, p ≤ 0.001; adjusted model; Table 6). The final overall model predicted 11% of variance in Hcy levels.

4. Discussion

The presented findings demonstrate independent and interactive influences of environmental UVR, 25(OH)D levels, and folate variant MTHFD1-rs2236225 on Hcy levels in an elderly Australian cohort. These findings are novel, with this being the first study to consider the independent and interactive influences of environmental UVR levels, vitamin D and folate levels and related genetic factors on Hcy levels.
Environmental UVR levels accumulated over six weeks and four months before clinic appointments (i.e., 6W-EDR and 4M-EDR) were independent predictors of Hcy levels. Hcy levels were inversely associated with 4M-EDR previously [11], with Hcy levels further shown to be related to UVR levels over a shorter six week period (6W-EDR) in the current study. Notably, an interaction between 25(OH)D levels and 6W-EDR, but not 4M-EDR, in predicting Hcy levels was shown, indicating an interaction between 25(OH)D levels and a UVR time frame representative of the time-lag between UVR exposure and 25(OH)D changes [52]. When examining the directionality of this interaction, the association between 6W-EDR and Hcy levels was found to be only significant in subjects within the lower quartiles of 25(OH)D status (<33.26 ng/mL). This finding suggests that a relationship between 6W-EDR and 25(OH)D level exists until adequate vitamin D status is achieved. However, the direct effect of 6W-EDR remained a significant independent predictor of Hcy levels in models considering 6W-EDR/25(OH)D interactions, indicating that the relationship between 6W-EDR and Hcy levels is not solely explained by changes in 25(OH)D levels.
Carriage of the MTHFD1-rs2236225 polymorphic allele was associated with increased Hcy levels. This finding is supported by a previous investigation finding homozygosity for this variant to be related to increases in plasma Hcy in healthy Mexican American women (n = 43) [63]. MTHFD1 is the gene for the methylenetetrahydrofolate dehydrogenase, cyclohydrolase and formyltetrahydrofolate synthetase 1 (MTHFD1) trifunctional enzyme, which is required for the interconversion of active folate forms 5,10-methylenetetrahydrofolate (THF), 5,10-methenylTHF, and 10-formylTHF [64]. 5,10-methyleneTHF is required for Hcy regulation, being reduced to 5-methylTHF via action of 5,10-methylenetetrahydrofolate reductase (MTHFR), with 5-methylTHF then used as a methyl donor in the remethylation of Hcy into methionine [3]. Murine models of MTHFD1 insufficiency display increased Hcy levels, as well as decreases in methionine and cystathionine, indicative of impaired Hcy regulation [65,66]. The MTHFD1-rs2236225 is a missense variant which results in a thermolabile protein with a reduced half-life and synthetase enzyme activity [67]. This variant causes a reduction in 5,10-methyleneTHF formation, which may adversely affect Hcy regulation [68].
The MTHFD1-rs2236225 variant was found to interact with 4M-EDR to predict Hcy levels, with the association between 4M-EDR and Hcy levels dependent on the carriage of the polymorphic MTHFD1-rs2236225 allele. This association may be explained in part by the thermolability of the MTHFD1 variant. A previous investigation by the current authors found the association between environmental UVR levels and RBC folate levels to be dependent on the presence of another thermolabile variant: MTHFR-rs1801133 [26]. It is possible that the presence of thermolabile enzymes paired with increased UVR exposure leads to compensatory shifts in folate metabolism, which results in Hcy levels being influenced by UVR through a currently unknown mechanism.
Investigations examining the potential relationship between UVR and Hcy status are limited. Prior to studies undertaken by the current group, only two studies had examined Hcy in this context, with potential seasonal variations in Hcy examined, but with no significant findings reported [69,70]. However, these studies examined < 100 subjects residing in Northern European locations and therefore findings may reflect differences in statistical power and environmental UVR levels [69,70]. It is possible that a UVR-related decline in Hcy may only be evident in populations exposed to high levels of environmental UVR, such as Australians [29], with further investigations examining cohorts located in high UVR areas needed. The reported inverse relationship between UVR and Hcy levels should be interpreted with caution, as it may not indicate a benefit of UVR exposure on Hcy metabolism, but an increase in Hcy turnover and/or use in response to increased oxidative stress following UVR exposure. Increased oxidative stress may promote Hcy auto-oxidization and the formation of Hcy into other oxidants not detected by the Hcy assay used here [71,72].
An apparent limitation of this study is that Hcy levels were examined against environmental UVR levels rather than personal UVR exposures. Assessment of personal exposures was not possible due to this being a secondary retrospective cohort analysis. However, given that original subject recruitments and clinic appointments in the RHLS spanned > 18 months (> 80 clinic dates), collective environmental UVR data captured a wide range of UVR levels, with reported findings providing a justification for future investigations using more precise methods. Further strengths of this investigation were that multiple vitamin D and folate variants were examined, with the RHLS cohort well-characterised in respect to several known Hcy determinants, which were corrected for in subsequent analyses. 25(OH)D levels were examined via a HPLC method [33], with HPLC methods shown to be a more precise and sensitive method for 25(OH)D measurement compared to other methods [73]. However, a limitation of this method is that external quality assessment was not obtained via external quality assessment schemes (e.g., via DEQAS). The focus of the current study was to examine these influences in an elderly Australian cohort, given that this sub-population is at heightened risk of UVR-related cellular damage [29,30]. Several novel independent and interactive influences of environmental UVR levels and 25(OH)D levels and folate variant MTHFD1-rs2236225 on Hcy levels are presented. However, findings may be age-specific, given age-related decreases in UVR-protection [30], and increases in oxidative stress [74], with further investigation required before findings can be generalised to other populations.

5. Conclusions

Accumulated environmental UVR levels and MTHFD1-rs2236225 were found to be independent predictors of Hcy status in an elderly Australian cohort. A nutrient–environment interaction between 25(OH)D levels and 6W-EDR was demonstrated, with the association between 6W-EDR and Hcy only significant for subjects in lower quartiles of 25(OH)D status. Furthermore, a gene–environment interaction between MTHFR-rs1801133 and 4M-EDR in predicting Hcy was found, with an association between 4M-EDR and Hcy only observed in subjects carrying the polymorphic MTHFD1-rs2236225 allele. Presented findings are novel, with further investigation needed to examine if these independent and interactive effects are population or age-specific, and whether they influence the risk of Hcy-related health outcomes.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2072-6643/12/5/1455/s1, Figure S1: Binary associations between biochemical variables of interest (Hcy, 25(OH)D and RBC folate levels), Table S1: Allelic and genotypic frequencies for folate and vitamin D variants and assessment for deviation from Hardy–Weinberg equilibrium, Table S2: Folate variants identified by stepwise regression for inclusion in models for Hcy prediction, with and without adjustments, Table S3: Assessment for gene-nutrient interactions in predicting Hcy, with and without adjustments for determinants of Hcy levels.

Author Contributions

Conceptualization by P.J., E.B., and M.L.; methodology developed by all authors; genetics examined by P.J., C.M., and E.B.; clinical analysis and oversight by R.T., M.G., and M.V.; software, P.J., E.B., and C.J.S.; formal analysis, P.J., and E.B.; resources, E.B., M.L., M.V., Z.Y., and C.J.S.; data curation, M.V., M.L., and E.B.; supervision, M.L., C.J.S., M.V., and E.B.; writing—original draft preparation, P.J. and E.B.; writing—review and editing, all authors; project administration, M.V., M.L., and E.B.; funding acquisition, M.V., M.L., and E.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research conducted as part of the Retirement Health and Lifestyle Study, with initial and ongoing funding provided by the Australian Research Council (G0188386), Central Coast Local Health District Public Health Unit (G0190658/G1700259), UnitingCare Ageing NSW/ACT (G0189230), Urbis Pty Ltd. (G0189232), Valhalla Village Pty Ltd. (G1000936), and Hunter Valley Research Foundation.

Acknowledgments

Patrice Jones is supported by an Australian Government Research Training Program scholarship and a Hunter Medical Research Institute (HMRI) Greaves Family Scholarship. Emma Beckett is supported by a National Health and Medical Research Council (NHMRC) Early Career Fellowship. The authors would also like to thank the RHLS participants, and researchers and clinic staff involved in the broader RHLS study including Katrina King, Suzanne Niblett, John Furst, Paul Roach, Peter Lewis, Jenny Marriott, Marie Mazaroli, Elizabeth Death, Jodi Humphreys, and Louise Lambeth.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Škovierová, H.; Vidomanová, E.; Mahmood, S.; Sopková, J.; Drgová, A.; Červeňová, T.; Halašová, E.; Lehotský, J. The Molecular and Cellular Effect of Homocysteine Metabolism Imbalance on Human Health. Int. J. Mol. Sci. 2016, 17, 1733. [Google Scholar] [CrossRef] [PubMed]
  2. Smith, A.D.; Refsum, H.; Bottiglieri, T.; Fenech, M.; Hooshmand, B.; McCaddon, A.; Miller, J.W.; Rosenberg, I.H.; Obeid, R. Homocysteine and Dementia: An International Consensus Statement. J. Alzheimer’s Dis. 2018, 62, 561–570. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Blom, H.J.; Smulders, Y. Overview of homocysteine and folate metabolism. With special references to cardiovascular disease and neural tube defects. J. Inherit. Metab. Dis. 2010, 34, 75–811. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Stover, P.J.; Field, M.S. Trafficking of Intracellular Folates. Advances in Nutrition. Int. Rev. J. 2011, 2, 325–331. [Google Scholar]
  5. Juzeniene, A.; Thu Tam, T.T.; Iani, V.; Moan, J. 5-Methyltetrahydrofolate can be photodegraded by endogenous photosensitizers. Free. Radic. Boil. Med. 2009, 47, 1199–1204. [Google Scholar] [CrossRef]
  6. Fukuwatari, T.; Fujita, M.; Shibata, K. Effects of UVA irradiation on the concentration of folate in human blood. Biosci. Biotechnol. Biochem. 2009, 73, 322–327. [Google Scholar] [CrossRef] [Green Version]
  7. El-Saie, L.T.; Rabie, A.R.; Kamel, M.I.; Seddeik, A.K.; Elsaie, M.L. Effect of narrowband ultraviolet B phototherapy on serum folic acid levels in patients with psoriasis. Lasers Med. Sci. 2011, 26, 481–485. [Google Scholar] [CrossRef]
  8. Borradale, D.; Isenring, E.; Hacker, E.; Kimlin, M.G. Exposure to solar ultraviolet radiation is associated with a decreased folate status in women of childbearing age. J. Photochem. Photobiol. B Boil. 2014, 131, 90–95. [Google Scholar] [CrossRef] [Green Version]
  9. Ma, S.; Fattah, A.; El-Borhamy, M.I. Analysis of Serum Folate Levels after Narrow Band UVB Exposure. Egypt. Derm. Online J. 2006, 2, 1–7. [Google Scholar]
  10. Murase, J.E.; Koo, J.Y.M.; Berger, T.G. Narrowband ultraviolet B phototherapy influences serum folate levels in patients with vitiligo. J. Am. Acad. Dermatol. 2010, 62, 710–711. [Google Scholar] [CrossRef]
  11. Jones, P.; Lucock, M.; Scarlett, C.J.; Veysey, M.; Beckett, E. Environmental UVR Levels and Skin Pigmentation Gene Variants Associated with Folate and Homocysteine Levels in an Elderly Cohort. Int. J. Environ. Res. Public Health 2020, 17, 1545. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Jones, P.; Beckett, E.; Yates, Z.; Veysey, M.; Lucock, M. Converging Evolutionary, Environmental and Clinical Ideas on Folate Metabolism. Explor. Res. Hypothesis Med. 2016, 1, 34–41. [Google Scholar]
  13. Jones, P.; Lucock, M.; Veysey, M.; Jablonski, N.; Chaplin, G.; Beckett, E. Frequency of folate-related polymorphisms varies by skin pigmentation. Am. J. Hum. Boil. 2017, 30, e23079. [Google Scholar] [CrossRef] [PubMed]
  14. Amer, M.; Qayyum, R. The Relationship Between 25-Hydroxyvitamin D and Homocysteine in Asymptomatic Adults. J. Clin. Endocrinol. Metab. 2013, 99, 633–638. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Pham, T.M.; Ekwaru, J.P.; Mastroeni, S.S.; Mastroeni, M.F.; Loehr, S.A.; Veugelers, P.J. The Effect of Serum 25-Hydroxyvitamin D on Elevated Homocysteine Concentrations in Participants of a Preventive Health Program. PLoS ONE 2016, 11, e0161368. [Google Scholar] [CrossRef] [PubMed]
  16. Glueck, C.J.; Jetty, V.; Rothschild, M.; Duhon, G.; Shah, P.; Prince, M.; Lee, K.; Goldenberg, M.; Kumar, A.; Goldenberg, N.; et al. Associations between Serum 25-hydroxyvitamin D and Lipids, Lipoprotein Cholesterols, and Homocysteine. N. Am. J. Med Sci. 2016, 8, 284–290. [Google Scholar] [CrossRef] [Green Version]
  17. Mao, X.; Xing, X.; Xu, R.; Gong, Q.; He, Y.; Li, S.; Wang, H.; Liu, C.; Ding, X.; Na, R.; et al. Folic Acid and Vitamins D and B12 Correlate With Homocysteine in Chinese Patients With Type-2 Diabetes Mellitus, Hypertension, or Cardiovascular Disease. Medicine 2016, 95, e2652. [Google Scholar] [CrossRef]
  18. Moretti, R.; Caruso, P.; Dal Ben, M.; Conti, C.; Gazzin, S.; Tiribelli, C. Vitamin D, Homocysteine, and Folate in Subcortical Vascular Dementia and Alzheimer Dementia. Front. Aging Neurosci. 2017, 9, 169. [Google Scholar] [CrossRef] [Green Version]
  19. Kriebitzsch, C.; Verlinden, L.; Eelen, G.; van Schoor, N.M.; Swart, K.; Lips, P.; Meyer, M.B.; Pike, J.W.; Boonen, S.; Carlberg, C.; et al. 1,25-dihydroxyvitamin D(3) influences cellular homocysteine levels in murine pre-osteoblastic MC3T3-E1 cells by direct regulation of cystathionine β-synthase. Journal of bone and mineral research. Off. J. Am.Soc. Bone Miner. Res. 2011, 26, 2991–3000. [Google Scholar] [CrossRef] [Green Version]
  20. Lucock, M.; Jones, P.; Martin, C.; Beckett, E.; Yates, Z.; Furst, J.; Veysey, M. Vitamin D: Beyond Metabolism. J. Evid. Based Complementary Altern. Med. 2015, 20, 310–322. [Google Scholar] [CrossRef] [Green Version]
  21. Tiosano, D.; Audi, L.; Climer, S.; Zhang, W.; Templeton, A.R.; Fernández-Cancio, M.; Gershoni-Baruch, R.; Sánchez-Muro, J.M.; El Kholy, M.; Hochberg, Z. Latitudinal Clines of the Human Vitamin D Receptor and Skin Color Genes. G3 Genes Genomes Genet. 2016, 6, 1251–1266. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Amato, R.; Pinelli, M.; Monticelli, A.; Miele, G.; Cocozza, S. Schizophrenia and vitamin D related genes could have been subject to latitude-driven adaptation. BMC Evol. Boil. 2010, 10, 351. [Google Scholar] [CrossRef] [PubMed]
  23. Kuan, V.; Martineau, A.R.; Griffiths, C.J.; Hyppönen, E.; Walton, R. DHCR7 mutations linked to higher vitamin D status allowed early human migration to northern latitudes. BMC Evol. Boil. 2013, 13, 144. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Constans, J.; Hazout, S.; Garruto, R.M.; Gajdusek, D.C.; Spees, E.K. Population distribution of the human vitamin D binding protein: Anthropological considerations. Am. J. Phys. Anthr. 1985, 68, 107–122. [Google Scholar] [CrossRef] [PubMed]
  25. Jones, P.; Lucock, M.; Chaplin, G.; Jablonski, N.G.; Veysey, M.; Scarlett, C.; Beckett, E. Distribution of variants in multiple vitamin D-related loci (DHCR7/NADSYN1, GC, CYP2R1, CYP11A1, CYP24A1, VDR, RXRα and RXRγ) vary between European, East-Asian and Sub-Saharan African-ancestry populations. Genes. Nutri. 2020, 15, 5. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Lucock, M.; Beckett, E.; Martin, C.; Jones, P.; Furst, J.; Yates, Z.; Jablonski, N.G.; Chaplin, G.; Veysey, M. UV-associated decline in systemic folate: Implications for human nutrigenetics, health, and evolutionary processes. Am. J. Hum. Boil. 2016, 29, e22929. [Google Scholar] [CrossRef]
  27. Abbott, K.A.; Veysey, M.; Lucock, M.; Niblett, S.; King, K.; Burrows, T.; Garg, M.L. Sex-dependent association between erythrocyte n-3 PUFA and type 2 diabetes in older overweight people. Br. J. Nutr. 2016, 115, 1379–1386. [Google Scholar] [CrossRef] [Green Version]
  28. Beckett, E.L.; Martin, C.; Boyd, L.; Porter, T.; King, K.; Niblett, S.; Yates, Z.; Veysey, M.; Lucock, M. Reduced plasma homocysteine levels in elderly Australians following mandatory folic acid fortification—A comparison of two cross-sectional cohorts. J. Nutr. Intermed. Metab. 2017, 8, 14–20. [Google Scholar] [CrossRef]
  29. Godar, D.E. UV doses worldwide. Photochem. Photobiol. 2005, 81, 736–749. [Google Scholar] [CrossRef]
  30. Baumann, L. Skin ageing and its treatment. J. Pathol. 2007, 211, 241–251. [Google Scholar] [CrossRef]
  31. Lucock, M.; Yates, Z.; Martin, C.; Choi, J.-H.; Beckett, E.; Boyd, L.; LeGras, K.; Ng, X.; Skinner, V.; Wai, R.; et al. Methylation diet and methyl group genetics in risk for adenomatous polyp occurrence. BBA Clin. 2015, 3, 107–112. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Tan, Y.; Hoffman, R.M. A highly sensitive single-enzyme homocysteine assay. Nat. Protoc. 2008, 3, 1388–1394. [Google Scholar] [CrossRef] [PubMed]
  33. Turpeinen, U.; Hohenthal, U.; Stenman, U.H. Determination of 25-hydroxyvitamin D in serum by HPLC and immunoassay. Clin. Chem. 2003, 49, 1521–1524. [Google Scholar] [CrossRef] [PubMed]
  34. Ambrosini, G.L.; Mackerras, D.; de Klerk, N.H.; Musk, A.W. Comparison of an Australian food-frequency questionnaire with diet records: Implications for nutrition surveillance. Public Heal. Nutr. 2003, 6, 415–422. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Murthy, J.; Gurramkonda, V.B.; Bhaskar, L. Genetic variant in MTRR A66G, but not MTR A2756G, is associated with risk of non-syndromic cleft lip and palate in Indian population. J. Oral Maxillofac. Surg. Med. Pathol. 2015, 27, 782–785. [Google Scholar] [CrossRef]
  36. Van der Put, N.M.J.; Trijbels, F.J.M.; van den Heuvel, L.P.; Blom, H.J.; Steegers-Theunissen, R.P.M.; Eskes, T.K.A.B.; Mariman, E.C.M.; den Heyer, M.; Frosst, P.; Rozen, R. Mutated methylenetetrahydrofolate reductase as a risk factor for spina bifida. Lancet 1995, 346, 1070–1071. [Google Scholar] [CrossRef] [Green Version]
  37. Van der Put, N.M.J.; Gabreëls, F.; Stevens, E.M.B.; Smeitink, J.A.M.; Trijbels, F.J.M.; Eskes, T.K.A.B.; van den Heuvel, L.P.; Blom, H.J. A Second Common Mutation in the Methylenetetrahydrofolate Reductase Gene: An Additional Risk Factor for Neural-Tube Defects? Am. J. Hum. Genet. 1998, 62, 1044–1051. [Google Scholar] [CrossRef] [Green Version]
  38. Heil, S.G.; Van der Put, N.M.J.; Waas, E.T.; den Heijer, M.; Trijbels, F.J.M.; Blom, H.J. Is Mutated Serine Hydroxymethyltransferase (SHMT) Involved in the Etiology of Neural Tube Defects? Mol. Genet. Metab. 2001, 73, 164–172. [Google Scholar] [CrossRef]
  39. Wang, L.; Ke, Q.; Chen, W.; Wang, J.; Tan, Y.; Zhou, Y.; Hua, Z.; Ding, W.; Niu, J.; Shen, J.; et al. Polymorphisms of MTHFD, plasma homocysteine levels, and risk of gastric cancer in a high-risk Chinese population. Clin. Cancer Res. 2007, 13, 2526–2532. [Google Scholar] [CrossRef] [Green Version]
  40. Dufficy, L.; Naumovski, N.; Ng, X.; Blades, B.; Yates, Z.; Travers, C.; Lewis, P.; Sturm, J.; Veysey, M.; Roach, P.D.; et al. G80A reduced folate carrier SNP influences the absorption and cellular translocation of dietary folate and its association with blood pressure in an elderly population. Life Sci. 2006, 79, 957–966. [Google Scholar] [CrossRef]
  41. Ulrich, C.M.; Bigler, J.; Velicer, C.M.; Greene, E.A.; Farin, F.M.; Potter, J.D. Searching expressed sequence tag databases: Discovery and confirmation of a common polymorphism in the thymidylate synthase gene. Cancer Epidemiol. Biomark. Prev. 2000, 9, 1381–1385. [Google Scholar]
  42. Lafi, Z.M.; Irshaid, Y.M.; El-Khateeb, M.; Ajlouni, K.M.; Hyassat, D. Association of rs7041 and rs4588 Polymorphisms of the Vitamin D Binding Protein and the rs10741657 Polymorphism of CYP2R1 with Vitamin D Status Among Jordanian Patients. Genet. Test. Mol. Biomark. 2015, 19, 629–636. [Google Scholar] [CrossRef] [PubMed]
  43. Penna-Martinez, M.; Ramos-Lopez, E.; Stern, J.; Kahles, H.; Hinsch, N.; Hansmann, M.-L.; Selkinski, I.; Grünwald, F.; Vorländer, C.; Bechstein, W.O.; et al. Impaired vitamin D activation and association with CYP24A1 haplotypes in differentiated thyroid carcinoma. Thyroid 2012, 22, 709–716. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Halsall, J.A.; Osborne, J.E.; Epstein, M.P.; Pringle, J.H.; Hutchinson, P.E. The unfavorable effect of the A allele of the vitamin D receptor promoter polymorphism A-1012G has different mechanisms related to susceptibility and outcome of malignant melanoma. Dermato-Endocrinology 2009, 1, 54–57. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Gong, Y.-L.; Xie, D.-W.; Deng, Z.-L.; Bostick, R.M.; Miao, X.-J.; Zhang, J.-H.; Gong, Z.-H. Vitamin D receptor gene Tru9I polymorphism and risk for incidental sporadic colorectal adenomas. World J. Gastroenterol. 2005, 11, 4794–4799. [Google Scholar] [CrossRef]
  46. Bhanushali, A.A.; Lajpal, N.; Kulkarni, S.S.; Chavan, S.S.; Bagadi, S.S.; Das, B.R. Frequency of fokI and taqI polymorphism of vitamin D receptor gene in Indian population and its association with 25-hydroxyvitamin D levels. Indian J. Hum. Genet. 2009, 15, 108–113. [Google Scholar]
  47. Zmuda, J.M.; Cauley, J.A.; Danielson, M.E.; Wolf, R.L.; Ferrell, R.E. Vitamin D Receptor Gene Polymorphisms, Bone Turnover, and Rates of Bone Loss in Older African-American Women. J. Bone Miner. Res. 1997, 12, 1446–1452. [Google Scholar] [CrossRef]
  48. Johnson, W.G.; Stenroos, E.S.; Spychala, J.R.; Chatkupt, S.; Ming, S.X.; Buyske, S. New 19 bp deletion polymorphism in intron-1 of dihydrofolate reductase (DHFR): A risk factor for spina bifida acting in mothers during pregnancy? Am. J. Med Genet. 2004, 124, 339–345. [Google Scholar] [CrossRef] [PubMed]
  49. Selvaraj, P.; Alagarasu, K.; Harishankar, M.; Vidyarani, M.; Narayanan, P.R. Regulatory region polymorphisms of vitamin D receptor gene in pulmonary tuberculosis patients and normal healthy subjects of south India. Int. J. Immunogenet. 2008, 35, 251–254. [Google Scholar] [CrossRef]
  50. Horie, N.; Aiba, H.; Oguro, K.; Hojo, H.; Takeishi, K. Functional Analysis and DNA Polymorphism of the Tandemly Repeated Sequences in the 5’-terminal Regulatory Region of the Human Gene for Thymidylate Synthase. Cell Struct. Funct. 1995, 20, 191–197. [Google Scholar] [CrossRef] [Green Version]
  51. Standard, C. Erythema reference action spectrum and standard erythema dose. Cie S 1998, 7, E1998. [Google Scholar]
  52. Lucas, J.A.; Bolland, M.J.; Grey, A.B.; Ames, R.W.; Mason, B.H.; Horne, A.M.; Gamble, G.D.; Reid, I.R. Determinants of vitamin D status in older women living in a subtropical climate. Osteoporos. Int. 2005, 16, 1641–1648. [Google Scholar] [CrossRef] [PubMed]
  53. Beckett, E.L.; Jones, P.; Veysey, M.; Duesing, K.; Martin, C.; Furst, J.; Yates, Z.; Jablonski, N.G.; Chaplin, G.; Lucock, M. VDR gene methylation as a molecular adaption to light exposure: Historic, recent and genetic influences. Am. J. Hum. Boil. 2017, 29, e23010. [Google Scholar] [CrossRef] [PubMed]
  54. Lucock, M.; Yates, Z.; Martin, C.; Choi, J.-H.; Boyd, L.; Tang, S.; Naumovski, N.; Furst, J.; Roach, P.; Jablonski, N.; et al. Vitamin D, folate, and potential early lifecycle environmental origin of significant adult phenotypes. Evol. Med. Public Heal. 2014, 2014, 69–91. [Google Scholar] [CrossRef] [Green Version]
  55. Beckett, E.L.; Martin, C.; Duesing, K.; Jones, P.; Furst, J.; Yates, Z.; Veysey, M.; Lucock, M. Vitamin D Receptor Genotype Modulates the Correlation between Vitamin D and Circulating Levels of let-7a/b and Vitamin D Intake in an Elderly Cohort. J. Nutr. Nutr. 2014, 7, 264–273. [Google Scholar] [CrossRef]
  56. Rasmussen, L.B.; Ovesen, L.; Bulow, I.; Knudsen, N.; Laurberg, P.; Perrild, H. Folate intake, lifestyle factors, and homocysteine concentrations in younger and older women. Am. J. Clin. Nutr. 2000, 72, 1156–1163. [Google Scholar] [CrossRef] [Green Version]
  57. Ganji, V.; Kafai, M.R. Demographic, health, lifestyle, and blood vitamin determinants of serum total homocysteine concentrations in the third National Health and Nutrition Examination Survey, 1988–1994. Am. J. Clin. Nutr. 2003, 77, 826–833. [Google Scholar] [CrossRef] [Green Version]
  58. Jacques, P.F.; Bostom, A.G.; Wilson, P.W.; Rich, S.; Rosenberg, I.H.; Selhub, J. Determinants of plasma total homocysteine concentration in the Framingham Offspring cohort. Am. J. Clin. Nutr. 2001, 73, 613–621. [Google Scholar] [CrossRef] [Green Version]
  59. Bland, J.M.; Altman, D.G. Multiple significance tests: The Bonferroni method. BMJ 1995, 310, 170. [Google Scholar] [CrossRef] [Green Version]
  60. Department of Health, Folate Testing. 2014. Available online: https://www.health.gov.au/internet/main/publishing.nsf/Content/F05A8741F610EDB7CA257EB30026794E/$File/Folate%20testing%20Review%20Report.docx (accessed on 1 November 2019).
  61. Maron, B.A.; Loscalzo, J. The Treatment of Hyperhomocysteinemia. Annu. Rev. Med. 2009, 60, 39–54. [Google Scholar] [CrossRef] [Green Version]
  62. Nowson, C.A.; McGrath, J.J.; Ebeling, P.R.; Haikerwal, A.; Daly, R.M.; Sanders, K.M.; Seibel, M.J.; Mason, R.S. Vitamin D and health in adults in Australia and New Zealand: A position statement. Med. J. Aust. 2012, 196, 686–687. [Google Scholar] [CrossRef] [PubMed]
  63. Ivanov, A.; Nash-Barboza, S.; Hinkis, S.; Caudill, M.A. Genetic Variants in Phosphatidylethanolamine N-methyltransferase and Methylenetetrahydrofolate Dehydrogenase Influence Biomarkers of Choline Metabolism When Folate Intake Is Restricted. J. Am. Diet. Assoc. 2009, 109, 313–318. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Anguera, M.C.; Field, M.S.; Perry, C.; Ghandour, H.; Chiang, E.P.; Selhub, J.; Shane, B.; Stover, P.J. Regulation of folate-mediated one-carbon metabolism by 10-formyltetrahydrofolate dehydrogenase. J. Boil. Chem. 2006, 281, 18335–18342. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Field, M.S.; Shields, K.S.; Abarinov, E.V.; Malysheva, O.V.; Allen, R.H.; Stabler, S.P.; Ash, J.A.; Strupp, B.J.; Stover, P.J.; Caudill, M.A. Reduced MTHFD1 activity in male mice perturbs folate- and choline-dependent one-carbon metabolism as well as transsulfuration. J. Nutr. 2012, 143, 41–45. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. MacFarlane, A.J.; Perry, C.A.; Girnary, H.H.; Gao, D.; Allen, R.H.; Stabler, S.P.; Shane, B.; Stover, P.J. Mthfd1 is an essential gene in mice and alters biomarkers of impaired one-carbon metabolism. J. Boil. Chem. 2008, 284, 1533–1539. [Google Scholar] [CrossRef] [Green Version]
  67. Sutherland, H.G.; Hermile, H.; Sanche, R.; Menon, S.; Lea, R.A.; Haupt, L.M.; Griffiths, L.R. Association study of MTHFD1 coding polymorphisms R134K and R653Q with migraine susceptibility. Headache 2014, 54, 1506–1514. [Google Scholar] [CrossRef] [Green Version]
  68. Christensen, K.E.; Rohlicek, C.V.; Andelfinger, G.U.; Michaud, J.; Bigras, J.L.; Richter, A.; Mackenzie, R.E.; Rozen, R. The MTHFD1 p.Arg653Gln variant alters enzyme function and increases risk for congenital heart defects. Hum. Mutat. 2009, 30, 212–220. [Google Scholar] [CrossRef]
  69. McKinley, M.C.; Strain, J.J.; McPartlin, J.; Scott, J.M.; McNulty, H. Plasma Homocysteine Is Not Subject to Seasonal Variation. Clin. Chem. 2001, 47, 1430–1436. [Google Scholar] [CrossRef]
  70. Clarke, R.; Woodhouse, P.; Ulvik, A.; Frost, C.; Sherliker, P.; Refsum, H.; Ueland, P.M.; Khaw, K.T. Variability and determinants of total homocysteine concentrations in plasma in an elderly population. Clin. Chem. 1998, 44, 102–107. [Google Scholar] [CrossRef] [Green Version]
  71. Hayden, M.R.; Tyagi, S.C. Homocysteine and reactive oxygen species in metabolic syndrome, type 2 diabetes mellitus, and atheroscleropathy: The pleiotropic effects of folate supplementation. Nutr. J. 2004, 3, 4. [Google Scholar] [CrossRef] [Green Version]
  72. Zou, C.G.; Banerjee, R. Homocysteine and Redox Signaling. Antioxidants Redox Signal. 2005, 7, 547–559. [Google Scholar] [CrossRef] [PubMed]
  73. Rezayi, M.; Ghayour-Mobarhan, M.; Tavakoly Sany, S.B.; Fani, M.; Avan, A.; Pasdar, Z.; Ferns, G.A.; Abouzari-Lotf, E.; Amiri, I.S. A comparison of analytical methods for measuring concentrations of 25-hydroxy vitamin D in biological samples. Anal. Methods 2018, 10, 5599–5612. [Google Scholar] [CrossRef]
  74. Cui, H.; Kong, Y.; Zhang, H. Oxidative stress, mitochondrial dysfunction, and aging. J. Signal Transduct. 2011, 2012, 1–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Nutrients 12 01455 g001 550
Figure 1. Examined genetic variants involved in folate and vitamin D metabolic pathways. Examined folate-related genetic variants reside within key genes coding for the enzymes (blue) involved in three interconnected pathways which support activities of folate in the methionine cycle, which allows for regulation homocysteine, and pathways related to the synthesis of thymidylate (thymine precursor) and purines. These pathways are interconnected as they rely on the same pool of reduced folates within the cell, so disruption to any of these pathways via genetic variation can lead to changes in Hcy levels [4]. RFC1 is a carrier protein which transports circulating folates into target cells. Hcy may be regulated through an alternative folate-independent process—the transsulfuration pathway. Hcy is converted to cysteine in this pathway via CSE and CBS. Vitamin D is previously shown to influence CBS activity [19], and this is the suggested mechanism, explaining previous reports of Hcy and vitamin D levels being inversely related [15,16]. Examined vitamin D-related genetic variants reside within genes coding for the key enzymes (red) involved in vitamin D synthesis and metabolism, which may influence vitamin D levels and/or activity and therefore have the potential to influence vitamin D’s role in influencing Hcy levels. Abbreviations: Folate metabolism; AICAR, 5-aminoimidazole-4-carboxamide ribonucleotide; CBS, cystathionine beta-synthase; CSE, Cystathionine gamma-lyase; DHF, dihydrofolate; DHFR, dihydrofolate reductase; dTMP, deoxyuridine monophosphate; dUMP, deoxyuridine monophosphate; GAR, glycinamide ribonucleotide; Hcy, homocysteine; MAT, adenosylmethionine synthetase; Met, methionine; MTHFD1, methylenetetrahydrofolate dehydrogenase 1; MTHFR, 5,10-methyleneTHF reductase; MTR, methionine synthase; MTRR, methionine synthase reductase; RFC1, reduced folate carrier gene; SAH, S-adenosyl hcy; SAHH, S-adenosylhomocysteine hydrolase; SAM, S-adenosyl methionine; SHMT, serinehydroxymethyl transferase; THF, tetrahydrofolate; TYMS, thymidylate synthase. Vitamin D metabolism; CYP24A1, 24-hydroxylase; CYP27B1, 1α-hydroxylase; CYP2R1, 25-hydroxylase; DBP/GC, vitamin D-binding protein; DHCR7, 7-dehydrocholesterol reductase; VDR, vitamin D receptor. Figure created with BioRender.
Figure 1. Examined genetic variants involved in folate and vitamin D metabolic pathways. Examined folate-related genetic variants reside within key genes coding for the enzymes (blue) involved in three interconnected pathways which support activities of folate in the methionine cycle, which allows for regulation homocysteine, and pathways related to the synthesis of thymidylate (thymine precursor) and purines. These pathways are interconnected as they rely on the same pool of reduced folates within the cell, so disruption to any of these pathways via genetic variation can lead to changes in Hcy levels [4]. RFC1 is a carrier protein which transports circulating folates into target cells. Hcy may be regulated through an alternative folate-independent process—the transsulfuration pathway. Hcy is converted to cysteine in this pathway via CSE and CBS. Vitamin D is previously shown to influence CBS activity [19], and this is the suggested mechanism, explaining previous reports of Hcy and vitamin D levels being inversely related [15,16]. Examined vitamin D-related genetic variants reside within genes coding for the key enzymes (red) involved in vitamin D synthesis and metabolism, which may influence vitamin D levels and/or activity and therefore have the potential to influence vitamin D’s role in influencing Hcy levels. Abbreviations: Folate metabolism; AICAR, 5-aminoimidazole-4-carboxamide ribonucleotide; CBS, cystathionine beta-synthase; CSE, Cystathionine gamma-lyase; DHF, dihydrofolate; DHFR, dihydrofolate reductase; dTMP, deoxyuridine monophosphate; dUMP, deoxyuridine monophosphate; GAR, glycinamide ribonucleotide; Hcy, homocysteine; MAT, adenosylmethionine synthetase; Met, methionine; MTHFD1, methylenetetrahydrofolate dehydrogenase 1; MTHFR, 5,10-methyleneTHF reductase; MTR, methionine synthase; MTRR, methionine synthase reductase; RFC1, reduced folate carrier gene; SAH, S-adenosyl hcy; SAHH, S-adenosylhomocysteine hydrolase; SAM, S-adenosyl methionine; SHMT, serinehydroxymethyl transferase; THF, tetrahydrofolate; TYMS, thymidylate synthase. Vitamin D metabolism; CYP24A1, 24-hydroxylase; CYP27B1, 1α-hydroxylase; CYP2R1, 25-hydroxylase; DBP/GC, vitamin D-binding protein; DHCR7, 7-dehydrocholesterol reductase; VDR, vitamin D receptor. Figure created with BioRender.
Nutrients 12 01455 g001
Table
Table 1. Subject Characteristics.
Table 1. Subject Characteristics.
Male (n = 273)Female (n = 346)All (n = 619)
Continuous VariablesMean95% CIRangeMean95% CIRangeMean95% CIRange
Age7776–7865–937776–7865–957776–7865–95
RBC folate (nmol/L)1625.71311.5–1939.8931.0–2539.01342.21292.0–1392.5385.0–2695.01340.91304.8–1377.1381.0–2695.0
25(OH)D (ng/mL)31.524.6–38.49.7–44.735.033.3–36.79.3–94.435.834.6–37.09.3–94.4
Hcy (μmol/L)11.310.5–12.20.1–47.410.19.5–10.60.1–31.210.49.9–10.90.1–47.4
Serum vitamin B12 (pmol/L)224.6207.4–241.812.0–1116.0251.0235.6–266.268.0–1500.0239.3227.8–251.012.0–1500.0
Creatinine (μmol/L)10.610.0–11.21.4–51.68.47.8–8.90.5–29.59.49.0–9.80.5–51.6
Vitamin B6 intake (mg/d)9.04.9–13.10.0–220.79.26.9–11.50.0–203.08.66.9–10.30.0–220.7
Alcohol intake (g/day)10.49.0–11.70.0–40.44.43.7–5.10.0–37.77.06.3–7.70.0–40.4
6W-EDR5920.85577.9–6263.72140.5–11057.45959.95662.1–6257.62140.5–11057.46014.45718.4–6166.92149.5–11057.4
4M-EDR16795.515888.5–17702.57788.2–28258.016739.015939.6–17538.47788.2–29160.616764.716165.8–17361.97788.2–29160.6
Categorical Variablesn% n% n%
Tea serves/day
<17830 10834 19332
1–211846 14145 27446
>26124 6721 13222
Coffee serves/day
<19035 11336 20835
1–211244 13543 26244
>25521 6822 12922
Smoking status
Current or ex-smoker18568 12636 30850
Never smoked8832 22064 31150
BMI category
Underweight or normal4920 8226 12223
Overweight12148 12840 25743
Obese8032 10734 19934
Hcy: Homocysteine, 4M-EDR: 4 month accumulated erythemal dose rate, 6W-EDR: 6 week accumulated erythemal dose rate, 2D-EDR: 2 day accumulated erythemal dose rate, BMI: body mass index.
Table
Table 2. Associations between accumulated erythemal dose rate (EDR) measures, and levels of 25(OH)D levels, and RBC folate on Hcy levels, with and without adjustments for Hcy confounders.
Table 2. Associations between accumulated erythemal dose rate (EDR) measures, and levels of 25(OH)D levels, and RBC folate on Hcy levels, with and without adjustments for Hcy confounders.
Hcy Levels
Unadjusted (n = 618)Model 1 (n = 579)Model 2 (n = 464)
βpβpβp
6W-EDR−0.25<0.001−0.24<0.0010.24<0.001
25(OH)D levels- 0.040.3−0.030.5
RBC folate levels- −0.020.7−0.030.6
4M-EDR−0.29<0.001−0.29<0.001−0.30<0.001
25(OH)D levels- −0.050.2−0.050.3
RBC folate levels- −0.040.3−0.050.2
Italics and bold indicate results that are statistically significant. Adjustments: Model 1 = 25(OH)D and RBC folate levels. Model 2 = Model 1 and sex, age, creatinine and vitamin B12 levels, reported dietary intake of alcohol, vitamin B6, tea and coffee, smoking status and BMI category. Total number of participants in each model vary due to missing data.
Table
Table 3. Assessment of interactions between EDR measures (6W-EDR and 4M-EDR) and levels of 25(OH)D or RBC folate in determining Hcy levels, with and without adjustments.
Table 3. Assessment of interactions between EDR measures (6W-EDR and 4M-EDR) and levels of 25(OH)D or RBC folate in determining Hcy levels, with and without adjustments.
Hcy Levels
UnadjustedAdjusted
βpβp
6W-EDR−0.24<0.001−0.24<0.001
25(OH)D levels−0.04 0.3−0.03 0.5
6W-EDR x 25(OH)D levels0.100.010.150.002
n = 582/467
4M-EDR−0.27<0.001−0.27 <0.001
25(OH)D levels−0.05 0.2−0.03 0.4
4M-EDR x 25(OH)D levels0.04 0.30.09 0.05
n = 582/467
6W-EDR−0.24 <0.001−0.24 <0.001
RBC folate levels0.00 0.9−0.01 0.9
6W-EDR x RBC folate levels0.01 0.90.01 0.9
n = 612/491
4M-EDR−0.29 <0.001−0.29 <0.001
RBC folate levels−0.02 0.6−0.04 0.4
4M-EDR x RBC folate levels0.01 0.90.01 0.9
n = 612/491
Italics and bold indicate results that are statistically significant. P-values for interactions were compared against Bonferroni-adjusted p thresholds of p < 0.008 for 6W-EDR models and p < 0.0125 for 4M-EDR models, to account for multiple testing. Tabled variables were entered into models as continuous variables. Adjustments: sex, age, creatinine and vitamin B12 levels, reported dietary intake of alcohol, vitamin B6, tea and coffee, smoking status, and BMI category. Total number of participants in each model varies due to missing data.
Table
Table 4. Association between 6W-EDR and Hcy levels stratified by quartiles of 25(OH)D levels.
Table 4. Association between 6W-EDR and Hcy levels stratified by quartiles of 25(OH)D levels.
Hcy Levels
6W-EDR
by 25(OH)D quartiles (ng/mL)nβpmean95% CI
Q1 (<23.95)
Mean: 18.60		107−0.350.00110.68.6–12.6
Q2 (23.96–33.25)
Mean: 28.24		120−0.320.00110.28.2–12.1
Q3 (33.26–45.60)
Mean: 39.51		108−0.180.099.87.9–11.7
Q4 (>45.60)
Mean: 56.88		1320.010.910.28.3–12.2
6W-EDR x 25(OH)D quartiles (as categorical)—pinteraction = 0.004
Italics and bold indicate results that are statistically significant. The p-value for interaction between 6W-EDR and 25(OH)D was compared against a Bonferroni-adjusted p threshold of p < 0.008 to account for multiple testing between 6W-EDR and variables of interest. Values shown are adjusted with sex, age, and Hcy determinants, creatinine and vitamin B12 levels, reported dietary intake of alcohol, vitamin B6, tea and coffee, smoking status, and BMI category.
Table
Table 5. Assessment of gene–environment interactions between DHFR-rs70991108 and MTHFD1-rs2236225 variants and 2D-EDR or 6W-EDR in determining Hcy levels.
Table 5. Assessment of gene–environment interactions between DHFR-rs70991108 and MTHFD1-rs2236225 variants and 2D-EDR or 6W-EDR in determining Hcy levels.
Hcy Levels
UnadjustedAdjusted
βpβp
6W-EDR−0.23<0.001−0.24<0.001
DHFR-rs70991108−0.060.1−0.030.6
6W-EDR x DHFR-rs709911080.020.7−0.010.9
n = 611/462
6W-EDR−0.20<0.001−0.170.001
MTHFD1-rs2236225−0.050.2−0.090.05
6W-EDR x MTHFD1-rs22362250.060.20.110.05
n = 611/461
4M-EDR−0.28<0.001−0.30<0.001
DHFR-rs70991108−0.060.1−0.030.5
4M-EDR x DHFR-rs709911080.020.7−0.010.9
n = 611/462
4M-EDR−0.23<0.001−0.22<0.001
MTHFD1-rs2236225−0.050.2−0.090.04
4M-EDR x MTHFD1-rs22362250.100.030.150.006
n = 611/461
Italics and bold indicate results that are statistically significant. P-values for interactions were compared against Bonferroni-adjusted p thresholds of p < 0.008 for 6W-EDR models and p < 0.0125 for 4M-EDR models, to account for multiple testing. Adjustments: RBC folate levels and Hcy determinants, sex, age, creatinine and vitamin B12 levels, reported dietary intake of alcohol, vitamin B6, tea and coffee, smoking status and BMI category. Totals shown are for unadjusted and adjusted models respectively. Total number of participants in each model varies due to missing data.
Table
Table 6. Independent and interactive influences of EDR measures, 25(OH)D levels and MTHFD1-rs2236225 on Hcy levels.
Table 6. Independent and interactive influences of EDR measures, 25(OH)D levels and MTHFD1-rs2236225 on Hcy levels.
Hcy Levels
Unadjusted (n = 578)Adjusted (n = 461)
βpβp
6W-EDR−0.20<0.001−0.190.001
25(OH)D levels−0.040.3−0.030.6
6W-EDR x 25(OH)D levels0.100.020.140.002
MTHFD1-rs2236225−0.040.3−0.090.05
6W-EDR x MTHFD1-rs22362250.060.20.110.05
Model—R2 (p)0.07<0.0010.09<0.001
4M-EDR−0.22<0.001−0.21<0.001
25(OH)D levels−0.050.2−0.030.4
4M-EDR x 25(OH)D levels0.030.40.090.06
MTHFD1-rs2236225−0.040.3−0.090.04
4M-EDR x MTHFD1-rs22362250.100.030.140.006
Model—R2 (p)0.08<0.0010.11<0.001
Italics and bold indicate results that are statistically significant. Adjustments: RBC folate levels and Hcy determinants, sex, age, creatinine and vitamin B12 levels, reported dietary intake of alcohol, vitamin B6, tea and coffee, smoking status and BMI category. Totals shown are for unadjusted and adjusted models, respectively. The total number of participants in each model varies due to missing data.

Share and Cite

MDPI and ACS Style

Jones, P.; Lucock, M.; Martin, C.; Thota, R.; Garg, M.; Yates, Z.; Scarlett, C.J.; Veysey, M.; Beckett, E. Independent and Interactive Influences of Environmental UVR, Vitamin D Levels, and Folate Variant MTHFD1-rs2236225 on Homocysteine Levels. Nutrients 2020, 12, 1455. https://0-doi-org.brum.beds.ac.uk/10.3390/nu12051455

AMA Style

Jones P, Lucock M, Martin C, Thota R, Garg M, Yates Z, Scarlett CJ, Veysey M, Beckett E. Independent and Interactive Influences of Environmental UVR, Vitamin D Levels, and Folate Variant MTHFD1-rs2236225 on Homocysteine Levels. Nutrients. 2020; 12(5):1455. https://0-doi-org.brum.beds.ac.uk/10.3390/nu12051455

Chicago/Turabian Style

Jones, Patrice, Mark Lucock, Charlotte Martin, Rohith Thota, Manohar Garg, Zoe Yates, Christopher J. Scarlett, Martin Veysey, and Emma Beckett. 2020. "Independent and Interactive Influences of Environmental UVR, Vitamin D Levels, and Folate Variant MTHFD1-rs2236225 on Homocysteine Levels" Nutrients 12, no. 5: 1455. https://0-doi-org.brum.beds.ac.uk/10.3390/nu12051455

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