Genome-wide analyses reveal novel opioid use disorder loci and genetic overlap with schizophrenia, bipolar disorder, and major depression
=========================================================================================================================================

* Børge Holen
* Alexey A. Shadrin
* Romain Icick
* Guy Hindley
* Linn Rødevand
* Kevin S. O’Connell
* Oleksandr Frei
* Shahram Bahrami
* Weiqiu Cheng
* Nadine Parker
* Markos Tesfaye
* Piotr Jahołkowski
* Naz Karadag
* Anders M. Dale
* Srdjan Djurovic
* Olav B. Smeland
* Ole A. Andreassen

## Abstract

Opioid use disorder (OUD) and mental disorders are often comorbid, with increased morbidity and mortality. The causes underlying this relationship are poorly understood. Although these conditions are highly heritable, their shared genetic vulnerabilities remain unaccounted for. We applied the conditional/conjunctional false discovery rate (cond/conjFDR) approach to analyse summary statistics from independent genome wide association studies of OUD, SCZ, BD and MD. Next, we characterized the identified shared loci using biological annotation resources. OUD data was obtained from the Million Veteran Program (15,756 cases 99,039 controls). SCZ (53,386 cases 77,258 controls), BD (41,917 cases 371,549 controls) and MD (170,756 cases 329,443 controls) data was provided by the Psychiatric Genomics Consortium. We discovered genetic enrichment for OUD conditional on associations with SCZ, BD, MD and vice versa, indicating polygenic overlap with identification of 14 novel OUD loci at condFDR<0.05 and 7 unique loci shared between OUD and SCZ (n=2), BD (n=2) and MD (n=7) at conjFDR<0.05 with concordant effect directions, in line with estimated positive genetic correlations. Two loci were novel for OUD, one for BD and one for MD. Three OUD risk loci were shared with more than one psychiatric disorder, at *DRD2* on chromosome 11 (BD and MD), at *FURIN* on chromosome 15 (SCZ, BD and MD), and at the major histocompatibility complex region (SCZ and MD). Our findings provide new insights into the shared genetic architecture between OUD and SCZ, BD, and MD, indicating a complex genetic relationship, suggesting overlapping neurobiological pathways.

Keywords
*   Bipolar disorder
*   genetic overlap
*   genome-wide association study (GWAS)
*   major depression
*   opioid use disorder
*   schizophrenia

## 1 INTRODUCTION

Opioid use disorder (OUD) causes substantial morbidity and mortality worldwide 1. Some countries are more affected, such as the USA with the opioid crisis 2. Although OUD is less prevalent than other substance use disorders, OUD is an enormous public health burden with high rates of overdose deaths, which increased during the COVID-19 pandemic 3.

OUD is relatively prevalent in patients with severe mental disorders, such as schizophrenia (SCZ), bipolar disorder (BD) and major depression (MD) 4-6, all of which are associated with affective and psychotic symptoms to various degrees. Comorbidities lead to greater suffering compared to having either disorder alone 7, and often impede treatment 4. In patients with severe mental disorders the estimated prevalence of co-occurring OUD was 2.6-5.3% 4. Among adults with OUD the prevalence of a co-occurring severe mental disorder was estimated to be 27 % 8. At a sub-diagnostic threshold, an increased risk of nonmedical opioid use was reported for patients who had already been diagnosed with a severe mental disorder 9. Similarly, an increased risk of developing severe mental disorders for patients with existing nonmedical opioid use was also found 9. Further, there are several lines of evidence suggesting an overlapping neurobiological substrate related to the reward system and dopamine across mental illness in general 10, OUD 11, mood disorders 12, and psychotic disorders 13. E.g., can negative symptoms in SCZ, a category of internal heterogeneity, be conceptualized as a reward processing impairment 14. Better understanding of the mechanisms underlying these comorbidities is crucial for improving treatment and quality of life of affected individuals.

OUD and severe mental disorders can increase the risk of developing one another, have common environmental and genetic risk factors 15. Both OUD and severe mental disorders have moderate to high heritability estimates from twin and family studies, with OUD at 0.50 16 SCZ at 0.80 17, BD at 0.70 18 and MDD at 0.35 19. Recent progress in genotyping technology and international collaborations assembling large genome-wide association studies (GWAS) have provided novel insights into the genetic architecture of these complex disorders. A key discovery is that these disorders are highly polygenic 20, i.e., associated with many genetic risk variants, each with a small effect. Genetic studies have also suggested shared genetic aetiologies across severe mental disorders 15, including recent research which has found significant concordant genetic correlations (rg) between substance use disorders and severe mental disorders 16 with genetic correlation estimates between OUD and SCZ at 0.29, for BD at 0.16 and MD at 0.35 21, suggesting common genetic risk factors. However, genetic correlation is a genome-wide measure, which includes the effects of all Single Nucleotide Polymorphisms (SNPs), and thus does not identify overlap at the individual locus level. Moreover, genetic correlation is unable to capture genetic overlap in the presence of shared genetic variants with a mixture of same and opposite effect directions as they “cancel each other out” 22. In recent years, genome-wide analyses demonstrated genetic overlap with mixed effect directions among a wide range of human traits and disorders, regardless of their genetic correlations 23. Therefore, investigating genetic overlap beyond genetic correlation is required to further elucidate the genetic architecture of complex disorders and their genetic relationships. Here, we aimed to reveal more of the shared genetic architecture between OUD and the mental disorders SCZ, BD and MD by applying analytical tools designed for polygenic architectures 24. We applied the conditional/conjunctional false discovery rate (cond/conjFDR) tools to boost discovery of genetic variants and to identify overlapping genetic loci between OUD, SCZ, BD and MD, beyond genetic correlation 24.

