- Research article
- Open Access
- Open Peer Review
Association between prodynorphin gene polymorphisms and opioid dependence susceptibility: a meta-analysis
BMC Psychiatryvolume 19, Article number: 281 (2019)
Prodynorphin (PDYN) gene polymorphisms have been linked with opioid dependence (OD) with conflicting outcomes, the aim of this study is to synthesize the existing evidence of the association between PDYN polymorphisms and OD susceptibility.
Four databases including PubMed, EMBASE, Web of Science, and Wanfang were retrieved for relevant studies before August, 2018. All identified studies were evaluated using predetermined inclusion and exclusion criteria. Summary odds ratio (OR) and 95% confidence interval (95%CI) were calculated to appraise the association. Statistical analysis was performed using RevMan 5.3 software.
A total of seven case-control studies with 3129 cases and 3289 controls were recruited in the meta-analysis. For rs910080, rs1997794, rs1022563, and rs2235749 polymorphisms of PDYN gene, there were six, four, five, and four studies eventually included, respectively. The findings indicated that rs910080 polymorphism was significantly correlated with OD among Asian population under allelic model (A vs. G, OR = 1.30, 95% CI 1.04–1.62, P = 0.02, FDR = 0.05) and dominant model (AA+AG vs. GG, OR = 1.25, 95% CI 1.04–1.51, P = 0.02, FDR = 0.05). However, rs1022563, rs1997794 and rs2235749 polymorphisms did not appear to associate with OD susceptibility.
There existed a significant association between rs1022563 polymorphism and OD among Asian population. As the included studies were not adequate to guarantee a robust and convincing conclusion, future studies with larger sample size among more ethnicities are recommended.
Opioid dependence (OD), mainly characterized by persistent drug-taking and drug-seeking behavior, is a chronic brain disorder that has seriously negative social and health consequences . A range of environmental factors like severe stressors, family separation, divorce, as well as death in the family have been considered to be associated with OD for a long time . However, researchers observe that only a few individuals that were exposed to drugs will finally develop specific addictions . Thus, there are some factors beyond environmental factors that can contribute to OD.
Actually, early in 1990s, the isolation of genes that encoded opioid peptide precursors had opened an era of molecular and genetic investigations of OD . Family-based studies in different countries and regions revealed the heritability estimates for substance dependence of twins ranged from 30 to 70% . Over the past years, numerous genes and single-nucleotide polymorphisms have been reported to be contributors of OD [6, 7].
Prodynorphin (PDYN), also known as preprodynorphin, is the precursor of neoendorphins and dynorphins . The human PDYN gene is mapped to chromosome 20p13 . PDYN gene knockout mice presented an increasingly explorative behavior in anxiety tests, which demonstrated the role of prodynorphin-derived peptides .
Since Zimprich et al.  initially reported the correlation between PDYN gene polymorphisms and OD, several researchers have subsequently replicated with inconsistent results [12,13,14,15,16,17,18]. The lack of convergent findings might ascribe to weak statistical power to ascertain a significant effect, population structure or between-study heterogeneity. Thus, this meta-analysis was conducted to provide a more precise estimation on the relationship of PDYN gene polymorphisms and OD.
The procedures of the systematic review and meta-analysis were conducted in accordance with the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines .
Search strategy of literature
A systematic search of potentially relevant literature was carried out in PubMed, Web of Science, EMBASE, and Wanfang databases. The following search string was employed: (diamorphine OR opioids OR narcotic OR opium OR opiate OR heroin OR “Heroin”[Mesh]) AND (PDYN OR prodynorphin) AND (“Polymorphism, Single Nucleotide” [Mesh] OR SNP OR mutation OR polymorphism OR variant). No restriction on language was imposed. The search period was set from inception to August, 2018.
Inclusion and exclusion criteria
Publications that were included in the present study had to satisfy the following criteria: (1) case-control studies looking into the association of PDYN gene polymorphisms and OD; (2) case participants with confirmed diagnosis of OD; (3) studies with sufficient genotype counts to estimate odds ratio (OR) as well as 95% confidence interval (95%CI). Accordingly, case-case studies, reviews, conference abstracts, as well as animal studies were excluded.
The quality of eligible studies was evaluated separately by two investigators (CW and MM) according to Newcastle-Ottawa Scale (NOS). A ‘star system’ was used to judge the quality of each study on the basis of selection, comparability, as well as exposure. The NOS ranges from zero to nine stars. Studies with five or more stars indicated a high quality.
The following data were extracted independently by two investigators (CW and MM) from each qualified article: author, year of publication, country, ethnicity and gender of participants, sample size, and genotype distribution. Once encountering discrepancies, two authors re-checked the original articles until an agreement was achieved.
Hardy-Weinberg Equilibrium (HWE)  of control subjects was calculated. The strength of association of PDYN gene polymorphisms and OD was presented as OR and 95%CI in allelic model (M vs. W), homozygous model (MM vs. WW), heterozygous model (MW vs. WW), dominant model (MM + MW vs. WW), as well as recessive model (MM vs. MW + WW). To evaluate the heterogeneity across included studies, Q-statistical test alongside I2 test were utilized. P < 0.1 and I2 > 50% suggested the heterogeneity was not significant, and the fixed-effect model could be applied, otherwise the random-effects model should be employed. Subgroup-analysis by ethnicity was performed to test if there was an ethnicity-specific effect. For multiple testing corrections, the Benjamini-Hochberg method was applied to control false discovery rate (FDR) . Forest plots were made by RevMan 5.3 software. Funnel plots were exploited to assess the publication bias.
An initial search of four databanks yielded 151 records, of which 44 from Web of Science, 78 from EMBASE, 21 from PubMed, and 8 from Wanfang. Forty-seven records were deleted because of duplication. Ninety-three irrelevant citations were removed from the remaining 104 records after title and abstract screening. Another 4 ineligible articles were erased after comprehensive assessment according to inclusion and exclusion criteria. Ultimately, 7 articles [12,13,14,15,16,17,18] were included in qualitative synthesis and meta-analysis. The process was displayed in Fig. 1.
The general characteristics of included studies were displayed in Table 1. A total of seven case-control studies [12,13,14,15,16,17,18] with 3129 cases and 3289 controls were included in the meta-analysis. Of the included studies, six [12,13,14,15,16, 18] were published in English and one  was in Chinese. It should be noted that Clarke et al.’s study  published in 2012 consisted of two different cohorts. The sample size of an individual cohort ranged from 119 to 1030 for cases, and 90 to 668 for controls. For rs910080 polymorphism, Hashemi et al.’s study  was not in HWE. When it came to rs1022563 polymorphism, three studies [16,17,18] were out of HWE. As to rs2235749 polymorphism, three studies [14, 16, 17] were not in HWE. According to NOS, each study received five or more stars (Table 2).
Results of meta-analysis
The correlations of PDYN gene polymorphisms (rs910080, rs1997794, rs1022563, and rs2235749) and OD vulnerability were summarized in Table 3.
Data synthesis for rs910080 was available on six studies [12,13,14,15,16,17] with seven independent cohorts including 3092 cases and 3136 controls. Due to considerable heterogeneity, random-effect model was employed with the exception of heterozygous model. The pooled data indicated a null correlation between rs910080 polymorphism and OD susceptibility under five genetic models. Subgroup analysis stratified by ethnicity suggested a significant association between rs910080 polymorphism and OD across Asian population in allelic model (A vs. G, OR = 1.30, 95% CI 1.04–1.62, P = 0.02, FDR = 0.05) and dominant models (AA+AG vs. GG, OR = 1.25, 95% CI 1.04–1.51, P = 0.02, FDR = 0.05). While no correlation was detected under other models.
As for rs1997794, four studies [12, 13, 15, 18] with five cohorts involving 2942 cases and 2605 controls went into final data synthesis. Because heterogeneity between studies was significant, the random-effect model was applied. The result indicated that rs1997794 polymorphism did not significantly associate with OD in five genetic models. Subgroup analysis stratified by ethnicity also yielded a null association.
In respect of rs1022563, five studies with six cohorts containing 2378 cases and 2400 controls satisfied the eventual data combination. Because of huge heterogeneity, random-effect model was utilized. Rs1022563 polymorphism appeared to have no significant correlation with OD in five genetic models. Subgroup analysis by ethnicity as well indicated a non-significant association.
When it came to rs2235749, four studies [12, 14, 16, 17] with 1273 cases and 1285 controls carried out among Asian population were analyzed. OR and 95%CI were estimated by the random-effect approach due to significant heterogeneity. The pooled results indicated a null association between rs2235749 polymorphism and OD susceptibility.
Sensitivity analysis and publication bias
After studies that deviated from HWE were excluded, the corresponding OR did not reverse. Consequently, they were not removed from our meta-analysis. Sensitivity analysis confirmed the stability of the outcomes, where exclusion of individual dataset did not significantly change the overall effects. Furthermore, visual inspection of funnel plots did not identify substantial asymmetry for all the four loci under any genetic model.
OD affects the lives of millions of people worldwide and imposes a heavy burden on society . Hence, seeking new strategies to prevent and treat OD is necessary. Investigation into the etiology of OD has focused on multiple areas. Genetic factors and environmental factors are both considered as contributors to OD. As a genetic factor, PDYN gene has received more and more attention .
Several studies across different populations have been performed on the association of PDYN gene polymorphisms and OD vulnerability, while it is quite regrettable that no convincing evidence is available. This inconsistency could be interpreted by several factors, including severity of diagnosis, limited sample size, inadequate power, ethnic heterogeneity, population stratification, as well as large phenotype range. The present study, with pooled data from Asian (Malaysia, Iran, and China), Caucasian (European American), and African (African American) populations, indicated that rs910080 polymorphism was significantly correlated with OD among Asian population. However, rs1022563, rs1997794 and rs2235749 polymorphisms did not appear to associate with OD susceptibility.
PDYN gene codes ligand of opioid receptor. Functional variant in PDYN gene hampers the downstream signaling pathway and affects the overall biological networks, thus implicating in the pathophysiology of OD. Dynorphins are a cluster of posttranslational products of PDYN gene. They bind to all three types of opiate receptors, but exhibit considerable affinity to kappa opioid receptor . Products of PDYN gene have been reported to inhibit neurotransmission and take part in mood regulation, stress responses , reward , and motor functions. Yuferov et al.  reported the mRNA of PDYN was rich in brain areas that were associated with drug-taking and drug-withdrawal. Wang et al.  found that acute intermittent morphine administration could raise PDYN and kappa opioid receptor mRNA level in the brains of rat models. Solecki et al.  observed that chronic exposure to heroin could lead to an increased PDYN mRNA in the nucleus accumbens shell and core of rats. All this evidence suggested that PDYN gene played a crucial role in the occurrence and development of OD.
Heterogeneity is a crucial issue when elucidating the outcomes of a meta-analysis. We attempted to minimize this issue by developing strict criteria for study inclusion, information extraction, as well as data analysis. Considerable heterogeneity between studies still existed despite subgroup analysis was stratified. Several potential reasons might give rise to the presence of heterogeneity. For demographics, age of individual study ranges greatly. For gender, although majority studies are mixed genders, some are focusing on only men or women. For genetic background, included studies are multi-national and many ethnic groups were represented within them. For environmental factor, participants from both developing countries and developed countries were incorporated into the present study.
Although the current study has collected available data on the association between PDYN gene polymorphisms and OD predisposition, several limitations could not be neglected. First, the present outcomes were based on unadjusted estimates. Because raw data for confounding factors like age , gender , socioeconomic status  of the participants were unavailable, we could not perform the corresponding stratification analysis. Second, the number of participants of each included study was limited, thus lacking enough statistic power to guarantee the association. Third, we only investigated the role of four independent polymorphisms in PDYN gene. However, finding susceptibility polymorphisms and genes was a difficult task owing to the intricate mode of inheritance and multiple polymorphisms in various genes contributing to increased disease risk. OD, without exception, is a multifactor disorder determined by the synthetic effects and interactions of a series of polymorphisms, genes, and environmental factors. These interactions might conceal or enlarge the function of the included polymorphisms. Lastly, only seven articles from four databanks were searched for this study, potentially relevant studies published in other databanks might be missed.
Taken together, the current study indicated the rs910080 polymorphism was significantly correlated with OD among Asian population. However, rs1022563, rs1997794 and rs2235749 polymorphisms did not appear to associate with OD susceptibility. Due to limitations demonstrated above, the associations between PDYN gene polymorphisms and the risk of OD could not be entirely concluded. As the included studies were not adequate to guarantee a robust and convincing conclusion, studies with larger sample size among more populations are recommended.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author Ru-qin Luo, firstname.lastname@example.org.
- PDYN :
Cami J, Farre M. Drug addiction. N Engl J Med. 2003;349(10):975–86.
Tyrfingsson T, Thorgeirsson TE, Geller F, Runarsdottir V, Hansdottir I, Bjornsdottir G, et al. Addictions and their familiality in Iceland. Ann N Y Acad Sci. 2010;1187:208–17.
Trigo JM, Martin-Garcia E, Berrendero F, Robledo P, Maldonado R. The endogenous opioid system: a common substrate in drug addiction. Drug Alcohol Depend. 2010;108(3):183–94.
Charbogne P, Kieffer BL, Befort K. 15 years of genetic approaches in vivo for addiction research: Opioid receptor and peptide gene knockout in mouse models of drug abuse. Neuropharmacology. 2014;76(Pt B):204–17.
Agrawal A, Lynskey MT. Are there genetic influences on addiction: evidence from family, adoption and twin studies. Addiction. 2008;103(7):1069–81.
Prom-Wormley EC, Ebejer J, Dick DM, Bowers MS. The genetic epidemiology of substance use disorder: a review. Drug Alcohol Depend. 2017;180:241–59.
Goldman D, Oroszi G, Ducci F. The genetics of addictions: uncovering the genes. Nat Rev Genet. 2005;6(7):521–32.
Yuferov V, Levran O, Proudnikov D, Nielsen DA, Kreek MJ. Search for genetic markers and functional variants involved in the development of opiate and cocaine addiction and treatment. Ann N Y Acad Sci. 2010;1187:184–207.
Litt M, Buroker NE, Kondoleon S, Douglass J, Liston D, Sheehy R, et al. Chromosomal localization of the human proenkephalin and prodynorphin genes. Am J Hum Genet. 1988;42(2):327–34.
Wittmann W, Schunk E, Rosskothen I, Gaburro S, Singewald N, Herzog H, et al. Prodynorphin-derived peptides are critical modulators of anxiety and regulate neurochemistry and corticosterone. Neuropsychopharmacology. 2009;34(3):775–85.
Zimprich A, Kraus J, Woltje M, Mayer P, Rauch E, Hollt V. An allelic variation in the human prodynorphin gene promoter alters stimulus-induced expression. J Neurochem. 2000;74(2):472–7.
Yuanyuan J, Rui S, Hua T, Jingjing C, Cuola D, Yuhui S, et al. Genetic association analyses and meta-analysis of Dynorphin-kappa opioid system potential functional variants with heroin dependence. Neurosci Lett. 2018;685:75–82.
Nagaya D, Zahari Z, Saleem M, Yahaya BH, Tan SC, Yusoff NM. An analysis of genetic association in opioid dependence susceptibility. J Clin Pharm Ther. 2018;43(1):80–6.
Hashemi M, Shakiba M, Sanaei S, Shahkar G, Rezaei M, Mojahed A, et al. Evaluation of prodynorphin gene polymorphisms and their association with heroin addiction in a sample of the southeast Iranian population. Mol Biol Res Commun. 2018;7(1):1–6.
Clarke TK, Ambrose-Lanci L, Ferraro TN, Berrettini WH, Kampman KM, Dackis CA, et al. Genetic association analyses of PDYN polymorphisms with heroin and cocaine addiction. Genes Brain Behav. 2012;11(4):415–23.
Wei SG, Zhu YS, Lai JH, Xue HX, Chai ZQ, Li SB. Association between heroin dependence and prodynorphin gene polymorphisms. Brain Res Bull. 2011;85(3–4):238–42.
Jia W, Sh JG, Ao L, Lei A, Zhang R. Association between polymorphism in 3′ untranslated regions of prodynorphin gene and heroin dependence. Progress Modern Biomed. 2010;10(13):2439–41.
Clarke TK, Krause K, Li T, Schumann G. An association of prodynorphin polymorphisms and opioid dependence in females in a Chinese population. Addict Biol. 2009;14(3):366–70.
Moher D, Liberati A, Tetzlaff J, Altman DG. Preferred reporting items for systematic reviews and meta-analyses: the PRISMA statement. BMJ. 2009;339:b2535.
Ryckman K, Williams SM. Calculation and use of the Hardy-Weinberg model in association studies. Curr Protoc Hum Genet. 2008;Chapter 1:Unit 1.18.
Benjamini Y, Yekutieli D. Quantitative trait loci analysis using the false discovery rate. Genetics. 2005;171(2):783–90.
Darcq E, Kieffer BL. Opioid receptors: drivers to addiction? Nat Rev Neurosci. 2018;19(8):499–514.
Reed B, Butelman ER, Yuferov V, Randesi M, Kreek MJ. Genetics of opiate addiction. Curr Psychiatry Rep. 2014;16(11):504.
Schwarzer C. 30 years of dynorphins--new insights on their functions in neuropsychiatric diseases. Pharmacol Ther. 2009;123(3):353–70.
Drolet G, Dumont EC, Gosselin I, Kinkead R, Laforest S, Trottier JF. Role of endogenous opioid system in the regulation of the stress response. Prog Neuro-Psychopharmacol Biol Psychiatry. 2001;25(4):729–41.
Kelley AE, Bakshi VP, Haber SN, Steininger TL, Will MJ, Zhang M. Opioid modulation of taste hedonics within the ventral striatum. Physiol Behav. 2002;76(3):365–77.
Yuferov VP, Laforge KS, Spangler R, Maggos CE, Kreek MJ. Guinea pig preprodynorphin mRNA: primary structure and regional quantitation in the brain. DNA Cell Biol. 1996;15(12):1105–12.
Wang XM, Zhou Y, Spangler R, Ho A, Han JS, Kreek MJ. Acute intermittent morphine increases preprodynorphin and kappa opioid receptor mRNA levels in the rat brain. Brain Res Mol Brain Res. 1999;66(1–2):184–7.
Solecki W, Ziolkowska B, Krowka T, Gieryk A, Filip M, Przewlocki R. Alterations of prodynorphin gene expression in the rat mesocorticolimbic system during heroin self-administration. Brain Res. 2009;1255:113–21.
Yuferov V, Randesi M, Butelman ER, van den Brink W, Blanken P, van Ree JM, et al. Association of variants of prodynorphin promoter 68-bp repeats in caucasians with opioid dependence diagnosis: effect on age trajectory of heroin use. Neurosci Lett. 2019;704:100–5.
Pear VA, Ponicki WR, Gaidus A, Keyes KM, Martins SS, Fink DS, et al. Urban-rural variation in the socioeconomic determinants of opioid overdose. Drug Alcohol Depend. 2019;195:66–73.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.