Skip to main content

Association of maternal phthalates exposure and metabolic gene polymorphisms with congenital heart diseases: a multicenter case-control study

Abstract

Background

The majority of congenital heart diseases (CHDs) are thought to result from the interactions of genetics and the environment factors. This study aimed to assess the association of maternal non-occupational phthalates exposure, metabolic gene polymorphisms and their interactions with risk of CHDs in offspring.

Methods

A multicenter case-control study of 245 mothers with CHDs infants and 268 control mothers of health infant was conducted from six hospitals. Maternal urinary concentrations of eight phthalate metabolites were measured by ultra-high performance liquid chromatography coupled with tandem mass spectrometry (UHPLC-MS/MS). Twenty single nucleotide polymorphisms (SNPs) in cytochrome P450 family 2 subfamily C member 9 (CYP2C9) and 19 (CYP2C19), uridine diphosphate (UDP) glucuronosyl transferase family 1 member A7 (UGT1A7), family 2 member B7 (UGT2B7) and B15(UGT2B15) genes were genotyped. The multivariate logistic regressions were used to estimate the association between maternal phthalates exposure or gene polymorphisms and risk of CHDs. Generalized multifactor dimensionality reduction (GMDR) was used to analyze the gene–gene and gene–phthalates exposure interactions.

Results

There was no significant difference in phthalate metabolites concentrations between the cases and controls. No significant positive associations were observed between maternal exposure to phthalates and CHDs. The SNPs of UGT1A7 gene at rs4124874 (under three models, log-additive: aOR = 1.74, 95% CI:1.28–2.37; dominant: aOR = 1.86, 95% CI:1.25–2.78; recessive: aOR = 2.50, 95% CI: 1.26–4.94) and rs887829 (under the recessive model: aOR = 13.66, 95% CI: 1.54–121) were significantly associated with an increased risk of CHDs. Furthermore, the associations between rs4124874 (under log-additive and dominant models) of UGT1A7 were statistically significant after the false discovery rate correction. No significant gene-gene or gene-phthalate metabolites interactions were observed.

Conclusions

The polymorphisms of maternal UGT1A7 gene at rs4124874 and rs887829 were significantly associated with an increased risk of CHDs. More large-scale studies or prospective study designs are needed to confirm or refute our findings in the future.

Peer Review reports

Introduction

Congenital heart diseases (CHDs), typically defined as structural and functional abnormalities of the heart and great vessels during embryonic development, are the most common type of birth defects. The reported birth prevalence of CHDs varies widely across countries and continents. A recent meta-analysis has estimate that the worldwide prevalence of CHDs in live births increased to 9.41/‰ in the period 2010–2017 [1]. In China, total CHDs birth prevalence increased continuously, from 0.201‰ in 1980–1984 to 4.905‰ in 2015–2019 [2]. CHDs are the leading cause of infant morbidity and mortality, and account for approximately 40% of prenatal deaths and 20% of mortality in the first year of life [3]. Therefore, CHDs place a heavy burden on the healthcare system, and have become a sizable public health concern. The etiology of CHDs is complex and multifactorial. To date, \(\sim\)20–30% of the CHDs cases could be identified with clear environmental or genetic factors, most CHDs cases are considered to be caused by the interaction of genetic and environmental factors [4,5,6].

Phthalates, are a family of synthetic chemicals widely used as plasticizers in many consumer products. Low-molecular-weight phthalates, such as di-n-butyl phthalate (DBP) and di-isobutyl phthalate (DiBP), are present in fragrances, pharmaceuticals, cosmetics, personal-care products, and packaging materials, whereas high-molecular-weight phthalates, such as di-2-ethylhexyl phthalate (DEHP) and benzyl-n-butyl ortho-phthalate (BBP), are used to soften plastics, particularly PVC building materials [7, 8]. Phthalates are not covalently bound to other substances and can easily permeate and migrate into environmental media, as result, the general population is ubiquitously exposed to phthalates via oral, dermal, inhalation, and intravenous routes [9, 10]. Phthalate metabolites have been consistently detected in urine from the general population including pregnant women worldwide [7,8,9, 11]. Between 2000 and 2017, indoor phthalate exposure in China has led to 3.32 million disability-adjusted life years (DALYs) per year, accounting for 0.90% of total DALYs across China, with DBP, DiBP and DEHP being the most abundant phthalates in indoor environments of residences, offices, and schools [12].

Phthalates can easily cross the placental barrier and impact the developing fetus in utero [7]. Growing literatures have reported that prenatal exposure to phthalates were associated with the risk of multiple adverse pregnancy outcomes, including preterm birth [13, 14], spontaneous pregnancy loss [15, 16], low birth weight [17], male newborn genital anomalies [18], and fetal retarded growth [19]. Two animal experiment have observed that BBP or DBP exposure could cause abnormalities in zebrafish embryo morphology, including cardiac structure deformities [20, 21]. To date, population epidemiological studies on the association between maternal occupational phthalates exposure and CHDs are limited, and the results are inconsistent. For example, two case-control studies have shown that maternal occupational phthalates exposures are associated with increased risk of total CHDs [22] or some CHDs subtypes [23]. However, three case-control studies have not found significant association between maternal occupational phthalates exposure to and CHDs [24,25,26]. Besides, phthalates exposure is assessed through subjective reporting, which may result in recall bias and not fully objectively reflect the exact level in vivo. To our knowledge, there have no studies on the association between maternal non-occupational phthalates exposure, identified by the biomarker phthalates, and the risk of CHDs.

The metabolism of environmental exposure is the first step to affect the occurrence of disease, and the individual metabolism level can significantly affect the risk of disease occurrence. The cytochrome P450 enzymes (CYPs), such as CYP2C9, and CYP2C19, and uridine diphosphate (UDP) -glucuronyl transferases (UGTs), such as UGT1A7, UGT2B7 and UGT2B15, have been proposed to play important roles in the phase I and phase II biotransformation of phthalates, respectively [7, 27, 28]. Common single nucleotide polymorphisms (SNPs) genetic polymorphisms in these genes could affect individual susceptibility to adverse effects of exposure to phthalates. One study suggested that the SNPs of UGT, such as rs7439366 of UGT2B7, rs1902023 of UGT2B15, were associated with the clearance of bisphenol A and phthalates in in patients with polycystic ovary syndrome [29]. One study observed that rs1799853 and rs1057910 of CYP2C9 could reduce DEHP biotransformation, and rs1799853 and rs1057910 of CYP2C9, rs12248560 of CYP2C19, and rs11692021of UGT1A7 might represent biomarkers of susceptibility or resilience in phthalates exposure [28]. However, few studies have investigated maternal genetic susceptibility to CHDs related to phthalates. Additionally, few studies have explored possible gene-environment interactions.

In the present study, we first evaluated the association between maternal exposure to phthalates by measuring urinary levels of phthalate metabolites during the second pregnancy and the risk of fetal CHDs. Then, we investigated the association between maternal genetic polymorphisms and the risk of fetal CHDs. Finally, we explored the potential interaction between maternal genetic variants and exposure to phthalates on the risk of CHDs.

Methods

Study population

This study was based on the project of a multicenter hospital-based case-control study which was performed from February 2010 to July 2015. The subjects were recruited from six tertiary maternal and child health hospitals with the qualification of performing prenatal diagnoses, located in Fuzhou, Nanning, Chengdu, Zhengzhou, Wuhan and Shenzhen city, respectively.

The cases group were pregnant women whose fetuses were diagnosed with CHDs by echocardiography and no any extracardiac abnormalities and with gestational age more than 12 weeks. CHDs were further confirmed by humanitarian examination of the pathological anatomy for aborted fetuses or ultrasound examination performed within 30 postnatal days and telephone follow-up performed within 60 days for liveborn fetuses. The exclusion criteria were as follows: (1) fetuses with syndromic diseases and chromosomal aberrations, (2) pregnant women with multiple pregnancies, (3) a mother with a family history of CHDs.

The control group were pregnant women with a singleton pregnancy without major congenital malformations diagnosed by echocardiography in the same hospital and with gestational age more than 12 weeks. A further ultrasound examination was performed within 30 postnatal days, and a telephone follow-up was performed within 60 days.

Each participant participated in face-to-face interview. The structured questionnaire was composed of parental social demographics, living environment and habits, working environment, maternal reproductive history, maternal illness and drug use history, maternal diet and nutrition, and maternal life events and mental state.

Four millilitres of blood was collected in EDTA from each participant during the second trimester by venepuncture and stored at − 70℃until genotyping. Ten millilitres single-spot urine samples of each participant during the second trimester was collected in a nonsterile clean polypropylene container, then divided into aliquots and stored at − 80 °C until further analysis.

Measurement of urinary concentrations of phthalate metabolites

Phthalate metabolites concentration was measured at the West China School of Public Health, Sichuan University by ultra-high performance liquid chromatography coupled with tandem mass spectrometry (UHPLC-MS/MS). Eight phthalate metabolites from four parent compounds were quantified, including mono-n-butyl phthalate (MnBP, metabolite of di-n-butyl phthalate (DnBP)); mono-isobutyl phthalate (MiBP, metabolite of diisobutyl phthalate (DiBP)); mono-benzyl phthalate (MBzP, metabolite of butylbenzyl phthalate (BBzP)); five metabolites of di (2-ethylhexyl) phthalate (DEHP):[mono(2-ethyl-hexyl) phthalate (MEHP),mono(2-ethyl-5-hydroxyhexyl) phthalate (MEHHP),mono(2-ethyl-5-oxohexyl) phthalate (MEOHP), mono (2-ethyl-5-carboxypentyl)phthalate (MECCP), and mono-2- carboxymethyl hexyl phthalate (MCMHP)].

The analysis was performed using an ultra-high performance liquid chromatography ACQUITY UPLC I-Class coupled to an Xevo TQ-XS triple stage quadrupole mass spectrometer (Waters, USA). The details regarding the preparation and analysis of the samples are available in Supplementary Appendix A.

To adjust the dilution of urine, creatinine-adjusted concentrations of urinary phthalate metabolites were calculated. The concentrations below the limit of detection (LOD) were replaced by LOD/\(\sqrt{2}\).

DNA extraction and genotyping

Genomic DNA was extracted with magnetic bead method (BioTeke, Wuxi, China) according to the recommended protocol. Twenty SNPs in the UGT1A7, UGT2B7, UGT2B15, CYP2C19, CYP2C9 were selected based on the following principal criteria: (1) an association with diseases in previous studies or the metabolic level of phthalates [28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44], (2) a minor allele frequency > 0.05 in Han Chinese. In total, 20 SNPs were selected. The genotypes of these SNPs were detected using multiple-polymerase chain reaction amplification and next generation sequencing (iGeneTech Bioscience Co., Ltd, Beijing, China). More detailed information about the studied genetic variants and genotyping is presented in Supplementary Appendix, Table S2.

For quality-control assessment, genotyping was repeated in 10% of samples, and the consistency rate was 100%.

Statistical analyses

The composition ratio of baseline characteristics between case and control groups was compared by χ2 test using Statistical Package for Social Sciences (SPSS) version 16.0 software (SPSS Inc., IBM, Chicago, USA).These characteristics included maternal age (at the time of the last menstrual period), maternal ethnicity, maternal education level, parental smoking or environmental tobacco smoke (ETS) exposure, maternal alcohol consumption, parity, maternal pre-pregnancy body mass index (ppBMI) (BMI = weight/(height × height) (kg/m2) before known pregnancy), maternal folic acid supplement.

As the distributions of phthalate metabolites concentrations did not meet the normality assumption, they were described as median (interquartile range) and compared with Mann–Whitney U test. The concentrations of phthalate metabolites were performed natural log10 transformation and further divided into three categories according to the tertile which they fell into. Logistic regression analysis was used to estimate the associations of the levels of phthalate metabolites with CHDs using SPSS version 16.0 software (SPSS Inc., IBM, Chicago, USA).

Hardy–Weinberg equilibrium was assessed in the controls using Plink software (http://zzz.bwh.harvard.edu/plink/ld.shtml). Logistic regression analysis was performed to investigate the association between individual genetic polymorphisms and CHDs using Plink software.

The effects of the gene–gene and gene-phthalates exposure interactions on CHDs occurrence were evaluated by logistic models using generalized multifactor dimensionality reduction (GMDR, version 0.7, University of Virginia, Charlottesville, VA).

All analyses were adjusted for covariates or potential confounders. False discovery rate (FDR) correction of multiple hypothesis testing was performed. Two-sided P < 0.05 was considered statistically significant.

Results

Characteristics of the study participants

In this study, a total of 513 subjects were analyzed, including 245 cases with CHDs fetuses and 268 controls. The baseline characteristics of participants were listed in Table 1.

Table 1 Descriptive characteristics of the participants

There were significant differences in the distributions of maternal ethnicity (96.64% Han in controls and 77.14% Han in cases, P < 0.001), maternal education level (70.52% controls vs. 37.55% cases had College or higher education, P < 0.001), parental smoking or ETS exposure (53.73% controls vs. 42.45% cases, P = 0.011).

Compared with control mothers, CHDs case mothers were also more likely to be multipara (46.64 controls vs.64.08% cases, P < 0.001). Maternal age, maternal alcohol consumption, pre-pregnancy BMI and folic acid supplements did not show statistical differences between the two groups (Pall>0.05).

Distribution of phthalate metabolite concentrations in maternal urine

The LODs, detection rates, and distributions of phthalate metabolites were summarized in Table 2. The LOD for MiBP, MnBP, MBzP, MEHP, MEHHP, MEOHP, MECPP and MCMHP was 0.50 ng/mL, 0.50 ng/mL, 0.05 ng/mL, 0.50 ng/mL, 0.05 ng/mL, 0.05 ng/mL, 0.10 ng/mL and 0.20 ng/mL, respectively. In addition to relatively low detection rate of 84.60% for MBzP, the detection rates of other seven metabolites were nearly or equal to 100%.

Table 2 Urinary concentrations of unadjusted (ng/mL) and adjusted PAEs (µg/g Cr)

The median concentrations of MnBP, MEHHP, MEOHP and MECPP were higher in CHDs than in controls (41.11 vs. 34.08, 2.51 vs. 2.35, 1.48 vs. 1.40, 5.49 vs. 5.43 µg/g creatinine, respectively), whereas the median concentrations of MiBP, MBzP, MEHP and MCMHP were lower in CHDs than in controls (10.01 vs. 10.43, 0.08 vs. 0.09, 2.37 vs. 2.71, 1.55 vs. 1.60 µg/g creatinine, respectively). However, these eight phthalate metabolites concentrations did not show significant differences between two groups (Pall>0.05).

Associations between maternal phthalates exposure and CHDs risk

The relation between maternal phthalates exposure and the risk of CHDs was displayed in Table 3. Among the eight metabolites, including MiBP, MnBP, MBzP, MEHP, MEHHP, MEOHP, MECPP, MCMHP, the first-tertile log10-transformed concentration of each was used as a reference, the second- and third-tertile concentrations were not associated with the risks of CHDs. Stratified analysis in the Han Chinese maternal population also did not observe significant positive associations between maternal exposure to phthalates and CHDs (Supplementary Appendix Table S3).

Table 3 Logistic regression analyses of the association between phthalate metabolites in maternal urinary samples and the risk of CHDs

Association between maternal gene polymorphisms and CHDs risk

The genotype distributions for polymorphisms of UGT1A7, UGT2B15, UGT2B7, CYP2C19 and CYP2C9 in the controls were consistent with Hardy-Weinberg equilibrium (see Supplementary Appendix, Table S4).

The association between single gene loci polymorphisms and the risk of CHDs when assuming various genetic models was shown in Table 4. In the UGT1A7 gene, the SNP rs4124874 was associated with an increased risk of CHDs (under the log-additive model: aOR = 1.74, 95% CI:1.28–2.37; under the dominant model: aOR = 1.86, 95% CI:1.25–2.78; under the recessive model: aOR = 2.50, 95% CI: 1.26–4.94), and the SNP rs887829 was associated with an increased risk of CHDs under the recessive model (aOR = 13.66, 95% CI: 1.54–121). However, for the recessive model, the associations were not statistically significant after the false discovery rate (FDR) correction. No significant association was found between any of the remaining 18 selected loci and the risk of CHDs.

Table 4 Association between maternal genotypes and the risk of CHDs

GMDR analyses for gene–gene and gene–environment interactions on CHDs

The gene-gene and gene-environment interaction model by GMDR was presented in Table 5. The P value was determined using the permutation test with 1000 replications. For gene-gene interaction, rs4124874 was a susceptibility locus for CHDs risk. Two-locus to five-locus interaction models were observed, but there were no statistical differences. In addition, for gene-phthalate metabolites interaction, five interaction combinations with no statistical significance were observed.

Table 5 GMDR analysis for gene–gene and gene–phthalates exposure interaction models in CHDs

Discussion

In this case-control study, we observed no difference in the concentration of maternal phthalate metabolites between the case and control groups, and did not find significant association between the concentrations of maternal phthalate metabolites and the risk of CHDs, while observed that the SNPs rs4124874 and rs887829 of UGT1A7 gene were associated with an increased risk of CHDs, but did not find significant gene-gene or gene-phthalate metabolite interactions on CHDs.

Our findings suggested that maternal phthalates exposure was not associated with the risk of CHDs, which was similar to the results of two Netherlands studies showing no association of maternal exposure to phthalates with risk of CHDs [24, 26]. Another case-control study in Hungary also reported that maternal job-exposure matrix (JEM)-assessed and self-reported exposures to phthalates were not associated with the risk of CHDs or subtypes [25]. Moreover, our analysis in the Han Chinese maternal population found no significant positive associations between phthalates exposure to and CHDs. The results of stratified analysis were inconsistent with one previous studies reported in China, which found that maternal occupational exposure to phthalates was associated with a higher incidence of total congenital heart defects, with aOR 1.6 (95% CI: 1.0-2.6) [22]. Similarly, a later Chinese study confirmed that maternal occupation exposure to phthalates was associated with CHDs subtypes, including ventricular septal defect (VSD), atrial septal defect (ASD), patent ductus arteriosus (PDA), pulmonary valve stenosis (PVS) [23]. However, due to the small sample size, we did not perform the association analysis between phthalates exposure with CHDs subtypes. In addition to CHDs, maternal occupational exposure to phthalate obtained by a job exposure matrix was associated with increased risks of hypospadias (OR = 3.12; 95% CI, 1.04–11.46) [45]. First trimester urinary DEHP metabolite concentrations were associated with increased odds of neonatal genital anomaly (OR = 2.54, 95% CI, 1.09–5.92) [18]. One study found that the detection ratio of positive BBP and its metabolites in maternal urine was obviously higher in neural tube defects population than that in normal controls [46]. Overall, there are only a few limited reports on phthalate exposure during pregnancy and birth defects.

It is important to note that studies on the association between maternal phthalates exposure and the risk of CHDs in the Chinese population have shown inconsistent results, which may be due to differences in the way exposure was evaluated. Phthalates exposure is ubiquitous. Two previous China studies focused occupational exposures, which only considered a single source of exposure at work, ignoring comprehensive exposure in the living environment. Moreover, it was difficult to estimate accurately the amount of exposure, the time and the frequency due to recall bias. In addition, the discrepancies across studies may be due to the heterogeneity of populations. For example, two European studies found that maternal occupational phthalates exposure was not associated with CHDs, however, the inverse associations were observed between maternal occupational phthalates exposure and the risk of CHDs in China studies. The conflicting results can likely be explained by differences in sensitivity to the biochemical and toxic effects of phthalates due to genetic polymorphisms.

Accumulating evidence have suggested that inter-individual differences in the ability of the xenobiotic metabolism due to high variability in certain metabolic enzyme activities can influence the effects of environmental exposure on birth defects (e.g., oral clefts, neural tube defects, CHDs) [47,48,49,50]. The uridine diphosphate (UDP)-glucuronosyl transferase (UGTs), UGT1A7, UGT2B7 and UGT2B15, catalyze the glucuronidation of multiple substrates. As for the UGT1A7 gene, rs11692021 polymorphism was related to higher risk of chronic pancreatitis (aOR = 1.76, 95% CI: 1.26–2.46) [30]. In addition, rs4124874 and rs4148323 polymorphism were associated with in increased risk of hyperbilirubinemia [31, 32]. However, rs887829 genotype was related to a decreased risk of hyperbilirubinemia (aOR = 0.55, 95% CI: 0.34–0.89) [32]. As for the UGT2B7 gene, rs7662029 polymorphism was associated with the withdrawal symptoms in methadone maintenance patients [33]. SNP rs4587017 polymorphism might influence the analgesic effects of fentanyl in the cold pressor-induced pain test [34]. SNP rs12233719 polymorphism was associated with never-smoking female lung cancer risk [35]. As for the UGT2B15 gene, rs2045100 locus, six SNPs (rs3100, rs4148269, rs9994887, rs13112099, rs7686914 and rs7696472), and rs1902023 polymorphism, showed significant associations with increased risk for prostate cancer [36,37,38]. In the present study, among 17 SNPs in UGT1A7, UGT2B7 and UGT2B15, the polymorphisms of maternal UGT1A7 gene at rs4124874 (under additive, dominant and recessive models), and rs887829 (under recessive model) were associated with increased CHDs risk. The molecular mechanism remains unclear. It is possible that these two intronic loci affect the alternative spicing of the gene products, or they might be in linkage disequilibrium with other causal loci or genes, thereby affecting the metabolism of phthalates.

Two members of the CYP family, CYP2C9 and CYP2C19 are the main Phase I metabolizing enzymes mediating the toxicity of phthalates, their polymorphisms are associated with the risk of many types of diseases. It has been reported that rs1057910 genotype of CYP2C9 was related to increased risk of adenoma recurrence (aRR = 1.47, 95% CI: 1.19–1.83) [39], or developing sporadic colorectal carcinoma (aOR = 2.77, 95% CI: 1.1653–4.643) [40]. In addition, for the two common polymorphic loci rs12248560 and rs4244285 of CYP2C19, the genotype of rs12248560 was associated with decreased breast cancer risk (aOR = 0.77, 95% CI: 0.65–0.93) [41], while rs4244285 polymorphism was related to higher risk of epilepsy (aOR = 4.24, 95% CI: 2.52–7.15) [42], long-term ischemic stroke events (hazard ratio: 1.64, 95% CI: 1.06–2.53) [43], hypertension (aOR = 2.433, 95% CI: 1.797–3.293) [44]. Meanwhile, several studies have observed that CYP2C9 rs1057910 polymorphism was associated with increased risk for adenoma recurrence (aRR = 1.47, 95% CI 1.19–1.83) [39], or developing sporadic colorectal carcinoma (aOR = 2.589, 95% CI: 1.549–4.330) [40]. In our study, no significant associations between CHDs risk and genotype were seen for the rs12248560 and rs4244285 of CYP2C19, or rs1057910 of CYP2C9 polymorphisms.

It is now widely believed that most structural birth defects including CHDs are caused by a complex combination of genetic and environmental factors that interact to interfere with morphogenetic processes. More and more studies have reported significant interaction effects between gene-environment interactions for the development of CHDs. One study found maternal dietary factors and cystathionine beta synthase (CBS) gene variants (rs2851391, rs234714) interactions were significantly associated with risk of CHDs [51]. One study observed the interaction between maternal tobacco exposure and polymorphisms of the MTHFD1 gene including rs1950902, rs2236222, rs1256142, rs11849530 and rs2236225, was significantly associated with the risk of CHDs in offspring [52]. One study reported a significantly positive interaction between maternal folic acid supplementation and genetic variation at rs828858 of MTHFD2 for the risk of CHDs [53]. In our previous study, we found that polymorphisms of maternal GST genes (GSTM1, GSTT1, GSTP1) might modify the association of maternal smoke exposure with CHDs [54]. In addition, we also observed that the polymorphisms of maternal AHR rs2158041 and rs7811989, or UGT1A1 rs4148323 might modify the association of PAHs exposure with CHDs, CYP1A2 rs4646425 or CYP2E1 rs915908 polymorphisms significantly interacted with PAHs exposure on CHDs [55, 56]. In the present study, no significantly positive interaction of gene-gene or gene-phthalate metabolites for the risk of CHDs was observed, more large-scale studies or prospective study designs are needed to explore the interactions in the future.

This study has several strengths. First, to the best of our knowledge, this is the first study to evaluate the effect of the interaction between maternal phthalates exposure and maternal gene polymorphism on the risk of CHDs. Second, compared with previous phthalates exposure assessment that relied solely on expert industrial hygienist consensus or the self-reported questionnaire, we used a urinary bio-markers based approach to evaluate the association between maternal phthalates exposure and the risk of CHDs, offering an objective measure of exposure. Finally, the subjects of our study were non-occupational, low-dose phthalates exposed pregnant women, thus, the results can be generalized to all women because the environmental factors can be assessed at the individual level.

However, our study still had several limitations. First, small sample size limited the statistical power; besides, due to a small number of cases, our study could not perform the analysis for type-specific CHDs; future studies with larger sample sizes are warranted to confirm or refute our findings. Second, we only measured phthalate metabolites from a single spot urine sample taken at the second trimester which could not precisely estimate the mother’s long-term exposure level. Thus, future studies are needed to collect multiple urine samples. Third, only maternal phthalates exposure and genetic susceptibilities were considered; future studies are needed to investigate the effects of fetal exposure, fetal genotypes, and the interaction between them on the risk of CHDs.

In conclusion, our analysis results indicated that maternal phthalates exposure was not associated with the risk of CHDs. The polymorphisms of maternal UGT1A7 gene at rs4124874 and rs887829 were significantly associated with an increased risk of CHDs. No significant gene-gene or gene-phthalate metabolites interactions on CHDs was observed. However, due to the complex pathogenesis of CHDs and the limitation of small sample size, more large-scale studies or prospective study designs are needed to confirm or refute our findings in the future.

Data availability

The variant data for this study have been deposited in the European Variation Archive (EVA) at EMBL-EBI under accession number PRJEB70402.

Abbreviations

CHDs:

congenital heart diseases

aOR:

adjusted odds ratio

95% CI:

95% confidence interval

UHPLC-MS/MS:

ultra-high performance liquid chromatography coupled with tandem mass spectrometry

SNPs:

single nucleotide polymorphisms

CYP2C9:

cytochrome P450 family 2 subfamily C member 9

CYP2C9:

cytochrome P450 family 2 subfamily C member 19

UGT1A7:

uridine diphosphate (UDP) glucuronosyl transferase family 1 member A7

UGT2B7:

uridine diphosphate (UDP) glucuronosyl transferase family 2 member B7

UGT2B15:

uridine diphosphate (UDP) glucuronosyl transferase family 2 member B15

GMDR:

generalized multifactor dimensionality reduction

FDR:

false discovery rate

MnBP:

mono-n-butyl phthalate

DnBP:

di-n-butyl phthalate

MiBP:

mono-isobutyl phthalate

DiBP:

diisobutyl phthalate

MBzP:

mono-benzyl phthalate

BBzP:

butylbenzyl phthalate

DEHP:

di (2-ethylhexyl) phthalate

MEHP:

mono(2-ethyl-hexyl) phthalate

MEHHP:

mono(2-ethyl-5-hydroxyhexyl) phthalate

MEOHP:

mono(2-ethyl-5-oxohexyl) phthalate

MECCP:

mono (2-ethyl-5-carboxypentyl)phthalate

MCMHP:

mono-2- carboxymethyl hexyl phthalate

LOD:

limit of detection

ETS:

environmental tobacco smoke

References

  1. Liu Y, Chen S, Zuhlke L, Black GC, Choy MK, Li N, Keavney BD. Global birth prevalence of congenital heart defects 1970–2017: updated systematic review and meta-analysis of 260 studies. Int J Epidemiol. 2019;48(2):455–63.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Zhao L, Chen L, Yang T, Wang T, Zhang S, Ye Z, Luo L, Qin J. Birth prevalence of congenital heart disease in China, 1980–2019: a systematic review and meta-analysis of 617 studies. Eur J Epidemiol. 2020;35(7):631–42.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Nagy O, Barath S, Ujfalusi A. The role of microRNAs in congenital heart disease. EJIFCC. 2019;30(2):165–78.

    CAS  PubMed  PubMed Central  Google Scholar 

  4. Pierpont ME, Brueckner M, Chung WK, Garg V, Lacro RV, McGuire AL, Mital S, Priest JR, Pu WT, Roberts A, et al. Genetic basis for congenital Heart Disease: Revisited: A Scientific Statement from the American Heart Association. Circulation. 2018;138(21):e653–e711.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Krauss RS, Hong M. Gene-environment interactions and the etiology of birth defects. Curr Top Dev Biol. 2016;116:569–80.

    Article  CAS  PubMed  Google Scholar 

  6. Moreau JLM, Kesteven S, Martin E, Lau KS, Yam MX, O’Reilly VC, Del Monte-Nieto G, Baldini A, Feneley MP, Moon AM et al. Gene-environment interaction impacts on heart development and embryo survival. Development 2019, 146(4).

  7. Benjamin S, Masai E, Kamimura N, Takahashi K, Anderson RC, Faisal PA. Phthalates impact human health: epidemiological evidences and plausible mechanism of action. J Hazard Mater. 2017;340:360–83.

    Article  CAS  PubMed  Google Scholar 

  8. Wang Y, Zhu H, Kannan K. A review of Biomonitoring of Phthalate exposures. Toxics 2019, 7(2).

  9. Johns LE, Cooper GS, Galizia A, Meeker JD. Exposure assessment issues in epidemiology studies of phthalates. Environ Int. 2015;85:27–39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Wang Z, Ma J, Wang T, Qin C, Hu X, Mosa A, Ling W. Environmental health risks induced by interaction between phthalic acid esters (PAEs) and biological macromolecules: a review. Chemosphere. 2023;328:138578.

    Article  CAS  PubMed  Google Scholar 

  11. Gao H, Zhu YD, Xu YY, Zhang YW, Yao HY, Sheng J, Jin ZX, Ren LL, Huang K, Hao JH, et al. Season-dependent concentrations of urinary phthalate metabolites among Chinese pregnant women: repeated measures analysis. Environ Int. 2017;104:110–7.

    Article  CAS  PubMed  Google Scholar 

  12. Liu W, Sun Y, Liu N, Hou J, Huo X, Zhao Y, Zhang Y, Deng F, Kan H, Zhao Z, et al. Indoor exposure to phthalates and its burden of disease in China. Indoor Air. 2022;32(4):e13030.

    Article  CAS  PubMed  Google Scholar 

  13. Gao H, Wang YF, Huang K, Han Y, Zhu YD, Zhang QF, Xiang HY, Qi J, Feng LL, Zhu P, et al. Prenatal phthalate exposure in relation to gestational age and preterm birth in a prospective cohort study. Environ Res. 2019;176:108530.

    Article  CAS  PubMed  Google Scholar 

  14. Zhang Y, Mustieles V, Yland J, Braun JM, Williams PL, Attaman JA, Ford JB, Calafat AM, Hauser R, Messerlian C. Association of Parental Preconception Exposure to phthalates and phthalate substitutes with Preterm Birth. JAMA Netw Open. 2020;3(4):e202159.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Zhang H, Gao F, Ben Y, Su Y. Association between phthalate exposure and risk of spontaneous pregnancy loss: a systematic review and meta-analysis. Environ Pollut. 2020;267:115446.

    Article  CAS  PubMed  Google Scholar 

  16. Ji H, Wu Z, Chen D, Miao M, Chen H, Shuai W, Liang H, Yuan W. Individual and joint effects of phthalates exposure on the risk of early miscarriage. J Expo Sci Environ Epidemiol 2023.

  17. Jin S, Cui S, Xu J, Zhang X. Associations between prenatal exposure to phthalates and birth weight: a meta-analysis study. Ecotoxicol Environ Saf. 2023;262:115207.

    Article  CAS  PubMed  Google Scholar 

  18. Sathyanarayana S, Grady R, Barrett ES, Redmon B, Nguyen RHN, Barthold JS, Bush NR, Swan SH. First trimester phthalate exposure and male newborn genital anomalies. Environ Res. 2016;151:777–82.

    Article  CAS  PubMed  Google Scholar 

  19. Casas M, Valvi D, Ballesteros-Gomez A, Gascon M, Fernandez MF, Garcia-Esteban R, Iniguez C, Martinez D, Murcia M, Monfort N, et al. Exposure to Bisphenol A and phthalates during pregnancy and ultrasound measures of fetal growth in the INMA-Sabadell cohort. Environ Health Perspect. 2016;124(4):521–8.

    Article  CAS  PubMed  Google Scholar 

  20. Sun G, Liu K. Developmental toxicity and cardiac effects of butyl benzyl phthalate in zebrafish embryos. Aquat Toxicol. 2017;192:165–70.

    Article  CAS  PubMed  Google Scholar 

  21. Sun G, Li Y. Exposure to DBP induces the toxicity in early development and adverse effects on cardiac development in zebrafish (Danio rerio). Chemosphere. 2019;218:76–82.

    Article  ADS  CAS  PubMed  Google Scholar 

  22. Wang C, Xie L, Zhou K, Zhan Y, Li Y, Li H, Qiao L, Wang F, Hua Y. Increased risk for congenital heart defects in children carrying the ABCB1 gene C3435T polymorphism and maternal periconceptional toxicants exposure. PLoS ONE. 2013;8(7):e68807.

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  23. Wang C, Zhan Y, Wang F, Li H, Xie L, Liu B, Li Y, Mu D, Zheng H, Zhou K, et al. Parental occupational exposures to endocrine disruptors and the risk of simple isolated congenital heart defects. Pediatr Cardiol. 2015;36(5):1024–37.

    Article  CAS  PubMed  Google Scholar 

  24. Snijder CA, Vlot IJ, Burdorf A, Obermann-Borst SA, Helbing WA, Wildhagen MF, Steegers EA, Steegers-Theunissen RP. Congenital heart defects and parental occupational exposure to chemicals. Hum Reprod. 2012;27(5):1510–7.

    Article  CAS  PubMed  Google Scholar 

  25. Fazekas-Pongor V, Fekete M, Csaky-Szunyogh M, Cseh K, Penzes M. Parental occupational exposure and congenital heart diseases in a Hungarian case-control study. Int Arch Occup Environ Health 2021.

  26. Wijnands KP, Zeilmaker GA, Meijer WM, Helbing WA, Steegers-Theunissen RP. Periconceptional parental conditions and perimembranous ventricular septal defects in the offspring. Birth Defects Res Clin Mol Teratol. 2014;100(12):944–50.

    Article  CAS  Google Scholar 

  27. Dominguez-Romero E, Scheringer M. A review of phthalate pharmacokinetics in human and rat: what factors drive phthalate distribution and partitioning? Drug Metab Rev. 2019;51(3):314–29.

    Article  CAS  PubMed  Google Scholar 

  28. Stajnko A, Runkel AA, Kosjek T, Snoj Tratnik J, Mazej D, Falnoga I, Horvat M. Assessment of susceptibility to phthalate and DINCH exposure through CYP and UGT single nucleotide polymorphisms. Environ Int. 2022;159:107046.

    Article  CAS  PubMed  Google Scholar 

  29. Luo Y, Nie Y, Tang L, Xu CC, Xu L. The correlation between UDP-glucuronosyltransferase polymorphisms and environmental endocrine disruptors levels in polycystic ovary syndrome patients. Med (Baltim). 2020;99(11):e19444.

    Article  CAS  Google Scholar 

  30. Ockenga J, Vogel A, Teich N, Keim V, Manns MP, Strassburg CP. UDP glucuronosyltransferase (UGT1A7) gene polymorphisms increase the risk of chronic pancreatitis and pancreatic cancer. Gastroenterology. 2003;124(7):1802–8.

    Article  CAS  PubMed  Google Scholar 

  31. Li Z, Song L, Hao L. The role of UGT1A1 (c.-3279 T > G) gene polymorphisms in neonatal hyperbilirubinemia susceptibility. BMC Med Genet. 2020;21(1):218.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Zhou Y, Wang SN, Li H, Zha W, Wang X, Liu Y, Sun J, Peng Q, Li S, Chen Y, et al. Association of UGT1A1 variants and hyperbilirubinemia in breast-fed full-term Chinese infants. PLoS ONE. 2014;9(8):e104251.

    Article  ADS  PubMed  PubMed Central  Google Scholar 

  33. Tian JN, Ho IK, Tsou HH, Fang CP, Hsiao CF, Chen CH, Tan HK, Lin L, Wu CS, Su LW, et al. UGT2B7 genetic polymorphisms are associated with the withdrawal symptoms in methadone maintenance patients. Pharmacogenomics. 2012;13(8):879–88.

    Article  CAS  PubMed  Google Scholar 

  34. Muraoka W, Nishizawa D, Fukuda K, Kasai S, Hasegawa J, Wajima K, Nakagawa T, Ikeda K. Association between UGT2B7 gene polymorphisms and fentanyl sensitivity in patients undergoing painful orthognathic surgery. Mol Pain. 2016;12:1744806916683182.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Qian Y, Xie L, Li L, Feng T, Zhu T, Wang R, Yang Y, Zhou B, Yu H, Qian B. Association between sex hormones regulation-related SNP rs12233719 and lung cancer risk among never-smoking Chinese women. Cancer Med. 2021;10(5):1880–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Sun T, Oh WK, Jacobus S, Regan M, Pomerantz M, Freedman ML, Lee GS, Kantoff PW. The impact of common genetic variations in genes of the sex hormone metabolic pathways on steroid hormone levels and prostate cancer aggressiveness. Cancer Prev Res (Phila). 2011;4(12):2044–50.

    Article  CAS  PubMed  Google Scholar 

  37. Vidal AC, Tucker C, Schildkraut JM, Richardson RM, McPhail M, Freedland SJ, Hoyo C, Grant DJ. Novel associations of UDP-glucuronosyltransferase 2B gene variants with prostate cancer risk in a multiethnic study. BMC Cancer. 2013;13:556.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Grant DJ, Hoyo C, Oliver SD, Gerber L, Shuler K, Calloway E, Gaines AR, McPhail M, Livingston JN, Richardson RM, et al. Association of uridine diphosphate-glucuronosyltransferase 2B gene variants with serum glucuronide levels and prostate cancer risk. Genet Test Mol Biomarkers. 2013;17(1):3–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Barry EL, Poole EM, Baron JA, Makar KW, Mott LA, Sandler RS, Ahnen DJ, Bresalier RS, McKeown-Eyssen GE, Ulrich CM. CYP2C9 variants increase risk of colorectal adenoma recurrence and modify associations with smoking but not aspirin treatment. Cancer Causes Control. 2013;24(1):47–54.

    Article  PubMed  Google Scholar 

  40. Liao LH, Zhang H, Lai MP, Lau KW, Lai AK, Zhang JH, Wang Q, Wei W, Chai JH, Lung ML, et al. The association of CYP2C9 gene polymorphisms with colorectal carcinoma in Han Chinese. Clin Chim Acta. 2007;380(1–2):191–6.

    Article  CAS  PubMed  Google Scholar 

  41. Justenhoven C, Hamann U, Pierl CB, Baisch C, Harth V, Rabstein S, Spickenheuer A, Pesch B, Bruning T, Winter S, et al. CYP2C19*17 is associated with decreased breast cancer risk. Breast Cancer Res Treat. 2009;115(2):391–6.

    Article  CAS  PubMed  Google Scholar 

  42. Makowska M, Smolarz B, Brys M, Forma E, Romanowicz H. An association between the rs1799853 and rs1057910 polymorphisms of CYP2C9, the rs4244285 polymorphism of CYP2C19 and the prevalence rates of drug-resistant epilepsy in children. Int J Neurosci. 2021;131(12):1147–54.

    Article  CAS  PubMed  Google Scholar 

  43. Wu P, Liu Z, Tian Z, Wu B, Shao J, Li Q, Geng Z, Pan Y, Lu K, Wang Q et al. CYP2C19 loss-of-function variants Associated with Long-Term ischemic stroke events during Clopidogrel Treatment in the Chinese Population. Clin Pharmacol Ther 2023.

  44. Cai N, Li C, Gu X, Zeng W, Zhong J, Liu J, Zeng G, Zhu J, Hong H. CYP2C19 loss-of-function is associated with increased risk of hypertension in a Hakka population: a case-control study. BMC Cardiovasc Disord. 2023;23(1):185.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Ormond G, Nieuwenhuijsen MJ, Nelson P, Toledano MB, Iszatt N, Geneletti S, Elliott P. Endocrine disruptors in the workplace, hair spray, folate supplementation, and risk of hypospadias: case-control study. Environ Health Perspect. 2009;117(2):303–7.

    Article  PubMed  Google Scholar 

  46. Song G, Wang R, Cui Y, Hao CJ, Xia HF, Ma X. Oxidative stress response associates with the teratogenic effects of benzyl butyl phthalate (BBP). Toxicol Res (Camb). 2020;9(3):222–9.

    Article  PubMed  Google Scholar 

  47. Shi M, Christensen K, Weinberg CR, Romitti P, Bathum L, Lozada A, Morris RW, Lovett M, Murray JC. Orofacial cleft risk is increased with maternal smoking and specific detoxification-gene variants. Am J Hum Genet. 2007;80(1):76–90.

    Article  CAS  PubMed  Google Scholar 

  48. Wang L, Jin L, Liu J, Zhang Y, Yuan Y, Yi D, Ren A. Maternal genetic polymorphisms of phase II metabolic enzymes and the risk of fetal neural tube defects. Birth Defects Res Clin Mol Teratol. 2014;100(1):13–21.

    Article  CAS  Google Scholar 

  49. Padula AM, Yang W, Schultz K, Lurmann F, Hammond SK, Shaw GM. Genetic variation in biotransformation enzymes, air pollution exposures, and risk of spina bifida. Am J Med Genet A. 2018;176(5):1055–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Liu Z, Wang M, Yu P, Li X, Lin Y, Duan Y, Tian Y, Zhu J, Deng Y, Li N. Maternal trichloroethylene exposure and metabolic gene polymorphisms may interact during fetal cardiovascular malformation. Reprod Toxicol. 2021;106:1–8.

    Article  CAS  PubMed  Google Scholar 

  51. Li Y, Diao J, Li J, Luo L, Zhao L, Zhang S, Wang T, Chen L, Yang T, Zhu P, et al. Association of maternal dietary intakes and CBS gene polymorphisms with congenital heart disease in offspring. Int J Cardiol. 2021;322:121–8.

    Article  PubMed  Google Scholar 

  52. Song X, Li Q, Diao J, Li J, Li Y, Zhang S, Zhao L, Chen L, Wei J, Shu J, et al. Association of MTHFD1 gene polymorphisms and maternal smoking with risk of congenital heart disease: a hospital-based case-control study. BMC Pregnancy Childbirth. 2022;22(1):88.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Liu H, Ou J, Chen Y, Chen Q, Luo M, Wang T, Qin J. Association of Maternal Folate Intake and offspring MTHFD1 and MTHFD2 genes with congenital heart disease. Nutrients 2023, 15(16).

  54. Li X, Liu Z, Deng Y, Li S, Mu D, Tian X, Lin Y, Yang J, Li J, Li N, et al. Modification of the association between maternal smoke exposure and congenital heart defects by polymorphisms in glutathione S-transferase genes. Sci Rep. 2015;5:14915.

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  55. Li N, Mu Y, Liu Z, Deng Y, Guo Y, Zhang X, Li X, Yu P, Wang Y, Zhu J. Assessment of interaction between maternal polycyclic aromatic hydrocarbons exposure and genetic polymorphisms on the risk of congenital heart diseases. Sci Rep. 2018;8(1):3075.

    Article  ADS  PubMed  PubMed Central  Google Scholar 

  56. Tao J, Li N, Liu Z, Deng Y, Li X, Luo F, Yu P, Zhu J. Polymorphisms in gene UGT1A1 modify the association of prenatal exposure to polycyclic aromatic hydrocarbons with congenital heart diseases risk. J Matern Fetal Neonatal Med. 2023;36(1):2183743.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

We are indebted to the pediatric cardiologists, geneticists, and epidemiologists who collaborated in this program and made the study possible. We thank the obstetricians, pediatricians, pathologists, experimental technician and other participants involved in the project for recruiting the participants and collecting the data. We thank all participating families for their cooperation and for providing personal information. We also thank the reviewers for their helpful comments.

Funding

This work was supported by the National Key Research and Development Program of China (grant numbers 2022YFC2703300, 2022YFC2703302, 2022YFC2703700, 2022YFC2703702); Sichuan Science and Technology Program (grant numbers 2022NSFSC0661, 2022NSFSC0656).

Author information

Authors and Affiliations

Authors

Contributions

P. Y. and N. L. developed the study design and drafted the manuscript. H. K. and J. T. assisted in conducting the experiment and analyzing the data. L. L., M. W., Y. L. and W. X. assisted in preparing samples and extracting DNA. L. Z., D. Y., X. L., Y. W. and J. Z. participated in reviewing, editing, and revising the manuscript. All authors have read and approved the final manuscript.

Corresponding authors

Correspondence to Jing Tao or Ping Yu.

Ethics declarations

Ethics approval and consent to participate

All the participants signed an informed consent form. This research was approved by the Ethics Committee of Sichuan University (No. 2010004) and followed the tenets of the Declaration of Helsinki.

Consent for publication

Not Applicable.

Competing interests

The authors declare no competing interests.

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 4.0 International License, which permits use, sharing, adaptation, 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 changes were made. 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/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, N., Kang, H., Liu, Z. et al. Association of maternal phthalates exposure and metabolic gene polymorphisms with congenital heart diseases: a multicenter case-control study. BMC Pregnancy Childbirth 24, 167 (2024). https://doi.org/10.1186/s12884-024-06343-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12884-024-06343-z

Keywords