Long-term mortality in mothers of infants with neonatal abstinence syndrome: A population-based parallel-cohort study in England and Ontario, Canada

Summary

Background

Opioid addiction is a major public health threat to healthy life expectancy; however, little is known of long-term mortality for mothers with opioid use in pregnancy. Pregnancy and delivery care are opportunities to improve access to addiction and supportive services. Treating neonatal abstinence syndrome (NAS) as a marker of opioid use during pregnancy, this study reports long-term maternal mortality among mothers with a birth affected by NAS in relation to that of mothers without a NAS-affected birth in 2 high-prevalence jurisdictions, England and Ontario, Canada.

Methods and findings

We conducted a population-based study using linked administrative health data to develop parallel cohorts of mother–infant dyads in England and Ontario between 2002 and 2012. The study population comprised 13,577 and 4,966 mothers of infants with NAS and 4,205,675 and 929,985 control mothers in England and Ontario, respectively. Death records captured all-cause maternal mortality after delivery through March 31, 2016, and cause-specific maternal mortality to December 31, 2014. The primary exposure was a live birth of an infant with NAS, and the main outcome was all deaths among mothers following their date of delivery. We modelled the association between NAS and all-cause maternal mortality using Cox regression, and the cumulative incidence of cause-specific mortality within a competing risks framework. All-cause mortality rates, 10-year cumulative incidence risk of death, and crude and age-adjusted hazard ratios were calculated. Estimated crude 10-year mortality based on Kaplan–Meier curves in mothers of infants with NAS was 5.1% (95% CI 4.7%–5.6%) in England and 4.6% (95% CI 3.8%–5.5%) in Ontario versus 0.4% (95% CI 0.41%–0.42%) in England and 0.4% (95% CI 0.38%–0.41%) in Ontario for controls (p < 0.001 for all comparisons). Survival curves showed no clear inflection point or period of heightened risk. The crude hazard ratio for all-cause mortality was 12.1 (95% 11.1–13.2; p < 0.001) in England and 11.4 (9.7–13.4; p < 0.001) in Ontario; age adjustment did not reduce the hazard ratios. The cumulative incidence of death was higher among NAS mothers than controls for almost all causes of death. The majority of deaths were by avoidable causes, defined as those that are preventable, amenable to care, or both. Limitations included lack of direct measures of maternal opioid use, other substance misuse, and treatments or supports received.

Conclusions

In this study, we found that approximately 1 in 20 mothers of infants with NAS died within 10 years of delivery in both England and Canada—a mortality risk 11–12 times higher than for control mothers. Risk of death was not limited to the early postpartum period targeted by most public health programs. Policy responses to the current opioid epidemic require effective strategies for long-term support to improve the health and welfare of opioid-using mothers and their children.

Keywords:

Death rates – Neonates – Labor and delivery – pregnancy – opioids – Infants – Ontario – England

Introduction

Opioid use is responsible for an important increase in premature mortality in young and middle-aged adults in the US [1] and Canada [2], 2 of the countries with the highest per capita prescription opioid consumption in Western industrialized nations [3]. Other countries such as England have seen similar rates of increase in prescription opioid use but not concomitant increases in mortality rates, likely related in part to better access to addiction treatment and more oversight of prescription opioids [4]. Across all of these jurisdictions, there is increasing opioid use by pregnant women, and while little is known about associated maternal mortality, a recent study using data from 22 US states and the District of Columbia reports a higher than 3-fold increase from 2007 to 2017 in opioid-related deaths in women during or within the first year after pregnancy [5].

Population-based surveillance of opioid use during pregnancy is difficult given the lack of prescription medication data in many jurisdictions and the challenges in measuring illicit use. Neonatal abstinence syndrome (NAS) is coded in the infant birth hospitalization record and offers a widely used but imperfect proxy measure of maternal opioid use during pregnancy. NAS manifests typically within hours to 1 to 2 days of delivery with autonomic, gastrointestinal, and neurologic symptoms of drug withdrawal in the infant that often require prolonged postnatal care [6]. Not all infants exposed to opioids in utero will experience NAS [7,8]. In trial and observational settings, approximately half of women receiving methadone or buprenorphine maintenance therapy in pregnancy gave birth to an infant with signs of NAS [7], and estimates of up to 91% have been reported in other groups of women with chronic opioid use [6]. The incidence of NAS rose dramatically between the early 2000s and 2014 in the US (from 28 to 144 per 10,000 births) [9] and Canada (from 18 to 54 per 10,000 live births) [10], but remained relatively stable in Australia and England [11,12].