## 2 MATERIALS AND METHODS

### 2.1 GWAS Samples

Independent GWAS data were obtained in the form of summary statistics (p-values and effect sizes) (Table 1.). For OUD, data from the Million Veteran Program (MVP), Yale-Penn and the Study of Addiction: Genetics and Environment (SAGE) 21 were used. In the MVP cohort, case and control status were defined using electronic health record–data. There were 8,529 affected European American (EA) individuals and 71,200 opioid-exposed EA controls, 4,032 affected African American (AA) individuals and 26,029 opioid-exposed AA controls. The Yale-Penn and SAGE datasets included 2,015 EA cases and 963 controls and 1,180 AA cases and 847 controls from a previous GWAS 25. To ensure compatibility of linkage disequilibrium (LD) pattern, our main analysis on OUD focused on individuals of EA ancestry and included a total number of 10,544 cases and 72,163 opioid-exposed controls. GWAS data on SCZ (53,386 cases and 77,258 controls) 26, BD (41,917 cases and 371,549 controls) 27 and MD (170,756 cases and 329,443 controls) 28 were obtained from the Psychiatric Genetics Consortium (PGC) and were of European ancestry. See Supplementary methods for details.

View this table:
[Table 1.](http://medrxiv.org/content/early/2022/09/09/2022.09.09.22279755/T1)

Table 1. 
Genome-wide association studies used in the present analyses. OUD = opioid use disorder; SCZ = schizophrenia; BD = bipolar disorder; MD = major depression; MVP = Million Veteran Program; PGC = Psychiatric Genomics Consortium.

### 2.2 Statistical Analyses

We visualized cross-trait enrichment using conditional quantile-quantile (Q-Q) plots, where distribution of p-values of all SNPs of a primary phenotype and for strata defined by the p-values for association with a secondary phenotype are represented. Cross-trait enrichment is evident if enrichment of statistical associations in the primary phenotype increases with increased significance of association with the secondary phenotype 29. We then applied the conditional FDR (condFDR), which leverages cross-trait enrichment between two phenotypes to improve genetic discovery. CondFDR readjusts the test-statistics in the primary phenotype by conditioning on SNP associations with the secondary phenotype, returning a condFDR value. OUD loci were identified at condFDR<0.05. We then applied the conjunctional FDR (conjFDR) tool to increase discovery of shared genetic loci jointly associated with OUD and in turn SCZ, BD and MD. The conjFDR approach is an extension of condFDR 29. Reversing the order of the phenotypes gives the condFDR value for the second phenotype conditioned on the first phenotype. ConjFDR then identifies SNPs which are significantly associated with both phenotypes 24. Shared loci were found at conjFDR<0.05, in line with prior literature 30. To control for spurious enrichment, random pruning was averaged over 500 iterations, and one SNP in each LD block (r2>0.1) was randomly selected for each iteration. We excluded SNPs within the major histocompatibility complex (MHC) (chr6:25000000-33000000), the chromosomal region 8p23.1 (location 7200000-12500000) and the gene *MAPT* (chr17:40000000-47000000) before fitting the FDR model given their complex LD structures that may bias FDR estimation 31.

We performed conjFDR analyses in independent samples for OUD (5,212 cases and 26,876 controls) 21, SCZ (22,778 cases and 35,362 controls) 32, BD (4,501 cases and 192,220 controls) [[https://r5.finngen.fi/](https://r5.finngen.fi/)], and MD (75,607 cases and 231,747 controls) 33 (Table 1). See Supplementary methods for details. We then conducted a version of the binomial test, i.e., sign concordance test 34. Finally, we tested for variants nominally significant at p<0.05 in each independent sample. All analyses were conducted in Oslo, Norway.

### 2.3 Genetic locus definition

We defined independent genetic loci according to the FUMA protocol 35. Briefly, independent significant genetic variants were identified as variants with conjFDR<0.05 and LD r2<0.6 with each other. A subset of these independent significant variants with LD r2<0.1 were then selected as lead variants. For each lead variant all candidate variants were identified as variants with LD r2≥0.6 with the lead variant. For a given lead variant the borders of the genetic locus were defined as min/max positional coordinates over all corresponding candidate variants. Loci were then merged if they were separated by less than 250kb. LD information was calculated from the 1000 Genomes Project European-ancestry reference panel 36. Directional effects of loci were evaluated by comparing z-scores of the lead SNPs. We investigated all discovered loci for overlap with previously identified loci using our internal database which includes updated GWAS on all the traits investigated in this study.

### 2.4 Functional annotation of loci shared between OUD and severe mental disorders

The FUMA application SNP2GENE was run to identify genes mapped to the candidate SNPs with a conjFDR value <0.1 35. All analyses were corrected for multiple comparisons. To minimize false positive and negative findings, we implemented a gene mapping strategy of qualification of at least two out of three of the following criteria: physical proximity of SNPs to genes, expression Quantitative Trait Loci (eQTL) and chromatin interaction mapping. SNP deleteriousness was measured by a Combined Annotation-Dependent Depletion (CADD)-score above 12.37 37. For complementarity we also used the open-source Open Targets Genetics application, Variant to Gene (V2G) 38 to map lead SNPs to genes. V2G utilizes physical proximity of SNPs to genes, molecular phenotype quantitative trait loci investigations (QTL) which maps SNPs to genes where expression level is influenced by allelic variation at the SNP level (eQTL), and protein Quantitative Trait Loci (pQTL) which maps SNPs to genes where protein level is influenced by allelic variation at the SNP level, and chromatin interaction where 3D DNA-DNA interactions are considered. V2G leverages this information in machine learning algorithms on the input of lead SNPs.

## 3 RESULTS

### 3.1 Cross-trait enrichment

To visualize cross-trait enrichment, we constructed stratified quantile-quantile (Q-Q) plots that show successive increments of SNP enrichment for OUD conditional on increasing levels of SNP associations with SCZ, BD and MD (Figure 1 A-C). We then reversed the stratified Q-Q plots displaying the SNP associations for SCZ, BD and MD conditional on SNP associations with OUD (Figure 1 D-F). In all cases we observed successive upward and leftward deflection of the Q-Q plots after conditioning the primary phenotype on the conditional phenotype, indicative of genetic overlap between the phenotypes.

![Figure 1.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/09/09/2022.09.09.22279755/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2022/09/09/2022.09.09.22279755/F1)

Figure 1. 
a-c: Conditional Q-Q plots of observed versus expected opioid use disorder (OUD) −log10 p-values (corrected for inflation) below the standard GWAS threshold of p < 5×10−8 as a function of the significance of the association with schizophrenia (SCZ), bipolar disorder (BD) and major depression (MD), at the level of p ≤ 0.1, p ≤ 0.01, p ≤ 0.001, respectively. d-f: Conditional Q-Q plots of observed versus expected SCZ, BD & MD −log10 p-values (corrected for inflation) below the standard GWAS threshold of p < 5×10−8 as a function of significance of association with OUD, at the level of p ≤ 0.1, p ≤ 0.01, p ≤ 0.001, respectively. The blue lines illustrate the standard enrichment for all SNPs irrespective of their association p-value in the second phenotype. The dashed line shows the null hypothesis. Successive leftward deflection for declining nominal p-values for a given disorder from the dashed line of no association, indicate that the proportion of non-null SNPs increase with higher levels of association with the other disorders.

### 3.2 CondFDR OUD loci

We applied condFDR to boost discovery of OUD loci conditional on mental disorders, and identified a total of 20 loci for OUD, including 14 novel (condFDR<0.05) (Table 2). This includes 9 loci for OUD conditioned on SCZ (6 novel), 8 loci conditioned on BD (5 novel) and 13 loci conditioned on MD (8 novel), including the locus involving the SNP rs1799971 (*OPRM1*) which is reported in the original OUD GWAS 21. The condFDR results are presented in Supplementary Tables 1-3.

View this table:
[Table 2.](http://medrxiv.org/content/early/2022/09/09/2022.09.09.22279755/T2)

Table 2. 
The most strongly associated single nucleotide polymorphisms (SNPs) in genomic loci associated with opioid use disorder (OUD) at conditional false discovery rate (condFDR) < 0.05 given association with schizophrenia (SCZ), bipolar disorder (BD) or major depression (MD) after merging regions < 250 kb apart into a single locus are shown. The table presents chromosomal position (Chr.), conditional false discovery rate (FDR) value, conditioned trait and OUD novelty status. For more details including genes mapped to these loci and the full list of all loci associated with OUD at condFDR < 0.05, see Supplemental Tables 1-3.

### 3.3 Shared loci between OUD and mental disorders (conjFDR)

To identify shared loci between OUD and the psychiatric disorders we performed conjFDR analyses with a threshold at <0.05 (Figure 2). We identified in total seven loci, two of which are novel for OUD and shared with MD at *RN7SKP157* and *B3GALTL* (Table 3). One locus was novel for BD at *DRD2* and one was novel for MD at *PPPC6* (Table 3). For all lead SNPs in the shared loci, the effect directions in OUD and each psychiatric disorder were concordant, i.e., risk of OUD was linked to higher risk of mental illness (Supplementary Tables 4-6). As expected, we found a locus shared between OUD, SCZ and MD located to the MHC region which has previously been implicated in these disorders at the genome wide significant level 39-41. This overlapping genetic signal suggests involvement of the immune system in the shared genetic risk of these disorders. However, the extended MHC region has highly complex LD structure spanning a large number of genes. We cannot reliably infer whether this overlapping signal reflects separate or shared loci, or causal genes.

![Figure 2.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/09/09/2022.09.09.22279755/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2022/09/09/2022.09.09.22279755/F2)

Figure 2. 
ConjFDR Manhattan plot. Common genetic variants jointly associated between opioid use disorder (OUD) and schizophrenia (SCZ), bipolar disorder (BD) and major depression (MD) at conjunctional false discovery rate (conjFDR) < 0.05. Showing the –log10 transformed conjFDR values for each single nucleotide polymorphism (SNP) on the y-axis and chromosomal positions along the x-axis. The dotted horizontal line represents the threshold for significant shared associations (conjFDR=0.05 i.e. −log10[conjFDR]□>1.3). Independent lead SNPs are encircled in black.

View this table:
[Table 3.](http://medrxiv.org/content/early/2022/09/09/2022.09.09.22279755/T3)

Table 3. 
Summary data of shared loci identified between opioid use disorder (OUD) and severe mental disorders (SMDs); schizophrenia (SCZ), bipolar disorder (BD) and major depression (MD) at conjunctional false discovery rate (conjFDR) < 0.05. Bp = base pair; SNP = single nucleotide polymorphism; Chr = chromosome, Bold p-value refers to the listed SMD z-score.

### 3.4 Functional characterization

Functional analyses of candidate SNPs in the loci conditioned on mental disorders with condFDR revealed that most were intergenic and intronic (Supplementary Tables 1-3). After testing for enrichment (e>1) according to functional category exonic (e=2.27, Fisher’s p=0.03) and intronic variants (e=1.15, Fisher’s p=0.03) were significantly enriched for OUD|SCZ. For OUD|BD intergenic (e=1.17, Fisher’s p=5.02E-07) and intronic (e=1.08, Fisher’s p=0.05) variants were significantly enriched. For OUD|MD intronic (e=1.50, Fisher’s p=2.42E-25) variant were significantly enriched. Functional analyses of candidate SNPs in the shared loci identified by conjFDR revealed that most were intergenic and intronic (Supplementary Tables 7-9). After testing for enrichment according to functional category UTR3 (e=5.43, Fisher’s p=0.01) and exonic variants (e=4.99, Fisher’s p=0.01) were significantly enriched for OUD/SCZ. For OUD/BD exonic (e=8.21, Fisher’s p= 9.66E-05) and intronic (e=2.21, Fisher’s p=2.92E-14) variants were significantly enriched. For OUD/MD intronic (e=1.77, Fisher’s p= 3.38E-35) variant were significantly enriched. Mapping of shared loci to genes 35 resulted in 98 identified genes shared between OUD and SCZ in two loci (Supplementary Table 10). Thirty of these were left after filtering on two out of three mapping criteria (positional, eQTL, and chromatin interaction information). We identified 19 genes shared between OUD and BD in two loci (Supplementary Table 11). Five of these passed filtering with two out of three mapping strategies. Finally, we identified 145 genes shared between OUD and MD in seven loci (Supplementary Table 12). Thirty-six genes passed the filtering strategy.

### 3.5 Validation analyses

We used the sign concordance test 34 to validate the shared loci which produced one of two concordant variants (50%; binomial p=0.75) for SCZ, two of two (100%; binomial p=0.25) concordant variants for BD, six of seven (86%; binomial p=0.06) concordant variants for MD, and four of eight (50%; binomial p=0.64) concordant variants for OUD. We further assessed validity of the lead SNPs in loci identified at conjFDR<0.05 by testing for presence in independent datasets at p<0.05. One lead variant (rs4702 at *FURIN*) was present in the SCZ GWAS of East Asian ancestry 32. No lead variants were present in the BD GWAS of European ancestry ([https://r5.finngen.fi/](https://r5.finngen.fi/)). For MD, three of seven lead SNPs (rs28524171 at *RN7SKP157*, rs6929812 within MHC, and rs1800498 at *DRD2*) was validated in the independent MD GWAS of European ancestry 33. Finally, no lead variants in the shared loci for OUD and psychiatric disorders, had a p-value <0.05 in the independent AA GWAS dataset on OUD 21. A total of four of the OUD risk loci identified by conjFDR were identified in recent larger GWAS on OUD, specifically the *PPP6C* 42, *NCAM1* 39, *DRD2* 39 and *FURIN* 39,42,43 loci, supporting the validity of the results.

## 4 DISCUSSION

In the current study, we used the cond/conjFDR analytical method and leveraged cross-trait enrichment between OUD, and SCZ, BD and MD to improve the statistical power for discovery of shared genetic loci to shed light on their genetic relationships. We identified polygenic overlap beyond genetic correlation. Specifically, we discovered 14 novel OUD loci and 7 novel loci jointly associated with OUD and psychiatric disorders (Figure 2, Table 2 and 3). There was a concordant allelic direction of effects in all the shared loci, indicating that these genetic variants increase the risk of both OUD and the respective psychiatric disorders. Several lines of evidence suggest a pleiotropic nature of mental traits and disorders 44. Our findings are in line with the observation that multiple genetic variants with small effect sizes influence many traits to different degrees 17. The currently identified shared loci explain only a small fraction of the genetic architecture of OUD and severe mental disorders and, thus, the liability to these disorders. We identified 14 novel loci for OUD, but there are still numerous undetected common variants, which will be identified with larger GWAS samples 17. In addition, environmental factors, copy number and rare genetic variants play an important role in the aetiology of OUD 45, SCZ 46, BD 47 and MD 48. The present results increase our understanding of the underlying genetic variants influencing OUD and the comorbidity between OUD and the severe psychiatric disorders SCZ, BD and MD. These common genetic variants may help generate novel hypotheses about the underlying molecular relationships and their role in the development of psychopathology, whether overlapping loci represent components of discrete clusters of risk (e.g., for substance use disorders) or a general vulnerability to psychiatric symptoms in general 49, is still an open question.

We identified a pleiotropic locus on chromosome 15 (lead SNP rs4702) which was shared between OUD, SCZ, BD and MD at conjFDR<0.05. It has previously been associated with SCZ 40, BD 27 and MD 41. This locus was not identified in the original OUD GWAS 21, but reached genome-wide significance in subsequent GWAS on OUD, validating this finding 39,42,43. The lead SNP rs4702 is shown to influence the expression of its nearest gene *FURIN* in neurons derived from human induced pluripotent stem cells, and influence synaptic function in synergy with other SCZ risk variants 50. *FURIN* encodes a cleaving and activating enzyme involved in cleavage of the endogenous opioid proenkephalin 51, which is found in most parts of the brain and nervous system 52. We discovered a novel OUD locus at chromosome 5 shared with MD, which has previously been associated with both SCZ 32 and MD 53. The nearest gene to the lead SNP is pseudogene *RN7SKP157*, while the most likely causal genes according to the V2G algorithm were *KIF2A* and *DIMT1*. We also discovered a novel OUD locus at chromosome 13, implicating three genes (*HSPH1, RXFP2* and *B3GALTL*). This locus is shared with MD and has previously been implicated in SCZ 32 and MD 28. FUMA and V2G both identified the nearest gene and most likely implicated gene as *B3GALTL*, respectively. We identified another OUD risk locus on chromosome 11 which was shared with both BD and MD. The locus is novel for BD and has previously been implicated in OUD 43, MD 41, as well as for alcohol use disorder 54 and SCZ 40.The nearest gene to the lead SNP is *DRD2*, while the most likely causal genes according to the V2G algorithm was *ANKK1* and *TTC12*. The gene cluster *TTC12-ANKK1-DRD2-NCAM1* (NTAD) have previously been widely linked to nicotine 55, alcohol 56 and other drug dependencies 57. *NCAM1* has previously been found to be associated with SCZ, BD, MD, alcohol use disorder and cannabis use disorder 54. Sequence variants in the *DRD2* gene has been associated with SCZ 58, Parkinson’s disease 59 and opioid addiction 60.

The identified genetic variants seem to support neurobiological hypotheses about the reward system in substance use disorders (SUD) and severe mental disorders (SMDs). The identified loci include variants that may result in changed function of the mesocorticolimbic circuit, influenced by altered metabolism of opioids (*FURIN* variants) and altered dopamine receptor function (*DRD2* variants). Such functional modification of the reward system may lead to an increased risk of SUDs and SMDs, for example via adapted dopamine transmission leading to reduced reward processing 14 which may underlie negative psychotic symptoms and depressive mood, and increased susceptibility for exogenous opioids due to elevated dopamine release 61 resulting in higher risk of opioid use disorder. Our findings warrant experimental validation to elucidate the pathophysiological mechanisms of the reward circuitry underlying comorbid SUDs and SMDs.

The current study has some limitations. We cannot exclude sporadic, non-systematic sample overlap; however, we calculated linkage disequilibrium score regression genetic covariance intercepts and neither is significant (all p-values>0.05), suggesting no significant effect of overlap between the samples. We included datasets from a diverse set of populations, but larger datasets for other cohorts than Europeans are needed to uncover more of the genetic underpinnings of OUD and SMDs.

To conclude, our study suggests overlapping genetic architecture between OUD and SCZ, BD, and MD, indicating similarities in the genetic architectures of these debilitating disorders. We highlight loci previously associated with SUDs and suggest underlying molecular pathways for OUD. Uncovering the underlying genetics of these disorders can lead to improvements in patient stratification and identification of novel targets for drug development.

## Supporting information

Supplemental information [[supplements/279755_file02.docx]](pending:yes)

Supplemental tables [[supplements/279755_file03.xlsx]](pending:yes)

## Data Availability

All data produced in the present study are available upon reasonable request to the authors

## DECLARATIONS OF INTERESTS

OAA has received speaker’s honorarium from Sunovion and Lundbeck and is a consultant for Healthlytix. AMD is a founder of and holds equity interest in CorTechs Labs and serves on its scientific advisory board. He is also a member of the Scientific Advisory Board of Healthlytix and receives research funding from General Electric Healthcare (GEHC). The terms of these arrangements have been reviewed and approved by the University of California, San Diego in accordance with its conflict-of-interest policies. Remaining authors have no conflicts of interest to declare.

## AUTHOR CONTRIBUTIONS

BH, AAS, OBS, and OAA were responsible for the study concept and design. AAS, KSOC and OAA contributed to the acquisition of summary data. BH, OBS, OAA, and AAS assisted with data analysis and interpretation of findings. AAS prepared figures. BH drafted the manuscript. OBS and OAA provided critical revision of the manuscript for important intellectual content. All authors critically reviewed content and approved final version for publication.

## ACKNOWLEDGMENTS

We thank the Psychiatric Genomics Consortium, the Million Veteran Program, FinnGen biobank and 23andMe for access to data, and the many people who provided DNA samples. The work was supported by the Research Council of Norway (324499, 324252, 273291, 223273), South-Eastern Norway Regional Health Authority (2019-108), European Union’s Horizon 2020 Research and Innovation Programme grant # 847776 (CoMorMent), and KG Jebsen Stiftelsen (SKGJ□MED□ 021). A.M.D. was supported by NIH grant U24DA041123. This work was performed on the TSD (Tjeneste for Sensitive Data) facilities, owned by the University of Oslo, operated and developed by the TSD service group at the University of Oslo, IT-Department (USIT). (tsd-drift{at}usit.uio.no).

*   Received September 9, 2022.
*   Revision received September 9, 2022.
*   Accepted September 9, 2022.


*   © 2022, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/)

## REFERENCES

1.  1.Strang J, Volkow ND, Degenhardt L, et al. Opioid use disorder. Nat Rev Dis Primers. 2020;6(1):3.
    
    

2.  2.Williams AR, Bisaga A. From AIDS to Opioids — How to Combat an Epidemic. N Engl J Med. 2016;375(9):813–815.
    
    

3.  3.The L. A time of crisis for the opioid epidemic in the USA. Lancet. 2021;398(10297):277.
    
    

4.  4.Toftdahl NG, Nordentoft M, Hjorthøj C. Prevalence of substance use disorders in psychiatric patients: a nationwide Danish population-based study. Soc Psychiatry Psychiatr Epidemiol. 2016;51(1):129–140.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00127-015-1104-4&link_type=DOI) 

5.  5.Chiappelli J, Chen S, Hackman A,  Elliot Hong L. Evidence for differential opioid use disorder in schizophrenia in an addiction treatment population. Schizophr Res. 2018;194:26–31.
    
    

6.  6.Martins SS, Fenton MC, Keyes KM, Blanco C, Zhu H, Storr CL. Mood and anxiety disorders and their association with non-medical prescription opioid use and prescription opioid-use disorder: longitudinal evidence from the National Epidemiologic Study on Alcohol and Related Conditions. Psychol Med. 2012;42(6):1261–1272.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0033291711002145&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21999943&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F09%2F2022.09.09.22279755.atom) 

7.  7.Drake RE, Brunette MF. Complications of severe mental illness related to alcohol and drug use disorders. Recent Dev Alcohol. 1998;14:285–299.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9751950&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F09%2F2022.09.09.22279755.atom) 

8.  8.Jones CM, McCance-Katz EF. Co-occurring substance use and mental disorders among adults with opioid use disorder. Drug Alcohol Depend. 2019;197:78–82.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F09%2F2022.09.09.22279755.atom) 

9.  9.Martins SS, Keyes KM, Storr CL, Zhu H, Chilcoat HD. Pathways between nonmedical opioid use/dependence and psychiatric disorders: Results from the National Epidemiologic Survey on Alcohol and Related Conditions. Drug Alcohol Depend. 2009;103(1):16–24.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.drugalcdep.2009.01.019&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19414225&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F09%2F2022.09.09.22279755.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000266760700003&link_type=ISI) 

10. 10.Pujara M, Koenigs M. Mechanisms of reward circuit dysfunction in psychiatric illness: prefrontal-striatal interactions. Neuroscientist. 2014;20(1):82–95.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/1073858413499407&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23924665&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F09%2F2022.09.09.22279755.atom) 

