Association between Maternal Drug Use and Cytochrome P450 Genetic Polymorphisms and the Risk of Congenital Heart Defects in Offspring＊

QIN Jia Bi ,LUO Liu ,SUN Meng Ting ,HUANG Peng ,WANG Ting Ting ,ZHANG Sen Mao ,LI Jin Qi,LI Yi Huan,CHEN Le Tao,DIAO Jing Yi,#,and ZHU Ping

1.Guangdong Cardiovascular Institute,Guangdong Provincial People’s Hospital,Guangdong Academy of Medical Sciences,Guangzhou 510100,Guangdong,China;2.Department of Epidemiology and Health Statistics,Xiangya School of Public Health,Central South University,Changsha 410078,Hunan,China;3.NHC Key Laboratory of Birth Defect for Research and Prevention,Hunan Provincial Maternal and Child Health Care Hospital,Changsha 410008,Hunan,China;4.Hunan Children's Hospital,Changsha 410007,Hunan,China

Abstract Objective This study aimed to assess the associations between maternal drug use,cytochrome P450(CYP450) genetic polymorphisms,and their interactions with the risk of congenital heart defects (CHDs)in offspring.Methods A case-control study involving 569 mothers of CHD cases and 652 controls was conducted from November 2017 to January 2020.Results After adjusting for potential confounding factors,the results show that mothers who used ovulatory drugs (adjusted odds ratio [aOR]=2.12;95% confidence interval [CI]:1.08–4.16),antidepressants (aOR=2.56;95% CI:1.36–4.82),antiabortifacients (aOR=1.55;95% CI:1.00–2.40),or traditional Chinese drugs (aOR=1.97;95% CI:1.26–3.09) during pregnancy were at a significantly higher risk of CHDs in offspring.Maternal CYP450 genetic polymorphisms at rs1065852 (A/T vs.A/A:OR=1.53,95% CI:1.10–2.14;T/T vs. A/A:OR=1.57,95% CI:1.07–2.31) and rs16947 (G/G vs. C/C:OR=3.41,95%CI:1.82–6.39) were also significantly associated with the risk of CHDs in offspring.Additionally,significant interactions were observed between the CYP450 genetic variants and drug use on the development of CHDs.Conclusions In those of Chinese descent,ovulatory drugs,antidepressants,antiabortifacients,and traditional Chinese medicines may be associated with the risk of CHDs in offspring.Maternal CYP450 genes may regulate the effects of maternal drug exposure on fetal heart development.

Key words:Congenital heart defect;Maternal drug use;Cytochrome P450 genes;Case-control study

INTRODUCTION

Congenital heart defects (CHDs) are the most common type of congenital malformation,representing a major global health problem.The global prevalence of CHDs is approximately 8.22 per 1,000 live births[1].The highest prevalence has been reported in Asia(9.3‰),followed by Europe (8.2‰),and the lowest occurs in Africa (1.9‰)[2].Presently,the prevalence of CHD in China is 8.98‰[3].As CHDs are the main cause of non-infectious death for infants aged＜ 5-year-old,all cardiac defects were responsible for 46% of infant deaths caused by congenital malformations[4].Additionally,children with a CHD are at significantly higher risk of cardiovascular diseases compared with healthy controls,and they achieve significantly lower height,weight,and linguistic competence in later life[5].Thus,further etiological studies are needed for more targeted planning and decision-making in controlling and managing these diseases.

It is generally accepted that most CHD cases result from a combination of environmental and genetic factors,indicating that the complex interactions between environmental teratogens and genetic factors should be the direction of etiological research on CHDs.Many environmental factors have been identified as risk factors of CHDs[6-8],but the association between maternal drug use and the risk of CHDs remains controversial.It is difficult for pregnant women to completely avoid exposure to some drugs such as oral contraceptives,ovulatory drugs,antibiotics,antidepressants,and antiabortifacients,indicating that more attention should be paid to the potential effects of these drugs on embryonic heart development.Particularly in China,few studies have considered the relationship between maternal drug intake and CHDs.

