Maternal prescribed opioid analgesic use during pregnancy and associations with adverse birth outcomes: A population-based study
Authors:
Ayesha C. Sujan aff001; Patrick D. Quinn aff002; Martin E. Rickert aff001; Kelsey K. Wiggs aff001; Paul Lichtenstein aff003; Henrik Larsson aff003; Catarina Almqvist aff003; A. Sara Öberg aff003; Brian M. D’Onofrio aff001
Authors place of work:
Department of Psychological & Brain Sciences, Indiana University Bloomington, Bloomington, Indiana, United States of America
aff001; Department of Applied Health Science, School of Public Health, Indiana University Bloomington, Bloomington, Indiana, United States of America
aff002; Department of Medical Epidemiology and Biostatistics, Karolinska Institutet, Stockholm, Sweden
aff003; School of Medical Sciences, Örebro University, Örebro, Sweden
aff004; Astrid Lindgren Children’s Hospital, Karolinska University Hospital, Stockholm, Sweden
aff005; Department of Epidemiology, T.H. Chan School of Public Health, Harvard, Boston, United States of America
aff006
Published in the journal:
Maternal prescribed opioid analgesic use during pregnancy and associations with adverse birth outcomes: A population-based study. PLoS Med 16(12): e32767. doi:10.1371/journal.pmed.1002980
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pmed.1002980
Summary
Background
Published research on prescribed opioid analgesic (POA) use during pregnancy and birth outcomes is limited in scope and has not adequately adjusted for potential confounding factors. To help address these gaps, we estimated associations between maternal POAs during pregnancy and two adverse birth outcomes using a large population-based dataset, multiple definitions of POA exposure, and several methods to evaluate the influence of both measured and unmeasured confounding factors.
Methods and findings
We obtained data by linking information from several Swedish registers and conducted a retrospective cohort study on a population-based sample of 620,458 Swedish births occurring between 2007 and 2013 (48.6% female; 44.4% firstborn). We evaluated associations between prenatal POA exposure and risk for preterm birth (PTB; <37 gestational weeks) and small for gestational age (SGA; birth weight 2 standard deviations below the expected weight for gestational age or lower). We evaluated the influence of confounding by adjusting for a wide range of measured covariates while comparing exposed and unexposed infants. Additionally, we adjusted for unmeasured confounding factors by using several advanced epidemiological designs. Infants exposed to POAs anytime during pregnancy were at increased risk for PTB compared with unexposed infants (6.4% exposed versus 4.4% unexposed; adjusted odds ratio [OR] = 1.38, 95% confidence interval [CI] 1.31–1.45, p < 0.001). This association was attenuated when we compared POA-exposed infants with acetaminophen-exposed infants (OR = 1.18, 95% CI 1.07–1.30, p < 0.001), infants born to women who used POAs before pregnancy only (OR = 1.05, 95% CI 0.96–1.14, p = 0.27), and unexposed siblings (OR = 0.99, 95% CI 0.85–1.14, p = 0.92). We also evaluated associations with short-term versus persistent POA use during pregnancy and observed a similar pattern of results, although the magnitudes of associations with persistent exposure were larger than associations with any use or short-term use. Although short-term use was not associated with SGA (adjusted ORsingle-trimester = 0.95, 95% CI 0.87–1.04, p = 0.29), persistent use was associated with increased risk for SGA (adjusted ORmultiple-trimester = 1.40, 95% CI 1.17–1.67, p < 0.001) compared with unexposed infants. The association with persistent exposure was attenuated when we used alternative comparison groups (e.g., sibling comparison OR = 1.22, 95% CI 0.60–2.48, p = 0.58). Of note, our study had limitations, including potential bias from exposure misclassification, an inability to adjust for all sources of confounding, and uncertainty regarding generalizability to countries outside of Sweden.
Conclusions
Our results suggested that observed associations between POA use during pregnancy and risk of PTB and SGA were largely due to unmeasured confounding factors, although we could not rule out small independent associations, particularly for persistent POA use during pregnancy.
Keywords:
Birth – pregnancy – Drug therapy – Analgesics – opioids – Infants – Sweden – Preterm birth
Introduction
While estimates of the prevalence of opioid prescriptions among pregnant women vary across studies, research suggests that a substantial proportion of women are treated with opioids during pregnancy [1–5]. For example, a study of pregnant women enrolled in Medicare in the United States reported that 30% of pregnant women fill prescriptions for opioids [5]. Given that opioids cross the human placental barrier, maternal use of these medications results in fetal exposure [6]. However, effects of prenatal exposure to both illicit and prescription opioids on child development are unclear [7–9]. The existing research on opioid use during pregnancy has largely focused on illicit opioids (e.g., heroin) or buprenorphine and methadone in the context of medication-assisted treatment for opioid use disorder [10–14]. However, the use of prescribed opioid analgesics (POAs) for treatment of pain is much more common than illicit use of opioids [15] or medication-assisted treatment [16]. For example, approximately 5% of women use an illicit substance during pregnancy in the US [15], and less than 3% of women prescribed opioids are prescribed these medications to treat opioid use disorder [16].
The limited research on POA use during pregnancy has focused on birth outcomes, particularly preterm birth (PTB), fetal growth, and structural birth defects [7–9], all of which have potential consequences for future morbidity and mortality [17–19]. Some observational studies have reported statistically significant associations between prenatal POA exposure and birth outcomes, whereas others have not [7–9]. Moreover, it is unclear if observed associations are due to causal effects of POA exposure or confounding from indications for maternal analgesic use (e.g., traumatic injury, acute or chronic inflammation, musculoskeletal or neuropathic pain) or other patient characteristics (e.g., psychiatric disorders, concurrent use of other psychiatric medication) [16,20–25]. To help address these gaps, we estimated associations between maternal POAs during pregnancy and PTB and small for gestational age (SGA; a proxy for fetal growth restriction) in a large, population-based dataset using multiple definitions of POA exposure and several methods to evaluate the influence of both measured and unmeasured confounding factors. To our knowledge, no previous study has used methods to evaluate the influence of unmeasured confounding, which is critical because solely adjusting for measured characteristics is unlikely to adequately capture all sources of confounding [26,27]. Because no single observational method can account for all plausible confounding factors, we sought converging evidence from multiple comparisons to rigorously test causal hypotheses [28–30]. Specifically, we compared POA-exposed infants with infants exposed to acetaminophen (i.e., paracetamol), infants born to mothers with POA prescriptions before but not during pregnancy, and their unexposed siblings.
Methods
We conducted a retrospective cohort, complete-case analysis study on a population-based sample of 620,458 Swedish births from July 1, 2007, to December 31, 2013, who were not exposed to opioids for the treatment of opioid use disorder. We processed and analyzed data using SAS 9.4 and STATA 15.1. A STROBE checklist (S1 Appendix) and a description of our planned analyses (S2 Appendix) are included in the supplemental materials.
The institutional review board at Indiana University and the regional ethical review board in Stockholm, Sweden, approved this study. The study used data available from national registers, and informed consent was not necessary.
We obtained data by linking information from several Swedish registers [31]. The Multi-Generation Register includes biological relationships for all individuals born from 1932 and residing in Sweden since 1961 [32]. The Swedish Prescribed Drug Register includes records of filled medication prescriptions since July 2005 [33,34]. The Medical Birth Register includes information on 96%–99% of births and pregnancy characteristics since 1973 [35–37]. The National Patient Register includes diagnostic codes from all hospital admissions since 1987 and 80% of specialist outpatient care since 2001 [38,39]. The National Crime Register includes criminal convictions since 1973 [40,41]. The Education Register includes highest level of completed formal education [42]. The integrated database for labor market research includes annual socioeconomic data for all individuals since 1990 [43].
Measures
POA exposures
S3 Appendix provides detailed information on the included POAs, including frequencies of prescriptions of specific POA medications and information on the type of clinic prescribing the medications. For example, while many prescriptions originated from obstetrics/gynecology and maternity care clinics, the majority of prescriptions originated from other types of clinics.
We included all prescriptions with Anatomical Therapeutic Chemical (ATC) codes beginning with N02A, acetaminophen/codeine combination medications, with ATC code N02BE51, and buprenorphine (N02AE01, N07BC01, N07BC51) and methadone (N07BC02) prescribed as POAs. We distinguished buprenorphine and methadone for pain treatment versus opioid use disorder treatment with criteria used in previous publications [44,45].
We considered several exposure definitions using maternal POA prescriptions during nonoverlapping windows of time (Fig A in S3 Appendix). We defined these windows relative to approximated last menstrual period (LMP) and conception dates. Consistent with prior research [46], we first estimated LMP by subtracting gestational age (predominantly based on ultrasound measurements occurring the 18th to 20th week of pregnancy) from birth date and then estimated conception date by adding 14 days to LMP. Gestational age was predominantly based on ultrasound measurements given that approximately 95% of all pregnant women in Sweden undergo at least one ultrasound [47].
Before-pregnancy-only exposure included those with a prescription 360 days before conception to 91 days before conception and no prescriptions 90 days before conception to birth.
Exposure during a washout period only (i.e., the period in which there would be ambiguity regarding whether a filled prescription would lead to use during pregnancy) included those with a prescription 90 days before conception to the day before conception and no prescriptions from conception to birth. We chose a 90-day window for the washout period because in Sweden the maximum amount of opioid medication dispensed is for a three-month period [48]. Including the washout-period-only variable as a predictor in models prevented individuals who only had filled prescriptions that occurred shortly before pregnancy from being classified as either exposed or unexposed during pregnancy.
We defined during-pregnancy exposure as filled prescriptions anytime from conception to birth. We also used the timing of maternal prescriptions across trimesters to create proxies for short-term use (occurring in a single trimester) versus persistent use (occurring in multiple trimesters). Any use was the primary exposure, and secondary analyses examined whether associations differed by duration of use. We did not consider timing of exposure as defined by use in specific trimesters because preliminary analyses did not provide support for sensitive periods of exposure during pregnancy (S4 Appendix). Specifically, these preliminary analyses showed that while the magnitude of the associations with first-trimester exposure only (conception to 89 days after LMP) were slightly larger than the magnitude of associations with second-trimester (90 to 179 days after LMP) or third-trimester (180 days after LMP to birth) exposure, the associations were not statistically significantly different for both PTB (p = 0.25) and SGA (p = 0.54).
Active comparator medication exposure
We considered exposures to maternal acetaminophen (N02BE01) as prescriptions (1) anytime during pregnancy, (2) in a single trimester, and (3) in multiple trimesters. For all three exposures, we excluded individuals with maternal POA prescriptions anytime during pregnancy in order to capture acetaminophen-only exposure.
Covariates
Pregnancy-related characteristics included birth order, year of birth, maternal-reported smoking during the first trimester, and maternal prescriptions of other psychiatric medications during pregnancy (Table A in S3 Appendix includes ATC codes).
We included maternal and paternal characteristics that were associated with subsequent receipt of POAs [20,49]. Specifically, maternal and paternal characteristics included pre-conception inpatient and outpatient diagnoses of opioid use disorder, non-opioid substance use disorder, schizophrenia or bipolar disorder, and definite or uncertain suicide attempts made according to International Classification of Disease (ICD) criteria; any pre-conception conviction of violent, nonviolent, drug or alcohol, or driving-related crimes; age at childbirth; highest level of education at year of birth; and country of origin.
Other familial and socioeconomic characteristics included maternal reports of parents cohabitating at birth, family-level income quintile at year of birth relative to the Swedish population in that year, and neighborhood deprivation score quintile at year of birth relative to the Swedish population in that year. Neighborhood deprivation score was based on a principal component analysis of several yearly indicators for small geographical areas constructed to delineate socially homogenous areas. The deprivation score for a given area incorporated the proportion of welfare recipients, unemployed individuals, immigrants, divorced individuals, and individuals with low educational attainment, as well as measures of residential mobility, crime rates, and neighborhood disposable income [50].
Outcomes
Outcomes were PTB (birth before 37 gestational weeks) and SGA (birth weight 2 standard deviations below the expected fetal weight for gestational age or lower).
Associations between POA exposure and birth outcomes
We fit five logistic regression models estimating associations between POA exposure and birth outcomes in the analytic sample using robust standard errors to account for clustering of individuals within nuclear families bound by the same biological mothers (i.e., siblings). In all models we included washout-period-only POA exposure as a dummy code.
We first fit all models using POA exposure anytime during pregnancy as the predictor. Then, we re-estimated the models using single-trimester and multiple-trimester exposure as predictors. S5 Appendix lists predictors and comparison groups for the main analytic models.
Population-wide associations comparing exposed with all unexposed infants
Models 1 and 2 estimated population-wide associations with POA exposure compared with unexposed infants. Model 1 (unadjusted) did not include any covariates. Model 2 (adjusted) included all measured characteristics as covariates in the regression models. We did not use propensity scores to adjust for measured characteristics because with a relatively limited number of covariates, propensity score methods make exactly as much adjustment as more traditional, outcome regression methods [51], and prominent researchers have expressed serious concerns about the causal inferences researchers have drawn based on propensity score methods (e.g., [51,52]). For the purpose of transparency, all parameter estimates from model 2 are listed in S6 Appendix.
Associations using alternative comparison groups
Models 3 through 5 used alternative comparison groups consisting of infants likely to share some maternal characteristics with exposed infants. Model 3 (comparative safety) [53] compared POA-exposed infants with acetaminophen-only-exposed infants to assess the relative safety of the two medications. Given that both POAs and acetaminophen are prescribed for the treatment of pain, a null difference between POA-exposed and acetaminophen-exposed infants would indicate that POA use is as safe as a commonly used pain medication and provide support for confounding by indication.
Model 4 (before-pregnancy use) [54] compared infants exposed to POAs during pregnancy with infants born to women with POA prescriptions before but not during pregnancy. By design, model 4 accounted for all unmeasured factors shared by women who were prescribed POAs around the time of pregnancy. No association between during-pregnancy exposure compared with exposure before pregnancy only would suggest that shared factors rather than a causal influence of intrauterine POA exposure is responsible for an observed increased occurrence in adverse birth outcomes among infants born to women who filled POA prescriptions during pregnancy.
Model 5 (sibling comparison) [55] compared POA-exposed infants with their siblings who were unexposed to POAs during pregnancy. By design, this within-family comparisons accounts for all unmeasured genetic and environmental factors that make siblings similar, including all familial factors that were stable across pregnancies. A lack of association in a sibling comparison would suggest that familial risk factors rather than a causal effect of POA exposure during pregnancy explain observed population-wide associations. In addition to adjusting for potential unmeasured confounding factors by design (i.e., factors shared by the comparison groups), models 3 through 5 also included all measured covariates that varied between comparison groups.
Sensitivity analyses
We also conducted several sensitivity analyses to evaluate our modeling assumptions and the extent to which the exposure and outcome definitions influenced the main analyses results (S7 to S13 Appendix).
Bias from exposure definitions
First, to explore potential bias from exposure misclassification, we estimated population associations while adjusting for all the covariates using several different exposure definitions: (1) an expanded exposure window that included the 90 days before conception in case a prescription that was filled shortly before pregnancy was used during pregnancy, (2) a restricted exposure window that excluded the three days before birth in case a woman was prescribed POAs at the end of her pregnancy to use after delivery, and (3) an exposure defined according to filled prescriptions or maternal-reported use in order to capture women who filled POA prescription before pregnancy but used them during pregnancy (S7 Appendix).
Second, to evaluate whether a specific type of POA medication was largely driving our findings, we re-estimated adjusted associations in (1) a subsample excluding births occurring to women with during-pregnancy prescriptions of dextropropoxyphene (N02AC04 and N02AC54) because this medication is no longer prescribed in Sweden and (2) a subsample excluding births occurring to women with during-pregnancy methadone or buprenorphine prescriptions in case these infants were exposed to prescribed opioids for the treatment of opioid use disorder rather than pain (S8 Appendix).
Third, in order to evaluate whether exposure to medications other than POAs that are included in POA combination medications influenced the main analyses result, we estimated adjusted associations in a subsample excluding births occurring to mothers with filled prescriptions of combination POA medications (i.e., oxycodone/naloxone, buprenorphine/naloxone, morphine/antispasmodics, ketobemidone/antispasmodics, hydromorphone/antispasmodics, codeine combinations, and dextropropoxyphene combinations; S9 Appendix).
Fourth, to further assess whether exposure to polypharmacy was responsible for observed associations with POA exposure, we re-estimated adjusted associations in a subsample excluding births occurring to women prescribed other psychiatric medications during pregnancy (S10 Appendix).
Sibling comparison assumptions
We conducted two sets of analyses to evaluate the assumptions of sibling comparisons (S11 Appendix). To evaluate whether individuals with siblings differed from the full population, we re-estimated adjusted associations in a subsample comprised of the women that contributed more than once to the full population sample (i.e., women that had more than one singleton birth in the study period). To evaluate if the sibling comparisons were biased by carryover effects (i.e., exposure in an earlier pregnancy affecting subsequent pregnancies), we compared differentially exposed pairs of firstborn cousins.
Bias from outcome definitions
To evaluate whether any potential influence of POA exposure was not reflected in the clinical cutoff values of the birth outcomes, we fit all of the main models predicting birth outcomes on continuous scales. Specifically, we predicated (1) gestational age measured in days and (2) birth weight measured in grams adjusted for gestational age (S12 Appendix).
Influence of missing data
We conducted two sensitivity analyses to evaluate whether missing data influenced our findings (S13 Appendix). First, we estimated the association between absence of data and POA exposure in the target sample.
Then, we assessed the potential confounding influence of the covariates with missing data by removing these covariates from the fully adjusted model in the analytic sample and fitting this alternative model in both the analytic and target samples to evaluate whether the findings in the restricted (analytic) sample appear to represent the findings of the full (target) sample.
Results
Participants
We started with a population-based sample of 711,986 births occurring in Sweden between July 1, 2007, and December 31, 2013. To create the target sample, we sequentially excluded infants with invalid child identifiers (2,648), maternal identifiers (107), and sex (6); multiple births (19,844); births missing gestational age (173); and births exposed to buprenorphine or methadone for opioid use disorder treatment rather than pain (276). To create the analytic sample from the target sample, we excluded 68,474 infants with missing covariate data. The final analytic cohort of 620,458 births represented approximately 90% of the target sample. We also identified 288,995 births occurring to the same mother in the analytic sample for inclusion in sibling comparison analyses. This sample of siblings that shared mothers included 9,201 unique mothers.
Demographics
Table 1 compares background characteristics among infants exposed to any POA during pregnancy and all the unexposed infants in the target sample. Additionally, S14 Appendix shows the prevalence of background characteristics stratified by all exposure groups, including infants with single- and multiple-trimester exposure. We provided demographic data from the target sample in order to document the prevalence of missing data on covariates, which ranged from none to approximately four percent.
In the final analytic sample, 4.4% of infants were exposed to POAs anytime during pregnancy, with 3.7% exposed in a single trimester and 0.7% exposed in multiple trimesters (Table 2).
PTB occurred more often among exposed (6.4%) than unexposed (4.4%) infants, particularly among infants with multiple-trimester exposure (5.8% of single-trimester-exposed, 9.6% of multiple-trimester-exposed); and about 2% of infants were SGA regardless of exposure status (2.1% of unexposed, 2.2% of exposed, 2.0% of single-trimester-exposed, and 3.04% of multiple-trimester-exposed; S15 Appendix).
Associations with PTB
Exposure anytime during pregnancy
The unadjusted association between any exposure to POAs and PTB (odds ratio [OR] = 1.48, 95% confidence interval [CI] 1.41–1.56, p < 0.001; Table 3) was attenuated but remained robust after adjustment for measured covariates (OR = 1.38, 95% CI 1.31–1.45, p < 0.001) and was further attenuated when we used the alternative comparison groups (comparative safety OR = 1.18, 95% CI 1.07–1.30, p < 0.001; before-pregnancy-use OR = 1.05, 95% CI 0.96–1.14, p = 0.27; siblings OR = 0.99, 95% CI 0.85–1.14, p = 0.92).
Short-term versus persistent exposure during pregnancy
We observed a similar pattern of attenuated associations with increasing confounder control also when separating single- and multiple-trimester exposure, although this also revealed that associations for multiple-trimester exposure were consistently more pronounced than for single-trimester exposure. The unadjusted associations (single-trimester OR = 1.34, 95% CI 1.26–1.41, p < 0.001 and multiple-trimester OR = 2.31, 95% CI 2.08–2.55, p < 0.001) suggested a relationship with duration of use that remained despite attenuation in adjusted models (single-trimester OR = 1.27, 95% CI 1.20–1.34, p < 0.001 and multiple-trimester OR = 1.97, 95% CI 1.77–2.18, p < 0.001); however, the associations were partially attenuated when using some of the alternative comparison groups (comparative safety: single-trimester OR = 0.91, 95% CI 0.83–0.99, p = 0.03 and multiple-trimester OR = 1.53, 95% CI 1.36–1.72, p < 0.001; before-pregnancy-use: single-trimester OR = 1.00, 95% CI 0.92–1.09, p = 0.99 and multiple-trimester OR = 1.52, 95% CI 1.28–1.80, p < 0.001) and fully attenuated in the sibling comparison (single-trimester OR = 0.99, 95% CI 0.85–1.15, p = 0.87 and multiple-trimester OR = 1.04, 95% CI 0.70–1.55, p = 0.83).
Associations with SGA
Exposure anytime during pregnancy
Across models 1 through 5, associations were negligible (Table 3).
Short-term versus persistent exposure during pregnancy
Similar to the findings for PTB, the association between POAs exposure and SGA differed according to the duration of use (Table 3). While there was no association with single-trimester exposure (unadjusted OR = 0.96, 95% CI 0.88–1.06, p = 0.43), a moderate association with multiple-trimester exposure (unadjusted OR = 1.45, 95% CI 1.22–1.73, p < 0.001) remained in the adjusted (OR = 1.40, 95% CI 1.17–1.67, p < 0.001) and comparative safety (OR = 1.41, 95% CI 1.15–1.73, p = 0.001) models but was attenuated in the before-pregnancy-use (OR = 1.20, 95% CI 0.89–1.60, p = 0.23) and siblings (OR = 1.22, 95% CI 0.60–2.48, p = 0.58) models.
Sensitivity analyses
Bias from exposure definitions
S7 Appendix examined potential bias from exposure misclassification and found the same pattern of results with a number of alternative exposure definitions (i.e., the expanded prescription window definition, the restricted prescription window definition, and the definition according to maternal-reported POA use or during-pregnancy prescriptions). S8 Appendix suggested that findings were not largely driven by dextropropoxyphene, an opioid that is no longer prescribed, or by opioids that are also prescribed for the treatment of opioid use disorder. S9 Appendix showed commensurate adjusted associations in a subsample excluding those with during-pregnancy prescriptions of combination POA medications, suggesting that the results were not driven by inclusion of combination POA medications. Similarly, S10 Appendix showed comparable adjusted associations in a subsample excluding those with during-pregnancy polypharmacy, suggesting that associations with POA exposure observed in main analyses were not driven by exposure to polypharmacy.
Sibling comparison assumptions
Two analyses in S11 Appendix examined assumptions of sibling comparisons. The first analysis suggested that attenuation in sibling comparisons was not due to restricting the sample to infants with siblings, because population-wide associations were comparable in samples of infants with and without siblings. The second analysis suggested that sibling comparison findings were not due to carryover effects (i.e., exposure in a prior pregnancy affecting subsequent pregnancies) because comparisons of firstborn cousins showed the same pattern of results as the sibling comparison findings.
Bias from outcome definitions
S12 Appendix demonstrated that using clinical cutoff values for the outcomes did not cause a failure to detect an influence of POA exposure by estimating associations with two continuously measured birth outcomes. Any POA exposure was associated with reduced gestational age in unadjusted models, and the associations attenuated in subsequent models. POA exposure was not associated with reduced birth weight adjusted for gestational age across any of the models.
Influence of missing data
S13 Appendix suggests that excluding infants with missing data did not bias our results. Absence of data was not associated with prescriptions anytime during pregnancy after adjusting for all measured covariates in our target sample (OR = 1.03, 95% CI 0.98–1.08), indicating that covariate adjustment helped minimize potential bias from missing data. Moreover, removing the covariates with missing data (smoking during pregnancy, maternal and paternal education, maternal and paternal country of origin, paternal age at childbearing, parental cohabitation status at birth, family income, and neighborhood deprivation) from the fully adjusted model did not meaningfully change the association in the analytic sample, and this alternative model yielded similar estimates of association in the target sample as the analytic sample.
Discussion
In a population-based sample of Swedish births occurring between 2007 and 2013, we used multiple observational designs to evaluate the consequences of prenatal POA exposure on the risk of two adverse birth outcomes—PTB and SGA. Compared with unexposed infants, infants exposed to POAs anytime during pregnancy, in a single trimester, and in multiple trimesters were all at increased risk for PTB, although infants exposed in multiple trimesters had the greatest risk. However, when we used comparison groups consisting of unexposed infants that shared characteristics with exposed infants, the associations with PTB were attenuated, suggesting that observed associations were largely due to unmeasured confounding factors. For example, siblings exposed anytime during pregnancy, in a single trimester, and in multiple trimesters were not at increased PTB risk compared with their unexposed siblings, which suggests that unmeasured genetic and environmental factors shared by siblings account for observed population-wide associations.
For SGA, we only observed associations with persistent POA exposure (i.e., in multiple trimesters). We observed the higher risk of SGA among exposed infants when we made comparisons to unexposed infants, as well as acetaminophen-exposed infants. However, the association was attenuated when we used infants born to mothers with POA prescriptions before pregnancy only and unexposed siblings as the comparison groups, which again suggests that the observed associations were largely due to confounding factors.
Our study had several noteworthy strengths that distinguish it from previous studies. First, we used multiple methods that were able to account for both measured and unmeasured sources of confounding and found converging evidence across these methods that suggested that the observed associations with birth outcomes were at least partially due to confounding. This converging evidence suggests that our conclusions are not due to violations of the assumptions of one design. Second, we reduced the potential influence from exposure misclassification by including a 90-day prepregnancy washout period and conducting multiple sensitivity analyses to evaluate the influence of potential exposure misclassification. Third, unlike most previous studies that examined any POA use anytime during pregnancy, we defined exposure by the number of trimesters with prescriptions and the number of prescriptions throughout pregnancy, as well as evaluated for sensitive periods of exposure by estimating associations with exposure earlier and later in pregnancy. Fourth, we excluded opioids used to treat opioid use disorder. Most previous studies focused on two opioids—methadone and buprenorphine—in the context of medication-assisted treatment of opioid use disorder or did not differentiate between POAs and opioids prescribed for opioid use disorder treatment. This is an important distinction because women who are prescribed opioids to manage pain presumably differ on several important background characteristics compared with women who are prescribed opioids for the treatment of opioid use disorder.
Our study also had several limitations. First, our exposure was subject to misclassification because mothers may not have taken their prescribed medications. Second, because of the rarity of the exposures, we did not evaluate associations with specific POA medications, and we did not examine pill dosages, number of pills prescribed, and number days covered. Therefore, future research should utilize finer-grained measurements, such as morphine equivalent units, which take into account the number of pills prescribed, numbers of days supplied, strength per pill, and type of POA [56]. Third, the comparative safety model evaluated the relative safety of POA exposure compared with acetaminophen exposure, and a null difference could have been due to adverse effects of both medications rather than confounding by common indications for the medications. Fourth, we did not have measures of indications for POA use, although we used several methods to assess the influence of unmeasured confounding factors. Nonetheless, future research should evaluate associations with POA treatment during pregnancy among women with conditions causing chronic pain (e.g., back pain, abdominal pain, fibromyalgia, and rheumatoid arthritis) [2] in order to compare exposed to unexposed with similar conditions for which POA treatment is indicated. Fifth, we do not know if unmeasured confounding biased results. In particular, we cannot rule out the possibility that further unmeasured confounding factors account for the remaining observed associations, particularly with multiple-trimester exposure. Indeed, women who use POAs in multiple trimesters are likely to have more severe conditions requiring chronic POA treatment, which could result in greater unmeasured confounding by indication. Sixth, we do not know if our findings will generalize to countries outside of Sweden given between-country differences (e.g., higher prescribing rates in the US compared to Sweden) [2,57]; however, in theory these differences would not impact the ability to detect a causal effect of POA exposure. Nonetheless, future research should evaluate whether our findings apply to populations outside of Sweden. Seventh, we excluded infants that did not have information on all covariates. However, sensitivity analyses suggested that excluding these infants from the main analyses did not influence our conclusions (S13 Appendix). Eighth, we only explored two birth outcomes and were not able to explore more rare birth outcomes, such as extreme PTB. However, sensitivity analyses showed the same pattern of results with continuous measures of gestational age and fetal growth (S12 Appendix). Future research should examine the consequences of prenatal POA exposure on other important adverse birth outcomes (e.g., birth defects), as well as outcomes presenting later in development (e.g., autism and attention-deficit hyperactivity disorder).
Despite the limitations, our study has important clinical implications. The vast majority of infants exposed during pregnancy did not have PTB and were not SGA. The absolute risk was low even among infants with persistent exposure, with 90% not having PTB and 97% not being SGA. Although we could not rule out small independent associations, particularly for persistent exposure during pregnancy, our results suggested that associations between prenatal POA exposure and PTB and SGA were largely due to confounding factors associated with maternal POA use during pregnancy rather than causal effects of POA exposure. We believe our findings are valuable as they may help doctors and patients better weigh the risks and benefits of POA use in women of childbearing age and pregnant women, although decision-making must consider a wide range of potential adverse outcomes. Our results also indicate that women of childbearing years should be assessed for a broad range of risk factors, and interventions aimed at reducing the incidence of adverse birth outcomes associated with maternal POA use during pregnancy should target co-occurring risk factors.
Supporting information
S1 Appendix [docx]
STROBE checklist.
S2 Appendix [docx]
Planned analyses.
S3 Appendix [docx]
Additional information on medication exposures.
S4 Appendix [docx]
Preliminary analyses exploring for sensitive periods of exposure.
S5 Appendix [docx]
Review of models used to evaluate associations between POA use and birth outcomes.
S6 Appendix [docx]
All parameter estimates from adjusted models estimating associations with maternal prescribed opioid filled prescriptions anytime during pregnancy.
S7 Appendix [docx]
Sensitivity analyses evaluating for potential bias from exposure misclassification.
S8 Appendix [docx]
Sensitivity analyses evaluating the influence of type of opioid.
S9 Appendix [docx]
Sensitivity analyses evaluating the influence of the inclusion of combination POA medications.
S10 Appendix [docx]
Sensitivity analyses evaluating the role of polypharmacy.
S11 Appendix [docx]
Sensitivity analyses evaluating assumptions of sibling comparisons results.
S12 Appendix [docx]
Sensitivity analyses evaluating associations with continuous outcomes.
S13 Appendix [docx]
Sensitivity analyses evaluating the influence of missing data.
S14 Appendix [docx]
Prevalence of the covariates stratified by exposure status in the target sample.
S15 Appendix [docx]
Prevalence of PTB and SGA among all exposure and comparison groups.
Zdroje
1. Desai RJ, Hernandez-Diaz S, Bateman BT, Huybrechts KF. Increase in prescription opioid use during pregnancy among Medicaid-enrolled women. Obstet Gynecol. 2014;123(5):997–1002. Epub 2014/05/03. doi: 10.1097/AOG.0000000000000208 PubMed Central PMCID: PMC4020039. 24785852
2. Bateman BT, Hernandez-Diaz S, Rathmell JP, Seeger JD, Doherty M, Fischer MA, et al. Patterns of opioid utilization in pregnancy in a large cohort of commercial insurance beneficiaries in the United States. Anesthesiology. 2014;120(5):1216–24. Epub 2014/02/15. doi: 10.1097/ALN.0000000000000172 24525628.
3. Kroll-Desrosiers AR, Skanderson M, Bastian LA, Brandt CA, Haskell S, Kerns RD, et al. Receipt of Prescription Opioids in a National Sample of Pregnant Veterans Receiving Veterans Health Administration Care. Womens Health Issues. 2016;26(2):240–6. doi: 10.1016/j.whi.2015.09.010 WOS:000375058000016. 26549242
4. Engeland A, Bramness JG, Daltveit AK, Ronning M, Skurtveit S, Furu K. Prescription drug use among fathers and mothers before and during pregnancy. A population-based cohort study of 106,000 pregnancies in Norway 2004–2006. Br J Clin Pharmacol. 2008;65(5):653–60. Epub 2008/02/26. doi: 10.1111/j.1365-2125.2008.03102.x 18294334; PubMed Central PMCID: PMC2432474.
5. Epstein RA, Bobo WV, Martin PR, Morrow JA, Wang W, Chandrasekhar R, et al. Increasing pregnancy-related use of prescribed opioid analgesics. Annals of Epidemiology. 2013;23(8):498–503. doi: 10.1016/j.annepidem.2013.05.017 23889859
6. Malek A, Mattison DR. Drugs and Medicines in Pregnancy: The Placental Disposition of Opioids. Current Pharmaceutical Biotechnology. 2011;12(5):797–803. WOS:000290665000012. doi: 10.2174/138920111795470859 21480827
7. Lind JN, Interrante JD, Ailes EC, Gilboa SM, Khan S, Frey MT, et al. Maternal Use of Opioids During Pregnancy and Congenital Malformations: A Systematic Review. Pediatrics. 2017;139(6). Epub 2017/06/01. doi: 10.1542/peds.2016-4131 28562278; PubMed Central PMCID: PMC5561453.
8. Yazdy MM, Desai RJ, Brogly SB. Prescription Opioids in Pregnancy and Birth Outcomes: A Review of the Literature. J Pediatr Genet. 2015;4(2):56–70. Epub 2016/03/22. doi: 10.1055/s-0035-1556740 26998394; PubMed Central PMCID: PMC4795985.
9. Brennan MC, Rayburn WF. Counseling about risks of congenital anomalies from prescription opioids. Birth Defects Res A Clin Mol Teratol. 2012;94(8):620–5. Epub 2012/08/01. doi: 10.1002/bdra.23064 22847955.
10. Disher T, Gullickson C, Singh B, Cameron C, Boulos L, Beaubien L, et al. Pharmacological Treatments for Neonatal Abstinence Syndrome: A Systematic Review and Network Meta-analysisPharmacological Treatments for Neonatal Abstinence SyndromePharmacological Treatments for Neonatal Abstinence Syndrome. JAMA Pediatrics. 2019;173(3):234–43. doi: 10.1001/jamapediatrics.2018.5044 30667476
11. Reddy UM, Davis JM, Ren Z, Greene MF. Opioid Use in Pregnancy, Neonatal Abstinence Syndrome, and Childhood Outcomes: Executive Summary of a Joint Workshop by the Eunice Kennedy Shriver National Institute of Child Health and Human Development, American College of Obstetricians and Gynecologists, American Academy of Pediatrics, Society for Maternal-Fetal Medicine, Centers for Disease Control and Prevention, and the March of Dimes Foundation. Obstet Gynecol. 2017;130(1):10–28. Epub 2017/06/09. doi: 10.1097/AOG.0000000000002054 28594753; PubMed Central PMCID: PMC5486414.
12. Wachman EM, Schiff DM, Silverstein M. Neonatal Abstinence Syndrome Advances in Diagnosis and Treatment. JAMA-J Am Med Assoc. 2018;319(13):1362–74. doi: 10.1001/jama.2018.2640 WOS:000429035900020. 29614184
13. Brogly SB, Saia KA, Walley AY, Du HM, Sebastiani P. Prenatal buprenorphine versus methadone exposure and neonatal outcomes: systematic review and meta-analysis. Am J Epidemiol. 2014;180(7):673–86. Epub 2014/08/26. doi: 10.1093/aje/kwu190 25150272.
14. Grossman M, Seashore C, Holmes AV. Neonatal Abstinence Syndrome Management: A Review of Recent Evidence. Rev Recent Clin Trials. 2017;12(4):226–32. doi: 10.2174/1574887112666170816144818 WOS:000415749200002. 28814260
15. Creanga AA, Berg CJ, Ko JY, Farr SL, Tong VT, Bruce FC, et al. Maternal Mortality and Morbidity in the United States: Where Are We Now? Journal of Women's Health. 2014;23(1):3–9. doi: 10.1089/jwh.2013.4617 24383493
16. Patrick SW, Dudley J, Martin PR, Harrell FE, Warren MD, Hartmann KE, et al. Prescription opioid epidemic and infant outcomes. Pediatrics. 2015;135(5):842–50. doi: 10.1542/peds.2014-3299 25869370
17. Class QA, Rickert ME, Langstrom N, Lichtenstein P, D'Onofrio BM. Birth weight, physical morbidity, and mortality: A population-based sibling-comparison study. American Journal of Epidemiology. 2014;179:550–8. PubMed Central PMCID: PMCJournal in progress. doi: 10.1093/aje/kwt304 24355331
18. Class QA, Rickert ME, Larsson H, Lichtenstein P, D'Onofrio BM. Fetal growth and psychiatric and socioeconomic problems: population-based sibling comparison British Journal of Psychiatry. 2014;205(5):355–61. PubMed Central PMCID: PMCID in progress: NIHMSID: 602974. doi: 10.1192/bjp.bp.113.143693 25257067
19. D'Onofrio BM, Class QA, Rickert ME, Larsson H, Langstrom N, Lichtenstein P. Preterm birth and mortality and morbidity: A population-based quasi-experimental study. JAMA Psychiatry. 2013;70:1231–40. PubMed Central PMCID: PMC3823714. doi: 10.1001/jamapsychiatry.2013.2107 24068297
20. Quinn PD, Hur K, Chang Z, Krebs EE, Bair MJ, Scott EL, et al. Incident and long-term opioid therapy among patients with psychiatric conditions and medications: A national study of commercial health care claims. Pain. 2017;158(1):140–8. doi: 10.1097/j.pain.0000000000000730 27984526
21. Gustavsson A, Bjorkman J, Ljungcrantz C, Rhodin A, Rivano-Fischer M, Sjolund KF, et al. Pharmaceutical treatment patterns for patients with a diagnosis related to chronic pain initiating a slow-release strong opioid treatment in Sweden. Pain. 2012;153(12):2325–31. Epub 2012/09/05. doi: 10.1016/j.pain.2012.07.011 22944610
22. Sun EC, Dixit A, Humphreys K, Darnall BD, Baker LC, Mackey S. Association between concurrent use of prescription opioids and benzodiazepines and overdose: Retrospective analysis. BMJ. 2017;356:j760. doi: 10.1136/bmj.j760 28292769
23. Fredheim OMS, Borchgrevink PC, Mahic M, Skurtveit S. A pharmacoepidemiological cohort study of subjects starting strong opioids for nonmalignant pain: A study from the Norwegian Prescription Database. Pain. 2013;154(11):2487–93. doi: 10.1016/j.pain.2013.07.033 24075311
24. Larochelle MR, Zhang F, Ross-Degnan D, Wharam JF. Trends in opioid prescribing and co-prescribing of sedative hypnotics for acute and chronic musculoskeletal pain: 2001–2010. Pharmacoepidemiology and Drug Safety. 2015;24(8):885–92. doi: 10.1002/pds.3776 25906971
25. Kobus AM, Smith DH, Morasco BJ, Johnson ES, Yang X, Petrik AF, et al. Correlates of higher-dose opioid medication use for low back pain in primary care. Journal of Pain. 2012;13(11):1131–8. doi: 10.1016/j.jpain.2012.09.003 23117108
26. Morales DR, Slattery J, Evans S, Kurz X. Antidepressant use during pregnancy and risk of autism spectrum disorder and attention deficit hyperactivity disorder: systematic review of observational studies and methodological considerations. BMC Medicine. 2018;16(1):6. doi: 10.1186/s12916-017-0993-3 29332605
27. Pourhoseingholi MA, Baghestani AR, Vahedi M. How to control confounding effects by statistical analysis. Gastroenterology and Hepatology From Bed to Bench. 2012;5(2):79–83. PMC4017459. 24834204
28. Rutter M. Psychosocial influences: Critiques, findings, and research needs. Development and Psychopathology. 2000;12:375–405. doi: 10.1017/s0954579400003072 11014744
29. Academy of Medical Sciences Working Group. Identifying the environmental causes of disease: How should we decide what to believe and when to take action? London: Academy of Medical Sciences; 2007.
30. Shadish WR, Cook TD, Campbell DT. Experimental and quasi-experimental designs for generalized causal inference. New York: Houghton Mifflin; 2002.
31. Statistics Sweden [cited 2019 Sep 23]. Available from: https://www.scb.se/en/.
32. Statistics Sweden. Multi-Generation Register 2009—A Description of Contents and Quality. Örebro, Sweden: 2010.
33. National Board of Health and Welfare. 4 years with the Prescribed Drug Register [in Swedish]. Stockholm: National Board of Health and Welfare, 2010.
34. Wettermark B, Hammar N, Fored CM, Leimanis A, Otterblad Olausson P, Bergman U, et al. The new Swedish Prescribed Drug Register—opportunities for pharmacoepidemiological research and experience from the first six months. Pharmacoepidemiology and drug safety. 2007;16(7):726–35. Epub 2006/08/10. doi: 10.1002/pds.1294 16897791.
35. National Board of Health and Welfare. The Swedish Medical Birth Register—A Summary of Content and Quality. Stockholm: National Board of Health and Welfare, 2003.
36. Cnattingius S, Ericson A, Gunnarskog J, Kallen B. A quality study of a medical birth registry. Scand J Soc Med. 1990;18(2):143–8. doi: 10.1177/140349489001800209 2367825.
37. National Board of Health and Welfare. Pregnancies, Deliveries and Newborn Infants—The Swedish Medical Birth Register 1973–2011, Assisted Reproduction, treatment 1991–2010 [in Swedish]. Stockholm: 2013.
38. National Board of Health and Welfare. Quality and Content in the Patient Register—Inpatient care discharges 1964–2007 and outpatient (specialized care) visits 1997–2007 [in Swedish]. Stockholm: National Board of Health and Welfare, 2009 2009-125–15.
39. WHO. International Classification of Diseases, World Health Organization WHO. Geneva1992.
40. BRÅ e. Kriminalstatistik 2010 [Criminal statistics]. In: Västerås, editor. 2011.
41. Fazel S, Grann M. The population impact of severe mental illness on violent crime. American Journal of Psychiatry. 2006;163(8):1397–403. doi: 10.1176/ajp.2006.163.8.1397 WOS:000239358800018. 16877653
42. Statistics Sweden. Education in Sweden 2009 2009 [updated February 21, 2014] [cited 2019 Sep 23]. Available from: http://www.scb.se/statistik/_publikationer/uf0527_2009a01_br_uf08br0901.pdf
43. LISA database [cited 2019 Sep 23]. Available from: http://www.scb.se/Pages/List____257743.aspx
44. WHO Collaborating Centre for Drug Statistics Methodology. Guidelines for ATC classification and DDD assignment 2014. Oslo: Author, 2013.
45. Abrahamsson T, Berge J, Ojehagen A, Hakansson A. Benzodiazepine, z-drug and pregabalin prescriptions and mortality among patients in opioid maintenance treatment-A nation-wide register-based open cohort study. Drug and alcohol dependence. 2017;174:58–64. Epub 2017/03/21. doi: 10.1016/j.drugalcdep.2017.01.013 28315808.
46. Hampp C, Gelperin K, Taylor L, Bird ST, Toh S, Sahin L, et al. Misclassification in Assessment of First Trimester In-utero Exposure to Drugs Used Proximally to Conception: the Example of Letrozole Utilization for Infertility Treatment. American Journal of Epidemiology. 2018;188(2):418–25. doi: 10.1093/aje/kwy237 30321259
47. Hagenfeldt K, Alton V, Axelsson O, Blennow M, Bojö F, Bygdeman M, et al. Routine ultrasound examination during pregnancy—summary in English. Stockholm: The Swedish Council on Health Technology Assessment (SBU). 1998.
48. Soderberg KC, Laflamme L, Moller J. Newly initiated opioid treatment and the risk of fall-related injuries. A nationwide, register-based, case-crossover study in Sweden. Cns Drugs. 2013;27(2):155–61. doi: 10.1007/s40263-013-0038-1 23345030
49. Quinn PD, Rickert ME, Franck J, Sariaslan A, Boersma K, Lichtenstein P, et al. Associations of mental health and family background with opioid analgesic therapy: a nationwide Swedish register-based study. Pain. 2019. Epub 2019/07/25. doi: 10.1097/j.pain.0000000000001643 31339870.
50. Sariaslan A, Langstrom N, D'Onofrio B, Hallqvist J, Franck J, Lichtenstein P. The impact of neighbourhood deprivation on adolescent violent criminality and substance misuse: a longitudinal, quasi-experimental study of the total Swedish population. International journal of epidemiology. 2013;42(4):1057–66. Epub 2013/09/26. doi: 10.1093/ije/dyt066 24062294; PubMed Central PMCID: PMC3780994.
51. Pearl J. Causality: Models, Reasoning and Inference. Cambridge: Cambridge University Press; 2000.
52. Luellen JK, Shadish WR, Clark MH. Propensity scores: an introduction and experimental test. Eval Rev. 2005;29(6):530–58. doi: 10.1177/0193841X05275596 16244051
53. Lund JL, Richardson DB, Stürmer T. The active comparator, new user study design in pharmacoepidemiology: Historical foundations and contemporary application. Current Epidemiology Reports. 2015;2(4):221–8. doi: 10.1007/s40471-015-0053-5 26954351
54. Smith GD. Assessing intrauterine influences on offspring health outcomes: can epidemiological studies yield robust findings? Basic Clin Pharmacol Toxicol. 2008;102(2):245–56. doi: 10.1111/j.1742-7843.2007.00191.x 18226080
55. D'Onofrio BM, Class QA, Rickert ME, Sujan AC, Larsson H, Kuja-Halkola R, et al. Translational Epidemiologic Approaches to Understanding the Consequences of Early-Life Exposures. Behav Genet. 2016;46(3):315–28. Epub 2015/11/23. doi: 10.1007/s10519-015-9769-8 26590988; PubMed Central PMCID: PMC4860044.
56. Svendsen K, Borchgrevink PC, Fredheim O, Hamunen K, Mellbye A, Dale O. Choosing the unit of measurement counts: The use of oral morphine equivalents in studies of opioid consumption is a useful addition to defined daily doses. Palliative Medicine. 2011;25(7):725–32. doi: 10.1177/0269216311398300 21378066
57. Handal M, Engeland A, Ronning M, Skurtveit S, Furu K. Use of prescribed opioid analgesics and co-medication with benzodiazepines in women before, during, and after pregnancy: a population-based cohort study. European Journal of Clinical Pharmacology. 2011;67(9):953–60. doi: 10.1007/s00228-011-1030-7 WOS:000293790000011. 21484468
Štítky
Interní lékařstvíČlánek vyšel v časopise
PLOS Medicine
2019 Číslo 12
- Příznivý vliv Armolipidu Plus na hladinu cholesterolu a zánětlivé parametry u pacientů s chronickým subklinickým zánětem
- Berberin: přírodní hypolipidemikum se slibnými výsledky
- Flexofytol® – přírodní revoluce v boji proti osteoartróze kloubů
- Vliv kombinace nutraceutik na remodelaci levé komory srdeční u osob s metabolickým syndromem
- Nutraceutika a jejich ovlivnění mírného kardiometabolického rizika
Nejčtenější v tomto čísle
- Ambient particulate matter pollution and adult hospital admissions for pneumonia in urban China: A national time series analysis for 2014 through 2017
- Association between gestational weight gain and severe adverse birth outcomes in Washington State, US: A population-based retrospective cohort study, 2004–2013
- Adherence to the 2017 French dietary guidelines and adult weight gain: A cohort study
- Acute kidney injury and adverse renal events in patients receiving SGLT2-inhibitors: A systematic review and meta-analysis