Skip to main content
Advertisement
  • Loading metrics

First-trimester exposure to benzodiazepines and risk of congenital malformations in offspring: A population-based cohort study in South Korea

  • Yunha Noh ,

    Contributed equally to this work with: Yunha Noh, Hyesung Lee

    Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Visualization, Writing – original draft, Writing – review & editing

    Affiliation School of Pharmacy, Sungkyunkwan University, Suwon, South Korea

  • Hyesung Lee ,

    Contributed equally to this work with: Yunha Noh, Hyesung Lee

    Roles Data curation, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing

    Affiliations School of Pharmacy, Sungkyunkwan University, Suwon, South Korea, Department of Biohealth Regulatory Science, Sungkyunkwan University, Suwon, South Korea

  • Ahhyung Choi,

    Roles Formal analysis, Investigation, Methodology, Visualization, Writing – review & editing

    Affiliation School of Pharmacy, Sungkyunkwan University, Suwon, South Korea

  • Jun Soo Kwon,

    Roles Conceptualization, Investigation, Methodology, Writing – review & editing

    Affiliations Department of Psychiatry, Seoul National University College of Medicine, South Korea, Department of Brain and Cognitive Sciences, Seoul National University College of Natural Sciences, South Korea, Institute of Human Behavioral Medicine, Seoul National University Medical Research Center, South Korea

  • Seung-Ah Choe,

    Roles Conceptualization, Investigation, Methodology, Writing – review & editing

    Affiliation Division of Life Sciences, Korea University, Seoul, South Korea

  • Jungmi Chae,

    Roles Data curation, Formal analysis, Investigation, Writing – review & editing

    Affiliation Department of Research, Health Insurance Review and Assessment Service, Wonju, South Korea

  • Dong-Sook Kim ,

    Roles Data curation, Formal analysis, Investigation, Project administration, Resources, Supervision, Validation, Writing – review & editing

    sttone@hira.or.kr (DSK); shin.jy@skku.edu (JYS)

    Affiliation Department of Research, Health Insurance Review and Assessment Service, Wonju, South Korea

  • Ju-Young Shin

    Roles Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Supervision, Validation, Writing – review & editing

    sttone@hira.or.kr (DSK); shin.jy@skku.edu (JYS)

    Affiliations School of Pharmacy, Sungkyunkwan University, Suwon, South Korea, Department of Biohealth Regulatory Science, Sungkyunkwan University, Suwon, South Korea, Department of Clinical Research Design & Evaluation, Samsung Advanced Institute for Health Sciences & Technology (SAIHST), Sungkyunkwan University, Seoul, South Korea

Abstract

Background

Benzodiazepines are frequently prescribed during pregnancy; however, evidence about possible teratogenicity is equivocal. We aimed to evaluate the association between first-trimester benzodiazepine use and the risk of major congenital malformations.

Methods and findings

Using Korea’s nationwide healthcare database, we conducted a population-based cohort study of women who gave birth during 2011 to 2018 and their live-born infants. The exposure was defined as one or more benzodiazepine prescriptions during the first trimester. We determined the relative risks (RRs) and confidence intervals (CIs) of overall congenital malformations and 12 types of organ-specific malformations. Infants were followed from birth to death or 31 December 2019, whichever came first (up to 8 years of age). Propensity score fine stratification was employed to control for 45 potential confounders. Among a total of 3,094,227 pregnancies, 40,846 (1.3%) were exposed to benzodiazepines during the first trimester (mean [SD] age, 32.4 [4.1] years). The absolute risk of overall malformations was 65.3 per 1,000 pregnancies exposed to benzodiazepines versus 51.4 per 1,000 unexposed pregnancies. The adjusted RR was 1.09 (95% CI 1.05 to 1.13, p < 0.001) for overall malformations and 1.15 (1.10 to 1.21, p < 0.001) for heart defects. Based on mean daily lorazepam-equivalent doses, the adjusted RRs for overall malformations and heart defects were 1.05 (0.99 to 1.12, p = 0.077) and 1.12 (1.04 to 1.21, p = 0.004) for <1 mg/day and 1.26 (1.17 to 1.36, p < 0.001) and 1.31 (1.19 to 1.45, p < 0.001) for >2.5 mg/day doses, respectively, suggesting a dose–response relationship. A small but significant increase in risk for overall and heart defects was detected with several specific agents (range of adjusted RRs: 1.08 to 2.43). The findings were robust across all sensitivity analyses, and negative control analyses revealed a null association. Study limitations include possible exposure misclassification, residual confounding, and restriction to live births.

Conclusions

