Skip to main content
  • Research article
  • Open access
  • Published:

Use of smoking cessation pharmacotherapies during pregnancy is not associated with increased risk of adverse pregnancy outcomes: a population-based cohort study



Varenicline, bupropion and nicotine replacement therapy (NRT) are three effective pharmacotherapies for smoking cessation, but data about their safety in pregnancy are limited. We assessed the risk of adverse perinatal outcomes and major congenital anomalies associated with the use of these therapies in pregnancy in Australia.


Perinatal data for 1,017,731 deliveries (2004 to 2012) in New South Wales and Western Australia were linked to pharmaceutical dispensing, hospital admission and death records. We identified 97,875 women who smoked during pregnancy; of those, 233, 330 and 1057 were exposed to bupropion, NRT and varenicline in pregnancy, respectively. Propensity scores were used to match exposed women to those who were unexposed to any smoking therapy (1:10 ratio). Propensity scores and gestational age at exposure were used to match varenicline-exposed to NRT-exposed women (1:1 ratio). Time-dependent Cox proportional hazards models estimated hazard ratios (HR) with 95% confidence intervals (95% CI) for any adverse perinatal event (a composite of 10 unfavourable maternal and neonatal outcomes) and any major congenital anomaly.


The risk of any adverse perinatal event was not significantly different between bupropion-exposed and unexposed women (39.2% versus 39.3%, HR 0.93, 95% CI 0.73–1.19) and between NRT-exposed and unexposed women (44.8% vs 46.3%, HR 1.02, 95% CI 0.84–1.23), but it was significantly lower in women exposed to varenicline (36.9% vs 40.1%, HR 0.86, 95% CI 0.77–0.97). Varenicline-exposed infants were less likely than unexposed infants to be born premature (6.5% vs 8.9%, HR 0.72, 95% CI 0.56–0.92), be small for gestational age (11.4% vs 15.4%, HR 0.68, 95% CI 0.56–0.83) and have severe neonatal complications (6.6% vs 8.2%, HR 0.74, 95% CI 0.57–0.96). Among infants exposed to varenicline in the first trimester, 2.9% had a major congenital anomaly (3.5% in unexposed infants, HR 0.91, 95% CI 0.72–1.15). Varenicline-exposed women were less likely than NRT-exposed women to have an adverse perinatal event (38.7% vs 51.4%, HR 0.58, 95% CI 0.33–1.05).


Pregnancy exposure to smoking cessation pharmacotherapies does not appear to be associated with an increased risk of adverse birth outcomes. Lower risk of adverse birth outcomes in varenicline-exposed pregnancies is inconsistent with recommendations that NRT be used in preference to varenicline.

Peer Review reports


Smoking cessation in pregnancy reduces the risk of perinatal adverse outcomes and has long-term maternal and child health benefits [1, 2]. However, a substantial proportion of women continue to smoke throughout gestation [3] and even in subsequent pregnancies despite having prior poor birth outcomes [4]. This illustrates the difficulty of quitting for women who have high levels of nicotine dependence [5], and highlights the need for effective assistance. Pharmacotherapies for smoking cessation including varenicline, bupropion and nicotine replacement therapy (NRT) are more effective than behavioural cessation interventions in non-pregnant adults [6, 7]; however, little is known about their safety and efficacy in pregnancy [1, 8,9,10].

Research evidence about pregnancy outcomes associated with varenicline exposure is limited to uncontrolled and small studies, including a case series for 24 women [11, 12] and 89 other cases identified from teratology counselling and surveillance services [13]. Available data regarding pregnancy outcomes of bupropion for smoking cessation (i.e. not for depression) are limited to two small pilot trials (n = 35) [14, 15] and two small observational studies (n < 140) [16, 17]. Given the lack of robust evidence regarding their safety, varenicline and bupropion have no therapeutic indications for smoking cessation in pregnancy [2, 9, 10]. For NRT on the other hand, there have been trials [1, 18] and large-scale observational studies [1, 19,20,21]; however, neither have been conclusive [1, 18,19,20,21]. This is most likely due to poor adherence to NRT [1] and heterogeneity in the NRT products (e.g. patches, gums, lozenges, spray) under investigation [19,20,21]. Transdermal NRT releases continuous doses of nicotine, thus is potentially associated with greater risk of harm than intermittent doses offered by oral forms. Clinical guidelines support the use of NRT during pregnancy only when the expected benefits outweigh the potential risks [2, 9, 10].

Although NRT clinical trials have been possible, based on the assumption that NRT is safer than continued smoking in pregnancy, this has not been the case for varenicline or bupropion. In the absence of sufficient evidence, we established the Smoking MUMS (Maternal Use of Medications and Safety) Study [22], a population-based cohort study focused on smoking cessation pharmacotherapies during pregnancy. In this Australian cohort of women who smoked during pregnancy, 3.6% used a smoking cessation therapy, mostly varenicline (1.8%) and NRT (1.7%) [23]. The current paper aimed to:

  1. 1)

    Compare the risk of adverse perinatal outcomes and major congenital anomalies between pregnancies exposed to these pharmacotherapies and pregnancies exposed to smoking but no pharmacotherapy and

  2. 2)

    Compare the risk of these outcomes between pregnancies that were exposed to different therapies.


Study data source

