Skip to main content

Association of maternal folic acid supplementation and offspring MTRR gene polymorphism with congenital heart disease: a hospital-based case-control study in Han population

Abstract

Background

Although many studies shown that the risk of congenital heart disease (CHD) was closely related to genetic and environmental factors, the exact mechanism was still unclear. This study was to assess the association of maternal folic acid supplementation (FAS), the 5-methyltetrahydrofolate-homocysteine methyltransferase reductase (MTRR) gene polymorphisms in offspring and their interaction effects with the risk of CHD and its subtypes.

Methods

A case-control study was conducted on 595 children with CHD and 605 healthy child controls. The multivariate logistic regression model was used to assess the association of maternal FAS, offspring MTRR gene polymorphisms and their interaction effects with CHD and its subtypes.

Results

This study shown that maternal FAS was significantly associated with a reduced risk of CHD (OR = 0.55, 95%CI: 0.36–0.83) and its subtypes including ASD (OR = 0.25, 95%CI: 0.14–0.45), VSD (OR = 0.42, 95%CI: 0.27–0.64), and CTD (OR = 0.23, 95%CI: 0.09–0.59) in offspring. Offspring MTRR gene polymorphisms at rs162048 (GG vs. AA: OR = 2.05, 95%CI: 1.35–3.13), rs1802059 (AA vs. GG: OR = 5.13, 95%CI: 2.15–12.23; GA vs. GG: OR = 1.81, 95%CI: 1.35–2.43), rs10380 (TT vs. CC: OR = 2.27, 95%CI: 1.20–4.31) and rs1801394 (GG vs. AA: OR = 1.58, 95%CI: 1.02–2.42) were significantly associated with the risk of CHD, and similar results were also found for three subtypes of CHD. Additionally, a statistically significant interaction effect between maternal FAS and offspring MTRR gene polymorphism at rs1802059 was observed (OR = 0.38, 95%CI: 0.15–0.94). Among children who had a variant genotype at rs1802059, the risk of CHD was significantly decreased when their mother used folate for this pregnancy compared with mothers not using folate.

Conclusions

In those of Chinese descent, maternal FAS and offspring MTRR gene polymorphisms are significantly associated with the risk of CHD and its three subtypes. Furthermore, maternal FAS may help to offset some of risks of CHD due to offspring MTRR genetic variants. However, more studies with prospective designs and larger samples are needed to confirm our findings.

Trial registration

Registration number: ChiCTR1800016635; Registration time: 14/06/2018.

Introduction

Congenital heart disease (CHD), as one of the most common birth defects, refers to structural abnormalities of the heart and/or great vessels present at birth due to abnormal cardiac development during embryonic development [1]. A recent meta-analysis showed that the global incidence of CHD was 3.81‰ [2]. Existing evidence suggests that both genetic and environmental factors contribute to the development of CHD [3].

Folic acid, also known as vitamin B9, is involved in carbon metabolism and affects the synthesis of nucleic acids and proteins, etc., and is an indispensable basic material for the human body to maintain normal physiological functions [4]. More and more studies have pointed to the important role of folic acid in preventing a variety of birth defects, including neural tube defects, cleft lip and palate, and CHD [5,6,7,8]. Previous epidemiological studies on the association between maternal folic acid supplementation (FAS) during pregnancy and the risk of CHD in offspring were conducted in the United States [9], China [10], Norway [11, 12], and Hungary [13], but the results were inconsistent. This could be caused by differences in study design, biases, racial differences, and in dietary patterns. In addition, it was proposed that the strength of the association between folic acid and CHD was heterogeneous in different subtypes [14].

As a metabolite in the folate pathway, homocysteine is an independent risk factor for CHD [15]. The 5-methyltetrahydrofolate-homocysteine methyltransferase reductase (MTRR) gene is one of the key regulatory enzymes involved in homocysteine metabolic pathway. The MTRR gene which is located on the short arm of chromosome 5 at position 15.31, maintains 5-methyltetrahydrofolate-homocysteine methyltransferase (MTR) in its functional state and assists in the folate/VB12 dependent conversion of Homocysteine (Hcy) to methionine. The MTRR gene plays an important role in folate metabolism and the maintenance of Hcy levels, and its mutation may alter the concentration of homocysteine, thereby affecting normal embryonic development [15]. A meta-analysis showed that the MTRR gene polymorphism at rs1801394 was significantly associated with the risk of CHD [16]. Some epidemiological studies and animal studies shown that MTRR gene mutations were associated with an increased risk of CHD, and even led to embryonic lethality [17,18,19]. However, some studies suggested no relationship [20]. Additionally, most of previous studies only focus on two genetic loci (i.e., rs1801394 and rs1532268) of MTRR gene [16, 21,22,23]. The MTRR gene contains other common variants within the coding sequence in addition to rs1801394 and rs1532268. Further investigation of other genetic loci in the MTRR gene is necessary to reveal clues related to the metabolism underlying the potential embryonic protective effects of maternal FAS during pregnancy.

Therefore, we performed a hospital-based case-control study to comprehensively assess the association of maternal FAS, multiple genetic variants of the offspring MTRR gene, and their interaction effects with the risk of CHD and its subtypes of Chinese descent, which will help to promote the development of primary prevention strategies to reduce the burden of CHD worldwide.

Methods

Ethics statement

