https://doi.org/10.22037/ijpr.2019.15097.12874

MDMA Abuse in Relation to MicroRNA Variation in Human Brain Ventral Tegmental Area and Nucleus Accumbens

authors:

avatar Goksun Demirel a , b , 1 , * , avatar Esra Guzel c , 1 , avatar Chad J Creighton d , e , avatar Yeter Erol Ozturk f , avatar Ceyhun Kucuk g , avatar Hızır Asliyuksek g , avatar Türkan Yurdun h

Department of Pharmaceutical Toxicology, Faculty of Pharmacy, Biruni University, Istanbul, Turkey.
Department of Pharmaceutical Toxicology, Faculty of Pharmacy, Çukurova University, Adana, Turkey.
Department of Molecular Biology and Genetics, Institute of Health Sciences, University of Health Sciences, Istanbul, Turkey.
Department of Medicine, Baylor College of Medicine, Houston, Texas.
Dan L. Duncan Cancer Center Division of Biostatistics, Houston, Texas.
Chemistry Department, Council of Forensic Medicine, Istanbul, Turkey.
Council of Forensic Medicine, Istanbul, Turkey.
Department of Pharmaceutical Toxicology, Faculty of Pharmacy, Marmara University, Istanbul, Turkey.
G D and E G. These authors contributed equally to this work.
Iranian Journal of Pharmaceutical Research : IJPR: Vol.18, issue 4; 1989-1999
article type: Original Article
received: October , 2018
accepted: April , 2019

How To Cite Demirel G, Guzel E, Creighton C J, Ozturk Y E, Kucuk C, et al. MDMA Abuse in Relation to MicroRNA Variation in Human Brain Ventral Tegmental Area and Nucleus Accumbens. Iran J Pharm Res. 2019;18(4):e124705. https://doi.org/10.22037/ijpr.2019.15097.12874.

Abstract

3,4-methylenedioxymethamphetamine (MDMA) is one of the most widespread illegal drugs, that have been used particularly by young people in the 15-34 age group. MicroRNAs (miRNAs) are endogenously synthesized, non-coding, and small RNAs that post-transcriptionally regulate their target genes’ expression by inhibiting protein translation or degradation. miRNAs are increasingly implicated in drug-related gene expressions and functions. Notably, there are no reports of miRNA variation in the human brain in MDMA abuse. We here present a miRNA profiling study – the first such study, to the best of our knowledge – into the post-mortem human brains of a sample of people with MDMA abuse, along with non-drug dependent controls. The miRNA profiling of nucleus accumbens (NAc) and ventral tegmental areas (VTA) was performed by microarray analysis. Subsequently, two candidate miRNA putative biomarkers were selected according to significant regional differential expression (miR-1202 and miR-7975), using quantitative reverse-transcription PCR (qRT-PCR). We showed that the expression level of miR-7975 was significantly lower in the VTA regions of the 30 MDMA users, as compared with the 30 control samples. Another significantly deregulated miR-1202 was down-regulated in the NAc regions of 30 MDMA samples in comparison to the control samples. Alteration of these miRNAs can potentially serve as novel biomarkers for MDMA abuse, and warrant further research in independent and larger samples of patients.

Introduction

MDM(3,4methylenedioxymethamphetamine), commonly known as ‘ecstasy’ or ‘molly’, is an illicit, synthetic drug which causes serotonergic neurotoxicity in the brain (1). Prosocial behavior and empathogenic feelings are boosted by MDMA’s effect on emotional and cognitive empathy processes. MDMA is considered to have therapeutic effects in psychotherapy; however, high-dose MDMA use has been associated with multiple complications such as persisting cognitive deficits and panic attacks (2-4), and can also lead to drug addiction, toxicity, or death (5, 6). MDMA, the most frequently-used amphetamine-type stimulant, causes the release of serotonin, dopamine, and norepinephrine in the parts of the brain that are targeted by other addictive drugs (7). Addictive drugs target the mesocorticolimbic dopamine (DA) system in the addiction center of the brain, including the ventral tegmental area (VTA), prefrontal cortex, and nucleus accumbens (NAc) (8, 9). The VTA, which contains the largest group of dopamine neurons, plays an important role in reward-related and goal-directed behaviors such as cognitive and emotional processes (10-12). The NAc receives dopaminergic inputs from the VTA and has been associated with control of behavior by drug-related stimuli and with relapse. (13). Hence, the VTA and NAc are central components of the brain’s reward system, and also play a major role in the response to drugs abuse (14, 15).