The Smoking MUMS Study is a cohort design based on all pregnancies that resulted in a birth (gestational age ≥ 20 weeks or birthweight ≥ 400 g) in two Australian states, New South Wales (NSW) and Western Australia (WA), between 2003 and 2012. The study protocol, data sources and data preparation have been described elsewhere [22, 24]. The current analyses used linked records from four data sources including perinatal data (deliveries in 2003–2012), dispensing data for pharmaceuticals subsidised through the Pharmaceutical Benefits Scheme (PBS, 2003–2013), hospital admissions (2001–2013) and deaths (2003–2014). Under Australian universal health care system, eligible residents had access to subsidised prescriptions of bupropion (for smoking cessation only) from February 2001 and varenicline from January 2008. Subsidy for NRT patches commenced in December 2008, initially only for Aboriginal and Torres Strait Islander people, and then for the general population since January 2011. At time of conducting this study, oral forms of NRT were not subsidised in Australia. All forms of NRT, including patches, were also available over the counter whilst bupropion and varenicline were only available on prescription.

Study population

We identified a base cohort of pregnant women (conception between 1 January 2004 and 9 April 2012) who smoked during pregnancy. Conception date was estimated using the equation: date of conception = date of delivery − gestational age at delivery × 7 + 14 days [22]. The beginning and end points of cohort entry ensure that dispensing data were available for at least 1 year prior to conception and deliveries in 2012 did not disproportionately include pregnancies with gestation shorter than 40 weeks.

Maternal smoking status in pregnancy was derived from either perinatal or maternal hospital admission data [25]. During the study period, there were changes in maternal smoking items in the perinatal data collections (2010 in WA and 2011 in NSW). Therefore, for those who delivered before 2010 in WA and before 2011 in NSW, a woman was identified as having smoked during pregnancy if the response was Yes to the item about whether a woman smoked during pregnancy. For those who delivered in 2010 onwards in WA, smoking status was defined as Yes if the reported quantity of tobacco cigarettes smoked each day during the first or the second half of pregnancy was greater than 0. For those who delivered in NSW in 2011 onwards, smoking status was based on a Yes response to items indicating whether a woman smoked in the first or the second half of the pregnancy. Based on the hospital admission corresponding to the delivery [25], a woman was identified as having smoked during pregnancy if a Z72.0 code (i.e. use of tobacco in the last month) [26] was recorded in any diagnosis field. In addition, some women who were not identified as having smoked in pregnancy but had dispensing of a smoking cessation therapy during pregnancy were reclassified as smokers according to a published algorithm [23].

From the base cohort above, we identified women who were exposed to bupropion, varenicline and NRT during pregnancy and women who had never been exposed to any of these medicines in pregnancy. As aforementioned, subsidy for these medicines commenced at different time points. Accordingly, for comparisons involving bupropion, varenicline and NRT, the study periods started from 1 January 2004, 1 January 2008 and 1 January 2009, respectively. When comparing outcomes of different therapies (aim 2), we required exposure to occur over the same calendar time to avoid confounding by the underlying temporal variation in pregnancy and labour care [27, 28]. As only 30 pregnancies were exposed to bupropion between 2009 and 2012 (when both varenicline and NRT were available), comparisons between bupropion and other therapies were not conducted. In accordance with ethical approvals for this study, congenital anomalies were examined among infants born in NSW only. Cell sizes were such that this outcome could be examined among only varenicline-exposed pregnancies relative to pregnancies exposed to smoking but no pharmacotherapy (Fig. 1).

Fig. 1
figure 1

Flowchart of pregnancies included in the analyses

Exposure to smoking cessation pharmacotherapy

Dispensing data included records of all subsidised varenicline (Anatomical Therapeutic Chemical [ATC] code N07BA03), NRT patches (N07BA01) and bupropion (code change from N07BA02 to N06AX12) [29]. The recommended course for non-pregnant adults is at least 7 weeks for bupropion (150 mg once daily for the first 3 days then 150 mg twice daily), at least 8 weeks for NRT (one patch daily) and 12 weeks for varenicline (0.5 mg once daily for the first 3 days then 1 mg once daily); one course generally involves two or more prescriptions [2]. Records of these therapies dispensed in the period from 100 days pre-conception to date of delivery were identified; the earliest dispensing in this period was referred to as the index dispensing. Days covered by each dispensing was estimated by dividing the quantity dispensed by recommended daily dosage. A pregnancy was identified as exposed if there were one or more dispensings of the therapy in the gestation period (i.e. conception to delivery) or pre-conception dispensings were sufficient to last into the gestation period (i.e. date of the index dispensing + days covered by pre-conception dispensings ≥ date of conception). Date of exposure was defined as either (i) date of the index dispensing if dispensed after conception or (ii) date of conception if pre-conception dispensings were sufficient to last into the gestation period or (iii) date of the first post-conception dispensing if pre-conception dispensings were insufficient to last into the gestation period. Gestational age at exposure (expressed as weeks, the whole number) was calculated using the formula [gestational age at delivery − (date of delivery − date of exposure)/7].

Study outcomes