11. 11.Zhu Y, Wienecke CFR, Nachtrab G, Chen X. A thalamic input to the nucleus accumbens mediates opiate dependence. Nature. 2016;530(7589):219–222.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature16954&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26840481&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F09%2F2022.09.09.22279755.atom) 

12. 12.Russo SJ, Nestler EJ. The brain reward circuitry in mood disorders. Nat Rev Neurosci. 2013;14(9):609–625.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrn3381&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23942470&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F09%2F2022.09.09.22279755.atom) 

13. 13.Nielsen M, Rostrup E, Broberg BV, Wulff S, Glenthøj B. Negative Symptoms and Reward Disturbances in Schizophrenia Before and After Antipsychotic Monotherapy. Clin EEG Neurosci. 2018;49(1):36–45.
    
    

14. 14.Galderisi S, Mucci A, Buchanan RW, Arango C. Negative symptoms of schizophrenia: new developments and unanswered research questions. Lancet Psychiatry. 2018;5(8):664–677.
    
    

15. 15.Smeland OB, Frei O, Fan C-C, Shadrin A, Dale AM, Andreassen OA. The emerging pattern of shared polygenic architecture of psychiatric disorders, conceptual and methodological challenges. Psychiatr Genet. 2019;29(5):152–159.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/YPG.0000000000000234&link_type=DOI) 