Mothers of infants with NAS represent a heterogeneous group including those using prescription opioids and opioid agonists for analgesia, those on medically supervised maintenance treatment for dependence, and those using illicit opioids [12]. Overall social disadvantage and opioid use are inextricably linked, and many women who misuse these drugs are at greater risk of adversity, including deprivation, violence and abuse, and use of other substances [13]. All of these factors impact negatively on maternal health and may diminish parenting capacity [14]. Improving maternal health and preventing premature maternal mortality in this population is therefore critical for both mothers and their children. Pregnancy and delivery care are opportunities to improve access to addiction and supportive services. While public health programs such as nurse home visits tend to focus on pregnancy and the early postpartum period [15], it is largely unknown whether this is the only period of risk for poor outcomes for mothers who use opioids in pregnancy.

Evidence on mortality for mothers who use opioids in pregnancy is limited but consistently shows increased rates around the time of delivery (Table 1) [5,16–18]. Very high rates of longer-term maternal mortality have been reported in 2 older studies (1 Australian and 1 Finnish) [19,20], but these studies may not reflect the current opioid epidemic in North America or the UK [4]. Recent estimates of mortality from a meta-analysis of people with substance misuse disorder and homeless and prison populations reported a standardized mortality ratio for women of 11.9 (95% CI 10.4–13.3), which was higher than the equivalent figure for men (7.9; 95% CI 7.0–8.7) [21]. There is a dearth of information about long-term health outcomes for women—particularly mothers—with opioid use, which is an important knowledge gap given rising rates of prescription and illicit opioid use. Pregnancy can be seen as a window of opportunity for identifying and managing substance misuse and the implications for parenting capacity.

In this study, we capitalize on linked population-based maternal–infant healthcare records and mortality files in 2 jurisdictions (England and Ontario, Canada). Both have similar healthcare systems, including universal access to healthcare and similar postnatal public health programs that focus predominantly on the year after birth. We hypothesized that NAS mothers would have significantly higher rates of mortality than control mothers. We aimed to quantify this excess mortality, investigating all-cause mortality as the primary outcome and cause-specific mortality as secondary outcomes. We present all-cause mortality in relation to some key maternal characteristics at the time of birth, to identify clinically useful sub-groups with poor prognostic outcomes and to guide opioid misuse policy and research.

Methods

An analytic plan was written and approved before starting statistical analyses (S1 Text). All analyses were undertaken as planned and reported, with the exception of the removal of comorbidity adjustment in the Cox proportional hazards model (see section on statistical analysis below) as a result of feedback from peer review. This study is reported as per the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guidelines (S2 Text).

Cohort

We derived whole-region population-based cohorts of mothers aged 12 to 49 years and their live born infants delivered between April 1, 2002, and December 31, 2012, in England and Ontario, Canada. This period coincides with a steep increase in opioid use, particularly in Canada. The cohorts used longitudinal hospital discharge records for mothers (back to April 1, 1997, for some covariates) linked to hospitalization records for the infant. Mortality registration records were linked to hospital records for the mother and infant. In England, de-identified data on inpatient admissions for all National Health Service (NHS) hospitals linked to Office for National Statistics mortality data were obtained from NHS Digital, linked between mother and infant using previously reported methods [30], and analyzed within the UCL Data Safe Haven, England [31]. In Ontario, hospital discharge records were obtained from the Canadian Institute for Health Information Discharge Abstract Database and the Ontario Mental Health Reporting System. Cause of death information was taken from the Office of the Registrar General’s vital statistics database (data available only until 2014), and demographic information captured from the 2006 Canadian Census and the Registered Persons Database (which includes date of death). Canadian datasets were linked using unique coded identifiers common across all aforementioned datasets and analyzed at ICES in Toronto, Ontario, Canada. We restricted the cohort to singleton births, and if a woman had more than 1 live birth delivery during the study period, 1 delivery was chosen at random as the focus of the study, i.e., a delivery date was selected at random and used as the entry point for the mother (referred to as the index delivery), and all other deliveries were ignored.