The two primary outcomes of the study were whether a woman or neonate experienced any adverse perinatal event at birth and whether the neonate had any major congenital anomaly. Any adverse perinatal event was a composite of 10 individual birth outcomes, including preterm birth (< 37 weeks, medically indicated or spontaneous), small for gestational age (SGA, birthweight < 10th percentile sex- and gestational age-specific) [30], Apgar score at 5 min < 7, admission to neonatal special care (NSC), severe neonatal morbidity complications [31], emergency caesarean section, severe maternal morbidity complications [32], preterm premature rupture of membranes (PPROM), placental abruption and perinatal death (stillbirth or 28-day neonatal death). These outcome variables were derived from the perinatal record, hospital record relating to the mother’s delivery, hospital record relating to the baby’s birth and mortality data. Outcomes including SGA, NSC admission, Apgar score and severe neonatal morbidity complications were assessed among live births only. The second primary outcome was a diagnosis of any major congenital anomaly recorded in hospital admissions occurring within 18 months from birth, among live-born babies in NSW. A detailed description of the study outcome variables is presented in Additional file 1.


Exclusion criteria included multiple births, conception within 6 months from the immediate preceding delivery [33, 34], interstate residents and overseas visitors (likely incomplete capture of hospital admission and dispensing data), use of multiple smoking therapies in the same pregnancy, use of potentially teratogenic medications during pregnancy (category D and X, according to the Australian Therapeutic Goods Administration classification system) [35], likely data errors (birthweight < 1000 g whilst gestational age > 38 weeks, birthweight > 5500 g whilst gestational age < 37 week, based on the Australian birthweight chart) [30] and missing data (mostly due to geographical area of residence and Apgar score not being recorded). In the analysis of the major congenital anomaly outcome, pregnancies linked to anomalies due to chromosomal malformations, viral infection (e.g. cytomegalovirus, rubella) and known exogenous causes were also excluded. The numbers of excluded pregnancies according to exclusion criteria are presented in Fig. 1 and Additional file 1: Table S1.4.

Statistical analyses

For aim 1 comparisons, among unexposed women, 24.5% had two or more pregnancies; their first pregnancy was selected. In the exposed groups, a small number of women (n < 5) had two pregnancies exposed to the same therapy; their first exposed pregnancy was selected. In the analyses for aim 2, there was no woman with a pregnancy exposed to both varenicline and NRT.

Propensity scores were calculated using logistic regression models in which the outcome variable was exposure to a smoking cessation therapy and explanatory variables included state of birth, year of conception, maternal age at conception, Aboriginal and/or Torres Strait Islander status, country of birth, marital status, quintiles of socio-economic disadvantage scores associated with the residential area [36], geographical remoteness of the residential area [37], private health insurance, parity, previous caesarean section, number of hospital admissions in the year prior to conception and pre-existing maternal morbidities (mental health, chronic airway conditions, diabetes, hypertension, cardiovascular diseases, epilepsy, chronic renal diseases, thyroid disorders, substance and alcohol diagnoses, anaemia and coagulation disorders, the use of steroids, non-steroid anti-inflammatory drugs and medications for gastro-oesophageal reflux diseases). For the analysis of congenital anomalies, the propensity score model included an additional variable indicating whether the mother had previously had a child with a major congenital anomaly. See Additional file 1 for a comprehensive description of variables included in propensity score models and how they were derived from the source data.

For analyses addressing aim 1, each exposed pregnancy was matched to ten unexposed pregnancies on propensity score, using a greedy five- to one-digit algorithm (matching on the five digits of propensity score first, then four digits if there was no five-digit match, and so on to one-digit matching, with no replacement) [38]. To address aim 2, the matching (1:1 ratio) used both propensity score and gestational age at exposure. Absolute standardised differences were calculated to assess balance in the characteristics of the comparison groups (balanced if the difference < 0.1) [39]. Cox proportional hazard modelling (discrete tier for matched data, gestational week as time scale) was used to estimate hazard ratios (HRs). In these Cox models, exposure was defined as a time-dependent variable, i.e. a pregnancy was considered unexposed until the gestational week at exposure. The window of exposure was the first trimester (gestational week at exposure < 13) when examining the congenital anomaly outcome, up to week 37 (gestational week at exposure < 37) for the composite adverse perinatal event and preterm birth, and until delivery for other individual perinatal outcomes.

For comparisons addressing aim 1, the matched samples were well-balanced on maternal characteristics (Table 1); thus, univariable Cox models were built. For aim 2 comparisons, multivariable Cox models were built to adjust for characteristics with a standardised difference > 0.1 [39]; however, only the multivariable model assessing the composite perinatal outcome converged. Therefore, for this composite perinatal outcome, both univariable and multivariable HRs were reported, and for individual perinatal outcomes, univariable HRs were presented.

Table 1 Demographic and clinical characteristics of exposed women and matched comparison group, stratified by therapy comparison. Data are numbers (percentages) and standardised differences

Given the concern that women might have ceased therapy when they became aware of their pregnancy, we conducted sensitivity analyses by restricting the analyses to women who initiated the therapy after week 4 of gestation. It was only possible to conduct sensitivity analyses for the composite adverse perinatal event outcome (see Additional file 2: Table S2.6) because of small sample sizes for the other individual perinatal outcomes. All analyses were carried out in SAS version 9.4.


In the Smoking MUMS Study, there were 1,017,731 pregnancies belonging to 686,884 women in NSW and WA (conception between 1 January 2004 and 9 April 2012). The base cohort included 103,753 women who smoked in 140,913 pregnancies; of those, 13,667 pregnancies were excluded. Following the selection of one pregnancy per woman, data analyses included 233, 1057 and 330 women who were exposed to bupropion, varenicline and NRT, respectively, and 96,255 unexposed women (Fig. 1). As presented in Fig. 2, most women had only one dispensing of a smoking cessation therapy. Therapy initiation before conception was more common for bupropion (65.1%) and varenicline (65.6%) than for NRT (20%).