The latest data suggests that MDMA has become one of the most popularly used illegal drugs in Europe (9.3 million males and 4.7 million females). The European Drug Report has estimated that the availability and use of high-dose MDMA tablets is increasing, and this trend is set to continue. The same report emphasized that the number of people dying from MDMA use has significantly increased in recent years (16). Therefore, it is necessary to understand the molecular mechanism underlying the progression of MDMA abuse in order to find specific biomarkers for therapeutic approaches. There is an increasing body of evidence indicating that microRNAs (miRNAs) are associated with the regulation of drug-induced pathways, including the improvement of drug-addiction behaviors (17).

miRNAs are approximately 18-22 bp long post-translational regulators that target mRNAs to repress post-transcriptional gene expression by binding to complementary sequences in the 3’UTR (18, 19). Studies have demonstrated that miRNAs play a significant role in the activity of the brain’s regions in relation to drug addiction, including the NAc, VTA and other regions of the brain’s reward pathways (20).

Many studies have demonstrated that miRNA’s expression levels are altered by drugs such as heroin (21), cocaine (22), nicotine (23), and opioids (24). However, miRNA expression profiles in postmortem human MDMA users have been never shown before. Therefore, it is of importance to study miRNA profiling specifically in the brain regions associated with MDMA seeking.

We aimed to determine the miRNA expression profiles of the NAc and VTA regions in MDMA users’ postmortem brain tissues and matched controls in this study. We found miR-1202 and miR-7975, which are significantly down-regulated in the brain samples with MDMA, can be used as candidate biomarkers for treatment approaches to MDMA-seeking behavior.

Experimental

Collection of Postmortem Brain Samples

This study has been reviewed and approved by an institutional review board of Marmara University Faculty of Medicine (IRB No: 09.2016.376). Postmortem human brain tissues consisting of the nucleus accumbens and ventral tegmental areas of 30 positive MDMA-users and 30 control samples were collected from the Morgue Department of the Council of Forensic Medicine in Istanbul, Turkey. Informed consent was obtained from family members of the postmortem individuals prior to autopsy, in accordance with the standards established by the Council of Forensic Medicine.

Blood tests were used to detect the presence of MDMA in those subjects included in this study. The drug-free control group was selected from postmortem individuals who had a negative result for all drugs, and no history of substance use prior to death. In both groups, the following cases were excluded: significant head injury, the presence of an infectious disease (e.g HIV/AIDS, hepatitis B or C), chronic disease, or a history of mental disorders. The brains were removed during autopsy, and tissue samples were immediately collected from the VTA and the NAc, then were flash-frozen in liquid nitrogen. The frozen tissues were individually kept in microcentrifuge tubes at -80 °C until RNA extraction.

Screening of Blood Samples in Postmortem Individuals

All blood samples were screened by general unknown screening. In this general method, the analysis of drugs was accomplished by liquid chromatography, coupled with tandem mass spectrometry (LC– MS/MS) and gas chromatography-mass spectrometry (GC–MS) techniques. GC– MS analyses were carried out using a Shimadzu 2010 QP GC–MS instrument and a Restek RTX-5-MS (30m 0.25mm 0.25mm; 5% phenylmethylsiloxane) column was used for chromatographic separations. The mass detector was operated in electronic impact mode at 70 eV and the mass spectra were scanned at m/z 60–550 amu. The injection volume was 1 µL in splitless mode. LC–MS/MS analyses were carried out using an Agilent 6460 triple quadruple tandem mass spectrometer instrument and Agilent Poroshell C18 (dimensions 150 mm 2.7 mm 4.6 mm) column for chromatographic separations. Tandem mass spectrometry (MS/MS) detection was performed using electrospray positive and negative ionization in MRM mode, with an injection volume of 10 µL. Many reports have shown that a toxic or even fatal dose entails blood levels ranging from 0.5 to 10 mg/L in most cases (25).

Total RNA Extraction

NAc and VTA brain tissues were blended in 1 mL TRIzol reagent (Invitrogen, USA) per 50-100 mg of tissue. Total RNA was extracted as per the manufacturer’s instructions. Total RNA concentrations and purities of the samples were quantified by spectrophotometry using a NanoDrop ND-2000 system (Thermo Fisher Scientific, Inc., Wilmington, DE). The RNA Integrity Number (RIN) was measured by the Bioanalyzer 2100 (Agilent Technologies, USA).

miRNA Microarray and Data Analysis