In both jurisdictions, infants with a diagnosis of NAS were identified using International Classification of Diseases and Related Health Problems—10^th Revision (ICD-10) codes P96.1 (neonatal withdrawal symptoms from maternal use of drugs of addiction) or P04.4 (newborn [suspected to be] affected by maternal use of drugs of addiction) recorded during the delivery admission or subsequent readmission within 14 days of birth [11]. The ICD-9 equivalent to P96.1 has been shown to have high sensitivity (88.1%; 95% CI 83.3%–91.7%) and specificity (97.0%; 95% CI 93.8%–98.5%) and a positive predictive value of 91.2% (95% CI 86.8%–94.2%) for measuring NAS [32] and is used by the US Agency for Healthcare Research and Quality [33]. We included P04.4 as it is often used for opioid withdrawal in both jurisdictions. In a sample of all Ontario women giving birth in hospital from 2014 to 2017 (n = 464,400), during which time all prescription opioids were registered, over 56% of women whose infant had a diagnosis of P04.4 had either a prescription for opioids or opioid agonists, or a healthcare encounter related to opioid use during pregnancy (personal communication, A. Camden, University of Toronto, September 4, 2019). In our study, NAS was ascertained by P96.1 in 83% of cases in England and 65.7% in Ontario. The English maternal cohort was identified using linked data for singleton babies and mothers, corresponding to 96% of all live births in England within the study period. Mothers and babies were matched deterministically using data on hospital, general practitioner practice, maternal age, birthweight, gestation, birth order, and sex, or probabilistically using additional data including admission dates, ethnicity, and partial postal code [34]. Ontario mothers were identified using a unique number linking all newborn and maternal hospital records that is assigned at the delivery hospitalization.

Outcomes and covariates

The main outcome was all-cause mortality measured from April 1, 2002, to March 31, 2016, derived from linked death registrations. Cause-specific mortality was available only up until December 31, 2014, for Ontario. Cause-specific deaths were classified as avoidable (defined as preventable, amenable to care, or both), unavoidable (as specified by the UK’s Office for National Statistics) [35], or cancer (avoidable and unavoidable).

Longitudinal hospital records were used to derive baseline characteristics at the index delivery including maternal age (12–19 years, 20–34 years, and 35+ years); time since last birth (defined as no previous births [with lookback to April 1, 1997, only], <2 years, 2–5 years, and 6+ years); neighbourhood income quintile; urban or rural residence; pregnancy-related outcomes: cesarean delivery, gestational age at delivery (<34 weeks, 34–36 weeks, and 37+ weeks), gestational hypertension, pre-eclampsia/eclampsia, and gestational diabetes; neonatal mortality (death within 28 days of birth); and infant discharge from hospital to social services. We excluded observations with extreme values, including women aged <12 years or >49 years at the date of delivery. Missing values for neighbourhood income quintile (n = 3,850 [0.41%] for Ontario; n = 24,067 [0.57%] for England) were categorized into the lowest income quintile, and missing area of residence information (n = 151 [0.02%] for Ontario; n = 22,970 [0.54%] for England) was categorized into the urban area of residence. In Ontario, these missing data are suppressed for neighbourhoods with high rates of residential instability, which are predominantly low income and urban.

Baseline maternal morbidities were measured for a 5-year lookback period using all diagnostic codes in all hospital admissions within this period to assign the Charlson comorbidity index (Deyo version) (0, 1, or 2+ comorbid conditions) and to identify hospitalizations related to (i) any psychiatric condition, (ii) addiction-related conditions, and (iii) other mental health conditions. We chose the Deyo version of the Charlson comorbidity index as it has been most widely used in maternal mortality studies and, unlike other indices, has been validated on longer-term mortality, although not in pregnant women [36,37]. S1 and S2 Tables list diagnostic codes and provide definitions for neighbourhood income quintile and urban or rural residence.

Statistical analysis

Descriptive statistics

We compared mothers with an infant affected by NAS and controls within each jurisdiction. We compared baseline characteristics using Pearson’s chi-squared test (categorical variables) and ANOVA (continuous variables).

All-cause mortality

Our primary outcome was all-cause mortality of mothers after the birth of an infant with NAS, relative to control mothers. Survival analysis for time to all-cause mortality was modelled using multivariable Cox regression, with proportionality of hazards assessed by Schoenfeld residuals and log–log plots, and included all deaths to the end of the study period. Crude and adjusted hazard ratios were produced that were adjusted only for maternal age group at delivery to describe the extent of the mortality gap between NAS mothers and controls. We explicitly did not attempt to adjust models for other covariates as we do not attempt to infer causality. This is primarily because our data reflect maternal opioid use at a single time point (delivery). Thus we are unable to examine the direction of effect for key factors such as mental illness or socioeconomic status that may either confound the association between opioid use and death or lie on the causal pathway. We initially adjusted the model for the Charlson comorbidity index but as a result of peer review do not report these estimates because of our aim to not make causal inferences. Crude and adjusted survival curves were plotted to estimate the absolute risk of mortality at 5 and 10 years after birth. We derived 95% confidence intervals for mortality at 5 and 10 years after delivery through log–log transformation of the survival function and computed p-values using the z test [38]. Mothers with missing values for maternal age at delivery were excluded from the adjusted models. Individuals surviving beyond the end of the study period were censored at March 31, 2016.