16. 16.Deak JD, Johnson EC. Genetics of substance use disorders: a review. Psychol Med. 2021:1–12.
    
    

17. 17.Smeland OB, Frei O, Dale AM, Andreassen OA. The polygenic architecture of schizophrenia — rethinking pathogenesis and nosology. Nat Rev Neurol. 2020;16(7):366–379.
    
    

18. 18.Gordovez FJA, McMahon FJ. The genetics of bipolar disorder. Mol Psychiatry. 2020;25(3):544–559.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41380-019-0634-7&link_type=DOI) 

19. 19.Coleman JRI, Peyrot WJ, Purves KL, et al. Genome-wide gene-environment analyses of major depressive disorder and reported lifetime traumatic experiences in UK Biobank. Mol Psychiatry. 2020;25(7):1430–1446.
    
    

20. 20.Sullivan PF, Geschwind DH. Defining the Genetic, Genomic, Cellular, and Diagnostic Architectures of Psychiatric Disorders. Cell. 2019;177(1):162–183.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2019.01.015&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30901538&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F09%2F2022.09.09.22279755.atom) 

21. 21.Zhou H, Rentsch CT, Cheng Z, et al. Association of OPRM1 Functional Coding Variant With Opioid Use Disorder: A Genome-Wide Association Study. JAMA Psychiatry. 2020;77(10):1072–1080.
    
    