miRNA expression profiling of the NAc and VTA regions of eight postmortem brain samples (four of them from MDMA users and four of them from the control group) were performed using Agilent 8X60K Human miRNA V21.0. An Agilent miRNA microarray slide contains 2,549 probes, including control probes, according to the latest version of Sanger miRBase (Release 21.0). All of the related reagents were purchased from Agilent Inc. 50 ng of total RNA from four MDMA users and four control brain tissues were labeled with Cy3, using the “Agilent miRNA Labeling and Hybridisation Kit”, following the protocol recommended by the manufacturers. Each labeled RNA sample was hybridized onto an Agilent 8X60K Human miRNA microarray V21 chip and incubated at 55 °C for 20 h. The chip was washed after hybridization and immediately scanned by an Agilent Microarray Scanner with Surescan High Resolution Technology (Agilent Technologies, Santa Clara, CA). Feature Extraction v10.7.3.1 (Agilent Technologies) software was used to analyze the scanned images. Bioconductor software, a normalization tool for microarray data, was used on the raw data via R version 2.15 quantile normalization. Log transformed data was used to calculate fold change and P value for each sample. Potential candidate miRNAs with p < 0.01 and 1.4-fold change were selected for further validation by quantitative reverse transcription-polymerase chain reaction (qRT-PCR). The microarray data were then deposited into the NCBI Gene Expression Omnibus (GEO, Accession number GSE104315).

cDNA Synthesis and qRT-PCR

To validate the differential expression of miR-1202 and miR-7975, RNA extracted from the 30 MDMA users and 30 control samples (including the eight RNA samples used in the microarray) were studied using qRT-PCR. Complementary DNA (cDNA) synthesis was studied by using a “TaqMan MicroRNA Reverse Transcription Kit”, following the manufacturer’s protocol (Applied Biosystems, Foster City, CA). For each selected miRNA, TaqMan assays were purchased from Applied Biosystems, Foster City, CA. small-nucleolar RNA (snoRNA) RNU43 was used for normalization of miRNA expression analyses. (mir7975 5’ AUCCUAGUCACGGCACCA-3’, mir1202 5’- GUGCCAGCUGCAGUGGGGGAG-3).

Quantitative RT-PCR was performed using Universal Master Mix (Applied Biosystems) in a Roche LightCycler 480-II real-time thermal cycler (Roche, Switzerland, Basel). The reactions were applied in the following conditions: one cycle of 95 °C for 5 min followed by 40 cycles of 95 °C for 10 s, 60 °C for 20 s, and 72 °C for 25 s. Each was performed in triplicate. The relative quantification analysis was performed by the delta-delta-Ct method (2-ΔΔCT) as described previously (26).

Statistical analysis

Data were plotted as mean ± standard error and two-sided Student’s tests were used for statistical analysis of the qRT-PCR. A two-sided p-value of 0.05 or below was considered statistically significant. Receiver operating characteristic (ROC) curves were plotted using SPSS (Statistical Package for the Social Sciences) to see the power of control and validated miRNAs to differentiate from the MDMA samples. For ROC analysis, logistic regression was conducted and the predicted probabilities were calculated for each miRNA for the VTA and NAc regions of brain. The area under the curve (AUC) was calculated with a 95% confidence interval and accepted as significantly different from 0.5 when the p-value was greater than 0.5, meaning that the logistic regression classifies the group significantly better than by chance.

Results

In this study, we conducted miRNA profiling of the NAc and VTA regions of subjects who had used MDMA and control subjects. A total of 120 postmortem human brain tissues (30 from VTA, 30 from NAc and their matched controls) were obtained from the Istanbul Morgue Department of the Council of Forensic Medicine (Turkey). All of the included cases had been part of the Turkish population. The samples were screened for the presence of MDMA by the Istanbul Toxicology Department of the Council of Forensic Medicine (Turkey). Clinical and toxicology information about the postmortem MDMA and control individuals are provided in Table 1.

In every case, toxicology reports were reviewed and the absence or presence of MDMA was documented. MDMA concentrations ranged from 197.6 to 2017.8 ng/mL, with a median range of 554.75 ng/mL in postmortem blood samples. In the 30 MDMA positive cases, no other drugs of abuse were detected. The mean age of the MDMA user group was 27.96 ± 9.45, whereas the control group had an average age of 29.76 ± 10.57. Brain pH level changed from 5.56 to 6.98 in the MDMA group and 6.1 to 6.93 in the control group. No significant difference was found in PMI and brain pH levels between the two groups. Likewise, there was no difference between the RIN values of RNA extracted from the VTA and NAc regions of the MDMA users and the control group. Information on the subjects, including gender, age, PMI, RIN, and brain pH levels are summarized in Table 2.

