- 1Institute for Medical Research and Occupational Health, Zagreb, Croatia
- 2Faculty of Food Technology and Biotechnology, University of Zagreb, Zagreb, Croatia
- 3School of Medicine, University of Zagreb, Zagreb, Croatia
- 4University Hospital Sveti Duh, Zagreb, Croatia
- 5University Hospital Center Rijeka, Rijeka, Croatia
Background: Previous research has identified that, in individuals with lower folate and/or vitamin B12 levels, homocysteine is associated with specific elements, while amino acids are associated with folate, thereby influencing folate availability during supplementation. This study investigates the impact of the MTHFR 677C/T polymorphism on levels of these compounds in individuals with differing meat consumption preferences, non-vegetarians and vegetarians.
Methods: The study was conducted on 162 Croatian subjects. MTHFR gene polymorphism was determined by PCR-RFLP, elements by ICP-QQQ, amino acids by GC–MS, and biochemical parameters by chemiluminescent immunoassay and enzymatic methods. Differences between the groups were tested by ANOVA, while associations between the parameters were examined by multiple regression analyses.
Results: The CT + TT genotype had lower folate levels in comparison to the CC genotype. The influence of MTHFR 677C/T polymorphism on element levels was limited to Al and Ca: individuals with the CT + TT genotype exhibited higher plasma Ca levels than the wild-type (CC) genotype, while lower plasma Al values were observed only in vegetarians with the CT + TT genotype. The CT + TT genotype was also associated with higher levels of glutamic acid, proline, glycine (only in vegetarians), and aspartic acid.
Conclusion: Although we observed an effect of the MTHFR 677C/T polymorphism on folate, certain elements and amino acid levels, further research is required to validate our findings and establish a more comprehensive understanding of the link between elements, amino acids and the MTHFR gene polymorphism.
1 Introduction
The enzyme methylenetetrahydrofolate reductase (MTHFR) plays a crucial role in the metabolism of folate and one-carbon transfer reactions, as well as in regulating homocysteine levels by reducing 5,10-methylenetetrahydrofolate to 5-methyltetrahydrofolate (5-MTHF), which serves as a donor of the methyl group for homocysteine to methionine conversion. These processes are crucial for DNA methylation. Several polymorphic variants have been identified on MTHFR gene, and the most extensively researched variant, located on chromosome 1 (1p36.3), is rs1801133, also known as 677C/T. This mutation leads to a reduction in the MTHFR enzyme activity (1–3), resulting in impaired folate metabolism and an accumulation of homocysteine in the body, causing hyperhomocysteinemia. This condition has been associated with various health issues, including cardiovascular disease, atherosclerosis, and Alzheimer’s disease (2, 4–7).
The influence of the MTHFR 677C/T polymorphism on metabolic pathways is complex and can be modulated by a variety of environmental factors, including diet. Adequate consumption of nutrients involved in one-carbon metabolism, particularly folate, vitamin B12, and vitamin B6, is essential to compensate for reduced enzyme activity in individuals carrying the polymorphism (8–10). The availability of amino acids such as methionine, serine and glycine, which are also part of the one-carbon metabolism, can also be affected by dietary choices (11, 12). Given their different dietary profiles, vegetarianism and non-vegetarianism represent two distinct dietary approaches that require special consideration when analyzing these interactions and their potential consequences on the metabolic implications of the MTHFR 677C/T variation. Vegetarian diets are usually rich in plant-based foods, providing abundant amounts of vitamins B1, C, E, and folate, along with fibre, antioxidants, polyunsaturated fatty acids and phytochemicals (13–16). However, based on the literature data and our earlier studies, vegetarian diet tends to be lower in proteins, essential amino acids, essential elements such as zinc (Zn), iron (Fe), selenium (Se), iodine (I), and calcium (Ca), as well as vitamins A, D3, and particularly B12 (17–22). These differences in diet are important, because individuals carrying the MTHFR 677C/T mutation may have different metabolic responses depending on their dietary habits, which can have a significant impact on metabolic function and overall metabolic health.
Recent studies, including data from the National Health and Nutrition Examination Survey (NHANES), suggest an association between metals, homocysteine, and folate (23, 24). Toxic elements, particularly lead (Pb), cadmium (Cd) and mercury (Hg), can interfere with homocysteine and folate metabolism by disrupting the normal function of enzymes involved in the one-carbon cycle, a pathway that relies on folate and vitamin B12 to regulate homocysteine levels (25). Folate deficiency may heighten the susceptibility of individuals to adverse effects of toxic elements (24, 26–28). In addition, in individuals with reduced folate, Pb was positively associated with homocysteine levels (26, 27), while in boys with lower B12 and folate, Hg was negatively correlated with homocysteine levels (29). Other studies have highlighted associations of amino acids with elements and folate, such as their beneficial effects on Zn absorption (30, 31) and the solubility of folic supplements (32). Folate deficiency may lead to higher histidine excretion which can be prevented by folic acid supplementation (33). Conversely, high-dose of folate supplementation can mask anemia, worsen neurological symptoms, and increase cognitive impairment (34). Moreover, vitamin B12 deficiency can interfere with methylation-dependent processes, given that methionine synthase uses 5-methyltetrahydrofolate as a one-carbon donor (35).
Prior research examining the influence of polymorphism on health has predominantly concentrated on the specific interaction between polymorphism and individual metabolic components, primarily lipids and folate (36–42), and less frequently amino acids (43–45) or particular essential/toxic substances (46–49). To date, no study has systematically explored the intricate interplay between genetic polymorphisms and the wide range of metabolic markers within a single, integrative framework.
In light of the established relationships between elements, amino acids, and folate, this research aims to investigate the interaction between the MTHFR 677C/T polymorphism and dietary habits (vegetarian vs. non-vegetarian), focusing on their combined influence on key biochemical markers. Specifically, it focuses on amino acids that are either essential–must be obtained from diet–or involved in one-carbon metabolism; essential and toxic elements that are influenced by dietary intake and can be reliably measured in blood and/or plasma; and biochemical parameters previously associated with trace elements or MTHFR genotype. This integrative approach could provide new insights into the role of personalized nutrition in managing the health risks associated with genetic variations in the MTHFR gene.
2 Materials and methods
2.1 Study population and sample collection
This cross-sectional study included 162 healthy adult participants recruited in Croatia in 2020 and 2021. Each participant was required to sign an informed consent form after they were fully informed about the purpose and methodology of the study. The confidentiality of participants was preserved through the use of unique codes assigned to each sample and study data. The Ethics committee of the Institute for Medical Research and Occupational Health, Zagreb reviewed and approved the study protocols, ensuring compliance with the principles of the Helsinki Declaration. A questionnaire was used to collect general information about age, sex, education, height, weight and socioeconomic status, as well as dietary preferences, smoking status, and other lifestyle habits. The inclusion criteria in this study were as follows: adult vegetarians and non-vegetarians, no clinical signs of acute infection, and for women, an absence of pregnancy or breastfeeding during the study period. The study participants were divided into two groups based on their dietary habits: the vegetarian group, consisting of subjects who were vegan or lacto/ovo-vegetarian (17), and the non-vegetarian group, which included omnivores and pescatarians (those who consumed fish and/or seafood). Pescatarians are classified as non-vegetarians since they consume animal meat, such as fish and seafood, and hence cannot be considered vegetarians (17, 50). Plant-based diets inside the green box should be considered vegetarian, since they avoid meat and seafood.
Blood samples from each participant were collected between 7 and 10 a.m. using special vacutainer tubes (Becton-Dickinson, New Jersey, USA; containing K2-EDTA anticoagulant) for trace element analysis. The samples were then aliquoted into clean micro-tubes (Eppendorf, Hamburg, Germany) and stored at +4 °C for elemental analysis, or at −80 °C for DNA isolation. The remaining blood samples were centrifuged for 20 min at 1100 g (Eppendorf, Hamburg, Germany) to separate the plasma, which was then transferred into clean micro-tubes (Eppendorf, Hamburg, Germany), and stored at −80 °C until the element, amino acid, and biochemical parameters analysis. This is especially important in terms of amino acid analysis, because it is well known that pre-analytical factors, such as freeze–thaw cycles, can affect the stability and concentration of amino acids in plasma (51).
2.2 Element analysis
The inductively coupled plasma–mass spectrometry (ICP-MS) on an Agilent 8800 (Agilent Technologies, Santa Clara, CA, USA) was used for the analysis of toxic and essential elements in blood/plasma samples. Prior to element analysis, the blood and plasma samples were diluted (1:70 for blood and 1:20 for plasma) using a solution that contained 3 μg L−1 of internal standards, 0.01 mM EDTA, 0.7 mM NH3, and 0.07% (v/v) Triton X-100 (52). The working conditions of the ICP-MS, which were optimized with a tune solution of 1 μg L−1 7Li, 59Co, 89Y, 140Ce, and 205Tl, are detailed in Supplementary Table S1. Sample preparation and analysis were done in the laboratory with a Heating, Ventilating, and Air Conditioning (HVAC) system combined with HEPA filters. To identify potential contamination during the sampling, aliquoting, and dilution of blood and plasma samples, blank samples were also analyzed alongside the real blood and plasma samples. The accuracy of measurement, which ranged from 93 to 113% (Supplementary Tables S2, S3), was verified through the use of commercially available reference materials [Seronorm™ Trace Elements Whole Blood Level I, II and III and Serum Level I and II (Sero AS, Billingstad, Norway), ClinChek® Whole Blood Control Level I, II and III, Serum Level I and II, and Plasma Level I and II (Recipe, Munich, Germany)] and regular participation in the Interlaboratory Comparison Programme UK NEQAS (Birmingham, UK) for Trace Elements in blood and serum.
2.3 Biochemical parameters analysis
Biochemical parameters (glucose, triglycerides, total cholesterol, HDL cholesterol, LDL cholesterol, folate, and B12 vitamin) in plasma samples were measured by standardized and harmonized methods (Supplementary Table S4) on the Atellica Solution Analyser (Siemens, Erlangen, Germany). The reagents and calibrators were provided by the analyser manufacturer (Siemens, Erlangen, Germany), while quality control samples were obtained from Bio-Rad Laboratories (Hercules, California, USA) for glucose, triglycerides, and total, HDL, and LDL cholesterol, and from Sero (Billingstad, Norway) for B12 vitamin and folates. Verification of biochemical parameters analyzed on the Atelica Solution analyzer was performed according to protocols that are compliant with Clinical and Laboratory Standards Institute (CLSI) guidelines and daily imprecision was followed using commercial quality control samples.
2.4 Amino acids analysis
Concentrations of 25 amino acids (alanine, sarcosine, glycine, α-aminobutyric acid, valine, β-aminoisobutyric acid, leucine, allo-iso-leucine, iso-leucine, threonine, serine, proline, asparagine, aspartic acid, methionine, 4-hydroxyproline, phenylalanine, glutamic acid, α-aminoadipic acid, cysteine, ornithine, lysine, histidine, tyrosine, tryptophan) were measured in the plasma samples by Trace 1,300 gas chromatograph (Thermo Scientific, Milan, Italy) coupled to an ITQ 700 ion trap mass spectrometer (Thermo Scientific, Austin, TX, USA). Before GC–MS analysis, samples were processed using the Phenomenex EZ:faast™ amino acid analysis kit as described in Badawy et al. (53). Norvaline was used as an internal standard. Chromatographic separation was performed on a ZB-AAA GC column (10 m × 0.25 mm ID; 0.25 μm film thickness). The injection port was operated in split mode (1:10) at 250 °C. The carrier gas (helium) flow was constant at 1.1 mL/min. The oven temperature program was set at 30 °C/min starting from 100 °C until reaching the final temperature of 320 °C. The operating conditions for the MS system were as follows: electron ionization mode at the energy of 70 eV, transfer line temperature 180 °C, and ion source temperature 240 °C. Full scan mode with a scan range of m/z 45–450 was used for data acquisition.
2.5 DNA isolation and genotyping of MTHFR 677C/T polymorphisms
DNA was isolated from blood using QIAamp DNA Blood Mini Kit (Qiagen, Hilden, Germany), while BioSpec-nano (Shimadzu, Japan) was used to check the concentration and purity of isolated DNA. The amount of DNA in the eluate, measured by absorbance at 260 and 280 nm, is used to determine the DNA yield. The ratio of absorbance at 260 nm to 280 nm (A260/A280) is used to assess DNA purity. The A260/A280 ratio for pure DNA was in the range from 1.7 to 1.9. After isolation, DNA was stored at −20 °C until polymorphism analysis.
Polymerase chain reaction (PCR) on a T100 thermal cycler (Bio-Rad Laboratories, Hercules, California, USA) was used to amplify the fragment of 198-bp with the following primers: 5′-TGA AGG AGA AGG TGT CTG CGG GA-3′ (forward) and 5′-AGG ACG GTG CGG TGA GAG TG-3′ (reverse). Amplification was conducted under the following conditions: initial denaturation at 95 °C for 5 min, followed by 35 cycles of denaturation at 95 °C for 30 s, annealing at 60 °C for 30 s, and extension at 72 °C for 30 s, with a final extension at 72 °C for 5 min. PCR products were digested with HinfI (New England BioLabs, Hitchin, UK) as a restriction enzyme (39, 54) at 37 °C for 4 h. After digestion, the fragments were separated by electrophoresis on 2.5% agarose gel with Lonza GelStar strain (Lonza Group Ltd., Basel, Switzerland) and visualized under an ultraviolet illuminator. The characteristic fragments for the recessive homozygotes genotype (TT) were 175-bp and 23-bp, for the heterozygotes genotype (CT) were 198-bp, 175-bp, and 23-bp, while for the dominant homozygotes genotype (CC) were 198-bp. To confirm the quality of the analysis, 15% of the randomly selected samples were repeated.
2.6 Statistical analysis
The concordance between observed allele frequency and Hardy–Weinberg equilibrium was calculated using the χ2 test, while the differences in genotype frequencies between non-vegetarians and vegetarians were tested by Fisher’s exact test. Study subjects were divided into two main groups based on the MTHFR 677C/T polymorphism: the wild-type genotype (CC genotype, n = 73) and the CT + TT genotype (n = 87). Considering the low frequency of the TT genotype (6.9%), assessing the effects associated with this group would have led to unreliable estimates and greatly reduced statistical power. Therefore, a dominant genetic model (CT + TT vs. CC) was used, acknowledging that TT-specific effects might be obscured. Participants were further stratified based on their dietary habits into non-vegetarians (CC genotype n = 43; CT + TT genotype n = 49) and vegetarians (CC genotype n = 30; CT + TT genotype n = 38).
Differences between groups based on dietary habits (non-vegetarians vs. vegetarians) and genotypes (CT + TT vs. CC genotype) were tested by the analysis of variance (ANOVA) and Duncan’s post-hoc test. Additionally, the associations between the studied parameters and dietary habit and genotype, adjusted for age, sex, smoking, BMI, education, and intake of dietary supplements, were analysed using multiple regression analyses. Key demographic parameters included age and gender; prospective lifestyle factors involved BMI and smoking status; education level served as a proxy for socioeconomic status; and dietary supplement intake was considered to account for potential side effects. Variables like smoking, sex, education, and supplement intake were considered as categorical. Sex was categorized as male/female, smoking as yes/no, supplement intake was specified by type of supplement (folic acid, vitamin B12, minerals, and amino acids), while education was divided into two levels: (1) secondary school vs. (2) university degree and post-graduate degree (PhD). Initial models also included the effect of interaction between dietary habit and genotype. In cases when this interaction was not statistically significant, it was omitted from the final model.
For variables with <10% of values below the limit of detection (LOD), values below LOD were set to 0.5 LOD for the subsequent analyses. In cases when more than 10% of values were below LOD (As in blood (63%), Al (41%), α-aminoadipic acid (12%), and cysteine (38%) in plasma), differences between groups were tested by ANOVA for censored data, followed by Tukey’s multiple comparison test (55). In the case of cysteine and α-aminoadipic acid, differences were tested by the Kruskal–Wallis test (values <LOD are considered as ties) as censored data ANOVA is not recommended in cases when the distribution of residuals deviates from lognormal distribution. For multiple regression analyses, censored likelihood multiple imputation was employed (56). Most of the dependent variables were log-transformed to approach the normal distribution of residuals in ANOVA and regression analysis. The relationships between measured parameters were analysed using Spearman’s correlation coefficients, calculated separately for each combination of genotype (CC and CT + TT genotype) and dietary habit (non-vegetarian and vegetarian). For meaningful interpretation, parameters were grouped based on their nutritional characteristics and biological functions as follows: B vitamins (folates and vitamin B12); lipid profile (triglycerides, total cholesterol, HDL, and LDL cholesterol); toxic elements (As, Cd, Hg, Pb, Al, Sn, Tl, Cs); essential elements (Ca, Cu, Fe, I, Mg, Se, Zn, Ba, Sr); trace essential elements (Mo, Mn, Co, Cr); essential amino acids (valine, leucine, isoleucine, threonine, methionine, phenylalanine, lysine, histidine, tryptophan); non-essential amino acids (alanine, glycine, serine, proline, asparagine, aspartic acid, 4-hydroxyproline, glutamic acid, cysteine, ornithine, tyrosine); and non-proteinogenic amino acids (sarcosine, α-aminobutyric acid, β-aminoisobutyric acid, allo-iso-leucine, α-aminoadipic acid). Further examining the differences between different types of diets, sparse discriminant analysis (57) was performed. Sparse discriminant analysis (SDA) is a variant of linear discriminant analysis that attempts to find sparse linear combinations of variables that discriminate between different subgroups best. It is especially suitable for limited sample sizes, including situations when the number of variables is larger than the sample size. Furthermore, compared to linear discriminant analysis, the sparseness criterion reduces the number of variables included in each discriminative vector, allowing easier interpretation of results. Differences between three types of dietary habits with respect to measured elements in blood and plasma and amino acids were sought: (1) omnivores (2) lacto-ovo vegetarians, and (3) vegans, while the smallest subgroup (pescatarians; N = 10) was excluded from the analysis. The number of variables per each discriminative vector was selected by leave-one-out cross-validation. Before the analysis, all variables were log-transformed and normalized to have equal mean and variance. Variables with more than 10% of values below the LOD (i.e., blood As, plasma Al, α-aminoadipic acid, and cysteine) were excluded from the analysis. After selecting the number of variables per discriminative vector by leave-one-out cross-validation, the model was fitted and evaluated on the whole dataset. Thus, the classification accuracy for the whole dataset would correspond to training set accuracy. Post-hoc power for detecting gene-diet interactions was calculated based on the F-test for ANOVA model with subgroup sizes of 30 (minimum subgroup size in our sample) and 40 (average subgroup size). Power was calculated for the effect sizes of 0.1, 0.25, and 0.4, corresponding to small, medium, and large effect size, respectively (58). The level of significance was set to 0.05. Power was calculated in G*Power v. 3.1.9.7 (59). Power for a medium effect size of 0.25 was 0.78 for a sample size of 30 per subgroup, and 0.88 for a sample size of 40 per subgroup. For a large effect size of 0.40, the power was larger than 0.99 for both scenarios. On the other hand, the power for detecting interactions with small effect sizes was low (0.19 and 0.24 for the two scenarios).
Statistical analysis was conducted using the software package TIBCO Statistica™, version 14.0.0.15 (TIBCO Software, Inc., USA), and the statistical software R, version 4.3.3 (R Foundation for Statistical Computing, Austria). Multiple regression and ANOVA for censored data were performed in R, as well as sparse discriminant analysis. Package lodi was used for censored likelihood multiple imputation in multiple regression analysis, package NADA2 for censored data ANOVA, while package sparseLDA was used for discriminant analysis. Any results with a p-value below 0.05 were considered significant.
3 Results
3.1 General characteristics of the study group
Table 1 displays the general characteristics [median and interquartile range (IQR)] of the study cohort, as well as the MTHFR 677C/T frequency information, while flowchart of participants is presented in Supplementary Figure S1. The median age of the study group was 35 years (range: 18–65 years), and 48.1% of the participants held a university degree, while 32.7% obtained a PhD. Three vegetarian participants did not disclose their educational background. In the study group, 28.4% of subjects were male, while 26.2% smoked cigarettes, with median value of 7 (4–11) cigarettes per day. In terms of dietary habits, 44.3% of vegetarian participants identified as vegans, while 55.7% identified as lacto-ovo vegetarians. They adhered to this diet for 9 years (4–21), with their primary motivations being animal welfare (86% of participants), planet sustainability (74%), health concerns (63%), and to a lesser extent, meat aversion (30%). Out of the participants in the non-vegetarian group, a small percentage (10.9%) followed a pescatarian diet, while the majority (89.1%) were omnivores. Vegetarians had significantly lower (p < 0.001) BMI [21.6 (20.2–23.7)] compared to those who consume meat [23.2 (21.7–24.9)]. However, they consumed significantly higher (p = 0.011) amounts of supplements (78.6%) in comparison to non-vegetarians (58.7%). Out of these 109 participants taking any type of dietary supplementation, 93 participants were taking vitamin B12, minerals, or folic acid, with 13 participants taking all three types of supplements and 27 participants taking two of these three supplements (N = 25 for simultaneous intake of minerals and B12 and N = 2 for minerals and folic acid). No significant differences were found in the allele frequencies between non-vegetarians and vegetarians (p = 0.730). Vegetarians exhibited a frequency prevalence of MTHFR 677C/T of 44.1% homozygous dominant (CC), 47.1% heterozygous (CT), and 8.8% homozygous recessive (TT), compared to non-vegetarians who had a frequency prevalence of 46.7% CC, 47.8% CT, and 5.5% TT. Allele frequencies were in accordance with the Hardy–Weinberg equilibrium for all participants combined (χ2 = 2.22, p = 0.136), as well as for non-vegetarians (χ2 = 2.16, p = 0.142) and vegetarians (χ2 = 0.384, p = 0.536) considered separately.
3.2 The influence of gene-dietary interaction on the levels of toxic and essential elements
Table 2 shows element concentrations in whole blood and plasma of the non-vegetarians and vegetarians with regard to the MTHFR 677C/T polymorphism obtained by ANOVA. Results showed that vegetarians had significantly higher levels of Cd and Co in blood, as well as Sr and Mo in plasma, than non-vegetarians. On the other hand, the levels of several elements in their blood (Cs, Hg, and As) and plasma (Sn, Se, Zn, I, Ca, and Cu) were significantly lower compared to non-vegetarians. Plasma Ca levels were increased in participants with the CT + TT genotype. Additionally, plasma Al levels of vegetarians with the CT + TT genotype were significantly lower than those with the CC genotype (2.00 (1.54–2.96) μg L−1 vs. 2.68 (1.78–3.91) μg L−1, respectively), whereas the plasma Al levels in non-vegetarians with the CT + TT genotype (2.92 (2.32–3.59) μg L−1) were comparable to those with the CC genotype (3.08 (2.53–3.43) μg L−1). To better understand the associations between the studied parameters and dietary habits and genotype adjusted for potential confounders, we conducted a multiple regression analysis (Figure 1; Supplementary Table S5). While ANOVA analysis indicated a significant interaction between dietary habits and genotype for blood Mn, this effect was not confirmed in multiple regression analysis adjusted for age, sex, smoking, BMI, education, and supplement intake. Interaction between genotype and dietary habits remained statistically significant for Al in plasma (β = −0.34 [−0.62, −0.06], p = 0.019). Furthermore, multiple regression analysis showed that participants with the CT + TT genotype among both non-vegetarians and vegetarians exhibited a lower blood As (β = −0.47 [−0.84, −0.09], p = 0.017) and an increased plasma Ca (β = 0.01 [0.002, 0.03], p = 0.015) compared to those with the CC genotype (Supplementary Table S5). All differences in dietary habits (vegetarians vs. non-vegetarians) observed in the ANOVA analysis, except for plasma Ca, remained significant in the multiple regression analysis, adjusted for age, sex, smoking, BMI, education, and supplement intake, across all measured elements in blood and plasma.
Table 2. The concentration of elements in whole blood and plasma of non-vegetarians and vegetarians with regard to the MTHFR 677C/T gene polymorphism.
Figure 1. Associations between dietary habits and elements in blood and plasma, represented as regression coefficients and 95% confidence intervals for each element, obtained using multiple linear regression (reference group: non-vegetarians) after adjusting for age, sex, genotype, smoking status, body mass index, level of education, and intake of dietary supplements.
In addition, results of the multiple linear regression analysis (Supplementary Table S5) revealed that the concentrations of Pb, Cs, and As in blood, as well as Cu and Sr in plasma increased with age, with p-values ranging from 0.007 to 0.033. The blood Pb (β = 0.29 [0.08, 0.51], p = 0.007) and plasma Fe (β = 0.43 [0.27, 0.58], p < 0.001), Mo (β = 0.16 [0.008, 0.31], p = 0.039), Se (β = 0.09 [0.03, 0.14], p = 0.003), Ca (β = 0.03 [0.01, 0.04], p < 0.001), and Sn (β = 0.15 [0.02, 0.29], p = 0.027) levels in male participants were significantly higher than those in female participants, whereas their blood levels of Mn (β = −0.12 [−0.23, −0.02], p = 0.023) and Co (β = −0.41 [−0.57, −0.24], p < 0.001) were significantly lower. Smoking was associated with lower plasma Cu (β = −0.07 [−0.13, −0.01], p = 0.023) and higher blood Cd (β = 0.90 [0.68, 1.10], p < 0.001) levels. Participants with increased BMI had higher blood Cs (β = 0.02 [0.002, 0.04], p = 0.028), while those with higher education levels exhibited lower Cd (β = −0.25 [−0.49, −0.005], p = 0.045) and higher Cs (β = 0.14 [0.004, 0.28], p = 0.043) in the blood. The only two elements in plasma that exhibited a negative correlation with supplement intake were Sn (β = −0.14 [−0.26, −0.02], p = 0.021) and Al (β = −0.13 [−0.25, −0.001], p = 0.048).
3.3 The influence of gene-dietary interaction on biochemical parameters
The concentrations of biochemical parameters in plasma (total, HDL, and LDL cholesterol, glucose, triglycerides, folate, and vitamin B12) in non-vegetarians and vegetarians in relation to the MTHFR 677C/T polymorphism are illustrated in Figures 2, 3. The results of ANOVA (Figures 2, 3) showed statistically significant differences in total cholesterol (p = 0.0002) and vitamin B12 (p = 0.012) levels between vegetarians and non-vegetarians in the CT + TT genotype group, with lower levels observed in vegetarians. Reduced levels of LDL cholesterol were found in vegetarians in both genotype groups compared to non-vegetarians (Figure 2). Furthermore, non-vegetarians with the CT + TT genotype had higher levels of LDL cholesterol (p = 0.049) and lower levels of folate (p = 0.039) than non-vegetarians with the CC genotype. The ANOVA analysis also indicated a statistically significant interaction between diet and genotype only for LDL cholesterol (F(1,156) = 5.64, p = 0.019), suggesting that non-vegetarian CT + TT genotype had higher LDL cholesterol levels (1.97 (1.54–2.32) mmol L−1) than the CC genotype (1.63 (1.4–2.06) mmol L−1), whereas LDL cholesterol levels in vegetarians were nearly identical for both genotypes (1.37 (0.99–1.62) mmol L−1 for CT + TT genotype and 1.35 (1.18–1.65) mmol L−1 for the CC genotype). On the other hand, the results of multiple regression analysis (Supplementary Table S6), adjusted for age, sex, genotype, smoking, BMI, education, and supplement intake, showed a significant genotype-diet interaction only for HDL cholesterol (β = 0.16 [0.003, 0.33], p = 0.046). Findings from the multiple regression analysis confirm associations between genotype and folate (CT + TT genotype had lower folate levels than CC genotype), although a significant association with LDL cholesterol was not confirmed. Significant associations between diet and total and LDL cholesterol, as well as vitamin B12 were confirmed in the multiple regression analysis, while also revealing a negative association with HDL cholesterol (β = −0.14 [−0.25, −0.01], p = 0.028).
Figure 2. Concentrations of total cholesterol, HDL, and LDL cholesterol in plasma of non-vegetarians (CC genotype n = 43; CT + TT genotype n = 49) and vegetarians (CC genotype n = 30; CT + TT genotype n = 38) in relation to MTHFR 677C/T polymorphism (CC genotype n = 73; CT + TT genotype n = 87), presented by box and whisker plots. Only significant differences between genotypes (tested separately within each diet group) and diet (tested separately within each genotype) obtained by ANOVA post-hoc comparisons are indicated.
Figure 3. Concentrations of vitamin B12, folate, glucose, and triglycerides in plasma of non-vegetarians (CC genotype n = 43; CT + TT genotype n = 49) and vegetarians (CC genotype n = 30; CT + TT genotype n = 38) in relation to MTHFR 677C/T polymorphism (CC genotype n = 73; CT + TT genotype n = 87), presented by box and whisker plots. Only significant differences between genotypes (tested separately within each diet group) and diet (tested separately within each genotype) obtained by ANOVA post-hoc comparisons are indicated.
Furthermore, the findings of the multiple regression analysis showed that total and LDL cholesterol increased with age (β = 0.01 [0.003, 0.009], p = 0.001 for total cholesterol; β = 0.01 [0.004, 0.01], p = 0.001 for LDL cholesterol), and BMI (β = 0.01 [0.004, 0.02], p = 0.006 for total cholesterol; β = 0.03 [0.01, 0.04], p < 0.001 for LDL cholesterol). Increased levels of triglycerides were also associated with higher BMI (β = 0.03 [0.009, 0.04], p = 0.003), age (β = 0.006 [0.0002, 0.01], p = 0.043), and smoking (β = 0.18 [0.05, 0.31], p = 0.007), while glucose was associated only with BMI (β = 0.009 [0.005, 0.01], p < 0.001). Additionally, vitamin B12 was positively associated with supplementary intake (β = 0.22 [0.10, 0.34], p = 0.001).
3.4 The influence of gene-dietary interactions on amino acid levels
Table 3 shows the concentrations of amino acids in the plasma of vegetarians and non-vegetarians with regard to the MTHFR polymorphism. The results of ANOVA (Table 3) showed that vegetarians had reduced levels of sarcosine, α-aminobutyric acid, lysine, leucine, 4-hydroxyproline, and higher levels of alanine, proline, glycine, ornithine, allo-iso-leucine, asparagine, and glutamic acid in comparison to non-vegetarians. Furthermore, CT + TT genotype had significantly higher levels of proline and glutamic acid than CC genotype. Results of ANOVA (Table 3) also indicated a significant interaction between dietary habits and genotype for asparagine (F(1,155) = 4.11, p = 0.044), glycine (F(1,155) = 4.59, p = 0.034), and alanine (F(1,155) = 4.91, p = 0.028), with latter two further confirmed in regression analysis (interaction β = 0.16 [0.005, 0.32], p = 0.043 and β = 0.21 [0.04, 0.37], p = 0.013, for glycine and alanine, respectively; Supplementary Table S7). Precisely, plasma levels of alanine, asparagine, and glycine in the CT + TT genotype were higher in vegetarians than in non-vegetarians, while in the CC genotype these differences were not observed. Additionally, multiple regression analysis adjusted for age, sex, smoking, BMI, education, and supplement intake (Figure 4; Supplementary Table S7) confirmed significant effects of dietary habits on sarcosine, α-aminobutyric acid, allo-iso-leucine, proline, asparagine, 4-hydroxyproline, glutamic acid, ornithine, and lysine, but not on leucine. Although both analyses indicated lower values of 4-hydroxyproline among vegetarians, residuals from the parametric regression model deviated significantly from a normal distribution. Multiple regression analysis confirmed the effects of genotype on glutamic acid and proline, and additionally revealed increased aspartic acid among participants with CT + TT genotype (β = 0.20 [0.02, 0.38], p = 0.027).
Table 3. The concentration of amino acids in plasma of non-vegetarians and vegetarians with regard to the MTHFR 677C/T gene polymorphism.
Figure 4. Associations between dietary habits and amino acid concentrations in plasma, represented as regression coefficients and 95% confidence intervals for each amino acid, obtained using multiple linear regression analysis (reference group: non-vegetarians) after adjusting for age, sex, genotype, smoking status, body mass index, level of education, and intake of dietary supplements.
Furthermore, regression analysis indicated that the sex of the participants had a notable impact on the levels of measured amino acids. Thus, male participants exhibited higher levels of over half of the measured amino acids (proline, methionine, phenylalanine, glutamic acid, α-aminoadipic acid, valine, leucine, allo-iso-leucine, iso-leucine, β-aminoisobutyric acid, ornithine, lysine, tyrosine, and tryptophan) than their female counterparts, with p-values ranging from <0.001 to 0.029 (Supplementary Table S7). Increases in serine (β = 0.14 [0.01, 0.27], p = 0.031), proline (β = 0.13 [0.02, 0.25], p = 0.020), aspartic acid (β = 0.25 [0.04, 0.45], p = 0.018), α-aminoadipic acid (β = 1.10 [0.23, 1.97], p = 0.014), glutamic acid (β = 0.17 [0.07, 0.28], p = 0.002), and ornithine (β = 0.11 [0.004, 0.22], p = 0.042) were significantly correlated with smoking. The data also revealed a negative association between age and concentrations of leucine (β = −0.47 [−0.88, −0.06], p = 0.026), while showing positive associations with glutamic acid (β = 0.006 [0.0004, 0.01], p = 0.036) and ornithine (β = 0.005 [0.0001, 0.01], p = 0.046). A positive association was observed between BMI and alanine (β = 0.02 [0.008, 0.03], p = 0.001), α-aminobutyric acid (β = 0.02 [0.004, 0.04], p = 0.019), sarcosine (β = 0.03 [0.006, 0.06], p = 0.015), and valine (β = 0.02 [0.005, 0.03], p = 0.003). Additionally, individuals with a high level of education had marginally higher levels of α-aminobutyric acid (β = 0.17 [0.007, 0.34], p = 0.041).
3.5 Results of Spearman’s correlation and sparse discriminant analysis
Results of the Spearman’s correlation analysis of all measured parameters are presented as heatmaps of correlation coefficients in Supplementary Figures S2–S5. Overall, positive correlations were more frequent and stronger than negative ones across the entire biomarker panel. The strongest associations were found within biologically related groups, whereas correlations between major biomarker classes (lipids, elements and amino acids) were generally weaker, suggesting largely independent variation across these domains. Within the lipid profile, triglycerides, total cholesterol, LDL and HDL cholesterol showed the expected positive inter-correlations, with the exception of associations of HDL cholesterol with triglycerides and LDL cholesterol. Folate and vitamin B12 were weakly related to most other markers, with a positive correlation with each other. Amino acids showed the most pronounced clustering. Essential amino acids, particularly valine, leucine, and iso-leucine, were strongly and positively inter-correlated. Non-essential amino acids also displayed positive correlations within their group (e.g., glutamic acid, aspartic acid, serine, and proline) and with essential amino acids. Tyrosine, in particular, displayed high correlations with several essential amino acids (threonine, methionine, phenylalanine, lysine, histidine, and tryptophan). Among elements, several toxic elements were significantly positively correlated (e.g., As with Hg, and Tl with Cd and Pb), whereas negative correlations were observed between selected essential and toxic elements, such as Cd with I, Mg, Se, and Cu.
Based on the cross-validation procedure, six variables per discriminant vector were selected in order to retain sparseness while keeping the classification accuracy sufficiently high (67%) in comparison to the maximum achieved accuracy (75%), which was achieved when the number of variables per discriminant vector exceeded 38. Figure 5. shows the projection of data on two resulting SDA discriminant vectors, after applying the SDA to the whole dataset. The classification accuracy for the whole dataset was 77%. The first discriminant vector (SD1), effectively discriminates between (1) omnivores and (2) lacto-ovo vegetarians and vegans. Namely, lacto-ovo vegetarians and vegans, predominantly positioned on the right side of the graph, had higher SD1 values (coefficients shown in the figure caption), suggesting that participants with higher levels of proline and asparagine and lower levels of Hg, Se, sarcosine and α-aminobutyric acid were more likely to be lacto-ovo vegetarians or vegans. With respect to SD2, the largest difference was noted between omnivores (on average highest values of SD2) and vegans (lowest values of SD2), while lacto-ovo vegetarians were positioned between these two groups (Figure 5). However, the overlap between the groups seemed to be more pronounced compared to SD1. In comparison to an omnivore diet, vegans generally had higher plasma Sr and Mo levels and lower plasma Cu, lysine, and α-aminobutyric acid levels. Regarding β-aminoisobutyric acid levels, vegans had, on average, higher values than lacto-ovo vegetarians but lower values than omnivores.
Figure 5. Outcomes of sparse discriminant analysis. The figure shows the projection of data onto the first (SD1) and second (SD2) discriminant vectors. Coefficients for SD1 were as follows: −1.56 (blood Hg), −0.65 (plasma Se), −3.16 (sarcosine), −0.46 (α-aminobutyric acid), 0.94 (proline), and 0.45 (asparagine), while coefficients for SD2 were: −0.06 (plasma Mo), 0.84 (plasma Cu), −0.67 (plasma Sr), 0.49 (α-aminobutyric acid), −0.04 (β-aminoisobutyric acid), and 0.42 (lysine).
4 Discussion
This study investigated the effects of MTHFR 677C/T gene polymorphism and dietary habits on the levels of elements and amino acids, and biochemical parameters in a Croatian adult population. According to our findings, vegetarians had better lipid profiles than non-vegetarians, as evidenced by their lower levels of total and LDL-cholesterol (usually referred as “bad cholesterol”). This suggests that vegetarians may be at a lower risk of developing atherosclerotic plaque, which results in the narrowing of blood vessels (60–62). We also found that vegetarians had lower levels of certain trace elements and minerals, such as iodine, selenium, copper, and zinc, than non-vegetarians. This is likely due to the low presence or reduced bioavailability of these nutrients in plant-based diets (63). The patterns observed for lipid profiles and trace elements are consistent with trends observed in our previous study exploring a vegetarian diet (18). Compared to non-vegetarians, vegetarians had higher levels of non-essential amino acids, such as alanine and glycine, while essential amino acid lysine was significantly lower. The variations among the study groups are consistent with the known nutritional makeup of plant-based foods, which typically have lower lysine content but higher glutamate and glycine contents (64). Furthermore, the factors identified by the SDA analysis confirmed that lacto-ovo vegetarians and vegans have higher levels of proline and asparagine, indicating greater intake of plant-based protein, which may be linked to metabolic and oxidative profiles (65, 66). Although there is limited evidence in human populations, reduced levels of sarcosine and α-aminobutyric acid may suggest different patterns of amino acid metabolism. Both metabolites have been linked to metabolic dysregulation, gut-brain communication, and oxidative stress (67, 68). At the same time, they have lower levels of blood Hg and plasma Se, indicating that they consume less fish and seafood, which reduces their exposure to Hg but also decreases Se intake, essential for antioxidant enzymes (69, 70). Vegans, in particular, tend to have lower lysine levels than omnivores, a key amino acid linked to impaired protein synthesis (71), and β-aminoisobutyric acid, which is important for lipid and glucose regulation, as well as exercise response (72). Their lower levels of plasma Cu indicate reduced bioavailability of Cu from plant-based foods (73), and elevated levels of plasma Sr and Mo likely reflect higher intake of plant-based foods and water (74, 75). In light of the observed differences in the biochemical data, element, and amino acid profiles between vegetarians and non-vegetarians, the study participants had to be categorized according to their eating preferences. This approach reduced the possibility of bias from dietary differences while enabling us to more thoroughly investigate the effects of genetic polymorphism.
4.1 Genotypic distribution and allelic frequencies
It is well established that both geographical location and ethnicity have a significant impact on the frequency of the homozygous recessive (TT) genotype of MTHFR 677C/T gene polymorphism. The prevalence of this genotype worldwide is the highest among Europeans (34.1%), followed by populations in Australia (20.5%), Asia (19.7%), South America (27.8%), North America (31.2%), Africa (10.3%), and Middle East (7.5%). On the other hand, the lowest prevalence of TT genotype was reported for people of African descent (0–1.2%) (2, 10, 76). A study conducted on Mexican Amerindians reported that geographical and ethnic factors influence the prevalence of the MTHFR 677C/T polymorphism. The prevalence in the Southeastern region ranged from 23 to 71%, followed by the Southern region (0.08 to 68%) and the Central Eastern region (0 to 57%), while the lowest frequencies were observed in Mexican Amerindians from the Northern region (0 to 19%) (77). In Europe, the lowest frequency of TT genotype was recorded in Germany (7.8%), followed by the Netherlands (8.9%), Norway (9.4%), France (9.8%), Sweden (10.3%), and Ireland (10.8%), while the highest prevalence was observed in Britain (13.2%) and Italy (18%) (76). The genotype frequencies obtained in this study (45.6% of CC, 47.5% of CT, and 6.9% of TT) are in agreement with the global occurrence of the TT genotype (10), and are comparable to those found in earlier studies conducted in the Croatian population (78–81).
4.2 The effects of the MTHFR gene polymorphism on folate and B12 levels
MTHFR 677C/T polymorphism is responsible for a 50–60% decrease of MTHFR activity in individuals with the homozygous recessive (TT) genotype compared to the wild-type (CC) genotype, leading to a decreased production of 5-MTHF from folate. This, in turn, impairs the conversion of homocysteine to methionine, resulting in elevated homocysteine levels and hyperhomocysteinemia (2). In the absence of sufficient levels of vitamin B12 and folate, which are crucial cofactors in this reaction, the body will accumulate homocysteine. A study of the German population with varying dietary habits found that vegetarians with homozygous recessive genotype (TT) had slightly lower (p = 0.06) serum folate levels than vegetarians with wild-type genotype (CC) (82). The association between the TT genotype and lower folate levels in comparison to those with the CC genotype has also been observed in study of the Chinese hypertensive population (83). In our study, CT + TT genotype exhibited significantly lower (ANOVA F(1,156) = 7.8, p = 0.006) plasma folate levels than those with the CC genotype. These results were also supported by the results of multiple regression analysis (β for genotype: −0.19 [−0.32, −0.06], p = 0.005; Supplementary Table S6). Furthermore, with a deeper post-hoc analysis we observed that non-vegetarians with the CT + TT genotype exhibited significantly lower (p = 0.039) plasma folate levels (11.9 (9.06–16.2) nmol L−1) than those with the CC genotype (14.0 (12.1–18.4) nmol L−1). Similarly, vegetarians with CT + TT genotype had lower folate levels (13.5 (10.6–17.5) nmol L−1) in comparison to those with the CC genotype (15.8 (12.7–20.1) nmol L−1), although the difference was not statistically significant (p = 0.094), which might be explained by lower statistical power in post-hoc tests. Interestingly, when obtained values were compared to the cut-off value for folate deficiency defined by the manufacturer [<7.63 nmol L−1 (84)], we discovered that a lower percentage of individuals with the CC genotype were folate deficient (6.9% of non-vegetarians and 3.3% of vegetarians) compared to those with CT and TT genotypes (14.3% of non-vegetarians and 7.9% of vegetarians, respectively). These results suggest that individuals with the CT + TT genotype may be more susceptible to hyperhomocysteinemia caused by folate deficiency than those with the wild-type (CC) genotype. However, the design of our study does not allow us to draw firm conclusions regarding the impact of the MTHFR 677C/T polymorphism on folate metabolism, as homocysteine levels, a key downstream biomarker in the 5-MTHF metabolic pathway, were not measured. Therefore, we cannot determine whether the observed differences are caused directly by disruptions in folate and homocysteine metabolism or result from other mechanisms linked to the MTHFR variant, such as epigenetic, neurochemical, oxidative, vascular, or gene–nutrient pathways. According to the literature, the only known mechanism responsible for intracellular activation of all forms of folate to 5-MTHF is MTHFR (2, 85). 5-MTHF serves as a methyl donor for the conversion of homocysteine to methionine and as a precursor to S-adenosylmethionine (SAM), which is further a cofactor or donor of methyl groups for DNA methyltransferases (86). A polymorphism in the MTHFR gene causes instability and reduces enzymatic activity, which leads to decreased methyl folate synthesis. This, in turn, diminishes one-carbon metabolism required for DNA synthesis and methylation (85, 87). Additionally, a meta-analysis of folic acid intervention studies revealed that individuals with the TT genotype had significantly lower folate status than CC/CT and CC genotypes (41). Surprisingly, these differences in folate status between genotypes persisted even after a minimum of 8 weeks of folate supplementation, despite an increase in mean serum folate levels, suggesting that individuals with TT genotype did not respond as expected to folic acid supplementation (41). According to a German population study, TT individuals with low folate or vitamin B12 status had a higher risk of hyperhomocysteinemia than those with the CC genotype, who have high levels of both folate and B12 (82). Lower plasma folate and vitamin B12 levels in TT genotype with the mild folate-deficiency was reported by Pereira et al. (88). Opposite to this finding, we did not detect any differences between the genotypes for vitamin B12 and this disparity between our and Pereira’ study can likely be attributed to variations in individuals’ supplementation to prevent B12 deficiency. For instance, 42% of our participants took vitamin B12 supplements. However, we are unable to make this assertion with certainty, since Pereira et al. (88) did not provide any information regarding vitamin B12 supplementation. Vitamin B12 deficiency results in lower remethylation of homocysteine and production of S-adenosylmethionine, which in turn affects the methyltransferase activity. Additionally, 5-MTHF becomes trapped between MTHFR and inhibited methionine synthase, rendering it incapable of participating in the remethylation pathway (89, 90).
4.3 The effects of the MTHFR gene polymorphism on lipid status
Numerous studies have explored the association between MTHFR 677C/T polymorphism and lipid levels, but the findings have been inconsistent (36, 38, 39, 42). Some authors (91, 92) suggested that these inconsistencies could be an indirect consequence of homocysteine mediation or modifications in methylation status, given that MTHFR is an important enzyme in the one-carbon metabolism, wherein 5-MTHF serves as an important methyl donor for various molecules, including homocysteine, DNA, and proteins (93). However, other factors like ethnicity, age, sex, and health status, may also modulate the relationship between lipids and polymorphism, contributing to the variability in findings (39, 42, 94, 95). For instance, while some studies have been unable to establish a correlation between the MTHFR 677C/T polymorphism and lipid status in the overall study population, differences became apparent when data were stratified by sex. In several studies, only female with the CT + TT genotype had higher levels of triglycerides and LDL cholesterol in comparison to those with the CC genotype (39, 94). In contrast, Zhang et al. (95) found that, in their study of two distinct Chinese ethnic groups, only males with the CT + TT genotype had higher levels of both LDL cholesterol and triglycerides when data were stratified by sex, although HDL cholesterol levels did not differ by genotypes. In our study, which focused on population stratified by diet, we found a significantly higher (p = 0.049) LDL cholesterol level of non-vegetarians with CT + TT genotype (1.97 (1.54–2.32) nmol L−1) in comparison to those with the CC genotype (1.63 (1.41–2.06) nmol L−1). These ANOVA results indicated higher LDL cholesterol in non-vegetarians with the MTHFR 677C/T polymorphism, suggesting that MTHFR 677C/T could be a risk factor for increasing LDL cholesterol in plasma. However, these results did not remain significant in the regression analysis after adjusting for potential confounders, implying that the apparent effect is context-dependent and likely influenced by confounding. Nevertheless, several biologically plausible pathways link dyslipidemia to MTHFR-related hyperhomocysteinemia. Elevated homocysteine can activate transcription factors [nuclear factor Y (NF-Y), cyclic adenosine monophosphate response element-binding protein (CREB), sterol regulatory element-binding protein 2 (SREBP-2)], increasing mRNA expression and activity of HMG-CoA reductase, promoting hepatic lipid accumulation, and ultimately contributing to hypercholesterolemia (96). In another pathway, elevated homocysteine can induce endoplasmic reticulum stress in hepatocytes, endothelial cells, and arterial smooth muscle cells, triggering the unfolded protein response and SREBP pathways and up-regulating genes involved in cholesterol synthesis and uptake, as well as the intracellular cholesterol acumulation (97). The loss of statistical significance after adjustment emphasizes the complexity of the relationships between MTHFR 677C/T and lipids, which are probably influenced by other metabolic modifiers and nutritional status. However, given the mechanisms described above, our initial ANOVA findings have biological plausibility. A limitation of the current study is that plasma homocysteine levels were not measured, which prevents a proper and logical explanation of the relationship between the examined polymorphism and serum lipid levels. Additionally, significant interaction between diet and genotype was observed for HDL cholesterol after adjustment for age, sex, smoking, BMI and education, indicating opposing trends among vegetarians (lower HDL values for CC genotype) and non-vegetarians (higher HDL for CC genotype). Diet and eating habits are well-recognized important factors that influence gene-lipid interactions. While our study design does not allow us to draw conclusions about mechanisms or pathways, the presence of this interaction suggests that the MTHFR 677C/T polymorphism can influence HDL cholesterol levels depending on dietary habits, especially among vegetarians and non-vegetarians. A plausible explanation is that MTHFR activity influences homocysteine metabolism and methylation processes, which in turn regulates lipid metabolism (98, 99). Carriers of the MTHFR 677C/T polymorphism exhibit significantly reduced MTHFR enzyme activity and elevated homocysteine levels (10, 100). Associated increased oxidative stress has been linked to elevated homocysteine and decreased expression of transport proteins like ATP-binding cassette transporter A1 (ABCA1), impairing HDL biogenesis (101, 102). Reduced levels of SAM in individuals with the MTHFR 677C/T polymorphism may also reduce methylation of genes involved in lipoprotein metabolism, such as apolipoprotein A-1 (APOA1) and lecithin-cholesterol acyltransferase (LCAT) (10, 103, 104). Moreover, it has been reported that many nutraceuticals, such as dietary supplements, can influence lipid profiles and dyslipidemia, potentially exerting a greater effect than individual genetic differences (105, 106).
Although we did not perform analyses stratified by sex, our regression models were adjusted for this variable and showed that male sex was not a significant predictor of plasma lipid status. A recent systematic review and meta-analysis by Luo et al. (42) reported that individuals with the CT and TT genotypes had higher levels of total and LDL cholesterol compared to those with the CC genotype. However, no association was found with triglycerides and HDL cholesterol. Interestingly, when data were stratified by sex and ethnicity, it became evident that the obtained associations were specific to female individuals and were predominantly observed in the Asian population. Gender differences may partly result from women consuming more folic acid than men. This, in turn, affects the increase in homocysteine levels and the differences among genotypes, leading to higher homocysteine levels and increased lipid production in those with the MTHFR 677C/T polymorphism (107). Sex-specific associations between lipid plasma levels and MTHFR 677C/T polymorphism have also been observed in studies of a long-lived Chinese cohort (39) and an overweight Chinese Han population (38). However, caution is needed when interpreting results of latter study, as Zhi et al. (38) found no significant correlation between the MTHFR 677C/T polymorphism and lipid levels in healthy Han individuals with normal BMI. Conversely, obesity and dyslipidaemia are known to be associated with elevated levels of triglycerides, total cholesterol, and LDL cholesterol, as well as a reduction in HDL cholesterol (108). This suggests that the effect attributed to the polymorphism on lipid levels may have been caused by obesity or dyslipidaemia, rather than the polymorphism itself, as these effects were not present in the healthy control group, indicating that the impact of MTHFR polymorphisms on lipid levels also depends on the health status of the population studied. In another study conducted by Yin et al. (109), both males and females with the CT + TT genotype and normal weight had higher levels of total, HDL, and LDL cholesterol in comparison to those with the CC genotype. On the other hand, overweight individuals with the CT + TT genotype had higher HDL cholesterol and triglyceride levels, but lower total and LDL cholesterol levels, than the CC genotype. Interestingly, a study of hypertensive individuals (110) found no significant association between the MTHFR 677C/T polymorphism and any blood lipid levels. The interaction between the health status of the population and lipid levels was also confirmed in studies that exclusively examined a single sex with specific health conditions, leading to variable results. Therefore, in a study of obese female adolescents, no significant differences in lipid parameters were observed between genotypes (111), whereas women with polycystic ovarian syndrome (PCOS) carrying the CT genotype exhibited higher triglyceride and cholesterol levels compared to those with the CC genotype (112).
4.4 The effects of the MTHFR gene polymorphism on element levels
Unlike biochemical markers, there is a scarcity of research that investigates the relationship between the MTHFR 677C/T polymorphism and element levels in the body. To the best of our knowledge, this study is among the few to examine the effects of this polymorphism on both toxic and essential elements (see Joneidi et al. (25) for a review).
The only literature currently available that investigates the link between the MTHFR 677C/T polymorphism and essential elements is limited to Fe, as no evidence was found linking other essential elements to this gene polymorphism. Bortolus et al. (46) investigated the effects of four one-carbon metabolism-related polymorphisms, including MTHFR 677C/T, on the Fe and ferritin status of healthy Italian men and women. They found no correlation between the investigated parameters, which is in accordance with our findings regarding the total Fe content in plasma. In addition, we observed significantly higher plasma Ca levels in individuals with the CT + TT genotype (β = 0.01 [0.003, 0.03], p = 0.015; Supplementary Table S5). This association raises an interesting question about the possible link between this genetic variant and the regulation of Ca homeostasis, although direct evidence connecting MTHFR 677C/T polymorphism to circulating Ca is, to our knowledge, lacking. Prior meta-analysis has reported associations between the MTHFR 677C/T and bone mineral density (113) as well as fracture risk (114). Because abnormalities in Ca metabolism are associated with the risk of osteopenia and osteoporosis (115) and excessive vascular Ca accumulation is associated with elevated cardiovascular risk (116), the observed genotype-related difference in plasma Ca may have clinical relevance. As mentioned above, MTHFR 677C/T polymorphism affects folate and homocysteine metabolism. Mendelian-randomization analyses by Casas et al. (117), who observed that the MTHFR T allele carriers have an increased risk of stroke, consistent with allele-related elevations in homocysteine concentration, supports a causal link between elevated homocysteine and stroke risk. It is therefore possible that the T allele, through mechanisms that still need to be clarified, influences various aspects of Ca metabolism with potential impacts on bone and vascular health. Given these considerations, our finding of higher plasma Ca in CT + TT genotype emphasize the need for long-term studies that can evaluate the clinical effects of these changes on the risk of osteoporosis and cardiovascular disease, as well as for functional studies examining the expression of Ca transport proteins and ion channels in relation to the MTHFR genotype.
In the present study, we observed that participants with the CT + TT genotype had lower levels of As in their blood than individuals with the CC genotype (β = −0.47 [−0.84, −0.09], p = 0.017; Supplementary Table S5), suggesting reduced susceptibility to As exposure through food. Nevertheless, we were unable to conclusively confirm this finding, as we employed blood as a biological indicator of exposure to arsenic from food, primarily from fish and seafood, which predominantly contain organic arsenic compounds that are less toxic than inorganic As (iAs) (118, 119). In contrast, most previous studies examining As metabolism in relation to MTHFR 677C/T polymorphism focused on exposure to iAs by analyzing its metabolites in urine samples. It is widely recognized that the MTHFR 677C/T polymorphism may lead to a decrease in methylation capacity. In the event of human exposure to inorganic As, this may lead to the incomplete methylation of monomethylarsonic acid (MMA) to dimethylarsinic acid (DMA), thereby increasing the risk of As-induced cancer (120, 121). According to a study conducted on populations from Hungary, Romania, and Slovakia (122), which investigated polymorphism and other factors influencing the metabolism of iAs at low environmental levels, males with the CT + TT genotype had a lower proportion of dimethylarsinic acid (%DMA) and a higher percentage of monomethylarsonic acid (%MMA) in their urine than those who were homozygous for the wild-type allele (CC). This suggests that the CT + TT genotype is more susceptible to the accumulation of iAs in the body and an increased occurrence of toxic effects caused by iAs in comparison to the CC genotype. Similar results were observed in Indigenous women from Argentina who were also exposed to iAs through water (123).
The MTHFR 677C/T polymorphism has been associated with several other potentially toxic elements. Saad-Hussein et al. (49) found that occupationally exposed workers with the TT genotype had higher urinary Pb concentrations than those with other genotypes. They hypothesized that this might be due to decreased bone mineral density, a condition that has been linked to MTHFR gene polymorphisms (124). However, Golbahar et al. (125) failed to establish a correlation between the MTHFR 677C/T polymorphism and bone mineral density in postmenopausal women, despite directly measuring bone density. In this study, both vegetarians and non-vegetarians with CT + TT genotype had slightly higher blood Pb levels in comparison to those with the CC genotype, although the difference was not statistically significant. Liu et al. (47) reported that CT and TT with higher Cd levels had an increased risk of neural tube defects. Although the authors did not discuss the relationship between Cd levels and MTHFR polymorphism, the presented proportions of new-borns with high Cd levels did not show a statistically significant difference between the CC and CT/TT genotypes. In a study of Inuit population, Parajuli et al. (48) reported that individuals with the CT + TT genotype had lower Hg levels compared to those with the CC genotype. By contrast, our findings indicated slightly higher Hg levels in participants with the CT + TT genotype compared to the CC genotype, although this difference was not statistically significant. The observed discrepancy between studies could stem from the traditionally higher exposure of the Inuit population to Hg through fish and seafood consumption, leading to 2–5 times higher blood Hg levels than those recorded in the Croatian population [this study; (52, 126, 127)] and 16 times higher than in the average Canadian adult (128). However, these conclusions are limited by the lack of supplementary information from Parajuli et al. (48) regarding potential exposure sources, including fish and seafood intake, which may influence Hg levels (129–131). Nevertheless, it is important to recognize that, in addition to the potential harmful impact of Hg from seafood on health, a moderate fish consumption may also have a beneficial role in maintaining DNA integrity (132).
Our results of analysis of variance revealed that plasma Al levels in vegetarians with CT + TT genotype were significantly lower in comparison to the CC genotype. This suggests that vegetarians with the CT + TT genotype might accumulate lower Al levels in the body in comparison to the homozygous CC genotype. This is a novel observation not previously documented in the literature, but it should be interpreted with caution since it was not observed in non-vegetarians. This finding may be of interest in light of recent reports linking the presence of Al in the body (urine and/or hair) with symptoms of autism spectrum disorder (ASD) in children (133–135). In parallel, a recent systematic review and meta-analysis (136, 137) indicated a significant association between the MTHFR 677C/T polymorphism and elevated ASD risk; nevertheless, these results should be interpreted cautiously due to small sample sizes and heterogeneity in study characteristics (e.g., ethnicity, genotyping methods, etc.). Furthermore, our study did not assess ASD outcomes and we cannot infer any relationship between plasma Al levels, dietary habits, MTHFR 677C/T polymorphism, and ASD. We also did not examine the expression of metallothioneins or other metal-binding proteins that could influence Al bioavailability. Although our results suggest that a vegetarian diet potentially reduces the risk of neurotoxic effects of Al through higher intake of antioxidants that reduce oxidative stress and bind metals, inadequate intake of vitamin B12 (crucial for methylation and neutralization of homocysteine) as well as suboptimal intake of other key nutrients, such as Zn, Fe, and omega-3 fatty acids, could also increase the risk in individuals with the MTHFR polymorphism. Further studies are needed to elucidate the molecular mechanisms behind this association. It is plausible that different mechanisms are involved – for example, MTHFR 677C/T polymorphism may increase the risk for autism through methylation disorders and folate metabolism, independent of Al, while Al act as an additional or separate risk factor that accumulates differently across different genotypic backgrounds. To better understand these relationships, further research should focus on the relevant molecular pathways. Since there is growing evidence that both genetic and environmental factors contribute to ASD, these fascinating findings encourage combined research in this area, which could improve our understanding of the mechanisms that may link exposure to Al and the MTHFR 677C/T polymorphism to ASD (138). In this regard, our findings could be helpful in the development of epidemiologic studies that examine the interaction between MTHFR polymorphism and the level of Al in the body and their influence on human health. More broadly, our findings raise the question of potential interactions between the MTHFR 677C/T polymorphism and exposure to environmental toxins, such as toxic metals that cause oxidative stress and interfere with enzyme activity. In T allele carriers, exposure to heavy metals may worsen metabolic imbalances, as the MTHFR 677C/T polymorphism already reduces the availability of 5-methyltetrahydrofolate and elevates levels of homocysteine. Individual differences in sensitivity to metal poisoning may partly reflect gene–environment interactions. To clarify mechanism and clinical relevance, future studies should integrate genetic predisposition with detailed heavy metal exposure assessment. If carriers of specific genetic variants such as MTHFR 677C/T are confirmed to have increased vulnerability to environmental toxins, including metals, this could guide exposure-prevention strategies and support the development of personalized public health approaches.
4.5 The effects of MTHFR gene polymorphism on amino acids
For many years, it has been demonstrated that small variations in the activity of certain enzymes can lead to minor modifications in metabolic fluxes, but can have substantial effects on metabolite concentrations (139). Merging metabolomics and genetic data is necessary to understand the interplay of genes, environment, and diet (140). This investigation involved a metabolomic analysis of plasma samples to identify the amino acids that could distinguish between individuals with and without the MTHFR 677C/T polymorphisms. We observed that participants with CT + TT genotype had elevated levels of proline, glutamic acid, and aspartic acid in comparison to the CC genotype, while glycine was elevated only in vegetarians with the CT + TT genotype. The literature on metabolomic changes associated with the MTHFR polymorphism is scarce, particularly with respect to variations in amino acids. In the Swedish National Study on Aging and Care, the MTHFR 677C/T polymorphism was associated with the homocysteine-to-methionine ratio, with lower levels in the TT genotype compared to the CC genotype. However, no differences were observed in the levels of methionine. Additionally, the research showed that the CT and TT genotypes with low methionine levels were more susceptible to developing cardiovascular multimorbidity than individuals with the CC genotype (45). In this study, methionine levels did not differ across the genotypes. However, the levels of glycine in vegetarians with the CT + TT genotype were significantly higher than those in the CC genotype, suggesting that vegetarians with CT + TT genotype may be more prone to mild hyperhomocysteinemia than those with the CC genotype. The results of our research on the adult population can be partially compared to those of a Spanish study conducted on cerebral creatine transporter deficiency patients, who underwent a regimen of oral glycine and L-arginine therapy administered twice daily over a period of 9 months. The authors observed that patients with the CT + TT genotype exhibited varying degrees of hyperhomocysteinemia, with the severity of the condition increasing with increasing dosages of glycine and L-arginine (43). It is well-established that the protein glycine N-methyltransferase (GNMT) is linked with MTHFR by S-adenosylmethionine (SAM). When SAM levels are elevated, MTHFR is inhibited, which in turn reduces 5-MTHF, thereby enabling GNMT to reduce SAM levels. Conversely, reduced SAM levels lead to higher production of 5-MTHF, which inhibits GNMT and spares methyl groups for DNA methylation (141–145). A recent study conducted on pregnant women reported that women with MTHFR polymorphism had lower levels of 1-monohexadecanoylglycerol, pyrophosphate, benzoin, and linoleic acid, and higher levels of glyceric acid, L-tryptophan, L-alanine, L-proline, norvaline, L-threonine, and myo-inositol than women without polymorphism (44). Our finding of higher proline levels in the CT + TT genotype agrees with their findings. Nevertheless, it is important to take caution when interpreting these data, as the authors did not specify whether they investigated MTHFR 677C/T or 1298A/C polymorphism, and their data was based on a small sample size of only 18 women (10 without and 8 with MTHFR polymorphism). Proline is essential for the tendons and joints, but it can also act as a neurotoxin and metabotoxin in certain circumstances (65).
As previously mentioned, the literature contains a limited number of studies that examine the relationship between human amino acid profiles and MTHFR polymorphisms. To the best of our knowledge, this study is the first to report the significant relationship between this polymorphism and glutamic and aspartic acid in humans. The sole study in the literature that addresses a similar subject is a pilot study conducted on 67 individuals diagnosed with schizophrenia (146), which utilized the homeostatic model assessment of insulin resistance (HOMA-IR) as a marker of cardiovascular disease. Unfortunately, while the authors identified the link between HOMA-IR and two combined MTHFR polymorphisms (677C/T and 1298A/C), as well as the associations between HOMA-IR and specific amino acids (a negative correlation with glycine, serine, and betaine, and a positive correlation with glutamate and threonine), they did not explore the relationship between amino acid profiles and MTHFR polymorphisms. We also observed significant gene-diet interactions for asparagine, alanine, and glycine, with the latter two remaining significant after multivariate regression analysis adjustment. Although the observational design of our study does not permit causal inference about possible biological mechanisms underlying this interaction, these patterns may be partly explained by their involvement in one-carbon metabolism and amino acid interconversion pathways. For instance, serine hydroxymethyltransferase (SHMT) uses 5,10-methylenetetrahydrofolate for the conversion of glycine into serine, which is generated by the cytosolic nicotinamide adenine dinucleotide phosphate (NADPH)-dependent activity of MTHFR (11, 147). These pathways may be significantly impacted by dietary habits (vegetarian vs. non-vegetarian), which differ in the quantity and quality of micronutrients and protein. These considerations provide testable hypotheses for the observed interactions; their confirmation will require genotype-stratified studies integrating amino acid intakes and one-carbon metabolism biomarkers to clarify the link between MTHFR polymorphism and amino-acid metabolism.
4.6 Study strengths and limitations
This study is one of the few that offers new insights into the impact of MTHFR genetic variants on an individual’s sensitivity to metals and the amino acid profile in human plasma. This is of relevance when assessing how one-carbon-related genetic variations might affect both element and amino acid status. An important strength of this study is the similarity in lifestyle and anthropometric traits between the two dietary groups, which helps to minimize any bias caused by differences in population characteristics. Additionally, the data linking dietary habits with element and amino acid levels, alongside biochemical parameters, may be valuable for designing future studies targeting specific populations, particularly those with individual sensitivities.
However, several limitations should be addressed when interpreting these findings. First, due to the low prevalence rate of the TT genotype (~7%) in our study group (n = 162), we were unable to assess its independent effects on the measured variables reliably. Analyzing the TT genotype separately would have substantially reduced the power of statistical tests and increased the risk of generating potentially misleading results. While adequate statistical power for detecting interactions of medium to large effect sizes was maintained in the dominant-model approach, grouping TT and CT genotypes might have obscured the TT-specific effects and introduced bias in interpreting the role of the TT variant. Furthermore, the study was underpowered for detecting interactions of small effect sizes. Second, although the study aimed to include a sample broadly representative of Croatian adults (vegetarians and non-vegetarians), full representativeness across all sociodemographic strata could not be ensured. The use of convenience sampling approach introduces the risk of selection bias. Furthermore, the high percentage of participants with doctoral degrees (particularly in the non-vegetarian group) may indicate a higher level of health awareness and health-oriented behavior compared to the general population, which could have affected their dietary preferences, supplement use and metabolic profiles. Third, our study group was unbalanced in terms of sex (72% females vs. 28% males), which may limit the generalizability of the results. However, this distribution aligns with global trends showing higher prevalence of vegetarian dietary patterns among women (148), and a greater willingness of women to participate in research (149, 150). Although regression models were adjusted for sex, the smaller sample size for male participants reduced power to detect sex-specific effects. Furthermore, the interaction effects for sex were not assessed. Future studies should include larger, more diverse, and sex-balanced populations to improve the external validity of the findings. Additionally, since the study was conducted within a single cohort, the results might not be generalizable to other populations with different dietary habits, lifestyle factors/behaviors. Despite these limitations, our research offers new integrative data linking dietary practices, the MTHFR gene, and several biochemical markers within a real-life cohort. Fourth, no urine samples were obtained, which prevented us to quantify the levels of iAs metabolites (monomethylarsonic acid, MMA, and dimethylarsinic acid, DMA). Consequently, we could not assess methylation capacity of the CT + TT genotype in comparison to those with the CC genotype. Fifth, we were unable to determine plasma homocysteine levels due to limitations in the study’s design and inadequate sample preparation for homocysteine analysis, thereby preventing us to assess well established relationship between homocysteine and the MTHFR 677C/T polymorphism (41, 151). Although we assessed homocysteine levels in a previous study, in which we found that vegetarians had higher levels than non-vegetarian subjects (12.30 ± 3.80 μmol L−1 vs. 10.70 ± 2.60 μmol L−1) (18), these data could not be integrated with the current dataset or directly compared due to differences in study design and study population. Consequently, we cannot conclusively determine whether the observed differences between the analyzed indicators are mediated by homocysteine or by mechanisms independent of homocysteine. Thus, our understanding of the mechanistic connection between the MTHFR genotype and observed effects remains limited. To examine the direct and indirect effects of MTHFR on biochemical parameters, elements, and amino acid profiles, future studies should include homocysteine in the analytical panel.
Sixth, aside from the use of dietary supplements, complete information on dietary habits, such as total calorie intake, distribution of macronutrients, and specific nutrient intake, was not included in the analysis. Instead, a broad classification (vegetarian vs. non-vegetarian) was used to reflect dietary habits. This limited our ability to capture the full range of individual dietary patterns and nutrient intake and may have constrained interpretation of interactions between MTHFR polymorphisms, diet, and biochemical outcomes. A detailed dietary assessment within this cohort is ongoing and will be addressed in future analysis. In addition to diet, various external factors, including water quality, environmental contamination, fish and seafood intake, and occupational exposure, can also influence the levels of both toxic and essential elements in blood and/or plasma. However, as the current study did not analyse these data on these non-dietary exposures, our ability to fully distinguish non-dietary environmental influences from nutritional determinants of elemental status is limited. Seventh, the classification accuracy for sparse discriminant analysis was evaluated on the same dataset that was used for fitting the model. Thus, the classification accuracy on an independent dataset is expected to be lower. Independent evaluation on a new sample is recommended in future studies. Lastly, potential sources of bias and confounding factors must be considered. The high proportion of individuals with PhD, the sex imbalance, and the convenience sampling could have influenced the results. Unmeasured factors such as physical activity, environmental exposures, and micronutrient intake may also distort the observed associations. Additionally, since the results are based on a Croatian population, the findings might not be widely applicable to other groups due to differences in dietary habits, lifestyle choices, and environmental factors.
5 Conclusion
In conclusion, this study demonstrates the significant impact of dietary habits on most of the analysed parameters and offers new insights into the relationship between the MTHFR 677C/T polymorphism and folate levels, showing that individuals with this polymorphism had lower folate levels. The study also found associations between the MTHFR 677C/T polymorphism and plasma Al and Ca levels, and identified novel relationships with the amino acids proline, glycine, aspartic acid, and glutamic acid. Overall, our results primarily highlight the significance of gene–diet interactions associated with trace elements, biochemical parameters, and amino acid profiles, underscoring the necessity for additional study. Although these findings may contribute to developing more personalized dietary plans over time, replication in larger, genotype-stratified cohorts and clinical validation will be required before translation into practice. To better distinguish between the respective contributions of gene-diet and gene–environment interactions, future studies should include more detailed dietary intake data and relevant environmental exposure data.
Data availability statement
The data that support the findings of this study are available from the corresponding authors, upon reasonable request.
Ethics statement
The studies involving humans were approved by the Ethics Committee of the Institute for Medical Research and Occupational Health, Zagreb, Croatia. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
TO: Data curation, Formal analysis, Investigation, Methodology, Supervision, Validation, Visualization, Writing – original draft. JK: Formal analysis, Investigation, Methodology, Software, Visualization, Writing – original draft. IP: Data curation, Investigation, Writing – review & editing. ZK-G: Formal analysis, Investigation, Methodology, Writing – review & editing. DP: Formal analysis, Investigation, Methodology, Validation, Writing – review & editing. AS: Conceptualization, Formal analysis, Investigation, Methodology, Supervision, Visualization, Writing – original draft, Writing – review & editing. AD: Formal analysis, Investigation, Methodology, Validation, Writing – review & editing. BTL: Formal analysis, Investigation, Methodology, Validation, Writing – review & editing. AJ: Formal analysis, Investigation, Methodology, Validation, Writing – review & editing. AP: Investigation, Writing – review & editing. FC: Formal analysis, Investigation, Methodology, Writing – review & editing. LD: Formal analysis, Investigation, Methodology, Validation, Writing – review & editing. MG: Data curation, Investigation, Writing – review & editing. IJ: Data curation, Funding acquisition, Investigation, Project administration, Resources, Supervision, Writing – review & editing. GG: Conceptualization, Data curation, Funding acquisition, Investigation, Project administration, Resources, Supervision, Writing – review & editing.
Funding
The author(s) declared that financial support was received for this work and/or its publication. This study is an integral part of the research projects “Interactions between Human and Environmental Health: Determinants of Health Preservation” (HumEnHealth), and “Assessment of Efficacy and Toxicity of Biologically Active Substances” (BioMolTox), funded by the European Union -Next Generation EU, Program Contract of 8 December 2023, Class: 643-02/23-01/00016, Reg. no. 533-03-23-0006.
Acknowledgments
This study was performed using the facilities and equipment funded within the European Regional Development Fund project KK.01.1.1.02.0007 “Research and Education Centre of Environmental Health and Radiation Protection – Reconstruction and Expansion of the Institute for Medical Research and Occupational Health.” All authors are non-vegetarians.
Conflict of interest
The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that Generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnut.2025.1710613/full#supplementary-material
References
1. Bhaskar, LVKS, and Saikrishna, L. Meta-analysis of MTHFR polymorphisms and pancreatic cancer susceptibility In: GP Nagaraju, Ahmad SBT-TA for PC, editors. Theranostic approach for pancreatic cancer. Cambridge, Massachusetts, USA: Academic Press (2019). 263–74.
2. Liew, SC, and Das, GE. Methylenetetrahydrofolate reductase (MTHFR) C677T polymorphism: epidemiology, metabolism and the associated diseases. Eur J Med Genet. (2015) 58:1–10. doi: 10.1016/j.ejmg.2014.10.004,
3. Leclerc, D, Sibani, S, and Rozen, R. Molecular biology of methylenetetrahydrofolate reductase (MTHFR) and overview of mutations/polymorphisms. Austin, TX, USA: Landes Bioscience (2015).
4. Ganguly, P, and Alam, SF. Role of homocysteine in the development of cardiovascular disease. Nutr J. (2015) 14:1–10. doi: 10.1186/1475-2891-14-6,
5. Huang, T, Ren, JJ, Huang, J, and Li, D. Association of homocysteine with type 2 diabetes: a meta-analysis implementing mendelian randomization approach. BMC Genomics. (2013) 14:867. doi: 10.1186/1471-2164-14-867,
6. Refsum, H, Nurk, E, Smith, AD, Ueland, PM, Gjesdal, CG, Bjelland, I, et al. The Hordaland homocysteine study: a community-based study of homocysteine, its determinants, and associations with disease. J Nutr. (2006) 136:1731S–40S. doi: 10.1093/jn/136.6.1731s,
7. Obersby, D, Chappell, DC, Dunnett, A, and Tsiami, AA. Plasma total homocysteine status of vegetarians compared with omnivores: a systematic review and meta-analysis. Br J Nutr. (2013) 109:785–94. doi: 10.1017/S000711451200520X,
8. Lyon, P, Strippoli, V, Fang, B, and Cimmino, L. B vitamins and one-carbon metabolism: implications in human health and disease. Nutrients. (2020) 12:1–24. doi: 10.3390/nu12092867,
9. Troesch, B, Weber, P, and Mohajeri, MH. Potential links between impaired one-carbon metabolism due to polymorphisms, inadequate B-vitamin status, and the development of Alzheimer’s disease. Nutrients. (2016) 8:803. doi: 10.3390/nu8120803,
10. Zarembska, E, Ślusarczyk, K, and Wrzosek, M. The implication of a polymorphism in the methylenetetrahydrofolate reductase gene in homocysteine metabolism and related civilisation diseases. Int J Mol Sci. (2024) 25:1–31. doi: 10.3390/ijms25010193,
11. Ducker, G, and Rabinowitz, J. One-carbon metabolism in health and disease. Cell Metab. (2017) 25:27–42. doi: 10.1016/B978-0-12-804274-8.00054-0,
12. Dai, Z, Zheng, W, and Locasale, JW. Amino acid variability, tradeoffs and optimality in human diet. Nat Commun. (2022) 13:6683. doi: 10.1038/s41467-022-34486-0,
13. Parker, HW, and Vadiveloo, MK. Diet quality of vegetarian diets compared with nonvegetarian diets: a systematic review. Nutr Rev. (2019) 77:144–60. doi: 10.1093/nutrit/nuy067,
14. Craig, WJ. Nutrition concerns and health effects of vegetarian diets. Nutr Clin Pract. (2010) 25:613–20. doi: 10.1177/0884533610385707,
15. Key, TJ, Appleby, PN, Crowe, FL, Bradbury, KE, Schmidt, JA, and Travis, RC. Cancer in British vegetarians: updated analyses of 4998 incident cancers in a cohort of 32,491 meat eaters, 8612 fish eaters, 18,298 vegetarians, and 2246 vegans. Am J Clin Nutr. (2014) 100:378S–85S. doi: 10.3945/ajcn.113.071266,
16. Lousuebsakul-Matthews, V, Thorpe, DL, Knutsen, R, Beeson, WL, Fraser, GE, and Knutsen, SF. Legumes and meat analogues consumption are associated with hip fracture risk independently of meat intake among Caucasian men and women: the Adventist health Study-2. Public Health Nutr. (2013) 17:2333–43. doi: 10.1017/S1368980013002693,
17. Gajski, G, Gerić, M, Jakaša, I, Peremin, I, Domijan, AM, Vučić Lovrenčić, M, et al. Inflammatory, oxidative and DNA damage status in vegetarians: is the future of human diet green? Crit Rev Food Sci Nutr. (2023) 63:3189–221. doi: 10.1080/10408398.2021.1986464,
18. Gajski, G, Gerić, M, Vučić Lovrenčić, M, Božičević, S, Rubelj, I, Nanić, L, et al. Analysis of health-related biomarkers between vegetarians and non-vegetarians: a multi-biomarker approach. J Funct Foods. (2018) 48:643–53. doi: 10.1016/j.jff.2018.07.054
19. Elorinne, AL, Alfthan, G, Erlund, I, Kivimäki, H, Paju, A, Salminen, I, et al. Food and nutrient intake and nutritional status of Finnish vegans and non-vegetarians. PLoS One. (2016) 11:1–14. doi: 10.1371/journal.pone.0148235,
20. Li, D. Chemistry behind vegetarianism. J Agric Food Chem. (2011) 59:777–84. doi: 10.1021/jf103846u,
21. Kristensen, NB, Madsen, ML, Hansen, TH, Allin, KH, Hoppe, C, Fagt, S, et al. Intake of macro- and micronutrients in Danish vegans. Nutr J. (2015) 14:115–21. doi: 10.1186/s12937-015-0103-3,
22. Schmidt, JA, Rinaldi, S, Scalbert, A, Ferrari, P, Achaintre, D, Gunter, MJ, et al. Plasma concentrations and intakes of amino acids in male meat-eaters, fish-eaters, vegetarians and vegans: a cross-sectional analysis in the EPIC-Oxford cohort. Eur J Clin Nutr. (2016) 70:306–12. doi: 10.1038/ejcn.2015.144,
23. Ledda, C, Cannizzaro, E, Lovreglio, P, Vitale, E, Stufano, A, Montana, A, et al. Exposure to toxic heavy metals can influence homocysteine metabolism? Antioxidants. (2020) 9:1–14. doi: 10.3390/antiox9010030,
24. Wu, Z, Hu, H, Wang, C, Wu, J, Xiong, Y, Fu, Y, et al. Association between serum folate levels and blood concentrations of cadmium and lead in US adults. Environ Sci Pollut Res. (2022) 29:3565–74. doi: 10.1007/s11356-021-15767-5,
25. Joneidi, Z, Mortazavi, Y, Memari, F, Roointan, A, Chahardouli, B, and Rostami, S. The impact of genetic variation on metabolism of heavy metals: genetic predisposition? Biomed Pharmacother. (2019) 113:108642. doi: 10.1016/j.biopha.2019.108642,
26. Pollack, AZ, Mumford, SL, Sjaarda, L, Perkins, NJ, Malik, F, Wactawski-Wende, J, et al. Blood lead, cadmium and mercury in relation to homocysteine and C-reactive protein in women of reproductive age: a panel study. Environ Heal A Glob Access Sci Source. (2017) 16:1–9. doi: 10.1186/s12940-017-0293-6,
27. Bakulski, KM, Park, SK, Weisskop, MG, Tucker, KL, Sparrow, D, Spiro, A, et al. Lead exposure, B vitamins, and plasma homocysteine in men 55 years of age and older: the VA normative aging study. Environ Health Perspect. (2014) 122:1066–74. doi: 10.1289/ehp.1306931,
28. He, H, Zhang, Z, and Li, M. Association between serum folate concentrations and blood lead levels in adolescents: a cross-sectional study. Front Pediatr. (2022) 10:1–9. doi: 10.3389/fped.2022.941651,
29. Gallagher, CM, and Meliker, JR. Total blood mercury, plasma homocysteine, methylmalonic acid and folate in US children aged 3-5years, NHANES 1999-2004. Sci Total Environ. (2011) 409:1399–405. doi: 10.1016/j.scitotenv.2011.01.021
30. Lönnerdal, B. Dietary factors influencing zinc absorption. J Nutr. (2000) 130:1378S–83S. doi: 10.1093/jn/130.5.1378S,
31. Xiaoming, XC, Zafra-Stone, S, Bagchi, M, and Bagchi, D. Bioavailability, antioxidant and immune-enhancing properties of zinc methionine. Biofactors. (2006) 27:231–44. doi: 10.1002/biof.5520270120,
32. Pérez-Carreón, K, Martínez, LM, Videa, M, Cruz-Angeles, J, Gómez, J, and Ramírez, E. Effect of basic amino acids on folic acid solubility. Pharmaceutics. (2023) 15:2544. doi: 10.3390/pharmaceutics15112544,
33. Cooperman, JM, and Lopez, R. The role of histidine in the anemia of folate deficiency. Exp Biol Med. (2002) 227:998–1000. doi: 10.1177/153537020222701107,
34. Smith, AD. Folic acid fortification: the good, the bad, and the puzzle of vitamin B-12. Am J Clin Nutr. (2007) 85:3–5. doi: 10.1093/ajcn/85.1.3,
35. Froese, DS, Fowler, B, and Baumgartner, MR. Vitamin B12, folate, and the methionine remethylation cycle—biochemistry, pathways, and regulation. J Inherit Metab Dis. (2019) 42:673–85. doi: 10.1002/jimd.12009,
36. Kucukhuseyin, O, Kurnaz, O, Akadam-Teker, AB, Isbir, T, Bugra, Z, Ozturk, O, et al. The association of MTHFR C677T gene variants and lipid profiles or body mass index in patients with diabetic and nondiabetic coronary heart disease. J Clin Lab Anal. (2013) 27:427–34. doi: 10.1002/jcla.21623,
37. Kurzawski, M, Wajda, A, Malinowski, D, Kazienko, A, Kurzawa, R, and Drozdzik, M. Association study of folate-related enzymes (MTHFR, MTR, MTRR) genetic variants with non-obstructive male infertility in a polish population. Genet Mol Biol. (2015) 38:42–7. doi: 10.1590/S1415-475738120140179,
38. Zhi, X, Yang, B, Fan, S, Wang, Y, Wei, J, Zheng, Q, et al. Gender-specific interactions of MTHFR C677T and MTRR A66G polymorphisms with overweight/obesity on serum lipid levels in a Chinese Han population. Lipids Health Dis. (2016) 15:185–9. doi: 10.1186/s12944-016-0354-9,
39. Chen, N. Y., Liu, C. W., Du, L. L., Xiao, L. P., Ge, L., Wang, Y. Y., et al. Enrichment of MTHFR 677 T in a Chinese long-lived cohort and its association with lipid modulation. Lipids Health Dis (2014) 13:1–12. doi: 10.1186/1476-511X-13-104, 24968810
40. Golja, MV, Šmid, A, Kuželički, NK, Trontelj, J, Geršak, K, and Mlinarič-Raščan, I. Folate insufficiency due to MTHFR deficiency is bypassed by 5-methyltetrahydrofolate. J Clin Med. (2020) 9:1–18. doi: 10.3390/jcm9092836,
41. Colson, NJ, Naug, HL, Nikbakht, E, Zhang, P, and McCormack, J. The impact of MTHFR 677 C/T genotypes on folate status markers: a meta-analysis of folic acid intervention studies. Eur J Nutr. (2017) 56:247–60. doi: 10.1007/s00394-015-1076-x,
42. Luo, Z, Lu, Z, Muhammad, I, Chen, Y, Chen, Q, Zhang, J, et al. Associations of the MTHFR rs1801133 polymorphism with coronary artery disease and lipid levels: a systematic review and updated meta-analysis. Lipids Health Dis. (2018) 17:1–15. doi: 10.1186/s12944-018-0837-y,
43. Villar, C, Campistol, J, Fons, C, Armstrong, J, Mas, A, Ormazabal, A, et al. Glycine and l-arginine treatment causes hyperhomocysteinemia in cerebral creatine transporter deficiency patients. JIMD Rep. (2011) 4:13–6. doi: 10.1007/8904_2011_41
44. Recber, T, Orgul, G, Aydın, E, Tanacan, A, Nemutlu, E, Kır, S, et al. Metabolic infrastructure of pregnant women with methylenetetrahydrofolate reductase polymorphisms: a metabolomic analysis. Biomed Chromatogr. (2020) 34:e4842–7. doi: 10.1002/bmc.4842,
45. Calderón-Larrañaga, A, Saadeh, M, Hooshmand, B, Refsum, H, Smith, AD, Marengoni, A, et al. Association of homocysteine, methionine, and MTHFR 677C>T polymorphism with rate of cardiovascular multimorbidity development in older adults in Sweden. JAMA Netw Open. (2020) 3:e205316. doi: 10.1001/jamanetworkopen.2020.5316,
46. Bortolus, R, Filippini, F, Chiaffarino, F, Udali, S, Rinaldi, M, Gandini, G, et al. Lifestyle and dietary factors, iron status and one-carbon metabolism polymorphisms in a sample of Italian women and men attending a transfusion medicine unit: a cross-sectional study. Br J Nutr. (2023) 130:65–70. doi: 10.1017/S0007114522003245,
47. Liu, M, Yu, J, Su, Z, Sun, Y, Liu, Y, Xie, Q, et al. Associations between prenatal exposure to cadmium and lead with neural tube defect risks are modified by single-nucleotide polymorphisms of fetal MTHFR and SOD2: a case–control study. Environ Health. (2021) 20:1–10. doi: 10.1186/s12940-021-00752-9,
48. Parajuli, RP, Goodrich, JM, Chan, HM, Lemire, M, Ayotte, P, Hegele, RA, et al. Variation in biomarker levels of metals, persistent organic pollutants, and omega-3 fatty acids in association with genetic polymorphisms among Inuit in Nunavik, Canada. Environ Res. (2021) 200:111393. doi: 10.1016/j.envres.2021.111393,
49. Saad-Hussein, A, Shousha, WG, Al-Sadek, SYM, and Ramadan, SS. Role of MTHFR 677C>T and 1298A>C gene polymorphisms on renal toxicity caused by lead exposure in wastewater treatment plant workers. Environ Sci Pollut Res. (2023) 30:84758–64. doi: 10.1007/s11356-023-28309-y,
50. EVU (Europen Vegetarian Union). EVU position paper. Definitions of “vegan” and “vegetarian” in accordance with the EU food information regulation. (2019). Available online at: https://www.euroveg.eu/wp-content/uploads/2018/08/EVU_PP_Definition.pdf (Accessed May 12, 2025).
51. Breier, M, Wahl, S, Prehn, C, Fugmann, M, Ferrari, U, Weise, M, et al. Targeted metabolomics identifies reliable and stable metabolites in human serum and plasma samples. PLoS One. (2014) 9:1–11. doi: 10.1371/journal.pone.0089728,
52. Sekovanić, A, Jurasović, J, Piasek, M, Pašalić, D, Orct, T, Grgec, AS, et al. Metallothionein 2A gene polymorphism and trace elements in mother-newborn pairs in the Croatian population. J Trace Elem Med Biol. (2018) 45:163–70. doi: 10.1016/j.jtemb.2017.10.011,
53. Badawy, AAB, Morgan, CJ, and Turner, JA. Application of the Phenomenex EZ:faast™ amino acid analysis kit for rapid gas-chromatographic determination of concentrations of plasma tryptophan and its brain uptake competitors. Amino Acids. (2008) 34:587–96. doi: 10.1007/s00726-007-0012-7,
54. Zappacosta, B, Romano, L, Persichilli, S, Cutrone, LA, Graziano, M, Vitrani, A, et al. Genotype prevalence and allele frequencies of 5,10-methylenetetrahydrofolate reductase (MTHFR) C677T and A1298C polymorphisms in Italian newborns. Lab Med. (2009) 40:732–6. doi: 10.1309/LMBES080OMMANDOZ
55. Helsel, DR. Statistics for censored environmental data using Minitab® and R. 2nd ed. Hoboken, NJ, USA: John Wiley & Sons, Inc. (2011).
56. Boss, J, Mukherjee, B, Ferguson, KK, Aker, A, Alshawabkeh, AN, Cordero, JF, et al. Estimating outcome-exposure associations when exposure biomarker detection limits vary across batches. Epidemiology. (2019) 30:746–55. doi: 10.1097/EDE.0000000000001052,
57. Clemmensen, L, Hastie, T, Witten, D, and Ersbøll, B. Sparse discriminant analysis. Technometrics. (2011) 53:406–13. doi: 10.1198/TECH.2011.08118
58. Cohen, J. Statistical power analysis for the behavioural sciences. New York: Academic Press (1969).
59. Faul, F, Erdfelder, E, Lang, A-G, and Buchner, A. G*power 3: a flexible statistical power analysis program for the social, behavioral, and biomedical sciences. Behav Res Methods. (2007) 39:175–91. doi: 10.3758/BF03193146,
60. Koch, CA, Kjeldsen, EW, and Frikke-Schmidt, R. Vegetarian or vegan diets and blood lipids: a meta-analysis of randomized trials. Eur Heart J. (2023) 44:2609–22. doi: 10.1093/eurheartj/ehad211,
61. Borén, J, John Chapman, M, Krauss, RM, Packard, CJ, Bentzon, JF, Binder, CJ, et al. Low-density lipoproteins cause atherosclerotic cardiovascular disease: pathophysiological, genetic, and therapeutic insights: a consensus statement from the European atherosclerosis society consensus panel. Eur Heart J. (2020) 41:2313–30. doi: 10.1093/eurheartj/ehz962,
62. Schade, DS, Helitzer, D, and Eaton, P. Evidence that low density lipoprotein is the primary cause of atherosclerotic cardiovascular disease: a Bradford-hill approach. World J Cardiovasc Dis. (2017) 7:271–84. doi: 10.4236/wjcd.2017.79025
63. Hunt, JR. Bioavailability of iron, zinc, and other trace minerals from vegetarian diets. Am J Clin Nutr. (2003) 78:633S–9S. doi: 10.1093/ajcn/78.3.633s,
64. Mattila, P, Mäkinen, S, Eurola, M, Jalava, T, Pihlava, JM, Hellström, J, et al. Nutritional value of commercial protein-Rich Plant products. Plant Foods Hum Nutr. (2018) 73:108–15. doi: 10.1007/s11130-018-0660-7,
65. Phang, JM. Proline metabolism in cell regulation and cancer biology: recent advances and hypotheses. Antioxidants Redox Signal. (2019) 30:635–49. doi: 10.1089/ars.2017.7350,
66. Zhang, J, Fan, J, Venneti, S, Cross, JR, Takagi, T, Bhinder, B, et al. Asparagine plays a critical role in regulating cellular adaptation to glutamine depletion. Mol Cell. (2014) 56:205–18. doi: 10.1016/j.molcel.2014.08.018,
67. James, SC, Fraser, K, Cooney, J, Günther, CS, Young, W, Gearry, RB, et al. Concentrations of plasma amino acids and neurotransmitters in participants with functional gut disorders and healthy controls. Meta. (2023) 13:313. doi: 10.3390/metabo13020313,
68. Tanas, A, Tozlu, ÖÖ, Gezmiş, T, Hacimüftüoǧlu, A, Abd El-Aty, AM, Ceylan, O, et al. In vitro and in vivo neuroprotective effects of sarcosine. Biomed Res Int. (2022) 2022:7498. doi: 10.1155/2022/5467498,
69. Bates, CJ, Prentice, A, Birch, MC, and Delves, HT. Dependence of blood indices of selenium and mercury on estimated fish intake in a national survey of British adults. Public Health Nutr. (2007) 10:508–17. doi: 10.1017/S1368980007246683,
70. Bai, S, Zhang, M, Tang, S, Li, M, Wu, R, Wan, S, et al. Effects and impact of selenium on human health, a review. Molecules. (2025) 30:50. doi: 10.3390/molecules30010050
71. Lin, X, Li, S, Zou, Y, Zhao, FQ, Liu, J, and Liu, H. Lysine stimulates protein synthesis by promoting the expression of ATB0,+ and activating the mTOR pathway in bovine mammary epithelial cells. J Nutr. (2018) 148:1426–33. doi: 10.1093/jn/nxy140,
72. Yi, X, Yang, Y, Li, T, Li, M, Yao, T, Hu, G, et al. Signaling metabolite β-aminoisobutyric acid as a metabolic regulator, biomarker, and potential exercise pill. Front Endocrinol (Lausanne). (2023) 14:1–12. doi: 10.3389/fendo.2023.1192458,
73. Aoyagi, S, Baker, DH, and Wedekind, KJ. Estimates of copper bioavailability from liver of different animal species and from feed ingredients derived from plants and animals. Poult Sci. (1993) 72:1746–55. doi: 10.3382/ps.0721746,
74. European Food Safety Authority (EFSA). Scientific opinion on dietary reference values for molybdenum. EFSA J. (2013) 11:1–35. doi: 10.2903/j.efsa.2013.3333
75. Melnyk, LJ, Donohue, MJ, Pham, M, and Donohue, J. Absorption of strontium by foods prepared in drinking water. J Trace Elem Med Biol. (2019) 53:22–6. doi: 10.1016/j.jtemb.2019.01.001,
76. Botto, LD, and Yang, Q. 5,10-methylenetetrahydrofolate reductase gene variants and congenital anomalies: a HuGE review. Am J Epidemiol. (2000) 151:862–77. doi: 10.1093/oxfordjournals.aje.a010290,
77. Garcı, H, Sa, E, Contreras-cubas, C, Barajas-olmos, F, Cid, M, Centeno-cruz, F, et al. Heterogenous distribution of MTHFR gene variants among mestizos and diverse Amerindian Groups from Mexico. PLoS One. (2016) 11:1–13. doi: 10.1371/journal.pone.0163248
78. Lovričević, I, Franjić, BD, Tomičić, M, Vrkić, N, De Syo, D, Hudorović, N, et al. 5,10-methylenetetrahydrofolate reductase (MTHFR) 677 C→T genetic polymorphism in 228 croatian volunteers. Coll Antropol. (2004) 28:647–54.
79. Zec Sambol, S, Obrovac Glišić, M, Marković Bašić, N, Dvornik, Š, Grahovac, B, Skočić Mihić, S, et al. Influence of dietary pattern and methylentetrahydrofolate reductase C677T polymorphism on the plasma homocysteine level among healthy vegetarians and omnivores. Hrana u Zdr i Boles Znan časopis za Nutr i dijetetiku. (2017) 6:54–62.
80. Žuntar, I, Topić, E, Vukosavić, D, Vuković, V, Demarin, V, Begonja, A, et al. Croatian population data for the C677T polymorphism in methylenetetrahydrofolate reductase: frequencies in healthy and atherosclerotic study groups. Clin Chim Acta. (2003) 335:95–100. doi: 10.1016/S0009-8981(03)00283-3
81. Reljic, A, Simundic, AM, Topic, E, Nikolac, N, Justinic, D, and Stefanovic, M. The methylenetetrahydrofolate reductase (MTHFR) C677T polymorphism and cancer risk: the Croatian case-control study. Clin Biochem. (2007) 40:981–5. doi: 10.1016/j.clinbiochem.2007.05.005,
82. Herrmann, W, Obeid, R, Schorr, H, Zarzour, W, and Geisel, J. Homocysteine, methylenetetrahydrofolate reductase C677T polymorphism and the B-vitamins: a facet of nature-nurture interplay. Clin Chem Lab Med. (2003) 41:547–53. doi: 10.1515/CCLM.2003.083,
83. Li, WX, Dai, SX, Zheng, JJ, Liu, JQ, and Huang, JF. Homocysteine metabolism gene polymorphisms (MTHFR C677T, MTHFR A1298C, MTR A2756G and MTRR A66G) jointly elevate the risk of folate deficiency. Nutrients. (2015) 7:6670–87. doi: 10.3390/nu7085303,
85. Ledowsky, CJ, Schloss, J, and Steel, A. Variations in folate prescriptions for patients with the MTHFR genetic polymorphisms: a case series study. Explor Res Clin Soc Pharm. (2023) 10:100277. doi: 10.1016/j.rcsop.2023.100277,
86. Rozen, R. Molecular biology of methylenetetrahydrofolate reductase (MTHFR): interrelationships with folic acid, homocysteine and vascular disease BT - homocysteine and vascular disease In: K Robinson, editor. Homocysteine and vascular disease. Developments in cardiovascular medicine. Dordrecht: Springer Netherlands (2000). 271–89.
87. Friso, S, Choi, SW, Girelli, D, Mason, JB, Dolnikowski, GG, Bagley, PJ, et al. A common mutation in the 5,10-methylenetetrahydrofolate reductase gene affects genomic DNA methylation through an interaction with folate status. Proc Natl Acad Sci USA. (2002) 99:5606–11. doi: 10.1073/pnas.062066299,
88. Pereira, AC, Schettert, IT, Morandini Filho, AAF, Guerra-Shinohara, EM, and Krieger, JE. Methylenetetrahydrofolate reductase (MTHFR) c677t gene variant modulates the homocysteine folate correlation in a mild folate-deficient population. Clin Chim Acta. (2004) 340:99–105. doi: 10.1016/j.cccn.2003.09.016,
89. Stanger, O. Physiology of folic acid in health and disease. Curr Drug Metab. (2002) 3:211–23. doi: 10.2174/1389200024605163,
90. Tisman, G, and Herbert, V. B 12 dependence of cell uptake of serum folate: an explanation for high serum folate and cell folate depletion in B 12 deficiency. Blood. (1973) 41:465–9. doi: 10.1182/blood.v41.3.465.465
91. Baszczuk, A, Thielemann, A, Musialik, K, Kopczyński, J, Bielawska, L, Dżumak, A, et al. The impact of supplementation with folic acid on homocysteine concentration and selected lipoprotein parameters in patients with primary hypertension. J Nutr Sci Vitaminol (Tokyo). (2017) 63:96–103. doi: 10.3177/jnsv.63.96,
92. Visram, M, Radulovic, M, Steiner, S, Malanovic, N, Eichmann, TO, Wolinski, H, et al. Homocysteine regulates fatty acid and lipid metabolism in yeast. J Biol Chem. (2018) 293:5544–55. doi: 10.1074/jbc.M117.809236,
93. Llanos, AAM, Marian, C, Brasky, TM, Dumitrescu, RG, Liu, Z, Mason, JB, et al. Associations between genetic variation in onecarbon metabolism and LINE-1 DNA methylation in histologically normal breast tissues. Epigenetics. (2015) 10:727–35. doi: 10.1080/15592294.2015.1062205,
94. Jiang, S, Zhao, R, Pan, M, Venners, SA, Zhong, G, and Hsu, YH. Associations of MTHFR and MTRR polymorphisms with serum lipid levels in Chinese hypertensive patients. Clin Appl Thromb Hemost. (2014) 20:400–10. doi: 10.1177/1076029612467226,
95. Zhang, L, Yin, R-X, Liu, W-Y, Miao, L, Wu, D-F, Aung, LHH, et al. Association of methylenetetrahydrofolate reductase C677T polymorphism and serum lipid levels in the Guangxi Bai Ku Yao and Han populations. Lipids Health Dis. (2010) 9:123. doi: 10.1186/1476-511X-9-123,
96. Woo, CWH, Siow, YL, Pierce, GN, Choy, PC, Minuk, GY, Mymin, D, et al. Hyperhomocysteinemia induces hepatic cholesterol biosynthesis and lipid accumulation via activation of transcription factors. Am J Physiol Endocrinol Metab. (2005) 288:1002–10. doi: 10.1152/ajpendo.00518.2004,
97. Werstuck, GH, Lentz, SR, Dayal, S, Hossain, GS, Sood, SK, Shi, YY, et al. Homocysteine-induced endoplasmic reticulum stress causes dysregulation of the cholesterol and triglyceride biosynthetic pathways. J Clin Invest. (2001) 107:1263–73. doi: 10.1172/JCI11596,
98. Guan, D, Ji, Y, Lu, X, Feng, W, and Ma, W. Associations of MTHFR gene polymorphism with lipid metabolism and risk of cerebral infarction in the Northwest Han Chinese population. Front Neurol. (2023) 14:2351. doi: 10.3389/fneur.2023.1152351,
99. Geisel, J, Hübner, U, Bodis, M, Schorr, H, Knapp, J, Obeid, R, et al. The role of genetic factors in the development of hyperhomocysteinemia. Clin Chem Lab Med. (2003) 41:1427–34. doi: 10.1515/CCLM.2003.219,
100. Steluti, J, Carvalho, AM, Carioca, AAF, Miranda, A, Gattás, GJF, Fisberg, RM, et al. Genetic variants involved in one-carbon metabolism: polymorphism frequencies and differences in homocysteine concentrations in the folic acid fortification era. Nutrients. (2017) 9:1–12. doi: 10.3390/nu9060539,
101. Jin, P, Bian, Y, Wang, K, Cong, G, Yan, R, Sha, Y, et al. Homocysteine accelerates atherosclerosis via inhibiting LXRα–mediated ABCA1/ABCG1–dependent cholesterol efflux from macrophages. Life Sci. (2018) 214:41–50. doi: 10.1016/j.lfs.2018.10.060,
102. Koklesova, L, Mazurakova, A, Samec, M, Biringer, K, Samuel, SM, Büsselberg, D, et al. Homocysteine metabolism as the target for predictive medical approach, disease prevention, prognosis, and treatments tailored to the person. EPMA J. (2021) 12:477–505. doi: 10.1007/s13167-021-00263-0,
103. Obeid, R, and Herrmann, W. Homocysteine and lipids: S-adenosyl methionine as a key intermediate. FEBS Lett. (2009) 583:1215–25. doi: 10.1016/j.febslet.2009.03.038,
104. Kluck, GEG, Durham, KK, Yoo, JA, and Trigatti, BL. High density lipoprotein and its precursor protein apolipoprotein A1 as potential therapeutics to prevent anthracycline associated cardiotoxicity. Front Cardiovasc Med. (2020) 7:1–16. doi: 10.3389/fcvm.2020.00065,
105. Cicero, AFG, Fogacci, F, Stoian, AP, Vrablik, M, Al Rasadi, K, Banach, M, et al. Nutraceuticals in the management of dyslipidemia: which, when, and for whom? Could nutraceuticals help low-risk individuals with non-optimal lipid levels? Curr Atheroscler Rep. (2021) 23:57. doi: 10.1007/s11883-021-00955-y,
106. Garoufi, A, Papadaki, M, Kalogiannis, M, Zerva, U, Charakida, M, Marmarinos, A, et al. Beneficial effects of a nutraceutical combination on lipid profiles in children with moderate and severe hypercholesterolemia. Biomolecules. (2024) 14:1–13. doi: 10.3390/biom14121608,
107. Ford, ES, and Bowman, BA. Serum and red blood cell folate concentrations, race, and education: findings from the third National Health and nutrition examination survey. Am J Clin Nutr. (1999) 69:476–81. doi: 10.1093/ajcn/69.3.476,
108. Hedayatnia, M, Asadi, Z, Zare-Feyzabadi, R, Yaghooti-Khorasani, M, Ghazizadeh, H, Ghaffarian-Zirak, R, et al. Dyslipidemia and cardiovascular disease risk among the MASHAD study population. Lipids Health Dis. (2020) 19:1–11. doi: 10.1186/s12944-020-01204-y
109. Yin, RX, Wu, DF, Miao, L, Aung, LHH, Cao, XL, Yan, TT, et al. Several genetic polymorphisms interact with overweight/obesity to influence serum lipid levels. Cardiovasc Diabetol. (2012) 11:123–18. doi: 10.1186/1475-2840-11-123,
110. Li, WX, Lv, WW, Dai, SX, Pan, ML, and Huang, JF. Joint associations of folate, homocysteine and MTHFR, MTR and MTRR gene polymorphisms with dyslipidemia in a Chinese hypertensive population: a cross-sectional study. Lipids Health Dis. (2015) 14:101–11. doi: 10.1186/s12944-015-0099-x,
111. Frelut, ML, Emery-Fillon, N, Guilland, JC, Dao, HH, and De Courcy, GP. Alanine amino transferase concentrations are linked to folate intakes and methylenetetrahydrofolate reductase polymorphism in obese adolescent girls. J Pediatr Gastroenterol Nutr. (2006) 43:234–9. doi: 10.1097/01.mpg.0000228110.83616.92,
112. Jain, M, Pandey, P, Tiwary, N, and Jain, S. MTHFR C677T polymorphism is associated with hyperlipidemia in women with polycystic ovary syndrome. J Hum Reprod Sci. (2012) 5:52–6. doi: 10.4103/0974-1208.97802,
113. Riancho, JA, Valero, C, and Zarrabeitia, MT. MTHFR polymorphism and bone mineral density: Meta-analysis of published studies. Calcif Tissue Int. (2006) 79:289–93. doi: 10.1007/s00223-006-0143-y,
114. Wang, H, and Liu, C. Association of MTHFR C667T polymorphism with bone mineral density and fracture risk: an updated meta-analysis. Osteoporos Int. (2012) 23:2625–34. doi: 10.1007/s00198-011-1885-6,
115. Dalemo, S, Eggertsen, R, Hjerpe, P, Almqvist, EG, and Boström, KB. Bone mineral density in primary care patients related to serum calcium concentrations: a longitudinal cohort study from Sweden. Scand J Prim Health Care. (2018) 36:198–206. doi: 10.1080/02813432.2018.1459430,
116. Reid, IR, Bristow, SM, and Bolland, MJ. Calcium and cardiovascular disease. Endocrinol Metab. (2017) 32:339–49. doi: 10.3803/EnM.2017.32.3.339,
117. Casas, JP, Bautista, LE, Smeeth, L, Sharma, P, and Hingorani, AD. Homocysteine and stroke: evidence on a causal link from mendelian randomisation. Lancet. (2005) 365:224–32. doi: 10.1016/s0140-6736(05)17742-3,
118. Sadiku, OO, and Rodríguez-Seijo, A. Metabolic and genetic derangement: a review of mechanisms involved in arsenic and lead toxicity and genotoxicity. Arh Hig Rada Toksikol. (2022) 73:244–55. doi: 10.2478/aiht-2022-73-3669,
119. EFSA (CONTAM). Scientific opinion on arsenic in food. EFSA J. (2009) 7:1351. doi: 10.2903/j.efsa.2009.1351
120. Martinez, VD, Vucic, EA, Becker-Santos, DD, Gil, L, and Lam, WL. Arsenic exposure and the induction of human cancers. J Toxicol. (2011) 2011:1–13. doi: 10.1155/2011/431287,
121. Wei, S, Zhang, H, and Tao, S. A review of arsenic exposure and lung cancer. Toxicol Res (Camb). (2019) 8:319–27. doi: 10.1039/c8tx000298c,
122. Lindberg, AL, Kumar, R, Goessler, W, Thirumaran, R, Gurzau, E, Koppova, K, et al. Metabolism of low-dose inorganic arsenic in a central European population: influence of sex and genetic polymorphisms. Environ Health Perspect. (2007) 115:1081–6. doi: 10.1289/ehp.10026,
123. Engström, KS, Broberg, K, Concha, G, Nermell, B, Warholm, M, and Vahter, M. Genetic polymorphisms influencing arsenic metabolism: evidence from Argentina. Environ Health Perspect. (2007) 115:599–605. doi: 10.1289/ehp.9734,
124. Shen, W, Zhang, B, Liu, S, Wu, H, Gu, X, Qin, L, et al. Association of blood lead levels with methylenetetrahydrofolate reductase polymorphisms among Chinese pregnant women in Wuhan City. PLoS One. (2015) 10:1–13. doi: 10.1371/journal.pone.0117366,
125. Golbahar, J, Hamidi, A, Aminzadeh, MA, and Omrani, GR. Association of plasma folate, plasma total homocysteine, but not methylenetetrahydrofolate reductase C667T polymorphism, with bone mineral density in postmenopausal Iranian women: a cross-sectional study. Bone. (2004) 35:760–5. doi: 10.1016/j.bone.2004.04.018,
126. Sekovanić, A, Piasek, M, Orct, T, Grgec, AS, Sarić, MM, Stasenko, S, et al. Mercury exposure assessment in mother–infant pairs from continental and coastal Croatia. Biomolecules. (2020) 10:19–22. doi: 10.3390/biom10060821,
127. Grgec, AS, Jurasović, J, Kljaković-Gašpić, Z, Orct, T, Samarin, IR, Janči, T, et al. Potential risks and health benefits of fish in the diet during the childbearing period: focus on trace elements and n-3 fatty acid content in commonly consumed fish species from the Adriatic Sea. Environ Adv. (2022) 8:100226. doi: 10.1016/j.envadv.2022.100226
128. Canada H. Second report on human biomonitoring of environmental Chemicals in Canada: results of the Canadian health measures survey cycle 1 (2007–2009). Ontario, Canada: Health Canada (2013). 71 p.
129. De Oliveira, RB, Da Silva, DM, Franco, TSBS, Vasconcelos, CRS, De Sousa, DJDA, Sarrazin, SLF, et al. Fish consumption habits of pregnant women in Itaituba, Tapajós River basin, Brazil: risks of mercury contamination as assessed by measuring total mercury in highly consumed piscivore fish species and in hair of pregnant women. Arh Hig Rada Toksikol. (2022) 73:131–42. doi: 10.2478/aiht-2022-73-3611,
130. EFSA NDA Panel. Scientific opinion on health benefits of seafood (fish and shellfish) consumption in relation to health risks associated with exposure to methylmercury. EFSA J. (2014) 12:3761. doi: 10.2903/j.efsa.2014.3761
131. FDA. Mercury in food. (2024). Available online at: https://www.fda.gov/food/environmental-contaminants-food/mercury-food (accessed March 7, 2025).
132. Gajski, G, Matković, K, Delić, L, and Gerić, M. Evaluation of primary DNA damage in Young healthy females based on their dietary preferences. Nutrients. (2023) 15:2218. doi: 10.3390/nu15092218,
133. Fiore, M, Barone, R, Copat, C, Grasso, A, Cristaldi, A, Rizzo, R, et al. Metal and essential element levels in hair and association with autism severity. J Trace Elem Med Biol. (2020) 57:126409. doi: 10.1016/j.jtemb.2019.126409,
134. Wahil, MSA, Ja’afar, MH, and Isa, ZM. Urinary aluminium and its association with autism spectrum disorder in urban preschool children in Malaysia. PeerJ. (2023) 11:1–16. doi: 10.7717/PEERJ.15132
135. Rashaid, AHB, Nusair, SD, Alqhazo, MT, Adams, JB, Abu-Dalo, MA, and Bashtawi, MA. Heavy metals and trace elements in scalp hair samples of children with severe autism spectrum disorder: a case-control study on Jordanian children. J Trace Elem Med Biol. (2021) 67:126790. doi: 10.1016/j.jtemb.2021.126790,
136. Rai, V. Association of methylenetetrahydrofolate reductase (MTHFR) gene C677T polymorphism with autism: evidence of genetic susceptibility. Metab Brain Dis. (2016) 31:727–35. doi: 10.1007/s11011-016-9815-0,
137. Sadeghiyeh, T, Dastgheib, SA, Mirzaee-Khoramabadi, K, Morovati-Sharifabad, M, Akbarian-Bafghi, MJ, Poursharif, Z, et al. Association of MTHFR 677C>T and 1298A>C polymorphisms with susceptibility to autism: a systematic review and meta-analysis. Asian J Psychiatr. (2019) 46:54–61. doi: 10.1016/j.ajp.2019.09.016
138. Bai, D, Yip, BHK, Windham, GC, Sourander, A, Francis, R, Yoffe, R, et al. Association of genetic and environmental factors with autism in a 5-country cohort. JAMA Psychiatry. (2019) 76:1035–43. doi: 10.1001/jamapsychiatry.2019.1411,
139. Fell, DA. Enzymes, metabolites and fluxes. J Exp Bot. (2005) 56:267–72. doi: 10.1093/jxb/eri011,
140. Krumsiek, J, Suhre, K, Evans, AM, Mitchell, MW, Mohney, RP, Milburn, MV, et al. Mining the unknown: a systems approach to metabolite identification combining genetic and metabolic information. PLoS Genet. (2012) 8:e1003005. doi: 10.1371/journal.pgen.1003005,
141. Beagle, B, Yang, TL, Hung, J, Cogger, EA, Moriarty, DJ, and Caudill, MA. The glycine N-methyltransferase (GNMT) 1289 C→T variant influences plasma total homocysteine concentrations in young women after restricting folate intake. J Nutr. (2005) 135:2780–5. doi: 10.1093/jn/135.12.2780,
142. Petr, M, Šteffl, M, and Kohlíková, E. Effect of the MTHFR 677c/t polymorphism on homocysteinemia in response to creatine supplementation: a case study. Physiol Res. (2013) 62:721–9. doi: 10.33549/physiolres.932542
143. Wagner, C, Briggs, WT, and Cook, RJ. Inhibition of glycine n-methyltransferase activity by folate derivatives: implications for regulation of methyl group metabolism. Biochem Biophys Res Commun. (1985) 127:746–52. doi: 10.1016/S0006-291X(85)80006-1,
144. Simile, MM, Latte, G, Feo, CF, Feo, F, Calvisi, DF, and Pascale, RM. Alterations of methionine metabolism in hepatocarcinogenesis: the emergent role of glycine N-methyltransferase in liver injury. Ann Gastroenterol. (2018) 31:552–60. doi: 10.20524/aog.2018.0288,
145. Wang, YC, Wu, MT, Lin, YJ, Tang, FY, Ko, HA, and Chiang, EP. Regulation of folate-mediated one-carbon metabolism by glycine N-methyltransferase (GNMT) and methylenetetrahydrofolate reductase (MTHFR). J Nutr Sci Vitaminol (Tokyo). (2015) 61:S148–50. doi: 10.3177/jnsv.61.S148,
146. Ward, KM, Burghardt, K, Kraal, AZ, Jaeger, A, Yeomans, L, McHugh, C, et al. Genetic and metabolite variability in one-carbon metabolism applied to an insulin resistance model in patients with schizophrenia receiving atypical antipsychotics. Front Psych. (2021) 12:1–8. doi: 10.3389/fpsyt.2021.623143,
147. Chen, J, Cui, L, Lu, S, and Xu, S. Amino acid metabolism in tumor biology and therapy. Cell Death Dis. (2024) 15:42. doi: 10.1038/s41419-024-06435-w,
148. Phillips, F. Vegetarian nutrition. Nutr Bull. (2005) 30:132–67. doi: 10.1111/j.1467-3010.2005.00467.x
149. Smith, W. Does gender influence online survey participation? A record-linkage analysis of university faculty online survey response behavior. (2008). Available online at: https://files.eric.ed.gov/fulltext/ED501717.pdf (accessed June 2, 2025).
150. Lobato, L, Bethony, JM, Pereira, FB, Grahek, SL, Diemert, D, and Gazzinelli, MF. Impact of gender on the decision to participate in a clinical trial: a cross-sectional study. BMC Public Health. (2014) 14:1156. doi: 10.1186/1471-2458-14-1156
151. Kumar, J, Garg, G, Kumar, A, Sundaramoorthy, E, Sanapala, KR, Ghosh, S, et al. Single nucleotide polymorphisms in homocysteine metabolism pathway genes: association of CHDH A119C and MTHFR C677T with hyperhomocysteinemia. Circ Cardiovasc Genet. (2009) 2:599–606. doi: 10.1161/CIRCGENETICS.108.841411,
Keywords: single nucleotide polymorphism, recessive genotype, toxic and essential elements, folate, proline, glutamic acid
Citation: Orct T, Kovačić J, Peremin I, Kljaković-Gašpić Z, Pašalić D, Sekovanić A, Dorotić A, Tariba Lovaković B, Jurič A, Pizent A, Crnjac F, Dukić L, Gerić M, Jakaša I and Gajski G (2026) MTHFR 677C/T gene polymorphism and dietary habits: effects on trace element levels, amino acids, and biochemical parameters. Front. Nutr. 12:1710613. doi: 10.3389/fnut.2025.1710613
Edited by:
L. Joseph Su, University of Texas Southwestern Medical Center, United StatesReviewed by:
Humberto Garcia-Ortiz, National Institute of Genomic Medicine (INMEGEN), MexicoHuoxiang Zhou, Zhengzhou University, China
Copyright © 2026 Orct, Kovačić, Peremin, Kljaković-Gašpić, Pašalić, Sekovanić, Dorotić, Tariba Lovaković, Jurič, Pizent, Crnjac, Dukić, Gerić, Jakaša and Gajski. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ankica Sekovanić, YXNla292YW5pY0BpbWkuaHI=; Goran Gajski, Z2dhanNraUBpbWkuaHI=
†These authors have contributed equally to this work and share first authorship
Daria Pašalić3