Skip to main content

Association studies of dopamine synthesis and metabolism genes with multiple phenotypes of heroin dependence

Abstract

Background

Heroin dependence is a complex disease with multiple phenotypes. Classification of heroin users into more homogeneous subgroups on the basis of these phenotypes could help to identify the involved genetic factors and precise treatments. This study aimed to identify the association between genetic polymorphisms of DA synthesis and metabolism genes, including tyrosine hydroxylase (TH), DOPA decarboxylase (DDC), solute carrier family 6 member 3 (SLC6A3) and DA beta-hydroxylase (DBH), with six important phenotypes of heroin dependence.

Methods

A total of 801 heroin dependent patients were recruited and fourteen potential functional single nucleotide polymorphisms (SNPs) were genotyped by SNaPshot. Associations between SNPs with six phenotypes were mainly assessed by binary logistic regression. Generalized multifactor dimensionality reduction was used to analyze the gene-by-gene and gene-by-environment interactions.

Results

We found that DBH rs1611114 TT genotype had a protective effect on memory impairment after heroin dependence (P = 0.002, OR = 0.610). We also found that the income-rs12666409-rs129915-rs1611114 interaction yielded the highest testing balance accuracy and cross-validation consistency for memory change after heroin dependence.

Conclusions

Our results suggest that the memory change after heroin dependence was a result of a combination of genetics and environment. This finding could lead to a better understanding of heroin dependence and further improve personalized treatment.

Peer Review reports

Background

Drug dependence is a chronic relapsing brain disease that results from a combination of genetic and environmental factors [1], and the genetic factor accounts for approximately 30–60% of the total risk [2, 3]. Although multiple genes for the susceptibility to heroin dependence have been identified, a lot of these susceptibility genes could not be confirmed in subsequent studies, which further reminded us of the complexity of the genetic mechanism of heroin dependence. Heroin dependence is a complex disease with multiple phenotypes [4]. Classification of heroin users into more homogeneous subgroups on the basis of clinical and/or pathophysiological features could help to identify the involved genetic factors and precise treatments [4]. A previous study showed no significant association between the dopamine (DA) receptor D1 gene and opioid dependence in the overall analysis, but the association became apparent after stratifying based on the duration of the transition from first use to dependence (DTFUD) [5]. The DTFUD determined the addictive liability of the dependent patients [5, 6]. Early opioid use was associated with higher rates of impulsivity and relapse of drug seeking behaviour [7, 8]. Individuals have different tolerances to heroin [9] or methadone [10] in terms of the amount required to achieve the same efficacy. Therefore, individuals require different doses of heroin or methadone. Heroin dependence leads to pathological memory formation and seriously affects normal working memory [11]. These phenotypes of heroin dependence vary widely among heroin dependent individuals and may be affected by genetic factors.

DA system modulates motivation, cognition and reward in the context of dependence [12,13,14]. Extracellular DA concentrations are the net result of release (exocytosis) and clearance (uptake) from the extracellular space. Tyrosine hydroxylase (TH) converts L-tyrosine to L-DOPA, and then DOPA decarboxylase (DDC) catalyzes the decarboxylation of L-DOPA to generate DA [15]. The solute carrier family 6 member 3 (SLC6A3) encodes a DA transporter (DAT) that can re-uptake DA released into the synaptic cleft to the presynaptic terminal [16]. The DA beta-hydroxylase (DBH) protein catalyses the oxidative hydroxylation of DA to norepinephrine [17]. These factors regulate extracellular DA concentrations at different levels.

The A allele of SLC6A3 rs27072 was associated with an early-onset age of nicotine dependence [18]. The frequency of the G allele of TH rs6356 was significantly lower in the late-onset Parkinson’s disease group [19]. Male ever-smokers with the AA genotype of DBH rs5320 smoked fewer cigarettes per day [20]. Haplotype-based association analysis revealed a protective T-G-T-G haplotype of DDC rs921451-rs3735273-rs1451371-rs2060762, which was significantly associated with nicotine dependence after correction for multiple testing [21]. However, the association of these genes with other important phenotypes of heroin dependence remained unclear.

