Abstract
The synaptosomal-associated protein SNAP25 is a key player in synaptic vesicle docking and fusion and has been associated with multiple psychiatric conditions, including schizophrenia, bipolar disorder, and attention-deficit/hyperactivity disorder. We recently identified a promoter variant in SNAP25, rs6039769, that is associated with early-onset bipolar disorder and a higher gene expression level in human prefrontal cortex. In the current study, we showed that this variant was associated both in males and females with schizophrenia in two independent cohorts. We then combined in vitro and in vivo approaches in humans to understand the functional impact of the at-risk allele. Thus, we showed in vitro that the rs6039769 C allele was sufficient to increase the SNAP25 transcription level. In a postmortem expression analysis of 33 individuals affected with schizophrenia and 30 unaffected control subjects, we showed that the SNAP25b/SNAP25a ratio was increased in schizophrenic patients carrying the rs6039769 at-risk allele. Last, using genetics imaging in a cohort of 71 subjects, we showed that male risk carriers had an increased amygdala–ventromedial prefrontal cortex functional connectivity and a larger amygdala than non-risk carriers. The latter association has been replicated in an independent cohort of 121 independent subjects. Altogether, results from these multilevel functional studies are bringing strong evidence for the functional consequences of this allelic variation of SNAP25 on modulating the development and plasticity of the prefrontal–limbic network, which therefore may increase the vulnerability to both early-onset bipolar disorder and schizophrenia.
SIGNIFICANCE STATEMENT Functional characterization of disease-associated variants is a key challenge in understanding neuropsychiatric disorders and will open an avenue in the development of personalized treatments. Recent studies have accumulated evidence that the SNARE complex, and more specifically the SNAP25 protein, may be involved in psychiatric disorders. Here, our multilevel functional studies are bringing strong evidence for the functional consequences of an allelic variation of SNAP25 on modulating the development and plasticity of the prefrontal–limbic network. These results demonstrate a common genetically driven functional alteration of a synaptic mechanism both in schizophrenia and early-onset bipolar disorder and confirm the shared genetic vulnerability between these two disorders.
Introduction
Bipolar disorder (BD) and schizophrenia (SZ) share many risk factors, cognitive features, neuroanatomical characteristics, clinical symptomatology, and medication. Shared risk factors include genetic risk factors with shared heritability and some common single nucleotide polymorphisms (SNPs), increasing the risk for both conditions (Lichtenstein et al., 2009; Purcell et al., 2009). Recent genome-wide association studies (GWASs) and whole exome sequence analyses conducted in SZ and BD revealed enrichment in genes encoding molecules involved in synaptic neurotransmission (Psychiatric GWAS Consortium Bipolar Disorder Working Group, 2011; Purcell et al., 2014; Schizophrenia Working Group of the Psychiatric Genomics, 2014). One of the key steps of neurotransmission is the Ca2+-dependent assembly of the synaptosomal-associated protein SNAP25 with synaptobrevin and syntaxin to form the soluble N-ethylmaleimide-sensitive factor attachment protein receptor (SNARE) complex. Multiple SNPs located in SNAP25 have already been associated with SZ (Carroll et al., 2009; Wang et al., 2015) and other psychiatric disorders (Barr et al., 2000; Terracciano et al., 2010; Braida et al., 2015). In addition, several studies have reported alteration of the SNAP25 protein level in the brains of patients with either SZ or BD (Thompson et al., 1998, 2003a; Fatemi et al., 2001; Honer et al., 2002; Scarr et al., 2006; Gray et al., 2010). In particular, alterations at the hippocampal level have been reported in both conditions (Young et al., 1998; Fatemi et al., 2001; Thompson et al., 2003a; Scherk et al., 2008).
At a neural level, common alterations in SZ and BD are mainly frontolimbic (d'Albis and Houenou, 2015). Hippocampal volumes are generally described as decreased in both conditions (Hibar et al., 2016; van Erp et al., 2016), while findings are more heterogeneous for amygdala, with an age-dependent directionality (Hajek et al., 2009; Welch et al., 2010; Ganzola et al., 2014). Prefrontal–limbic connectivity has also repeatedly been reported to be altered in both conditions, probably leading to emotional dysregulation in both conditions (Hart et al., 2013; Cannon, 2015; d'Albis and Houenou, 2015).
We recently identified one polymorphism located in the promoter region of SNAP25 that was associated with early-onset BD (EOBD) as well as with a higher expression level in the prefrontal cortex of subjects carrying the at-risk allele (Etain et al., 2010). Among the various forms of BD, EOBD is probably a more homogeneous subtype, more severe and more frequently associated with psychotic features, genetic heritability, and neurodevelopmental abnormalities than forms with a later age of onset. Thus, EOBD may be a neurodevelopmental subtype of BD, closer to SZ than BD. Given this, we hypothesized that the rs6039769 SNAP25 at-risk variant for EOBD may be associated with an increased risk of SZ. However, one usual limitation of genetic association studies is the absence of information on the plausibility and nature of the relation between a genotype and a phenotype. Functional studies are thus necessary to bring support and shed light on such links. We thus aimed at understanding the functional impact of this DNA at-risk variant combining multiple approaches (i.e., gene expression analyses both in vitro and in brain postmortem tissues as well as in vivo using a clinical human genetics imaging approach), demonstrating that this at-risk variant is associated with common neuroanatomical features of EOBD and SZ, namely alterations in the prefrontal limbic networks.
Materials and Methods
Subjects
As false-positive results are a recurrent concern in both the genetics and neuroimaging fields (Meyer-Lindenberg et al., 2008), we used the following two independent samples: a “discovery” sample and a “replication” sample for both genetic analyses and neuroimaging genetic analyses for replication.
Genetic analyses
Discovery sample.
Four hundred sixty-one patients (67.9% males) meeting the Diagnostic and Statistical Manual of Mental Disorders, Fourth Edition, Text Revision (DSM-IV-TR) criteria (American Psychiatric Association, 2000) for SZ and 315 unaffected control subjects (54.6% males) have been recruited as described previously (Dumaine et al., 2011; Schürhoff et al., 2015). All subjects were interviewed by trained investigators (psychiatrists or psychologists) in 10 Expert Centres for SZ of the Fondation FondaMental (https://www.fondation-fondamental.org), using the diagnostic interview for genetic studies (DIGS version 3), yielding DSM-IV diagnoses (Nurnberger et al., 1994).
Replication sample.
The latest Psychiatric Genomics Consortium GWAS summary statistics for rs6039769 have been downloaded from the Broad Institute web site (https://data.broadinstitute.org/mpg/ricopili/). To match with the geographical origin of our discovery sample, only data from the 49 Caucasian populations, including 34,241 patients with SZ and 45,604 control subjects, have been considered for the genetic analysis (Schizophrenia Working Group of the Psychiatric Genomics, 2014).
Neuroimaging genetics analyses
Discovery sample.
Twenty-five adult inpatients and outpatients with BD type I (by DSM-IV-R criteria) were recruited from two university-affiliated participating centers (AP-HP Henri Mondor and Lariboisière Hospitals, Paris, France; mean age, 38.9 years; SD, 12.1; 10 females); 46 control subjects were recruited from media announcements and registry offices, and had no personal or family history of Axis I mood disorder, SZ, or schizoaffective disorder (mean age, 38.7 years; SD, 11.4; age comparison with patients: t(69) = 0.07, p = 0.93; 27 females χ2 = 2.27, p = 0.13). All participants underwent clinical assessment by trained raters using the DIGS. Exclusion criteria for all participants consisted of a history of neurological disease or head trauma with loss of consciousness and contraindications for magnetic resonance imaging (MRI).
Replication sample.
An independent sample of 121 healthy control subjects (mean age, 29.7 years; SD. 10.8; 65 women) has been recruited at the Central Institute of Mental Health in Mannheim, Germany. A comparison with the discovery sample yielded a significant difference in mean age (the replication sample comprising younger subjects: t(190) = 5.46, p = 1.5 × 10−7).
Local ethics committees (CPP Ile de France IX and the ethics committee at the University of Heidelberg, Heidelberg, Germany) approved the protocols and procedures, and written informed consent was obtained from all subjects before study participation.
Genetic association
Genomic DNA has been extracted from blood lymphocyte, B-lymphoblastoid cell line, or from saliva for each subject. We performed the rs6039769 genotyping using TaqMan SNP genotyping assays (C_29497348_10) on a 7900 Real-Time PCR System (Applied Biosystems).
Functional in vitro study
The Neuro-2a (N2a) mouse neuroblastoma cell line (catalog #400394/p451_Neuro-2A, CLS; RRID:CVCL_0470) was grown at 37°C, with 5% CO2 in DMEM containing GlutaMAX (Thermo Fisher Scientific) and was supplemented with 10% fetal bovine serum, 100 U/ml penicillin, and 100 μg/ml streptomycin. The plasmid containing the 784 bp upstream of the SNAP25 transcription start site has been cloned in the pGL3-Basic vector (Promega) upstream of the luciferase reporter gene as described previously (Cai et al., 2008) and was provided by Professor Weihong Song (The University of British Columbia, Vancouver, Canada). The rs6039769 polymorphism was generated by in vitro mutagenesis using the QuikChange II XL Site-Directed Mutagenesis Kit (Agilent Technologies). N2a cells were seeded onto 24-well plates 1 d before transfection and were transfected with 50 ng of plasmid DNA per well using 1 μl of Lipofectamine 2000 (Thermo Fisher Scientific) and 1 ng of Renilla luciferase vector to normalize the transfection efficiency. Cells were harvested at 24 h after transfection and lysed with 100 μl of passive lysis buffer (Promega) per well. Luciferase activities were measured sequentially using the Dual-Luciferase Reporter Assay System (Promega).
SNAP25 postmortem expression study
RNA, cDNA, and DNA from 33 individuals affected with SZ (25 males and 8 females) and 30 unaffected control subjects (23 males and 7 females) were donated by the Stanley Medical Research Institute (Chevy Chase, MD), as part of the Array Collection that consisted of samples from the dorsolateral prefrontal cortex (Brodmann's area 46; Torrey et al., 2000). Diagnoses were made according to the DSM-IV. The samples were coded, and genotypes and disease status were known only after expression analyses. Expression levels of the two isoforms of SNAP25 were determined using TaqMan gene expression assays as described previously (Etain et al., 2010). Normalization was performed using an endogenous housekeeping gene encoding the human β-actin (ACTB) with limited primers (Applied Biosystems). As the ratio between the two isoforms of SNAP25 has been shown to play a critical role in the synaptic maturation at puberty (Bark et al., 2004), a key step for SZ onset, we compared the ratio between SNAP25b and SNAP25a depending on the genotype.
Experimental design and statistical analysis
Statistical testing for allelic and genotypic associations was performed using the PLINK version 1.07 software (RRID:SCR_001757; http://zzz.bwh.harvard.edu/plink/). We performed statistical analyses of expression using the statistical software package R version 2.9.2 (RRID:SCR_003005; http://cran.r-project.org). Groupwise comparisons included parametric tests (Student's t test/ANOVA) and nonparametric tests (Mann–Whitney U test/Kruskal–Wallis test), according to the normality of distribution, and were tested using the Shapiro–Wilk method. A potential correlation between SNAP25a and SNAP25b expression levels and age, postmortem interval, refrigerator interval, and brain pH has been tested using Spearman's rank correlation test in case patients and control subjects.
Neuroimaging genetics study
Structural and resting-state fMRI (rs-fMRI) scans were collected at Neurospin (Commissariat à l'Énergie Atomique, Saclay, France) for the discovery sample, using a 3 T MRI (12-channel head coil; Tim Trio, Siemens) with the following acquisition parameters: 3DT1 sequence: TE = 2.98 ms; TR = 2.3 s; 160 sections; voxel size = 1 × 1 × 1.1 mm; rs-fMRI: echoplanar imaging; voxel size = 3 × 3 × 3 mm; 35 slices; TR = 2000 ms; 360 time points (12 min). For rs-fMRI, subjects were instructed to close their eyes and try not to think about anything specific without sleeping.
Structural scans were acquired at the Mannheim Central Institute for Mental Health (Mannheim, Germany) for the replication sample, using the same 3 T MRI as in Paris (12-channel head coil; Tim Trio, Siemens) with the same acquisition parameters.
Automated volumetry, as implemented in FIRST 1.2, part of FSL 4.1 software (Patenaude et al., 2011), was applied on T1 images to segment the hippocampi and amygdala. Registrations and segmentations were visually inspected as recommended by the authors of the software. The volumes of hippocampi and amygdalae were measured on segmented images. As our results with the discovery sample concerned only the amygdala, we took advantage of the larger replication sample to test for lateralized associations of the amygdala volumes with the genotype.
The rs-fMRI images were preprocessed using Data Processing Assistant for Resting-State fMRI (DPARSF) toolbox (RRID:SCR_002372; http://www.restfmri.net) with Statistical Parametric Mapping (SPM8; RRID:SCR_007037; http://www.fil.ion.ucl.ac.uk/spm/). Briefly, the 10 first volumes were discarded for the stabilization of the signal. Functional volumes were corrected for slice timing acquisition and head motion (>2 mm), then were realigned on T1 images and segmented into six tissue types (white matter, gray matter, CSF, bone, and air). We then spatially normalized these images to the Montreal Neurological Institute space using Dartel and smoothed them by a 4 mm full-width at half-maximum. We performed a linear trend using a temporal bandpass filter (0.001–0.08 Hz). Finally, we used signal in gray matter, CSF, and head motion parameters (Friston et al., 1996) as regressors of nuisance.
For the analysis of frontolimbic functional connectivity, we used a seed-based analysis approach. We defined three bilateral regions of interest (ROIs) using the Wake Forest University PickAtlas software (RRID:SCR_007378; http://www.fmri.wfubmc.edu/download.htm): amygdala, hippocampus, and the ventromedial prefrontal cortex (VMPFC), defined as the addition of Brodmann's areas 11, 25, and 47 (Ladouceur et al., 2011; Fig. 1). Using the REST toolbox (http://www.restfmri.net) for SPM8, we computed correlation coefficients between the times series for four pairs of ROIs (bilaterally amygdala/VMPFC; bilaterally hippocampus/VMPFC).
Illustration of the amygdala and VMPFC on an axial view of a single subject.
To explore the effect of rs6039769, we conducted multiple linear regressions using age, sex, diagnosis, and lateralization as covariates, with hippocampal volumes (left + right/2), amygdala volumes, and correlation coefficients as dependent variables in independent analyses using R 2.14.1. The rs6039769 risk variant was modeled in a model comparing AA and AC genotype versus CC as the AA group was relatively small. For analyses of volumes, we used the total intracranial volume as an additional covariate. Total intracranial volumes were computed as the sum of gray matter, white matter, and CSF volumes obtained after segmentation of the T1 images with the “new segment” procedure implemented in SPM8 software. We checked for the normality of residuals using Q–Q plots. For significant results, we additionally checked for an interaction between genotype and gender, as such interactions have previously been described on emotional neuroimaging variables; when significant, we performed post hoc stratified analyses by gender. We also tested for an interaction between genotype and status (control/BD) and could not find any interaction.
Results
The SNAP25 promoter variant rs6039769 is associated with SZ
We genotyped rs6039769 in 461 patients with SZ and 315 unaffected control subjects. Both populations were in Hardy–Weinberg equilibrium (patients: χSZ2 = 1.81, 2 df, pSZ = 0.18; control subjects: χctrl2 = 0.83, 2 df, pctrl = 0.36). We showed a significant increase in the C allele frequency in patients compared with control subjects [χ2 = 5.33; 1 df; p = 0.02; odds ratio (OR) = 1.29; 95% CI = 1.04–1.60; Table 1). The genotypic association study also showed a significant difference in genotypic distribution between patients and control subjects (χ2 = 8.03, 2 df, p = 0.02). The most significant difference was observed for CC genotypes under a recessive model (χ2 = 7.76, 1 df, p = 0.006). As more males were present in the patient group than in the control subject group (χ2 = 14.11, 1 df, p = 0.0002), we checked for gender-specific difference in genotype frequencies. However, no difference was observed between males and females, either in patients (χ2 = 0.19, 1 df, p = 0.91) or in control subjects (χ2 = 0.76, 1 df, p = 0.68).
rs6039769 genotype and allele frequency in the population of patients with SZ and unaffected control subjects
As a replication sample, we used results from the latest Psychiatric Genomics Consortium GWAS on schizophrenia (Schizophrenia Working Group of the Psychiatric Genomics, 2014). Similar to what we observed in the discovery sample, a significant increase in allele frequency was observed for the rs6039769 C at-risk allele in patients with SZ when compared with control subjects (χ2 = 5.03; 1 df; p = 0.02; OR = 1.03; 95% CI = 1.01–1.04).
The rs6039769 C allele is sufficient to increase the SNAP25 transcription level in vitro
We previously showed that rs6039769 was associated with an increased expression level of SNAP25 in prefrontal cortex (Etain et al., 2010). To determine whether this effect was due to this rs6039769 or another SNP in linkage disequilibrium, we cloned 784 bp of the SNAP25 promoter region upstream to the luciferase reporter gene in the pGL3-basic vector. Plasmids with either the C or the A alleles were thus transfected in mouse neuroblastoma cells. Both plasmids showed a significant increase in the expression level when compared with the empty vector (Fig. 2). The presence of the at-risk C allele resulted in a significant increase of the luciferase activity of almost 20% (Student's t test, t(40) = 4.76, p = 7.7 × 10−5).
In vitro analysis of rs6039769 in N2a cells. Quantification of the SNAP25 expression level was estimated using a luciferase reporter gene under control of the human SNAP25 promoter region containing either the C or the A allele of rs6039769. Both constructs showed a higher reporter expression level than the empty vector. The C allele showed a 20% increase of the reporter expression level than the A allele. The values are reported as the mean ± SD. ***p < 0.0001.
The SNAP25b/SNAP25a ratio is increased in schizophrenic patients carrying the rs6039769 C at-risk allele
Using postmortem studies, the expression level of SNAP25 has been previously reported to be modified in patients with SZ (Thompson et al., 1998, 2003a; Fatemi et al., 2001; Honer et al., 2002; Scarr et al., 2006; Gray et al., 2010). To determine whether such modification in brains of patients was dependent on rs6039769, we analyzed the transcription level of SNAP25 using quantitative real-time PCR on postmortem prefrontal cortex from patients with SZ and control subjects. Two isoforms exist for SNAP25 (SNAP25a and SNAP25b), which differ only by 9 aa (Bark et al., 1995), but have different function and expression pattern (Sørensen et al., 2003; Bark et al., 2004; Yamamori et al., 2011). We thus analyzed both isoforms independently as well as their combination. No significant difference was observed according to the rs6039769 genotype for SNAP25a or for SNAP25b, either in patients with SZ or in control subjects (Fig. 3A–D). While no difference in the SNAP25b/SNAP25a ratio was observed according to rs6039769 genotype in control subjects (Kruskal–Wallis test, K = 0.59, p = 0.75; Fig. 3E), we showed a significant increase of this ratio when patients had the at-risk allele C (Kruskal–Wallis test, K = 6.23, p = 0.04; Fig. 3F). We checked that no significant interaction was found between the SNAP25b:SNAP25a ratio and postmortem interval, refrigerator interval, brain pH, and gender (data not shown).
A–F, SNAP25a and SNAP25b mRNA relative expression level in postmortem prefrontal cortex of patients with schizophrenia (B, D, F) and healthy control subjects (A, C, E) according to rs6039769 genotype. Data are expressed as a mean value of relative mRNA expression level of SNAP25a (A, B), SNAP25b (C, D), and SNAP25b/SNAP25a ratio (E, F). No significant difference was observed for patients and control subjects for either SNAP25a or SNAP25b. A significant difference was observed for SNAP25b/SNAP25a ratio only for patients with schizophrenia according to the rs6039769 genotype, with an increased ratio for individuals homozygous for the C allele (Kruskal–Wallis test, K = 6.23, p = 0.04). Subjects homozygous for the rs6039769 C at-risk allele are shown with a black symbol, whereas others genotypes are shown with a clear symbol. Individual results and median values are reported. NS, Not significant.
SNAP25 rs6039769 C at-risk allele and the prefrontal–limbic network in humans
We combined genetic and brain imaging data in two independent cohorts to determine the in vivo functional consequences of the SNAP25 at-risk allele. No difference in genotype frequencies was observed between the discovery and replication cohorts (χ2 = 1.39, 2 df, pdisc = 0.50) and both were in Hardy–Weinberg equilibrium (discovery cohort: χdisc2 = 1.96, 2 df, pdisc = 0.16; replication cohort: χrep2 = 0.03, 2 df, prep = 0.86; Table 2). Our analyses revealed a gender-dependent significant correlation between rs6039769 genotype and the amygdala volume (multiple linear regression genotype × gender interaction, t = −2.8; p = 0.008; Fig. 4A). Both in patients with BD and healthy control subjects, males with an at-risk genotype (CC) had a larger amygdala than non-risk carriers (AA or AC; Student's t test, t(31) = 2.73, p = 0.01), but this difference was not observed in females (Student's t test, t(33) = 1.31, p = 0.20). We did not find any effect on the hippocampal volume. In the “replication” sample, we also detected a significant association between the risk genotype and a larger amygdala volume, localized on the left side (ANCOVA, F(1,117) = 8.05; p = 0.005).
rs6039769 genotype and allele frequency in the two populations used for imaging genetic studies
Structural and functional prefrontal limbic network characterization in regard to SNAP25 genotype. A, Amygdala volume in patients with schizophrenia and unaffected control subjects. Males homozygous for the SNAP25 rs6039769 C at-risk allele have a larger amygdala than non-risk carriers (AA or AC) or than females regardless of their genotype. B, Amygdala–VMPFC connectivity in patients with schizophrenia and unaffected control subjects. Again, males homozygous for the SNAP25 rs6039769 C at-risk allele have a higher connectivity than non-risk carriers (AA or AC) or than females regardless of their genotype. Squares show males, and circles show females. Black symbols show individuals with CC genotype and clear symbols show those with AA or AC genotype. Individual results are reported with mean values ± SD.
For prefrontal–limbic functional connectivity, we identified a gender-dependent significant relationship between rs6039769 genotype and functional connectivity measures between the amygdala and the VMPFC (multiple linear regression genotype × gender interaction, t = −2.3, p = 0.03; Fig. 4B). In both patients with BD and healthy control subjects, male risk carriers had a greater correlation coefficient than non-risk carriers or than females regardless of their genotype.
Discussion
SNAP25 genetic variants are associated with an increased risk of multiple psychiatric disorders
In this study, we showed that the rs6039769 SNAP25 at-risk allele for EOBD was also associated with SZ. We then explored the functional impact of this variant both in vitro and in vivo using postmortem studies and brain imaging. We showed that this variation was associated with an increased transcription level, which was consistent with the increased SNAP25b/SNAP25a ratio observed in patients with SZ. This was also consistent with a larger amygdala volume and an increased amygdala–VMPFC connectivity reported in men carrying the CC at-risk genotype.
The association between rs6039769 and SZ is consistent with the abundant literature identifying SNAP25 genetic variations as risk factors for many psychiatric conditions, such as SZ (Wang et al., 2015), BD (Etain et al., 2010) or attention-deficit/hyperactivity disorder (ADHD; Liu et al., 2017), and strongly supports the hypothesis of a shared genetic vulnerability that has been suggested for these mental disorders (Lichtenstein et al., 2009; Lee et al., 2013). The association that we previously reported between rs6039769 and BD was specific to early-onset forms of the disorder, which is highly comorbid with ADHD (West et al., 1995; Wozniak et al., 1995). This common vulnerability allele may thus be associated with a specific comorbid condition shared among several psychiatric disorders, such as psychotic symptoms that are more frequently observed in early-onset forms of BD (Leboyer et al., 2005).
Patients with psychiatric disorders exhibit alterations in the expression level of SNAP25
The functional consequences of this promoter variation help us in understanding how it may contribute to a higher liability for SZ and EOBD. The increased transcription level we observed is consistent with the increased SNAP25 levels reported in CSF of patients with SZ (Thompson et al., 1999, 2003b) and with altered SNAP25 levels in postmortem brain samples from patients with SZ or BD. These reports have highlighted a decreased level of SNAP25 in the hippocampus (Thompson et al., 2003a) but increased densities in prefrontal cortex area 9 in patients with SZ or BD (Thompson et al., 1998; Scarr et al., 2006). Animal models of SZ, such as mutant human DISC1 mice, also exhibit altered SNAP25 expression (Pletnikov et al., 2008), and mutant mice for SNAP25 exhibit several phenotypic traits similar to those present in psychosis, such as impaired sensorimotor gating (Jeans et al., 2007).
The SNARE complex, including SNAP25, is crucial for vesicular membrane fusion and neurotransmitter release in synapses (Graham et al., 2002). SNAP25 is also necessary for axonal growth (Osen-Sand et al., 1993; Shimojo et al., 2015), dendrite formation (Kimura et al., 2003), and calcium homeostasis (Verderio et al., 2004; Pozzi et al., 2008). Therefore, an altered expression level of SNAP25 may impact the brain neurotransmission. More specifically, we observed an altered balance in the SNAP25b/SNAP25a isoforms ratio, which probably results from the SNAP25 expression increase in subjects with rs6039769 CC genotype. The two isoforms are differentially expressed during the development, with the SNAP25a being associated with development and synaptogenesis, while SNAP25b is the major form in mature brains (Prescott and Chamberlain, 2011; Yamamori et al., 2011). Note that during postnatal development both isoforms have a similar expression level and SNAP25b is dramatically upregulated in early adulthood (Bark et al., 1995), a critical period for SZ and BD onset. A modified SNAP25b/SNAP25a balance may alter the brain maturation or the balance among neurite outgrowth, synaptogenesis, and neurotransmitter release (Yamamori et al., 2011). This increase of SNAP25b compared with SNAP25a may modify the stability of vesicles in the releasable vesicle pool (Sørensen et al., 2003) and thus increase neurotransmitter exocytosis.
Synaptic transmission is impaired in SZ and BD and is also the target of their common medication (i.e., antipsychotic agents). Haloperidol has been shown to increase SNAP25 levels in the striatum (Barakauskas et al., 2010). The mouse mutant coloboma, which bears a mutation on chromosome 2 including SNAP25, displays a decreased glutamate level in the neocortex and a low dopamine level in the ventral striatum (Raber et al., 1997). In another mouse model, carrying a single amino acid change in SNAP25, Kataoka et al. (2011) reported abnormal serotonin and dopamine release in the amygdala, along with anxiety-related behavior. This rodent model exemplifies the potential impact of SNAP25 on serotonin, amygdala, and, in turn, behaviors. Though we did not have human behavioral dimensional measures in this study, a previous study has reported a link between a SNAP25 genotype and emotional regulation in a population of children between the ages of 6 and 20 months (Sheese et al., 2009). Another group identified an association between neuroticism and SNAP25 in a GWAS that included 3972 individuals (Terracciano et al., 2010).
SNAP25 genetic variants are associated with structural and functional brain alterations
Consistent with what was observed in animal models, we found, through our imaging genetics approach in humans, that the rs6039769 CC genotype was associated with a larger amygdala volume and an increased functional connectivity between VMPFC and amygdala. The amygdala is a core component of emotion production and emotional processing, mostly implicated in automatic emotional regulation (Phillips et al., 2008). Alterations in amygdala volume, function, and/or connectivity have been repeatedly reported in individuals with BD and SZ (Hajek et al., 2009; van Erp et al., 2016). Amygdala volume has been found to be enlarged in at-risk offspring of patients with BD (Bauer et al., 2014) and in adolescents at risk for SZ (Welch et al., 2010; Ganzola et al., 2014). Intriguingly, a meta-analysis (Hajek et al., 2009) identified a decreased volume in children and adolescents with BD compared with control subjects and a trend toward a larger left amygdala volume in adult patients with BD. In our study, we only studied adult subjects in both of our samples (discovery and replication) and identified an increase in the amygdala volume associated with the risk variant. In our replication sample, our effect was also left sided. One possible interpretation is that both BD and the SNAP25 risk variant are associated with an abnormal neurodevelopment of the amygdala, leading to a smaller amygdala in children and adolescents and a larger volume in adults. This is speculative and should be tested in children and adolescents or through longitudinal studies. A higher SNAP25 expression in patients with psychiatric disorders may lead to a chronic hyperactivity, which has been proposed as a putative mechanism for an enlargement of amygdala volumes, resulting in emotional and social dysregulation and increasing the risk for BD, SZ, and ADHD (Rosenfeld et al., 2011). Prefrontal limbic connectivity is also a key network implicated in emotion regulation, allowing a prefrontal modulation over the amygdala activation (Phillips et al., 2008). An increased functional connectivity between VMPFC and amygdala is thus a crucial feature of BD (d'Albis and Houenou, 2015). This altered connectivity is also present in healthy relatives of patients with BD, suggesting it as a possible endophenotype, a marker of a genetic liability to BD. This altered prefrontal–limbic connectivity may explain the inability of patients with BD to efficiently regulate emotional states, leading to emotional hyperreactivity (Henry et al., 2012). Altered prefrontal–amygdala functional connectivity has also been reported to be altered in patients with SZ (Vai et al., 2015). Structural connectivity of the uncinate fasciculus that connects VMPFC and amygdala has consistently been found altered in first-episode patients with SZ (Wheeler and Voineskos, 2014) and in healthy relatives (Muñoz Maniega et al., 2008). Prefrontal–limbic dysconnectivity is assumed to be one of the neural networks abnormalities common to SZ and BD (d'Albis and Houenou, 2015). Overall, our results suggest that the increased liability to SZ and EOBD associated with the rs6039769 risk allele may be mediated by its action on limbic networks (amygdala volume and prefrontal–limbic connectivity). Our associations between the SNP variation and neuroimaging variables were only present in male, suggesting a gender-dependent impact of SNAP25 on limbic networks. Gender differences in emotion processing are largely described in the literature (Nolen-Hoeksema, 2012; Lungu et al., 2015). In parallel, the regional SNAP25 density varies according to gender with a higher density in men than in women in frontal regions, with an inverse relationship in temporal areas (Downes et al., 2008). In addition, hormones have been found to regulate the expression of SNAP25 in animal models (Jacobsson et al., 1998). One study (Cagliani et al., 2012) also identified a gender-specific correlation between a SNAP25 allelic variation and verbal performance. All these studies, along with the fact that SNAP25 has been associated with several male-predominant psychiatric conditions (ADHD, SZ, and autism spectrum disorder) shed light on our gender-specific finding. However, we cannot fully explain why this result was gender specific while the other ones were not. Note that all of our samples were enriched in males and may thus not be fully appropriate to observe significant differences between genders.
We selected different populations for the samples of the imaging genetics study: patients and control subjects for the discovery sample, and control subjects only for the replication group. Selecting the better population for imaging genetics is highly debated in the literature (Bigos and Weinberger, 2010). Further, the replication sample was significantly younger than the discovery sample. However, we found similar results in those two different groups. These points strengthen the validity and generalizability of our results but may also explain why we did not identify gender effects in the replication sample.
Despite bringing evidence from multiple levels, we recognize some limitations in our study. First, the sample size for the MRI study is relatively modest, and these results should be taken with caution. We thus did not correct our analyses for multiple comparisons. However, we were able to replicate the association between the risk genotype and a larger amygdala volume in an independent sample. Independent samples with a larger sample size would also allow conducting whole-brain structural and connectivity analyses. Second, the sample size of the genetic study was also modest, but the same genetic effect was observed in the Psychiatric Genomics Consortium cohorts. Third, the limited number of available brain samples prevents any replication of our data and more specifically of the importance of the SNAP25b/SNAP25a ratio. Because of the crucial importance of the amygdala and hippocampus in SZ and BD, studying SNAP25 expression levels in these regions depending on the genotype would be of high interest. However, alteration of the SNAP25 level and its binding partners of the SNARE complex have been consistently reported in humans with psychiatric disorders as well as in mouse models, demonstrating the importance of the regulation of this complex. Further, study of animal models modulating the ratio between the two SNAP25 isoforms and its interaction with environmental factors, such as pre-perinatal stress, should help in understanding the physiological impact of such deregulation and its impact on behaviors.
In conclusion, through our multimodal approach, from in vitro construct to imaging genetics and postmortem brain expression, we brought strong evidence for the functional consequences of this allelic variation of SNAP25. Our findings suggest that rs6039769, as a functional polymorphism located in the promoter region of SNAP25, may modulate the development and plasticity of the prefrontal–limbic network and therefore increase vulnerability to EOBD and SZ.
Footnotes
This work was supported by the Investissements d'Avenir program managed by the Agence Nationale pour la Recherche (ANR) under reference ANR-11-IDEX-0004-02 (Labex BioPsy). This work also received financial support from the ANR (Grant ANR-08-MNPS-041), the Institut National Pour La Sante et la Recherche Medicale, the Réseau Thématique de Recherche et de Soins en Santé Mentale (Fondation FondaMental), and the Institut Thématique Multi Organismes Neurosciences, Sciences Cognitives, Neurologie et Psychiatrie. The Jamain group is affiliated with the Paris School of Neuroscience and the Bio-Psy Laboratory of Excellence. We thank the patients and control subjects for their participation. We also thank Pr. W. Song for SNAP25 promoter constructs. In addition, we thank E. Abadie, the Cochin Hospital Cell Library (Professor J. Chelly), the Clinical Investigation Centre (Dr. P. Lecorvoisier). and the Biological Resources Platform (Professor B. Ghaleh) of Mondor Hospital for technical assistance. Finally, We thank the personnel of the UNIACT, NeuroSpin Platform for their help in the MRI scanning of patients and control subjects.
The authors declare no competing financial interests.
- Correspondence should be addressed to Dr. Stéphane Jamain, INSERM U955, Psychiatrie Translationnelle, Faculté de Médecine, 8 rue du Général Sarrail, Créteil, 94000, France. stephane.jamain{at}inserm.fr
