Genome-wide association study in individuals of European and African ancestry and multi-trait analysis of opioid use disorder identifies 19 independent genome-wide significant risk loci

Despite the large toll of opioid use disorder (OUD), genome-wide association studies (GWAS) of OUD to date have yielded few susceptibility loci. We performed a large-scale GWAS of OUD in individuals of European (EUR) and African (AFR) ancestry, optimizing genetic informativeness by performing MTAG (Multi-trait analysis of GWAS) with genetically correlated substance use disorders (SUDs). Meta-analysis included seven cohorts: the Million Veteran Program, Psychiatric Genomics Consortium, iPSYCH, FinnGen, Partners Biobank, BioVU, and Yale-Penn 3, resulting in a total N = 639,063 (N(cases) = 20,686;N(effective) = 77,026) across ancestries. OUD cases were defined as having a lifetime OUD diagnosis, and controls as anyone not known to meet OUD criteria. We estimated SNP-heritability (h(2)(SNP)) and genetic correlations (r(g)). Based on genetic correlation, we performed MTAG on OUD, alcohol use disorder (AUD), and cannabis use disorder (CanUD). A leave-one-out polygenic risk score (PRS) analysis was performed to compare OUD and OUD-MTAG PRS as predictors of OUD case status in Yale-Penn 3. The EUR meta-analysis identified three genome-wide significant (GWS; p ≤ 5 × 10(−)(8)) lead SNPs—one at FURIN (rs11372849; p = 9.54 × 10(−)(10)) and two OPRM1 variants (rs1799971, p = 4.92 × 10(−)(09); rs79704991, p = 1.11 × 10(−)(08); r(2) = 0.02). Rs1799971 (p = 4.91 × 10(−)(08)) and another OPRM1 variant (rs9478500; p = 1.95 × 10(−)(08); r(2) = 0.03) were identified in the cross-ancestry meta-analysis. Estimated h(2)(SNP) was 12.75%, with strong r(g) with CanUD (r(g) = 0.82; p = 1.14 × 10(−)(47)) and AUD (r(g) = 0.77; p = 6.36 × 10(−)(78)). The OUD-MTAG resulted in a GWAS N(equivalent) = 128,748 and 18 independent GWS loci, some mapping to genes or gene regions that have previously been associated with psychiatric or addiction phenotypes. The OUD-MTAG PRS accounted for 3.81% of OUD variance (beta = 0.61;s.e. = 0.066; p = 2.00 × 10(−)(16)) compared to 2.41% (beta = 0.45; s.e. = 0.058; p = 2.90 × 10(−)(13)) explained by the OUD PRS. The current study identified OUD variant associations at OPRM1, single variant associations with FURIN, and 18 GWS associations in the OUD-MTAG. The genetic architecture of OUD is likely influenced by both OUD-specific loci and loci shared across SUDs.