- About this Journal
- Abstracting and Indexing
- Aims and Scope
- Annual Issues
- Article Processing Charges
- Articles in Press
- Author Guidelines
- Bibliographic Information
- Citations to this Journal
- Contact Information
- Editorial Board
- Editorial Workflow
- Free eTOC Alerts
- Publication Ethics
- Reviewers Acknowledgment
- Submit a Manuscript
- Subscription Information
- Table of Contents
BioMed Research International
Volume 2014 (2014), Article ID 560183, 9 pages
The Effect of Multiple Single Nucleotide Polymorphisms in the Folic Acid Pathway Genes on Homocysteine Metabolism
1Obstetrics and Gynecology Hospital, Fudan University, Shanghai 200011, China
2Department of Obstetrics and Gynecology, Shanghai Medical College, Fudan University, Shanghai 200032, China
3Shanghai Key Laboratory of Female Reproductive Endocrine Related Diseases, Shanghai 200011, China
4State Key Laboratory of Genetic Engineering, Institute of Biostatistics, School of Life Science, Fudan University, Shanghai 200433, China
5Children’s Hospital of Fudan University, Shanghai 201102, China
6Institute of Biomedical Sciences, Fudan University, Shanghai 201508, China
7Department of Pharmaceutical Sciences, College of Pharmacy, University of South Florida, Tampa, FL 33612, USA
8Central Lab, Shanghai Xuhui Central Hospital, Shanghai 200031, China
Received 5 April 2013; Revised 29 September 2013; Accepted 17 October 2013; Published 12 January 2014
Academic Editor: Wen-Hwa Lee
Copyright © 2014 Shuang Liang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Objective. To investigate the joint effects of the single nucleotide polymorphisms (SNPs) of genes in the folic acid pathway on homocysteine (Hcy) metabolism. Methods. Four hundred women with normal pregnancies were enrolled in this study. SNPs were identified by MassARRAY. Serum folic acid and Hcy concentration were measured. Analysis of variance (ANOVA) and support vector machine (SVM) regressions were used to analyze the joint effects of SNPs on the Hcy level. Results. SNPs of MTHFR (rs1801133 and rs3733965) were significantly associated with maternal serum Hcy level. In the different genotypes of MTHFR (rs1801133), SNPs of RFC1 (rs1051266), TCN2 (rs9606756), BHMT (rs3733890), and CBS (rs234713 and rs2851391) were linked with the Hcy level adjusted for folic acid concentration. The integrated SNPs scores were significantly associated with the residual Hcy concentration (RHC) (). The Hcy level was significantly higher in the group with high SNP scores than that in other groups with SNP scores of less than 0.2 (). Moreover, this difference was even more significant in moderate and high levels of folic acid. Conclusion. SNPs of genes in the folic acid pathway possibly affect the Hcy metabolism in the presence of moderate and high levels of folic acid.
The folic acid pathway is essential for hundreds of intracellular transmethylation reactions including DNA methylation and DNA synthesis, processes that are closely related to homocysteine (Hcy) metabolism [1, 2]. Folic acid deficiency and abnormal metabolism of folic acid and Hcy not only play an important role in neural tube defects (NTDs) , but also are key factors for congenital heart disease, cleft lip and palate, late pregnancy complications, premature labor different kinds of neurodegenerative and psychiatric diseases, and cancer [3–10]. It was reported that folic acid preconceptional supplementation is effective for NTD prevention. However, it remains unclear whether 30–50% of cases are still unpreventable with folic acid supplementation .
Variations in genes that play key roles in the folic acid cycle have been widely investigated, where single nucleotide polymorphisms (SNPs) have been found to be associated with folic acid and Hcy metabolism. The variants of MTHFR and RFC1 were found to interact with Hcy levels [11, 12], while the combined effects of the 2 MTHFR polymorphisms (rs1801133 and rs1801131) were found to be associated with Hcy concentration . However, studies have been limited to one or a few SNPs with joint effects [11, 14–16].
With the rapid development of ongoing high-throughput human gene sequencing and bioinformatics, abundant SNPs can be identified. Support vector machines (SVMs) are a classic supervised machine learning algorithm typically used for classification and regression analysis. SVM has been widely used in solving biological problems, including gene function prediction, gene expression data analysis, and even cancer diagnosis [12, 17, 18]. It has also been applied to the analysis of SNPs.
The aim of the study is to identify the genotypes of 18 SNPs using MassARRAY and to investigate an association between the cumulative effects of 18 SNPs in the 9 genes of the folic acid pathway and Hcy metabolism.
2. Material and Method
2.1. Study Population
Four hundred pregnant women at 11–25 gestational weeks were enrolled in the study conducted at the Obstetrics and Gynecology Hospital of Fudan University, China, from April to May 2011. Women included in the study were not smokers, did not drink alcohol, had no chronic diseases, and were not taking prescription medications. Blood samples were taken from fasting subjects. The serum was separated for the measurement of the concentrations of folic acid and Hcy and the remaining blood clots were used for DNA extraction and genotyping. All individuals enrolled in this study signed the informed consent and the study was approved by the Ethics Committee of the Obstetrics and Gynecology hospital of Fudan University, China.
2.2. Biochemical Measurement of Serum Folic Acid and Hcy Concentrations
Serum folic acid concentration was measured by a chemiluminescent microparticle immunoassay (Architect Folic acid Reagent; Abbott, Lisnamuck, Longford) using the ARCHITECT I Systems following the manufacturer’s recommended protocols.
Serum Hcy measurements were carried out by Liquid Chromatography Coupled to Tandem Mass Spectrometry (LC/MS/MS) using an API 3000 LC/MS/MS system (Applied Biosystem) equipped with an electrospray ionization interface that was used in the positive ion mode ([M+H]+) according to the manufacturer’s instructions.
2.3. Identification of SNPs Using MassARRAY
18 SNPs in 9 folic acid pathway (Figure 1) related genes that have been documented in the literature to be associated with Hcy related disease were selected in our population [19–22] based on the CHB data using the following criteria: MAF > 0.1 by the Haploview program (version 4.0) (Table 1).
Genomic DNA was prepared from peripheral leukocytes using Relax Gene blood DNA System (Relax Gene; TIANGEN, Beijing, China). The genotypes were determined by the Sequenom MassARRAY MALDI-TOF system. Primer sequences of the 18 SNPs were shown in Supplementary Table 1 (see Supplementary Table 1 in the Supplementary Material available online at http://dx.doi.org/10.1155/2014/560183).
2.4. SVM Regression Model
SVM regression model was used to analyze the relationship between the SNPs of genes in the folic acid pathway and the changes in serum Hcy concentration. To minimize the effects of folic acid concentration on Hcy concentration, the residual homocysteine concentration (RHC, RHC = actual Hcy concentration − predicted value with a liner regression function of folic acid) was used as the dependent variable. The predicted RHC was defined as the SNP scores. For the independent variables (features) of the SVM regression model, we coded each SNP as two independent variables. For example, if a SNP has three types, AA, Aa, and aa, its coding features would be only two independent variables, which are feature_AA and feature_Aa, following Table 2. Thus, the input space initially has 30 independent variables for 15 SNPs.
Basically, the SVM regression model was implemented using “SMOReg” algorithm of Weka software package with default parameters, where C equals 1.0; the kernel is polynomial kernel with exponent value equals 1.0. Using all of the 15 SNPs to construct a data model would result in suboptimal accuracy because the different variables may contain overlapping information that disturbs the model-constructing process and thus variable selection was needed. Variable selection out of the total 15 SNP variables was conducted using recursive addition and stepwise addition of the input variables. The basic idea of the method was that, beginning with two features, we tried to add features to the SVM regression model which would improve the correlation coefficient most and then add another feature until none of the added features could improve the model or the improvement was less than the 0.01 correlation coefficient (see van Looy et al.’s paper  for more details).
After features selection, a tenfold cross-validation model was used to assess the SVM regression model. In this method, we divided the dataset into 10 subsets of approximately equal size and built the model ten times, each time leaving out one of the subsets as testing set and the others as training sets.
2.5. Statistical Analysis
Statistical analyses were performed using SPSS (SPSS Inc., Chicago, IL, USA), version 16.0 for Windows. A value of <0.05 was considered statistically significant. The Hardy-Weinberg equilibrium constant was assessed using the chi-squared () test. Pairwise linkage disequilibrium of SNPs was estimated using Haploview. The square of the correlation coefficient between markers was used to define linkage using the data from the study population. Linear regression was used for detecting the association between folic acid and Hcy concentration, while analysis of variance (ANOVA) was used for the association analysis between SNP and Hcy concentration with serum Hcy concentration as dependent variable, serum SNP genotypes as the fix factor, and folic acid concentration as covariance.
3.1. Demographic Characteristics of the Participants
Genotyping results were available for 386 of the 400 subjects (96.5%). A total of 14 subjects were excluded: 5 samples due to sample hemolysis and 9 samples due to sample DNA degradation. Data on serum folic acid and Hcy were available for all subjects. The regression equation used was (, , Supplementary Figure 1).
All subjects were genotyped for 18 SNPs. Three SNPs were eliminated from further analysis: TCN2 (rs1801198) showed significant () deviation from the Hardy-Weinberg proportions and CBS (rs5742905) and MTR (rs74767314) were monomorphic. The remaining 15 SNPs are shown in Table 1. The genotyping call rate for each SNP ranged from 96% to 100%. There was no evidence for linkage disequilibrium in our database (, Supplementary Figure 2). There was no significant difference in age, gestational weeks, parities, and pregnancies among each genotype of the 15 SNPs (, Supplementary Tables 2–5).
3.2. The Effects of SNPs on Serum Hcy Concentration Adjusted for Folic Acid Concentration
MTHFR SNPs (rs1801133 and rs3737965) were associated with serum Hcy concentrations which were adjusted for folic acid concentration (Table 3). There were only two cases with homozygous MTHFR SNP (rs3737965). Due to the low frequency of variants for MTHFR SNP (rs3737965) polymorphism, it was not included in the data analysis.
Figure 2 shows the effects of SNPs on Hcy concentration in the different genotypes of MTHFR (rs1801133) after the Hcy concentration was adjusted for folic acid concentration. The SNPs MTHFR (rs1801133) CC, RFC1 (rs1051266), and TCN2 (rs9606756) were significantly associated with Hcy concentration (Figures 2(c) and 2(k)). A similar association was observed with the SNPs CBS (rs2851391) in MTHFR (rs1801133) CT and MTHFR (rs3733890) and CBS (rs234713) in MTHFR (rs1801133) TT (Figures 2(m), 2(i), and 2(n)).
3.3. SVM Model of Multiple SNPs and the Residual Hcy Concentration
In the SVM regression, five SNPs were selected: MTHFR (rs1801133, rs1801131, and rs3737965), CBS (rs234713), and BHMT (rs3733890). The weights of the five SNP variables are shown in Table 4 and the relationship between SNP scores and residual Hcy concentration is shown in Figure 3. The correlation coefficient between RHC and SNP scores was 0.275 in training sets and only 0.247 in the cross-validation combined test sets (Supplementary Table 6).
All subjects were divided into four groups according to the 25%, 50% and 75% of the SNP scores (−0.26, 0, and 0.2, resp.). The Hcy concentration was significantly higher in the group with SNP scores of more than 0.2 than that in groups with SNP scores less than 0.2 (, Figure 4(a)). For those with folic acid levels more than 25% (13.1 ng/mL), a possible interaction between Hcy concentration and SNP scores was detected () (Figures 4(c) and 4(d)). However, for those subjects with folic acid levels less than 25%, SNP scores appeared not to be associated with Hcy concentration (Figure 4(b)).
We first used SVM regression to predict Hcy concentration from the SNPs of genes in the folic acid pathway after analysis of the joint effect between MTHFR (rs1801133) and other genes related to Hcy metabolism. The results revealed that the integrate SNPs scores of SVM were significantly associated with Hcy concentration, especially at moderate and high levels of folic acid. This finding suggests that the variations in the genes of the folic acid pathway may be an important contributor to Hcy related diseases in women with moderate and high folic acid levels from folic acid supplementation.
Our finding that the SNPs in MTHFR (rs1801133, rs1801131, and rs3737965), RFC1 (rs1051266), CBS (rs2851391 and rs234713), TCN2 (rs9606756), and BHMT (rs3733890) were associated with the Hcy level adjusted for folic acid level is partly consistent with previous studies. Moreover, it is worth noting that the variation of MHTFR (rs1801133) was included, along with other SNPs. The thermolabile protein MTHFR is of great importance for the regulation of available 5-MTHF, which serves as the main circulating folate necessary for Hcy remethylation. MTHFR (rs1801133) can result in 50–60% reduced enzyme activity, while MTHFR (rs1801131) can also decrease MTHFR activity [24, 25]. MTHFR (rs3737965) was moderately associated with plasma folic acid concentration by genome-wide association studies . RFC1 is necessary for the uptake of folic acids such as 5-MTHF and RFC1 (rs1051266) was reported to have an impact on folic acid and Hcy concentrations [19, 22]. BHMT (rs3733890), TCN2 (rs9606756), and CBS (rs2851391 and rs234713) have been reported to be associated with Hcy related diseases [27–29].
High-throughput techniques have spawned a mass of complex biological data. However, analysis of these data creates a bottleneck seen in current studies .
In our study, the joint effects of SNPs generated nonlinear, noisy, and complex data sets that also contained a great deal of irrelevant information. Despite the effects of serum folic acid concentration, we have successfully presented a SVM regression model that can evaluate the RHC by five SNPs of genes in the folic acid pathway as inputs. This model translated the complex SNP patterns into a simple output of SNP scores which was significantly related to the changes in Hcy concentration adjusted by serum folic acid concentration, and it was found that 5 out of 15 SNPs were useful as inputs. The RHC was constructed to eliminate the effects of folic acid, for folic acid itself can affect Hcy level. It was suggested that the SVM model could be a potential algorithm for predicting Hcy related diseases.
Furthermore, we found that the Hcy concentration was significantly higher in the group with SNP scores of more than 0.2 than that in groups with SNP scores of less than 0.2, especially for those with folic acid level more than 25%. However, there were no changes in Hcy concentration detected for those subjects with folic acid levels less than 25%. Abnormal metabolism of Hcy is related to many diseases, such as congenital heart disease , cleft lip with or without cleft palate , and NTDs . The causes of these diseases have not been identified under normal concentrations of folic acid although folic acid periconceptional supplementation can effectively prevent many diseases related to folic acid deficiency. Our study provides the evidence that the joint effects of SNPs in the folic acid pathway may play an important role in Hcy related diseases, especially under sufficient support of folic acid.
In conclusion, the joint effects of SNPs in genes that belong to the folic acid pathway can affect Hcy metabolism especially under normal and high levels of folic acid. Further research that includes a bigger sample size is needed to test this SVM model.
|ANOVA:||Analysis of variance|
|CBS:||Cystathionine β synthase|
|LC/MS/MS:||Liquid Chromatography Coupled to Tandem Mass Spectrometry|
|MTRR:||Methionine synthase reductase|
|NTDs:||Neural tube defects|
|RFC1:||Reduced folate carrier 1|
|RHC:||Residual homocysteine concentration|
|SNPs:||Single nucleotide polymorphisms|
|SPSS:||Statistical packages for social sciences|
|SVM:||Support vector machine|
Conflict of Interests
There is no conflict of interests about this paper among authors.
Shuang Liang, Yuanpeng Zhou, and Huijun Wang are cofirst authors, equally contributing to the paper.
This work was supported by National Basic Research Program of China (2010CB529500), National Science Fund of China (81270712, 81300506, and 81200449), National Academic-Specific Program of Health Care, the Public Health Care Program of Shanghai (12GWZX0301) from Shanghai Municipal Health Bureau, and the Training Program of Shanghai Academic Leading Talent by Shanghai Committee of Science and Technology and Shanghai Bureau of Human Resources. The authors would like to thank all of the subjects who contributed their time to this study. They would also like to thank Weili Yan and Naiqing Zhao for their assistance with data management.
- A. J. Levine, W. Lee, J. C. Figueiredo et al., “Variation in folate pathway genes and distal colorectal adenoma risk: a sigmoidoscopy-based case-control study,” Cancer Causes and Control, vol. 22, no. 4, pp. 541–552, 2011.
- I. J. M. van der Linden, L. A. Afman, S. G. Heil, and H. J. Blom, “Genetic variation in genes of folate metabolism and neural-tube defect risk,” Proceedings of the Nutrition Society, vol. 65, no. 2, pp. 204–215, 2006.
- O. Stanger, B. Fowler, K. Pietrzik et al., “Homocysteine, folate and vitamin B12 in neuropsychiatric diseases: review and treatment recommendations,” Expert Review of Neurotherapeutics, vol. 9, no. 9, pp. 1393–1412, 2009.
- I. A. L. M. van Rooij, M. C. Ocké, H. Straatman, G. A. Zielhuis, H. M. W. M. Merkus, and R. P. M. Steegers-Theunissen, “Periconceptional folate intake by supplement and food reduces the risk of nonsyndromic cleft lip with or without cleft palate,” Preventive Medicine, vol. 39, no. 4, pp. 689–694, 2004.
- E. Goldmuntz, S. Woyciechowski, D. Renstrom, P. J. Lupo, and L. E. Mitchell, “Variants of folate metabolism genes and the risk of conotruncal cardiac defects,” Circulation: Cardiovascular Genetics, vol. 1, no. 2, pp. 126–132, 2008.
- D. Taruscio, P. Carbone, O. Granata, F. Baldi, and A. Mantovani, “Folic acid and primary prevention of birth defects,” BioFactors, vol. 37, no. 4, pp. 280–284, 2011.
- D. L. Furness, N. Yasin, G. A. Dekker, S. D. Thompson, and C. T. Roberts, “Maternal red blood cell folate concentration at 10–12 weeks gestation and pregnancy outcome,” Journal of Maternal-Fetal and Neonatal Medicine, vol. 25, no. 8, pp. 1423–1427, 2012.
- L. Zhang, W. Liu, Q. Hao, L. Bao, and K. Wang, “Folate intake and methylenetetrahydrofolate reductase gene polymorphisms as predictive and prognostic biomarkers for ovarian cancer risk,” International Journal of Molecular Sciences, vol. 13, no. 4, pp. 4009–4020, 2012.
- J. Song, K. Sohn, A. Medline, C. Ash, S. Gallinger, and Y. Kim, “Chemopreventive effects of dietary folate on intestinal polyps in Apc+/- Msh2-/- Mice,” Cancer Research, vol. 60, no. 12, pp. 3191–3199, 2000.
- J. Song, A. Medline, J. B. Mason, S. Gallinger, and Y. I. Kim, “Effects of dietary folate on intestinal tumorigenesis in the Apc(Min) mouse,” Cancer Research, vol. 60, no. 19, pp. 5434–5440, 2000.
- A. M. Devlin, R. Clarke, J. Birks, J. G. Evans, and C. H. Halsted, “Interactions among polymorphisms in folate-metabolizing genes and serum total homocysteine concentrations in a healthy elderly population,” The American Journal of Clinical Nutrition, vol. 83, no. 3, pp. 708–713, 2006.
- M. P. S. Brown, W. N. Grundy, D. Lin et al., “Knowledge-based analysis of microarray gene expression data by using support vector machines,” Proceedings of the National Academy of Sciences of the United States of America, vol. 97, no. 1, pp. 262–267, 2000.
- T. Huang, K. L. Tucker, Y. Lee et al., “Methylenetetrahydrofolate reductase variants associated with hypertension and cardiovascular disease interact with dietary polyunsaturated fatty acids to modulate plasma homocysteine in Puerto Rican adults,” Journal of Nutrition, vol. 141, no. 4, pp. 654–659, 2011.
- K. N. Saraswathy, M. Asghar, R. Samtani et al., “Spectrum of MTHFR gene SNPs C677T and A1298C: a study among 23 population groups of India,” Molecular Biology Reports, vol. 39, no. 4, pp. 5025–5031, 2012.
- M. Y. Tsai, M. Bignell, F. Yang, B. G. Welge, K. J. Graham, and N. Q. Hanson, “Polygenic influence on plasma homocysteine: association of two prevalent mutations, the 844ins68 of cystathionine β-synthase and A2756G of methionine synthase, with lowered plasma homocysteine levels,” Atherosclerosis, vol. 149, no. 1, pp. 131–137, 2000.
- V. Komlósi, E. Hitre, E. Pap et al., “SHMT1 1420 and MTHFR 677 variants are associated with rectal but not colon cancer,” BMC Cancer, vol. 10, article 525, 2010.
- L. Pena-Castillo, M. Tasan, C. L. Myers, et al., “A critical assessment of Mus musculus gene function prediction using integrated genomic evidence,” Genome Biology, vol. 9, supplement 1, article S2, 2008.
- S. Ramaswamy, P. Tamayo, R. Rifkin et al., “Multiclass cancer diagnosis using tumor gene expression signatures,” Proceedings of the National Academy of Sciences of the United States of America, vol. 98, no. 26, pp. 15149–15154, 2001.
- A. Stanisławska-Sachadyn, L. E. Mitchell, J. V. Woodside et al., “The reduced folate carrier (SLC19A1) c.80G>A polymorphism is associated with red cell folate concentrations among women,” Annals of Human Genetics, vol. 73, no. 5, pp. 484–491, 2009.
- J. Chen, M. J. Stampfer, J. Ma et al., “Influence of a methionine synthase (D919G) polymorphism on plasma homocysteine and folate levels and relation to risk of myocardial infarction,” Atherosclerosis, vol. 154, no. 3, pp. 667–672, 2001.
- A. Parle-McDermott, F. Pangilinan, K. K. O'Brien et al., “A common variant in MTHFD1L is associated with neural tube defects and mRNA splicing efficiency,” Human Mutation, vol. 30, no. 12, pp. 1650–1656, 2009.
- A. Chango, N. Emery-Fillon, G. P. de Courcy et al., “A polymorphism (80G->A) in the reduced folate carrier gene and its associations with folate status and homocysteinemia,” Molecular Genetics and Metabolism, vol. 70, no. 4, pp. 310–315, 2000.
- S. van Looy, T. Verplancke, D. Benoit et al., “A novel approach for prediction of tacrolimus blood concentration in liver transplantation patients in the intensive care unit through support vector regression,” Critical Care, vol. 11, no. 4, article R83, 2007.
- R. H. Finnell, G. M. Shaw, E. J. Lammer, and K. A. Volcik, “Does prenatal screening for 5,10-methylenetetrahydrofolate reductase (MTHFR) mutations in high-risk neural tube defect pregnancies make sense?” Genetic Testing, vol. 6, no. 1, pp. 47–52, 2002.
- N. M. J. van der Put, F. Gabreëls, E. M. B. Stevens et al., “A second common mutation in the methylenetetrahydrofolate reductase gene: an additional risk factor for neural-tube defects?” The American Journal of Human Genetics, vol. 62, no. 5, pp. 1044–1051, 1998.
- A. Hazra, P. Kraft, R. Lazarus et al., “Genome-wide significant predictors of metabolites in the one-carbon metabolism pathway,” Human Molecular Genetics, vol. 18, no. 23, pp. 4677–4687, 2009.
- G. M. Shaw, W. Lu, H. Zhu et al., “118 SNPs of folate-related genes and risks of spina bifida and conotruncal heart defects,” BMC Medical Genetics, vol. 10, article 49, 2009.
- B. Hooshmand, A. Solomon, I. Kåreholt et al., “Homocysteine and holotranscobalamin and the risk of Alzheimer disease: a longitudinal study,” Neurology, vol. 75, no. 16, pp. 1408–1414, 2010.
- B. M. Riedel, A. M. Molloy, K. Meyer et al., “Transcobalamin polymorphism 67A->G, but not 776C->G, affects serum holotranscobalamin in a cohort of healthy middle-aged men and women,” Journal of Nutrition, vol. 141, no. 10, pp. 1784–1790, 2011.
- J. Listgarten, S. Damaraju, B. Poulin et al., “Predictive models for breast cancer susceptibility from multiple single nucleotide polymorphisms,” Clinical Cancer Research, vol. 10, no. 8, pp. 2725–2737, 2004.