In this large nationwide cohort study, we found that first-trimester benzodiazepine exposure was associated with a small increased risk of overall malformations and heart defects, particularly at the higher daily dose. The absolute risks and population attributable fractions were modest. The benefits of benzodiazepines for their major indications must be considered despite the potential risks; if their use is necessary, the lowest effective dosage should be prescribed to minimize the risk.

Trial registration

ClinicalTrials.gov NCT04856436.

Author summary

Why was this study done?

  • Anxiety and insomnia are common during pregnancy, and benzodiazepines are frequently prescribed for managing these conditions.
  • The safety of benzodiazepines during pregnancy remains uncertain, as their evidence from epidemiological studies is limited and conflicting.

What did the researchers do and find?

  • In this large nationwide cohort study of more than 3 million pregnancies, we found a small increased risk of overall and heart defects associated with first-trimester benzodiazepine use.
  • The risk of overall and heart defects was slightly increased at the high daily dose group, suggesting the dose–response relationship.
  • A small but significant increased risk for overall and heart defects was detected with several specific benzodiazepines.

What do these findings mean?

  • The findings suggest that, although small, the potential risks should be evaluated against the efficacy of benzodiazepines and the lowest effective dosage should be recommended when prescribed in early pregnancy.
  • The increased risk for congenital malformations observed with several specific agents should be carefully monitored in future research as a potential safety signal.

Introduction

Anxiety, insomnia, and mood disorders are common during pregnancy [1,2], and benzodiazepines are frequently prescribed to pregnant women to manage these conditions [3]. The worldwide prevalence of benzodiazepine use during pregnancy is approximately 2% [4], and more than 1% of pregnant women in South Korea are prescribed these agents during the first trimester (S1 Fig). Despite their regular use, over the last decades, with a parallel growing focus on newer antidepressants or antipsychotics, benzodiazepines have received limited consideration, resulting in a lack of evidence on the safety of their use [5].

Benzodiazepines readily cross the human placenta and may accumulate in fetal tissues at concentrations higher than those detected in maternal serum [6,7]. The teratogenicity of benzodiazepines is biologically plausible, as they bind to receptors in peripheral tissues as well as the brain and are involved in cell proliferation and differentiation [8]. Given their potential to harm the fetus, regulatory agencies recommend that benzodiazepines should be avoided during pregnancy [9,10]. However, evidence indicating the teratogenicity of benzodiazepines is uncertain as pregnant women are usually excluded from clinical trials. Although different meta-analyses of epidemiological studies have reported no association between benzodiazepines and congenital malformations [11,12], a majority of these studies had significant methodological limitations, including low statistical power due to a small number of exposed women; potential recall and selection bias, as most previous studies were case–control designs; no detailed information regarding benzodiazepine prescriptions (e.g., dosage and indication); and lacked control for important confounders (e.g., psychiatric comorbidities, concomitant medications) [1120]. Moreover, although pharmacokinetic and pharmacodynamic profiles substantially differ among individual benzodiazepines [21], most studies have interpreted them simply as a class effect, and few studies have demonstrated a dose–response relationship.

Accordingly, to address these limitations and encourage optimal therapeutic decisions for pregnant women, additional research is needed in an adequately large-scale pregnancy cohort. Thus, we conducted a nationwide cohort study in South Korea to examine the association between maternal exposure to benzodiazepines during the first trimester and the risk of major congenital malformations in their offspring.

Methods

Data source and study cohort

We conducted a nationwide retrospective cohort study using healthcare data retrieved from the Health Insurance Review and Assessment Service (HIRA) database, which covers 50 million people (approximately 99% of the South Korean population), from 1 July 2009 to 31 December 2019. These data comprise individual-level demographics and all records of diagnosis and healthcare utilization (e.g., drug prescription and medical procedure), provided through inpatient, outpatient, and emergency department visits. In a validation study comparing our database and electronic medical records, the overall positive predictive value of diagnosis records was 82% [22].

Our cohort included all pregnancies resulting in live births from 1 January 2011 to 31 December 2018, identified with procedure codes of delivery (S1 Table). We included all live-born infants who were linked with their mothers and restricted the pregnancy cohort to women aged 20 to 45 years at delivery. We excluded pregnancies diagnosed with a chromosomal abnormality; those with exposure to known teratogenic drugs during the first trimester; and those unexposed to benzodiazepines during the first trimester, but exposed at least once during the 3 months preceding the last menstrual period, to avoid contaminating the unexposed group with women who could have taken benzodiazepines post-pregnancy (S2 Table). Infants were followed from birth to 8 years, death, or the end of the study period (December 2019), whichever came first. We calculated the date of the last menstrual period using an algorithm to estimate the gestational age in administrative databases [23].

