Association of PICK1 and BDNF variations with increased risk of methamphetamine dependence among Iranian population: a case–control study

Abstract

Background

Genetic factors play an important role in susceptibility to methamphetamine dependency. In this line, protein that interact with C-kinase-1 (PICK1) and brain-derived neurotrophic factor (BDNF) genes are linked to methamphetamine dependence (substance use disorder). Thus, in a case–control study, we investigated the association between polymorphisms of PICK1 and BDNF genes and methamphetamine dependence in an Iranian population.

Methods

Total of 235 cases and 204 controls were recruited in a period between 2015 to 2018. The PICK1-rs713729, -rs2076369 and BDNF-rs6265 genotypes were determined via ARMS-PCR assay. Statistical analysis was performed, using SPSS 20.0, PHASE 2.1.1 program as well as SNP Analyzer 2.0.

Results

In the present study, two polymorphisms including PICK1-rs713729 (OR 1.38 (CI 1.08–1.52; P-_value 0.004) in multiplicative and dominant models, and PICK1-rs2076369 (OR 1.31 (CI 1.10–1.56; P-_value 0.002) in multiplicative, dominant and co-dominant models were associated with the risk of methamphetamine abuse. Moreover, haplotype analysis showed a significant association of haplotype AG (OR 2.50 (CI 1.50–4.16; P-_value 0.0002) in dominant, recessive and co-dominant models, and haplotype TT (OR 0.67 (CI 0.50–0.91; P-_value 0.009) in dominant and co-dominant models with the risk of methamphetamine abuse. None of the polymorphisms in this study had a high level of linkage disequilibrium.

Conclusion

Our findings indicate that the PICK1 gene polymorphism might affect the risk of methamphetamine dependency in our population.

Background

Methamphetamine is the most extensively abused illegal drug and this is a growing global problem [1,2,3]. Many changes occur in the human brain as a result of methamphetamine abuse [4]. Methamphetamine leads to the release of synaptic dopamine, which might be a reason for the increased addiction [4]. In this context, methamphetamine can have neurotoxic effects on dopaminergic neurons and also contribute to changes in another neurotransmitter system, particularly glutamate [5, 6]. Glutamate-mediated excitotoxicity is the major mechanism by which methamphetamine damages the central nervous system (CNS) [7, 8]. Moreover, the primary target of methamphetamine is the dopamine transporters, which remove dopamine from the extracellular space at the synapse and controls dopamine signals [9]. The function and surface availability of the dopamine transporter are regulated via different cellular mechanisms [10, 11]. Additionally, it was stated that the density of dopamine transporter is low in the caudate/putamen of methamphetamine abusers, suggesting that the long-term use of methamphetamine results in damage to dopaminergic neurons [12, 13].

Several factors are involved in methamphetamine abuse. In this line, genetic factors play an important role in susceptibility to the use of methamphetamine [1,2,3]. The dopamine transporter, polymorphisms have been shown to be a risk factor for prognosis of prolonged-type methamphetamine psychosis, such as a single nucleotide polymorphism (SNP) [14, 15]. Moreover, in a systematic review, Bousman et al. showed that several genes are involved in methamphetamine-related disorders. It has been shown that three genes including catechol-O-methyltransferase (COMT), gamma-aminobutyric acid type A receptor subunit alpha1 (GABRA1) and dopamine receptor D4 (DRD4), nine genes including, brain-derived neurotrophic factor (BDNF), arrestin beta 2 (ARRB2), cytochrome P450 2D6 (CYP2D6), glutathione S-transferase mu 1 (GSTM1), glycine transporter-1 (GLYT1), glutathione S-transferase P1 (GSTP1), solute carrier family 22 member 3 (SLC22A3), prodynorphin (PDYN), and protein interacting with C kinase (PICK1), two genes including gamma-aminobutyric acid receptor subunit gamma-2 (GABRG2) and v-akt murine thymoma viral oncogene homolog 1 (AKT1) seem to be associated with methamphetamine abuse or dependence in Japanese and Han-Chinese populations [15]. It is indicated that the genetic epidemiology of methamphetamine abuse is very complex. There is an association between PICK1 gene and drug addiction, such as methamphetamine abuse [16]. In this respect, the PICK1 gene is mapped to chromosome 22q13.1 [17]. Interestingly, PICK1 interacts with dopamine transporter, which leads to the clustering of dopamine transporter on the cell surface and subsequently the improvement of dopamine transporter uptake activity [16, 18].

Furthermore, BDNF-rs6265 is a functional SNP, which affects drug dependency [19]. BDNF gene is mapped to chromosome 11p. It was found that BDNF-rs6265 (Val66Met) SNP is linked to susceptibility to methamphetamine dependence in a Thai and Malaysian population [20,21,22]. In this respect, BDNF proteins are involved in the regulation of synaptic transmission as well as the process underlying substance use disorder (SUD) [23]. BDNF was shown to support the survival and protection of dopaminergic neurons following methamphetamine administration in mice [6]. The functional effects of SNPs, which are linked to SUD, are often unclear; hence, more investigations are required on different populations to define how these variants influence gene expression and function. In this case, worldwide participatory attempts are warranted to promote the accessibility of large population-based datasets/samples and to increase the feasibility of genetic association studies [15]. On the other hand, a systematic review by Alam-MehrJerdi et al. indicated that the state of methamphetamine abuse requires further research on the epidemiology and health-related implications in Iran [24]. Therefore, we aimed to investigate the association between PICK1 and BDNF main SNPs and methamphetamine dependence (substance use disorder (SUD)) in a cohort of Iranian individuals.

Methods

The aim, design and setting of the study

This case–control study involved 235 cases (with methamphetamine dependence (SUD)) and 204 gender-matched controls (healthy individuals). In our study, SUD was recognized by the 11 criteria. The 11 criteria are divided into four categories of behavior, such as impaired control, social impairment, risky use and pharmacological indicators (tolerance and withdrawal) related to the substance use (Table 1) [25, 26].

Table 1 Criteria for substance use disorders (SUD)

The PICK1-rs713729, PICK1-rs2076369 and BDNF-rs6265 genotypes were analyzed by amplification refractory mutation system-polymerase chain reaction (ARMS-PCR) assay. Subsequently, statistical analysis was performed, using SPSS 20.0 (IBM Inc., Chicago, IL, USA), PHASE 2.1.1 program as well as SNP Analyzer 2.0.

The characteristics of participants

Study population and clinical data

This study was performed in accordance with the Declaration of Helsinki (1964) and its subsequent amendments. Moreover, approval was obtained from the local Ethics committee of Mashhad University of Medical Sciences (IR.MUMS.fm.REC.1394.421), 439 blood samples were collected from 204 controls and 235 cases from Mashhad, Iran from 2015 to 2018. After explaining the study objectives, a written informed consent was obtained from all participants. A questionnaire was used to collect demographic and other essential information from all participants (Table 2). The selection procedure included confirmed urine test (addiction test) and also the availability of complete patient’s follow-up data. Moreover, healthy participants, individually matched on age, were recruited from the Health Examination Centre who were receiving routine medical examinations.

Table 2 Demographic and clinical characteristics of controls and methamphetamine dependence (SUD)

Description of materials

Blood collection and DNA extraction

In order to extract the DNA, approximately 10 millilitres (ml) of peripheral blood was obtained from all participants and immediately subdivided into tubes containing sterile ethylene diamine tetra acetic acid (EDTA) [27]. DNA extraction from whole blood was extracted using salting-out technique. Then, the extracted DNA was quantified by the ratio of absorbance at 260 nm (nm) and 280 nm (A_260/280) via BioTek™ Epoch™ Microplate Spectrophotometer (Winooski, VT, USA,) as well as via gel electrophoresis and finally stored at -20 °C until used.

Target single nucleotide polymorphisms (SNPs) determinations (Marker selection)

In this study, the SNPs were selected using available SNPs databases and published articles. Such articles examined intron and exon SNPs, which might alter the affinity of PICK1-rs713729, PICK1-rs2076369 and BDNF-rs6265 to methamphetamine dependence (SUD) (Table 3). Moreover, potential functional SNPs were included in order to meet the following criteria: minor allele frequency (MAF) > 0.05 (5%), heterozygosity > 0.15 (15%) and also validated SNPs in articles and databases. Furthermore, in order to inhibit redundancy in SNPs genotyping, SNPs that are not located in strong linkage disequilibrium (LD) were chosen.

Table 3 Characteristics of the investigated polymorphisms in this study

Genotyping

To determine the genotype frequency of PICK1-rs713729, PICK1-rs2076369 and BDNF-rs6265 an ARMS-PCR method was used. Specific primers for PCR amplification were designed via web tools, such as Primer1 and also WASP (web-based allele-specific primer designing tool) [28].

PCR amplifications for PICK1-rs713729, PICK1-rs2076369 and BDNF-rs6265 were conducted in a 10–15 µl (μl) volume per reaction, containing 3 µl Taq 2 × master mix (Ampliqon, Germany), 10 µM of each primer and 100 nanogram (ng) DNA. Moreover, the specific primers used to detect PICK1-rs713729, PICK1-rs2076369 and BDNF-rs6265 SNPs are listed in Table 4. For PICK1-rs2076369, we also used betaine (Ampliqon, Germany) as an enhancer in PCR.

Table 4 Primer sequences used for genotyping in ARMS-PCR

The ARMS-PCRs condition for each primer is as follows, Table 5. In general, initial denaturation at temperature 94 °C for five minutes, then 35 cycles including denaturation at 94 °C for 25 s, annealing at alternative °C for 25 s (based on each primer), an elongation at 72 °C for 30 s followed by 72 °C for seven minutes as the final elongation step (Table 5).

Table 5 The ARMS-PCRs condition for targeted SNPs

Absence or presence of mutant or normal PCR products were detected via gel electrophoresis in 3% agarose gel by ultraviolet (UV) trans illuminator (Gel Doc; U:Genius).

Statistical analysis

A Hardy–Weinberg equilibrium (HWE) method was used to evaluate the differences in data for statistical significance. HWE assumption was investigated by the Pearson χ² distribution with 1 degree of freedom. Allele and genotype frequencies were calculated, and the differences between groups were evaluated by Chi-squared tests. Then, the association between methamphetamine, risk factors and alleles/genotypes was evaluated by binary logistic regression, estimating Odds ratios (ORs) and also 95% confidence intervals (CIs). Three logistic regression models were used to analyse the SNPs, using different genetic models (additive, dominant, and recessive). For the analysis of SNP-SNP interactions, an adjusted logistic regression model was used to estimate the multiplicative interaction effect of the SNPs, located on the same haplotype. P-_value = 0 < 0.05 was considered to be statistically significant. SPSS 20.0 (Inc., Chicago, IL, USA), PHASE program as well as SNP Analyser 2 software were used for further statistical analysis [29].

Haplotype analysis

Haplotypes were generated and assembled from the genotyped data by PHASE program, to reconstruct haplotypes, and SNP Analyzer 2 software [29, 30]. In the present study, P-_values of less than 0.05 were considered to be statistically significant. Moreover, Bonferroni correction was also used to account for multiple testing; thus, a two-tailed P-_value < 0.016 (₌0.05/3 SNPs) was considered to be statistically significant in the present study.

Results

Identification of single nucleotide polymorphisms (SNPs) and association studies

There were no significant associations between BDNF-rs6265 and the risk of methamphetamine dependence (SUD). On the contrary, a significant association was observed between two SNPs and the risk of methamphetamine abuse including PICK1-rs713729 (OR 1.38 (CI 1.08–1.52; P-_value 0.004) in multiplicative and dominant models, and PICK1-rs2076369 (OR 1.31 (CI 1.10–1.56; P-_value 0.002) in multiplicative, dominant and co-dominant models (Table 6). Moreover, haplotype analysis showed that specific haplotypes related to these SNPs were associated with methamphetamine dependence (SUD). In this line, analysis of PICK1-rs713729 and PICK1-rs2076369 haplotypes in our population showed that the haplotype AG (OR 2.50 (CI 1.50–4.16; P-_value 0.0002) in dominant, recessive and co-dominant models and haplotype TT (OR 0.67 (CI 0.50–0.91; P-_value 0.009) in dominant model and co-dominant model had a significant association with the risk of methamphetamine dependence SUD) (Table 7).

Table 6 Analysis based on different genetic models

Table 7 Haplotype analysis based on different genetic models

Distribution of single nucleotide polymorphisms (SNPs)

The frequencies of the genotypes with high-quality genotype call were as follows: PICK1-rs713729 “TT” 174 (77.67%), “AT” 33 (14.73%), and “AA” 17 (7.58%) among cases and “TT” 168 (87.95%), “AT” 17 (8.90%) and “AA” 6 (3.14%) among controls; PICK1-rs2076369 “GG” 117 (52.23%), “GT” 80 (35.71%) and “TT” 27 (12.05%) among cases and 70 (36.08%), 100 (51.54%), 24 (12.37%) and among controls, respectively. Moreover, BDNF-rs6265 genotype frequencies were: “CC” 161 (72.52%), “CT” 55 (24.77%) and “TT” 6 (2.70%) among cases and 142 (73.57%), 48 (24.87%) and 3 (1.55%) among controls, respectively. Table 7 shows the conferred risk by each haplotype. In this study, two haplotypes, AG and TT were significantly different between the normal and methamphetamine dependence (SUD) individuals. The prevalence of TT haplotype in the case group (27%) was lower than that of the normal individuals (34%), and the frequency of the AG haplotype in the methamphetamine dependence (SUD) individuals was higher (12%) than that (3%) of the normal individuals (Additional file 1: Table 1). We also examined whether the three SNPs were in LD. In this line, there was no strong LD between these SNPs (Additional file 2: Table 1).

Discussion

The SUD is a persistent relapsing disorder with harmful consequences [31]. In this respect, the brain reward system such as mesocortical dopamine system is the common feature of this disorder [32, 33]. Furthermore, abnormal dopamine and glutamate systems are associated with the pathophysiology of SUD and dependency. However, the development of drug dependency is affected via factors, such as pharmacological effects on mental status, environmental and individual factors, such as genetics. In this respect, genetic factors are suggested to have a greater effect on drug dependency. It was shown that drugs abuse and addiction cause dopamine and opioid peptides to be released into the ventral striatum, which causes “high” sensation in abusers [33, 34].

The use of methamphetamine has become a serious health concern in Iranian methadone dependent patients [24, 35]. In Iran, it is used to improve sexual performance and promote physical energy due to its stimulating effects [36]. It is shown that methamphetamine abuse has surged from 3.9% amongst both genders in 2007 to 60.3% for men in 2014 and 89.5% for women in 2015–2016. Remarkably, the frequency of methamphetamine dependence amongst female methadone dependent patients in the Iranian population was higher. This dependency has been linked to multiple health issues in the social and health contexts of both genders, particularly women [35].

Showing a possible association between PICK1-rs713729 and PICK1-rs2076369 in the PICK1 gene promoter with methamphetamine dependence were the key findings of the present study. Consistent with our study, Matsuzawa et al. showed that PICK1-rs713729 and PICK1-rs2076369 were significantly associated with methamphetamine abuse in a Japanese population (Table 8). Additionally, they revealed that PICK1-rs713729 was linked to those with spontaneous relapse of psychosis [16]. Moreover, PICK1 gene was reported as the beginning of methamphetamine addiction, worse prognosis as well as spontaneous relapse [37]. PICK1 is involved in the targeting and localization of synaptic membranes proteins and also in the surface of dopamine transporter clustering on the cell surface [18, 38]. This clustering leads to an increase of dopamine transporter uptake activity. Thus, PICK1 expression levels can alter the role of the dopamine transporter and can affect the pathogenesis of methamphetamine abuse/dependency [16]. Moreover, the interactions of PICK1 with α-amino-3-hydroxy-5-methyl-4-isoxazole propionate (AMPA) receptors, metabotropic glutamate receptors and D-serine synthesizing enzyme have been shown and might be implicated in the pathogenesis of drug-related disorders [39,40,41,42]. Further studies are still required to show the direct association of these proteins in the pathogenesis of methamphetamine abuse.

Table 8 The association studies regarding methamphetamine and BDNF/PICK1 gene polymorphisms in the methamphetamine abusers and controls

Based on previous studies, it was reported that gene variations, which are associated with glutamatergic and serotonin systems explain differences in SUD and/or dependency risk between individuals. It is indicated that glutamate-related genes influence the risk of SUD and/or dependency. The glutamate receptor genes were suggested to interact with BDNF by BDNF- tropomyosin kinase B (TrKB) transduction signaling cascade. BDNF is a neurotrophic factor, which is involved in the expansion, maintenance and survival of dopaminergic neurons in CNS [43].

In the present study, no association was found between BDNF-rs6265 and the risk of methamphetamine dependence (SUD) in our population. This finding is consistent with the result of methamphetamine dependent male Caucasian individuals, which showed no association in 193 non-psychotic males (117 methamphetamine-dependent cases and 76 controls) [44]. However, Cheng et al. showed that there was an association between methamphetamine-dependency and BDNF gene in 103 methamphetamine abusers and 122 normal controls. They noted that the lower 66Met carriers were linked to substance abuse [21]. Furthermore, Sim et al. described an increase in a Chinese subgroup of Malaysian methamphetamine-dependent subjects (n = 24), which was not found among other Malaysian ethnic groups [22]. The various results for this SNP might be due to different sample sizes and also different ethnicities as well as genetic diversity (Table 8). It is indicated that ethnic differences can affect the frequency of BDNF Val66Met [45, 46].

Long-term methamphetamine-induced brain changes are significantly dependent on BDNF genetic variation [47]. Val66Met might enhance the risk of suicide behavior [45]. Moreover, BDNF variations might be involved in methamphetamine withdrawal. BDNF levels, equal or less than 1,251.0 pg/ml has been stated to be linked to depression symptoms during methamphetamine withdrawal [48]. BDNF Val66Met had an important effect on the Treatment Effectiveness Score (TES), methamphetamine-negative urine drug screens during treatment, which was higher among Val/Val Caucasian carriers [49]. Interestingly, there were more pretreatment days with methamphetamine use in females with Val/Val genotype than males [10]. Estrogen improves BDNF expression, which has implications on the release of striatal dopamine caused by methamphetamine and protects against neurotoxicity caused by methamphetamine [50, 51]. Thus, females can be able to use methamphetamine more frequently [10].

Furthermore, there was no association between PICK1 and BDNF gene polymorphisms in our study. That may be attributed to the different positions of the gene polymorphisms on different chromosomes [PICK1 (22q13.1) and BDNF (11p13–15)].

There were some limitations in our study. We did not investigate all the SNPs in these genes. Although, the sample size was large enough to detect an association between these SNPs and SUD in comparison with previous studies, larger sample size would strengthen the results. We therefore recommend further replication studies with larger sample size in order to validate and explain this association with conclusive findings.

Conclusions

Collectively, the variation in the PICK1 gene was associated with methamphetamine dependence (SUD), reflecting the underlying biological mechanisms, which can make a bridge between pathways and methamphetamine dependence (SUD). Our findings suggest that the PICK1 gene might be involved in susceptibility to SUD in our population. These findings can be helpful in rehabilitation programs and psycho-education for those who have substance dependency. In this case, potential genetic predictors can be used for individuals susceptible to SUD. All in all, understanding genetic variations might help to understand the biological mechanisms of progression, and suppression of methamphetamine. In addition, our findings provide the basis for future genetic research on the use of methamphetamine dependency and related neurological side effects.