22. 22.Bulik-Sullivan B, Finucane HK, Anttila V, et al. An atlas of genetic correlations across human diseases and traits. Nat Genet. 2015;47(11):1236.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3406&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26414676&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F09%2F2022.09.09.22279755.atom) 

23. 23.Smeland OB, Bahrami S, Frei O, et al. Genome-wide analysis reveals extensive genetic overlap between schizophrenia, bipolar disorder, and intelligence. Mol Psychiatry. 2020;25(4):844–853.
    
    

24. 24.Andreassen OA, Djurovic S, Thompson WK, et al. Improved detection of common variants associated with schizophrenia by leveraging pleiotropy with cardiovascular-disease risk factors. Am J Hum Genet. 2013;92(2):197–209.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2013.01.001&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23375658&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F09%2F2022.09.09.22279755.atom) 

25. 25.Cheng Z, Zhou H, Sherva R, Farrer LA, Kranzler HR, Gelernter J. Genome-wide Association Study Identifies a Regulatory Variant of RGMA Associated With Opioid Dependence in European Americans. Biol Psychiatry. 2018;84(10):762–770.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2017.12.016&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F09%2F2022.09.09.22279755.atom) 

26. 26.Ripke S, Walters JT, O’Donovan MC. Mapping genomic loci prioritises genes and implicates synaptic biology in schizophrenia. medRxiv. 2020:2020.2009.2012.20192922.
    
    