The need for informed consent was waived, as this study was conducted using anonymized claims data. This study was approved by the Institutional Review Board of Sungkyunkwan University, South Korea (No. 2021-04-005). We registered the study protocol on ClinicalTrials.gov (NCT04856436).

Exposure

The exposed group was composed of women who filled at least one benzodiazepine prescription during the first trimester (first 90 days of pregnancy), known as the etiologically relevant period for congenital malformations. The unexposed group comprised women who were not prescribed any benzodiazepine from 3 months before the last menstrual period to the end of the first trimester. Furthermore, we evaluated short-acting (half-life ≤24 h) and long-acting (half-life >24 h) benzodiazepines based on the duration of action [24], as well as the individual medications. To assess a dose–response relationship, we calculated all benzodiazepine doses by converting them to lorazepam-equivalent doses [25,26], and we classified them into 3 mean daily dose groups: <1 mg/day, 1 to 2.5 mg/day, and >2.5 mg/day.

Outcomes

Major congenital malformations were identified by diagnostic records, according to the ICD-10 codes defined by the European Surveillance of Congenital Anomalies classification (S1 Table) [27]. Major congenital malformations were further categorized into 12 types of organ-specific malformations: (1) nervous system; (2) eye; (3) ear, face, and neck; (4) heart; (5) respiratory system; (6) oral cleft; (7) digestive system; (8) abdominal wall; (9) urinary system; (10) genital organs; (11) limb; and (12) other malformations. The primary outcomes of interest were overall major congenital malformations and heart defects (accounting for most malformations); other types of organ-specific malformations were considered secondary outcomes, owing to the anticipated small number of events. For exploratory analyses, we evaluated the risks of individual categories of heart defects and digestive system, as previous studies have reported that these specific malformations may be correlated with benzodiazepines [15,1719].

Covariates

We considered a broad range of covariates as potential confounders or proxies of potential confounders: maternal demographics status (e.g., age and the type of insurance), psychiatric conditions (e.g., bipolar disorder, depression/mood disorder, anxiety, and sleep disorder), maternal conditions (e.g., epilepsy/seizures, headache/migraine, diabetes, hypertension), obstetric conditions (parity, plurality), concomitant medications (e.g., antidepressants, anticonvulsants, antipsychotics, and opioid analgesics), and measures of healthcare utilization (e.g., obstetric comorbidity index [28,29], number of distinct diagnoses, and hospital admission) (S1 Table). Maternal comorbidities and concomitant medication use were measured from 6 months before the last menstrual period to the end of the first trimester. Measures of healthcare utilization were measured during the 6 months before, but not during, pregnancy to avoid these variables being affected by early awareness of pregnancy, except for obstetric comorbidity index [28,29].

Statistical analysis

Baseline characteristics of women exposed and those unexposed to benzodiazepines were compared using absolute standardized differences (aSDs; ≥0.1 indicates a significant imbalance between the 2 groups). The absolute risks (per 1,000 pregnancies), risk differences, and unadjusted relative risks (RRs) with 95% confidence intervals (CIs) were calculated for each outcome, stratified by exposure to benzodiazepines. Moreover, we calculated the population-attributable fractions (PAFs) to estimate the impact on an outcome of exposure in the study population [30]. In this study, PAF means the fractions of overall or individual congenital malformations in pregnant women attributable to benzodiazepines. We calculated PAF as (O-E)/O, where O is the observed number of outcomes and E is the expected number of outcomes under no exposure [31]. We used the propensity score (PS) fine stratification method to control for potential confounders [32]. Accordingly, the PS for exposure to benzodiazepines versus nonexposure was derived using a logistic regression model, which included all covariates without additional selection. After trimming the observations in nonoverlapping regions of the PS distribution, the exposed women were divided into 50 equally sized strata based on the PS distribution; then, unexposed women were weighted using the distribution of the exposed women for each stratum. We estimated the adjusted RR with 95% CI using a generalized linear model (log-binomial model).

We conducted 7 prespecified subgroup analyses to determine whether the risk of congenital malformations varied across exposure and maternal characteristics. Stratified analyses were performed according to the duration of action, individual benzodiazepines, and mean daily dose group with the 3 categories. Furthermore, we conducted stratified analyses by maternal age at delivery (≤35 and >35 years), multifetal pregnancy, history of epilepsy, and concomitant use of antidepressants.