Age-standardized all-cause mortality rates were estimated for England and Ontario using the direct method of standardization and the Canadian 2006 Census as the standard population. We present age-standardized all-cause mortality rates stratified by maternal age at delivery (excluding those with missing age information), neighbourhood income quintile (Q1 [most deprived] versus Q2–Q5 due to low numbers of deaths in the most affluent quintiles), urban/rural residence, previous addiction-related or other mental-health-related hospitalization, Charlson comorbidity index (0 versus 1+ due to low numbers of deaths in the 2+ category for NAS mothers), and infant discharge to social service out-of-home care.

Cause-specific mortality

Ten-year cumulative incidence of death was calculated for each cause of death category with consideration of other causes of death as competing risk events using Gray’s test for the homogeneity of 2 or more cumulative incidence functions. Individuals surviving beyond the end of the study period were censored at December 31, 2014, for this particular analysis.

Statistical analyses were performed using Stata version 15 for England data and SAS version 9.4 statistical software for Ontario data in a Unix environment. p-Values for age-standardized all-cause mortality and cause-specific mortality were derived using the z test.

ICES is a prescribed entity under section 45 of Ontario’s Personal Health Information Protection Act. Section 45 authorizes ICES to collect personal health information, without consent, for the purpose of analysis or compiling statistical information with respect to the management of, evaluation or monitoring of, the allocation of resources to, or planning for all or part of the health system. Projects conducted under section 45, by definition, do not require review by a research ethics board. This project was conducted under section 45, and approved by the ICES Privacy and Legal Office. Research ethics approval was granted by the Hospital for Sick Children Research Ethics Board for Ontario analyses. The English analyses were exempt from UK NHS Research Ethics Committee approval because it involved the analysis of de-identified administrative data.

Results

Cohort characteristics of mothers in England and Ontario

After applying our inclusion criteria, there were 13,577 and 4,966 mothers of infants with NAS and 4,205,675 and 929,985 control mothers in England and Ontario, respectively. Baseline characteristics are described in Table 2. In both jurisdictions, the majority of mothers had no previous recorded hospital birth (lookback to April 1, 1997) and lived in urban areas. Compared to controls, a larger proportion of mothers of infants with NAS lived in neighbourhoods in the lowest income quintile (37.3% versus 23.5% in England and 44.9% versus 22.9% in Ontario; p < 0.001 for both jurisdictions); mothers of infants with NAS were on average younger than controls, and, in Ontario, a higher proportion were teenage mothers (9.8% versus 3.7%; p < 0.001). Mothers of infants with NAS were more likely to have higher comorbidity scores (1 or 2+ on the Charlson comorbidity index) than controls (21.2% versus 6.9% in England and 7.6% versus 1.7% in Ontario; p < 0.001 for both jurisdictions), and a greater proportion also had a previous psychiatric hospitalization (13.5% versus 0.6% in England and 15.3% versus 1.1% in Ontario; p < 0.001 for both jurisdictions). Infant discharge to care by social services was much more common among infants with NAS than among control infants in both England (9.7% versus 0.1%; p < 0.001) and Ontario (15.2% versus 0.1%; p < 0.001). We found no differences in neonatal mortality between infants with NAS and controls in Ontario, but a marginally higher rate among infants with NAS in England (0.3% versus 0.2% p = 0.01).

Risk of death among mothers of infants born with NAS