Fig. 2
figure 2

Timing and duration of exposure to a smoking cessation pharmacotherapy in pregnancy. The beginning of each horizontal line indicates the date of dispensing and the length of the line represents the number of days covered by the dispensing. *LMP last menstrual period (estimated)

Table 1 shows that the matching procedures resulted in well-balanced baseline characteristics between exposed and unexposed pregnancies (see Additional file 2: Tables S2.1-S2.5 for before-matching characteristics). Nevertheless, in comparison to NRT-exposed women, varenicline-exposed women were more likely to be born overseas, had a partner and a private health insurance but were less likely to have a mental health condition and drug and alcohol disorders.

The analyses of the first primary outcome (Table 2) showed that the risk of having any adverse perinatal event was not significantly different between bupropion-exposed and unexposed pregnancies (hazard ratio [HR] 0.93, 95% confidence interval [CI] 0.73 to 1.19) or between NRT-exposed and unexposed pregnancies (HR 1.02, 95% CI 0.84 to 1.23). Compared to unexposed women, there was a lower risk of any adverse perinatal event in those who were exposed to varenicline (HR 0.86, 95% CI 0.77 to 0.97). When exposure to NRT was the reference group, exposure to varenicline was not associated with higher risk of any adverse event (multivariable HR 0.58, 95% CI 0.33 to 1.05, univariable HR 0.59, 95% CI 0.38 to 0.91). As presented in Table 3, 2.9% of infants exposed to varenicline in the first trimester had a major congenital anomaly compared to 3.5% of infants not exposed to any smoking cessation therapy (HR 0.91, 95% CI 0.72 to 1.15).

Table 2 Adverse perinatal outcomes in exposed women and the matched comparison group, stratified by therapy comparison. Data presented are numbers (percentage) of outcomes and hazard ratios (95% confidence interval)
Table 3 Major congenital anomalies in infants exposed to varenicline in the first trimester and matched comparison group. Data presented are numbers (percentage) of outcomes and hazard ratios (95% confidence interval)

Our analyses of individual birth outcomes (Table 2) also showed that varenicline-exposed women were less likely than unexposed women to have a baby who was born preterm (HR 0.72, 95% CI 0.56 to 0.92), was small for gestational age (HR 0.68, 95% CI 0.56 to 0.83), had an Apgar score at 5 min less than 7 (HR 0.47, 95% CI 0.27 to 0.81) and had severe morbidity complications (HR 0.74, 95% CI 0.57 to 0.96). Varenicline-exposed women were also less likely than NRT-exposed women to have a preterm birth (HR 0.29, 95% CI 0.12 to 0.71).

In the sensitivity analyses restricted to women who initiated therapy after week 4 of gestation, there were large reductions in the number of exposed pregnancies. Hazard ratios for the composite adverse perinatal event were similar to those of the main analyses; however, the effects of varenicline (compared to no therapy) became statistically non-significant (see Additional file 2: Table S2.6).


Principal findings

Overall, our population-based study found significant reduction in the risk of any adverse perinatal event associated with pregnancy exposure to varenicline and no increased risk of major congenital anomalies associated with first trimester exposure to varenicline. There was also no significant increase in risk of any adverse perinatal event associated with exposure to bupropion and NRT. When individual perinatal outcomes were examined, findings were most encouraging for varenicline such as a significant reduction in the risk of preterm birth, SGA and severe neonatal morbidity complications.

Comparison with other studies

This is the first controlled comparison of varenicline exposure in pregnancy; prior studies were not able to draw a causal relationship [11,12,13]. Protective effects of varenicline in this study could be potentially explained by it being the most efficacious pharmacotherapy for smoking cessation as reported by studies in non-pregnant adults [40], and the timing of its use. In our study, the majority of varenicline-exposed pregnancies had the therapy dispensed prior to conception or in the first few weeks of the first trimester, a pattern that is consistent with the recommendation that varenicline not be used in pregnancy [2, 9, 10]. Early use of varenicline in pregnancy could have resulted in early quitting in pregnancy. Prior studies have shown that earlier smoking cessation in pregnancy is associated with a greater reduction in the risk of preterm birth [41, 42] and SGA [41, 43]. Whilst the early use of varenicline may explain why pregnancies exposed to varenicline had better outcomes than pregnancies not exposed to any therapy, superior efficacy of varenicline in non-pregnant adults compared to NRT [40] may explain why pregnancies exposed to varenicline had more favourable outcomes than those exposed to NRT, because in the varenicline-NRT comparisons, pregnancies were matched on gestational age at exposure. Within the Smoking MUMS Study, another analysis is underway to compare the effectiveness of varenicline and NRT when used in pregnancy.

Our study found that exposure to varenicline in the first trimester was not associated with increased risk of any major congenital anomaly. A previous study of varenicline was limited to two cases with congenital anomalies [13]. In our study, the proportion of infants with a major congenital anomaly (exposed to varenicline or exposed to smoking but not pharmacotherapy) are consistent with the underlying Australian and international figures [44,45,46]. To date, there were no human studies that have examined the pharmacokinetics of varenicline in pregnancy. Animal studies found that maternal exposure to high dosage of varenicline was associated with low foetal weight and development toxicity in offspring [47]. Studies in rats and rabbits did not find a teratogenic effect of varenicline, even with administered dosages 23 and 50 times higher, respectively, than the maximum recommended human daily dose [47]. However, it has been long established that animal studies are seriously limited in their ability to predict human teratogenesis [48, 49].