Additionally, we performed diverse sensitivity analyses for all primary and secondary outcomes to evaluate the robustness of the main findings. First, we redefined the use of exposure as having filled at least 2 benzodiazepine prescriptions during the first trimester. Second, we redefined the outcome definition as the presence of ≥2 diagnoses of congenital malformations. Third, we restricted the study cohort to those who had underlying comorbidities related to the indication for benzodiazepines (bipolar disorder, depression/mood disorder, anxiety, sleep disorder, and gastrointestinal disease) to mitigate confounding by indication. Fourth, we restricted the cohort to nulliparous women to account for intraindividual correlations that might arise from repeated measurements of the same women. Fifth, we conducted a negative control analysis by comparing negative exposure control (defined as pregnancies exposed to benzodiazepines between 180 days and 90 days before the last menstrual period, which is not an etiologically relevant window for congenital malformations) with the reference group in the main analysis (pregnancies not exposed to benzodiazepines in the first trimester). If the main finding is subject to residual confounding, we can expect a non-null result from the negative control analysis. Sixth, for outcomes that presented an increased risk, we used a rule-out approach to explore the impact of unmeasured confounders (e.g., maternal smoking status) (S1 Appendix). Lastly, we conducted a quantitative bias analysis based on the probabilistic method to address the impact of selection bias (S2 Appendix). All analyses were conducted according to the prespecified analysis plan (S1 Protocol). Statistical analyses were performed using SAS version 9.4 (SAS Institute, Cary, NC, USA).

Results

Our cohort consisted of 3,094,227 pregnancies, of which 40,846 (1.3%) were prescribed at least one benzodiazepine during the first trimester (mean [standard deviation (SD)] age, 32.4 [4.1] years) (Table 1). The benzodiazepine-exposed group had higher medical conditions and concomitant medication use for both psychiatric-related and unrelated conditions (e.g., 12.6% exposed versus 0.7% unexposed had an anxiety disorder; 16.9% exposed versus 1.4% unexposed received an antidepressant) than the unexposed group. All cohort characteristics were well balanced between the exposed and unexposed groups after PS adjustment, with an aSD <0.1.

thumbnail
Table 1. Cohort characteristics of pregnancies with and without benzodiazepine exposure during the first trimester.

https://doi.org/10.1371/journal.pmed.1003945.t001

The absolute risk difference for overall malformations was 13.9 per 1,000 pregnancies (65.3 versus 51.4 per 1,000 in the exposed and unexposed groups, respectively) and that for heart defects was 11.5 per 1,000 (38.9 versus 27.4 per 1,000) (Fig 1). The PAFs of overall malformations and heart defects were 0.36% and 0.55%, respectively. The unadjusted RRs increased for overall malformations, heart defects, digestive system defects, abdominal wall defects, urinary system defects, genital defects, and other malformations. After adjustment for potential confounders, the RR estimates shifted substantially toward a null value; however, the risk for overall malformations and heart defects, although small, remained significantly elevated (adjusted RR 1.09 [95% CI 1.05 to 1.13, p < 0.001] and 1.15 [1.10 to 1.21, p < 0.001], respectively). Of individual heart defects, significant associations were found in cardiac septal defects (adjusted RR 1.13, 95% CI 1.07 to 1.20, p < 0.001) and defects of the great arteries (1.28, 95% CI 1.16 to 1.42, p < 0.001) (Fig 2).

thumbnail
Fig 1. Absolute and relative risks of congenital malformations in infants following maternal exposure to benzodiazepines during the first trimester.

CI, confidence interval; PAF, population attributable fraction; PS, propensity score; *RD1000, risk difference per 1,000 births.

https://doi.org/10.1371/journal.pmed.1003945.g001

thumbnail
Fig 2. Risks of individual categories of heart and digestive system defects in infants following maternal exposure to benzodiazepines during the first trimester.

CI, confidence interval; PAF, population attributable fraction; PS, propensity score; *RD1000, risk difference per 1,000 births.

https://doi.org/10.1371/journal.pmed.1003945.g002

We observed that the risks for overall malformations were comparable between short- and long-acting benzodiazepines (adjusted RR 1.09 [95% CI 1.03 to 1.14, p < 0.001] versus 1.07 [95% CI 1.01 to 1.13, p = 0.014]) (Figs 3, S2 and S3). The risk of overall malformations, although small but significantly, increased with midazolam, diazepam, flunitrazepam, and chlordiazepoxide administration; the risk of heart defects increased with midazolam, alprazolam, diazepam, flunitrazepam, and ethyl loflazepate administration (range of adjusted RRs for the primary outcomes: 1.08 to 2.43). The increased risk for overall malformations and heart defects was highest in the group of a mean daily lorazepam-equivalent dose of >2.5 mg/day (adjusted RR 1.26 [95% CI 1.17 to 1.36, p < 0.001] and 1.31 [95% CI 1.19 to 1.45, p < 0.001], respectively) (Fig 4). In addition, the RR for both primary outcomes were higher among women aged >35 years, those with multifetal pregnancy, and those without a history of epilepsy than in their counterparts. The RR for heart defects was higher among women who used both benzodiazepines and antidepressants than in women unexposed to these agents.