From 2002 to 2016, there were 112,890 total person-years of follow-up for mothers of infants with NAS in England (34.9 million for controls) and 35,740 total person-years of follow-up for mothers of infants with NAS in Ontario (7.9 million for controls). The mean duration of follow-up was 8.4 years (England) and 7.2 years (Ontario) for NAS mothers and 8.2 years (England) and 8.6 years (Ontario) for controls. Crude mortality rate for mothers of infants with NAS was 5.01 (95% CI 4.62–5.44) per 1,000 person-years in England and 4.28 (95% CI 3.63–5.02) per 1,000 person-years in Ontario. In both jurisdictions, the crude cumulative mortality incidence (superimposed on the survival curves in Fig 1) was significantly different between mothers of infants with NAS and controls. In England, 5- and 10-year mortality (95% CI) was 1.81% (1.59%–2.05%) and 5.13% (4.69–5.62%), respectively, for mothers of infants with NAS and 0.15% (0.15%–0.16%) and 0.42% (0.41–0.43%), respectively, for controls (p < 0.001 for 5- and 10-year mortality for NAS mothers versus controls); in Ontario, it was 1.85% (1.51%–2.28%) and 4.58% (3.81%–5.49%), respectively, for mothers of infants with NAS and 0.15% (0.15%–0.16%) and 0.40% (0.38%–0.41%), respectively, for controls (p < 0.001 for 5- and 10-year mortality for NAS mothers versus controls). The decline in survival of mothers of infants with NAS over time was steady in both jurisdictions, with no clear inflection point or distinct period of risk.

Table 3 outlines age-standardized mortality rates overall and stratified by clinical and demographic characteristics for the whole study period. In both jurisdictions, rates of death were higher among mothers of infants with NAS compared to controls across most risk groups. Differences in rates were not significant (95% confidence intervals overlapped) for NAS mothers and controls for those with a history of addiction-related admission (England and Ontario) or who had an infant discharged to social services (Ontario only) or who were 19 years old or under at delivery (Ontario only).

The crude hazard ratio for all-cause mortality among mothers of infants with NAS was 12.1 (95% CI 11.1–13.2; p < 0.001) in England and 11.4 (95% CI 9.7–13.4; p < 0.001) in Ontario (Table 4). After adjustment for age, the hazard ratio increased for England and was unchanged for Ontario (England adjusted HR 13.0; 95% CI 11.9–14.1; p < 0.001; and Ontario adjusted HR 11.4; 95% CI 9.7–13.4; p < 0.001).

Table 5 presents results for cause-specific mortality between 2002 and 2014 (cause-specific data for 2015–2016 were not available in Ontario). Avoidable deaths were the most common cause of death among mothers of infants with NAS in both jurisdictions (accounting for >85% in England and 75% in Ontario), with a 10-year cumulative incidence risk of 42.9 deaths per 1,000 population (95% CI 38.4–47.9) among English mothers and 30.8 deaths per 1,000 population (95% CI 24.1–38.8) among Ontario mothers. Intentional and unintentional injuries (e.g., transport injuries, unintentional falls) made up the majority of avoidable mortality in the mothers of infants with NAS in both jurisdictions.

Discussion

In this large population-based study across 2 countries, 1 in 20 mothers of infants with NAS died within 10 years of delivery—a mortality risk that was 11–12 times higher than for control mothers. Findings were consistent across both jurisdictions. For virtually all causes of death, mortality rates were substantially higher for mothers of infants with NAS than for controls, with the majority of deaths attributable to avoidable causes such as intentional and unintentional injuries. We also identified universally high mortality rates among mothers who had a history of hospitalization for addiction, irrespective of whether or not their infant had NAS. We found no evidence of a high risk period in the 1–2 years after birth (corresponding to the period typically targeted by public health nursing or other support for high-risk families) for maternal deaths in the NAS group.

Other population-based studies report increased perinatal maternal mortality in mothers using opioids [5,16,17] and longer-term risk in mothers with alcohol or drug misuse during pregnancy [20,22] (Table 1). None of these studies address long-term all-cause mortality for mothers of NAS-affected infants as we have done. High mortality rates have been reported for non-pregnant populations of female opioid users. Notably, 3 large-scale studies present crude mortality rates per 1,000 person-years of 6.5 (95% CI 6.1–6.9) in New South Wales [39], between 7.5 and 13.9 in opioid-using women aged 15–44 years in California [29], and 12.2 (95% CI 10.3–14.4) and 19.7 (95% CI 15–25.8), respectively, for female users of heroin and other opioids in Denmark [40]. Our crude mortality rates per 1,000 person-years for mothers with infants with NAS of 5.01 (95% CI 4.62–5.44) in England and 4.28 (95% CI 3.63–5.02) in Ontario are comparatively lower, which is likely driven in part by a relatively younger age distribution in our study population. For example, in the New South Wales study, Degenhardt et al. included all ages and reported increasing mortality rates with older age. The higher risk of all-cause mortality we report is not surprising and has been shown in other marginalized groups or those with mental health problems [21,41]. In particular, unintentional injury deaths (which include those related to victimization) predominate and may result from social vulnerability or misclassification of injuries that are intentional [42–46].