Comparison of the miRNA expression profiles of the NAc and VTA brain regions of the postmortem human brain tissues from four of the MDMA users and four of the control samples was carried out using microarray analysis. According to the quantile normalized data, we found 41 significant deregulated miRNAs in the VTA regions, and 46 significant deregulated miRNAS in the NAc regions of both the MDMA-using group and the control group, with fold >1.4 between the two groups in both regions. A heat map illustration of differentially expressed miRNAs is presented in Figure 1.

Among the significantly differentially regulated miRNAs, miR-7975 was found to be down-regulated in only the VTA regions of the MDMA samples and miR-1202 to be down-regulated in only the NAc regions of the MDMA samples. Hence, miR-7975 and miR-1202 were selected, following a detailed literature search, for further qRT-PCR validation experiments. The qRT-PCR results confirmed that the expression of miR-7975 was significantly lower in the VTA region of the MDMA users than in that of the control samples. miR-1202 was found to have significantly reduced expression in the NAc region of the MDMA users compared to the control samples (Figure 2).

To test the power of miR-1202 and miR-7975 for distinguishing MDMA specimens from control samples, ROC were plotted. The values showed that miR-1202 (MDMA NAc vs control NAc) and miR-7975 (MDMA VTA vs control VTA) had an area under the curve (AUC) of 0.822, and 0.778, respectively (Figure 3). The results demonstrate their sufficiency to distinguish MDMA specimens from control specimens.

We also used three different miRNA target prediction databases (TargetScan 7.0 (27), mirDB (28) and miRWalk 2.0 (29)) to identify putative target genes and how they are affected by miR-1202 and miR-7975. Our bioinformatics studies have demonstrated that among the putative targets of mir-1202, GRM4 (30), GLP1 (31), CTNN1 (32), GRIN1, GRIN2A, GRIN2B, SLC38, and GRIK3 (33) were considerable candidates because their functions are associated with drug addiction or glutamatergic neurotransmission. In addition, human miRNA expression databases show that miR-1202 is highly expressed in the brain compared to other tissues (34, 35). This result indicates that miR-1202 may play an important regulatory role in its target genes in drug addiction. The analysis revealed that potential targets of miR-7975 – PDE5A (36), SLC6A6, SLC6A9 (37), CALN1 (38), GDNF (39), TMOD2, CHRNB2 (40) – play a role in several drugs of abuse which affect neurotransmitter processes.

Table 1

Clinical information for the 60 postmortem brains included in this study

Clinical DiagnosisAgeGenderCause of DeathMDMA in Blood [ng/mL]PMI[hours]pH
1Control31MRoad traffic accidentN/A176.13
2Control49MCardiac arrestN/A166.46
3Control38FAcute coronary thrombosisN/A206.32
4Control32FWork accidentN/A176.12
5Control32MWork accidentN/A146.67
6Control25MCardiac arrestN/A186.84
7Control23MGunshot woundsN/A236.9
8Control18MRoad traffic accidentN/A196.41
9Control17FRoad traffic accidentN/A96.65
10Control19MRoad traffic accidentN/A236.48
11Control48MGunshot woundN/A186.59
12Control38MRoad traffic accidentN/A246.64
13Control52FAtherosclerotic heart diseaseN/A196.4
14Control19MRoad traffic accidentN/A86.7
15Control21MRoad traffic accidentN/A156.54
16Control24MRoad traffic accidentN/A196.43
17Control26MGunshot woundN/A236.1
18Control30FRoad traffic accidentN/A136.14
19Control37MHypertensive heart diseaseN/A136.32
20Control18MRoad traffic accidentN/A106.8
21Control35FHypertensive heart diseaseN/A176.46
22Control41MCardiac arrestN/A126.93
23Control28FRoad traffic accidentN/A86.43
24Control15MRoad traffic accidentN/A226.41
25Control36FGunshot woundsN/A146.52
26Control32MWork accidentN/A226.36
27Control47MCardiac arrestN/A206.84
28Control18FWork accidentN/A186.48
29Control24FGunshot woundsN/A156.76
30Control20FWork accidentN/A116.85
31MDMA user32MDeath due to MDMA intoxication1523.6236.64
32MDMA user28MDeath due to MDMA intoxication606.8156.88
33MDMA user33MDeath due to MDMA intoxication1119.2126.66
34MDMA user19MRoad traffic accident318.696.36
35MDMA user18FHyperpyrexia402.3156.66
36MDMA user32FGunshot wounds207.1126.48
37MDMA user47MRoad traffic accident276.9226.68
38MDMA user22FGunshot wound401.3186.98
39MDMA user30MDeath due to MDMA intoxication1023.6126.87
40MDMA user31FDeath due to MDMA intoxication756.3136.63
41MDMA user46MDeath due to MDMA intoxication2017.8156.55
42MDMA user16MGunshot wounds269.3226.7
43MDMA user20FRoad traffic accident252.8216.8
44MDMA user17FFall from high197.6176.94
45MDMA user35FDeath due to MDMA intoxication569.3236.4
46MDMA user36MGunshot wounds262.7145.92
47MDMA user20FDeath due to MDMA intoxication669.7196.62
48MDMA user35MDeath due to MDMA intoxication1324.6236.9
49MDMA user27MDeath due to MDMA intoxication540.2246.14
50MDMA user20FDeath due to MDMA intoxication524.6226.42
51MDMA user16MDeath due to MDMA intoxication825.6146.62
52MDMA user47MCardiac arrest450.6146.31
53MDMA user33MDeath due to MDMA intoxication895.3166.7
54MDMA user21FDeath due to MDMA intoxication908.7156.8
55MDMA user17MDeath due to MDMA intoxication654.4126.43
56MDMA user20FRoad traffic accident306.8136.6
57MDMA user24MGunshot wound379.5185.56
58MDMA user42FDeath due to MDMA intoxication648.7246.86
59MDMA user31MRoad traffic accident295.3176.6
60MDMA user24MDeath due to MDMA intoxication852.7196.39
Table 2