Overall, this is the most rigorous study to assess safety of bupropion for smoking cessation in pregnancy. Previous smaller studies [14,15,16,17] were based on self-reported exposure [16, 17] and were not able to separate bupropion use for smoking cessation from its use for depression [16, 17]. Our study found a 15% increase in the risk of emergency caesarean section although this was not statistically significant. This raised questions about the biological pathway of bupropion, given that bupropion and its metabolites can cross the placenta into foetal circulation [50].

The medicinal form of nicotine prevents a foetus from being exposed to a multitude of toxic substances in cigarette smoke [2, 9]. As nicotine and cotinine (nicotine’s metabolite) pass through the placenta [21, 51], nicotine obtained through the patches could pose health risk to the foetus, e.g. disrupted development of foetal cholinergic system in the first trimester [51]. However, due to faster clearance of nicotine and cotinine in pregnancy [52], effects of NRT would be less profound than cigarette smoke. Our study found no differences in perinatal outcomes between NRT-exposed and unexposed pregnancies which are generally consistent with reports from clinical trials [1, 18] and observational studies [19, 20]. A previous study which reported an increased risk of preterm birth and low birthweight [21] failed to control for the effect of maternal smoking. Nevertheless, a clinical trial has showed that high-dose NRT patches may increase diastolic blood pressure in late pregnancy, which may potentially lead to unfavourable pregnancy outcomes [18].

Strengths and limitations of the study

To our knowledge, this is the largest study to date about the safety of varenicline in pregnancy. We used different study designs (i.e. non-user and active comparator) to examine outcomes associated with varenicline use; this addressed concerns about confounding by indication and health-seeking behaviours. The use of propensity score matching also addresses confounding by maternal characteristics [53]. Exposure was analysed in a time-dependant manner, eliminating immortal time bias [54]. Reliable recording of date of delivery and gestational age in perinatal data [55, 56] enabled accurate ascertainment of pregnancy exposure.

This study has some limitations. The study included pregnancies delivered at least 20 weeks of gestation; thus, outcomes such as pregnancy loss or termination before week 20 were not examined. Our study had inadequate statistical power to address effects of exposure to bupropion and NRT on rare perinatal adverse outcomes such as severe maternal morbidity complications, PPROM, Apgar score at 5 min < 7, placental abruption and perinatal mortality. It was not possible for the study to examine specific major congenital anomalies associated with first trimester varenicline exposure. Outcomes such as changes in mood, behaviour or thinking were not examined in this study. Although the risk of these mental health side effects associated with varenicline or bupropion use has been found to be less profound than previously thought, a risk remains in people with a history of mental illness [57]. In Australia, NRT can be purchased over the counter, and in geographically remote regions, clients of Aboriginal Health Services may receive NRT free of charge without the need for a prescription [58]. Such data were not captured in PBS dispensing data; thus, some exposed pregnancies would have been classified as unexposed. Nevertheless, we believe the extent of misclassification of exposure status would be negligible due to financial barriers such as low income among a large proportion of smoking women [59] and high out-of-pocket costs of NRT purchased over the counter [60]. Although the main analyses in this study assumed that exposed women took all of the dispensed medicine as per recommended dosages, this may not be the case and actual exposure status and periods of exposure may deviate from this assumption. A particular concern was that women might have ceased therapy when they became aware of their pregnancy, but sensitivity analyses focused on women who initiated therapy after week 4 of gestation suggest that the study findings are robust.


Given the uncertainty about the safety of varenicline and bupropion during pregnancy, these therapies are not recommended in pregnant women. As evidence regarding the safety of NRT during pregnancy is inconsistent, clinical guidelines place the onus on the physician to decide whether the potential benefits of NRT use during pregnancy outweigh the risk of harm [2, 9, 10]. As a result, women and healthcare providers are currently in a bind when deciding whether a smoking cessation pharmacotherapy should be used, and if so, which one. Our study showed that varenicline is being used in pregnancy and its benefit is likely to outweigh the harm. This could help assess whether experimental studies might be ethical. Further multi-jurisdictional collaborations are needed for more robust evidence, which could allow investigations of outcomes such as miscarriage, terminations, specific malformations, maternal cardiovascular and neuropsychiatric events, and long-term outcomes for the babies given concerns about neurological development [47].

Availability of data and materials

The data sets were constructed with the permission of each of the source data custodians and with specific ethical approvals. The authors do not have permission to share patient-level data because of the highly confidential nature of the data. Permission to access to the data is restricted to researchers named and approved by relevant Human Research Ethics Committees.



Australian Classification of Health Interventions, Eighth Edition


Anatomical Therapeutic Chemical


Hazard ratio


International Statistical Classification of Diseases and Related Health Problems, Tenth Revision, Australian Modification


Nicotine replacement therapy


Neonatal special care


New South Wales


Pharmaceutical Benefits Scheme


Preterm premature rupture of membranes


Small for gestational age


