ABSTRACT
Background Intrauterine exposure to antiseizure and antidepressant drugs is associated with adverse pregnancy outcomes, including low birth weight. Gene variants in placental efflux transporter genes may alter foetal exposure to these drugs.
Methods We investigated whether genetic variants in placental efflux transporters modified the impact of maternal antiseizure and antidepressant drug use on offspring birth weight, using data from the Norwegian Mother, Father and Child Cohort Study and the Medical Birth Registry of Norway (69,828 offspring with their genotype data, 81,189 with maternal genotype data). We systematically searched for gene variants in placental efflux transporters influencing drug exposure [MDR1-ABCB1 (7 alleles), MRP1-ABCC1 (2), MRP2-ABCC2 (3) and, BCRP-ABCG2 (2)] and calculated genetic scores (sum of alleles potentially related to low transporter activity). We assessed the interaction between prenatal drug use and genetic scores on birth weight.
Findings Prenatal antiseizure medication exposure was associated with lower birth weight (-95.5 grams, 95% CI -190 to -0.78). This relationship depended on the offspring’s MRP2-ABCC2 genetic score. Birth weight differences between exposed vs. unexposed were 70.3 g (95% CI -494 to 634) in the lowest genetic score category and -306 g (95% CI -361 to -31.8) in the highest (interaction p-value = 0.019; main variant: rs3740066). The antiseizure medication-lower birth weight association also depended on the maternal MDR1-ABCB1 genetic score. Birth weight differences between exposed vs. unexposed were -66.8 g (95% CI -225 to 91.2) in the lowest genetic score category and -317 g (95% CI -517 to -117) in the highest (interaction p-value = 0.037; main variants: rs10248420 and rs2235015). Prenatal antidepressant exposure was associated with low birth weight, but no gene-drug interactions were observed.
Conclusions MRP2-ABCC2 and MDR1-ABCB1 variants may influence prenatal antiseizure medication’s impact on neonatal birth weight.
What is already known on this topic
- Low birth weight increases the risk of neonatal death, morbidity, and lifelong susceptibility to various health conditions, including neurodevelopmental, cardiometabolic, respiratory, and psychiatric disorders.
- Previous studies have reported an association between maternal use of centrally acting medication during pregnancy and low birth weight in the offspring, but the underlying causes remain poorly understood.
- It is of great concern to identify strategies that can modulate foetal exposure to these medications and enable individualised, safe use since discontinuation may also increase the risk of harm to the mother.
What this study adds
- Our study uncovered genetic variants in both offspring (rs3740066) and mother (rs10248420 and rs2235015) linked to placental efflux transporters that modulated the association between maternal antiseizure drug use during pregnancy and the risk of low birth weight in the offspring, marking the first exploration of such interactions.
- If these findings are confirmed and information about specific drugs is incorporated, this could potentially be used in the development of personalized recommendations for pregnant women undergoing these treatments.
INTRODUCTION
The use of medication during pregnancy must carefully balance benefits to the mother and potential harm to the offspring.1 Centrally acting drugs (such as antiseizure drugs [mainly lamotrigine, 0.3% of pregnancies] and antidepressants [1.5% of pregnancies]) have been increasingly used during pregnancy.2,3 However, intrauterine exposure to these medications may be associated with adverse pregnancy outcomes such as impaired foetal growth and reduced birth weight in offspring.4–6 Furthermore, these conditions are linked to a multitude of morbidities, resulting in both immediate and long-term adverse consequences for newborns, their families, and society.7 The placenta plays an essential role in substance exchange between mother and foetus (e.g., nutrients, metabolic by-products). Small, non-polar, lipophilic xenobiotics such as antiseizure drugs and antidepressants can pass through the placenta via passive diffusion.8 However, several active proteins can efflux these compounds from the foetus back to the mother.9 Active efflux proteins from the adenosine triphosphate-binding cassette superfamily (such as P-glycoprotein, breast cancer resistance protein, and others multidrug resistance proteins) are involved in the transplacental passage of drugs particularly affecting the passage of centrally-acting drugs and their metabolites.10,11 Several antiseizure drugs and antidepressants are substrates for these transporters.12,13 Efflux transporters may therefore have a foetoprotective effect by lowering the concentrations of potentially toxic drugs on the foetal side. Genetic variants in the transporters may modify their activity and could be located in the offspring or maternal genotype. 14,15 Little is known about the role of genetic variants in placental efflux transporters on the modulation of the association of antiseizure drugs and antidepressants during pregnancy with offspring outcomes such as low birth weight.
Our aim was to explore whether there is an interaction between the use of antiseizure drugs or antidepressants during pregnancy and genetic variants related to placental efflux transporters (in P-glycoprotein, breast cancer resistance protein, and others multidrug resistance proteins) on offspring birth weight.
MATERIALS AND METHODS
Study participants
We included participants in the Norwegian Mother, Father, and Child Cohort Study (MoBa). MoBa is a prospective, population-based pregnancy cohort conducted by the Norwegian Institute of Public Health. Pregnant women and their partners were recruited across Norway between 1999-2008 at the time of routine ultrasound screening (∼17th gestational week). The cohort includes approximately 114,000 children, 95,000 mothers, and 75,000 fathers.
This work used a subsample of offspring from singleton pregnancies with available information on genotype and birth weight (from version #12 of the quality-assured data files released on May 11, 2022). Genotype data was obtained from blood samples provided during pregnancy and at birth (mothers) and umbilical cord blood at birth (offspring).16 A total of 238,001 samples have been genotyped in 24 genotyping batches with varying selection criteria, genotyping platforms, and genotyping centres.17 Quality control, phasing, and imputation were performed using the MoBaPsychGen pipeline as previously described.17
Our work is described according to the Strengthening the Reporting of Observational studies in Epidemiology guidelines (Supplemental Table 3).18
Use of centrally acting drugs
We used self-reported information on the use of centrally acting drugs from MoBa questionnaires in the 18th and 30th gestational weeks.19 Specifically, women were asked whether they had epilepsy or depression, and if they answered yes, whether any medications were used to treat the condition(s). This project has no information on the specific drug(s) the woman used. Any reported use of medication(s) for epilepsy (yes/no) or depression (yes/no) were categorized and evaluated separately.
Birth weight
Data on offspring birth weight in grams was obtained from the Medical Birth Registry of Norway, a national health registry which contains information about all births in Norway since 1967.20
Covariates
Information on maternal age at delivery (continuous), parity (1, 2, 3, ≥ 4) and sex of the offspring (female/male) was obtained from the birth registry. Furthermore, information on maternal years of education (continuous), pre-pregnancy body mass index (continuous), having ever smoked (yes/no), and self-reported depression or epilepsy was gathered in the questionnaire administered at gestational week 18.
Genetic variants related to placental transporters and genetic scores
We performed a systematic review in PubMed, Web of Science, pharmacogenetic websites (https://www.pharmgkb.org) and GWAS catalog (https://www.ebi.ac.uk/gwas/) to identify gene variants on the adenosine triphosphate-binding cassette superfamily of placental efflux transporters (MDR1-ABCB1, MRP1-ABCC1, MRP2-ABCC2, and BCRP-ABCG2). We followed the recommendations of the Preferred Reporting Items for Systematic reviews and Meta-Analyses statement.21 A detailed description of the search is available in Supplemental Materials.
After the systematic review, we only conserved genetic variants that were not in linkage disequilibrium according to linkage disequilibrium block analysis in the 1000 Genomes CEU and GBR populations (given the MoBa participants passing the post-imputation quality control cluster with these populations; R2 < 0.8) and with a minor allele frequency > 1%.22 Genetic scores were calculated as the sum of the number of risk alleles for all genetic variants related to each transporter in a participant. The genetic score was represented in the results as a genetic score if the number of score values were ≤5. Otherwise, we represented them in quartiles. We considered a risk allele to be the allele associated with any of these five theoretical situations potentially related to an increased exposure of the drug to the foetus: decreased function of the transporter, declined expression of the transporter, non-resistance to the treatment, higher concentrations of the drug, or more adverse effects. Information related to the studies is available in Supplemental Table 2.
Ethical approval
The MoBa study is conducted according to the Declaration of Helsinki for medical research involving human subjects. The establishment of MoBa and initial data collection was based on a license from the Norwegian Data Protection Agency. It is now based on regulations related to the Norwegian Health Registry Act. Participants provided written informed consent before joining the cohort. This project was approved by the Regional Committee for Medical and Health Research Ethics of South/East Norway (reference: 2017/1362).
Statistical analyses
Normally distributed continuous variables were described by means and standard deviations, nonnormally distributed continuous variables by medians and 25th-75th percentiles, and categorical variables by proportions.
We first assessed the relationship between maternal use of centrally acting drugs during pregnancy and birth weight using multivariable linear regressions adjusted for offspring sex and maternal factors (age at delivery, years of education, pre-pregnancy body mass index, parity and having ever smoked). Clustered standard errors were computed to account for dependence among births to mothers contributing with more than one pregnancy in MoBa. To determine whether there were significant interactions between maternal prenatal use of centrally acting drugs and the genetic scores on the offspring birth weight, we applied likelihood ratio tests between nested linear regression models with and without an interaction product-term of “exposure group × genetic score”. The nested models were further adjusted for the first twenty ancestry-informative genetic principal components and genotyping batch. We considered any interaction with a p-value for the overall interaction test < 0.05 and a differential association between medication use and offspring birth weight in the two extreme categories according to the genetic score. Whether an interaction was found, we further explored interactions with the individual variants in the genetic score as sensitivity analyses using the same strategy (as these were exploratory sensitivity analyses, we here considered any interaction with a p-value < 0.1).
Statistical analyses were performed in R Software, version 4.1.0. Our analysis code is available in GitHub: https://github.com/alvarohernaez/Gene_psychoactivedrug_BW_MoBa/.
RESULTS
Study population
In our study, we included 69,828 singleton pregnancies with offspring genotype information and offspring birth weight, and 81,189 singleton pregnancies with maternal genotype data and offspring birth weight (Figure 1). Mean birth weight was 3,639 grams (SD 522). A total of 174 children (0.25%) were exposed to maternal use of antiseizure drugs during pregnancy and 766 (1.10%) to antidepressants. Women using centrally acting drugs were younger (only for antiseizure drugs), had lower educational attainment, and were more likely to have ever smoked (Table 1).
Search of genetic variants on placental transporters
The systematic search identified 26 genetic variants in the placental transporter genes associated with differences in antiseizure drugs and/or antidepressant outcomes (Supplemental Figure 1), of which 14 were independent and available in the MoBa genotype database (Supplemental Table 1). Seven variants were found for MDR1-ABCB1, two for MRP1-ABCC1, three for MRP2-ABCC2, and two for BCRP-ABCG2. We used these genetic variants to calculate risk scores. The ranges of values for the four genetic scores calculated were: 0 to 11 (MDR1-ABCB1; we presented the stratified associations according to quartiles), 0 to 4 (MRP1-ABCC1), 2 to 6 (MRP2-ABCC2), and 0 to 4 (BCRP-ABCG2).
Interaction between genetic variants of placental transporters and exposure to antiseizure drugs during pregnancy on birth weight
Antiseizure medication use during pregnancy was associated with lower birth weight (-95.5 g, 95% confidence interval [CI] -190 to -0.78). Greater values of the MRP2-ABCC2 genetic score in the offspring were associated with lower birth weight in offspring of mothers who used antiseizure drugs during pregnancy. The difference in birth weight between exposed vs. unexposed was 70.3 g (95% CI -494 to 634) in the lowest genetic score category (2 risk alleles) and -306 g (95% CI -361 to -31.8) in the highest (6 risk alleles) (p-value for interaction = 0.019; Figure 2). The variant with the strongest evidence of an interaction in the sensitivity analyses of the MRP2-ABCC2 genetic score was rs3740066 (p-value for interaction = 0.023; Supplemental Figure 2).
Greater values of the MDRP1-ABCB1 genetic score in the mother were also linked to lower birth weight among offspring of mothers who used antiseizure drugs during pregnancy. The difference in birth weight between exposed vs. unexposed were -66.8 g (95% CI -225 to 91.2) in the lowest genetic score category (first quartile) and - 317 g (95% CI -517 to -117) in the highest (fourth quartile) (p-value for interaction = 0.037; Figure 3). Variants rs10248420 and rs2235015 included in the MDRP1-ABCB1 genetic score had some evidence of an interaction in the sensitivity analyses (p-values for interactions = 0.042 and 0.058, respectively; Supplemental Figure 3).
Interaction between genetic variants of placental transporters and exposure to antidepressants during pregnancy on birth weight
Use of antidepressant drugs during pregnancy was associated with lower birth weight (-60.5 g, 95% CI -97.6 to -23.3). However, we observed no interactions between prenatal use of antidepressants and any genetic score of placental efflux transporters (in offspring or mothers) on birth weight in the offspring (Supplemental Figures 4 and 5).
DISCUSSION
In this large population-based study, involving offspring of mothers who used antiseizure medication during pregnancy, a high genetic score for MRP2-ABCC2 in the offspring and MDR1-ABCB1 in the mother, suggesting a greater exposure to these drugs, was associated with lower birth weight in the offspring.
MDR1-ABCB1 and MRP2-ABCC2 are the only transporters located in the apical membrane of the trophoblast (maternal side) and their function is to return xenobiotics to maternal circulation.23 Thus, it seems plausible that genetic variants linked to lower transporter activity are associated with higher foetal exposure and greater toxicity of antiseizure drugs. In addition, the main variant in the interaction between the MRP2-ABCC2 genetic score and prenatal antiseizure drug use (rs3740066) has been previously linked to a poorer response to antiseizure medication24 and a higher risk of other adverse events during pregnancy (intrahepatic cholestasis of pregnancy).25 The differences in the interactions found within genetic variants in the offspring (MRP2-ABCC2) and maternal genome (MDR1-ABCB1) can be explained by physiological changes in the placenta structure. Early in pregnancy, the separation between the maternal circulation and the foetal circulation by the trophoblasts villous barrier is 50-100 μm (second month). It becomes progressively thinner as the pregnancy proceeds (only 4-5 μm at term).26 In addition, the total placenta surface area increases from about 5 m2 at 28 weeks of gestation to 12 m2 at term. Thus, passive diffusion of medications increases with gestational age, which means that medication can reach the foetal circulation more easily later in pregnancy.26 At the same time, the expression of MRP2-ABCC2 increases with advancing gestational age, whereas the expression of MDR1-ABCB1 declines.27,28 As a consequence of all these phenomena, the presence of risk alleles for MRP2-ABCC2 may play a more significant role in determining the offspring’s exposure. However, an alternative hypothesis could also explain our findings. A greater foetal exposure to medication due to lower activity of the efflux transporters may also be linked to lower medication levels in the mother and a poorer control of the disease. Greater disease effects in mothers may be associated with alternative mechanisms that may also explain an increased risk of low birth weight in the offspring.
There could be several reasons for the lack of gene-drug interactions with antidepressants. First, the association between exposure to antidepressants during pregnancy and birth weight is of a smaller magnitude in our data, which can make the search of gene-drug interactions more challenging. Second, the accuracy of our data regarding the usage of antidepressants may be less precise, as women may have varying criteria for identifying such treatments. This discrepancy may not be present regarding antiseizure drugs. Finally, a larger proportion of women may discontinue the use of antidepressants during pregnancy compared to antiseizure medications.29
Our study has some limitations. First, due to the unavailability of information on the specific drug utilized, we conducted our analysis based on aggregated drug groups. Although the main medications in each group show similar pharmacodynamic behaviour and their usage during pregnancy reduces the number of medications under consideration,3 this may interfere with our findings because there may be variations in the medication groups’ interactions due to differences in transporter affinity or dose-dependent effects.26,30 Second, due to the lack of a genome-wide association study showing which genetic variants are associated with high or low activity of efflux transporters, we calculated risk scores using studies on genetic variants indirectly related to high or low activity (e.g., treatment resistance, treatment efficacy, presence of adverse effects). Third, the interactions found in our data are only suggestive and our findings require replication in independent populations to confirm their validity. Finally, the characteristics of our participants (individuals of northern European ancestry with moderate-high socioeconomic status) limit our capacity to generalize our conclusions to other populations.
CONCLUSIONS
In summary, genetic variants in MRP2-ABCC2 and MDR1-ABCB1 placental transporters may modulate the association between prenatal exposure to antiseizure drugs and low birth weight in the offspring. To the best of our knowledge, this is the largest gene-drug interaction study on offspring birth weight to date and it was performed in a well-characterized population with genome-wide genotype information. If our findings are confirmed, particularly in studies in which we know exactly which drug the mothers have used during pregnancy, the assessment of these genetic variants may be implemented in clinical settings to support safer use of antiseizure drugs during pregnancy.
DATA AVAILABILITY STATEMENT
Consent given by the participants does not open for storage of data on an individual level in repositories or journals. Researchers who want access to datasets for replication should apply to datatilgang{at}fhi.no. Access to datasets requires approval from the Regional Committee for Medical and Health Research Ethics in Norway and an agreement with MoBa.
FUNDING
The MoBa Cohort Study is supported by the Norwegian Ministry of Health and Care Services and the Norwegian Ministry of Education and Research. This project received funding from the European Research Council under the European Union’s Horizon 2020 research and innovation program (grant agreement No 947684). This work was also partly supported by the Research Council of Norway through its Centres of Excellence funding scheme, project number 262700 and 223273, and the project “Women’s fertility – an essential component of health and well-being”, number 320656, and co-funded by the European Research Council (grant agreement No 101071773). M.H.H. was supported by the fellowship “Estancias de movilidad Jose Castillejo dentro del Programa Estatal de Promoción del Talento y su empleabilidad, en el marco del Plan Estatal de lnvestigación Científica y Técnica y de lnnovación 2017-2020”, from the Spanish government. P.R.N. was supported by the European Research Council (grant agreement No 293574), the Novo Nordic Foundation (#NNF19OC0054741) and the Trond Mohn Foundation (PRECISE-DIA). E.C. and A.Havdahl were supported by the Research Council of Norway (274611) and the South-Eastern Norway Regional Health Authority (project numbers 2020022, 2021045). O.A.A. was supported by the European Research Council (grant agreement No 964874) and the Research Council of Norway (300309, 273291). The funders had no role in the study design; the collection, analysis, and interpretation of data; in the writing of the report; or in the decision to submit the article for publication. Views and opinions expressed in this paper are those of the authors only and do not necessarily reflect those of the funders. Neither the European Union nor the granting authority can be held responsible for them.
CONFLICTS OF INTEREST
O.A.A. is a consultant for cortechs.ai and has received speaker’s honorarium from Janssen, Lundbeck, and Sunovion unrelated to the current work. The rest of the authors declare that no competing interests exist.
AUTHORS’ CONTRIBUTION
M.H.H. conceived and designed the study, was responsible for data curation, formal analyses, data interpretation, and drafting the article. J.M.C. critically revised the article and contributed to interpretation of the data. K.H.S. provided support in data analysis, software use, and visualization of results, and critically revised the article. T.K.W.G. was involved in the interpretation of data and critically revised the article. Y.L. contributed to data acquisition and critically revised the article. P.M. contributed to data acquisition and critically revised the article. P.R.N. contributed to data acquisition and critically revised the article. O.A.A. contributed to data acquisition and critically revised the article. E.C. contributed to data acquisition and critically revised the article. A.Havdahl contributed to data acquisition and critically revised the article. E. Molden critically revised the article. S.E.H. contributed to data acquisition, obtained funding, and critically revised the article. K.F. provided support in data analysis and critically revised the article. M.C.M. and A.Hernaez coordinated the project, contributed to data interpretation, and critically revised the article. M.H.H. and A.Hernaez are the guarantors of this study, accept full responsibility for the work and the conduct of the study, had access to the data, and controlled the decision to publish.
SUPPORTING INFORMATION
SUPPLEMENTAL MATERIALS
Search strategy of genetic variants
All studies that examined the association between the placental transporters of interest and epilepsy or depression were researched. We systematically searched PubMed, Web of Science, pharmacogenetic websites (https://www.pharmgkb.org) and GWAS catalog (https://www.ebi.ac.uk/gwas/home). All databases were searched using the Boolean method with the following terms (1 AND 2 AND 3): 1) “antiseizure” OR “antiepilepsy” OR “anticonvulsants” OR “antidepressant” OR “benzodiazepines”; 2) “polymorphism” OR “genetic polymorphism” OR “genetic variant” OR “pharmacogenetics”; 3) “ATP binding cassette” or “ABCB1” OR “ABCC1” OR “ABCC2” OR “ABCG2” OR “ABCC5” OR “ABCC3”. In addition, the reference and discussions of all pooled articles were carefully scanned for additional publications. We excluded studies: 1) not performed in humans; 2) not published in English; 3) not published in scientific journals; and 4) focused on other drug groups.
SUPPLEMENTAL TABLES
ACKNOWLEDGMENTS
The MoBa Cohort Study is supported by the Norwegian Ministry of Health and Care Services and the Ministry of Education and Research. We thank all the participating families in Norway who take part in this ongoing cohort study, and those who contributed to the recruitment and the infrastructure of the cohort.
We thank the Norwegian Institute of Public Health for generating high-quality genomic data. This research is part of the HARVEST collaboration, supported by the Research Council of Norway (project reference: #229624). We also thank the NORMENT Centre for providing genotype data, funded by the Research Council of Norway (project reference: #223273), South East Norway Health Authority, and Stiftelsen Kristian Gerhard Jebsen. We further thank the Center for Diabetes Research (University of Bergen) for providing genotype data funded by the European Research Council Advanced Grant project SELECTionPREDISPOSED, Stiftelsen Kristian Gerhard Jebsen, the Trond Mohn Foundation, the Research Council of Norway, the Novo Nordisk Foundation, the University of Bergen, and the Western Norway Health Authority.
This work was performed on the TSD (Tjeneste for Sensitive Data) facilities, owned by the University of Oslo, operated, and developed by the TSD service group at the University of Oslo, IT-Department (tsd-drift@usit.uio.no).
This paper does not necessarily reflect the position or policy of the Norwegian Research Council.
REFERENCES
SUPPLEMENTAL REFERENCES
- 1.
- 2.
- 3.
- 4.
- 5.
- 6.
- 7.
- 8.
- 9.
- 10.
- 11.
- 12.
- 13.
- 14.
- 15.
- 16.
- 17.
- 18.
- 19.
- 20.
- 21.
- 22.
- 23.
- 24.
- 25.
- 26.
- 27.
- 28.
- 29.
- 30.
- 31.
- 32.
- 33.
- 34.
- 35.
- 36.
- 37.
- 38.
- 39.
- 40.