27. 27.Mullins N, Forstner AJ, O’Connell KS, et al. Genome-wide association study of more than 40,000 bipolar disorder cases provides new insights into the underlying biology. Nat Genet. 2021;53(6):817–829.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-021-00857-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34002096&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F09%2F2022.09.09.22279755.atom) 

28. 28.Howard DM, Adams MJ, Shirali M, et al. Genome-wide association study of depression phenotypes in UK Biobank identifies variants in excitatory synaptic pathways. Nat Commun. 2018;9(1):1470.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-018-03819-3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F09%2F2022.09.09.22279755.atom) 

29. 29.Smeland OB, Frei O, Shadrin A, et al. Discovery of shared genomic loci using the conditional false discovery rate approach. Hum Genet. 2020;139(1):85–94.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00439-019-02060-2&link_type=DOI) 

30. 30.Smeland OB, Frei O, Kauppi K, et al. Identification of Genetic Loci Jointly Influencing Schizophrenia Risk and the Cognitive Traits of Verbal-Numerical Reasoning, Reaction Time, and General Cognitive Function. JAMA Psychiatry. 2017;74(10):1065–1075.
    
    

31. 31.Schwartzman A, Lin X. The effect of correlation in false discovery rate estimation. Biometrika. 2011;98(1):199–214.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/biomet/asq075&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23049127&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F09%2F2022.09.09.22279755.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000287759000015&link_type=ISI) 

