Abstract
Cocaine dependence is a complex psychiatric disorder involving both genetic and environmental factors. Several neurotransmitter systems mediate cocaine’s effects, dependence and relapse, being the components of the neurotransmitter release machinery good candidates for the disorder. Previously, we identified a risk haplotype for cocaine dependence in the NSF gene, encoding the protein N-Ethylmaleimide-Sensitive Factor essential for synaptic vesicle turnover. Here we examined the possible contribution to cocaine dependence of a large copy number variant (CNV) that encompasses part of the NSF gene. We performed a case-control association study in a discovery sample (359 cases and 356 controls) and identified an association between cocaine dependence and the CNV (P = 0.013), that was confirmed in the replication sample (508 cases and 569 controls, P = 7.1e-03) and in a pooled analysis (P = 1.8e-04), with an over-representation of low number of copies in cases. Subsequently, we studied the functional impact of the CNV on gene expression and found that the levels of two NSF transcripts were significantly increased in peripheral blood mononuclear cells (PBMC) along with the number of copies of the CNV. These results, together with a previous study from our group, support the role of NSF in the susceptibility to cocaine dependence.
Similar content being viewed by others
Introduction
Cocaine dependence is a complex psychiatric disorder that results from the interaction between genetic and environmental factors, with an estimated heritability of 65–79%1,2,3. Cocaine inhibits dopamine, serotonin and norepinephrine reuptake by binding to their transporters (DAT1, SERT and NET, respectively), increasing their levels at the synaptic cleft4,5. Also, other neurotransmitter pathways participate in cocaine dependence, withdrawal or relapse, including the glutamatergic, GABAergic and endocannabinoid systems6,7,8,9,10,11,12. Their involvement in mediating cocaine effects and dependence makes the molecular machinery that controls neurotransmitter release to the synaptic cleft a strong candidate for participating in the development of cocaine dependence. As part of this complex process that involves many different proteins, the N-ethylmaleimide sensitive factor (NSF) plays an essential role in the synaptic vesicle turnover. NSF is responsible for the recycling of Soluble NSF Attachment Protein Receptor (SNARE) complexes, preventing their accumulation and ensuring the maintenance of intracellular membrane trafficking13.
Previous genetic evidence supports the contribution of single-nucleotide polymorphisms (SNPs) in the NSF gene to cocaine addiction, where a two-SNPs risk haplotype was associated with a rapid development of cocaine dependence, shortly after the first cocaine use14. To date, association studies with SNPs have unraveled only a small fraction of the estimated heritability for cocaine dependence. Thus, exploring structural variation may help to capture part of the missing heritability underlying this complex disorder15. Interestingly, a large copy number variant (CNV) of about 734 kb encompasses the first 13 exons of the NSF gene (including two additional genes, LRRC37A2 and ARL17A) and ranges from 2 to 6 copies per individual (1 to 3 haploid copies)16. In the European population, about 72% of individuals have more than two copies. This CNV is located immediately distal from the MAPT locus at 17q21.31, which was previously found to be subjected to positive selection and rapid evolutionary turnover, and associated with neurological disorders17.
Since CNVs are common forms of genetic variation that often impact gene expression and may have a greater effect on transcriptional regulation than that attributed to SNPs18,19,20, we hypothesize that the CNV spanning NSF has an influence on the expression of the gene and may contribute to cocaine dependence. To our knowledge, association between CNVs and drug dependence has been explored mainly for alcohol or opioid dependence21,22,23,24, but no study has investigated cocaine dependence so far. To address this issue, we explored the role of the CNV encompassing NSF in cocaine addiction through a two stage case-control association study with cocaine dependence. Subsequently, we studied the linkage disequilibrium (LD) between the previous associated SNPs and the CNV, as well as the impact of this CNV on NSF expression.
Results
A CNV encompassing part of the NSF gene was previously reported16. According to published data and information in the Ensembl database (www.ensembl.org), it spans the first 13 exons of the gene. We refined the distal end of the CNV in our population by amplifying two NSF regions by quantitative real-time PCR (qRT-PCR), one in exon 13 and another one in intron 13 (Fig. 1a), and comparing the results with an internal CNV control in exon 2. All subjects were confirmed to bear only two copies of the DNA segment in intron 13, and a variable (and identical) number of copies in exons 13 and 2 (Fig. 1b), confirming that the CNV encompasses the first 13 exons of the gene in all 174 individuals assessed, ranging from 2 to 6 copies per diploid genome (Fig. 2).
In order to explore whether the NSF CNV was associated with cocaine dependence, we performed a case-control association study. The comparison of the number of copies of the CNV between 359 cocaine-dependent patients and 356 controls showed significant differences between groups (P = 0.036, Supplementary Table S1). When we subdivided subjects into ‘low number of copies of the CNV’ (2 and 3) and ‘high number of copies’ (4, 5 and 6), we identified a significant overrepresentation of low number of copies in cocaine-dependent patients (72% vs 63%; P = 0.013; OR = 1.51, CI = [1.10–2.07]). This overrepresentation was replicated in an independent sample (69% vs 63%; P = 7.1e-03; OR = 1.42, CI = [1.10–1.83]). The pooled analysis of both the discovery and replication samples also rendered statistically significant results (71% vs 62%; P = 1.8e-04; OR = 1.459, CI = [1.20 –1.78]), which remained significant after adjusting by age (Padj = 3.1e-03), with a higher frequency of subjects carrying 2 or 3 copies of the NSF CNV in the group of cocaine-dependent subjects when compared to controls (Table 1).
Then, we observed that individuals with low number of copies showed an earlier dependence onset than individuals with high number of copies (P = 0.042). The main differences were observed in the group of patients that develop cocaine dependence within a year after the first drug use (Fig. 3). When we subdivided subjects in two subgroups, early (<1 year) and late (≥1 year) dependence onset, we identified an overrepresentation of individuals with low number of copies in the group of patients showing early regular cocaine consumption (P = 5e-03; OR = 1.62, CI = [1.13–2.32]).
Previously, a SNP located in intron 13 of the NSF gene, rs183211, was associated with cocaine dependence either as a single marker or as part of a risk haplotype with rs17698176 (intron 17; Fig. 2a)14, also associated with an early onset of cocaine dependence. These two SNPs are located distal from the CNV according to our above experiments. Thus, we studied whether the associations between cocaine dependence and either the CNV or the two SNPs in the risk haplotype at NSF were related. To do so, we estimated SNP-CNV haplotype frequencies using an expectation-maximization (EM) algorithm, and subsequently assessed LD between the CNV and rs183211 (Supplementary Table S2) and the CNV and the rs183211-rs17698176 risk haplotype (Table 2). We found that the rs183211G risk allele was associated with low number of copies at the CNV (P = 6e-12) as it was in moderate LD with one copy (r2 = 0.30). Furthermore, similar results were obtained for the risk haplotype: the rs183211G-rs17698176T haplotype allele was associated with low copy number of the CNV (P < 2.2e-16) as this previously reported risk haplotype was in LD with one copy (r2 = 0.61).
The two intronic SNPs previously reported to be associated with cocaine dependence, rs183211 and rs17698176, had no predicted functional effects according to the software SNP Function Prediction (snpinfo.niehs.nih.gov/snpinfo/snpfunc.htm). As we found association between the risk haploallele rs183211G-rs17698176T and low number of CNV copies, we subsequently explored the potential impact of the CNV on the expression of the NSF gene. Expression levels for NSF_001 (full length transcript), NSF_002 (including the first four exons) and NSF_003 (including the last five exons) could be determined in PBMCs. Moreover, NSFP1, that spans the first 13 exons of the gene and seems to be a pseudogene that results from the CNV dynamics, encodes another transcript named NSFP1_001 (Fig. 2b). We analyzed the expression levels of this last transcript with an assay that also detects one of the gene products, NSF_001 (Fig. 4a, primer pair ‘a’). As shown in Fig. 4b, the expression levels of transcripts NSF_001 + NSFP1_001 and NSF_002 were significantly increased along with the number of copies of the CNV (P = 4.3e-06 and P = 8.9e-05, respectively). Similar results were observed when the individuals were grouped into low and high number of copies, with transcript levels increased by 2.1-fold and 1.5-fold for NSF_001 + NSFP1_001 and NSF_002, respectively, in the group of high number of copies as compared to the group of low number of copies (P = 8.5e-09 and P = 3.7e-04, respectively). Finally, we assessed the linearity of this correlation and found a significant association between the number of copies of the CNV and the level of transcripts NSF_001 + NSFP1_001 (ρ = 0.87, P = 9.4e-14) and NSF_002 (ρ = 0.68, P = 1.18e-06). As expected, no significant changes were observed for transcripts NSF_001 and NSF_003, as they are encoded by genomic regions outside the CNV (Fig. 4b).
Discussion
Many studies have been performed to identify SNP variants associated with drug dependence and, even though some of them have reported significant findings, only a small fraction of the estimated heritability of this complex phenotype can be explained by them. Methodological limitations may account in part for the missing heritability, but other sources of genetic variation, like copy number variants, could also contribute to the genetic landscape of the disorder15. In this study, we identified a significant association between cocaine dependence and a CNV spanning the first 13 exons the NSF gene. To our knowledge, this is the first study that associates a CNV with cocaine dependence.
In a previous study, our group identified an association between a quick transition from cocaine use to dependence and a risk haplotype (rs183211G-rs17698176T) in the NSF gene. Furthermore, a large polymorphic CNV partially overlapping the NSF gene has been described16. We thus hypothesized that the NSF CNV could contribute to the susceptibility to cocaine dependence by altering gene expression. First, we performed a case-control association study and found a significant association between this CNV and cocaine dependence, with an overrepresentation of individuals with low number of copies in cocaine patients. Then, we observed that a low number of CNV copies correlated with a rapid development of dependence, and that the CNV and the previously described risk haplotype were in LD. Furthermore, we observed an effect of the CNV on NSF expression, identifying a positive correlation between the number of copies of the CNV and the levels of two transcripts, NSF_002 and NSFP1_001, both encoded by the segment of the gene that is included in the CNV. Since we had no access to brain samples, gene expression was assessed in samples obtained from PBMC, which may differ from brain NSF expression. Some studies, however, reveal a correlation in gene expression levels between peripheral blood and brain25.
In addition, we have characterized the 3′ end of the CNV, located within the NSF gene, and confirmed that the two disease-associated SNPs are located outside the CNV but in moderate LD with it. In consequence, the association identified previously, which was detected under the assumption of diploidy for the SNP variants, was calculated correctly and suggests that the two intronic SNPs may tag the NSF CNV.
CNVs play an important role in the generation of genetic diversity and so in human evolution. They have also been shown to have phenotypic effects, including susceptibility to disease, depending on whether or not dosage-sensitive genes or regulatory sequences are affected by the rearrangement26,27. Many genes involved in brain development are enriched in CNVs28, the biological importance as well as the evolutionary conservation of these genes suggest that they could have phenotype implications and might contribute to the susceptibility to complex neuropsychiatric disorders29,30. Genome-wide studies performed in drug abuse have unraveled a number of CNVs possibly contributing to drug dependence. For example, a CNV at 11q14.2 was found associated with variations in brain volume and drinking behaviour21, and losses at 22q13.1 seem to increase the response to alcohol in the precuneus, a region previously associated with severity of alcohol dependence (AD)22. Furthermore, a GWAS assessing common CNVs in AD in a Korean cohort identified losses at 20q13.33 as protective factors for AD23. Recently, a GWAS performed in opioid dependence (OD) identified many CNVs which could contribute to resistance or susceptibility to OD, three of them genome-wide significant24. These findings suggest that CNVs may have an important role in drug dependence. However, to our knowledge, the present study is the first one performed in cocaine dependence considering CNVs, and obtaining significant results.
The NSF gene encodes the hexameric ATPase N-ethylmaleimide-sensitive factor, which plays an essential role in the synaptic vesicle turnover, modulating the kinetics of neurotransmitters and the integrative properties of synapses31,32,33. NSF operates in the recycling of the SNARE complexes by preventing the accumulation of cis-complexes and assuring a sufficient amount of free SNARE for the maintenance of the intracellular membrane trafficking13. Furthermore, several studies suggest that NSF might modulate the function of several neurotransmitters in the central nervous system by interacting with their receptors, such as the glutamatergic alpha-amino-3-hydroxy-5-methylisoxazole-4-propionic acid receptor (AMPA), the beta-2 adrenergic receptor (β2-AR), the dopamine receptors (more strongly D1 and D5), the adrenomedulin (AM) receptor, the γ-amino-butyric acid receptor (GABAA) and the serotonin transporter (SERT)34,35,36,37,38,39,40,41,42,43,44,45,46. Additionally, a recent study identified interactions between NSF and BLOC-1 (biogenesis of lysosome-related organelles complex 1), which participate in synaptic vesicle assembly, neurotransmission and plasticity, and this interaction modulates synaptic plasticity47.
Interestingly, a recent study showed that withdrawal from repeated cocaine exposure increased interactions between NSF and glutamate receptor 2 (GluR2), and this complex acts as a negative regulator of the expression of behavioral sensitization in rats48. Moreover, the NSF gene was found significantly down-regulated in postmortem hippocampus of cocaine addicts samples, which is consistent with our results showing an overrepresentation of low number of copies (which correlates with lower expression) in our sample of cocaine dependent patients49. Furthermore, NSF decrease correlates with a robust down-regulation of GABBR1 (Gamma-aminobutyric acid type B receptor subunit 1), a subunit of GABAB receptor, which could be used as a predictor for increased risk for addiction49. And in autistic patients, NSF expression was also found significantly reduced in lymphocytes and tended to be reduced in post-mortem brains, correlating with the severity of clinical traits46.
The CNV spanning the NSF gene was previously associated with elevated peak antibody response to the AVA-Biothrax vaccine, suggesting that this CNV has an effect on the expression of the genes in this region (NSF, ARL17 and LRRC37A)50. We have shown that a CNV spanning the NSF gene produces changes in the expression levels of some NSF transcripts, which could have an effect on the availability of neurotransmitter vesicles (loaded with e.g. dopamine, essential in the reward system) and its turnover, producing changes in the neurotransmitter release to the synaptic cleft. This would be in line with the “reward deficiency syndrome”, which hypothesizes that hypodopaminergic activity predisposes to cocaine dependence51. Although DA neurotransmission has been established as a central mediator in cocaine effects52,53,54,55, other neurotransmitters systems, such as the glutamatergic, GABAergic and endocannabinoid ones, also dependent on NSF activity, are very important and both association studies and animal models support their role in cocaine dependence56,57,58,59,60,61.
The present case-control association study raises several methodological issues: (i) sample heterogeneity might alter the results of association studies. To minimize the risks we considered sex-matched Caucasian patients and controls from Spain, all of them recruited in a small geographical area around Barcelona and evaluated following a common clinical assessment; (ii) cocaine dependence was not discarded in the control sample. However, this problem, which tends to balance allele frequencies in cases and controls, would increase the probability of false negative results, which is not the case in the present study; (iii) qRT-PCR assays are widely used to validate CNVs because of its specificity and reproducibility, even though in some cases the results may fall between integers (e.g. copy number measuring 1.4), making interpretation difficult. Ambiguous results were discarded from our study in order to ensure a correct assignation of the number of copies in all cases; (iv) Hardy-Weinberg equilibrium (HWE) could not be tested in our samples since our genotyping strategy did not allow establishing phases. However, the effects of CNVs are usually related to genetic dose, therefore what matters is the total number of copies and not the phased genotype; (v) considering the number of CNV copies as low (2 and 3) or high (4, 5 and 6) is somehow arbitrary (Supplementary Table S1). In any case, if we consider individuals with 4 number of copies into the ‘low’ group, the association remains statistically significant (P = 8.3e-04); (vi) we assessed the effect of the CNV on gene expression in RNA from PBMC, which may differ from expression in the brain; (vi) it is important to note that other genes in addition to NSF lie in the CNV and could also be contributing to cocaine dependence predisposition.
To conclude, we found a significant association between cocaine dependence and a CNV spanning the NSF gene. Functional studies revealed that the levels of two transcripts in blood, NSFP_001 and NSF_002, correlate with the number of copies of the CNV, suggesting that this CNV might have a functional impact on gene expression. These results, together with previous studies by our group14, suggest that NSF contributes to the genetic predisposition to cocaine dependence and, especially, to an early onset of dependence.
Materials and Methods
Delimitation of the copy number variation
The extension of the CNV spanning the NSF gene was ascertained in 174 individuals (89 controls and 85 cases). Detailed information about the subjects is provided in the ‘Association study’ section.
The number of copies at the CNV region and outside was assessed by qRT-PCR considering three regions of the NSF gene: exon 2, exon 13 and intron 13 (the latter located 102 bp 5′ from SNP rs183211, previously associated with cocaine dependence) to narrow down the CNV extension. The albumin gene (ALB) was used as a diploid control. Primer sequences are available in Supplementary Table S3. QRT-PCR assays were performed as follows: 7.5 μl of LightCycler® 480 SYBR Green I Master (Roche Life Sciences, Branford, CT, USA), 400 nM of each primer, 10 ng of DNA and PCR-grade water up to 15 μl. The qRT-PCR conditions were: 10 min at 95 °C, 40 cycles of 95 °C for 15 s and 60 °C for 1 min. Specificity was ascertained after completion of the amplification in the presence of SYBR Green by the melting procedure (ramping from +65 °C to +97 °C, rising the temperature by 2.2 °C at every step with 2-seconds intervals). QRT-PCR experiments were performed using the LightCycler® 480 System and results were analyzed with the LightCycler® 480 Software, Version 1.5 (Roche Life Science, Branford, CT, USA). Standard curves for all amplicons were constructed through triplicate qRT-PCR amplifications with serial dilutions (factor of 2.5) of the DNA template and detection with SYBR Green. The final amount of template in the curve ranged from 300 ng to 3.14e-02 ng for DNA. Standard curves were used to determine the linearity (R2) and the efficiency (E)62 of the qRT-PCR amplifications and all transcripts were amplified with an efficiency around 94% (E = 1.94) and R2 = 0.9964, without significant differences among them. Five HapMap subjects with known number of copies at the NSF CNV (NA18506, NA19239, NA19240, NA12156, NA18942) were used as a reference for 2, 3, 4, 5 and 6 copies, respectively16. Genomic samples of these individuals were obtained from the Coriell Institute for Medical Research (New Jersey, USA) and were included in all assays as an internal control. NSF gene copy number was calculated by the comparative Ct (ΔΔCt) method and all experiments displaying ambiguous results were discarded from the study.
Association study
Subjects
The discovery sample consisted of 359 cocaine-dependent patients (mean age 35.5 ± 8.1 years, 81.3% males (n = 292)) and 356 controls (52.7 ± 16.8 years, 81.5% males (n = 290)). The replication sample consisted of 508 cocaine-dependent patients (mean age 35.7 ± 7.8 years, 79.9% males (n = 406)) and 569 controls (52.1 ± 16.9 years, 79.8% males (n = 454)). All patients were recruited and evaluated at the Psychiatry Department of the Hospital Universitari Vall d’Hebron (Barcelona, Spain) according to DSM-IV TR criteria (Diagnostic and Statistical Manual of Mental Disorders, 4th ed., text revision) and the SCID (Structured Clinical Interview)63. A total of 75% of patients (n = 652) reported age at the initial consumption as well as age at dependence onset. All controls were recruited at the Blood and Tissues Bank of Barcelona. None of them had injected drugs intravenously. Patients and controls were Spanish, Caucasian and sex-matched. The study was approved by the ethics committee of our institution, the Institutional Review Board of the University of Barcelona (IRB00003099), and informed consent was obtained from all participants, in accordance with the Helsinki Declaration. Population stratification was previously discarded in our sample64.
DNA isolation and quantification
Genomic DNA was isolated from peripheral blood mononuclear cells (PBMC) using the salting-out method65. DNA concentrations were determined using NanoDrop ND-1000 Spectrophotometer (NanoDrop Technologies, Termo Fisher Scientific Inc., Wilmington, DE, USA).
Copy number assay
The number of copies of the CNV of each individual was assessed by qRT-PCR considering a region of exon 2 in the NSF gene. The albumin gene (ALB) was used as a diploid control. QRT-PCR conditions are described above. In all assays, samples with known number of copies at the NSF CNV were included, showing a concordance rate of 100%. All experiments displaying ambiguous results were not considered in the study.
Statistical analysis
The comparison of the number of copies of the CNV between cases and controls was performed using the Fisher’s exact test (two-tailed), considering the five possible copy number assignations per individual (from 2 to 6) and also a dichotomization into low (2 or 3 copies) and high (4, 5 or 6 copies) copy number. Age was considered as a covariate using binary logistic regression test. To analyze the effect of the number of copies of the CNV in the lapse between initial and regular consumption of cocaine, the comparison of medians was performed using the Mann-Whitney U non-parametric test as normality was rejected using Kolmgorov-Smirnov test. In addition, the time between initial consumption and onset of cocaine dependence was dichotomized intro early (<1 years) or late (≥1 years) dependence onset, and the comparison of the number of copies of the NSF CNV between these groups was performed with the Fisher’s exact test (one-tailed). Statistical analysis was performed using SPSS Statistics Version 22.0 software (IBM Corp. Released 2013. Armonk, NY).
Linkage disequilibrium analysis
We estimated the allele and haplotype frequencies of the CNV and two NSF SNPs found previously associated with cocaine dependence (rs183211 and rs17698176)14. Genotypes for both SNPs and CNV were available for 47% of the sample used in the CNV association study (n = 407, 36.07 ± 9.5 years, 85.7% males). Traditional approaches to assess LD between SNPs and CNVs are insufficient for higher order CNVs as the phase of the CNV cannot be directly inferred. To address this issue, we developed an EM algorithm-based method to estimate SNP-CNV haplotype frequencies under an assumption of Hardy-Weinberg equilibrium between haplotypes. The approach is similar to an existing method66. The algorithm begins by seeding initial SNP-CNV haplotype frequencies according to known SNP-CNV genotype counts (e.g., an individual with a homozygote AA SNP genotype and 2 copies of the CNV would likely have two copies of a haplotype which included the A SNP allele and a single copy of the CNV). In the E-step, haplotype counts are calculated based on the estimated haplotype frequencies; followed by the M-step, which recalculates haplotype frequencies based on the estimated haplotype counts. This process is repeated in an iterative fashion until the haplotype frequencies become constant from one iteration to the next. In this study, the frequency estimates were constant to five decimal places by the twentieth iteration. To estimate these haplotype frequencies, genotype data for rs183211 and rs17698176 from the 407 individuals previously described was used. From the resulting haplotype frequencies, we then assessed linkage disequilibrium (r2) between rs183211 and low number of copies of the CNV (individuals with 2 and 3 copies). Furthermore, we likewise assessed LD between the rs183211-rs17698176 haplotype and low number of copies of the CNV.
NSF expression assays
To correlate the number of copies of the CNV with NSF transcription levels, 45 Spanish, Caucasian, healthy volunteers (mean age 32.8 ± 11.8 years, 66.7% females (n = 30)) were recruited at the Genetics Department of the University of Barcelona (Barcelona, Spain). Genomic DNA was isolated as previously described. Total RNA was extracted from PBMC within 2 h after sampling using BD Vacutainer® CPTTM tubes (BD, Franklin Lakes, NJ, USA) and the High Pure RNA Isolation Kit (Roche Life Science, Branford, CT, USA). RNA was retro-transcribed to cDNA using the High Capacity cDNA Reverse Transcription Kit (Applied Biosystems, Foster City, CA, USA). Specific primers were designed to detect each of the eight NSF transcripts according to ENSEMBL, as well as NSFP1, a pseudogene that probably arose as a consequence of the CNV dynamics. Three transcripts could be amplified specifically: NSF_001 (primers set ‘b’ in Fig. 4a), NSF_002 (primers set ‘c’) and NSF_003 (primers set ‘d’) (ENST00000398238, ENST00000486366 and ENST00000465370, respectively, according to Ensembl release 83). For NSFP1_001 (ENST00000570034.1) we could not design specific primers since it overlaps with NSF_001; for this reason we designed a qRT-PCR assay to detect simultaneously both transcripts (primers set ‘a’). Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and hypoxanthine phosphoribosyltransferase (HPRT1) were used to normalize the relative amounts of mRNA. QRT-PCR assays were performed as described above and all transcripts were amplified with an efficiency around 94% (E = 1.94) and R2 = 0.9963, without significant differences among them.
Statistical analysis
Expression levels of each transcript were compared among the five groups generated according to the number of copies of the CNV per individual (from 2 to 6) or between low and high CNV copies using Kruskal-Wallis and Mann-Whitney U non-parametric tests, respectively. In both cases normality was previously rejected using the Kolmogorov-Smirnov test. To estimate the relationship among the number of copies of the CNV and the expression levels of each transcript, a linear regression analysis was performed and its correlation was confirmed using the Spearman’s rank correlation coefficient (ρ), threshold set at P < 0.05 in all tests. Statistical analysis was performed using SPSS Statistics Version 22.0 software (IBM Corp. Released 2013. Armonk, NY).
Ethics statement.
This study was approved by the local Ethics Committee and informed consent was obtained from all adult subjects, children, and their parents according to the Helsinki declaration.
Additional Information
How to cite this article: Cabana-Domínguez, J. et al. A highly polymorphic copy number variant in the NSF gene is associated with cocaine dependence. Sci. Rep. 6, 31033; doi: 10.1038/srep31033 (2016).
References
Kendler, K. S. & Prescott, C. A. Cocaine use, abuse and dependence in a population-based sample of female twins. Br. J. Psychiatry 173, 345–350 (1998).
Kendler, K. S., Karkowski, L. M., Neale, M. C. & Prescott, C. A. Illicit Psychoactive Substance Use, Heavy Use, Abuse, and Dependence in a US Population-Based Sample of Male Twins. Arch. Gen. Psychiatry 57, 261–269 (2000).
Wagner, F. A. & Anthony, J. C. From first drug use to drug dependence: Developmental periods of risk for dependence upon marijuana, cocaine, and alcohol. Neuropsychopharmacology 26, 479–488 (2002).
Kuhar, M. J., Ritz, M. C. & Boja, J. W. The dopamine hypothesis of the reinforcing properties of cocaine. Trends Neurosci. 14, 299–302 (1991).
Kalivas, P. W. Neurobiology of cocaine addiction: implications for new pharmacotherapy. Am. J. Addict. 16, 71–78 (2007).
Cornish, J. L. & Kalivas, P. W. Glutamate transmission in the nucleus accumbens mediates relapse in cocaine addiction. J. Neurosci. 20, RC89 (2000).
Gonzalez, G. et al. Tiagabine increases cocaine-free urines in cocaine-dependent methadone-treated patients: results of a randomized pilot study. Addiction 98, 1625–1632 (2003).
Morgan, A. E. & Dewey, S. L. Effects of pharmacologic increases in brain GABA levels on cocaine-induced changes in extracellular dopamine. Synapse 28, 60–65 (1998).
Dewey, S. L. et al. GABAergic attenuation of cocaine-induced dopamine release and locomotor activity. Synapse 25, 393–398 (1997).
De Vries, T. J. & Schoffelmeer, A. N. M. Cannabinoid CB1 receptors control conditioned drug seeking. Trends Pharmacol. Sci. 26, 420–426 (2005).
De Vries, T. J. et al. A cannabinoid mechanism in relapse to cocaine seeking. Nat. Med. 7, 1151–1154 (2001).
Dackis, C. A. et al. Modafinil and cocaine: a double-blind, placebo-controlled drug interaction study. Drug Alcohol Depend. 70, 29–37 (2003).
Barszczewski, M. et al. A novel site of action for alpha-SNAP in the SNARE conformational cycle controlling membrane fusion. Mol. Biol. Cell 19, 776–784 (2008).
Fernàndez-Castillo, N. et al. Candidate pathway association study in cocaine dependence: the control of neurotransmitter release. World J. Biol. Psychiatry 13, 126–134 (2012).
Manolio, T. A. et al. Finding the missing heritability of complex diseases. Nature 461, 747–753 (2009).
Sudmant, P. H. et al. Diversity of human copy number variation and multicopy genes. Science 330, 641–646 (2010).
Stefansson, H. et al. A common inversion under selection in Europeans. Nat. Genet. 37, 129–137 (2005).
Perry, G. H. et al. Diet and the evolution of human amylase gene copy number variation. Nat. Genet. 39, 1256–1260 (2007).
Henrichsen, C. N., Chaignat, E. & Reymond, A. Copy number variants, diseases and gene expression. Hum. Mol. Genet. 18, 1–8 (2009).
Schlattl, A., Anders, S., Waszak, S. M., Huber, W. & Korbel, J. O. Relating CNVs to transcriptome data at fine resolution: Assessment of the effect of variant size, type, and overlap with functional regions. Genome Res. 21, 2004–2013 (2011).
Boutte, D. et al. Association of genetic copy number variations at 11 q14.2 with brain regional volume differences in an alcohol use disorder population. Alcohol 46, 519–527 (2012).
Liu, J., Calhoun, V. D., Chen, J., Claus, E. D. & Hutchison, K. E. Effect of homozygous deletions at 22q13.1 on alcohol dependence severity and cue-elicited BOLD response in the precuneus. Addict. Biol. 18, 548–558 (2013).
Bae, J. S. et al. The genetic effect of copy number variations on the risk of alcoholism in a korean population. Alcohol. Clin. Exp. Res. 36, 35–42 (2012).
Li, D. et al. Genome-Wide Association Study of Copy Number Variations (CNVs) with Opioid Dependence. Neuropsychopharmacology 40, 1016–1026 (2015).
Sullivan, P. F., Fan, C. & Perou, C. M. Evaluating the comparability of gene expression in blood and brain. Am. J. Med. Genet. - Neuropsychiatr. Genet. 141B, 261–268 (2006).
Lupski, J. R. & Stankiewicz, P. Genomic disorders: Molecular mechanisms for rearrangements and conveyed phenotypes. PLoS Genet. 1, 627–633 (2005).
Stankiewicz, P. & Lupski, J. R. Structural variation in the human genome and its role in disease. Annu. Rev. Med. 61, 437–455 (2010).
de Smith, A. J. et al. Array CGH analysis of copy number variation identifies 1284 new genes variant in healthy white males: Implications for association studies of complex diseases. Hum. Mol. Genet. 16, 2783–2794 (2007).
de Smith, a J., Walters, R. G., Froguel, P. & Blakemore, a I. Human genes involved in copy number variation: mechanisms of origin, functional effects and implications for disease. Cytogenet. Genome Res. 123, 17–26 (2008).
Malhotra, D. & Sebat, J. CNVs: Harbingers of a rare variant revolution in psychiatric genetics. Cell 148, 1223–1241 (2012).
Schweizer, F. E. et al. Regulation of neurotransmitter release kinetics by NSF. Science 279, 1203–1206 (1998).
Malsam, J., Kreye, S. & Söllner, T. H. Membrane traffic in the secretory pathway: Membrane fusion: SNAREs and regulation. Cell. Mol. Life Sci. 65, 2814–2832 (2008).
Littleton, J. T. et al. SNARE-complex disassembly by NSF follows synaptic-vesicle fusion. Proc. Natl. Acad. Sci. USA 98, 12233–12238 (2001).
Bomberger, J. M., Parameswaran, N., Hall, C. S., Aiyar, N. & Spielman, W. S. Novel function for Receptor Activity-modifying Proteins (RAMPs) in post-endocytic receptor trafficking. J. Biol. Chem. 280, 9297–9307 (2005).
Chou, W.-H. et al. GABAA receptor trafficking is regulated by protein kinase C(epsilon) and the N-ethylmaleimide-sensitive factor. J. Neurosci. 30, 13955–13965 (2010).
Cong, M. et al. Binding of the β2 Adrenergic Receptor to N-Ethylmaleimide-sensitive Factor Regulates Receptor Recycling. J. Biol. Chem. 276, 45145–45152 (2001).
Evers, D. M. et al. Plk2 attachment to NSF induces homeostatic removal of GluA2 during chronic overexcitation. Nat. Neurosci. 13, 1199–1207 (2010).
Hanley, J. G., Khatri, L., Hanson, P. I. & Ziff, E. B. NSF ATPase and alpha-/beta-SNAPs disassemble the AMPA receptor-PICK1 complex. Neuron 34, 53–67 (2002).
Heydorn, A. et al. Distinct in vitro interaction pattern of dopamine receptor subtypes with adaptor proteins involved in post-endocytotic receptor targeting. FEBS Lett. 556, 276–280 (2004).
Kittler, J. T. et al. The subcellular distribution of GABARAP and its ability to interact with NSF suggest a role for this protein in the intracellular transport of GABA(A) receptors. Mol. Cell. Neurosci. 18, 13–25 (2001).
Leil, T. A., Chen, Z. W., Chang, C. S. S. & Olsen, R. W. GABAA receptor-associated protein traffics GABAA receptors to the plasma membrane in neurons. J. Neurosci. 24, 11429–11438 (2004).
Nishimune, A. et al. NSF Binding to GluR2 Regulates Synaptic Transmission. Neuron 21, 87–97 (1998).
Osten, P. et al. The AMPA receptor GluR2 C terminus can mediate a reversible, ATP-dependent interaction with NSF and alpha- and beta-SNAPs. Neuron 21, 99–110 (1998).
Pontier, S. M. et al. Coordinated action of NSF and PKC regulates GABAB receptor signaling efficacy. EMBO J. 25, 2698–2709 (2006).
Song, I. et al. Interaction of the N-ethylmaleimide-sensitive factor with AMPA receptors. Neuron 21, 393–400 (1998).
Iwata, K. et al. N-ethylmaleimide-sensitive factor interacts with the serotonin transporter and modulates its trafficking: implications for pathophysiology in autism. Mol. Autism 5, 33 (2014).
Gokhale, A. et al. The N-ethylmaleimide-sensitive factor and dysbindin interact to modulate synaptic plasticity. J. Neurosci. 35, 7643–76453 (2015).
Lu, H. et al. Interactions between N- Ethylmaleimide-Sensitive Factor and GluR2 in the Nucleus Accumbens Contribute to the Expression of Locomotor Sensitization to Cocaine. 34, 3493–3508 (2014).
Enoch, M. A. et al. GABAergic gene expression in postmortem hippocampus from alcoholics and cocaine addicts; corresponding findings in alcohol-naïve P and NP rats. PLoS One 7, 1–11 (2012).
Falola, M. I. et al. Genomic Copy Number Variants: Evidence for Association with Antibody Response to Anthrax Vaccine Adsorbed. PLoS One 8 (2013).
Comings, D. E. & Blum, K. Reward deficiency syndrome: genetic aspects of behavioral disorders. Prog. Brain Res. 126, 325–341 (2000).
Woolverton, W. L. & Johnson, K. M. Neurobiology of cocaine abuse. Trends Pharmacol. Sci. 13, 193–200 (1992).
Volkow, N. D. et al. Relationship between psychostimulant-induced ‘high’ and dopamine transporter occupancy. Proc. Natl. Acad. Sci. USA 93, 10388–10392 (1996).
Volkow, N. D., Fowler, J. S. & Wang, G. J. Imaging studies on the role of dopamine in cocaine reinforcement and addiction in humans. J. Psychopharmacol. 13, 337–345 (1999).
Volkow, N. D., Fowler, J. S., Wang, G. J. & Goldstein, R. Z. Role of dopamine, the frontal cortex and memory circuits in drug addiction: insight from imaging studies. Neurobiol. Learn. Mem. 78, 610–624 (2002).
Gass, J. T. & Olive, M. F. Glutamatergic substrates of drug addiction and alcoholism. Biochem. Pharmacol. 75, 218–265 (2008).
Kalivas, P. W. & O’Brien, C. Drug addiction as a pathology of staged neuroplasticity. Neuropsychopharmacology 33, 166–180 (2008).
Thomas, M. J., Kalivas, P. W. & Shaham, Y. Neuroplasticity in the mesolimbic dopamine system and cocaine addiction. Br. J. Pharmacol. 154, 327–342 (2008).
Uys, J. D. & LaLumiere, R. T. Glutamate: the new frontier in pharmacotherapy for cocaine addiction. CNS Neurol . Disord. Drug Targets 7, 482–491 (2008).
Soria, G. et al. Lack of CB1 cannabinoid receptor impairs cocaine self-administration. Neuropsychopharmacology 30, 1670–1680 (2005).
Xi, Z.-X. & Stein, E. A. Nucleus accumbens dopamine release modulation by mesolimbic GABAA receptors—an in vivo electrochemical study. Brain Res. 798, 156–165 (1998).
Livak, K. J. & Schmittgen, T. D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 25, 402–408 (2001).
First, M. B., Spitzer, S. R., Gibbon, M. & Williams, J. B. M. Structured Clinical Interview for DSM-IV Axis I Disorders (SCID). (American Psychiatric Press, 1997).
Fernàndez-Castillo, N. et al. Association study of 37 genes related to serotonin and dopamine neurotransmission and neurotrophic factors in cocaine dependence. Genes. Brain. Behav. 12, 39–46 (2013).
Miller, S. A., Dykes, D. D. & Polesky, H. F. A simple salting out procedure for extracting DNA from human nucleated cells. Nucleic Acids Res. 16, 1215 (1988).
Wineinger, N. E., Pajewski, N. M. & Tiwar, H. K. A method to assess linkage disequilibrium between CNVs and SNPs inside copy number variable regions. Front. Genet. 2, 1–7 (2011).
Acknowledgements
We are grateful to all the patients, their families and clinical collaborators who contributed with the patient recruitment (Nieves Martinez-Luna, R. Felipe Palma-Álvarez, Constanza Daigre, Nina Vela and Alejandra Herrero). We thank Mónica Cozar for her help with blood extractions and Gemma Marfany for helpful suggestions on gene expression analysis. N.F.-C., J.C.-D. and B.T. were supported by contracts of the ‘Centro de Investigación Biomédica en Red de Enfermedades Raras’ (CIBERER), and J.C.-D. was additionally supported by ‘Generalitat de Catalunya’ (2015 FI_B 00448). B.T. was supported by AGAUR (Generalitat de Catalunya) and by H2020 (grant agreement 667302, CoCA). Major financial support for this research was received by B.C. from the Spanish ‘Ministerio de Economía y Competitividad’ (SAF2015-68341-R) and AGAUR, ‘Generalitat de Catalunya’ (2014SGR932). M.R. is a recipient of a Miguel de Servet contract from the ‘Instituto de Salud Carlos III’, Spain (CP09/00119) and received financial support from ‘Instituto de Salud Carlos III’ (PI11/00571, PI11/01629, PI12/01139, PI14/01700), and by the European Regional Development Fund (ERDF) and AGAUR, ‘Generalitat de Catalunya’ (2014SGR1357). C.R. received financial support by ‘Instituto de Salud Carlos III’ (PI13/1911) and ‘Plan Nacional Sobre Drogas’ (2013/044). N.E.W. is supported by the NIH/NCATS grant 5UL1 TR001114. We also received funding from the European Community’s Seventh Framework Program (under grant agreement number 602805, Aggressotype) and from the European Community’s H2020 Program (under grant agreement number 667302, CoCA).
Author information
Authors and Affiliations
Contributions
J.C.-D. performed the genotyping assay of the CNV in the NSF gene and the association study, the CNV delimitation, the expression study of NSF and the statistical analysis; C.R., L.G.-L., L.R.-C., C.B., C.A. and M.C. participated in the recruitment of patients and clinical assessment and coordinated the clinical research; C.A. participated in the statistical analysis; G.E. and N.E.W. performed the linkage disequilibrium analysis; J.C.-D., B.T. and N.F.-C. isolated genomic DNA from samples; J.C.-D., B.C. and N.F.-C. designed and coordinated the study; J.C.-D. prepared the first draft of the manuscript and all figures and tables; M.R., B.C. and N.F.-C. supervised the manuscript preparation. All authors contributed to and approved the final manuscript.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Cabana-Domínguez, J., Roncero, C., Grau-López, L. et al. A Highly Polymorphic Copy Number Variant in the NSF Gene is Associated with Cocaine Dependence. Sci Rep 6, 31033 (2016). https://doi.org/10.1038/srep31033
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/srep31033
This article is cited by
-
Different MAPT haplotypes influence expression of total MAPT in postmortem brain tissue
Acta Neuropathologica Communications (2023)
-
Genome-wide meta-analysis of copy number variations with alcohol dependence
The Pharmacogenomics Journal (2018)