Western Australia


  1. Coleman T, Chamberlain C, Davey MA, Cooper SE, Leonardi-Bee J. Pharmacological interventions for promoting smoking cessation during pregnancy. Cochrane Database Syst Rev. 2015;12.

  2. The Royal Australian College of General Practitioners. Supporting smoking cessation: a guide for health professionals. Melbourne: The Royal Australian College of General Practitioners; 2014.

    Google Scholar 

  3. Schneider S, Huy C, Schutz J, Diehl K. Smoking cessation during pregnancy: a systematic literature review. Drug Alcohol Rev. 2010;29(1):81–90.

    Article  PubMed  Google Scholar 

  4. Tran D, Roberts C, Jorm L, Seeho S, Havard A. Change in smoking status during two consecutive pregnancies: a population-based cohort study. BJOG. 2014;121(13):1611–20.

    Article  CAS  Google Scholar 

  5. Koren G. Nicotine-replacement therapy in pregnancy-the end of the road? Ther Drug Monit. 2012;34(5):487–8.

    Article  CAS  PubMed  Google Scholar 

  6. Haddad A, Davis AM. Tobacco smoking cessation in adults and pregnant women: behavioral and pharmacotherapy interventions. JAMA. 2016;315(18):2011–2.

    Article  CAS  PubMed  Google Scholar 

  7. Cahill K, Stevens S, Perera R, Lancaster T. Pharmacological interventions for smoking cessation: an overview and network meta-analysis. Cochrane Database Syst Rev. 2013;5:CD009329.

    Google Scholar 

  8. Turner E, Jones M, Vaz LR, Coleman T. Systematic review and meta-analysis to assess the safety of bupropion and varenicline in pregnancy. Nicotine Tob Res. 2019;21(8):1001–10.

    Article  PubMed  Google Scholar 

  9. The American College of Obstetricians and Gynecologists. Committee opinion no. 721: smoking cessation during pregnancy. Obstet Gynecol. 2017;130(\):e200.

    Google Scholar 

  10. World Health Organization. WHO recommendations for the prevention and management of tobacco use and second-hand smoke exposure in pregnancy. Geneva: World Health Organization; 2013.

    Google Scholar 

  11. Kaplan YC, Olgac Dündar N, Kasap B, Karadas B. Pregnancy outcome after varenicline exposure in the first trimester. Case Rep Obstet Gynecol. 2014;2014:263981.

    PubMed  PubMed Central  Google Scholar 

  12. Harrison-Woolrych M, Paterson H, Tan M. Exposure to the smoking cessation medicine varenicline during pregnancy: a prospective nationwide cohort study. Pharmacoepidemiol Drug Saf. 2013;22(10):1086–92.

    CAS  PubMed  Google Scholar 

  13. Richardson JL, Stephens S, Yates LM, Diav-Citrin O, Arnon J, Beghin D, et al. Pregnancy outcomes after maternal varenicline use; analysis of surveillance data collected by the European Network of Teratology Information Services. Reprod Toxicol. 2017;67:26–34.

    Article  CAS  PubMed  Google Scholar 

  14. Stotts AL, Northrup TF, Cinciripini PM, Minnix JA, Blalock JA, Mullen PD, et al. Randomized, controlled pilot trial of bupropion for pregnant smokers: challenges and future directions. Am J Perinatol. 2015;32(4):351–6.

    PubMed  Google Scholar 

  15. Nanovskaya TN, Oncken C, Fokina VM, Feinn RS, Clark SM, West H, et al. Bupropion sustained release for pregnant smokers: a randomized, placebo-controlled trial. Am J Obstet Gynecol. 2017;216(4):420.e1–9.

    Article  CAS  Google Scholar 

  16. Chun-Fai-Chan B, Koren G, Fayez I, Kalra S, Voyer-Lavigne S, Boshier A, et al. Pregnancy outcome of women exposed to bupropion during pregnancy: a prospective comparative study. Am J Obstet Gynecol. 2005;192(3):932–6.

    Article  CAS  PubMed  Google Scholar 

  17. Bérard A, Zhao J-P, Sheehy O. Success of smoking cessation interventions during pregnancy. Am J Obstet Gynecol. 2016;215(5):611.e1–8.

    Article  Google Scholar 

  18. Berlin I, Grangé G, Jacob N, Tanguy M-L. Nicotine patches in pregnant smokers: randomised, placebo controlled, multicentre trial of efficacy. BMJ. 2014;348:g1622.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  19. Lassen TH, Madsen M, Skovgaard LT, Strandberg-Larsen K, Olsen J, Andersen AM. Maternal use of nicotine replacement therapy during pregnancy and offspring birthweight: a study within the Danish National Birth Cohort. Paediatr Perinat Epidemiol. 2010;24(3):272–81.

    Article  PubMed  Google Scholar 

  20. Strandberg-Larsen K, Tinggaard M, Nybo Andersen AM, Olsen J, Gronbaek M. Use of nicotine replacement therapy during pregnancy and stillbirth: a cohort study. BJOG. 2008;115(11):1405–10.

    Article  CAS  PubMed  Google Scholar 

  21. Gaither KH, Brunner Huber LR, Thompson ME, Huet-Hudson YM. Does the use of nicotine replacement therapy during pregnancy affect pregnancy outcomes? Matern Child Health J. 2009;13(4):497–504.

    Article  PubMed  Google Scholar 

  22. Havard A, Jorm LR, Preen D, Daube M, Kemp A, Einarsdóttir K, et al. The Smoking MUMS (Maternal Use of Medications and Safety) Study: protocol for a population-based cohort study using linked administrative data. BMJ Open. 2013;3(9).

    Article  PubMed  PubMed Central  Google Scholar 

  23. Roper L, Tran DT, Einarsdóttir K, Preen DB, Havard A. Algorithm for resolving discrepancies between claims for smoking cessation pharmacotherapies during pregnancy and smoking status in delivery records: the impact on estimates of utilisation. PLoS One. 2018;13(8):e0202999.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  24. Tran DT, Havard A, Jorm LR. Data cleaning and management protocols for linked perinatal research data: a good practice example from the Smoking MUMS (Maternal Use of Medications and Safety) Study. BMC Med Res Methodol. 2017;17(1):97.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Tran DT, Roberts CL, Havard A, Jorm LR. Linking birth records to hospital admission records enhances the identification of women who smoke during pregnancy. Aust N Z J Public Health. 2014;38(3):258–64.

    Article  PubMed  Google Scholar 

  26. National Centre for Classification in Health. The International Statistical Classification of Diseases and Related Health Problems, Tenth Revision, Australian Modification (ICD-10-AM)-Fifth Edition-Tabular list of diseases and Alphabetic index of diseases. Sydney: National Centre for Classification in Health, Faculty of Health Sciences, The University of Sydney; 2008.

    Google Scholar 

  27. New South Wales Department of Health. Maternity - Towards Normal Birth in NSW. Maternity and Child Health Unit - Strategic Development Division 2010.

    Google Scholar 

  28. Hutchinson M. Western Australia’s mothers and babies, 2012: 30th annual report of the Western Australian midwives notification system; 2015.

    Google Scholar 

  29. WHO Collaborating Centre for Drug Statistics Methodology. ATC alterations from 2005-2019. Geneva: World Health Organization; 2019. Accessed 7 July 2019

    Google Scholar 

  30. Dobbins TA, Sullivan EA, Roberts CL, Simpson JM. Australian national birthweight percentiles by sex and gestational age, 1998-2007. MJA. 2012;197(5):291–4.

    PubMed  Google Scholar 

  31. Lain SJ, Algert CS, Nassar N, Bowen JR, Roberts CL. Incidence of severe adverse neonatal outcomes: use of a composite indicator in a population cohort. Matern Child Health J. 2012;16(3):600–8.

    Article  PubMed  Google Scholar 

  32. Roberts CL, Cameron CA, Bell JC, Algert CS, Morris JM. Measuring maternal morbidity in routinely collected health data: development and validation of a maternal morbidity outcome indicator. Med Care. 2008;46(8):786–94.

    Article  PubMed  Google Scholar 

  33. Conde-Agudelo A, Rosas-Bermudez A, Kafury-Goeta AC. Birth spacing and risk of adverse perinatal outcomes: a meta-analysis. JAMA. 2006;295(15):1809–23.

    Article  CAS  PubMed  Google Scholar 

  34. Appareddy S, Pryor J, Bailey B. Inter-pregnancy interval and adverse outcomes: evidence for an additional risk in health disparate populations. J Matern Fetal Neonatal Med. 2017;30(21):2640–4.

    Article  PubMed  Google Scholar 

  35. Australian Therapeutic Goods Administration. Reporting adverse events: Australian Department of Health; 2018. Accessed 25 Oct 2019.

    Google Scholar 

  36. Australian Bureau of Statistics 2033.0.55.001. Census of population and housing: socio-economic indexes for areas (SEIFA), Australia 2006: Australian Bureau of Statistics; 2008. Accessed 5 Feb 2018.

  37. Australian Bureau of Statistics. Australian Standard Geographic Classification Remoteness Structure: Commonwealth of Australia; 2010. Accessed 15 Apr 2019.

  38. Parsons LS. Performing a 1:N case-control match on propensity sore. Proceedings of the Twenty -Sixth Annual SAS Users Group International Conference; 2004.

    Google Scholar 

  39. Nguyen T-L, Collins GS, Spence J, Daurès J-P, Devereaux PJ, Landais P, et al. Double-adjustment in propensity score matching analysis: choosing a threshold for considering residual imbalance. BMC Med Res Methodol. 2017;17(1):78.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Hartmann-Boyce J, Stead LF, Cahill K, Lancaster T. Efficacy of interventions to combat tobacco addiction: Cochrane update of 2013 reviews. Addiction. 2014;109(9):1414–25.

    Article  PubMed  Google Scholar 

  41. McCowan LME, Dekker GA, Chan E, Stewart A, Chappell LC, Hunter M, et al. Spontaneous preterm birth and small for gestational age infants in women who stop smoking early in pregnancy: prospective cohort study. BMJ. 2009;338:b1081.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Soneji S, Beltrán-Sánchez H. Association of maternal cigarette smoking and smoking cessation with preterm birth. JAMA Netw Open. 2019;2(4):e192514–e.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Blatt K, Moore E, Chen A, Van Hook J, DeFranco EA. Association of reported trimester-specific smoking cessation with fetal growth restriction. Obstet Gynecol. 2015;125(6):1452–9.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Centre for Epidemiology and Evidence. Quality and coverage of the NSW Register of Congenital Conditions using Admitted Patient Data: a record linkage study. Sydney: NSW Ministry of Health; 2016. Accessed 4 Apr 2018

    Google Scholar 

  45. Centers for Disease Control and Prevention. Update on overall prevalence of major birth defects-Atlanta, Georgia, 1978-2005. MMWR Morb Mortal Wkly Rep. 2008;57(1):1–5.

    Google Scholar 

  46. European Surveillance of Congenital Anomalies. EUROCAT prevalence data tables. Cases and prevalence (per 10,000 births) of all congenital anomaly subgroups for all registries, from 2012–2016 2018. Accessed 9 Apr 2019.

    Google Scholar 

  47. Pfizer Laboratories. Chantix (full prescribing information) 2019. Accessed 2 Apr 2019.

    Google Scholar 

  48. Scialli AR, Daston G, Chen C, Coder PS, Euling SY, Foreman J, et al. Rethinking developmental toxicity testing: evolution or revolution? Birth Defects Res. 2018;110(10):840–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Bracken MB. Why animal studies are often poor predictors of human reactions to exposure. J R Soc Med. 2009;102(3):120–2.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Fokina VM, West H, Oncken C, Clark S, Ahmed MS, Hankins G, et al. Bupropion therapy during pregnancy: concentrations of the drug and its major metabolites in umbilical cord plasma and amniotic fluid. Drug Alcohol Depend. 2017;171:e63.

    Article  Google Scholar 

  51. Ross EJ, Graham DL, Money KM, Stanwood GD. Developmental consequences of fetal exposure to drugs: what we know and what we still must learn. Neuropsychopharmacology. 2015;40(1):61–87.

    Article  CAS  PubMed  Google Scholar 

  52. Dempsey D, Jacob P, Benowitz NL. Accelerated metabolism of nicotine and cotinine in pregnant smokers. J Pharmacol Exp Ther. 2002;301(2):594.

    Article  CAS  PubMed  Google Scholar 

  53. Jackson JW, Schmid I, Stuart EA. Propensity scores in pharmacoepidemiology: beyond the horizon. Curr Epidemiol Rep. 2017;4(4):271–80.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Suissa S. Immortal time bias in pharmacoepidemiology. Am J Epidemiol. 2007;167(4):492–9.

    Article  PubMed  Google Scholar 

  55. Taylor L, Pym M, Bajuk B, Sutton L, Travis S, Banks C. Validation study NSW midwives data collection 1998. Sydney: NSW Health Department; 2000.

    Google Scholar 

  56. Downey F. Validation study of the Western Australian midwives notification system. 2005 data. Department of Health, Western Australia: Perth; 2007.

  57. US Food and Drug Administration. FDA drug safety communication: FDA revises description of mental health side effects of the stop-smoking medicines Chantix (varenicline) and Zyban (bupropion) to reflect clinical trial findings 2016. Accessed 25 Oct 2019.

    Google Scholar 

  58. Davey ME, Hunt JM, Foster R, Couzos S, van der Sterren AE, Sarin J, et al. Tobacco control policies and activities in aboriginal community-controlled health services. Med J Aust. 2015;202(10):S63–6.

    PubMed  Google Scholar 

  59. Greenhalgh EM, Scollo MM, Pearce M. 9.1 socio-economic position and disparities in tobacco exposure and use. In: Scollo MM, Winstanley MH, editors. Tobacco in Australia: facts and issues. Melbourne: Cancer Council Victoria; 2018.

    Google Scholar 

  60. Cancer Council New South Wales. NRT factsheet 2016. Accessed 25 Oct 2019.