The study was approved by Ethics Committee for Clinical Research of Xiangya School of Public Health of Central South University (No. XYGW-2018-36), and followed the principles of the Declaration of Helsinki. Before the participants’ personal information collection and biological samples, we have obtained all mothers’ written informed consent, to ensure their voluntary participation and fully understand the research purpose and procedures. The protocol of this study was registered at the Chinese Clinical Trial Registry (registration number: ChiCTR1800016635; registration time: 14/06/2018). It is available at http://www.chictr.org.cn/listbycreater.aspx.

Study design and recruitment of study participants

We performed a hospital-based case-control study from November 2017 to December 2019 using the Hunan Provincial Children’s Hospital as a recruitment site. Children with CHD and their mothers from the Department of Cardiothoracic Surgery were identified as case group, and the healthy individuals were recruited from the Department of Child Healthcare during the same period. To minimize potential maternal recall bias about preconception and first-trimester related exposures, we recruited only children younger than one year of age. All cases were simple CHD that excluded malformations from other systems and were diagnosed by echocardiography as well as confirmed by surgery. All controls were confirmed by physical examination and cardiac ultrasound that the structure and function of the heart were normal, and there were no chromosomal abnormalities and systemic malformations. In addition, to minimize predisposing factors due to genetic and cultural differences, only Han subjects were recruited, and there were no familial relationships between the case and control groups. Subjects who were unable to provide a sample or were unwilling to cooperate in filling out the questionnaire were excluded from the study.

Information and biological sample collection

In addition to CHD, this study also focused on three important CHD subtypes: arterial septal defect (ASD), ventricular septal defect (VSD), and conotruncal defects (CTD). Among them, CTD specifically includes the following four subtypes: tetralogy of Fallot, double outlet ventricle, persistent truncus arteriosus and various transposition of the great arteries. According to World Health Organization (WHO) guidelines, FAS is defined as folic acid intake of more than 0.4 mg per day on at least five days per week continuously before conception or throughout the first trimester of pregnancy. In this study, professionally trained investigators conducted face-to-face interviews to collect data. According to our previous study [24], the following covariates were used as confounders in subsequent analyses: socio-demographic characteristics (residence, maternal age, education level, child gender), family history (consanguineous marriage, birth defects), history of pregnancy complications (gestational diabetes, gestational hypertension), maternal perinatal lifestyle (active smoking, passive smoking, drinking alcohol).

After enrollment, a trained nurse in the department collected 5 ml of peripheral venous blood using an EDTA anticoagulant tube. Collected blood samples were immediately centrifuged with a low-speed centrifuge at 3500 r/min for 15 min. Subsequently, isolated plasma and blood cells were carefully packaged, labeled, and stored in a cryogenic freezer at − 80 ° C. In this study, blood cells were used for single nucleotide polymorphism (SNP) detection because of their high sensitivity and reliability.

Questionnaires were conducted face-to-face between professionally trained investigators and mothers to collect information. The questionnaire used has been described in our previous published articles [24]. One of the exposures of interest was maternal FAS, which was defined as a mother who had taken folic acid in 3 months before pregnancy and/or during the first-trimester pregnancy. Based on the previous literature [24], we selected the potentially confounding factors of significance as follows: demographic characteristics of pregnant women and children (sex of child, maternal age at pregnancy onset, maternal residence location and maternal educational level); family history (consanguineous marriage and history of birth defects in family); maternal exposure in 3 months before pregnancy and/or during the first-trimester pregnancy (drinking alcohol, active smoking and passive smoking); history of pregnancy complications (gestational diabetes and gestational hypertension). In China, each pregnant woman has a perinatal care manual, which contains pregnancy and personal information. The information collected was checked for accuracy against their health manuals and related medical records, which recorded their basic demographic characteristics, illness and lots of physical and laboratory examination results every time they went to a pregnancy test.

Sequencing of MTRR gene SNPs

The MTRR gene was a candidate gene in this study. In short, SNP markers were selected using the SNPBrowserTM program (version 3.0) from AppliedBiosystems, Inc. The program allows SNP markers to be selected from the HapMap database. Sites with r2 less than 0.8 and minimum allele frequency (MAF) less than 10% were excluded. The MTRR gene polymorphisms were detected using the matter-assisted laser desorption ionization time-of-flight mass spectrometry MassARRAY (Agena iPLEX assay, San Diego, CA, USA). Laboratory technicians were responsible for genotyping, retyping and double-checking samples, and recording genotype data. We specified a minimum call rate of 50% for SNP genotyping to ensure the integrity of the individual genotype data that has been called. In addition, the experimenters did not know whether the samples came from the case group or the control group to ensure the integrity of the experiment. The success rate of SNPlex detection for these eleven SNPs is greater than 90%. Finally, these loci (rs162036, rs326120, rs162048, rs1802059, rs2287779, rs10380, rs16879334, rs1801394, rs3776455, rs2303080 and rs1532268) were selected as candidate loci in this study.

Statistical analysis