thumbnail
Fig 3. Risks of congenital malformations in infants according to individual benzodiazepine exposure during the first trimester.

BZD, benzodiazepine; CI, confidence interval; RR, relative risk; PS, propensity score. Clorazepate was not analyzed owing to the small sample size (n = 14).

https://doi.org/10.1371/journal.pmed.1003945.g003

thumbnail
Fig 4. Risks of congenital malformations in infants following maternal exposure to benzodiazepines in the first trimester: subgroup analyses.

CI, confidence interval; PS, propensity score. *Mean daily dose based on the lorazepam-equivalent dose.

https://doi.org/10.1371/journal.pmed.1003945.g004

Our main findings remained largely consistent in all sensitivity analyses when redefining exposure and outcome and restricting to pregnancies with benzodiazepine-related underlying diseases and nulliparous pregnancies, as well as the negative control analysis (Figs 5, S4, S5 and S6).

thumbnail
Fig 5. Risks of congenital malformations in infants following maternal exposure to benzodiazepines in the first trimester: sensitivity analyses.

CI, confidence interval; PS, propensity score.

https://doi.org/10.1371/journal.pmed.1003945.g005

Discussion

Main findings

In this nationwide cohort study of approximately 3.1 million pregnancies, first-trimester benzodiazepine use was associated with a small increased risk of overall malformations, particularly heart defects. The risk of primary outcomes increased with a higher mean daily dose of benzodiazepines (>2.5 mg/day of lorazepam-equivalent dose), suggesting a dose–response relationship. Although these risks were similar between short- and long-acting benzodiazepines, a small but significant increase in risk was detected with several specific agents. For the other 11 organ-specific malformation types, we found no significant increase in the risk associated with first-trimester benzodiazepine exposure. Based on the upper limit of the 95% CI from adjusted estimates, the maximum observed risk was 63% for the abdominal wall defects (RR 1.08, 95% CI 0.72 to 1.63); thus, we could rule out the possibility of a more than 63% increase in the risk for all organ-specific malformations. Our findings were consistent across diverse sensitivity analyses, and the null finding in the negative control analysis strengthened the suggestion that the estimate is unlikely to be due to residual confounding.

A recent study has reported an increased risk of spontaneous abortion associated with benzodiazepine exposure in early pregnancy (adjusted odds ratio 1.85; 95% CI 1.61 to 2.12), indicating the gestational teratogenicity of benzodiazepines [33]. In the latest meta-analysis, which included 8 cohort studies, odds ratios of 1.13 (95% CI 0.99 to 1.30) and 1.27 (95% CI 0.98 to 1.65) for overall malformations and heart defects, respectively, were reported (5,195 exposed pregnancies) [11]; the magnitude of the observed effect was similar to that in our study. Although the meta-analysis concluded that no association existed between benzodiazepine use in pregnancy and congenital malformations with null findings, the lower limit of the 95% CI of estimates was close to 1. Moreover, most previous studies had limited power and were insufficient to adequately assess the potential risk of organ-specific malformations. Notably, our study had a much larger cohort of pregnancies (40,846 exposed pregnancies) than any study published to date, thus expanding on previous findings by providing more precise estimates, as well as controlling numerous potential confounders.

An increased risk of oral cleft, which was reported in early studies evaluating in utero benzodiazepine exposure [34,35], was not observed in our study and has also not been observed in recent studies [15,36]. Three previous studies have suggested an increased risk of digestive system defects associated with benzodiazepines [1719], but the risk was not confirmed in our study. Although a study using the Swedish birth registry has reported an increased risk of pyloric stenosis (odds ratio 3.31, 95% CI 1.53 to 7.84), the result was based on 8 events among 2,537 infants exposed to benzodiazepines [19]. In the present study, only one case of the upper alimentary tract defect was reported among 40,846 pregnancies exposed to benzodiazepines.

Furthermore, we observed a small but significantly increased risk of congenital malformations in the high-dose group (the mean daily lorazepam-equivalent dose >2.5 mg/day, which is higher than the daily dose defined by the World Health Organization) [37]. According to literature, neonates can slowly metabolize small doses of benzodiazepine; however, the drug persists at pharmacologically active concentrations for at least 1 week when high doses are administered to the mother [16,38]. The benzodiazepine residues that exceed the fetal metabolic capacity might impact the incidence of congenital malformations. Furthermore, we observed increased risks for overall or cardiac malformations with several specific benzodiazepines. The increased risk observed with some specific agents should be prudently interpreted, as no clear biological mechanism can explain these risks and chance finding cannot be excluded. Accordingly, this finding could be construed as a safety signal that should be carefully monitored in future studies. Moreover, the risk of heart defects associated with alprazolam use has been described previously [14,16], reported a nearly doubled risk (odds ratio 2.43, 95% CI 1.42 to 4.15); however, the study failed to consider the indication of use and psychiatric comorbidities, which may result in increased effect size in the exposure versus nonexposure groups [19,39].