The characteristics of MDMA and control samples that are involved in the study

Control (N=30)MDMA (N=30)p value
Sex19 M/11 F18 M/12 F-
Age29.76 ±10.5727.96±9.450.48
Brain pH6.52±0.236.57±0.300.50
RIN (VTA)6.38±0.286.47±0.340.25
RIN (NAc)6.76±0.356.51±0.720.34
PMI17.01±4.3416.57±4.650.64
Heat-map illustration of deregulated miRNAs in VTA with MDMA samples vs VTA control samples and NAc with MDMA vs NAc control postmortem human brain tissues

Heat-map illustration of deregulated miRNAs in VTA with MDMA samples vs VTA control samples and NAc with MDMA vs NAc control postmortem human brain tissues

Figure 1

Relative expression levels of miR-1202 in 30 NAc with MDMA samples compared to 30 NAc control samples (p = 0.008), miR-7975 in 30 VTA with MDMA samples compared to 30 VTA control samples (p = 0.01).

Relative expression levels of miR-1202 in 30 NAc with MDMA samples compared to 30 NAc control samples (p = 0.008), miR-7975 in 30 VTA with MDMA samples compared to 30 VTA control samples (p = 0.01).

Figure 2

Discussion

MDMA use has increased dramatically among adolescents despite growing evidence of its potential harmful effects (16). Although various therapies have been developed to treat MDMA abusers, most of the treatment options available have proved insufficiently effective. Thus, there is an urgent need to identify the underlying mechanisms of MDMA-seeking behavior and to develop effective therapeutic approaches to MDMA abuse.

miRNAs are among the most promising potential therapeutic targets for MDMA because of their ability to be ‘key organizers’’ of gene regulation and expression. There are now several reports that show evidence of the correlation among drug addiction and miRNA regulation (22, 41).

Latest studies have defined changes in aberrant miRNA expression profiles detected with drug consumption, along with drug induced neuroplasticity in the VTA and NAc, as well as other regions of the mesocorticolimbic DA system (21, 42, 43). Although a number of comprehensive analyses of several drug abuse miRNA expressions in animal models and human postmortem brains have been published over the years, profiling studies of differentially expressed miRNAs have not been conducted in connection with human MDMA seeking until now.

In this experiment, we researched the role of miRNAs in searching for therapeutic targets in MDMA users by examining post-mortem human brain tissues using microarray profiling of miRNA expression levels.

Acting on the miRNA microarray levels results, we elective miR-1202 and miR-7975 the best highly dysregulated miRNAs that were altered expression in postmortem brain tissues. We validated the down-regulation of miR-7975 in 30 VTA regions from MDMA users and the down-regulation of miR-1202 in 30 NAc regions, with MDMA users’ brain tissues compared with matched control postmortem brain tissues.

MiR-1202, is a neuron system-enriched miRNA that has been shown to have a key role in cognitive processes, and act through the modulation of the glutamatergic system (30, 44). Moreover, previous research has demonstrated that a decreased miR-1202 expression level was found to negatively correlate with GRM4, also detected in the postmortem brain samples. GRM4 is expressed through the brain and neuronal system which is involved in glutamatergic, dopaminergic, GABAergic, and serotonergic neurotransmission processes. Moreover, higher GRM4 protein expression has beforehand been reported to be incorporated with major depression and connectivity in brain activity (30).