Categorical variables were described using frequencies and percentages. In the univariate analysis, the Pearson chi-square test or Fisher’s exact test was used to compare categorical variables; The Wilcoxon rank sum test was used for ordinal categorical variables. SNPs of the MTRR gene were analyzed for deviations from Hardy-Weinberg equilibrium (HWE) in the control group. A locus with a P value of less than 0.1 indicates an imbalance in gene frequency and is therefore eliminated. Three common genetic models were used: dominant model (Aa + AA vs. aa), recessive model (AA vs. Aa + aa), and additive model (AA vs. Aa vs. aa). Odds ratios (ORs) and their 95% confidence intervals (CIs) were used to measure the level of association. Unadjusted ORs were calculated by univariate logistic regression. Adjusted ORs were calculated by multivariable logistic regression. Logistic regression was used to examine the main effects and interactive effects of the gene-environment interaction effect of maternal FAS and the offspring MTRR gene on CHD by controlling for other influencing factors. Bioinformatics analysis included: analysis of GTEx and based on STRING protein-protein interaction network. Statistical analysis was performed using SPSS (IBM SPSS Statistics for Macintosh, Version 26.0. Armonk, NY: IBM Corp) and R software, version 4.2.2 (R Foundation for Statistical Computing, Vienna, Austria.) 4.2.2. All tests were performed significantly for a two-sided P value not exceeding 0.05, except where otherwise specified.

Results

Comparison of baseline characteristics across groups

From November 2017 to December 2019, a total of 595 eligible CHD infants and their mothers were recruited into the case group. Among 595 CHD cases, 126 (21.2%) were diagnosed with ASD, 428 (71.9%) with VSD, 31 (5.2%) with CTD, 125 (21.0%) with other CHD subtypes. Of note, some cases have been diagnosed with multiple subtypes of CHD. Therefore, the sum of the various subtypes was not equal to 595. During the same period, 605 eligible healthy children and their mothers were identified as the control group. Baseline characteristics were summarized in Table 1.

Overall, there were no statistically significant differences across groups for maternal pregnancy age. However, our study showed that there were statistically significant differences between two groups for the following characteristics: Sex of child, residence location, maternal education level (years), history of consanguineous marriage, history of birth defects in family, gestational hypertension, gestational diabetes during this pregnancy, and personal lifestyle and habit in the periconceptional period including drinking alcohol, active smoking and passive smoking (all P values < 0.05). Therefore, these factors were adjusted when accessing the association of maternal FAS, offspring MTRR gene polymorphisms, and their interaction effect with the risk of CHD in offspring.

Table 1 Comparison of baseline characteristics between case and control groups

Maternal FAS and risk of CHD and its subtypes in offspring

Associations of maternal FAS with the risk of CHD and its subtypes were summarized in Table 2. The proportion of mothers of cases (total CHD:19.16%; ASD: 23.81%; VSD: 15.19%; CTD: 25.81%) who did not take FAS during pregnancy was significantly higher than that of the controls (7.93%). After adjustment for baseline data, mothers taking FAS for this pregnancy, compared with those not taking FAS, were at a significantly decreased risk of total CHD (aOR = 0.55; 95% CI: 0.36–0.83) and its subtypes including ASD (aOR = 0.25, 95%CI: 0.14–0.45), VSD (aOR = 0.42, 95%CI: 0.27–0.64), and CTD (aOR = 0.23, 95%CI: 0.09–0.59) in offspring.

Table 2 Association of maternal FAS with the risk of CHD and its subtypes in offspring

Offspring MTRR gene polymorphisms and risk of CHD and its subtypes

Genotype frequencies for each SNP of the offspring MTRR gene and P values of the HWE test in the control group were summarized in Table S1. The genotype distributions in the control group were within HWE (P > 0.1) except for rs162036, rs326120 and rs1532268. The results of intergroup distribution of alleles/genotypes of offspring MTRR locus were shown in Table S2. The genotype and allele distributions of rs162048, rs1802059, rs10380 and rs3776455 as well as the genotype distribution of rs1801394 showed statistically significant differences between the two groups (all P < 0.05), while the genotype and allele distributions of the other loci did not show statistically significant differences. As a consequence, the association of the MTRR gene polymorphism at rs162048, rs1802059, rs10380, rs1801394 and rs3776455 with the risk of CHD and its subtypes was assessed in the subsequent analyses.

Associations between offspring MTRR gene polymorphism and the risk of CHD and its subtypes in offspring based on logistic regression analyses were summarized in Table 3. Multivariable logistic regression showed that offspring MTRR gene polymorphisms at rs162048, rs1802059, rs10380 and rs1801394 were significantly associated with the risk of total CHD.

Table 3 Associations between polymorphisms of MTRR gene and the risk of CHD and its subtypes

For rs162048, offspring with the GG genotypes, compared with those with the AA genotype, had a significantly higher risk of total CHD (aOR = 2.05; 95% CI: 1.35–3.13); in addition, the dominant model (aOR = 1.53; 95% CI: 1.16–2.01), recessive model (aOR = 1.75; 95% CI: 1.18–2.59) and additive model (aOR = 1.40; 95% CI: 1.15–1.71) were significantly associated with an increased risk of total CHD. The results of the three subtypes were similar to those of the total CHD group. Specifically, the genotype distribution of rs162048 was statistically significant in all three subtypes, but the dominant/recessive/additive model of rs162048 in the ASD group was not statistically significant.

For rs1802059, offspring with the GA genotype (aOR = 1.81; 95% CI: 1.35–2.43), and AA (aOR = 5.13; 95% CI: 2.15–12.23) genotype were at a significantly increased risk of total CHD compared with those with the GG genotype; besides, the dominant model (aOR = 1.90; 95% CI: 1.43–2.51), recessive model (aOR = 4.29; 95% CI: 1.81–10.16) and additive model (aOR = 1.94; 95% CI: 1.06–2.49) all significantly increased the risk of total CHD. Results in the ASD and VSD groups were similar to total CHD group, while only the additive model was statistically significant in the CTD group.