Strengths and limitations

To our knowledge, this is the largest cohort study evaluating the association between first-trimester benzodiazepine exposure and congenital malformations. The large sample size allowed us to evaluate the risk of rare malformations stratified by individual agents and specific dose groups. Additionally, we used a nationwide database representing the entire population of South Korea, which allowed us to generate generalizable real-world evidence. Our study also had potential limitations. First, misclassification of exposure is possible. Thus, we redefined exposure as at least 2 filled benzodiazepine prescriptions during the first trimester on the assumption that if a woman refilled the prescription, she probably took them. Second, outcome misclassification is possible. Therefore, we conducted a sensitivity analysis requiring at least 2 diagnoses, which increased the likelihood that outcomes reflect the actual occurrence of congenital malformations, and the results were consistent with the main findings. Third, our results could be influenced by unmeasured confounders despite the adjustment for abundant confounders. To address this concern, we conducted a negative control analysis, and this analysis revealed no association, suggesting that our main finding was not attributed to residual confounding. In addition, we used the rule-out approach to explore the effect of unmeasured confounders; the result indicated that it was unlikely that the unmeasured confounder would explain the observed association (S1 Appendix). Fourth, our findings could be affected by confounding by indication, as we used the unexposed group as the reference group. However, the results were consistent with the main findings when restricting the study cohort to women who had underlying disease related to the indication for benzodiazepine. Fifth, our study cohort included live births only, which might lead to selection bias because severe malformations that resulted in pregnancy terminations would be missed. Therefore, we conducted quantitative bias analysis, and the corrected RR for the primary outcomes remained below 1.4, under the most extreme scenario, suggesting that the effect of such selection bias is minimal (S2 Appendix).

Conclusions

In this nationwide cohort study, benzodiazepine use during the first trimester was associated with a small increased risk of overall malformations and heart defects, particularly in the high daily dose group (at doses higher than the usual daily dose). However, the absolute risks and population attributable fractions were modest. Our findings suggest that the benefits of benzodiazepines for their major indications must be considered despite the potential risks. Nonetheless, to minimize the potential risk, alternative nonpharmacological strategies could be considered for managing anxiety and insomnia during pregnancy; if benzodiazepines are necessary, the lowest effective dosage should be prescribed during early pregnancy.

Supporting information

S1 STROBE Checklist. STROBE, Strengthening the Reporting of Observational Studies in Epidemiology.

https://doi.org/10.1371/journal.pmed.1003945.s001

(DOC)

S1 Table. Codes used to define the inclusion/exclusion criteria, exposures, outcomes of interest, maternal comorbidities, and concomitant medications.

https://doi.org/10.1371/journal.pmed.1003945.s002

(DOCX)

S2 Table. Selection of the study cohort using the HIRA database between 1 July 2009 and 31 December 2019.

https://doi.org/10.1371/journal.pmed.1003945.s003

(DOCX)

S1 Fig. Frequency of pregnancies exposed to benzodiazepines during the first trimester in South Korea between 2011 and 2018.

https://doi.org/10.1371/journal.pmed.1003945.s004

(DOCX)

S2 Fig. Absolute and relative risks of overall congenital malformations in infants according to individual benzodiazepine exposure during the first trimester.

https://doi.org/10.1371/journal.pmed.1003945.s005

(DOCX)

S3 Fig. Absolute and relative risks of congenital heart defects in infants according to individual benzodiazepine exposure during the first trimester.

https://doi.org/10.1371/journal.pmed.1003945.s006

(DOCX)

S4 Fig. Risks of congenital malformations in infants following maternal exposure to benzodiazepines during the first trimester: sensitivity analyses.

https://doi.org/10.1371/journal.pmed.1003945.s007

(DOCX)

S5 Fig. Risks of congenital malformations in infants following maternal exposure to benzodiazepines during the first trimester: sensitivity analyses II.

https://doi.org/10.1371/journal.pmed.1003945.s008

(DOCX)

S6 Fig. Risks of congenital malformations in infants following maternal exposure to benzodiazepines during the first trimester: sensitivity analyses III.

https://doi.org/10.1371/journal.pmed.1003945.s009

(DOCX)

S1 Appendix. Sensitivity analysis of unmeasured confounders (rule-out approach).

https://doi.org/10.1371/journal.pmed.1003945.s010

(DOCX)

S2 Appendix. Potential effect of including live births only.

https://doi.org/10.1371/journal.pmed.1003945.s011

(DOCX)

References

  1. 1. Gentile S. Anxiety and sleep disorders, psychopharmacology, and pregnancy. Psychopharmacology and Pregnancy. Springer; 2014. p. 87–102.
  2. 2. Goodman JH, Chenausky KL, Freeman MP. Anxiety disorders during pregnancy: a systematic review. J Clin Psychiatry. 2014;75(10):e1153–84. pmid:25373126.
  3. 3. Andrade SE, Gurwitz JH, Davis RL, Chan KA, Finkelstein JA, Fortman K, et al. Prescription drug use in pregnancy. Am J Obstet Gynecol. 2004;191(2):398–407. pmid:15343213.
  4. 4. Bais B, Molenaar NM, Bijma HH, Hoogendijk WJG, Mulder CL, Luik AI, et al. Prevalence of benzodiazepines and benzodiazepine-related drugs exposure before, during and after pregnancy: A systematic review and meta-analysis. J Affect Disord. 2020;269:18–27. pmid:32217339.
  5. 5. Dell’Osso B, Albert U, Atti AR, Carmassi C, Carrà G, Cosci F, et al. Bridging the gap between education and appropriate use of benzodiazepines in psychiatric clinical practice. Neuropsychiatr Dis Treat. 2015;11:1885–909. pmid:26257524.
  6. 6. Kanto JH. Use of benzodiazepines during pregnancy, labour and lactation, with particular reference to pharmacokinetic considerations. Drugs. 1982;23(5):354–80. pmid:6124415.
  7. 7. Kanto J, Erkkola R, Sellman R. Accumulation of diazepam and N-demethyldiazepam in the fetal blood during the labour. Ann Clin Res. 1973;5(6):375–9. pmid:4792202.
  8. 8. Lee DH, Kang SK, Lee RH, Ryu JM, Park HY, Choi HS, et al. Effects of peripheral benzodiazepine receptor ligands on proliferation and differentiation of human mesenchymal stem cells. J Cell Physiol. 2004;198(1):91–9. pmid:14584048.
  9. 9. National Institute for Health and Care Excellence: Clinical Guidelines. Antenatal and postnatal mental health: clinical management and service guidance. London: National Institute for Health and Care Excellence (UK); 2014.
  10. 10. Content and format of labeling for human prescription drug and biological products; requirements for pregnancy and lactation labeling. Final rule. Fed Regist. 2014;79(233):72063–103. pmid:25509060.
  11. 11. Grigoriadis S, Graves L, Peer M, Mamisashvili L, Dennis CL, Vigod SN, et al. Benzodiazepine Use During Pregnancy Alone or in Combination With an Antidepressant and Congenital Malformations: Systematic Review and Meta-Analysis. J Clin Psychiatry. 2019;80(4). pmid:31294935.
  12. 12. Enato E, Moretti M, Koren G. The fetal safety of benzodiazepines: an updated meta-analysis. J Obstet Gynaecol Can. 2011;33(1):46–8. pmid:21272436.
  13. 13. Bellantuono C, Martellini M, Orsolini L. Benzodiazepines and Z-Drugs in Pregnancy. Perinatal Psychopharmacology. Springer; 2019. p. 203–13.
  14. 14. Bellantuono C, Tofani S, Di Sciascio G, Santone G. Benzodiazepine exposure in pregnancy and risk of major malformations: a critical overview. Gen Hosp Psychiatry. 2013;35(1):3–8. pmid:23044244.
  15. 15. Dolovich LR, Addis A, Vaillancourt JM, Power JD, Koren G, Einarson TR. Benzodiazepine use in pregnancy and major malformations or oral cleft: meta-analysis of cohort and case-control studies. BMJ. 1998;317(7162):839–43. pmid:9748174
  16. 16. Iqbal MM, Sobhan T, Ryals T. Effects of commonly used benzodiazepines on the fetus, the neonate, and the nursing infant. Psychiatr Serv. 2002;53(1):39–49. pmid:11773648.
  17. 17. Bonnot O, Vollset SE, Godet PF, d’Amato T, Dalery J, Robert E. In utero exposure to benzodiazepine. Is there a risk for anal atresia with lorazepam? Encéphale. 2003;29(6):553–9. pmid:15029090.
  18. 18. Wikner BN, Stiller CO, Bergman U, Asker C, Källén B. Use of benzodiazepines and benzodiazepine receptor agonists during pregnancy: neonatal outcome and congenital malformations. Pharmacoepidemiol Drug Saf. 2007;16(11):1203–10. pmid:17894421.
  19. 19. Källén B, Borg N, Reis M. The use of central nervous system active drugs during pregnancy. Pharmaceuticals (Basel). 2013;6(10):1221–86. pmid:24275849.
  20. 20. Ban L, West J, Gibson JE, Fiaschi L, Sokal R, Doyle P, et al. First trimester exposure to anxiolytic and hypnotic drugs and the risks of major congenital anomalies: a United Kingdom population-based cohort study. PLoS ONE. 2014;9(6):e100996. pmid:24963627.
  21. 21. Vgontzas AN, Kales A, Bixler EO. Benzodiazepine side effects: role of pharmacokinetics and pharmacodynamics. Pharmacology. 1995;51(4):205–23. pmid:8577815.
  22. 22. Park E, Jang S, Jeon S, Lee S, Lee J, Choi D. Report of the evaluation for validity of discharged diagnoses in Korean Health Insurance database. Health Insurance Review and Assessment Service; 2017.
  23. 23. Margulis AV, Setoguchi S, Mittleman MA, Glynn RJ, Dormuth CR, Hernández-Díaz S. Algorithms to estimate the beginning of pregnancy in administrative databases. Pharmacoepidemiol Drug Saf. 2013;22(1):16–24. pmid:22550030.
  24. 24. Schaefer C, Peters PW, Miller RK. Drugs during pregnancy and lactation: treatment options and risk assessment: Academic Press; 2014.
  25. 25. Drug and Alcohol Services South Australia. Benzodiazepines: Information for GPs. In: Australia GoS, editor. Updated: 14 Mar 2021.
  26. 26. Ashton HC. Benzodiazepine equivalence table. 2007. Available from: https://www.benzo.org.uk/bzequiv.htm.
  27. 27. EUROCAT. EUROCAT Guide 1.4: Instruction for the registration of congenital anomalies. EUROCAT Central Registry, University of Ulster; 2013.
  28. 28. Bateman BT, Mhyre JM, Hernandez-Diaz S, Huybrechts KF, Fischer MA, Creanga AA, et al. Development of a comorbidity index for use in obstetric patients. Obstet Gynecol. 2013;122(5):957–65. pmid:24104771.
  29. 29. Metcalfe A, Lix LM, Johnson JA, Currie G, Lyon AW, Bernier F, et al. Validation of an obstetric comorbidity index in an external population. BJOG. 2015;122(13):1748–55. pmid:25559311.
  30. 30. Mansournia MA, Altman DG. Population attributable fraction. BMJ. 2018;360:k757. pmid:29472187.
  31. 31. Lin CK, Chen ST. Estimation and application of population attributable fraction in ecological studies. Environ Health. 2019;18(1):52. pmid:31196209.
  32. 32. Desai RJ, Rothman KJ, Bateman BT, Hernandez-Diaz S, Huybrechts KF. A Propensity-score-based Fine Stratification Approach for Confounding Adjustment When Exposure Is Infrequent. Epidemiology. 2017;28(2):249–57. pmid:27922533.
  33. 33. Sheehy O, Zhao JP, Bérard A. Association Between Incident Exposure to Benzodiazepines in Early Pregnancy and Risk of Spontaneous Abortion. JAMA Psychiat. 2019;76(9):948–57. pmid:31090881.
  34. 34. Saxén I. Associations between oral clefts and drugs taken during pregnancy. Int J Epidemiol. 1975;4(1):37–44. pmid:1116890.
  35. 35. Safra MJ, Oakley GP Jr. Association between cleft lip with or without cleft palate and prenatal exposure to diazepam. Lancet. 1975;2(7933):478–80. pmid:51287.
  36. 36. Rosenberg L, Mitchell AA, Parsells JL, Pashayan H, Louik C, Shapiro S. Lack of relation of oral clefts to diazepam use during pregnancy. N Engl J Med. 1983;309(21):1282–5. pmid:6633586.
  37. 37. World Health Organisation (WHO) Collaborating Centre for Drug Statistics Methodology. The anatomical therapeutic chemical (ATC) classification system and defined daily dose (DDD) index 2021. Available from: https://www.whocc.no/atc_ddd_index/.
  38. 38. Mandelli M, Morselli PL, Nordio S, Pardi G, Principi N, Sereni F, et al. Placental transfer to diazepam and its disposition in the newborn. Clin Pharmacol Ther. 1975;17(5):564–72. pmid:1126114.
  39. 39. Kyriacou DN, Lewis RJ. Confounding by Indication in Clinical Research. JAMA. 2016;316(17):1818–9. pmid:27802529.