32. 32.Lam M, Chen CY, Li Z, et al. Comparative genetic architectures of schizophrenia in East Asian and European populations. Nat Genet. 2019;51(12):1670–1678.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/S41588-019-0512-X&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F09%2F2022.09.09.22279755.atom) 

33. 33.Hyde CL, Nagle MW, Tian C, et al. Identification of 15 genetic loci associated with risk of major depression in individuals of European descent. Nat Genet. 2016;48(9):1031–1036.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3623&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27479909&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F09%2F2022.09.09.22279755.atom) 

34. 34.Lee JJ, Wedow R, Okbay A, et al. Gene discovery and polygenic prediction from a genome-wide association study of educational attainment in 1.1 million individuals. Nat Genet. 2018;50(8):1112–1121.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0147-3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30038396&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F09%2F2022.09.09.22279755.atom) 

35. 35.Watanabe K, Taskesen E, van Bochoven A, Posthuma D. Functional mapping and annotation of genetic associations with FUMA. Nat Commun. 2017;8(1):1826–1826.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-017-01261-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F09%2F2022.09.09.22279755.atom) 

36. 36.Auton A, Brooks LD, Durbin RM, et al. A global reference for human genetic variation. Nature. 2015;526(7571):68–74.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature15393&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26432245&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F09%2F2022.09.09.22279755.atom) 

37. 37.Kircher M, Witten DM, Jain P, O’Roak BJ, Cooper GM, Shendure J. A general framework for estimating the relative pathogenicity of human genetic variants. Nat Genet. 2014;46(3):310–315.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.2892&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24487276&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F09%2F2022.09.09.22279755.atom) 

38. 38.Ghoussaini M, Mountjoy E, Carmona M, et al. Open Targets Genetics: systematic identification of trait-associated genes using large-scale genetics and functional genomics. Nucleic Acids Res. 2020;49(D1):D1311–D1320.
    
    

39. 39.Kember RL, Vickers-Smith R, Xu H, et al. Cross-ancestry meta-analysis of opioid use disorder uncovers novel loci with predominant effects on brain. medRxiv. 2021:2021.2012.2013.21267480.
    
    

40. 40.Trubetskoy V, Pardiñas AF, Qi T, et al. Mapping genomic loci implicates genes and synaptic biology in schizophrenia. Nature. 2022.
    
    

41. 41.Bahrami S, Steen NE, Shadrin A, et al. Shared Genetic Loci Between Body Mass Index and Major Psychiatric Disorders: A Genome-wide Association Study. JAMA Psychiatry. 2020;77(5):503–512.
    
    

