Key PointsQuestion
Does prescription opioid medication have a potentially causal role in the risk for depression and anxiety disorder?
Findings
In this 2-sample mendelian randomization study using genetic instruments for common pain medications, the genetic liability for prescription opioid use was associated with increased risk for major depression.
Meaning
While further work is needed, this genetics-based study supports conventional observational literature suggesting prescription opioid use increases the risk for depression.
Importance
Growing evidence suggests that prescription opioid use affects depression and anxiety disorders; however, observational studies are subject to confounding, making causal inference and determining the direction of these associations difficult.
Objective
To investigate the potential bidirectional associations between the genetic liability for prescription opioid and other nonopioid pain medications and both major depressive disorder (MDD) and anxiety and stress-related disorders (ASRD) using genetically based methods.
Design, Setting, and Participants
We performed 2-sample mendelian randomization (MR) using summary statistics from genome-wide association studies (GWAS) to assess potential associations of self-reported prescription opioid and nonopioid analgesics, including nonsteroidal anti-inflammatories (NSAIDs) and acetaminophen-like derivatives use with MDD and ASRD. The GWAS data were derived from participants of predominantly European ancestry included in observational cohorts. Data were analyzed February 20, 2020, to May 4, 2020.
Main Outcomes and Measures
Major depressive disorder, ASRD, and self-reported pain medications (opioids, NSAIDs, anilides, and salicylic acid).
Results
The GWAS data were derived from participants of predominantly European ancestry included in the population-based UK Biobank and Lundbeck Foundation Initiative for Integrative Psychiatric Research studies: approximately 54% of the initial UK Biobank sample and 55.6% of the Lundbeck Foundation Initiative for Integrative Psychiatric Research sample selected for the ASRD GWAS were women. In a combined sample size of 737 473 study participants, single-variable MR showed that genetic liability for increased prescription opioid use was associated with increased risk of both MDD (odds ratio [OR] per unit increase in log odds opioid use, 1.14; 95% CI, 1.06-1.22; P < .001) and ASRD (OR, 1.24; 95% CI, 1.07-1.44; P = .004). Using multivariable MR, these opioid use estimates remained after accounting for other nonopioid pain medications (MDD OR, 1.14; 95% CI, 1.04-1.25; P = .005; ASRD OR, 1.30; 95% CI, 1.08-1.46; P = .006), and in separate models, accounting for comorbid pain conditions. Bidirectional analyses showed that genetic liability for MDD but not ASRD was associated with increased prescription opioid use risk (OR, 1.18; 95% CI, 1.08-1.30; P < .001). These estimates were generally consistent across single-variable and multivariable inverse variance–weighted (MV-IVW) and MR-Egger sensitivity analyses. Pleiotropy-robust methods did not indicate bias in any MV-IVW estimates.
Conclusions and Relevance
The findings of this mendelian randomization analysis suggest evidence for potential causal associations between the genetic liability for increased prescription opioid use and the risk for MDD and ASRD. While replication studies are necessary, these findings may inform prevention and intervention strategies directed toward the opioid epidemic and depression.
The United States is in the middle of an opioid epidemic,1 with an approximately 5-fold increase in opioid prescription use over the past 20 years resulting in large increases in opioid misuse.2 Opioid-related deaths are rapidly increasing, and approximately 68% of the 702 000 drug overdose deaths in the United States from 2013 to 2017 involved opioids.3,4 It is also reported that prescription opioids (vs illicit opioids) are the first opioids to be misused,5 with almost 30% of patients prescribed opioids for chronic pain misusing them6 and about 12% developing opioid use disorder (OUD).6
Informed prescribing practices require comprehensive understanding of the treatment’s risks and benefits, and while opioids alleviate pain, chronic use is associated with numerous adverse effects, including immunosuppression, natural reward processes dysregulation, and neurohormonal deficits.7,8 Observational studies have also found opioid use to be highly comorbid with both major depressive disorder (MDD) and anxiety and stress-related disorders (ASRD).9-12 It is estimated that the approximately 7.8 million adults with psychiatric disorders (primarily MDD and ASRD) receive more than half of the almost 200 million yearly opioid prescriptions.10 Compared with other patient populations, individuals with psychiatric disorders are also more likely to report long-term opioid use.13-15
Because MDD and ASRD are leading global causes of disability and death,16,17 elucidating the direction and potentially causal effect of these associations would be useful to inform prevention strategies. Literature suggests prescription opioid use increases MDD risk,18-21 and while observational findings suggest a potential association between prescription opioid use, MDD, and ASRD, observational data are subject to confounding and reverse causation, making causal inference difficult.22,23 Therefore it remains to be elucidated whether MDD is a cause or consequence of prescription opioid use.18 While prospective randomized clinical trials (RCTs) are the criterion standard of causal inference,24 performing RCTs to evaluate the effects of prescription opioid use or other nonopioid pain medications is often complicated by preexisting psychiatric comorbidities.
Mendelian randomization (MR), which uses single-nucleotide variants (SNVs) as unconfounded proxies for exposures to estimate their effect on outcomes of interest, minimizes the bias affecting observational epidemiologic studies.23,25-27 Conceptually, MR has analogies with RCTs, with randomization occurring at meiosis,23,26 and is an important strategy for strengthening causal inference when RCTs are impractical or unethical.23 Given the public health importance of assessing the possible bidirectional associations between prescription opioid use and neuropsychiatric disorders, gaining causal inference into these associations would be important to aid prevention strategies.
In the absence of RCTs, we used a 2-sample MR study design of summary-level data on self-reported prescription opioid use in the UK Biobank (UKB),28 the largest genome-wide association study (GWAS) to date on MDD29 and the largest GWAS to date on ASRD,30 to conduct bidirectional MR analyses investigating potential causal associations between the genetic liability for opioid pain medications and MDD and ASRD. Given prescription opioid medications are given for pain, and pain increases the risk for both MDD and ASRD,31,32 we aimed to provide context to our primary prescription opioid analysis by including other common pain medications, including anilides, salicylic acid and derivatives, and nonsteroidal anti-inflammatory drugs (NSAIDs) and chronic pain conditions to the study. We also leveraged multivariable MR (MVMR) methods developed in 201933 to account for potential confounding owing to these pain medications and chronic pain conditions that may affect the use of prescription opioid analgesics and the risk for these psychiatric disorders.
Study Design and Data Sources
A detailed description of the methods used in this study is provided in the eMethods in the Supplement. We used publicly available summary statistics from 3 GWAS sources of predominantly European ancestry (Figure 1; eTable 1 in the Supplement). All studies have existing ethical permissions from their respective institutional review boards and include participant written informed consent and rigorous quality control. Because all analyses herein are based on publicly available summary data, no ethical approval from an institutional review board was required for this study. Data for this study were analyzed from February 20, 2020, to May 4, 2020.
We used summary statistics from the first medication use case-control GWAS conducted among UKB study participants to generate genetic instruments for opioid and nonopioid pain medications.28 Pain medication categories were classified by active ingredient using the Anatomical Therapeutic Chemical Classification System and then assigned to 23 categories by active ingredient, including opioids (eg, morphine, oxycodone, codeine, fentanyl, pethidine, and tramadol), NSAIDs, anilides, and salicylic acid products (eTable 2 in the Supplement). Approximately 54% of UKB study participants are women, and the mean (SD) age of UKB study participants on the first visit UKB assessment was 56.5 (8.1) years.28 We also used summary statistics from MRC-IEU UKB GWASs34 for 6 PHEnome Scan Analysis Tool (PHESANT)35 phenotypes related to site-specific (back, knee, hip, neck/shoulder, headache, and abdominal/stomach) pain occurring for more than 3 months.
We used summary statistics from the largest publicly available GWASs for our MDD and ASRD gene associations.29,30 The MDD cases were required to meet international consensus criteria (DSM-IV, International Classification of Diseases, Ninth Revision, or International Statistical Classification of Diseases and Related Health Problems, Tenth Revision) for a lifetime diagnosis of MDD (45 591 cases and 97 674 controls).29 The MDD GWAS meta-analysis29 incorporated 29 cohorts. Cohorts had a wide range of ages (18 to >80 years), and approximately 56.0% of the participants were women.29 The ASRD cases and controls (12 665 cases and 19 220 controls) were selected from the Danish Lundbeck Foundation Initiative for Integrative Psychiatric Research (iPSYCH) cohort, a representative sample from the population of Denmark born between May 1, 1981, and December 31, 2005.30 Approximately 55.6% of ASRD sample participants were women. The ASRD cases were assigned by a psychiatrist during routine clinical care according to the international Statistical Classification of Diseases and Related Health Problems (10th Revision).30 Both MDD and ASRD cohort participants were of European ancestry, and all GWASs included age and sex as covariates in the association analysis; population stratification is taken into account in the principal components analysis.28-30
Participant overlap in samples used to estimate genetic associations between exposures and outcomes can increase weak instrument bias in MR analyses.36,37 We avoided overlap in our analyses of pain medication use on MDD and ASRD: we used meta-analytic results for MDD that excluded the UKB cohort (n = 29 740). To improve precision for the bidirectional analyses of MDD on pain medication use, we used meta-analytic results for MDD including additional cohorts from 23andMe (n = 307 354) and UKB (overall n = 480 359). Hence, sample overlap (for bidirectional analyses) was minimal (6.2%), and because MDD instrument strength was considered strong (F statistic of approximately 32, see subsequent text), considerable weak instrument bias is not expected.37
All analyses were performed using the TwoSampleMR and MendelianRandomization R packages (The R Foundation). We selected all relevant SNVs identified in each GWASs as having reached selection threshold P less than 5 × 10−6 and being uncorrelated (10 000 kilobase pairs apart and R2 ≤.001). We obtained SNV effects and corresponding standard errors from the exposure and outcome GWASs.38 We harmonized exposure and outcome data, removed palindromic SNVs with intermediate allele frequencies, and estimated the F parameter to evaluate instrument strength.37 We applied Steiger filtering to the harmonized data to identify and remove those SNVs exhibiting reverse causation by the test metric. The observed variance of the outcome exceeded the observed variance of the exposure explained by the SNVs (eTables 3-8 in the Supplement: harmonized data sets).39
For single-variable MR analysis, we used inverse variance–weighted (IVW) MR as the primary method. However, because this method gives consistent estimates only if all genetic variants are valid instrumental variables (IVs), we considered complementary MR-Egger and weighted median-based regression methods, which make different IV assumptions, as a sensitivity analysis to address the question of robustness of our IVW estimate.23,38,40 The MR-Egger regression gives consistent estimates when 100% of genetic variants are invalid IVs; weighted median requires 50% of the weight to come from valid IVs. However, regarding efficiency, weighted median estimates generally are nearly as precise as IVW estimates; both are substantially more precise than MR-Egger estimates, with MR-Egger regression estimates particularly imprecise if all IVs have similar magnitudes of association with the exposure.41 We used the MR-Egger intercept test,42 Cochran Q heterogeneity test,43 and MR pleiotropy residual sum and outlier (MR-PRESSO) test44 to evaluate potential IV violations. We also performed leave-one-out analyses to detect high influence points.39
For MVMR analyses, we constructed instruments using SNVs in each of the GWASs meeting our single-variable MR selection criteria, described previously. We used the MVMR extension of the inverse-variance weighted MR method36 and MR-Egger method to correct for both measured and unmeasured pleiotropy.45 We combined the SNVs from the relevant GWASs: prescription opioid plus nonprescription pain medications; prescription opioid plus each of the chronic pain conditions; and the bidirectional MDD plus ASRD instrument. We removed duplicate and correlated SNVs (within 10 000 kilobase pairs; R2 ≥0.001), then extracted the SNV effects and corresponding standard errors from the exposures and outcome GWASs (eTables 9-14 in the Supplement: harmonized data sets).
The MR estimates are reported as odds ratios (ORs) interpreted as MDD or ASRD risk per unit increase in log odds of the opioid or nonopioid medication use (or alternatively, for bidirectional results, the medication use risk per unit increase in log odds of MDD or ASRD). We report confidence intervals for the main results. We index the strength of evidence against the null hypotheses (no association) by the exact P value. While we caution against interpreting study findings solely on the basis of a P-value threshold,46 we use a 2-sided α of .00625 based on testing 4 pain medications categories against 2 psychiatric outcomes as a heuristic allowing for follow-up of a plausible number of findings. In assessing consistency and robustness, we looked for estimates substantially agreeing in direction and magnitude (overlapping confidence intervals) across complementary MR methods.
Single-variable conventional MR analysis showed only genetically determined prescription opioid use having an effect estimate consistent with increased risk for both MDD (IVW OR, 1.14; 95% CI, 1.06-1.22; P < .001) and ASRD (IVW OR, 1.24; 95% CI, 1.07-1.44; P = .004) (Table 1 and Figure 2; eTables 15-18 in the Supplement). These estimates were broadly consistent with estimates from the weighted median and MR-Egger sensitivity analyses; the MR-Egger estimates were substantially less precise (MDD WMOR, 1.18; 95% CI, 1.07-1.30; P < .001; ASRD WM OR, 1.22; 95% CI, 1.00-1.49; P = .05; MDD MR-Egger OR, 1.17; 95% CI, 0.92-1.47; P = .21; ASRD MR Egger OR, 1.26; 95%CI, 0.78-2.02; P = .36). The MR-Egger intercept analysis did not indicate horizontal pleiotropy. Conventional IVW leave-one-out analysis did not identify any high leverage points with high influence (Figure 2; eFigures 1-5 in the Supplement). The F statistics for the genetic instruments were consistent with an absence of weak instrument bias.
In MVMR, assessing the genetic liabilities for prescription opioid and nonopioid pain medications use jointly, opioid use retained a robust relationship with both MDD (IVW OR, 1.14; 95% CI, 1.04-1.25; P = .005) and ASRD (IVW OR, 1.30; 95% CI, 1.08-1.56; P = .006) (Table 2; eTable 19 in the Supplement). These estimates were broadly consistent with estimates from the MVMR-Egger sensitivity analyses, although the MVMR-Egger estimates were again substantially less precise for MDD (OR, 1.13; 95% CI, 1.02-1.26; P = .02), and similarly, ASRD (OR, 1.18; 95% CI, 0.95-1.46; P = .13). The MR-Egger intercept analysis again did not indicate horizontal pleiotropy. In additional MVMR analyses, assessing the genetic liability for prescription opioid use jointly with site-specific pain, opioid use retained that association with both MDD and ASRD (eg, controlling for back pain: MDD IVW OR, 1.12; 95% CI, 1.05-1.21; P = .001; ASRD OR, 1.26; 95% CI, 1.09-1.45; P = .002), again, the MVMR-Egger estimates consistent but substantially less precise (MDD MR-Egger OR, 1.12; 95% CI, 0.95-1.30; P = .19; ASRD MR-Egger OR, 1.27; 95% CI, 0.90-1.80; P = .18). The MR-Egger intercept analysis again did not indicate horizontal pleiotropy (eTable 20 in the Supplement).
Bidirectional single-variable analyses showed the genetic liability for MDD, but not ASRD, having effect size estimates consistent with increased risk of use for all opioid and nonopioid pain medications (opioids OR, 1.18; 95% CI, 1.08-1.30; P < .001; anilides OR, 1.20; 95% CI, 1.00-1.28; P < .001; NSAIDs OR, 1.15; 95% CI, 1.09-1.22; P < .001; salicylic acid/derivatives OR, 1.10; 95% CI, 1.03-1.17; P = .002) (Table 3; eTables 15-18 in the Supplement). These estimates were broadly consistent with estimates from the weighted median and MR-Egger sensitivity analyses; the MR-Egger estimates were substantially less precise. The MR-Egger intercept analysis did not indicate horizontal pleiotropy. The F statistics for the genetic instruments were consistent with an absence of weak instrument bias.
We evaluated potential bidirectional associations between the genetic liability for prescription opioid and nonopioid pain medication use and both MDD and ASRD and found genetic evidence that prescription opioid use was associated with increased MDD and ASRD risk, while strikingly, nonopioid analgesics had no direct association with the risk for MDD or ASRD.
Opioid MDD MR estimates were consistent in magnitude and direction across IVW, weighted median, and MR-Egger analyses, with the MR-Egger estimate substantially less precise, as is typically expected in MR genetic association studies, and the MR Egger intercept terms, being close to zero, consistent with absence of pleiotropy.47 Our findings extend observational literature suggesting prescription opioid use increases the risk for MDD,19-21 and unlike 2020 evidence for a causal effect of ICD-defined OUD on MDD,48 our opioid use instrument only included prescription opioids and not heroin, suggesting important neuropsychiatric implications, beyond the risk for opioid misuse, and OUD, that may be considered when prescribing opioids.
The underlying mechanisms of prescription opioids in the pathophysiology of MDD remain to be elucidated, but potentially include opioid-induced dysregulation of reward circuitry that results in reduced reward perception or pleasure and relief generation49 or other physical medical dysregulation (ie, endocrine and autonomic nervous system abnormalities50) that potentially contributes to the physical symptoms of MDD. Notably, it has been suggested that the endogenous opioid system is directly involved in the regulation of mood and the dysregulation of that system may factor into depression and anxiety11,51,52; μ-opioid receptors (MOR) are widely distributed within the brain, including regions involved in emotion regulation,11 and κ-opioid receptors (KOR) are expressed in the cortex, striatum, hippocampus, amygdala, and thalamus,53 suggesting a role in reward, pain, and emotion. Opioid receptor antagonism has also been shown to attenuate ketamine’s antidepressant effects,54 and subeuphoric doses of partial opioid agonists have been shown to improve MDD symptoms.11 However, while small opioid doses may improve mood by activating MORs, prolonged opioid use may saturate MORs and activate the KOR, which also modulate mood55,56 and are associated with depression54,57; in rodents, KOR agonists increase anxiety while deficiencies in the KOR system decrease anxiety.51 In humans, prolonged opioid use (30 days or more), which may saturate the opioid receptor system to affect mood,54 increased risk for developing treatment-resistant depression by more than 25%, compared with opioid use for less than 30 days.21 Elucidating the underlying mechanisms of the opioid system dysregulation potentially shared between MDD and OUD may be important to combat these crises.11
Our bidirectional analyses, with genetic liability for prescription opioid use as an outcome, point to genetic liability for MDD but not ASRD as a possible causal risk factor of opioid use in psychiatric populations. Depression severity has been shown to be associated with increasing likelihood of misusing opioid medications for nonpain symptoms and self-increasing opioid dosage58 and up to 30% of long-term opioid users who have MDD qualify for moderate-to-severe OUD,15,59 suggesting that targeting opioid use prevention for patients with MDD may help mitigate the US opioid epidemic.1,10 Given more than half of individuals with OUD have comorbid MDD,10 the increased use of opioids by individuals with MDD may be owing to self-medication of social or emotional pain,11 suggesting that the development of therapeutic interventions with minimal risk targeted at endogenous opioid dysregulation represents another important prevention opportunity.11 In sum, our findings support recommendations that caution is needed with prescribing in settings of mood disorders in favor of nonopioid alternatives, with screening for MDD prior to initiating opioid treatment.10 Further, our findings that the genetic liability for MDD increased the risk of the genetic liability to take NSAIDs, anilide, and salicylic acid and derivatives support the well-known association between depression and physical pain; comorbid depression and pain experience reduce physical, mental, and social functioning when beyond either depression or pain.60 Further still, our findings that genetic liability for NSAIDs use increased the risk for MDD and ASRD supports observational evidence suggesting NSAIDs have adverse neural effects, including neuropsychiatric symptoms,61,62 although the association direction is opposite to reported RCTs using NSAIDs as adjunct therapies for depression.63 The NSAID-related adverse neuropsychiatric symptoms are most often reported with indomethacin and selective cyclooxygenase-2 inhibitors, which may modulate neural processes and synaptic signaling processes where cyclooxygenase-2 is localized.61
Strengths and Limitations
This innovative 2-sample MR study investigating the association between genetic liability for opioid and other pain medication use and neuropsychiatric outcomes has several strengths. We use summary genetic associations from the largest available GWASs, important genetic analyses investigating small effect sizes; generally, larger sample sizes increase measurement precision. We also use complementary 2-sample MR methods for sensitivity analysis.64 The MVMR models are a major strength, enabling us to account for potential confounding owing to chronic pain, and also to estimate the direct association for genetic liability for prescription opioid use. Relatedly, including nonopioid pain medications provided additional context to the prescription opioid findings because none of the nonopioid pain medications retained a direct association with our neuropsychiatric outcomes in MVMR analyses. Finally, the list of prescribed prescription opioids comprising our opioid variable excludes heroin and other illicit opioid substances, enabling us to evaluate the specific association between the genetic liability for prescribed opioids and MDD.
There are also limitations to this study. The pain medication use and psychiatric disorder SNV effect estimates were obtained in mostly European studies, thus minimizing the possibility of population stratification bias and increasing the plausibility of the 2-sample MR assumption that summary associations derived from comparable populations.39,65 We performed sensitivity analyses to assess and minimize heterogeneity and pleiotropy. Nonetheless, we emphasize the importance of triangulating multiple lines of experimental evidence to strengthen causal inference.66
Regarding genetic instrument selection, following the example of prior MR studies,67,68 we used a stringent selection threshold (P < 5 × 10−6) for the pain medication use and ASRD risk instruments to compensate for lack of SNVs with effect P values less than conventional genome-wide significance (P < 5 × 10−8). We also used stringent LD clumping thresholds to ensure instrument independence. The biologic mechanisms of the selected SNVs are unknown; however, sensitivity analyses failed to find evidence for horizontal pleiotropy. Further, each pain medication phenotype was a binary variable (use/nonuse), so we could not assess potential dose-dependent changes in risk associated with prescription opioid use. Also, the genetic variants for pain medication use may be related to underlying disease, pain conditions, or even subclinical levels of the traits that influence MDD risk,28 which may affect results, although we note the association of the genetic liability of opioid use found in single-variable MR persisted in MVMR accounting for genetic liability for site-specific chronic pain. Nonetheless, future genetic studies using detailed hospital-based information about pain medication and including negative control populations would possibly allow to further strengthen causal inference. Finally, as has been noted previously, the UKB cohorts may not represent the general UK population,69 and analyses were limited to individuals of European ancestry; caution is warranted before generalizing findings to other populations.
We provide preliminary genetic evidence that prescription opioid use increases MDD and ASRD risk, suggesting important clinical consequences. We also find genetic evidence that MDD is a potential causal risk factor for increased prescription opioid use, which may help identify patient populations to aim prevention strategies to curb the ongoing opioid epidemic.
Corresponding Author: Falk W. Lohoff, MD, National Institute on Alcohol Abuse and Alcoholism, National Institutes of Health, 10 Center Dr (10CRC/2-2352), Bethesda, MD 20892 (falk.lohoff@nih.gov).
Published Online: November 11, 2020. doi:10.1001/jamapsychiatry.2020.3554
Open Access: This is an open access article distributed under the terms of the CC-BY License. © 2020 Rosoff DB et al. JAMA Psychiatry.
Author Contributions: Mr Rosoff had full access to all of the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.
Concept and design: All authors.
Acquisition, analysis, or interpretation of data: Rosoff, Lohoff.
Drafting of the manuscript: Rosoff, Lohoff.
Critical revision of the manuscript for important intellectual content: All authors.
Statistical analysis: Rosoff.
Obtained funding: Lohoff.
Administrative, technical, or material support: Lohoff.
Supervision: Lohoff.
Conflict of Interest Disclosures: None reported.
Funding/Support: This work was supported by the National Institutes of Health intramural funding (ZIA-AA000242 to Dr Lohoff); Division of Intramural Clinical and Biological Research of the National Institute on Alcohol Abuse and Alcoholism.
Role of the Funder/Sponsor: The funding sources had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.
Additional Contributions: This research was made possible by previous studies from the Substance Use Disorders Working Group of the Psychiatric Genomics Consortium (PGC-SUD) (supported by funds from National Institute on Drug Abuse and National Institute of Mental Health to MH109532 and, previously, with analyst support from National Institute on Alcohol Abuse and Alcoholism to U01AA008401 [COGA]); the Medical Research Council Integrative Epidemiology Unit (MRC-IEU, University of Bristol, England), especially the developers of the MRC-IEU UK Biobank Genome-Wide Association Study Pipeline; and the Lundbeck Foundation Initiative for Integrative Psychiatric Research (iPSYCH). We acknowledge their contributing studies and the participants in those studies, without whom this effort would not be possible.
Additional Information: All analyses were conducted using publicly available data. Complete genetic instrument list and harmonized datasets for every analysis included in this study are available in the eAppendix in the Supplement. The publicly available summary-level data for opioid and nonopioid pain medication use is available at https://www.ebi.ac.uk/gwas/downloads/summary-statistics; for major depression/major depressive disorder, at https://www.med.unc.edu/pgc/results-and-downloads/mdd/ and through MR Base at https://www.mrbase.org/; and for anxiety and stress-related disorders, at https://ipsych.dk/forskning/downloads/. Code Availability: the analysis code in R is available on request, and all data displayed in the figures are available in the eTables in the Supplement.
1.US Centers for Disease Control and Prevention. Understanding the epidemic: drug overdose deaths in the United States continue to increase in 2015. Published 2011. Updated December 19, 2018. Accessed November 5, 2019.
https://www.cdc.gov/drugoverdose/epidemic/ 2.Stoicea
N, Costa
A, Periel
L, Uribe
A, Weaver
T, Bergese
SD. Current perspectives on the opioid crisis in the US healthcare system: a comprehensive literature review.
Medicine (Baltimore). 2019;98(20):e15425. doi:
10.1097/MD.0000000000015425PubMedGoogle Scholar 13.Goesling
J, Henry
MJ, Moser
SE,
et al. Symptoms of depression are associated with opioid use regardless of pain severity and physical functioning among treatment-seeking patients with chronic pain.
J Pain. 2015;16(9):844-851. doi:
10.1016/j.jpain.2015.05.010PubMedGoogle ScholarCrossref 20.Wilson
L, Bekeris
J, Fiasconaro
M,
et al. Risk factors for new-onset depression or anxiety following total joint arthroplasty: the role of chronic opioid use.
Reg Anesth Pain Med. 2019;rapm-2019-100785.
PubMedGoogle Scholar 29.Wray
NR, Ripke
S, Mattheisen
M,
et al; eQTLGen; 23andMe; Major Depressive Disorder Working Group of the Psychiatric Genomics Consortium. Genome-wide association analyses identify 44 risk variants and refine the genetic architecture of major depression.
Nat Genet. 2018;50(5):668-681. doi:
10.1038/s41588-018-0090-3PubMedGoogle ScholarCrossref 33.Sanderson
E, Davey Smith
G, Windmeijer
F, Bowden
J. An examination of multivariable Mendelian randomization in the single-sample and two-sample summary data settings.
Int J Epidemiol. 2019;48(3):713-727. doi:
10.1093/ije/dyy262PubMedGoogle ScholarCrossref 35.Millard
LAC, Davies
NM, Gaunt
TR, Davey Smith
G, Tilling
K. Software application profile: PHESANT: a tool for performing automated phenome scans in UK Biobank.
Int J Epidemiol. 2018;47(1):29-35. doi:
10.1093/ije/dyx204PubMedGoogle ScholarCrossref 41.Bowden
J, Davey Smith
G, Haycock
PC, Burgess
S. Consistent estimation in mendelian randomization with some invalid instruments using a weighted median estimator.
Genet Epidemiol. 2016;40(4):304-314. doi:
10.1002/gepi.21965PubMedGoogle ScholarCrossref 42.Bowden
J, Del Greco M
F, Minelli
C, Davey Smith
G, Sheehan
N, Thompson
J. A framework for the investigation of pleiotropy in two-sample summary data Mendelian randomization.
Stat Med. 2017;36(11):1783-1802. doi:
10.1002/sim.7221PubMedGoogle ScholarCrossref 43.Bowden
J, Del Greco M
F, Minelli
C,
et al. Improving the accuracy of two-sample summary-data Mendelian randomization: moving beyond the NOME assumption.
Int J Epidemiol. 2019;48(3):728-742. doi:
10.1093/ije/dyy258PubMedGoogle ScholarCrossref 45.Rees
JMB, Wood
AM, Burgess
S. Extending the MR-Egger method for multivariable Mendelian randomization to correct for both measured and unmeasured pleiotropy.
Stat Med. 2017;36(29):4705-4718. doi:
10.1002/sim.7492PubMedGoogle ScholarCrossref 48.Zhou
H, Rentsch
CT, Cheng
Z,
et al; Veterans Affairs Million Veteran Program. Association of OPRM1 functional coding variant with opioid use disorder: a genome-wide association study.
JAMA Psychiatry. Published online June 3, 2020. doi:
10.1001/jamapsychiatry.2020.1206PubMedGoogle Scholar 58.Grattan
A, Sullivan
MD, Saunders
KW, Campbell
CI, Von Korff
MR. Depression and prescription opioid misuse among chronic opioid therapy recipients with no history of substance abuse.
Ann Fam Med. 2012;10(4):304-311. doi:
10.1370/afm.1371PubMedGoogle ScholarCrossref 59.Boscarino
JA, Hoffman
SN, Han
JJ. Opioid-use disorder among patients on long-term opioid therapy: impact of final DSM-5 diagnostic criteria on prevalence and correlates.
Subst Abuse Rehabil. 2015;6:83-91. doi:
10.2147/SAR.S85667PubMedGoogle ScholarCrossref 65.Bowden
J, Davey Smith
G, Haycock
PC, Burgess
S. Consistent estimation in Mendelian randomization with some invalid instruments using a weighted median estimator.
Genet Epidemiol. 2016;40(4):304-314. doi:
10.1002/gepi.21965PubMedGoogle ScholarCrossref 66.Lawlor
DA, Tilling
K, Davey Smith
G. Triangulation in aetiological epidemiology.
Int J Epidemiol. 2016;45(6):1866-1886.
PubMedGoogle Scholar 67.Gage
SH, Bowden
J, Davey Smith
G, Munafò
MR. Investigating causality in associations between education and smoking: a two-sample Mendelian randomization study.
Int J Epidemiol. 2018;47(4):1131-1140. doi:
10.1093/ije/dyy131PubMedGoogle ScholarCrossref 68.Choi
KW, Chen
CY, Stein
MB,
et al; Major Depressive Disorder Working Group of the Psychiatric Genomics Consortium. Assessment of bidirectional relationships between physical activity and depression among adults: a 2-sample mendelian randomization study.
JAMA Psychiatry. 2019;76(4):399-408. doi:
10.1001/jamapsychiatry.2018.4175PubMedGoogle ScholarCrossref 69.Fry
A, Littlejohns
TJ, Sudlow
C,
et al. Comparison of sociodemographic and health-related characteristics of UK Biobank participants with those of the general population.
Am J Epidemiol. 2017;186(9):1026-1034. doi:
10.1093/aje/kwx246PubMedGoogle ScholarCrossref