In our study, the high rate of premature mortality among the mothers of infants with NAS was mirrored by high rates for the mothers of control infants discharged to social services and for mothers with a history of hospitalization for addiction. Other studies have demonstrated maternal well-being declining in association with the loss of a child to foster care [42,47], whereas retaining care of the child may help facilitate treatment [48]. An estimated 7%–20% of NAS-affected infants do not return home with their mother at the time of postnatal discharge from hospital [49–52], which is similar to the percentage in our study (10%–15%). However, these figures may not reflect the much higher cumulative risk of foster care placements occurring later in childhood [53]. Our findings also mirror the growing body of literature describing the constellation of psychosocial risk factors linking mental illness, addiction, and social adversity [43,54] and suggest the need for multifaceted support for these mothers irrespective of whether their children are living with them.

The longitudinal and population-based nature of this study, its size, and comparison of similar universal healthcare systems are strengths. Limitations include potential linkage error, misclassification of mothers using opioids whose babies did not develop NAS, and lack of direct measures of maternal opioid use or treatment, or other substance misuse, which may underestimate the burden of mortality in mothers with opioid use within our study or make findings less generalizable to opioid-using mothers whose infants do not have NAS or who live in other jurisdictions with other approaches to treatment and available supportive services. In the extract of English hospital data, 96% of live births were matched to maternal records, but linkage was lower (88%) for infants with NAS. Comparison of the characteristics of these English infants with NAS by linkage status indicated an association between non-linkage and both longer hospital stays and greater risk of placement in out-of-home care [55], with the implication that our results may underestimate the risk of death among women with opioid use. Our results reflect the subset of mothers with prenatal opioid exposure who gave birth to an infant with NAS, and as such they likely underestimate the true risk of maternal mortality associated with prenatal opioid use. However, we also describe mortality in mothers without an infant with NAS, but who received care for mental illness and addictions, thus broadening the scope for generalizability. Our study cohort may also include the rare cases of NAS related to withdrawal from other substances [56] or from postnatal opioid use (“iatrogenic NAS”) [57], but data on these other exposures were not available.

Our study has implications for research, practice, and policy to improve maternal and, arguably, child outcomes related to prenatal opioid use. Enhanced treatment programs for opioid dependence that integrate maintenance therapy, psychotherapy, reproductive health, and obstetric care have been found to be effective in reducing substance misuse, unplanned pregnancies, and obstetric complications during the perinatal period [58]. Some evaluations of programs supporting mothers with opioid use and their children suggest that multifaceted services addressing health, addiction, housing, and parenting needs can improve parenting capacity and attachment and reduce child apprehension [56,59–64]. However, rigorous evidence on interventions promoting long-term support is limited and should be a research and policy priority. New funding for child welfare agencies in the US to provide services related to mental health, addictions, and parenting in response to the growing numbers of mothers using opioids is an opportunity for evaluation of different models of support [65]. Most current home-visiting programs target only families with children and only for a short period of time. Our findings suggest that interventions need to extend past the early postpartum period and include mothers whose children may not return home. Finally, the findings from our study also indicate that studies and surveillance focused only on deaths directly attributable to opioid overdose will miss the full extent of the problem, given the importance of deaths due to unintentional and intentional injuries, not all of which involve opioids.

In conclusion, while much attention and research on NAS has focused on infant and child outcomes in isolation, our study is the first population-based analysis to our knowledge of long-term maternal mortality following the birth of an infant with NAS. The findings provide a stark reminder of the vulnerability and sustained poor outcomes of these mothers. Policy responses to the current opioid epidemic will require effective strategies for risk mitigation and ongoing support for families affected by opioid use. Large-scale linkage of health and social care administrative data would facilitate ongoing research, program evaluation, and surveillance.

Supporting information

S1 Table [docx]
Study diagnostic codes and description of sociodemographic characteristics.

S2 Table [docx]
Description of baseline sociodemographic characteristics for neighbourhood income quintile and urban and rural area of residence.

S1 Text [docx]
Dataset creation plan and analytic plan: Long-term mortality in mothers of infants with neonatal abstinence syndrome: A population-based parallel-cohort study in England and Ontario, Canada.

S2 Text [docx]
STROBE Statement—Checklist of items that should be included in reports of cohort studies.