The effects of environmental factors may be modified by the genes that are responsible for activating and detoxifying exogenous toxins,which leads to increased resistance or sensitivity of the embryonic heart to these environmental factors[9,10].It is well known that phase I metabolic enzymes are involved in the metabolic transformation of exogenous compounds,and their activation effects are directly related to the sensitivity of the body to environmental toxins.Cytochrome P450 (CYP450)enzymes are important phase I metabolic enzymes with a strong catalytic effect that promotes the metabolism and activation of various drugs,lipids,carcinogens,and teratogens;thus,impacting their possible carcinogenic and teratogenic effects[11,12].Single nucleotide polymorphisms (SNPs) ofCYP450genes can cause variability in CYP450 activity[13].Previous studies have considered the role ofCYP450genes in the metabolism of drugs and the relationship betweenCYP450SNPs and cancer susceptibility in adults[14,15]. However,the relationship between maternalCYP450SNPs and CHDs in offspring has not been explored.Thus,we hypothesized that genetic variations inCYP450genes might increase the susceptibility to drug intakerelated fetal diseases,including CHDs.

We conducted this study with the following objectives:(i) to examine whether maternal drug exposure was significantly associated with risk of CHDs in offspring;(ii) to assess the association between genetic variations in maternalCYP450genes and the risk of CHDs in offspring;and (iii) to analyze the interactions between maternal drug exposure andCYP450genetic variants on CHDs.

METHODS

Recruitment of the Study Participants

A hospital-based case-control study was performed from November 2017 to January 2020 to assess the association between maternal drug use,CYP450gene polymorphisms,and their interactions with the risk of CHDs in offspring.This study was approved by the Ethics Committee of the Xiangya School of Public Health of Central South University,and written informed consent was obtained from all mothers.This study was registered at the Chinese Clinical Trial Registry Center (registration number:ChiCTR1800016635).The research procedure and design have been described in our previous studies[16,17].The study participants were recruited from two Hunan Children’s Hospital clinics.The case group was recruited from the Department of Cardiothoracic Surgery,which provides diagnosis,treatment,surgery,and management for CHDs,and the control group was recruited from the Department of Child Healthcare after health counseling or a medical examination.The controls were selected from the same hospital during the same study period as the cases.

Inclusion and Exclusion Criteria

The outcome of interest was CHDs.Children with CHD and their mothers were identified as the case group. Children without any congenital malformations after a medical examination and their mothers were identified as the control group.The CHD cases included the following subtypes:atrial septal defect (ASD),ventricular septal defect (VSD),atrioventricular septal defect (AVSD),patent ductus arteriosus (PDA),aortopulmonary septal defect(APSD),tetralogy of Fallot (TOF),and complete transposition of the great arteries (CTOGA).Nonsyndromic CHD was of interest,but patients with other organ malformations or known abnormalities were excluded.All CHD patients were diagnosed by ultrasonography.Similarly,the control children were confirmed to have no malformations by ultrasonography.The exposures of interest were maternal drug use and genetic variants in theCYP450genes.To minimize recall bias of exposure by mothers,all cases and controls were recruited when the children were ＜ 1-year-old.Additionally,this study was aimed at mothers of Han Chinese descent with singleton pregnancies. Eligible mothers completed a questionnaire and provided informed consent and blood samples.We excluded mothers who achieved pregnancy by assisted reproductive technology,includingin vitrofertilization or intracytoplasmic sperm injection.

Information Collection

A standardized questionnaire was used to collect the corresponding maternal drug use information by specially trained investigators.As the critical stage of fetal heart development is the early stage of pregnancy,we investigated the maternal drug use information from 3 months before pregnancy to the first trimester of pregnancy.In our study,the drugs of interest were oral contraceptives,ovulatory drugs(e.g.,clomiphene,luteinizing hormone-releasing hormone,and bromocriptine),macrolide antibiotics(e.g.,erythromycin,azithromycin,clarithromycin,roxithromycin,dierythromycin,and fluoroerythromycin),antidepressants (e.g.,clomipramine,citalopram,sympathomimetic decongestants,paroxetine,and phenylpropanolamine),traditional Chinese medicines,and antiabortifacients (e.g.,progesterone preparations,chorionic gonadotropin,human blood immunoglobulin,and indomethacin).Additionally,we also collected data regarding the mothers’sociodemographic characteristics (i.e.,child-bearing age,education level,annual family income in the past year,residential area);obstetric history (i.e.,adverse pregnancy outcomes and pregnancy-related complications);family history (i.e.,consanguineous marriage and congenital malformations);personal medical history (i.e.,congenital malformations,and cold or fever history in the 3 months before this pregnancy);personal lifestyle and habits (i.e.,history of active and passive smoking in the 3 months before this pregnancy,and drinking history in the 3 months before this pregnancy); exposure history to environmentally hazardous substances (i.e.,radioactive substances,house decorations,and harmful chemicals);and folate supplementation.After completing the questionnaire,we consulted the perinatal health care handbook (PHCH) and maternal medical records to further confirm the corresponding information.Notably,not all the information we needed was recorded by the PHCH.Therefore,we used face-to-face interviews and consulted the PHCH to ensure that the relevant information was accurately collected.When the information provided by these two methods was inconsistent,we relied on the information provided by the PHCH,considering that there may be recall bias during a face-to-face interview.These measures may have helped to reduce recall bias.

Genotyping

After completing the questionnaire,all mothers were required to provide blood samples after consent.Three to five milliliters of peripheral venous blood were collected by a qualified phlebotomist from all participants using a 21-gage needle.Plasma and blood cells were separated by centrifugation and stored at -80 °C for later laboratory testing.The genomic DNA was extracted from blood cells using the QIAamp DNA Mini Kit (Qiagen,Valencia,CA,USA) according to the manufacturer’s protocol and dissolved in sterile TBE buffer.

We chose four SNPs of theCYP450genes(CYP1A1at rs1048943 and rs4646903,CYP2D6at rs1065852 and rs16947),which have been associated with adverse pregnancy outcomes and metabolism of drugs of interest[18,19],as candidate loci.Because the related primers were unavailable,we used rs4646421 instead of rs4646903 (r2=1.000),rs5751210 instead of rs1065852 (r2=0.900),and rs4147641 instead of rs16947 (r2=0.965)according to linkage disequilibrium analysis.All SNPs(rs1048943,rs4646421,rs5751210,and rs4147641)were genotyped using the matrix-assisted laser desorption and ionization time-of-flight mass spectrometry Mass Array system (Agena iPLEXassay,San Diego,CA,USA).Different cycling conditions were used to optimally amplify the target sequences.The details of polymerase chain reaction primers are described in Supplementary Table S1 (available in www.besjournal.com).The laboratory technician,who performed the genotyping,retyped and doublechecked each sample,and recorded the genotype data,was blinded to the status of the control and case groups.The genotyping error rate was ＜ 5%.

Statistical Analysis

Unordered categorical variables are described using percentages and compared between two groups using Pearson’s chi-square test or Fisher’s exact test,as appropriate.The Wilcoxon rank-sum test was used for ordinal categorical variables.Hardy-Weinberg equilibrium (HWE) was tested in the control group (significance level atP＜ 0.05).Odds ratios (ORs) and their 95% confidence intervals(CIs) were used to assess the level of the associations between maternal drug exposure and CYP450 genetic variants with the risk of CHDs in offspring.The crudeORs were calculated by univariate logistic regression.AdjustedORs (aORs) were calculated by multivariate logistic regression.Moreover,all SNPs were estimated under the three genetic models(dominant,recessive,and additive).

Additionally,we used logistic regression analysis and controlled for other influencing factors to examine the effects of the environment-gene interaction of the maternal drug exposure and theCYP450genes for the risk of CHDs in offspring.We determined the patterns and implications of the environment-gene interactions according to a method introduced by Wallace[20].Interactions were determined using interaction coefficients (γ).Theγvalues were calculated using the regression coefficient (β) from logistic regression analysis (i.e.,γ1=βe*g/βeandγ2=βe*g/βgforthegeneenvironment interaction).γvalues ＞ 1 indicated a positive interaction;andγvalues ＜ 1 indicated a negative interaction;when theγvalue=1,there was no interaction.To minimize type I error,the false discovery rateP-value (FDR_P),which was adjusted for multiple testing,was estimated to obtain a more preciseP-value[16].All tests were performed for a two-sidedP-value not exceeding 0.05.All analyses were performed using SAS 9.1 software (SAS Institute,Cary,NC,USA).

RESULTS

Baseline Characteristics

A total of 569 eligible mothers were included in the case group,and 652 were included in the control group from November 2017 to January 2020.Among the 569 CHD cases,95 (16.7%) were diagnosed with ASD,353 (62.0%) with VSD,60 (10.5%) with AVSD,170 (29.9%) with PDA,8 (1.4%) with APSD,32 (5.6%)with TOF,and 2 (0.4%) with CTOGA.Some cases were diagnosed with multiple CHD subtypes.Therefore,the sum of the various subtypes was not equal to 569.Maternal baseline characteristics in the case and control groups are summarized in Table 1.Significant differences in education level,annual family income,residential area,history of adverse pregnancy outcomes,pregnancy-related complications,consanguineous marriage,and congenital malformation history,drinking history before pregnancy,active or passive smoking history before pregnancy,cold and fever history before this pregnancy,and folate supplementation were observed between the two groups (allP＜ 0.05).These factors were adjusted when assessing the associations between maternal drug use,theCYP450gene polymorphisms,and their interactions with the risk of CHDs in offspring.

Maternal Drug Use and the Risk of CHDs in Offspring

The associations between maternal drug use and the risk of CHDs in offspring are summarized in Table 2.After adjusting for other influencing factors,the results show that mothers who used ovulatory drugs(aOR=2.12,95%CI:1.08–4.16),antidepressants(aOR=2.56,95%CI:1.36–4.82),antiabortifacients(aOR=1.55,95%CI:1.00–2.40),or traditional Chinese medicines (aOR=1.97,95%CI:1.26–3.09)during this pregnancy were at a significantly higher risk of CHDs in offspring compared with the reference group.However,we did not observe a significant association between oral contraceptive or macrolide antibiotic use and the risk of CHDs.

Maternal CYP450 Genetic Polymorphisms and the Risk of CHDs in Offspring

The genotype frequencies for each SNP of the maternalCYP450genes andP-values for the HWE test in the control group are summarized in Table 3.The genotype distributions in the control group were within HWE (allP＞ 0.05).The associations between maternalCYP450genetic variants and the risk of CHDs in offspring based on logistic regression analyses are summarized in Table 4.After adjusting for influencing factors,the results suggested that maternalCYP450genetic polymorphisms at rs1065852 and rs16947 were significantly associated with the risk of CHDs in offspring.Mothers with the A/T (aOR=1.53,95%CI:1.10–2.14) and T/T (aOR=1.57,95%CI:1.07–2.31) genotype for rs1065852 compared with those with the A/A genotype experienced a significantly increased risk of CHD in offspring.The dominant (aOR=1.55,95%CI:1.13–2.12) and additive models (aOR=1.25;95%CI:1.03–1.52)were significantly associated with an increased risk of CHDs in offspring.

Mothers with the G/G genotype for rs16947 were at a significantly higher risk of CHDs in offspring than those with the C/C genotype (aOR=3.41,95%CI:1.82–6.39);and the dominant (aOR=1.47,95%CI:1.10–1.95),recessive (aOR=3.23,95%CI:1.73–6.02),and additive models (aOR=1.51,95%CI:1.20–1.90) were significantly associated with an increased risk of CHDs in offspring.

