Baseline characteristics of study population
After considering the inclusion criteria, we finally recruited 464 CHD cases and their parents into the case group and 504 health infants and their parents into the control group. Comparisons of baseline characteristics across groups were summarized in Table 1. Our study showed that there were statistically significant differences between two groups for the following characteristics: maternal education level (years), annual income in the past 1 year (RMB), history of adverse pregnancy outcomes, consanguineous marriage, history of congenital malformations in family, cold or fever in the periconceptional period, and personal lifestyle and habit in the periconceptional period including drinking alcohol, drinking tea, living near environmental pollution source, dyeing hair or perming and folate use (all P values < 0.05). Thus, these factors were adjusted when accessing the association of maternal tobacco exposure, the genetic variants of the maternal MTHFD gene, and their interactions with the risk of CHD in offspring.
Table 1. Baseline characteristics in case and control groupsa
Baseline characteristics
|
Control group (n=504)
|
Case group (n=464)
|
Univariate analysisd
|
Age at pregnancy onset (years)
|
|
|
χ2=0.191; P = 0.662
|
<35
|
434(86.1%)
|
404(87.1%)
|
|
≥35
|
70(13.9%)
|
60(12.9%)
|
|
Residence location
|
|
|
χ2 =39.390; P = 0.662
|
Rural areas
|
276(54.8%)
|
344(74.1%)
|
|
Urban areas
|
228(45.2%)
|
120(25.9%)
|
|
Education level (years)
|
|
|
Z =12.306; P <0.001b
|
≤9
|
6(1.2%)
|
66(14.2%)
|
|
9-12
|
100(19.8%)
|
190(40.9%)
|
|
12-16
|
168(33.3%)
|
130(28.0%)
|
|
>16
|
230(45.6%)
|
78(16.8%)
|
|
Annual income in the past 1 year (RMB)
|
|
Z = 15.946; P <0.001b
|
≤50,000
|
144(28.6%)
|
372(80.2%)
|
|
50,000-100,000
|
216(42.9%)
|
68(14.7%)
|
|
100,000-150,000
|
46(9.1%)
|
10(2.2%)
|
|
>150,000
|
98(19.4%)
|
14(3.0%)
|
|
History of adverse pregnancy outcomes
|
|
χ2 = 12.033; P = 0.001
|
No
|
280(55.6%)
|
206(44.4%)
|
|
Yes
|
224(44.4%)
|
258(55.6%)
|
|
Consanguineous marriage χ2 = 14.480; P<0.001
|
No
|
502(99.6%)
|
446(96.1%)
|
|
Yes
|
2(0.4%)
|
18(3.9%)
|
|
History of congenital heart disease in family
|
|
χ2 = 20.759; P <0.001
|
No
|
500(99.2%)
|
436(94.0%)
|
|
Yes
|
4(0.8%)
|
28(6.0%)
|
|
Cold or feverc
|
|
|
χ2 = 16.513; P <0.001
|
No
|
446(88.5%)
|
366(78.9%)
|
|
Yes
|
58(11.5%)
|
98(21.1%)
|
|
Drinking alcoholc
|
|
|
χ2 = 9.060; P = 0.003
|
No
|
468(92.9%)
|
404(87.1%)
|
|
Yes
|
36(7.1%)
|
60(12.9%)
|
|
Drinking teac
|
|
|
χ2 = 9.257; P = 0.002
|
No
|
402(79.8%)
|
404(87.1%)
|
|
Yes
|
102(20.2%)
|
60(12.9%)
|
|
Living near environmental pollution sourcec
|
|
χ2 = 38.443; P <0.001
|
No
|
470(93.3%)
|
370(79.7%)
|
|
Yes
|
34(6.7%)
|
94(20.3%)
|
|
Dyeing hair or permingc
|
|
|
χ2 = 12.532; P <0.001
|
No
|
474(94.0%)
|
406(87.5%)
|
|
Yes
|
30(6.0%)
|
58(12.5%)
|
|
Folatec
|
|
|
χ2 = 23.917; P <0.001
|
No
|
34(6.7%)
|
78(16.8%)
|
|
Yes
|
470(93.3%)
|
386(83.2%)
|
|
Decorating housingc
|
|
|
χ2 = 2.757; P = 0.097
|
No
|
462(91.7%)
|
438(94.4%)
|
|
Yes
|
42(8.3%)
|
26(5.6%)
|
|
aData presented as number (percentage) unless otherwise indicated.
bThe Wilcoxon rank-sum test method was used; otherwise, the χ2 test was used.
cThe exposure occurred in the periconceptional period.
dP < 0.05 was considered to indicate a statistically significant difference.
Maternal tobacco exposure and risk of CHD in offspring
Table 2 showed the association between maternal smoking and the risk of CHD in offspring. The prevalence rate of active smoking in 3 months before pregnancy in our controls (2.0%) was lower than the smoking rate among Chinese women in China Adult Tobacco Survey Report in 2015 (2.7%)37. None of the mothers in cases and controls reported active smoking in the first trimester of pregnancy. Mothers who reported active smoking in 3 months before pregnancy had an increased risk of CHD in offspring compared with the controls (P <0.001), but this association was not independent of potential confounders (P = 0.052). After adjustment for baseline data, mothers exposed to secondhand smoke in 3 months before pregnancy were observed an increased risk of CHD in offspring (aOR = 1.56; 95% CI: 1.13-2.15). Additionally, the risk of CHD in offspring was significantly higher among mothers who were exposed to secondhand smoke in the first trimester of pregnancy (aOR = 2.24; 95% CI: 1.57-3.20).
Table 2. Maternal smoking and risk of congenital heart defects in offspringa
Exposure
|
Control group
(n=504)
|
Case group
(n=464)
|
Unadjusted OR (95% CI)
|
P
|
Adjusted OR (95% CI)b
|
P c
|
Active smoking in 3 months before pregnancy
|
|
|
|
No
|
494 (98.0%)
|
432 (93.1%)
|
1.00 (reference)
|
-
|
1.00 (reference)
|
-
|
Yes
|
10 (2.0%)
|
32 (6.9%)
|
3.66 (1.78-7.53)
|
<0.001
|
2.37(0.99-5.65)
|
0.052
|
Active smoking in the first trimester
|
|
|
|
No
|
504 (100%)
|
464 (100%)
|
-
|
-
|
-
|
-
|
Yes
|
0 (0%)
|
0 (0%)
|
-
|
-
|
-
|
-
|
Passive smoking in 3 months before pregnancy
|
|
|
|
No
|
316 (62.7%)
|
222 (47.8%)
|
1.00 (reference)
|
-
|
1.00 (reference)
|
-
|
Yes
|
188 (37.3%)
|
242 (52.2%)
|
1.83 (1.42-2.37)
|
<0.001
|
1.56 (1.13-2.15)
|
0.007
|
Passive smoking in the first trimester
|
|
|
|
No
|
406 (80.6%)
|
274 (59.1%)
|
1.00 (reference)
|
-
|
1.00 (reference)
|
-
|
Yes
|
98 (19.4%)
|
190 (40.9%)
|
2.87 (2.154-3.83)
|
<0.001
|
2.24 (1.57-3.20)
|
<0.001
|
Abbreviations: CI = confidence interval.
aData presented as number (percentage) unless otherwise indicated.
bAdjusted for maternal education level (years), annual income in the past 1 year (RMB), history of adverse pregnancy outcomes, consanguineous marriage, history of congenital malformations in family, cold or fever in the periconceptional period, and personal lifestyle and habit in the periconceptional period including drinking alcohol, drinking tea, living near environmental pollution source, dyeing hair or perming and folate use.
cP < 0.05 was considered to indicate a statistically significant difference.
Genotypes Frequencies of SNPs and the results of HWE tests and GMS
The genotype frequencies for each SNP of the maternal MTHFD gene and the results of HWE tests were summarized in Table S2. The HWE tests showed that the genotype frequencies of the 5 SNPs of maternal MTHFD gene in the control group were all within HWE (all P values >0.01). The results of the GMS of each SNP were presented in Table S3. The genetic models of SNPs including rs1950902, rs2236225, rs2236222 and rs11849530 were all classified into the addictive model since ZHWDTT > -1.645 and ZHWDTT < 1.645. The genetic model of rs1256142 was classified into the dominant model because of ZHWDTT < -1.645. We initially ascertained the genetic models of overall SNPs of the maternal MTHFD gene, which was used for accessing the association between each SNP and risk of CHD in offspring based on the corresponding genetic model.
Genetic variants of maternal MTHFD gene and risk of CHD in offspring
The association between each maternal SNP of the MTHFD gene and the risk of CHD in the Han Chinese population was shown in Table 3. The univariate analyses suggested that there were statistically significant differences for the genetic variants at rs1950902 (AA vs. GG: P = 0.006; the addictive model: P = 0.002), rs2236222 (GG vs. AA: P = 0.001; the addictive model: P = 0.001) and rs1256142 (GA vs. GG: P = 0.035)
Table 3. MTHFD genes in mothers and risk of congenital heart disease in offspring
SNPs
|
Univariate logistic regression
|
|
Multivariable logistic regressionc
|
|
Unadjusted OR (95% CI)
|
P
|
|
Adjusted OR (95% CI)
|
P
|
FDR_Pd
|
rs1950902
|
|
|
|
|
|
|
G/G
|
1.00 (reference)
|
-
|
|
1.00 (reference)
|
-
|
-
|
G/A
|
1.29 (0.85-1.96)
|
0.224
|
|
1.38 (0.80-2.39)
|
0.247
|
0.309
|
A/A
|
1.80 (1.18-2.73)
|
0.006
|
|
1.73 (1.01-2.97)
|
0.046
|
0.090
|
Additivea
|
1.36 (1.12-1.64)
|
0.002
|
|
1.30 (1.02-1.65)
|
0.025
|
0.090
|
rs2236225
|
|
|
|
|
|
|
G/G
|
1.00 (reference)
|
-
|
|
1.00 (reference)
|
-
|
-
|
G/A
|
1.21 (0.92-1.59)
|
0.178
|
|
1.16 (0.81-1.65)
|
0.427
|
0.493
|
A/A
|
1.16 (0.61-2.20)
|
0.648
|
|
1.03 (0.47-2.27)
|
0.937
|
0.937
|
Additive
|
1.15 (0.92-1.44)
|
0.211
|
|
1.09 (0.82-1.45)
|
0.546
|
0.585
|
rs2236222
|
|
|
|
|
|
|
A/A
|
1.00 (reference)
|
-
|
|
1.00 (reference)
|
-
|
-
|
G/A
|
1.27 (0.97-1.66)
|
0.087
|
|
1.27 (0.96-1.67)
|
0.096
|
0.144
|
G/G
|
2.50 (1.46-4.29)
|
0.001
|
|
2.38 (1.38-4.12)
|
0.002
|
0.015
|
Additive
|
1.42 (1.16-1.75)
|
0.001
|
|
1.40 (1.14-1.73)
|
0.002
|
0.015
|
rs11849530
|
|
|
|
|
|
|
A/A
|
1.00 (reference)
|
-
|
|
1.00 (reference)
|
-
|
-
|
G/A
|
0.91 (0.69-1.20)
|
0.498
|
|
1.24 (0.87-1.77)
|
0.243
|
0.309
|
G/G
|
1.13 (0.77-1.65)
|
0.536
|
|
1.68 (1.02-2.77)
|
0.042
|
0.090
|
Additive
|
1.02 (0.86-1.22)
|
0.809
|
|
1.28 (1.02-1.62)
|
0.037
|
0.090
|
rs1256142
|
|
|
|
|
|
|
G/G
|
1.00 (reference)
|
-
|
|
1.00 (reference)
|
-
|
-
|
G/A
|
1.44 (1.03-2.02)
|
0.035
|
|
1.57 (1.01-2.45)
|
0.048
|
0.090
|
A/A
|
1.20 (0.83-1.74)
|
0.340
|
|
1.57 (0.97-2.56)
|
0.068
|
0.113
|
Dominantb
|
0.92 (0.70-1.21)
|
0.550
|
|
1.57 (1.03-2.40)
|
0.037
|
0.090
|
Abbreviations: CI = confidence interval; SNPs = single nucleotide polymorphisms; MTHFD = methylenetetrahydrofolate dehydrogenase; FDR_P = false discovery rate P value.
aAddictive means wild type homozygotes vs. heterozygote vs. mutant type homozygote.
bDominant means wild type homozygote vs. mutant type homozygotes and heterozygote.
cAdjusted for maternal education level (years), annual income in the past 1 year (RMB), history of adverse pregnancy outcomes, consanguineous marriage, history of congenital malformations in family, cold or fever in the periconceptional period, and personal lifestyle and habit in the periconceptional period including drinking alcohol, drinking tea, living near environmental pollution source, dyeing hair or perming and folate use.
dFDR_P < 0.1 was considered to indicate a statistically significant difference.
between the case and control groups. We further accessed the potential associations between genetic variants of maternal MTHFD gene and risk of CHD in offspring by aORs and their 95% CIs from logistic regression analysis. After adjustment for the potential confounders, the polymorphisms including rs1950902 (AA vs. GG: aOR = 1.73, 95% CI: 1.01-2.97; the addictive model: aOR =1.30, 95% CI: 1.02-1.65), rs2236222 (GG vs. AA: aOR = 2.38, 95% CI: 1.38-4.12; the addictive model: aOR = 1.40, 95% CI: 1.14-1.73) and rs1256142 (GA vs.GG: aOR = 1.57, 95% CI: 1.01-2.45; the dominant model: aOR = 1.57, 95% CI: 1.03-2.40) were still observed an increased risk for CHD, respectively. In addition, rs11849530 (GG vs. AA: aOR = 1.68, 95% CI: 1.02-2.77; the addictive model: aOR = 1.28, 95% CI: 1.02-1.62) was also observed a significant association with higher CHD risk.
Table 4. Interactions between SNPs of MTHFD gene and maternal smoking detected by logistic regression
SNPs
|
Active smoking before pregnancy
|
Passive smoking before pregnancy
|
Passive smoking in the first trimester
|
aOR (95%CI)a
|
FDR_Pb
|
aOR (95%CI)
|
FDR_P
|
aOR (95%CI)
|
FDR_P
|
rs1950902 (additive)
|
1.44 (1.02-2.04)
|
0.075
|
1.22(1.07-1.40)
|
0.008
|
1.38(1.18-1.60)
|
<0.001
|
rs2236225 (additive)
|
2.15 (1.14-4.05)
|
0.075
|
1.38(1.12-1.70)
|
0.008
|
1.62(1.29-2.04)
|
<0.001
|
rs2236222 (additive)
|
1.71 (0.87-3.36)
|
0.149
|
1.17(0.97-1.42)
|
0.129
|
1.38(1.11-1.72)
|
0.004
|
rs11849530 (additive)
|
2.10 (1.21-3.64)
|
0.075
|
1.35(1.13-1.61)
|
0.008
|
1.64(1.34-2.00)
|
<0.001
|
rs1256142 (dominant)
|
1.52 (0.96-2.41)
|
0.149
|
1.26(1.06-1.50)
|
0.023
|
1.53(1.26-1.86)
|
<0.001
|
Abbreviations: aOR = adjusted odds ratio; CI = confidence interval; SNPs = single nucleotide polymorphisms; MTHFD = methylenetetrahydrofolate dehydrogenase; FDR_P = false discovery rate P value.
aAdjusted for maternal education level (years), annual income in the past 1 year (RMB), history of adverse pregnancy outcomes, consanguineous marriage, history of congenital malformations in family, cold or fever in the periconceptional period, and personal lifestyle and habit in the periconceptional period including drinking alcohol, drinking tea, living near environmental pollution source, dyeing hair or perming and folate use.
bFDR_P < 0.1 was considered to indicate a statistically significant difference.
Interactions between maternal smoke exposure and MTHFD gene for risk of CHD
We modestly identified the four polymorphisms including rs1950902, rs2236222, rs1256142, and rs11849530 with significant main effects on CHD risk in the Han Chinese population. Moreover, though rs2236225 was not observed a significant main effect on CHD risk, this polymorphism was the most extensively studied one. Thus, we kept the 5 SNPs of the MTHFD gene including rs1950902, rs2236222, rs1256142, rs11849530 and rs2236225 for the interactions analysis. Interactions between maternal SNPs of MTHFD gene in the corresponding genetic model and maternal smoke exposure on CHD risk were summarized in Table 4. Our results showed there were statistically significant interaction effects between active smoking in 3 months before pregnancy and genetic variants of maternal MTHFD gene at rs1950902, rs2236225 and rs11849530. Specifically, mothers with GG/GA genotypes at rs1950902 (the addictive model: aOR = 1.44, 95% CI: 1.02-2.04), AA/GA genotypes at rs2236225 (the addictive model: aOR = 2.15, 95% CI: 1.14-4.05) and GG/GA genotypes at rs11849530 (the addictive model: aOR = 2.10, 95% CI: 1.21-3.64) generated a 1.44-fold, 2.10-fold and 2.15-fold increased CHD risk when they smoked in 3 months before pregnancy, respectively. Additionally, maternal passive smoking was also observed interaction effects with MTHFD gene at rs1950902, rs2236225, rs2236222, rs11849530 and rs1256142 in the Han Chinese population. To be specific, the mothers with GG/GA genotypes at rs1950902 (the addictive model: aOR = 1.22, 95% CI: 1.07-1.40), AA/GA genotypes at rs2236225 (the addictive model: aOR = 1.38, 95% CI: 1.12-1.70), GG/GA genotypes at rs11849530 (the addictive model: aOR = 1.35, 95% CI: 1.13-1.61) or AA genotype at rs1256142 (the dominant model: aOR = 1.26, 95% CI: 1.06-1.50) had significantly higher CHD risk in offspring when they were exposed to secondhand smoke in 3 months before pregnancy. Moreover, when mothers carried the AA/GA genotypes at rs1950902 (the addictive model: aOR = 1.38, 95% CI: 1.18-1.60), AA/GA genotypes at rs2236225 (the addictive model: aOR = 1.62, 95% CI: 1.29-2.04), GG/GA genotypes at rs2236222 (the addictive model: aOR = 1.38, 95% CI: 1.11-1.72), GG/GA genotypes at rs11849530 (the addictive model: aOR = 1.64, 95% CI: 1.34-2.00) or AA genotype at rs1256142 (the dominant model: aOR = 1.53, 95% CI: 1.26-1.86), exposure to secondhand smoke in the first trimester would increase the susceptibility to suffer from a CHD-affected delivery.