For rs10380, offspring with the TT genotypes, compared with those with the CC genotype, had a significantly higher risk of total CHD (aOR = 2.27; 95% CI: 1.20–4.31); and recessive model (aOR = 2.25; 95% CI: 1.19–4.24) were significantly associated with an increased risk of total CHD; MTRR gene polymorphisms at rs10380 were associated with the risk of disease in the ASD (CT vs. CC; TT vs. CC; dominant model; recessive model and additive model), VSD (TT vs. CC and additive model), and CTD (TT vs. CC and additive model) groups.

For rs1801394, offspring with the GG genotypes, compared with those with the AA genotype, had a significantly higher risk of total CHD (aOR = 1.58; 95% CI: 1.02–2.42); and additive model (aOR = 1.25; 95% CI: 1.03–1.53) were significantly associated with an increased risk of total CHD; MTRR gene polymorphisms at rs1801394 were associated with the risk of disease in the ASD (AG vs. AA; dominant model and recessive model) and VSD groups (GG vs. AA), but not in the CTD group.

Interaction effects between maternal FAS and offspring MTRR gene polymorphisms

Interaction effects between maternal FAS and offspring MTRR gene polymorphisms (rs162048, rs1802059, rs10380 and rs1801394) on the risk of total CHD were summarized in Table 4. Multivariable logistic regression suggested that there were statistically significant interaction effects between maternal FAS and offspring MTRR gene polymorphisms at rs1802059 (aOR = 0.38, 95%CI: 0.15–0.94; P = 0.035). Additionally, we further analyzed the independent effect of offspring MTRR gene polymorphisms at rs1802059 on the risk of total CHD by stratification of maternal FAS (Table 5). The results suggested that the risk of total CHD was significantly decreased among offspring who had variant genotypes but maternal FAS when compared with those who had variant genotypes and no maternal FAS. Without maternal FAS, the risk of CHD for offspring with the variant genotype was 4.41 (95% CI: 1.91–10.17) for rs1802059. However, with maternal FAS, the risk of CHD for offspring with the variant genotype was 1.66 (95% CI: 1.24–2.24) for rs1802059.

Table 4 Interaction effects between maternal FAS and offspring MTRR gene polymorphisms on the risk of total CHD
Table 5 Genotype of offspring MTRR gene polymorphisms at rs1802059 by stratification of maternal FAS and risk (odds ratio) of total CHD

Bioinformatics analysis