Gene-environment Interaction Associated with the Risk of CHDs in Offspring

The gene-environment interactions between the maternalCYP450genetic variants and drug use for the risk of CHDs in offspring are summarized in Table 5.In the interaction analyses,we only considered those drugs and SNPs that were significantly associated with the risk of CHDs in the above analyses.Significant interactive effects were observed between the maternal genetic variants at rs1065852 and drug use for the development of CHDs in offspring.Overall,when mothers had a risk genotype(A/T or T/T) at rs1065852,the risk of CHDs in offspring increased significantly if they used anantidepressant (aOR=8.32,95%CI:3.56–19.45) or an antiabortifacient (aOR=2.47,95%CI:1.48–4.10).

Table 1. Maternal baseline characteristics in the case and control groups

There were also statistically significant interactive effects between maternal genetic variants at rs16947 and drug use for the development of CHDs in offspring.Overall,mothers with a risk genotype (C/G or G/G) at rs16947 had a significantly higher risk of CHDs in offspring if they used ovulatory drugs (aOR= 6.49,95%CI:1.82–23.14) or traditional Chinese medicines (aOR=3.55,95%CI:1.72–7.29).

DISCUSSION

The past few years have seen an increased interest in exploring the etiology of CHDs.Although some genetic and environmental factors have been identified,the etiologies of most non-syndromic CHDs are still unknown.In the present study,we examined whether maternal drug use before or during early pregnancy would increase the risk of CHDs in offspring.We assessed the associations between maternal CYP450 genetic variants and the risk of CHDs in offspring and analyzed the interactions between maternal drug use and the CYP450 genetic variants on the risk of CHDs.To our knowledge,this is the first exhaustive exploration of the association between maternal drug use and CYP450 genetic polymorphisms,and their interactions with the risk of CHDs in offspring,which provides new insight into the etiology of CHDs,and guides pregnant women with a high risk of CHDs touse drugs reasonably during pregnancy.

Table 3. Maternal cytochrome P450 genotype frequencies and P values of HWE test

The findings from the present study indicate that using ovulatory drugs,antidepressants,antiabortifacients,or traditional Chinese medicines before or during the early stage of pregnancy was significantly associated with the risk of CHDs in offspring.However,our study showed that using oral contraceptives and macrolide antibiotics was not associated with the risk of CHDs in offspring.Mothers who used ovulatory drugs for this pregnancy compared with those who did not were at a significantly increased risk of CHDs in offspring(OR=2.12),which is generally consistent with previous studies[21,22].Clomiphene,which is often used as a first-line infertility drug,is associated with the risk of congenital defects,including CHDs.For example,one study that used data from the National Birth Defects Prevention Study 1997–2005 to assess the association between clomiphene use and birth defects showed that clomiphene use before pregnancy significantly increases the risk of several malformations,including anencephaly,Dandy-Walker malformation,septal heart defects,muscular ventricular septal defect,coarctation of the aorta,esophageal atresia,cloacal exstrophy,craniosynostosis,and omphalocele[21].K?llén et al.reported a similar result[22]. Additionally,our previous cohort study showed that preterm birth,low birth weight,and congenital defects increase significantly among mothers who use ovulatory drugs compared with a reference group[23].

Table 4. Maternal cytochrome P450 genetic polymorphisms and risk of congenital heart defect in offspring

One study showed no association between the use of most antibiotics and cardiovascular defects[22],which is consistent with our study.However,some studies have reported that using macrolides,such as erythromycin,during early pregnancy significantly increases the risk of cardiovascular defects in offspring[24-27].Of note,as maternal viral infections,such as rubella virus and cytomegalovirus,have been confirmed to be significantly associated with the risk of CHDs[6],we cannot rule out the possibility that a teratogenic effect of an antibiotic may be due to a viral infection.The relationship between exposure to antidepressants and birth defects has always been the focus of attention.The present study showed that antidepressants were significantly associated with the risk of CHDs in offspring,which was supported by previous studies[22,28,29].For example,the risk of CHDs increases significantly among mothers who use clomipramine during early pregnancy[22].Additionally,several meta-analyses have confirmed the association between antidepressant use and the risk of CHDs[28,29].

