Genetic susceptibility of opioid receptor genes polymorphism to drug addiction: A candidate-gene association study
BMC Psychiatry volume 21, Article number: 5 (2021)
Like other complex diseases including drug addiction, genetic factors can interfere with the disease. In this study, three opioid genes (OPRM1, OPRD1, and OPRK1) were examined for an association with drug addiction among Jordanian males.
The study involved 498 addicts, in addition to 496 healthy controls and all from Arab descent.
The findings in this study showed that rs1799971 of the OPRM1 gene was in association with drug addiction for both alleles and genotypes with P-values = 0.002 and 0.01, respectively. In addition, a significant association between the dominant model (A/A vs G/A-G/G) of rs1799971 (OPRM1) and drug addiction (P-value = 0.003, OR = 1.59 (1.17–2.15)) was detected. Moreover, a genetic haplotype (AGGGCGACCCC) of theOPRM1 gene revealed a significant association with drug addiction (P-value = 0.01, OR = 1.56 (1.15–2.12)). We also found that the age of addicts, smoking, and marital status with genetic variants within OPRM1, OPRD1, and OPRK1 genes may be implicated in drug addiction risk.
We propose that rs1799971 of the OPRM1gene is a genetic risk factor for drug addiction among Jordanian males.
Drug addictions are chronic complex diseases that featured by repetitive and irresistible consumption of certain drugs that lead to temporary euphoria. Subsequently, the addict will encounter tolerance, craving, and relapse as a result of neurological changes . Up to 5% of the adults worldwide have been reported as drug abusers in (2010), drug addiction burden both society and individual’s life in many aspects including physical, mental, and social aspects. Moreover, over 200 thousand deaths annually are recording due to drug addiction . Addiction has been linked to environmental factors together with the genetic one . Inspecting the role of genetic variations in the etiology of the addiction may improve response to the treatments and help in disease prevention [4, 5]. Both genome-wide methods and candidate gene studies, looking for the underlying role of the genetic factor in drug addiction by identification of genes that are crucial to neuroadaptation .
Opioid receptors (OPR) are part of the opioid system and responsible for mediating drug reward and addiction. These receptors in the nervous system attached to substances and induce a series of chemical reactions that lead to pain relief . The OPRM1 gene located on chromosome number 6 and encode for endogenous G protein-coupled mu-opioid receptors [8, 9]. The OPRM1 is one of the most studied genes regarding drug addiction and several genetic variants within the gene were significantly implicated with drug addiction. (rs1799971) SNP of the OPRM1 gene has been extensively investigated in association with drug addiction in different ethnic groups including Caucasian  European American, African American , German , and Japanese .
In addition, the OPRD1 and OPRK1 genes encode for the δ-opioid receptor (DOR) and κ-opioid receptor (KOR) respectively are also belong to the opioid system. The OPRD1gene located on chromosome 1, while the OPRK1 gene located on chromosome 8 . Activation of the μ-opioid receptor and δ-opioid receptor mediate drug rewarding DOR has been assigned as mood regulator for anxiety and depression . Within the OPRD1 gene, both rs1042114 (G80T) and rs2234918 (T921C) have been in the spotlight for association with drug addiction. The latter SNP found to be correlated with heroin addiction in a German population . Furthermore, the two intronic SNPs (rs2236857 and rs581111) within the OPRD1 gene found to be linked with heroin addiction risk .
Areas κ-opioid receptor activation is associated with increases in both physical and psychological stress responses. Genetic polymorphisms of the OPRK1 gene have been previously studied in association with substances . In this study, we investigated the genetic susceptibility of three OPR genes (OPRM1, OPRD1 and OPRK1) to drug addiction among Jordanian males by screening and analyzing 29 candidate SNPs.
The study cohort consisted of 498 drug addicts selected depending on drug addiction criteria according to the Manual of Mental Disorders (DSM-IV) criteria . All participants were males from Jordanian and were hospitalized in 2018 for 8 months in the National Centre for Rehabilitation of Addicts (NCRA) of the ministry of health in Jordanian and the Drug Rehabilitation Centre of the Jordanian Public Security Directorate (DRC-PSD). In addition, 496 healthy Jordanian male subjects with no history of drug addiction or psychiatric disorders were chosen as controls. Clinical data and structured questionnaires were collected according to the Human Ethics Committee of the Jordanian Ministry of Health (MOH/REC/180057), and the Institutional Review Board (IRB) at Jordan University of Science and Technology (43/114/2018). Patients who meet the inclusion criteria and had none of the exclusion criteria were interviewed, had the study objectives explained to them, and were asked to sign a consent form to participate in the study. Required written informed consent has been also obtained from all participants. All methods of the current study were performed in accordance with the IRB guidelines and regulations. Finally, all data was coded and no specific individual was identified.
In this study, the 29 SNPs of three OPR genes (OPRM1, OPRD1, and OPRK1) were selected based on their clinical significance and according to early publications that reveled the plausible implication of these variants with drug addiction among other populations which make them polymorphisms of interest to be investigated as most of them have robust functions in neuroadaptation [20, 21]. Moreover, these polymorphisms selected to declare their influence on the Jordanian population of Arab descent.
DNA extraction and SNP genotyping
Venous blood samples were used to purify genomic DNA according to Wizard® Genomic DNA Purification Kit (Promega Corporation, USA). Agarose gel electrophoresis and the Nano-Drop ND-1000 UV-Vis Spectrophotometer (BioDrop, UK) were used to select the quality and quantity of the extracted genetic material. Candidate SNPs within (OPRM1, OPRD1, and OPRK1) were screened using the Sequencing technique. Selected samples for genotyping were diluted using nuclease-free water with a final concentration of 20 ng/μl (50-500 μl) and shipped on wet ice to the Australian Genome Research Facility (AGRF) (Australia). At the AGRF, the samples were genotyped using the Agena Bioscience MassARRAY® on a Compact Spectrometer, iPLEX GOLD chemistry.
The P-Value of Hardy-Weinberg equilibrium (HWE), was used to detect if the selected variants fulfill the (HWE) equation. Genotypic and allelic frequencies, genetic association, multiple genetic models, and haplotyping analyses were performed using SNPStats software (InstitutCatalàd’Oncologia, 2006). In addition, multinomial logistic regression was done using the Statistical Package for the Social Sciences (SPSS), version 25.0 (SPSS, Inc., Chicago, IL). P-values less than 0.05 were considered to be statistically significant.
Multiple testing corrections
Clinical characteristic of participants
All participants in this study were Jordanian males of Arab descent. Four hundred ninety-eighr cases were identified as drug addicts with respect to different substances including; synthetic cannabinoids (47.5%), cannabinoids (19.6%), amphetamine (5.7%), alcohol (5.5%), benzodiazepines (4.6%), opiates (4.4%), cocaine (1.1%), and cannabis (0.4%). Among the cases 89% used only one substance. However, 11% of the addicts in this study were multi substance user. The routes for substances consumption were varying between substances. For instance, smoking was commonly used for cannabinoids use, while opiate was administrated by the cases using injection route. In addition, benzodiazepines, alcohol, and amphetamine were consumed orally. Cases were using inhalation for cocaine substance use.
The range of study cohort’s ages was from 18 to 40 years with mean age 29 ± 6.9 for controls and 28.7 ± 9.3 for patients. Age, smoking, work, and marital status for cases were considered also in this study; 88.2% of the patients were smokers, (70.4%) were single and (26.8%) were unemployed.
Table 1 summarizes the list of selected genetic variants of three OPR genes (OPRM1, OPRD1, and OPRK1). The table also shows the chromosomal positions of SNPs, minor alleles and their frequencies, and the P-value of HWE for cases and controls. Among the investigated SNPs, rs495491 SNP of the OPRM1 gene did not fulfill the HWE eq. (P-value < 0.05) and were excluded from this study.
Genetic association analyses
The correlation between the studied SNPs and drug addiction was conducted by performing different statistical tests. Table 2 illustrates the differences in the allelic and genotypic distribution between patients and controls. Our findings exposed that only rs1799971 of the OPRM1 gene was in association with drug addiction in this study (P-values were; 0.002, 0.01) for alleles and genotypes difference respectively. A significant variation in the variant allele (G) frequency distribution was observed between cases and controls; 14% of patients carried the (G) allele compared to only1% of healthy participants with the same variant (Table 2). In addition, the genotype (GG) among patients was twice it among control, this finding may propose the (G) allele of rs1799971 (OPRM1) as a risk genetic locus for drug addiction. Moreover, different genetic models analysis was done for further confirmation. Table S1 displays the included models in this study (dominant and recessive). According to Bonferroni correction, the P-value was set to 0.02. The result show a significant association between the dominant model (A/A vs G/A-G/G) of rs1799971 (OPRM1) and drug addiction (P-value = 0.003, OR = 1.59 (1.17–2.15)). Haplotyping was also investigated as a part of the genetic association’s analyses in this study. As Table S2 demonstrates, one block of the OPRM1 gene (AGGGCGACCCC) exhibits significant association with drug addiction (P-value = 0.01, OR = 1.56 (1.15–2.12)).
In this study, regression analysis was done to detect the relationship between four major features that were significant for drug addiction among cases and the selected SNPs. Tables 3, 4, and 5 depict the outcome of regression analysis. However, age of the addicts, smoking and marital status revealed a relationship with different genotypes of some candidate SNPs of the investigated genes was only significant for the risk of drug addiction among males Jordanians.
Drug addiction is a chronic disease that involves both the environmental and genetic factors . Despite the vague knowledge about the robust role of the genetic parameters in drug addiction risk, several approaches were extensively investigated to put the baseline of genetic predisposition to drug addiction such as family, twin, and adoption studies [25,26,27]. This study focus on investigating the impact of many Single nucleotide variants within critical genes on drug addiction among Jordanians. Several previous studies conducted on Jordanian Arabs have investigated the involvement of candidate genes in drug addiction including DRD4, 5-HTTLPR  SLC6A4  and OPRM1 [3, 30, 31].
Our findings revealed that the OPRM1 gene may be implicated in drug addiction among Jordanians, in particular the functioning SNP (rs1799971) found to be significantly associated with drug addiction which is in consistence with a Chinese study on336 Han heroin addicts . In this study, the distribution of the variant genotype GG of rs1799971 within cases was twice it within controls which propose that the G allele of (rs1799971) within OPRM1 gene may influence drug addiction increased risk. The aforementioned SNP is a missense variant located on exon 1 and encodes for a change from asparagine to aspartate amino acids. It has been suggested that this change results in removing an N-glycosylation site in the extracellular domain which alters the activity of endogenous opioid binding receptors . These actions end in pain harbor and reduction in the response to analgesic drugs [9, 32] furthermore, it tend to increase administration of opioids . In contrast to this study finding [11, 13] reported no evidence of an association between rs1799971 polymorphism and opioid addiction among African-Americans/ European-Americans (EA) and Japanese population.
rs648893 (C/T) polymorphism an Intronic SNP of the OPRM1 gene that implicated in DNA splicing has been recently nominated as candidate for substance dependence including drugs . In accordance to our study, [34, 35] found that the rs648893 of OPRM1 was not associated with drug addiction among Iranian population and European Americans respectively. After all, ethnic variations reflect the differences in the drug addiction etiology among populations. Another significant outcome of this recent study, we found that one block of the OPRM1 gene (AGGGCGACCCC) exhibits significant association with drug addiction among Jordanian males.
In this study, none of the investigated SNPs of both OPRD1 and OPRK1 genes were in correlation with drug addiction among Jordanians. However, in a review of several studies about genetic variants of addiction susceptibility summarizes that both rs1051660 and rs1051660 of OPRK1 are involved in drug addiction , while rs1042114 polymorphism of OPRD1 was significantly correlated with substance dependence disorder among EA .
On the other hand, our regression analysis declared that the relationship between the age of addicts, smoking and marital status with genetic variants within OPRM1, OPRD1 and OPRK1 genes may be implicated in drug addiction risk. Remarkably, it has been reported that suicide attempts can be induced by drug addiction disorder. In regard, no correlation between significant OPRM1 SNPs (rs1799971, rs609148 and rs648893) and suicide attempts behavior among 426 European-Americans addicts was found . However, the real relation between genetic factors and clinical features of drug addiction is not fully elucidated, more studies and analyses are needed to understand the biological etiology of the disorder.
In conclusion, we propose that rs1799971 of the OPRM1 gene is a genetic risk factor for drug addiction among Jordanian males. To avoid genotyping bias in this study, we excluded the genetic association outcome caused by population stratification. This study was the first to investigate a genetic association of opioid gene variations with drug addiction in males of Arab descent. These findings may provide crucial knowledge to understand the drug addiction mechanisms in Middle Eastern population of Arab descent. Overall, the current study has examined three genes of opioid receptors that appear to be involved in the vulnerability to and the treatment of drug addiction in a Middle Eastern population of Arab descent. Finally, the results of this study provided additional clinical, epidemiological and genetic knowledge that may be useful in the context of further genetic and pharmacogenetic analyses to reduce the severity of drug consumption and improve drug abstinence.
Availability of data and materials
The complete processed SNP genotypic data for the three OPR genes is available as a supplementary file.
Manual of Mental Disorders
Drug Rehabilitation Centre of the Jordanian Public Security Directorate
National Centre for Rehabilitation of Addicts
Ahmed M, UlHaq I, Faisal M, Waseem D, Taqi MM. Implication of OPRM1 A118G polymorphism in opioids addicts in Pakistan: in vitro and Insilico analysis. J Mol Neurosci. 2018;65:472.
Jadidi N, Nakhaee N. Etiology of drug abuse: a narrative analysis. J Addict. 2014;2014:352835.
Al-Eitan LN, Jaradat SA, Hulse GK, Tay GK. Custom genotyping for substance addiction susceptibility genes in Jordanians of Arab descent. BMC Rese Notes. 2012;5:497.
AL-Eitan L, Tarkhan A. Practical challenges and translational issues in pharmacogenomics and personalized medicine from 2010 onwards. Curr Pharmacogenom Personalized Med. 2016;14:7–17.
Al-Eitan LN, Jaradt SA, Hulse GK, Tay GK. Pharmacogenetic approach to treating drug dependence: serotonin transporter gene (SLC6A4) promoter polymorphisms as treatment predictors in Jordanian Arabs. Curr Pharmacogenomics Pers Med. 2012;10:293–305.
Bevilacqua L, Goldman D. Genes and addictions. Clin Pharmacol Ther. 2009;85:359–61.
Wang SC, Chen YC, Lee CH, Cheng CM. Opioid addiction, genetic susceptibility, and medical treatments: A review. Int J Mol Sci. 2019;20:4294.
Ninković J, Roy S. Role of the mu-opioid receptor in opioid modulation of immune function. Amino Acids. 2013;45:9–24.
Randesi M, Rotrosen J, Nunes EV, Lee JD, Novo P, Levran O, et al. Variants of opioid genes and response to treatment of opioid use disorder with buprenorphine-naloxone versus extended-release naltrexone in Caucasians. Am J Drug Alcohol Abuse. 2020:1–8.
Colloca L, Wang Y, Martinez PE, Chang Y-PC, Ryan KA, Hodgkinson C, et al. OPRM1 rs1799971, COMT rs4680, and FAAH rs324420 genes interact with placebo procedures to induce hypoalgesia. Pain. 2019;160:1824–34.
Crowley JJ, Oslin DW, Patkar AA, Gottheil E, DeMariaJr PA, O'Brien CP, et al. A genetic association study of the mu opioid receptor and severe opioid dependence. Psychiatr Genet. 2003;13:169–73.
Mistry CJ, Bawor M, Desai D, Marsh DC, Samaan Z. Genetics of opioid dependence: A review of the genetic contribution to opioid dependence. Curr Psychiatry Rev. 2014;10:156–67.
Gelernter J, Kranzler H, Cubells J. Genetics of two μ opioid receptor gene (OPRM1) exon 1 polymorphisms: population studies, and allele frequencies in alcohol-and drug-dependent subjects. Mol Psychiatry. 1999;4:476–83.
Yuferov V, Fussell D, LaForge KS, Nielsen DA, Gordon D, Ho A, et al. Redefinition of the human kappa opioid receptor gene (OPRK1) structure and association of haplotypes with opiate addiction. Pharmacogenetics. 2004;14:793–804.
Crist RC, Ambrose-Lanci LM, Vaswani M, Clarke TK, Zeng A, Yuan C, et al. Case-control association analysis of polymorphisms in the δ-opioid receptor, OPRD1, with cocaine and opioid addicted populations. Drug Alcohol Depend. 2013;127:122–8.
Mayer P, Rochlitz H, Rauch E, Rommelspacher H, Hasse HE, Schmidt S, et al. Association between a delta opioid receptor gene polymorphism and heroin dependence in man. Neuroreport. 1997;8:2547–50.
Nelson EC, Lynskey MT, Heat AC, Wray N, Agrawal A, Shand FL, et al. Association of OPRD1 polymorphisms with heroin dependence in a large case-control series. Addict Biol. 2014;19:111–21.
Xu K, Seo D, Hodgkinson C, Hu Y, Goldman D, Sinha R. A variant on the kappa opioid receptor gene (OPRK1) is associated with stress response and related drug craving, limbic brain activation and cocaine relapse risk. Transl Psychiatry. 2013;3:e292.
Corley RP, Zeiger JS, Crowley T, Ehringer MA, Hewitt JK, Hopfer CJ, et al. Association of candidate genes with antisocial drug dependence in adolescents. Drug Alcohol Depend. 2008;96:90–8.
Zhang H, Kranzler HR, Yang BZ, Luo X, Gelernter J. The OPRD1 and OPRK1 loci in alcohol or drug dependence: OPRD1 variation modulates substance dependence risk. Mol Psychiatry. 2008;13:531–43.
American Psychiatric Association, DSM-5 Task Force. Diagnostic and statistical manual of mental disorders: DSM-5™ (5th ed.). American Psychiatric Publishing, Inc; 2013.
Nyholt DR. A simple correction for multiple testing for single-nucleotide polymorphisms in linkage disequilibrium with each other. Am J Hum Genet. 2004;74:765–9.
Li J, Ji L. Adjusting multiple testing in multilocus analyses using the eigenvalues of a correlation matrix. Heredity. 2005;95:221–7.
Agrawal A, Lynskey MT. Are there genetic influences on addiction: evidence from family, adoption and twin studies. Addiction. 2008;103:1069–81.
Wang JC, Kapoor M, Goate AM. The genetics of substance dependence. Annu Rev Genomics Hum Genet. 2012;13:241–61.
Al-Eitan L, Haddad Y. Emergence of pharmacogenomics in academic medicine and public health in Jordan: history, present state and prospects. Curr Pharmacogenomics Pers Med. 2014;12:167–75.
Alblooshi H, Hulse G, Osman W, El Kashefm A, Shawky M, Ghaferi H, et al. The frequency of DRD2 rs1076560 and OPRM1 rs1799971 in substance use disorder patients from the United Arab Emirates. Ann General Psychiatry. 2018;17:22.
Al-Eitan LN, Alshudaifat KM, Anani JY. Association of the DRD4 exon III and 5-HTTLPR VNTR polymorphisms with substance abuse in Jordanian Arab population. Gene. 2020;733:144267.
Al-Eitan LN, Jaradat SA, Qin W, Wildenauer DMB, Wildenauer DD, Hulse GK, et al. Characterization of serotonin transporter gene (SLC6A4) polymorphisms and its association with drug dependence in a Jordanian Arab population. Toxicol Ind Health. 2014;30:598–610.
Al-Eitan LN, Jaradat SA, Su SY, Tay GK, Hulse GK. Mu opioid receptor (OPRM1) as a predictor of treatment outcome in opiate-dependent individuals of Arab descent. Pharmacogenomics Pers Med. 2012;5:99–111.
Chou WY, Wang CH, Liu PH, Liu CC, Tseng CC, Jawan B. Human opioid receptor A118G polymorphism affects intravenous patient-controlled analgesia morphine consumption after total abdominal hysterectomy. Anesthesiology. 2006;105:334–7.
Oertel BG, Kettner M, Scholich K, Renné C, Roskam B, Geisslinger G, et al. A common human μ-opioid receptor genetic variant diminishes the receptor signaling efficacy in brain regions processing the sensory information of pain. J BiolChem. 2009;284:6530–5.
Zhang D, Shao C, Shao M, Yan P, Wang Y, Liu Y, et al. Effect of mu-opioid receptor gene polymorphisms on heroin-induced subjective responses in a Chinese population. Biol Psychiatry. 2007;61:1244–51.
Rezaeifar A, Dahmardeh F. The effect of OPRM1 rs648893 gene polymorphism on opioid addiction in an Iranian population in Zabol: A case-control study. Int J Basic Sci Med. 2019;4:143–7.
Xuei X, Flury-Wetherill L, Bierut L, Dick D, Jr JN, Foroud T, et al. The opioid system in alcohol and drug dependence: family-based association study. Am J Med Genet B Neuropsychiatr Genet. 2007;144:877–.
Arias AJ, Chan G, Gelernter J, Farrer L, Kranzler HR. Variation in OPRM1 and risk of suicidal behavior in drug-dependent individuals. Am J Addict. 2012;21:5–10.
We gratefully acknowledge the contribution of participating patients whose cooperation made this study possible. The authors also would like to express their gratitude to Jordan University and Science and Technology (JUST, Irbid, Jordan) for providing administrative and technical support.
The Scientific Research Support Fund (SRSF) at Jordanian Ministry of Higher Education supported this work under grant number MPH/1/43/2017.
Ethics approval and consent to participate
This study was approved by the Human Ethics Committee of the Jordanian Ministry of Health and the Institutional Review Board (IRB) at Jordan University of Science and Technology. Written informed consent has been obtained from all participants. All methods of the current study were performed in accordance with the IRB guidelines and regulations.
Consent for publication
The authors declare no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
AL-Eitan, L.N., Rababa’h, D.M. & Alghamdi, M.A. Genetic susceptibility of opioid receptor genes polymorphism to drug addiction: A candidate-gene association study. BMC Psychiatry 21, 5 (2021). https://doi.org/10.1186/s12888-020-03006-z