We selected six phenotypes: memory change after heroin dependence, age of onset for heroin use, euphoria associated with first heroin use, DTFUD, daily heroin dose and daily methadone dose. We identified 14 potential functional single nucleotide polymorphisms (SNPs) by bioinformatics. We then analyzed the associations between the SNPs located in SLC6A3 (rs10064525, rs27072, rs1042098 and rs6347), TH (rs10770140, rs10770141, rs3842727 and rs6356), DDC (rs11575553 and rs12666409) and DBH (rs129882, rs129915, rs1611114 and rs5320) with these phenotypes. We performed survival analysis for DTFUD. Finally, we analyzed the gene-by-gene and gene-by-environment interactions in these significant phenotypes.

Methods

Subjects

A total of 801 heroin-addicted individuals (mean age 39.30 ± 8.291 years) participated in our study. All of these participants were Chinese Han population. They were all unrelated individuals that were neither first-degree relatives, second-degree relatives nor third-degree relatives. This study included 709 males and 92 females. These participants were recruited from Lantian Compulsory Isolated and Detoxification Center, Xi’an Methadone Maintenance Treatment (MMT) Center and Xin’an Central Hospital Drug Rehabilitation Center. Heroin dependence was evaluated according to the Diagnostic and Statistical Manual of Mental Disorders fifth edition (DSM-V) criteria and a urine test. A structured interview was performed to assist with the diagnosis and to ask questions regarding (1) basic information: age, gender, nationality, education level, employment status, individual income, marital status and family atmosphere; (2) heroin use: the age of onset for heroin use, DTFUD, the daily dose of heroin used, reasons for first use of heroin, the administration routes of heroin, the euphoria associated with the first heroin use and the euphoria associated with heroin use after dependence; (3) individual responses after heroin use: body weight, sleep quality, memory, appetite and sexual desire; (4) methadone maintenance therapy or withdrawal symptoms: the number of withdrawal and relapse for heroin (absence of heroin use over 7 days counted as one withdrawal), the daily methadone dose used, and the methadone maintenance treatment duration. The detailed content of the questionnaire was shown in Supplementary file 1. Exclusion criteria were as follows: age < 18 years old; dependence on other substances; presence of other neuropsychiatric diseases; and involvement in other clinical trials. A total of 737 (92%) subjects were tobacco smokers in this study. All smokers consumed< 20 cigarettes per day and did not meet the criteria for nicotine dependence according to the Fagerström Test for Nicotine Dependence (FTND, score of 4 or lower defined as no dependence).

Written informed consent was obtained from all the study participants. The Ethical Committee of Xi’an Mental Health Center approved our study. We followed the approved guidelines to perform the protocol.

SNP selection

