- Open Access
Associations of the MTHFR rs1801133 polymorphism with coronary artery disease and lipid levels: a systematic review and updated meta-analysis
Lipids in Health and Disease volume 17, Article number: 191 (2018)
The associations of the 5,10-methylenetetrahydrofolate reductase gene (MTHFR) rs1801133 polymorphism with coronary artery disease (CAD) and plasma lipid levels have been widely investigated, but the results were inconsistent and inconclusive. This meta-analysis aimed to clarify the relationships of the rs1801133 polymorphism with CAD and plasma lipid levels.
By searching in PubMed, Google Scholar, Web of Science, Cochrane Library, Wanfang, VIP and CNKI databases, 123 studies (87,020 subjects) and 65 studies (85,554 subjects) were identified for the CAD association analysis and the lipid association analysis, respectively. Odds ratio (OR) and standardized mean difference (SMD) were used to determine the effects of the rs1801133 polymorphism on CAD risk and lipid levels, respectively.
The variant T allele of the rs1801133 polymorphism was associated with increased risk of CAD under allelic model [OR = 1.11, 95% confidence interval (CI) = 1.06–1.17, P < 0.01], additive model (OR = 1.25, 95% CI = 1.14–1.37, P < 0.01), dominant model (OR = 1.11, 95% CI = 1.04–1.17, P < 0.01), and recessive model (OR = 1.22, 95% CI = 1.12–1.32, P < 0.01). The T carriers had higher levels of total cholesterol (TC) (SMD = 0.04, 95% CI = 0.01–0.07, P = 0.02) and low-density lipoprotein cholesterol (LDL-C) (SMD = 0.07, 95% CI = 0.01–0.12, P = 0.01) than the non-carriers.
The meta-analysis suggested that the T allele of the rs1801133 polymorphism is a risk factor for CAD, which is possibly and partly mediated by abnormal lipid levels.
Coronary artery disease (CAD) is currently the leading cause of death in developed countries, and in some developing countries like China . CAD is a multifactorial disease and a number of risk factors have been identified in the past few decades. Genetic polymorphism and dyslipidemia are the two most important risk factors for CAD [2, 3]. Genetic polymorphism is a term used to describe multiple forms of a single gene. Dyslipidemia is a state of abnormal amounts of lipids (e.g. triglycerides, cholesterol and/or phospholipids) in the blood, and is characterized by increased levels of triglycerides (TG), total cholesterol (TC) and low-density lipoprotein cholesterol (LDL-C), and/or decreased level of high-density lipoprotein cholesterol (HDL-C). Intensive efforts have been made in the scientific community to investigate the associations of the genetic polymorphisms in some specific genes with CAD risk and plasma lipid levels, but the results were inconsistent and inconclusive. It is difficult to identify the CAD- or dyslipidemia-related genetic polymorphisms successfully due to various reasons such as small sample sizes and ethnic differences.
5,10-methylenetetrahydrofolate reductase (MTHFR) is a rate-limiting enzyme in the one-carbon metabolism pathway and plays a key role in one-carbon metabolism by irreversibly catalyzing the conversion of 5,10-Methylenetetrahydrofolate (5,10-MTHF) to 5-Methyltetrahydrofolate (5-MTHF). 5-MTHF is a direct one-carbon donor (methyl group) for many substrates such as DNA , RNA  and proteins . More importantly, 5-MTHF is the only one-carbon donor for the remethylation of homocysteine which is produced in methionine cycle. In methionine cycle, methionine first reacts with adenosine triphosphate to form S-adenosylmethionine (SAM) under the catalysis of methionine adenosyltransferase. The methyl group in SAM is activated, and SAM is called activated methionine. The activated methyl group of SAM can be transferred to target substrates such as DNA under the catalysis of methyltransferase, and SAM itself is converted into S-adenosine homocysteine after demethylation. Homocysteine is produced after the removal of adenosine from S-adenosine homocysteine under the catalysis of S-adenosylhomocysteine hydrolase. In the last step, homocysteine accepts the methyl group from 5-MTHF and methionine is formed again.
Homocysteine a potential risk factor for CAD and MTHFR has been reported to be associated with CAD risk and abnormal lipid levels [7,8,9]. Mikael et al.  reported that MTHFR(+/−) mice had significantly higher levels of plasma TG and more lipid deposition in aortic sinus compared with MTHFR (+/+) mice. In another study, Christensen et al.  found that high folic acid consumption led to pseudo-MTHFR deficiency in mice, and this deficiency resulted in altered lipid metabolism and liver injury.
The rs1801133 polymorphism (also known as the 677C > T polymorphism) is located in exon 4 of the MTHFR gene and formed by a transition from cytosine (C) to thymine (T). The 222nd genetic code of the MTHFR gene is changed accordingly from GCC to GTC, resulting in the replacement of alanine (Ala) by valine (Val) in the MTHFR polypeptide. A large number of studies have investigated the associations of the rs1801133 polymorphism with CAD and lipid levels. In some of these studies, the T allele of the rs1801133 polymorphism was reported to be associated with an increased risk of CAD [10,11,12] and elevated levels of TG [13, 14], TC [13,14,15] and LDL-C [13,14,15,16], and reduced levels of HDL-C [13, 17]. However, the results obtained from other studies did not support these findings [18,19,20,21,22]. Hence, a meta-analysis is required to clarify the relationships of the rs1801133 polymorphism with CAD and lipid levels.
Although two meta-analyses [23, 24] have addressed the issue of the association between the rs1801133 polymorphism and CAD in 2002 and 2005, respectively, their sample sizes were relatively small, and blood lipid variables were not considered in the analyses. In this study, a systematic review and updated meta-analysis was performed based on previous publications to investigate the associations of the rs1801133 polymorphism with CAD and lipid levels. The results of this meta-analysis can provide an opportunity to unveil the interrelationships among the rs1801133 polymorphism, dyslipidemia and susceptibility to CAD.
Characteristics of the included studies
Initial search of the databases yielded 5197 articles. Four thousand nine hundred and eighty-one studies were excluded according to the titles and abstracts. Then full-text articles were retrieved and assessed on the basis of inclusion criteria. Thirty-seven studies were ineligible for the following reasons: 28 studies presented data for other polymorphisms; 5 studies had subjects overlapping with other publications; 3 studies were based on pedigree analysis; 1 study presented invalid data. In the end, 179 studies were selected for this meta-analysis (Fig. 1). One hundred twenty-three studies (87,020 subjects) of them were included in the CAD association analysis, and 65 studies (85,554 subjects) were included in the lipid association analysis. The references for the studies included in the present meta-analysis are listed in Additional file 1.
The characteristics of the studies included in the CAD association analysis are summarized in Additional file 2: Table S1. Seventy-two studies, 35 studies, 5 studies and 11 studies involved in Caucasians, Asians, Africans and the subjects of other ethnic origins, respectively. The characteristics of the studies included in the lipid association analysis are summarized in Additional file 2: Table S2. The plasma lipid levels according to the genotypes of the rs1801133 polymorphism are presented in Additional file 2: Table S3. Twenty-four studies, 25 studies, 3 studies and 13 studies involved in Caucasians, Asians, Africans and the subjects of other ethnic origins, respectively. Eleven studies, 5 studies, 6 studies and 33 studies involved in CAD, diabetes, hypertension and healthy subjects, respectively. Fifty-three studies, 61 studies, 49 studies and 58 studies presented the data for TG, TC, LDL-C and HDL-C, respectively.
Association of the r1801133 polymorphism with CAD
The variant T allele of the rs1801133 polymorphism was associated with increased risk of CAD under allelic model [odds ratio (OR) = 1.11, 95% confidence interval (CI) = 1.06–1.17, P < 0.01], additive model (OR = 1.25, 95% CI = 1.14–1.37, P < 0.01), dominant model (OR = 1.11, 95% CI = 1.04–1.17, P < 0.01) and recessive model (OR = 1.22, 95% CI = 1.12–1.32, P < 0.01) (Table 1). When the analyses were limited to the studies in Hardy-Weinberg equilibrium (HWE), the association between the rs1801133 polymorphism and CAD under allelic model (OR = 1.11, 95% CI = 1.06–1.16, P < 0.01), additive model (OR = 1.23, 95% CI = 1.12–1.35, P < 0.01), dominant model (OR = 1.11, 95% CI = 1.05–1.17, P < 0.01) and recessive model (OR = 1.20, 95% CI = 1.11–1.30, P < 0.01) was also significant (Table 1). Then the subgroup analyses stratified by the ethnicity of the subjects were performed, and the results showed that the rs1801133 polymorphism was strongly associated with CAD in Asians, but not in Africans or the subjects of other ethnic origins. In Caucasians, the association between the rs1801133 polymorphism and CAD was ambiguous, i.e. the rs1801133 polymorphism was marginally significantly associated with CAD under additive model (OR = 1.08, 95% CI = 1.00–1.17, P = 0.05) and recessive model (OR = 1.09, 95% CI = 1.01–1.17, P = 0.02), but not under allelic model and dominant model (Table 1).
Associations of the r1801133 polymorphism with plasma lipid levels
The outcomes of the analyses on all comparisons showed that the T allele carriers had higher levels of TC [standardized mean difference (SMD) = 0.04, 95% CI = 0.01–0.07, P = 0.02] and LDL-C (SMD = 0.07, 95% CI = 0.01–0.12, P = 0.01) than the non-carriers (Table 2, Figs. 2 and 3), and that there were no associations detected between the rs1801133 polymorphism and plasma levels of TG (SMD = 0.03, 95% CI = − 0.01-0.06, P = 0.11) and HDL-C (SMD = − 0.02, 95% CI = − 0.05-0.02, P = 0.30) (Table 2, Figs. 4 and 5). When the analyses were limited to the studies in HWE, the T allele carriers had higher levels of LDL-C (SMD = 0.04, 95% CI = 0.01–0.08, P = 0.03) and lower levels of HDL-C (SMD = − 0.02, 95% CI = − 0.03--0.00, P = 0.04) than the non-carriers (Table 2).
Then the subgroup analyses stratified by the characteristics of the subjects were performed. The significant associations of the rs1801133 polymorphism with higher levels of TC (SMD = 0.11, 95% CI = 0.02–0.19, P = 0.02) and LDL-C (SMD = 0.09, 95% CI = 0.01–0.16, P = 0.02) were detected in females, but not in males. The significant associations of the rs1801133 polymorphism with higher levels of TC (SMD = 0.10, 95% CI = 0.04–0.15, P < 0.01) and LDL-C (SMD = 0.10, 95% CI = 0.05–0.15, P < 0.01) were detected in Asians, but not in the subjects of other ethnic origins. However, the rs1801133 polymorphism was significantly associated with lower level of HDL-C (SMD = − 0.10, 95% CI = − 0.18--0.02, P = 0.01) only in the subjects of other ethnic origins. When health status was taken into account, the significant associations of the rs1801133 polymorphism with higher levels of TC (SMD = 0.07, 95% CI = 0.02–0.12, P < 0.01) and LDL-C (SMD = 0.07, 95% CI = 0.01–0.13, P = 0.02) were detected in healthy or control subjects, but not in the patients with CAD, type 2 diabetes mellitus (T2DM) and hypertension (Table 2).
In the CAD association analysis, there was significant heterogeneity in allelic model (I2 = 68.2%, Pheterogeneity < 0.01), additive model (I2 = 55.6%, Pheterogeneity < 0.01), dominant model (I2 = 64.0%, Pheterogeneity < 0.01) and recessive model (I2 = 51.5%, Pheterogeneity < 0.01). Twenty-seven studies, 22 studies, 19 studies and 17 studies were identified as the main contributors to the heterogeneity for allelic model, additive model, dominant model and recessive model, respectively, by using Galbraith plots. The heterogeneity was effectively removed or decreased after exclusion of these outlier studies, but OR values and 95% CIs did not change substantially (allelic model: OR = 1.04, 95% CI = 1.01–1.07, Pheterogeneity = 0.12, POR = 0.01; additive model: OR = 1.12, 95% CI = 1.05–1.19, Pheterogeneity = 0.09, POR < 0.01; dominant model: OR = 1.04, 95% CI = 1.00–1.07, Pheterogeneity = 0.14, POR = 0.05; recessive model: OR = 1.12, 95% CI = 1.06–1.18, Pheterogeneity = 0.13, POR < 0.01) (Table 3).
In the lipid association analysis, there was significant heterogeneity in the total comparisons for TG (I2 = 36.1%, Pheterogeneity < 0.01), TC (I2 = 54.1%, Pheterogeneity < 0.01), LDL-C (I2 = 69.6%, Pheterogeneity < 0.01) and HDL-C (I2 = 60.1%, Pheterogeneity < 0.01). Three comparisons, 9 comparisons, 3 comparisons and 6 comparisons were identified as the main contributors to the heterogeneity for TG, TC, LDL-C and HDL-C, respectively, by using Galbraith plots. SMD values and 95% CIs of TG and LDL-C did not change substantially after excluding the outlier comparisons. However, SMD values and 95% CIs of TC (SMD = 0.01, 95% CI = − 0.00-0.03, Pheterogeneity = 0.38, PSMD = 0.15) and HDL-C (SMD = − 0.02, 95% CI = − 0.03--0.00, Pheterogeneity = 0.48, PSMD = 0.05) changed significantly after excluding these outlier comparisons (Table 4).
Publication bias test
Begg’s test and Egger’s test were used to evaluate the publication bias of the included studies, and the results showed that there might be a publication bias in the analysis between the rs1801133 polymorphism and CAD (P < 0.05 for all genetic models). To clarify this problem, a trim-and-fill method was used to adjust the results, and no trimming was performed and the results were unchanged. It indicates that there is no publication bias in the literature. The significant P values of Begg’s test and Egger’s test were originated from other factors, e.g. mixed ethnicity in some studies. No publication bias was detected in the lipid association analysis.
In the present meta-analysis, the variant T allele of the rs1801133 polymorphism was associated with increased risk of CAD, and elevated levels of TC and LDL-C in the total population. It indicates that the T allele of the rs1801133 polymorphism is a risk factor for CAD, which is at least partly mediated by abnormal lipid levels.
A large number of studies have investigated the association between the rs1801133 polymorphism and CAD risk, as well as the underlying mechanisms, but most of them just focused on homosysteine. It was widely reported that the rs1801133 polymorphism influences the plasma levels of homosysteine in various populations such as Americans , Africans [17, 26, 27], Asians [18, 28, 29], Turkish [30, 31] and Brazilians . In methionine cycle, homocysteine is formed after the removal of adenosine from S-adenosine homocysteine. Under normal condition, homocysteine is remethylated to methionine by accepting a methyl group from 5-MTHF. 5-MTHF is formed from the reduction of 5,10-MTHF under the catalysis of MTHFR. In the T allele carriers of the rs1801133 polymorphism, the function of MTHFR may be affected since the normal alanine residue is replaced by a valine residue in the polypeptide, which in turn affects the production of 5-MTHF and the remethylation of homcysteine, resulting in elevated plasma homosysteine levels. Studies have shown that homocysteine is a risk factor for CAD, and oxidative stress , vascular inflammation  and endothelial injury  are involved in the underlying mechanisms in which homocysteine causes CAD. All these events are likely to trigger the development of atherosclerosis and arterial thrombosis.
However, the role of homocysteine in the pathogenesis of CAD is controversial. Several studies [36,37,38] demonstrated no association between homocysteine and CAD. It indicates that some other risk factors are involved in the association between the rs1801133 polymorphism and CAD. In this meta-analysis, the results showed that the variant T allele carriers of the rs1801133 polymorphism have higher levels of TC and LDL-C than the non-carriers, which indicates that the abnormal lipid levels caused by the T allele of the rs1801133 polymorphism might be one of the important reasons in the development of CAD since dyslipidemia is closely associated with the progression of coronary atherosclerosis, and it accounts for around 50% of the population-attributable risk for CAD . According to the 2013 ACC/AHA blood cholesterol guidelines  and the Adult Treatment Panel III (ATP III) Guidelines  of the United States, LDL-C was considered as a major cause of CAD and used as the primary target for therapy, and other lipid parameters were used as the secondary or supplementary targets.
The mechanisms in which the rs1801133 polymorphism is associated with plasma lipid levels have not been clarified yet. Several possible reasons can be proposed to explain the association between the rs1801133 polymorphism and plasma lipid levels. Firstly, the rs1801133 polymorphism may indirectly affect plasma lipid levels through the mediation of homocysteine [42,43,44]. Baszczuk et al.  reported that a daily administration of 15 mg of folic acid led to a considerable decrease in homocysteine levels, and a substantive increase in HDL-C levels in the patients with primary hypertension. In yeast cells, homocysteine supplementation increased cellular fatty acid and TG contents, induced a shift in fatty acid composition, and decreased the condensing enzymes involved in very long-chain fatty acid synthesis . Secondly, the rs1801133 polymorphism may modulate the lipid metabolism by affecting the methylation state of DNA or proteins. 5-MTHF is not only the methyl donor for homocysteine, but for many other target molecules, including DNA and proteins . Conceivably, the functions of the genes or proteins involved in lipid metabolism will be affected if their methylation state changes.
In most of the studies included in the lipid association analysis, a dominant model was used, i.e. TT + CT vs. CC. Therefore, a dominant model was also adopted in this meta-analysis. In subgroup analyses, we found that the differences in TC and LDL-C levels between the genotypes were mainly from Asian populations, whose SMD values were larger than those calculated in Caucasians, Africans and the subjects of other ethnicities (Table 2). The associations of the rs1801133 polymorphism with TC and LDL-C in Asians were consistently larger, which shows that there is a stronger association between the rs1801133 polymorphism and CAD in Asians as compared with other ethnicities (Table 1). Studies will be needed to elucidate the mechanisms that the rs1801133 polymorphism has different effects on blood lipid levels and CAD risk in different ethnicities.
In the lipid association analysis, subgroup analyses by gender and health status was performed since there might be important factors affecting the associations between the rs1801133 polymorphism and lipid levels. For example, the present meta-analysis indicates that gender might modulate the associations of the rs1801133 polymorphism with TC and LDL-C levels since there were significant associations existing only in females but not in males (Table 2). Health status might also modulate the associations of the rs1801133 polymorphism with TC and LDL-C levels. The significant associations of the rs1801133 polymorphism with TC and LDL-C levels only existed in healthy subjects, but not in the patients with CAD, T2DM and hypertension. The reason might be that the patients with CAD, T2DM and hypertension had serious metabolic disorders, which masked the effects of the rs1801133 polymorphism on TC and LDL-C levels. As compared with the diseases such as CAD, T2DM and hypertension, genetic polymorphisms usually have less impact on plasma lipid levels. In line with the results from the present study, several studies also reported that the rs1801133 polymorphism is associated TC and LDL-C levels in healthy subjects [13, 15, 46], but not in the patients with CAD [11, 20], T2DM [47, 48] and hypertension [49, 50]. Of the 179 studies included, 141 studies used polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method; 27 studies used real-time PCR method; 7 studies used DNA sequencing method; 1 study used gene chip method; and 3 studies did not report the genotyping method(s). Subgroup analyses stratified by the genotyping methods were conducted, and the results showed that there were differences in OR or SMD values among the studies with different genotyping methods (data not shown). In most cases, the results from the studies by PCR-RFLP method were in line with the results from all studies. The reason might be that most of the studies used PCR-RFLP method. The results from the studies by real-time PCR method or DNA sequencing method could have been affected by the small number of studies and small sample sizes.
Significant heterogeneity was detected in the total and subgroup analyses between the rs1801133 polymorphism and CAD, and the outlier studies were identified by using the Galbraith plots. No significant changes in OR values and 95% CIs were found after excluding the outlier studies (Table 3), which indicates that the association between the rs1801133 polymorphism and CAD is very strong. Significant heterogeneity was also detected in the total and some of the subgroup analyses between the rs1801133 polymorphism and plasma lipid levels. The outlier studies were identified and excluded, but SMD values and 95% CIs of LDL-C were not significantly changed in the total population or in Asians, which indicates that there is a strong association between the rs1801133 polymorphism and LDL-C levels, especially in Asians.
The associations of the rs1801133 polymorphism with CAD and plasma lipid levels are not likely to be type I errors (false-positive results). Firstly, the results from this meta-analysis are based on four different models for CAD association analysis, and on random effects model for lipid association analysis if the heterogeneity among the studies is significant (I2 > 50%). As compared with fixed effects model, the random effects model is a more conservative method and less likely to produce false-positive results. Secondly, 87,020 subjects and 85,554 subjects were included in the analysis for the CAD association analysis and the lipid association analysis, respectively. Among the subjects, 44.0% (CAD association analysis) and 52.3% (lipid association analysis) of them are the carriers of the variant T allele. Since the incidence of the T allele carriers is very high, type I errors could have been prevented in both the CAD and lipid association analyses. Thirdly, this meta-analysis only included the studies published in English and Chinese as it was very difficult to get the full papers published in various languages.
The current meta-analysis demonstrates that the rs1801133 polymorphism is associated with increased risk of CAD and elevated levels of TC and LDL-C. Further studies will be needed to elucidate the underlying mechanisms in which the rs1801133 polymorphism affects plasma lipid levels.
The articles published before September 2017 on the associations of the rs1801133 polymorphism with CAD and/or plasma lipid levels were identified. The languages of the articles were limited to English and Chinese. A comprehensive search was conducted and nine electronic databases were searched to identify all relevant articles. The databases are as follows: PubMed, Embase, Baidu Scholar, Google Scholar, Web of Science, Cochrane Library, Wanfang, CBM and CNKI. The following keywords were used: (“5,10-Methylenetetrahydrofolate reductase” or “Methylenetetrahydrofolate reductase” or “MTHFR”), (“polymorphism” or “mutation” or “variant” or “C677T” or “rs1801133” or “Ala222Val”), (“coronary artery disease” or “coronary heart disease” or “heart disease” or “coronary disease” or “cardiovascular disease” or “angina pectoris” or “acute coronary syndrome” or “myocardial infarction” or “CAD” or “CHD” or “HD” or “CD” or “AP” or “ACS” or “MI”), (“plasma lipid” or “blood lipid” or “serum lipid”).
Inclusion and exclusion criteria
The inclusion criteria for the association analysis between the rs1801133 polymorphism and CAD are as follows: 1) studies using a population-based case-control design; 2) CAD cases were angiographically defined; 3) number or frequency of cases according to the r1801133 genotypes was available. The inclusion criteria for the association analysis between the rs1801133 polymorphism and lipid levels are as follows: 1) studies in which mean lipids and standard deviations (SD) or standard errors (SE) by the rs1801133 genotypes were available; 2) studies which reported at least one of the four lipid variables, i.e. TG, TC, LDL-C and HDL-C; 3) baseline data were used for interventional studies. All references cited by the included articles were reviewed to check the published work which was not indexed by PubMed, Embase, Baidu Scholar, Google Scholar, Web of Science, Cochrane Library, Wanfang, CBM and CNKI. Reports with incomplete data, studies based on pedigree data, case reports, review articles, abstracts and animal studies were excluded from the meta-analysis.
Data were extracted from each study by using a structured data collection form and by two investigators independently according to the pre-specified selection criteria. Decisions were compared and disagreements about study selection were resolved by consensus or by involving a third investigator. For the overlapping articles, only the publications that presented the most detailed information were included. In this meta-analysis, the data extracted from each of the included studies are as follows: first author, year of publication, age, ethnicity, gender, health status, type of study, genotyping method, lipid assay method, sample size, and mean with SD or SE according to the r1801133 genotypes. If data in a study were unconvincing, we attempted to contact the corresponding or first author by e-mail and telephone.
Statistical analysis was performed by using STATA version 12.0 (Stata Corporation LP, College Station, TX, USA). All the tests were two-sided and a P-value of less than 0.05 for any test or model was considered to be statistically significant. OR with 95% CI was used to evaluate the strength of the association between the rs1801133 polymorphism and CAD. The pooled OR was performed for allelic model (T vs C), additive model (TT vs CC), dominant model (TT + CT vs CC) and recessive model (TT vs CT + CC). SMD with 95% CI was used to assess the strength of the associations between the rs1801133 polymorphism and plasma lipid levels. A fixed-effect model (Mantel-Haenszel method) was used to evaluate the results if heterogeneity among the included studies was not significant (I2 < 50%). Otherwise, the random-effect model (DerSimonian-Laird method) was used . Heterogeneity was investigated by Cochran’s χ2-based Q-statistic, and Galbraith plots were used to detect the potential sources of heterogeneity. OR and SMD values were recalculated after excluding the outlier studies. Subgroup analyses were performed according to ethnicity for CAD association analysis, and according to ethnicity, gender and health status for lipid association analysis. Ethnic subgroups were defined as Caucasian, Asian, African, and the subjects of other ethnic origin. Health status was defined as CAD, T2DM and hypertension. HWE was assessed by Fisher’s exact test. OR and SMD values were recalculated after excluding the studies which were not in HWE. Publication bias was tested by Begg’s funnel plots and Egger’s test .
- 95% CI:
95% confidence interval
coronary artery disease
high-density lipoprotein cholesterol
low-density lipoprotein cholesterol
standardized mean difference
type 2 diabetes mellitus
Mozaffarian D, Benjamin EJ, Go AS, Arnett DK, Blaha MJ, Cushman M, Das SR, de Ferranti S, Després JP, Fullerton HJ, Howard VJ, Huffman MD, Isasi CR, Jiménez MC, Judd SE, Kissela BM, Lichtman JH, Lisabeth LD, Liu S, Mackey RH, Magid DJ, DK MG, Mohler ER 3rd, Moy CS, Muntner P, Mussolino ME, Nasir K, Neumar RW, Nichol G, Palaniappan L, Pandey DK, Reeves MJ, Rodriguez CJ, Rosamond W, Sorlie PD, Stein J, Towfighi A, Turan TN, Virani SS, Woo D, Yeh RW, Turner MB. Heart Disease and Stroke Statistics-2016 Update: A Report From the American Heart Association. Circulation. 2016;133:e38–360.
Erdmann J, Kessler T, Munoz Venegas L, Schunkert H. A decade of genome-wide association studies for coronary artery disease: the challenges ahead. Cardiovasc Res. 2018; https://doi.org/10.1093/cvr/cvy084. [Epub ahead of print]
O'Donnell FL, Stahlman S, Oetting AA. Incidence rates of diagnoses of cardiovascular diseases and associated risk factors, active component, U.S. Armed Forces, 2007–2016. MSMR. 2018;25:12–8.
Wang L, Shangguan S, Chang S, Yu X, Wang Z, Lu X, Wu L, Zhang T. Determining the association between methylenetetrahydrofolate reductase (MTHFR) genepolymorphisms and genomic DNA methylation level: a meta-analysis. Birth Defects Res A Clin Mol Teratol. 2016;106:667–74.
Xhemalce B. From histones to RNA: role of methylation in cancer. Brief Funct Genomics. 2013;12:244–53.
Lanouette S, Mongeon V, Figeys D, Couture JF. The functional diversity of protein lysine methylation. Mol Syst Biol. 2014;10:724.
Kukrele P, Sharma RS. High homocysteine level- an important risk factor for coronary artery disease in vegetarians. J Assoc Physicians India. 2016;64:75.
Mikael LG, Wang XL, Wu Q, Jiang H, Maclean KN, Rozen R. Hyperhomocysteinemia is associated with hypertriglyceridemia in mice with methylenetetrahydrofolate reductase deficiency. Mol Genet Metab. 2009;98:187–94.
Christensen KE, Mikael LG, Leung KY, Lévesque N, Deng L, Wu Q, Malysheva OV, Best A, Caudill MA, Greene ND, Rozen R. High folic acid consumption leads to pseudo-MTHFR deficiency, altered lipid metabolism, and liver injury in mice. Am J Clin Nutr. 2015;101:646–58.
Husemoen LL, Skaaby T, Jørgensen T, Thuesen BH, Fenger M, Grarup N, Sandholt CH, Hansen T, Pedersen O, Linneberg A. MTHFR C677T genotype and cardiovascular risk in a general population without mandatory folic acid fortification. Eur J Nutr. 2014;53:1549–59.
Ramkaran P, Phulukdaree A, Khan S, Moodley D, Chuturgoon AA. Methylenetetrahydrofolate reductase C677T polymorphism is associated with increased risk of coronary artery disease in young south African Indians. Gene. 2015;571:28–32.
Chen W, Hua K, Gu H, Zhang J, Wang L. Methylenetetrahydrofolate reductase C667T polymorphism is associated with increased risk of coronary artery disease in a Chinese population. Scand J Immunol. 2014;80:346–53.
Zhi X, Yang B, Fan S, Wang Y, Wei J, Zheng Q, Sun G. 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.
Yin RX, Wu DF, Miao L, Htet Aung LH, Cao XL, Yan TT, Long XJ, Liu WY, Zhang L, Li M. Interactions of several single nucleotide polymorphisms and high body mass index on serum lipid traits. Biofactors. 2013;39:315–25.
Chen NY, Liu CW, Du LL, Xiao LP, Ge L, Wang YY, Wei Z, Wu HY, Luo CY, Liang L, Peng JH, Luo XQ, Yin RX, Nguyen CP, Pan SL. Enrichment of MTHFR 677 T in a Chinese long-lived cohort and its association with lipid modulation. Lipids Health Dis. 2014;13:104.
Jiang S, Zhao R, Pan M, Venners SA, Zhong G, Hsu YH. Associations of MTHFR and MTRR polymorphisms with serum lipid levels in Chinese hypertensive patients. Clin Appl Thromb Hemost. 2014;20:400–10.
Abd-Elmawla MA, Rizk SM, Youssry I, Shaheen AA. Impact of genetic polymorphism of methylenetetrahydrofolate reductase C677T on development of Hyperhomocysteinemia and related oxidative changes in Egyptian β-thalassemia major patients. PLoS One. 2016;11:e0155070.
Saffari B, Senemar S, Karimi M, Bahari M, Jooyan N, Yavarian M. An MTHFR variant, plasma homocysteine levels and late-onset coronary artery disease in subjects from southern Iran. Pak J Biol Sci. 2013;16:788–95.
Senemar S, Saffari B, Sharifkazemi MB, Bahari M, Jooyan N, Dehaghani ED, Yavarian M. 5,10-methylene tetrahydrofolate reductase C677T gene polymorphism, homocysteine concentration and the extent of premature coronary artery disease in southern Iran. EXCLI J. 2013;12:437–48.
Kucukhuseyin O, Kurnaz O, Akadam-Teker AB, Isbir T, Bugra Z, Ozturk O, Yilmaz-Aydogan H. 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.
Frederiksen J, Juul K, Grande P, Jensen GB, Schroeder TV, Tybjaerg-Hansen A, Nordestgaard BG. Methylenetetrahydrofolate reductase polymorphism (C677T), hyperhomocysteinemia, and risk of ischemic cardiovascular disease and venous thromboembolism: prospective and case-control studies from the Copenhagen City heart study. Blood. 2004;104:3046–51.
Vasisht S, Gulati R, Narang R, Srivastava N, Srivastava LM, Manchanda SC, Agarwal DP. Polymorphism (C677T) in the 5,10-methylenetetrahydrofolate reductase (MTHFR) gene: a preliminary study on north Indian men. Indian J Clin Biochem. 2002;17:99–107.
Klerk M, Verhoef P, Clarke R, Blom HJ, Kok FJ, Schouten EG, MTHFR Studies Collaboration Group. MTHFR 677C-->T polymorphism and risk of coronary heart disease: a meta-analysis. JAMA. 2002;288:2023–31.
Lewis SJ, Ebrahim S, Davey SG. Meta-analysis of MTHFR 677C->T polymorphism and coronary heart disease: does totality of evidence support causal role for homocysteine and preventive potential of folate? BMJ. 2005;331:1053.
Brilakis ES, Berger PB, Ballman KV, Rozen R. Methylenetetrahydrofolate reductase (MTHFR) 677C>T and methionine synthase reductase (MTRR) 66A>G polymorphisms: association with serum homocysteine and angiographic coronary artery disease in the era of flour products fortified with folic acid. Atherosclerosis. 2003;168:315–22.
Bennouar N, Allami A, Azeddoug H, Bendris A, Laraqui A, El Jaffali A, El Kadiri N, Benzidia R, Benomar A, Fellat S, Benomar M. Thermolabile methylenetetrahydrofolate reductase C677T polymorphism and homocysteine are risk factors for coronary artery disease in Moroccan population. J Biomed Biotechnol. 2007;2007:80687.
Ghazouani L, Abboud N, Mtiraoui N, Zammiti W, Addad F, Amin H, Almawi WY, Mahjoub T. Homocysteine and methylenetetrahydrofolate reductase C677T and A1298C polymorphisms in Tunisian patients with severe coronary artery disease. J Thromb Thrombolysis. 2009;27:191–7.
Yadav S, Hasan N, Marjot T, Khan MS, Prasad K, Bentley P, Sharma P. Detailed analysis of gene polymorphisms associated with ischemic stroke in south Asians. PLoS One. 2013;8:e57305.
Biswas A, Ranjan R, Meena A, Akhter MS, Yadav BK, Munisamy M, Subbiah V, Behari M, Saxena R. Homocystine levels, polymorphisms and the risk of ischemic stroke in young Asian Indians. J Stroke Cerebrovasc Dis. 2009;18:103–10.
Ilhan N, Kucuksu M, Kaman D, Ilhan N, Ozbay Y. The 677 C/T MTHFR polymorphism is associated with essential hypertension, coronary artery disease, and higher homocysteine levels. Arch Med Res. 2008;39:125–30.
Ozbek Z, Kucukali CI, Ozkok E, Orhan N, Aydin M, Kilic G, Sazci A, Kara I. Effect of the methylenetetrahydrofolate reductase gene polymorphisms on homocysteine, folate and vitamin B12 in patients with bipolar disorder and relatives. Prog Neuro-Psychopharmacol Biol Psychiatry. 2008;32:1331–7.
Pereira AC, Miyakawa AA, Lopes NH, Soares PR, de Oliveira SA, Cesar LA, Ramires JF, Hueb W, Krieger JE. Dynamic regulation of MTHFR mRNA expression and C677T genotype modulate mortality in coronary artery disease patients after revascularization. Thromb Res. 2007;121:25–32.
Ma SC, Hao YJ, Jiao Y, Wang YH, Xu LB, Mao CY, Yang XL, Yang AN, Tian J, Zhang MH, Jin SJ, Xu H, Jiang YD, Zhang HP. Homocysteine induced oxidative stress through TLR4/NF κB/DNMT1 mediated LOX 1 DNA methylation in endothelial cells. Mol Med Rep. 2017;16:9181–8.
Leng YP, Ma YS, Li XG, Chen RF, Zeng PY, Li XH, Qiu CF, Li YP, Zhang Z, Chen AF. l-Homocysteine-induced cathepsin V mediates the vascular endothelial inflammation inhyperhomocysteinaemia. Br J Pharmacol. 2018;175:1157–72.
Lian Z, Lv FF, Yu J, Wang JW. The anti-inflammatory effect of microRNA-383-3p interacting with IL1R2 against homocysteine-induced endothelial injury in rat coronary arteries. J Cell Biochem. 2018;119:6684–94.
Chen CJ, Yang TC, Chang C, Lu SC, Chang PY. Homocysteine is a bystander for ST-segment elevation myocardial infarction: a case-control study. BMC Cardiovasc Disord. 2018;18:33.
Clarke R, Bennett DA, Parish S, Verhoef P, Dötsch-Klerk M, Lathrop M, Xu P, Nordestgaard BG, Holm H, Hopewell JC, Saleheen D, Tanaka T, Anand SS, Chambers JC, Kleber ME, Ouwehand WH, Yamada Y, Elbers C, Peters B, Stewart AF, Reilly MM, Thorand B, Yusuf S, Engert JC, Assimes TL, Kooner J, Danesh J, Watkins H, Samani NJ, Collins R, Peto R, MTHFR Studies Collaborative Group. Homocysteine and coronary heart disease: meta-analysis of MTHFR case-control studies, avoiding publication bias. PLoS Med. 2012;9:e1001177.
Han L, Wu Q, Wang C, Hao Y, Zhao J, Zhang L, Fan R, Liu Y, Li R, Chen Z, Zhang T, Chen S, Ma J, Liu S, Peng X, Duan S. Homocysteine, ischemic stroke, and coronary heart disease in hypertensive patients: APopulation-based. Prospective Cohort Study Stroke. 2015;46:1777–86.
Yusuf S, Hawken S, Ounpuu S, Dans T, Avezum A, Lanas F, McQueen M, Budaj A, Pais P, Varigos J, Lisheng L, INTERHEART Study Investigators. Effect of potentially modifiable risk factors associated with myocardial infarction in 52 countries(the INTERHEART study): case-control study. Lancet. 2004;364:937–52.
Stone NJ, Robinson JG, Lichtenstein AH, Bairey Merz CN, Blum CB, Eckel RH, Goldberg AC, Gordon D, Levy D, Lloyd-Jones DM, McBride P, Schwartz JS, Shero ST, Smith SC Jr, Watson K, Wilson PW, Eddleman KM, Jarrett NM, LaBresh K, Nevo L, Wnek J, Anderson JL, Halperin JL, Albert NM, Bozkurt B, Brindis RG, Curtis LH, DeMets D, Hochman JS, Kovacs RJ, Ohman EM, Pressler SJ, Sellke FW, Shen WK, Smith SC Jr, Tomaselli GF, American College ofCardiology/American Heart Association Task Force on Practice Guidelines. 2013 ACC/AHA guideline on the treatment of blood cholesterol to reduce atherosclerotic cardiovascular risk in adults: a report of the American College of Cardiology/American Heart Association task force on practice guidelines. Circulation. 2014;129:S1–45.
National Cholesterol Education Program (NCEP) Expert Panel on Detection, Evaluation, and Treatment of High Blood Cholesterol in Adults (Adult Treatment Panel III). Third Report of the National Cholesterol Education Program (NCEP) Expert Panel on Detection, Evaluation, and Treatment of High Blood Cholesterol in Adults (Adult Treatment Panel III) final report. Circulation. 2002;106:3143–421.
Baszczuk A, Thielemann A, Musialik K, Kopczynski J, Bielawska L, Dzumak A, Kopczynski Z, Wysocka E. 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.
Visram M, Radulovic M, Steiner S, Malanovic N, Eichmann TO, Wolinski H, Rechberger GN, Tehlivets O. Homocysteine regulates fatty acid and lipid metabolism in yeast. J Biol Chem. 2018;293:5544–55.
Vijayakumar A, Kim EK, Kim H, Choi YJ, Huh KB, Chang N. Effects of folic acid supplementation on serum homocysteine levels, lipid profiles, and vascular parameters in post-menopausal Korean women with type 2 diabetes mellitus. Nutr Res Pract. 2017;11:327–33.
Llanos AA, Marian C, Brasky TM, Dumitrescu RG, Liu Z, Mason JB, Makambi KH, Spear SL, Kallakury BV, Freudenheim JL, Shields PG. Associations between genetic variation in one-carbon metabolism and LINE-1 DNA methylation inhistologically normal breast tissues. Epigenetics. 2015;10:727–35.
Zhang L, Yin RX, Liu WY, Miao L, Wu DF, Aung LH, Hu XJ, Cao XL, Wu JZ, Pan SL. 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.
Pollex RL, Mamakeesick M, Zinman B, Harris SB, Hanley AJ, Hegele RA. Methylenetetrahydrofolate reductase polymorphism 677C>T is associated with peripheral arterial disease in type 2 diabetes. Cardiovasc Diabetol. 2005;4:17.
Maeda M, Yamamoto I, Fukuda M, Motomura T, Nishida M, Nonen S, Fujio Y, Kasayama S, Azuma J. MTHFR gene polymorphism is susceptible to diabetic retinopathy but not to diabetic nephropathyin Japanese type 2 diabetic patients. J Diabetes Complicat. 2008;22:119–25.
Lv CF, Zhang T, Li LL, Wang CY, Liu K. Correlation of metabolism-related indicators with MTHFR gene polymorphism in H type hypertension patients. China Trop Med. 2014;14:86–8.
Fan GJ, Xu R, Zhang Q, Yun L, Cao YL, Zhang CM. Correlation between methylenetetrahydrofolate reductase gene C677T polymorphism and blood lipid abnormality among hypertensives. Chin J Arterioscler. 2017;25:153–8.
DerSimonian R, Kacker R. Random-effects model for meta-analysis of clinical trials: an update. Contemp Clin Trials. 2007;28:105–14.
Begg CB, Mazumdar M. Operating characteristics of a rank correlation test for publication bias. Biometrics. 1994;50:1088–101.
This research was supported by the grants from the Key Project of Education Department of Sichuan Province, P. R. China (17ZA0172) and the Cooperative Project on Scientific Research between Nanchong city and North Sichuan Medical College, P. R. China (NSMC20170403).
This research was supported by the grants from the Key Project of Education Department of Sichuan Province, P. R. China (17ZA0172) and the Cooperative Project on Scientific Research between Nanchong city and North Sichuan Medical College, P. R. China (NSMC20170403).
Availability to data and materials
All data generated or analysed during this study are included in this published article [and its supplementary information files].
The authors declare that they have no competing interests.
Ethics approval and consent to participate
Consent for publication
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The reference list for the studies included in the present meta-analysis. (DOC 106 kb)
Table S1. Characteristics of the individual studies included in the meta-analysis between the MTHFR rs1801133 polymorphism and CAD; Table S2. Characteristics of the individual studies included in the meta-analysis between the MTHFR rs1801133 polymorphism and plasma lipid levels; Table S3. Plasma lipid levels according to the genotypes of the MTHFR rs1801133 polymorphism. (DOC 544 kb)