For the MTRR gene polymorphisms at rs162048, rs1802059, rs10380, rs1801394 with statistical significance, the relevant expression quantitative trait locus (eQTL) data were searched and downloaded from the GTEx database (https://www.gtexportal.org/home/), and the results were shown in Table S3. Except for rs1802059, significant associations were observed between rs162048, rs10380, rs1801394 and MTRR gene expression in whole blood. In addition, these sites were also statistically correlated with MTRR gene expression levels in other human tissues, such as brain tissue, adipose tissue, muscle tissue, and skin tissue.

Figure 1 showed the interactions between MTRR proteins and other proteins using the STRING 11.0 protein-protein interaction network database. Table 6 listed the combined scores of MTRR proteins and predicted protein interactions. From the composite scores, MTRR showed the strongest evidence of interaction with MTR, followed by Methylenetetrahydrofolate reductase (MTHFR).

Fig. 1
figure 1

Diagram of protein-protein interaction networks associated with MTRR proteins

Table 6 Composite score of MTRR proteins and predicted protein interaction

Discussion

Our study showed that maternal FAS, and offspring MTRR gene polymorphisms were significantly associated with the risk of CHD and its subtypes of Chinese descent. Additionally, there were statistically significant interaction effects between maternal FAS and offspring MTRR gene polymorphisms on the risk of CHD in offspring. As far as we know, this is the first time that the association of maternal FAS, offspring MTRR gene polymorphisms, and their interaction effects with the risk of CHD has been exhaustively explored, which highlight the importance of periconception FAS and may help provide new clues for future etiological studies and interventions in CHD.

In this study, we observed that maternal FAS significantly reduced the risk of CHD (OR = 0.55) in the offspring. The results of previous studies on this topic were consistent with our findings. For example, Wang et al. [25] analyzed the data from the birth defect surveillance system in a district of Beijing and found that periconception maternal FAS could significantly reduce the risk of CHD in their offspring. Additionally, case-control studies from the United States [26], New Zealand [27], and a birth cohort study from Gansu Province, China [28] also observed a protective effect of maternal FAS against multiple types of CHD. It is well known that maternal FAS during pregnancy has been the main preventive measure to reduce the risk of neural tube malformations in offspring [5]. In recent years, more and more studies have shown that maternal FAS also has a preventive effect on CHD in offspring [8, 10]. Besides, there were studies claiming that folic acid can prevent birth defects such as cleft lip and palate, as well as cardiovascular diseases and cancers [6, 7]. The present study also found that the protective effect of maternal FAS on CHD subtypes such as ASD (OR = 0.25), VSD (OR = 0.42) and CTD (OR = 0.23) was stronger than that of total CHD, and the same results were observed in the studies of Morikawa et al. [14] and Qu et al. [29]. Folic acid played an irreplaceable role in one-carbon metabolism and was essential for the synthesis of many substances such as proteins, deoxyribonucleic acid, neurotransmitters and phospholipids [30,31,32]. Folic acid deficiency could affect this process, which could lead to the development of a variety of diseases. Therefore, our findings emphasized that women of childbearing age should raise their awareness of FAS and ensure adequate intake of folic acid during periconception to prevent CHD in offspring.

This study also assessed the association of MTRR gene polymorphisms with the risk of CHD. The results suggested that polymorphisms of the MTRR gene at rs162048, rs1802059, rs10380 and rs1801394 were significantly associated with the risk of CHD and its subtypes. Previous studies reached the same conclusion [33,34,35], and proposed that MTRR gene polymorphism might be associated with acyanotic CHD [36]. It was worth noting that there were also some studies suggesting that the association between MTRR gene polymorphism and CHD was not statistically significant. For example, a case-control study from the Netherlands did not find a significant association between offspring MTRR rs1801394 polymorphism and the development of heart defects in children [20]. The inconsistent results may be explained by racial differences, differences in study design, and the presence of bias. In addition, several studies confirmed that MTRR gene polymorphisms were also associated with a variety of cancers [37,38,39], reproductive disorders [40] and cleft lip and palate [41]. However, most of the previous studies only focused on two loci (i.e., rs1801394 and rs1532268). Our study was the first to explore the relationship between five polymorphisms (rs162048, rs1802059, rs10380, rs1801394 and rs3776455) of MTRR gene and the risk of CHD and its subtypes, which provided new clues for further research on MTRR gene variants associated with CHD. MTRR gene was one of the key regulatory enzymes involved in homocysteine metabolism pathway, and homocysteine was an independent risk factor for CHD [15], so there was a high possibility that MTRR gene polymorphisms might be associated with CHD. MTRR maintained sufficient levels of activated cobalamin, which served as a cofactor for MTR. In the process of MTR catalyzed re-methylation of homocysteine to methionine, cobalamin acted as an intermediated methyl carrier between methyltetrahydrofolate and homocysteine. The cobalamin cofactor cycled between the cob (I) alamin and methyl cob (III) alamin, but the cob (I) alamin could be oxidized to the unactivated cob (II) alamin form, and for it to regain activity, cob (II) needed to be converted to the methyl cob (III) alamin form by obtaining a methyl donor for S-adenosylmethionine catalyzed by MTRR. This cycle ensured the activity of MTR, and MTRR acted as a “companion” that played an important role in keeping MTR in an active state [42]. Therefore, MTRR gene mutation might change the concentration of homocysteine and thus affected the normal development of embryonic heart, and studies confirmed this speculation [15]. However, the specific mechanism will still require further studies to elucidate. In addition, bioinformatics analysis showed that the evidence for the interaction between MTRR and MTR was strongest, with a combined score of 0.999 for both, followed by sub-MTHFR with a combined score of 0.981. MTR catalyzes the transfer of a methyl group from methyl-cobalamin to homocysteine, yielding enzyme-bound cob(I)alamin and methionine, while MTRR is involved in the reductive regeneration of cob(I)alamin cofactor. The utilization of methylgroups in the folate cycle is necessary, and these proteins play an important role in methylation. In addition, MTHFR also plays an important role in the folate cycle by catalyzing the conversion of 5, 10-methyltetrahydrofolate to 5-methyltetrahydrofolate, a co-substrate for homocysteine remethylation to methionine. Therefore, we believed that MTRR had a strong interaction with MTR and MTHFR, which were closely related in folate metabolism.

In this study, we also found an interaction effect between maternal FAS and offspring MTRR gene polymorphisms at rs1802059 in the development of CHD. This study showed that children who carried the variant genotype and with maternal FAS had a significantly lower risk of developing CHD compared with children who carried the variant genotype but without maternal FAS. This suggested that maternal FAS might help to offset some of the risk of developing CHD due to variations in the MTRR gene in offspring. The interaction effect between maternal FAS and offspring MTRR gene polymorphisms in CHD had not been investigated before, and this study was the first to find this association. Previous literatures indicated that the insufficient activity of MTRR and MTR caused by MTRR gene variation will decrease the utilization rate of folic acid and increase the level of Hcy, which might eventually induce vascular endothelial injury and cardiovascular diseases [43, 44]. Additionally, insufficient maternal folic acid intake may affect the methylation of genes related to fetal growth and development [45,46,47]. We speculate that this might be a mechanism of interaction between maternal FAS and offspring MTRR gene polymorphisms. However, the exact mechanism remained unclear and required further study. More and more evidence show that the interaction of gene and environment could affect CHD, so the influence of interaction would be one of the future research directions of CHD.

There were some limitations in this study. First, selection bias was hard to avoid. This study was a hospital-based case-control study, and our target population only included those children with CHD who were born successfully and survived. We did not know the cases of termination of pregnancy due to CHD or death after birth due to severe CHD, which may lead to a series of problems, such as the representativeness of the sample. Second, this study classified the exposure to maternal FAS as “yes” and “no”, and did not investigate other ways of folic acid intake, such as dietary conditions, so that the true folic acid level of each pregnant woman during the perinatal period could not be obtained, which may affect the study results to some extent. Third, residual confounding was of concern, even though we used multivariate logistic regression to control for the effect of confounding. Fourth, previous studies shown MTRR gene polymorphisms at rs1532268 was associated with the risk of CHD [16]. However, in this study, the distribution of this locus in the control population did not meet the HWE test, so we did not analyze this locus. In addition, as this study included only Han Chinese, more studies in other populations need to generalize the findings. Moreover, there were no animal experimental studies on MTRR gene variants and the risk of congenital heart disease, and the functional evidence of the effects of these gene variants could not be directly obtained. These limitations suggest that future studies should conduct prospective cohort studies with larger sample sizes in different ethnic populations and experimental study.

Conclusion

In conclusion, the present study further confirms that maternal FAS is significantly associated with the risk of CHD and its subtypes in offspring of Chinese descent. This study also shows that MTRR gene polymorphisms at rs162048, rs1802059, rs10380 and rs1801394 are significantly associated with the risk of CHD and its subtypes. Additionally, interaction effect analyses indicate that maternal FAS can help to offset some of the risks of CHD due to genetic variants of offspring MTRR gene. However, the mechanism by which these factors relate to fetal heart development remained unknown, and more studies with prospective designs and larger samples as well as experimental studies need to confirm our findings.

Data availability

No datasets were generated or analysed during the current study.

Abbreviations

CHD:

Congenital heart disease

FAS:

folic acid supplementation

MTRR:

5-methyltetrahydrofolate-homocysteine methyltransferase reductase

VSD:

ventricular septal defect

MTR:

5-methyltetrahydrofolate-homocysteine methyltransferase

Hcy:

Homocysteine

ASD:

arterial septal defect

CTD:

conotruncal defects

SNPs:

Single nucleotide polymorphisms

MAF:

minimum allele frequency

HWE:

Hardy-Weinberg equilibrium

aOR:

Adjusted odds ratio

CI:

Confidence interval

MTHFR:

Methylenetetrahydrofolate reductase

eQTL:

expression quantitative trait locus

References

  1. Triedman JK, Newburger JW. Trends in congenital heart disease: the Next Decade. Circulation. 2016;133(25):2716–33.

    Article  PubMed  Google Scholar 

  2. Liu Y, Chen S, Zuhlke L, Babu-Narayan SV, Black GC, Choy MK, Li N, Keavney BD. Global prevalence of congenital heart disease in school-age children: a meta-analysis and systematic review. BMC Cardiovasc Disord. 2020;20(1):488.

    Article  PubMed  PubMed Central  Google Scholar 

  3. van der Bom T, Zomer AC, Zwinderman AH, Meijboom FJ, Bouma BJ, Mulder BJ. The changing epidemiology of congenital heart disease. Nat Rev Cardiol. 2011;8(1):50–60.

    Article  PubMed  Google Scholar 

  4. Molloy AM. Folate bioavailability and health. Int J Vitam Nutr Res. 2002;72(1):46–52.

    Article  CAS  PubMed  Google Scholar 

  5. Murray LK, Smith MJ, Jadavji NM. Maternal oversupplementation with folic acid and its impact on neurodevelopment of offspring. Nutr Rev. 2018;76(9):708–21.

    Article  PubMed  Google Scholar 

  6. Shulpekova Y, Nechaev V, Kardasheva S, Sedova A, Kurbatova A, Bueverova E, Kopylov A, Malsagova K, Dlamini JC, Ivashkin V. The Concept of Folic Acid in Health and Disease. Molecules 2021, 26(12).

  7. Rahimi S, Martel J, Karahan G, Angle C, Behan NA, Chan D, MacFarlane AJ, Trasler JM. Moderate maternal folic acid supplementation ameliorates adverse embryonic and epigenetic outcomes associated with assisted reproduction in a mouse model. Hum Reprod. 2019;34(5):851–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Rosenquist TH. Folate, Homocysteine and the cardiac neural crest. Dev Dynam. 2013;242(3):201–18.

    Article  CAS  Google Scholar 

  9. Shaw GM, O’Malley CD, Wasserman CR, Tolarova MM, Lammer EJ. Maternal periconceptional use of multivitamins and reduced risk for conotruncal heart defects and limb deficiencies among offspring. Am J Med Genet. 1995;59(4):536–45.

    Article  CAS  PubMed  Google Scholar 

  10. Qin JB, Li JQ, Li F, Sun MT, Wang TT, Diao JY, Zhang SM, Luo L, Li YH, Chen LT, et al. Association of maternal folate use and reduced folate carrier gene polymorphisms with the risk of congenital heart disease in offspring. Eur J Pediatr. 2021;180(10):3181–90.

    Article  CAS  PubMed  Google Scholar 

  11. Oyen N, Olsen SF, Basit S, Leirgul E, Strom M, Carstensen L, Granström C, Tell GS, Magnus P, Vollset SE et al. Association between maternal folic acid supplementation and congenital heart defects in offspring in birth cohorts from Denmark and Norway. J Am Heart Assoc 2019, 8(6).

  12. Leirgul E, Gildestad T, Nilsen RM, Fomina T, Brodwall K, Greve G, Vollset SE, Holmstrom H, Tell GS, Oyen N. Periconceptional Folic Acid supplementation and infant risk of congenital heart defects in Norway 1999–2009. Paediatr Perinat Ep. 2015;29(5):391–400.

    Article  Google Scholar 

  13. Csáky-Szunyogh M, Vereczkey A, Kósa Z, Gerencsér B, Czeizel AE. Risk and protective factors in the origin of Conotruncal defects of Heart-A Population-based case-control study. Am J Med Genet A. 2013;161(10):2444–52.

    Article  Google Scholar 

  14. Morikawa Y, Cserjesi P. Cardiac neural crest expression of Hand2 regulates Outflow and Second Heart Field Development. Circ Res. 2008;103(12):1422–U1168.

    Article  CAS  PubMed  Google Scholar 

  15. Zhao JY, Yang XY, Shi KH, Sun SN, Hou J, Ye ZZ, Wang J, Duan WY, Qiao B, Chen YJ, et al. A functional variant in the cystathionine β-synthase gene promoter significantly reduces congenital heart disease susceptibility in a Han Chinese population. Cell Res. 2013;23(2):242–53.

    Article  CAS  PubMed  Google Scholar 

  16. Xu AP, Wang WP, Jiang XL. The roles of MTRR and MTHFR gene polymorphisms in congenital heart diseases: a meta-analysis. Bioscience Rep 2018, 38.

  17. Zhao JY, Yang XY, Gong XH, Gu ZY, Duan WY, Wang J, Ye ZZ, Shen HB, Shi KH, Hou J, et al. Functional variant in Methionine Synthase Reductase Intron-1 significantly increases the risk of congenital heart disease in the Han Chinese Population. Circulation. 2012;125(3):482.

    Article  CAS  PubMed  Google Scholar 

  18. Elmore CL, Wu XC, Leclerc D, Watson ED, Bottiglieri T, Krupenko NI, Krupenko SA, Cross JC, Rozen R, Gravel RA, et al. Metabolic derangement of methionine and folate metabolism in mice deficient in methionine synthase reductase. Mol Genet Metab. 2007;91(1):85–97.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Swanson DA, Liu ML, Baker PJ, Garrett L, Stitzel M, Wu JM, Harris M, Banerjee R, Shane B, Brody LC. Targeted disruption of the methionine synthase gene in mice. Mol Cell Biol. 2001;21(4):1058–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Clin Chem Lab Med 2006, 44(11):1317–1323.

  21. Christensen KE, Zada YF, Rohlicek CV, Andelfinger GU, Michaud JL, Bigras JL, Richter A, Dube MP, Rozen R. Risk of congenital heart defects is influenced by genetic variation in folate metabolism. Cardiol Young. 2013;23(1):89–98.

    Article  PubMed  Google Scholar 

  22. Yu D, Yang L, Shen S, Fan C, Zhang W, Mo X. Association between methionine synthase reductase A66G polymorphism and the risk of congenital heart defects: evidence from eight case-control studies. Pediatr Cardiol. 2014;35(7):1091–8.

    Article  PubMed  Google Scholar 

  23. Karas Kuzelicki N, Doljak B. Congenital Heart Disease and Genetic Changes in Folate/Methionine cycles. Genes (Basel) 2024, 15(7).

  24. Song XL, Li QX, Diao JY, Li JQ, Li YH, Zhang SM, Zhao LJ, Chen LT, Wei JH, Shu J et al. Association of gene polymorphisms and maternal smoking with risk of congenital heart disease: a hospital-based case-control study. Bmc Pregnancy Childb 2022, 22(1).

  25. Wang D, Jin L, Zhang J, Meng WY, Ren AG. Maternal periconceptional folic acid supplementation and risk for fetal congenital heart defects. J Pediatr-Us. 2022;240:72–8.

    Article  CAS  Google Scholar 

  26. Botto LD, Mulinare J, Erickson JD. Occurrence of congenital heart defects in relation to maternal multivitamin use. Am J Epidemiol. 2000;151(9):878–84.

    Article  CAS  PubMed  Google Scholar 

  27. van Beynum IM, Kapusta L, Bakker MK, den Heijer M, Blom HJ, de Walle HEK. Protective effect of periconceptional folic acid supplements on the risk of congenital heart defects: a registry-based case-control study in the northern Netherlands. Eur Heart J. 2010;31(4):464–71.

    Article  PubMed  Google Scholar 

  28. Mao BH, Qiu J, Zhao N, Shao YW, Dai W, He XC, Cui HM, Lin XJ, Lv L, Tang ZF et al. Maternal folic acid supplementation and dietary folate intake and congenital heart defects. PLoS ONE 2017, 12(11).

  29. Qu YJ, Lin S, Zhuang J, Bloom MS, Smith M, Nie ZQ, Mai JZ, Ou YQ, Wu Y, Gao XM et al. First-trimester maternal folic acid supplementation reduced risks of severe and most congenital heart diseases in offspring: a large case-control study. J Am Heart Assoc 2020, 9(13).

  30. Blount BC, Mack MM, Wehr CM, MacGregor JT, Hiatt RA, Wang G, Wickramasinghe SN, Everson RB, Ames BN. Folate deficiency causes uracil misincorporation into human DNA and chromosome breakage: implications for cancer and neuronal damage. Proc Natl Acad Sci USA. 1997;94(7):3290–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Menezo Y, Elder K, Clement A, Clement P. Folic acid, Folinic Acid, 5 Methyl TetraHydroFolate supplementation for mutations that affect Epigenesis through the Folate and one-Carbon cycles. Biomolecules 2022, 12(2).

  32. Donnelly JG. Folic acid. Crit Rev Clin Lab Sci. 2001;38(3):183–223.

    Article  CAS  PubMed  Google Scholar 

  33. Zeng W, Liu L, Tong Y, Liu HM, Dai L, Mao M. A66G and C524T polymorphisms of the methionine synthase reductase gene are associated with congenital heart defects in the Chinese Han population. Genet Mol Res. 2011;10(4):2597–605.

    Article  CAS  PubMed  Google Scholar 

  34. Su J, Li ZZ. Analysis of and Gene Polymorphisms in Chinese patients with ventricular septal defect. Appl Immunohisto M M. 2018;26(10):769–74.

    CAS  Google Scholar 

  35. Pishva SR, Vasudevan R, Etemad A, Heidari F, Komara M, Ismail P, Othman F, Karimi A, Sabri MR. Analysis of and Gene Polymorphisms in Iranian ventricular septal defect subjects. Int J Mol Sci. 2013;14(2):2739–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Hassan FM, Khattab AA, Abo El Fotoh WMM, Zidan RS. And polymorphisms of methionine synthase reductase gene are linked to the development of acyanotic congenital heart diseases in Egyptian children. Gene. 2017;629:59–63.

    Article  CAS  PubMed  Google Scholar 

  37. Metayer C, Scélo G, Chokkalingam AP, Barcellos LF, Aldrich MC, Chang JS, Guha N, Urayama KY, Hansen HM, Block G, et al. Genetic variants in the folate pathway and risk of childhood acute lymphoblastic leukemia. Cancer Cause Control. 2011;22(9):1243–58.

    Article  Google Scholar 

  38. Zhong GP, Luo XJ, Li J, Liao YH, Gui G, Sheng JW. MTRR rs1532268 polymorphism and gastric cancer risk: evidence from a meta-analysis. J Int Med Res 2022, 50(5).

  39. Wu MH, Chen CH, Chen CP, Huang TL, Yueh TC, Wang ZH, Tsai CW, Pei JS, Mong MC, Yang YC, et al. Contribution of 5-Methyltetrahydrofolate-homocysteine methyltransferase reductase genotypes to Colorectal Cancer in Taiwan. Anticancer Res. 2022;42(5):2375–82.

    Article  CAS  PubMed  Google Scholar 

  40. Xu WH, Zhang L, Wu XL, Jin F. Association between Methionine Synthase Reductase A66G polymorphism and male infertility: a Meta-analysis. Crit Rev Eukar Gene. 2017;27(1):37–46.

    Article  CAS  Google Scholar 

  41. Li QY, Xu LD, Jia XY, Saleem K, Zaib T, Sun WJ, Fu SB. SNPs in folate pathway are associated with the risk of nonsyndromic cleft lip with or without cleft palate, a meta-analysis. Bioscience Rep 2020, 40.

  42. Mascarenhas R, Gouda H, Ruetz M, Banerjee R. Human B12-dependent enzymes: methionine synthase and Methylmalonyl-CoA mutase. Methods Enzymol. 2022;668:309–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Smith AD, Refsum H. Homocysteine - from disease biomarker to disease prevention. J Intern Med. 2021;290(4):826–54.

    Article  CAS  PubMed  Google Scholar 

  44. Hasan T, Arora R, Bansal AK, Bhattacharya R, Sharma GS, Singh LR. Disturbed homocysteine metabolism is associated with cancer. Exp Mol Med 2019, 51.

  45. Qian YY, Huang XL, Liang H, Zhang ZF, Xu JH, Chen JP, Yuan W, He L, Wang L, Miao MH, et al. Effects of maternal folic acid supplementation on gene methylation and being small for gestational age. J Hum Nutr Diet. 2016;29(5):643–51.

    Article  PubMed  Google Scholar 

  46. Kim JM, Hong K, Lee JH, Lee S, Chang N. Effect of folate deficiency on placental DNA methylation in hyperhomocysteinemic rats. J Nutr Biochem. 2009;20(3):172–6.

    Article  CAS  PubMed  Google Scholar 

  47. Calzada-Dávila M, Calvo-Anguiano G, Martínez-de-Villarreal LE, Lugo-Trampe JJ, González-Peña SM. Ancer-Rodríguez PR, Hernández-Almaguer MD, Campos-Acevedo LD: congenital heart diseases: genetic risk variants and their methylation status. Genes-Basel 2022, 13(11).

Download references

Acknowledgements

The authors thank all pediatricians and parents who participated in the study.

Author information

Authors and Affiliations

Authors

Contributions

L L, Y C and Q C participated in Investigation; L L, J O and T W participated in Methodology; L L, Q C and Y Z participated in Software; L L participated in Formal analysis and Writing – original draft; J O participated in Conceptualization; J O, Y G and T W participated in Data curation; J O and Y G participated in Writing – review & editing; Y C participated in Resources; Q C and M L participated in Validation; M L and T W participated in Visualization; Y Z and J B participated in Funding acquisition and Supervision; Jiabi Qin participated in Project administration; All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Yingfan Zhang or Jiabi Qin.

Ethics declarations

Ethics approval and consent to participate

This study was performed in line with the principles of the Declaration of Helsinki. Approval was granted by the Ethics Committee of Xiangya School of Public Health Central South University (No. XYGW-2018-36). All mothers provided their written informed consent to participate in this study.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Funding

Funding for this study was provided by the Hunan Provincial Natural Science Foundation Youth Fund Project (grant number: 2023JJ41022); the National Natural Science Foundation Program of China (grant numbers 82073653, 81803313); Hunan Outstanding Youth Fund Project (grant number 2022JJ10087).

Additional information

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, L., Ou, J., Chen, Y. et al. Association of maternal folic acid supplementation and offspring MTRR gene polymorphism with congenital heart disease: a hospital-based case-control study in Han population. J Health Popul Nutr 43, 220 (2024). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s41043-024-00699-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s41043-024-00699-w

Keywords