Other studies have demonstrated that GLP1, also known as EMHT1, targets the miR-1202 gene, and plays a role in memory consolidation and retrieval in maintenance mode (45). In addition, the findings showed that alcohol, amphetamine, cocaine, and nicotine are modulated by GLP1 in brain stimulation and reward processes. Researchers suggest that GLP1 analogues could be used in the treatment of drug dependence. It has been suggested that a down-regulated miR-1202 effect is mediated by the suppression of the target gene GLP1, which fuels drug addiction. (31). According to these results and those of our own labor, miR-1202 may be important as a potential therapy for MDMA abuse. In this study, we further highlighted miR-7975, whose functions have not yet been studied. Ron et al have shown that one target gene of miR-7975 is GDNF, which has a key role as a negative regulator of drug of abuse like cocaine and alcohol (39). Moreover, TMOD2, which is another target of mir-7975, has been previously shown to be involved in neurotransmitter discharge and synaptic plasticity in the adult brain (40). Every single microRNAs is able to regulate several hundred mRNAs; therefore, many other predicted mRNA targets may also subscribe to phenotypes consolidated with miR -1202 and miR-7975 deregulation by MDMA use. Therefore, altering the level of these miRNAs can regulate the differential expression of multiple genes, in replication to the MDMA toxicity, which can result in the molecular conformation that leads to abuse of drug.

Conclusion

To the best of our knowledge, this is the first study in literature of a consistent miRNA dysregulation in MDMA-affected postmortem human brain tissues. Overall, our findings show that miR-1202 and miR-7975 are differentially expressed in MDMA-positive individuals, and these miRNAs can serve as potential biomarkers for predicting MDMA abuse and treatment response. Predicted miRNAs such as mir-1202 and miR-7975, have been linked in the mechanisms of abuse of drug, and forward studies can be highly important to developing preventive strategies and fresh therapeutic interventions for MDMA abuse.

Acknowledgements