During pregnancy,antiabortifacients,such as dydrogesterone and chorionic gonadotropin,are often used to prevent early miscarriage and preterm labor.Our study supports a positive association between antiabortifacients and the risk of CHDs,which is consistent with a previous study[30].Dydrogesterone is a progestogen that is teratogenic,and some studies have demonstrated cardiac anomalies after progestin use during pregnancy[22,31].For example,Zaqout et al.showed that taking dydrogesterone during early pregnancy is significantly associated with the risk of CHDs[31].K?llén et al. reported a significant positive association between the administration of chorionic gonadotropin and CHDs[22].Although there are some different views,previous studies have generally found that maternal exposure to various drugs,including ovulatory drugs,antiabortifacients,and antidepressants,increases the risk of CHDs in offspring.However,whether these associations between maternal drug use and CHDs are due to the underlying diseases themselves,true drug effects,or combinations of underlying diseases and drugs remains uncertain and warrants further research.

In our study,we also assessed the SNPs of maternalCYP450genes (i.e.,CYP1A1gene at rs1048943 and rs4646903 as well as theCYP2D6gene at rs1065852 and rs16947) associated with the risk of CHDs in offspring.The present results indicate that genetic variants of maternalCYP450genes could play an important role in the development of the fetal heart.After adjusting for other influencing factors,the results suggested that maternalCYP450genetic polymorphisms at rs1065852 (A/Tvs.A/A:OR=1.53;T/Tvs.A/A:OR=1.57) and rs16947 (G/Gvs.C/C:OR=3.41) were significantly associated with an increased risk of CHDs in offspring. The importance of these results depends on the finding that theCYP450genes might be potential CHD susceptibility genes.A specific CHD gene has not been found,although much research has been conducted on this topic.

CYP450 enzymes,as a superfamily of activation enzymes,exist in many cell types.According to the sequence homology of two-fifths or more,CYP450 enzymes are divided into 18 families,and according to homologies ≥ 55%,they are categorized into 44 subfamilies.CYP1A1 and CYP2D6,the enzymes included in our study,are involved in drug metabolism.WhenCYP450genes mutate,the activity of the CYP450 enzymes is higher than that coded by the wild-type,which may change susceptibility to disease[32].Although many studies have assessed the associations betweenCYP450SNPs and the risk of various cancers[33-35]and specific adverse pregnancy outcomes[36,37],the results have often been inconsistent.For example,a metaanalysis showed that theCYP1A1rs1048943 and rs4646903 polymorphisms contribute to the risk of laryngeal cancer in the Asian population[38].However,a study published in 2016 showed that theCYP1A1rs1048943 polymorphism is not significantly associated with the risk of gastric cancer[35].Some differences in the results have been reported for other diseases[33,36-37].In the present study,no significant positive associations were observed between theCYP1A1rs1048943 and rs4646903 polymorphisms and the risk of CHDs.However,our study showed that genetic polymorphisms at rs1065852 and rs16947 were significantly associated with the risk of CHDs.No study has evaluated the association betweenCYP450genetic polymorphisms and the risk of CHDs,until now.Undoubtedly,our study will provide a new clue for screening candidate genes for CHDs.