42. 42.Gaddis N, Mathur R, Marks J, et al. Multi-trait genome-wide association study of opioid addiction: OPRM1 and Beyond. medRxiv. 2021:2021.2009.2013.21263503.
    
    

43. 43.Deak JD, Zhou H, Galimberti M, et al. 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. Mol Psychiatry. 2022.
    
    

44. 44.Watanabe K, Stringer S, Frei O, et al. A global overview of pleiotropy and genetic architecture in complex traits. Nat Genet. 2019;51(9):1339–1348.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-019-0481-0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F09%2F2022.09.09.22279755.atom) 

45. 45.Li D, Zhao H, Kranzler HR, et al. Genome-Wide Association Study of Copy Number Variations (CNVs) with Opioid Dependence. Neuropsychopharmacology. 2015;40(4):1016–1026.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/npp.2014.290&link_type=DOI) 

46. 46.Singh T, Poterba T, Curtis D, et al. Rare coding variants in ten genes confer substantial risk for schizophrenia. Nature. 2022;604(7906):509–516.
    
    

47. 47.Sul JH, Service SK, Huang AY, et al. Contribution of common and rare variants to bipolar disorder susceptibility in extended pedigrees from population isolates. Transl Psychiatry. 2020;10(1):74.
    
    

48. 48.Cheng S, Cheng B, Liu L, et al. Exome-wide screening identifies novel rare risk variants for major depression disorder. Mol Psychiatry. 2022;27(7):3069–3074.
    
    

49. 49.Caspi A, Houts RM, Belsky DW, et al. The p Factor: One General Psychopathology Factor in the Structure of Psychiatric Disorders? Clin Psychol Sci. 2014;2(2):119–137.
    
    

50. 50.Schrode N, Ho S-M, Yamamuro K, et al. Synergistic effects of common schizophrenia risk variants. Nat Genet. 2019;51(10):1475–1485.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-019-0497-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31548722&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F09%2F2022.09.09.22279755.atom) 

51. 51.Breslin MB, Lindberg I, Benjannet S, Mathis JP, Lazure C, Seidah NG. Differential processing of proenkephalin by prohormone convertases 1(3) and 2 and furin. J Biol Chem. 1993;268(36):27084–27093.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamJjIjtzOjU6InJlc2lkIjtzOjEyOiIyNjgvMzYvMjcwODQiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wOS8wOS8yMDIyLjA5LjA5LjIyMjc5NzU1LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

52. 52.Schafer MKH, Day R, Watson SJ, Akil H. Distribution of Opioids in Brain and Peripheral Tissues. In: Almeida OFX, Shippenberg TS, eds. Neurobiology of Opioids. Berlin, Heidelberg: Springer Berlin Heidelberg; 1991:53–71.
    
    

53. 53.Howard DM, Adams MJ, Clarke T-K, et al. Genome-wide meta-analysis of depression identifies 102 independent variants and highlights the importance of the prefrontal brain regions. Nat Neurosci. 2019;22(3):343–352.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F09%2F2022.09.09.22279755.atom) 

54. 54.Gutiérrez-Sacristán A, Grosdidier S, Valverde O, et al. PsyGeNET: a knowledge platform on psychiatric disorders and their genes. Bioinformatics. 2015;31(18):3075–3077.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btv301&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25964630&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F09%2F2022.09.09.22279755.atom) 

55. 55.Gelernter J, Yu Y, Weiss R, et al. Haplotype spanning TTC12 and ANKK1, flanked by the DRD2 and NCAM1 loci, is strongly associated to nicotine dependence in two distinct American populations. Hum Mol Genet. 2006;15(24):3498–3507.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/ddl426&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17085484&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F09%2F2022.09.09.22279755.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000242716200002&link_type=ISI) 

56. 56.Yang B-Z, Kranzler HR, Zhao H, Gruen JR, Luo X, Gelernter J. Association of haplotypic variants in DRD2, ANKK1, TTC12 and NCAM1 to alcohol dependence in independent case– control and family samples. Hum Mol Genet. 2007;16(23):2844–2853.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/ddm240&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17761687&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F09%2F2022.09.09.22279755.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000251036400006&link_type=ISI) 

57. 57.Nelson EC, Lynskey MT, Heath AC, et al. ANKK1, TTC12, and NCAM1 polymorphisms and heroin dependence: importance of considering drug exposure. JAMA psychiatry. 2013;70(3):325–333.
    
    

58. 58.Wang J, Su P, Yang J, et al. The D2R-DISC1 protein complex and associated proteins are altered in schizophrenia and normalized with antipsychotic treatment. J Psychiatry Neurosci. 2022;47(2):E134–e147.
    
    

59. 59.Pelzer EA, Stürmer S, Feis D-L, et al. Clustering of Parkinson subtypes reveals strong influence of DRD2 polymorphism and gender. Sci Rep. 2022;12(1):6038.
    
    

60. 60.Clarke TK, Weiss AR, Ferarro TN, et al. The dopamine receptor D2 (DRD2) SNP rs1076560 is associated with opioid addiction. Ann Hum Genet. 2014;78(1):33–39.
    
    

61. 61.Berridge KC. The debate over dopamine’s role in reward: the case for incentive salience. Psychopharmacology. 2007;191(3):391–431.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00213-006-0578-x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17072591&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F09%2F2022.09.09.22279755.atom)