SNP selection criteria were as follows: (1) located in the promoter, intron-exon border, exon, the 3’near region and untranslated regions (UTRs) of gene; (2) the minor allele frequency (MAF) greater than 0.05 based on dbSNP Han Chinese Beijing (HCB) and HapMap databases; (3) each SNP had three genotypes. Then 5 SNPs in SLC6A3 (rs10064525, rs27072, rs1042098 and rs6347), 4 SNPs in TH (rs10770140, rs10770141, rs3842727 and rs6356), 2 SNPs in DDC (rs11575553 and rs12666409) and 4 SNPs in DBH (rs129882, rs129915, rs1611114 and rs5320) were selected. Some important information of these SNPs was summarized in Table 1. Potential functional effects of these SNPs were predicted by SNP Function Prediction (FuncPred, https://snpinfo.niehs.nih.gov/snpinfo/snpfunc.html). Detailed effects of these SNPs were shown in Table S1 of Supplementary file 2.

Table 1 The information of selected SNPs

Genotyping

The genomic DNA was extracted by E.Z.N.A™ Blood DNA Midi Kit (Omega Bio-TeK, Norcross, GA, USA). Then, we used SNaPshot technology to genotype these SNPs. These specific experimental steps can be found in our previous research [22].

Statistical analysis

Hierarchical data on phenotype variables were used to identify subtypes, and associations between these subtypes with SNPs were first assessed by Pearson’s chi-square test or Fisher’s exact test. Then, these associations were assessed again by binary logistic regression, and the basic information (age, gender, education level, employment status, individual income, marital status and family atmosphere) were set as covariates.

We assessed the potential correlation of DTFUD with these genotypes using Kaplan–Meier survival analysis. Survival curves were used to estimate the probability that individuals had not experienced dependence over a period of time following their first heroin use. The DTFUD was regarded as survival time in months. The origin of the first use was specified as the time of first heroin use, and the outcome of interest was the first occurrence of dependence. We used three methods (Wilcoxon, log rank and Tarone–Ware tests) to compare the survival curves, and these methods give varying weight to different phases of follow-up time.

We used generalized multifactor dimensionality reduction (GMDR; http://sourceforge.net/projects/gmdr/) to analyze the gene-by-gene and gene-by-environment interactions. This method permits adjustment for discrete and quantitative covariates and is applicable to both dichotomous and continuous phenotypes in various population-based study designs [23]. The multi-locus model with the maximum testing balance accuracy (Testing Bal. Acc) and cross-validation consistency (CV Consistency) was the best interaction model.

P values were calculated based on our previous research [22] P values were presented as two-sided, and the statistically significant was defined as P < 0.05. We used Bonferroni’s correction to adjust the significance level, and the P value was divided by the total number of loci analyzed (15). Consequently the loci should be considered significant at P < 0.003. All statistical analyses were performed using SPSS 20.0 software (SPSS, USA).

Results

Association of DBH with memory change after dependence

For the phenotype of memory, a total of 586 subjects (73.2%) that memory impaired were assigned to the “reduction” group, and 215 subjects (26.8%) that memory no changed were assigned to the “no difference” group.

The frequency of DBH rs1611114 TT genotype was significantly lower in the reduction group (P = 0.002, OR = 0.610, 95%CI = 0.445–0.835, Table 2). This association was still significantly after adjustment of covariates (P = 0.002, OR = 0.607, 95%CI = 0.440–0.837, Table 2). These associations of other SNPs with memory change after dependence were listed in Table S2 of Supplementary file 2.

Table 2 The association between SNP and memory change after dependence

SNPs and the age of onset for heroin use

The age of onset for heroin use was converted to binary using median (28 years old), the sample was divided into earlier-onset (n = 440, 54.9%) and later-onset (n = 361, 45.1%) groups accordingly.

After Bonferroni’s correction, TH rs10770140 and rs10770141 were not associated with the age of onset for heroin use (P = 0.018 and P = 0.015, respectively, Table 3). These patients who carrying the rs10770140 TT genotype (P = 0.015, OR = 0.599, 95%CI = 0.395–0.908, Table 3) or the T allele (P = 0.013, OR = 0.604, 95%CI = 0.405–0.901, Table 3) were later-onset, but these findings were not observed after Bonferroni’s correction. These patients who carrying the rs10770141 GG genotype (P = 0.012, OR = 0.577, 95%CI = 0.374–0.890, Table 3) or the G allele (P = 0.010, OR = 0.580, 95%CI = 0.382–0.881, Table 3) were later-onset, but these findings were not observed after Bonferroni’s correction.

Table 3 The associations between SNPs and age of onset for heroin use

After Bonferroni’s correction, these associations of rs10770140 TT (P = 0.019, OR = 0.568, 95%CI = 0.354–0.912, Table 3) and rs10770141 GG (P = 0.014, OR = 0.538, 95%CI = 0.328–0.881, Table 3) genotypes with age of onset were still not significant after adjustment of covariates. These associations of other SNPs with age of onset for heroin use were listed in Table S3 of Supplementary file 2.

SNPs and the euphoria of the first heroin use

For the phenotype of euphoria, a total of 173 subjects (21.6%) that have obvious feelings of euphoria were assigned to the “strong euphoria” group, and 628 subjects (78.4%) that have minimal level of happiness, no euphoria or unclear were assigned to the “week euphoria” group.

After Bonferroni’s correction, TH rs10770140 (P = 0.047, Table 4) and rs10770141 (P = 0.035, Table 4) were not associated with the intensity of euphoria. The subjects carrying the rs10770140 CC (P = 0.046, OR = 1.012, 95%CI = 0.996–1.028, Table 4) or rs10770141 AA (P = 0.046, OR = 1.012, 95%CI = 0.996–1.028, Table 4) genotype had a stronger euphoric response to heroin, but this effect was not observed after Bonferroni’s correction. These associations were not significant after adjustment for the covariates. These associations of other SNPs with the intensity of euphoria were listed in Table S4 of Supplementary file 2.

Table 4 The associations between SNPs and euphoria associated with the first heroin use

SNPs and the daily dose of methadone or heroin

For the phenotype of methadone dose, the sample was divided into a high methadone dose group (n = 285, 44.9%) and a low methadone dose group (n = 350, 55.1%) based on the median daily methadone dose (45 mg). For the phenotype of heroin dose, the sample was divided into a high heroin dose group (n = 375, 46.8%) and a low heroin dose group (n = 426, 53.2%) based on the median daily heroin dose (0.5 g). No SNP was associated with the daily dose of methadone or heroin use (Table S5 and S6 of Supplementary file 2, respectively).

SNPs and DTFUD

For the phenotype of DTFUD, the sample was divided into a long duration group (n = 376, 46.9%) and a short duration group (n = 425, 53.1%) based on the median of DTFUD (1 month). No SNP was associated with DTFUD in either Pearson’s chi-square test (Table S7 of Supplementary file 2) or the Kaplan–Meier survival analysis (Fig. S1 of Supplementary file 3).

Gene-by-gene and gene-by-environment interactions

The one- to four-way gene-by-gene and gene-by-environment interaction models of the memory change after dependence were analyzed by GMDR.

For the phenotype of memory change after dependence, we found that these interactions in the one- to four-way models were significant. In two-way model (income-rs3842727), income and rs3842727 have interaction. In three-way model (income-rs129915-rs1611114), income, rs129915 and rs1611114 have interaction. In four-way model (income-rs12666409-rs129915-rs1611114), income, rs12666409, rs129915 and rs1611114 have interaction. Among these significant interactions, the CV Consistency and Testing Bal. Acc in the four-way model (income-rs12666409-rs129915-rs1611114) were the most significant, this model was the best model (Table 5).

Table 5 Gene-by-gene and gene-by-environment interactions in phenotype of memory change after dependence

These sequencing data of SNPs were listed in Supplementary file 4.

Discussion

Heroin dependence is a heterogeneous disease, and its clinical features and therapeutic outcomes vary substantially among patients. Regardless of linkage analysis or associative analysis, the results are obviously affected by disease grouping within samples, which finally result in poor stability. Features that correlated with dependence could be used as indicators to categorize heroin dependent patients into subgroups with uniform disease features.

In this study, we demonstrated that DBH rs1611114 was significantly associated with memory change after dependence. SNP rs1611114 was located in the 5’near region of DBH and was predicted to be located in the transcription factor binding site (Table S8 of Supplementary file 2). Research found that DBH − 1021C/T variant accounted for 35–52% of the variation in plasma DBH activity, and the T allele predicted very low DBH activity [24]. Studies have also found that lower DBH expression levels lead to higher levels of extracellular DA and therefore have better working memory [25, 26]. Hippocampal cannabinoid receptor 1, found exclusively in DBH-expressing cells, was both necessary and sufficient to impair memory consolidation produced by acute stress [27]. DBH also modulated working memory and verbal and visual memory and sustained attention during the cannabis intoxication period [28]. These studies shown that rs1611114 may affect memory phenotype by regulating DBH expression. Similar to these findings, we found that the frequency of the rs1611114 TT in the memory reduction group was significantly lower, which indicated that subjects in the memory reduction group may have lower levels of extracellular DA and therefore have worse memory. Furthermore, income and TH rs3842727 exhibited interactions in the memory change phenotype. Income, DBH rs129915, rs1611114 exhibited interactions in this phenotype. Rs1611114, DDC rs1266640, DBH rs129915 and income also exhibited interactions in this phenotype. These results showed that the memory change phenotype was a result of a combination of genetics and environment. TH converts L-tyrosine to L-DOPA, DDC catalyzes the decarboxylation of L-DOPA to generate DA, and then DBH catalyses the oxidative hydroxylation of DA to norepinephrine. SNPs in these genes may interact with each other in regulating the concentration of DA, and then affected the memory change after heroin dependence. Research showed that the higher-income group had greater memory capacity [29]. Children from economically disadvantaged background showed lower memory performance and had a smaller hippocampal volume [30]. Lower income was associated with a flatter cortisol awakening response, blunted reactivity and recovery to acute stress [31]. Therefore, higher income may have a protective effect on memory impairment.

TH rs10770140 and rs10770141 were promoter variants that differentially affected the transcriptional activity of TH [32] (Table S9 and Table S10 of Supplementary file 2, respectively). Under basal circumstances and nicotine stimulation circumstances, the rs10770141 A and rs10770140 C alleles exerted a greater effect on TH transcription [32, 33]. However, after Bonferroni’s correction, these associations of rs10770140 C and rs10770141 A alleles with age of onset for heroin use were not significant. These associations of rs10770140 CC and rs10770141 AA genotypes with the intensity of euphoria were still not significant after Bonferroni’s correction. The TH Val-81-Met polymorphism was found to be associated with early-onset alcoholism [34], as the P value of this article was not corrected.

Some studies have found that SLC6A3 rs27072 may be related to the age of onset for tobacco or alcohol [35,36,37]. In particular, rs27072-A allele carriers were more likely to initiate smoking onset before 18 years old [36]. However, our research did not find that SLC6A3 rs27072 was related to the age of onset for heroin use. We proposed three hypotheses to explain the difference. First, our study excluded alcohol and nicotine dependent subjects. Second, the proportion of our sample that use heroin before 18 years old was relatively small (Fig. S2 of Supplementary file 3), and age of onset in our study was converted to binary using median (28 years old). Third, the mechanism of heroin dependence may different from that of alcohol and nicotine dependence.

Conclusions

In conclusion, our results supported the important role of DBH rs1611114 in memory change of heroin dependence. This phenotype was a result of a combination of genetics and environment. These findings were valuable for obtaining a better understanding of heroin dependence and improving personalized treatment.

Availability of data and materials

All the data supporting the conclusions of this article are presented within the manuscript or are available in the additional supporting file containing the supplementary material.

Abbreviations

DA:

Dopamine

DAT:

DA transporter

DBH:

DA beta-hydroxylase

DDC:

DOPA decarboxylase

DSM-V:

Diagnostic and Statistical Manual of Mental Disorders fifth edition

DTFUD:

Duration of the transition from first use to dependence

FTND:

Fagerström Test for Nicotine Dependence

GMDR:

Generalized multifactor dimensionality reduction

HCB:

Han Chinese Beijing

MAF:

Minor allele frequency

MMT:

Methadone Maintenance Treatment

SLC6A3:

Solute carrier family 6 member 3

SNPs:

Single nucleotide polymorphisms

TH:

Tyrosine hydroxylase

UTRs:

Untranslated regions

References

  1. Vanyukov MM, Tarter RE. Genetic studies of substance abuse. Drug Alcohol Depend. 2000;59(2):101–23.

    CAS  PubMed  Google Scholar 

  2. Kendler KS, Jacobson KC, Prescott CA, Neale MC. Specificity of genetic and environmental risk factors for use and abuse/dependence of cannabis, cocaine, hallucinogens, sedatives, stimulants, and opiates in male twins. Am J Psychiatry. 2003;160(4):687–95.

    PubMed  Google Scholar 

  3. Tsuang MT, Lyons MJ, Meyer JM, Doyle T, Eisen SA, Goldberg J, et al. Co-occurrence of abuse of different drugs in men: the role of drug-specific and shared vulnerabilities. Arch Gen Psychiatry. 1998;55(11):967–72.

    CAS  PubMed  Google Scholar 

  4. Sun J, Bi J, Chan G, Oslin D, Farrer L, Gelernter J, et al. Improved methods to identify stable, highly heritable subtypes of opioid use and related behaviors. Addict Behav. 2012;37(10):1138–44.

    PubMed  PubMed Central  Google Scholar 

  5. Zhu F, Yan CX, Wen YC, Wang J, Bi J, Zhao YL, et al. Dopamine D1 receptor gene variation modulates opioid dependence risk by affecting transition to addiction. PLoS One. 2013;8(8):e70805.

    CAS  PubMed  PubMed Central  Google Scholar 

  6. Ridenour TA, Maldonado-Molina M, Compton WM, Spitznagel EL, Cottler LB. Factors associated with the transition from abuse to dependence among substance abusers: implications for a measure of addictive liability. Drug Alcohol Depend. 2005;80(1):1–14.

    PubMed  PubMed Central  Google Scholar 

  7. Basu D, Ghormode D, Madan R, Mattoo S, Nehra R, Prabhakar S. Age of onset of dependence: does it help our understanding of opioid dependence by generating meaningful categories or by acting as a useful dimension? A critical examination of the classic debate in psychiatry. Indian J Psychiatry. 2014;56(3):228–37.

    PubMed  PubMed Central  Google Scholar 

  8. Chahua M, Sordo L, Barrio G, Domingo-Salvany A, Brugal MT, Molist G, et al. Non-fatal opioid overdose and major depression among street-recruited young heroin users. Eur Addict Res. 2014;20(1):1–7.

    PubMed  Google Scholar 

  9. Warner-Smith M, Darke S, Lynskey M, Hall W. Heroin overdose: causes and consequences. Addiction. 2001;96(8):1113–25.

    CAS  PubMed  Google Scholar 

  10. Bell J, Burrell T, Indig D, Gilmour S. Cycling in and out of treatment; participation in methadone treatment in NSW, 1990-2002. Drug Alcohol Depend. 2006;81(1):55–61.

    PubMed  Google Scholar 

  11. Yan WS, Li YH, Xiao L, Zhu N, Bechara A, Sui N. Working memory and affective decision-making in addiction: a neurocognitive comparison between heroin addicts, pathological gamblers and healthy controls. Drug Alcohol Depend. 2014;134:194–200.

    PubMed  Google Scholar 

  12. Baik JH. Dopamine signaling in reward-related behaviors. Front Neural Circuits. 2013;7:152.

    PubMed  PubMed Central  Google Scholar 

  13. Chen W, Nong Z, Li Y, Huang J, Chen C, Huang L. Role of dopamine signaling in drug addiction. Curr Top Med Chem. 2017;17(21):2440–55.

    CAS  PubMed  Google Scholar 

  14. Zilverstand A, Parvaz MA, Moeller SJ, Goldstein RZ. Cognitive interventions for addiction medicine: understanding the underlying neurobiological mechanisms. Prog Brain Res. 2016;224:285–304.

    PubMed  Google Scholar 

  15. Chang YT, Mues G, Hyland K. Alternative splicing in the coding region of human aromatic L-amino acid decarboxylase mRNA. Neurosci Lett. 1996;202(3):157–60.

    CAS  PubMed  Google Scholar 

  16. Giros B, Caron MG. Molecular characterization of the dopamine transporter. Trends Pharmacol Sci. 1993;14(2):43–9.

    CAS  PubMed  Google Scholar 

  17. Gonzalez-Lopez E, Vrana KE. Dopamine beta-hydroxylase and its genetic variants in human health and disease. J Neurochem. 2020;152(2):157–81.

    CAS  PubMed  Google Scholar 

  18. Ling D, Niu T, Feng Y, Xing H, Xu X. Association between polymorphism of the dopamine transporter gene and early smoking onset: an interaction risk on nicotine dependence. J Hum Genet. 2003;49:35.

    PubMed  Google Scholar 

  19. Yan YP, Zhang B, Shen T, Si XL, Guo ZY, Tian J, et al. Study of GCH1 and TH genes in Chinese patients with Parkinson's disease. Neurobiol Aging. 2018;68:159.e3–6.

    CAS  Google Scholar 

  20. Ella E, Sato N, Nishizawa D, Kageyama S, Yamada H, Kurabe N, et al. Association between dopamine beta hydroxylase rs5320 polymorphism and smoking behaviour in elderly Japanese. J Hum Genet. 2012;57(6):385–90.

    CAS  PubMed  Google Scholar 

  21. Ma JZ, Beuten J, Payne TJ, Dupont RT, Elston RC, Li MD. Haplotype analysis indicates an association between the DOPA decarboxylase (DDC) gene and nicotine dependence. Hum Mol Genet. 2005;14(12):1691–8.

    CAS  PubMed  Google Scholar 

  22. Li Y, Qiao X, Yin F, Guo H, Huang X, Lai J, et al. A population-based study of four genes associated with heroin addiction in Han Chinese. PLoS One. 2016;11(9):e0163668.

    PubMed  PubMed Central  Google Scholar 

  23. Xu HM, Xu LF, Hou TT, Luo LF, Chen GB, Sun XW, et al. GMDR: versatile software for detecting Gene-gene and Gene-environ- ment interactions underlying complex traits. Curr Genomics. 2016;17(5):396–402.

    CAS  PubMed  PubMed Central  Google Scholar 

  24. Zabetian CP, Anderson GM, Buxbaum SG, Elston RC, Ichinose H, Nagatsu T, et al. A quantitative-trait analysis of human plasma-dopamine beta-hydroxylase activity: evidence for a major functional polymorphism at the DBH locus. Am J Hum Genet. 2001;68(2):515–22.

    CAS  PubMed  PubMed Central  Google Scholar 

  25. Greenwood PM, Lin MK, Sundararajan R, Fryxell KJ, Parasuraman R. Healthy aging increases the cognitive effects of two genes that influence extracellular dopamine. Psychol Aging. 2014;29(2):363–73.

    PubMed  Google Scholar 

  26. Parasuraman R, Greenwood PM, Kumar R, Fossella J. Beyond heritability: neurotransmitter genes differentially modulate visuospatial attention and working memory. Psychol Sci. 2005;16(3):200–7.

    PubMed  PubMed Central  Google Scholar 

  27. Busquets-Garcia A, Gomis-Gonzalez M, Srivastava RK, Cutando L, Ortega-Alvaro A, Ruehle S, et al. Peripheral and central CB1 cannabinoid receptors control stress-induced impairment of memory consolidation. Proc Natl Acad Sci U S A. 2016;113(35):9904–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  28. Cosker E, Schwitzer T, Ramoz N, Ligier F, Lalanne L, Gorwood P, et al. The effect of interactions between genetics and cannabis use on neurocognition. A review Prog Neuropsychopharmacol Biol Psychiatry. 2018;82:95–106.

    CAS  PubMed  Google Scholar 

  29. Lee S, Buring JE, Cook NR, Grodstein F. The relation of education and income to cognitive function among professional women. Neuroepidemiology. 2006;26(2):93–101.

  30. Raffington L, Czamara D, Mohn JJ, Falck J, Schmoll V, Heim C, et al. Stable longitudinal associations of family income with children's hippocampal volume and memory persist after controlling for polygenic scores of educational attainment. Dev Cogn Neurosci. 2019;40:100720.

    PubMed  PubMed Central  Google Scholar 

  31. Raffington L, Prindle J, Keresztes A, Binder J, Heim C, Shing YL. Blunted cortisol stress reactivity in low-income children relates to lower memory function. Psychoneuroendocrinology. 2018;90:110–21.

    CAS  PubMed  Google Scholar 

  32. Rao F, Zhang K, Zhang L, Rana BK, Wessel J, Fung MM, et al. Human tyrosine hydroxylase natural allelic variation: influence on autonomic function and hypertension. Cell Mol Neurobiol. 2010;30(8):1391–4.

    CAS  PubMed  PubMed Central  Google Scholar 

  33. Zhang K, Zhang L, Rao F, Brar B, Rodriguez-Flores JL, Taupenot L, et al. Human tyrosine hydroxylase natural genetic variation: delineation of functional transcriptional control motifs disrupted in the proximal promoter. Circ Cardiovasc Genet. 2010;3(2):187–98.

    CAS  PubMed  PubMed Central  Google Scholar 

  34. Dahmen N, Volp M, Singer P, Hiemke C, Szegedi A. Tyrosine hydroxylase Val-81-met polymorphism associated with early-onset alcoholism. Psychiatr Genet. 2005;15(1):13–6.

    PubMed  Google Scholar 

  35. Schmid B, Blomeyer D, Becker K, Treutlein J, Zimmermann US, Buchmann AF, et al. The interaction between the dopamine transporter gene and age at onset in relation to tobacco and alcohol use among 19-year-olds. Addict Biol. 2009;14(4):489–99.

    CAS  PubMed  Google Scholar 

  36. Ling D, Niu T, Feng Y, Xing H, Xu X. Association between polymorphism of the dopamine transporter gene and early smoking onset: an interaction risk on nicotine dependence. J Hum Genet. 2004;49(1):35–9.

    CAS  PubMed  Google Scholar 

  37. Genetic Study. The interaction between the dopamine transporter gene and age at onset in relation to tobacco and alcohol use among 19-year-olds. Addict Biol. 2009;14(4):489–99.

    Google Scholar 

Download references

Acknowledgments

We warmly thank the participants for their contribution.

Funding

This work was supported by the National Natural Science Foundation of China (No. 81571856). The funding bodies played no direct role in the design of the study, the collection, analysis, or interpretation of data, or the writing of the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

YL, YZ, SW and JL were responsible for the study concept and design; YL and YZ contributed to the acquisition of phenotype data and sample; YL, YZ, XS, YC and JZ performed the genotyping; YL and YZ assisted with data analysis and interpretation of findings; YL drafted the manuscript; YL, YZ, SW and JL provided critical revision of the manuscript for important intellectual content. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Shuguang Wei.

Ethics declarations

Ethics approval and consent to participate

Written informed consent was obtained from all the study participants. All participants read the informed consent, and we explained the relevant contents of the informed consent to them. Then all participants signed to agree to the relevant contents. This study protocol was approved by the Ethical Committee of the Medical College, Xi’an Jiaotong University.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1:.

Supplementary file 1: Questionnaire.

Additional file 2:

Supplementary file 2: Table S1 The functional effects of all selected SNPs. Table S2 Results of association analysis between other SNPs and phenotype of memory. Table S3 Results of association analysis between other SNPs and age of onset for heroin use. Table S4 Results of association analysis between other SNPs and phenotype of euphoria. Table S5 Results of association analysis between SNPs and phenotype of methadone dosage. Table S6 Results of association analysis between SNPs and phenotype of heroin dosage. Table S7 Results of association analysis between SNPs and phenotype of DTFUD. Table S8 Transcription factor binding sites of rs1611114. Table S9 Transcription factor binding sites of rs10770140. Table S10 Transcription factor binding sites of rs10770141.

Additional file 3:.

Supplementary file 3: Figure. S1. Kaplan–Meier survival analysis of DTFUD. Figure. S2. Frequency distribution of age of onset for heroin use.

Additional file 4:.

Supplementary file 4: Sequence data.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, Y., Zhu, Y., Lai, J. et al. Association studies of dopamine synthesis and metabolism genes with multiple phenotypes of heroin dependence. BMC Med Genet 21, 157 (2020). https://doi.org/10.1186/s12881-020-01092-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12881-020-01092-0

Keywords