The most impressive findings of this study were the interactions between maternal drug use and theCYP450genetic polymorphisms in the development of CHDs.The interaction analyses showed that mothers with aCYP450gene risky genotype at rs1065852 and rs16947 had a significantly increased risk of CHDs in offspring if they had used ovulatory drugs (OR=6.49),antidepressant (OR=8.32),antiabortifacients (OR=2.47),or traditional Chinese medicines (OR=3.55),indicating that maternalCYP450genes may regulate the effects of the corresponding drugs on embryonic heart development. However,the interactions between theCYP450genetic polymorphisms and drug use on the development of CHD were different,which was not surprising given the differences in the teratogenic effects of different drugs on fetal CHDs,as well as the differences in the degradation of different drugs by CYP450 enzymes.Teratogenesis is a process accompanied by complex steps and various factors,indicating different gene mutations and several biological pathways.Therefore,interactive effects of different factors exist in the development of CHDs.One view suggests that genetic polymorphisms in metabolic enzymes play a pathogenic role by altering the effects of environmental factors rather than directly causing disease[39].We supposed that the genetic variants in theCYP450genes promote metabolic activation of exogenous compounds,including drugs,by coding higher activity in the CYP450 enzymes.Then,teratogenic intermediate drug metabolites accumulate in the maternal body and influence the development of the embryo.Although previous studies have shown that genetic polymorphisms inCYP450genes may affect the metabolism of drugs[40,41],no study has assessed the interactions between maternalCYP450genes and drug use in the development of CHD until now.However,the specific mechanism of the combined effect between theCYP450genes and drug use on the risk of CHDs remains unclear and needs further research.There is increasing evidence that geneenvironment interactions affect CHDs,so the effects of these interactions will be the future research direction for CHDs.

Some limitations in our study should be discussed.First,our study may have suffered from selection bias.For example,a convenience sample,driven mainly by the number of respondents,was used for our study.The target population did not include mothers whose pregnancy was terminated because of a CHD.This limitation could lead to subsequent problems,including sample representativeness and generalization of the findings.Second,we replaced threeCYP450gene SNPs according to linkage disequilibrium,which may have overestimated the statistical correlation in our findings.Third,residual confounders are always a concern in observational studies.Although we adjusted for a wide range of factors influencing CHDs,we could not exclude the possibility that other unmeasured or inadequately measured factors confounded the true associations.Fourth,as this was a case-control study,we could not verify causality.Fifth,although some measures were adopted,recall bias was inevitable.Sixth,because of the limited sample size,we analyzed drugs in a wide range of categories rather than specific drugs,and we did not assess the risk of specific CHD phenotypes.Additionally,the frequency of drug use was not measured,and exposure to a drug was only dichotomized as yes or no independently.Finally,because the present study only included Han Chinese subjects,additional research in other populations is warranted to generalize the findings.These limitations highlight the need for larger samples,a prospective approach,and different ethnic populations.

In conclusion,our study indicates that using ovulatory drugs,antidepressants,antiabortifacients,or traditional Chinese medicines before or during the early stage of pregnancy may be significantly associated with an increased risk of CHDs in offspring of Chinese descent.Additionally,this study supports a significant positive association between maternalCYP450genetic variants at rs1065852 and rs16947 and the risk of fetal CHDs.Significant interactions were observed between the maternalCHY450genetic polymorphisms and drug use,suggesting that maternalCYP450genes may regulate the effects of maternal drug exposure on fetal heart development.However,the mechanism by which these factors affect the development of the fetal heart remains unknown,and more studies in different ethnic populations and with a larger sample are required to confirm these findings.

ACKNOWLEDGMENTS

The authors thank the editors and reviewers for their suggestions and all colleagues working in the Maternal and Child Health Promotion and Birth Defect Prevention Group.

DECLARATION OF INTEREST STATEMENT

The authors declare no conflicts of interest.

AUTHOR CONTRIBUTIONS

JBQ,JYD,and PZ designed and supervised the project.LL,MTS,PH,TTW,JYD,SMZ,JQL,YHL,and LTC collected and analyzed the data.JBQ,JYD,and PZ conducted the experiments.JBQ and JYD wrote the manuscript.JBQ and PZ edited the manuscript.

Supplementary Table S1.Primer sequences for cytochrome P450 genetic polymorphisms

Received:January 5,2021;

Accepted:August 17,2021