References

  • 1.

    White SR, Obradovic T, Imel KM, Wheaton MJ. The effects of methylenedioxymethamphetamine (MDMA, «Ecstasy») on monoaminergic neurotransmission in the central nervous system. Prog. Neurobiol. 1996;49:455-79. [PubMed ID: 8895996].

  • 2.

    Patel A, Moreland T, Haq F, Siddiqui F, Mikul M, Qadir H, Raza S. Persistent Psychosis After a Single Ingestion of «Ecstasy» (MDMA). Prim Care Companion CNS Disord. 2011;13:6.

  • 3.

    Robison AJ, Nestler EJ. Transcriptional and epigenetic mechanisms of addiction. Nat. Rev. Neurosci. 2011;12:623-37. [PubMed ID: 21989194].

  • 4.

    Wagner MT, Mithoefer MC, Mithoefer AT, MacAulay RK, Jerome L, Yazar-Klosinski B, Doblin R. Therapeutic effect of increased openness: Investigating mechanism of action in MDMA-assisted psychotherapy. J. Psychopharmacol. 2017;31:967-74. [PubMed ID: 28635375].

  • 5.

    Nadesan K, Kumari C, Afiq M. Dancing to death: A case of heat stroke. J. Forensic. Leg. Med. 2017;50:1-5. [PubMed ID: 28651196].

  • 6.

    Pence HH, Pence S, Kurtul N, Bakan E, Kok AN, Kocoglu H. The alterations in adenosine nucleotides and lactic acid levels in striated muscles following death with cervical dislocation or electric shock. Soud. Lek. 2003;48:8-11. [PubMed ID: 12673946].

  • 7.

    Verrico CD, Miller GM, Madras BK. MDMA (Ecstasy) and human dopamine, norepinephrine, and serotonin transporters: implications for MDMA-induced neurotoxicity and treatment. Psychopharmacology. 2007:489-503.

  • 8.

    Koob GF, Weiss F. Neuropharmacology of cocaine and ethanol dependence. Recent Dev. Alcohol. 1992;10:201-33. [PubMed ID: 1350359].

  • 9.

    Wise RA, Newton P, Leeb K, Burnette B, Pocock D, Justice JB Jr. Fluctuations in nucleus accumbens dopamine concentration during intravenous cocaine self-administration in rats. 1995. Psychopharmacology (Berl:10-20.

  • 10.

    Morales M, Margolis EB. Ventral tegmental area: cellular heterogeneity, connectivity and behaviour. Nat. Rev. Neurosci. 2017;18:73-85. [PubMed ID: 28053327].

  • 11.

    Ikemoto S. Brain reward circuitry beyond the mesolimbic dopamine system: a neurobiological theory. Neurosci. Biobehav. Rev. 2010;35:129-50. [PubMed ID: 20149820].

  • 12.

    Wise RA. Dopamine and reward: the anhedonia hypothesis 30 years on. Neurotox. Res. 2008;14:169-83. [PubMed ID: 19073424].

  • 13.

    Kelley AE, Berridge KC. The neuroscience of natural rewards: relevance to addictive drugs. J. Neurosci. 2002;22:3306-11. [PubMed ID: 11978804].

  • 14.

    Wise RA, Morales M. A ventral tegmental CRF-glutamate-dopamine interaction in addiction. Brain Res. 2010;1314:38-43. [PubMed ID: 19800323].

  • 15.

    Nestler EJ. Is there a common molecular pathway for addiction? Nat. Neurosci. 2005;8:1445-9.

  • 16.

    EMCDD. European Monitoring Centre for Drugs and Drug Addiction European, Drug Report (2017): Trends and Developments, Publications Office of the European Union, Luxembourg. 2017.

  • 17.

    Thomas MJ, Kalivas PW, Shaham Y. Neuroplasticity in the mesolimbic dopamine system and cocaine addiction. Br. J. Pharmacol. 2008;154:327-42. [PubMed ID: 18345022].

  • 18.

    Pashaei E, Guzel E, Ozgurses ME, Demirel G, Aydin N, Ozen M. A Meta-Analysis: Identification of Common Mir-145 Target Genes that have Similar Behavior in Different GEO Datasets. PLoS One. 2016;11:e0161491. [PubMed ID: 27655328].

  • 19.

    Guzel E, Karatas OF, Semercioz A, Ekici S, Aykan S, Yentur S, Creighton CJ, Ittmann M, Ozen M. Identification of microRNAs differentially expressed in prostatic secretions of patients with prostate cancer. Int. J. Cancer. 2015:136 875-9.

  • 20.

    Bastle RM, Oliver RJ, Gardiner AS, Pentkowski NS, Bolognani F, Allan AM, Chaudhury T, St Peter M, Galles N, Smith C, Neisewander JL, Perrone-Bizzozero NI. In silico identification and in-vivo validation of miR-495 as a novel regulator of motivation for cocaine that targets multiple addiction-related networks in the nucleus accumbens. Mol. Psychiatry. 2018;23:434. [PubMed ID: 28044061].

  • 21.

    Yan B, Hu Z, Yao W, Le Q, Xu B, Liu X, Ma L. MiR-218 targets MeCP2 and inhibits heroin seeking behavior. Sci. Rep. 2017;7:40413. [PubMed ID: 28074855].

  • 22.

    Chandrasekar V, Dreyer JL. Regulation of MiR-124, Let-7d, and MiR-181a in the accumbens affects the expression, extinction, and reinstatement of cocaine-induced conditioned place preference. Neuropsychopharmacology. 2011;36:1149-64. [PubMed ID: 21307844].

  • 23.

    Gomez AM, Altomare D, Sun WL, Midde NM, Ji H, Shtutman M, Turner JR, Creek KE, Zhu J. Prefrontal microRNA-221 Mediates Environmental Enrichment-Induced Increase of Locomotor Sensitivity to Nicotine. Int. J. Neuropsychopharmacol. 2016;19:pyv090.

  • 24.

    Zheng H, Zeng Y, Zhang X, Chu J, Loh HH, Law PY. mu-Opioid receptor agonists differentially regulate the expression of miR-190 and NeuroD. Mol. Pharmacol. 2010;77:102-9. [PubMed ID: 19854889].

  • 25.

    Garcia-Repetto R, Moreno E, Soriano T, Jurado C, Gimenez MP, Menendez M. Tissue concentrations of MDMA and its metabolite MDA in three fatal cases of overdose. Forensic. Sci. Int. 2003;135:110-4. [PubMed ID: 12927411].

  • 26.

    Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25:402-8. [PubMed ID: 11846609].

  • 27.

    Agarwal V, Bell GW, Nam JW, Bartel DP. Predicting effective microRNA target sites in mammalian mRNAs. Elife. 2015;4. doi:10.

  • 28.

    Wang X. miRDB: a microRNA target prediction and functional annotation database with a wiki interface. RNA. 2008;14:1012-7. [PubMed ID: 18426918].

  • 29.

    Dweep H, Gretz N. miRWalk20: a comprehensive atlas of microRNA-target interactions. Nat. Methods. 2015;12:697. [PubMed ID: 26226356].

  • 30.

    Lopez JP, Lim R, Cruceanu C, Crapper L, Fasano C, Labonte B, Maussion G, Yang JP, Yerko V, Vigneault E, Mestikawy S, Mechawar N, Pavlidis P, Turecki G. miR-1202 is a primate-specific and brain-enriched microRNA involved in major depression and antidepressant treatment. Nat Med. 2014;20:764-8. [PubMed ID: 24908571].

  • 31.

    Engel JA, Jerlhag E. Role of appetite-regulating peptides in the pathophysiology of addiction: implications for pharmacotherapy. CNS. Drugs. 2014;28:875-86. [PubMed ID: 24958205].

  • 32.

    Dias C, Feng J, Sun H, Shao NY, Mazei-Robison MS, Damez-Werno D, Scobie K, Bagot R, LaBonté B, Ribeiro E, Liu X, Kennedy P, Vialou V, Ferguson D, Peña C, Calipari ES, Koo JW, Mouzon E, Ghose S, Tamminga C, Neve R, Shen L, Nestler EJ. Beta-catenin mediates stress resilience through Dicer1/microRNA regulation. Nature. 2014;516:51-5. [PubMed ID: 25383518].

  • 33.

    Horzmann KA, Freeman JL. Zebrafish Get Connected: Investigating Neurotransmission Targets and Alterations in Chemical Toxicity. Toxics. 2016;4:3-19.

  • 34.

    Panwar B, Omenn GS, Guan Y. miRmine: a database of human miRNA expression profiles. Bioinformatics. 2017;33:1554-60. [PubMed ID: 28108447].

  • 35.

    Gong J, Wu Y, Zhang X, Liao Y, Sibanda VL, Liu W, Guo A. Comprehensive analysis of human small RNA sequencing data provides insights into expression profiles and miRNA editing. RNA Biol. 2014;11:1375-85. [PubMed ID: 25692236].

  • 36.

    Xie XQ, Wang L, Liu H, Ouyang Q, Fang C, Su W. Chemogenomics knowledgebased polypharmacology analyses of drug abuse related G-protein coupled receptors and their ligands. Front Pharmacol. 2014;5:3. [PubMed ID: 24567719].

  • 37.

    Rudnick G, Kramer R, Blakely RD, Murphy DL, Verrey F. The SLC6 transporters: perspectives on structure, functions, regulation, and models for transporter dysfunction. Pflugers. Arch. 2014;466:25-42. [PubMed ID: 24337881].

  • 38.

    Engmann O, Labonte B, Mitchell A, Bashtrykov P, Calipari ES, Rosenbluh C, Loh YE, Walker DM, Burek D, Hamilton PJ, Issler O, Neve RL, Turecki G, Hurd Y, Chess A, Shen L, Mansuy I, Jeltsch, Akbarian S, Nestler EJ. Cocaine-Induced Chromatin Modifications Associate With Increased Expression and Three-Dimensional Looping of Auts2. Biol. Psychiatry (2017.

  • 39.

    Ron D, Janak PH. GDNF and addiction. Rev. Neurosci. 2005;16:277-85. [PubMed ID: 16519005].

  • 40.

    Tapocik JD, Luu TV, Mayo CL, Wang BD, Doyle E, Lee AD, Lee NH, Elmer GI. Neuroplasticity, axonal guidance and micro-RNA genes are associated with morphine self-administration behavior. Addict. Biol. 2013;18:480-95. [PubMed ID: 22804800].

  • 41.

    Hollander JA, Im HI, Amelio AL, Kocerha J, Bali P, Lu Q, Willoughby D, Wahlestedt C, Conkright MD, Kenny PJ. Striatal microRNA controls cocaine intake through CREB signalling. Nature. 2010;466:197-202. [PubMed ID: 20613834].

  • 42.

    Zhu L, Li J, Dong N, Guan F, Liu Y, Ma D, Goh EL, Chen T1. mRNA changes in nucleus accumbens related to methamphetamine addiction in mice. Sci. Rep. 2016;6:36993. [PubMed ID: 27869204].

  • 43.

    Bosch PJ, Benton MC, Macartney-Coxson D, Kivell BM. mRNA and microRNA analysis reveals modulation of biochemical pathways related to addiction in the ventral tegmental area of methamphetamine self-administering rats. BMC. Neurosci. 2015;16:43.

  • 44.

    Lopez JP, Pereira F, Richard-Devantoy S, Berlim M, Chachamovich E, Fiori LM, Niola P, Turecki G1, Jollant F. Co-Variation of Peripheral Levels of miR-1202 and Brain Activity and Connectivity During Antidepressant Treatment. Neuropsychopharmacology. 2017;42:2043-51. [PubMed ID: 28079059].

  • 45.

    Stilling RM, Fischer A. A Drosophila model for the role of epigenetics in brain function and development. Genome Biol. 2011;12:103. [PubMed ID: 21320297].