Download references


The authors would like to thank the NSW Ministry of Health, the WA Department of Health, the Australian Government Department of Health, the Department of Human Services, the NSW Centre for Health Record Linkage, the Australian Institute for Health and Welfare, the Western Australia Data Linkage Branch and data custodians of the NSW Perinatal Data Collection, WA Midwife Notification System, NSW Admitted Patient Data Collection, WA Hospital Morbidity Data Collection, NSW Registry of Births, Deaths and Marriages, WA Registry of Births, Deaths and Marriages, and Pharmaceutical Benefits Scheme data for allowing access to the data and conducting the linkage of records. The authors would also like to thank the Centre for Big Data Research in Health Aboriginal and Torres Strait Islander Maternal and Child Reference Group for their input.


This study was funded by the Australian National Health and Medical Research Council (NHMRC #1028543). Author AH was supported by a National Heart Foundation of Australia Future Leader Fellowship (#100411). The funding source had no influence on the design and conduct of the study; collection, management, analysis and interpretation of the data; preparation, review or approval of the manuscript; or the decision to submit the manuscript for publication.

Author information

Authors and Affiliations



AH conceived the study, and all authors designed the study. AH acquired the data. DTT performed the analyses and wrote the manuscript. All authors participated in the interpretation of the results and critical review of the manuscript. All authors read and approved the final manuscript. The corresponding author attests that all listed authors meet authorship criteria and that no others meeting the criteria have been omitted.

Corresponding author

Correspondence to Duong Thuy Tran.

Ethics declarations

Ethics approval and consent to participate

The project was approved by the Australian Institute of Health and Welfare Ethics Committee, the NSW Population and Health Services Research Ethics Committee, the Department of Health WA Human Research Ethics Committee, the Aboriginal Health and Medical Research Council of NSW Ethics Committee and the Western Australian Aboriginal Health Ethics Committee. The study used routinely collected data that have been anonymised. Waiver of consent to participate was obtained.

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1:

Description of study variables, missing data and exclusion.

Additional file 2:

Maternal characteristics before and after matching, results of sensitivity analyses.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Tran, D.T., Preen, D.B., Einarsdottir, K. et al. Use of smoking cessation pharmacotherapies during pregnancy is not associated with increased risk of adverse pregnancy outcomes: a population-based cohort study. BMC Med 18, 15 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: