Skip to main content

SNAP25 is a potential target for early stage Alzheimer’s disease and Parkinson’s disease

Abstract

Background

Alzheimer’s disease (AD) and Parkinson’s disease (PD), two common irreversible neurodegenerative diseases, share similar early stage syndromes, such as olfaction dysfunction. Yet, the potential comorbidity mechanism of AD and PD was not fully elucidated.

Methods

The gene expression profiles of GSE5281 and GSE8397 were downloaded from the Gene Expression Omnibus (GEO) database. We utilized a series of bioinformatics analyses to screen the overlapped differentially expressed genes (DEGs). The hub genes were further identified by the plugin CytoHubba of Cytoscape and validated in the hippocampus (HIP) samples of APP/PS-1 transgenic mice and the substantial nigra (SN) samples of A53T transgenic mice by real-time quantitative polymerase chain reaction (RT-qPCR). Meanwhile, the expression of the target genes in the olfactory epithelium/bulb was detected by RT-qPCR. Finally, molecular docking was used to screen potential compounds for the target gene.

Results

One hundred seventy-four overlapped DEGs were identified in AD and PD. Five of the top ten enrichment pathways mainly focused on the synapse. Five hub genes were identified and further validated. As a common factor in AD and PD, the changes of synaptosomal-associated protein 25 (SNAP25) mRNA in olfactory epithelium/bulb were significantly decreased and had a strong association with those in the HIP and SN samples. Pazopanib was the optimal compound targeting SNAP25, with a binding energy of − 9.2 kcal/mol.

Conclusions

Our results provided a theoretical basis for understanding the comorbidity mechanism of AD and PD and highlighted that SNAP25 in the olfactory epithelium may serve as a potential target for early detection and intervention in both AD and PD.

Introduction

Neurodegenerative diseases (NDDs) are common irreversible diseases, which are the leading causes of disability in people over 65 years [1, 2]. NDDs mainly include Alzheimer’s disease (AD), Parkinson’s disease (PD), amyotrophic lateral sclerosis (ALS), Huntington’s disease (HD), and multiple sclerosis (MS). Among them, AD and PD have higher incidences and are usually sporadic in most cases. The characteristic pathological hallmarks of AD and PD are progressive loss of neurons and accumulation of pathological proteins such as extracellular amyloid-β (Aβ) and intraneuronal neurofibrillary tangles (NFTs) in AD and α-synuclein (α-Syn) and Lewy bodies/Lewy neurites (LBs/LNs) in PD. Accumulating evidence indicates that the pathologies of AD and PD begin in the asymptomatic stage, and there are no effective therapies available to prevent or delay the progressions of AD and PD [3, 4]. Meanwhile, several studies have found that LBs are deposited in the brains of AD patients [5, 6]. Studies also proved that α-Syn was the major non-Aβ component of Aβ plaques, while NFTs presented in vulnerable brain regions of PD patients [7]. The co-occurrence of different pathologies of the two diseases suggests that AD and PD have a common molecular basis, which makes the early diagnosis and treatments of AD and PD more challenging. However, the potential common molecular basis of AD and PD remains unclear.

Besides the common pathological characteristics, similar clinical symptoms, such as cognitive disorder and olfactory dysfunction, also coexist in AD and PD patients [8, 9]. More than 90% of AD and PD patients emerged with olfactory dysfunction in the early stages. Accumulating evidence demonstrated that olfactory dysfunction occurred at Braak stages 1 and 2 of sporadic PD, with damage to olfactory structures and pathological aggregation of α-syn in the olfactory bulb in the early stage, and then spread to the brain stem and related brain regions [10,11,12,13]. Similar deposition patterns of pathological proteins can be found in AD [14, 15]. More importantly, many studies have demonstrated that olfactory dysfunction is related to cognitive impairment and preceded cognitive decline, which can be used to predict the cognitive decline of AD or PD patients [16, 17]. The olfactory system consists of the olfactory sensory neurons, the olfactory bulb (OB), and other multiple olfactory cortices. Olfactory sensory neurons in the olfactory epithelium (OE) of the nasal cavity, receive chemical information from external sources and convert the chemical signals into neural signals to the OB. The olfactory system is a complex neural network characterized by long axonal projection, with nerve fibers having neural projection connections directly or indirectly to multiple brain regions via synapse connections [18, 19]. Such structural features facilitate the diffusion of α-Syn and Aβ in the olfactory bulb to the susceptible brain regions and cause neuronal degeneration. The OB had structural–functional connections with multiple brain regions, such as the hippocampus (HIP) [14, 20,21,22]. The HIP is a crucial brain region charging for learning, memory, and cognition. The HIP is severely damaged in AD patients, as well as the substantial nigra (SN) in PD patients. Olfactory dysfunction is characterized by not only an increase in the olfactory threshold but also a decrease in olfactory discrimination, especially in early stage AD and PD [23, 24]. Several studies have found that pathological proteins accumulated in the OB disturbed synaptic proteins’ synthesis and induced synaptic functional dysfunction and cognitive impairment [25]. However, the common mechanism of olfactory dysfunction in AD and PD has not been fully explained yet.

Transcriptomics sequencing integrated with bioinformatics analysis has been widely applied to assay the alteration of gene expression levels and predict their possible implications in given diseases. A large amount of data are available in the authoritative public GEO database (http://www.ncbi.nlm.nih.gov/geo/). In this paper, we utilized two original gene profiles to explore the common molecular changes in AD and PD. Through a series of bioinformatics analyses, we found that the common hub genes in AD and PD mainly focused on synapse dysfunction. The mRNA expression levels of the hub genes were validated by real-time quantitative polymerase chain reaction (RT-qPCR). The odor detection and discrimination tests were used to assess the olfactory function of mice. By further molecular docking, we predicted the potential compounds for the target gene.

The present study aimed to identify the common alternations of molecules and molecular networks in early stage AD and PD and try to screen the potential targets.

Results

Identification of the overlapping DEGs

The gene expression profiles of GSE5281 (AD) and GSE8397 (PD), were used to identify overlapping DEGs in AD and PD brain samples. twofold change and P < 0.05 were as set as the cutoff criteria. A total of 2159 genes (915 up-regulated genes and 1244 down-regulated genes) were extracted from the GSE5281 data set (Fig. 1A), and 334 genes (76 up-regulated genes and 258 down-regulated genes) from the GSE8397 data set (Fig. 1B). We compared the up-regulated and down-regulated DEGs of GSE5281 and GSE8397, respectively, and finally identified 28 overlapping up-regulated DEGs and 146 overlapping down-regulated DEGs as shown in the Venn diagram (Fig. 1C, D). The 174 overlapping DEGs were used in the subsequent analysis.

Fig. 1
figure 1

Gene expression profiles of GSE5281 and GSE8397 databases. A Volcano plot of gene expression profile data in GSE5281. B Volcano plot of gene expression profile data in GSE8397. Red dots (significantly up-regulated DEGs), blue dots (significantly down-regulated DEGs), and Grey dots (no significant changes of genes). |logFC| > 1 and adjusted P value < 0.05 were as cutoff threshold. C Venn diagram of the up-regulated overlapping DEGs. D Venn diagram of the down-regulated overlapping DEGs

Enrichment analysis of the overlapping DEGs

GO and KEGG enrichment terms of the overlapping DEGs were obtained with a P value < 0.01, using Metascape and the plugin ClueGo of Cytoscape. The top 10 enrichment terms were selected and shown in Fig. 2A. Notably, five terms out of the top 10 terms including synaptic signaling, synapse organization, transmission across chemical synapses, synaptic vesicle endocytosis, and synaptic vesicle cycle in the synapse were relative to synaptic structure and function. As shown in Fig. 2B, vesicle-mediated transport in synapses, spontaneous synaptic transmission, synaptic signaling, post-synapse organization, synapse maturation, and structural constituent of synapse had high correlations with synaptic structure and function. Taken together, these results revealed that overlapping DEGs mainly focused on synaptic dysfunction in the occurrence of AD and PD.

Fig. 2
figure 2

Enrichment analyses of the overlapping DEGs in GSE5281 and GSE8397 databases. A GO enrichment analysis of the overlapping DEGs. The X-axis represents the q value (− log10), and the Y-axis represents the GO terms. The GO terms were measured by the rich factor, q value, and number of genes enriched. The greater the Rich factor is, the greater the degree of enrichment and the greater the P value. The brighter the color of red is, the more significant the term. B ClueGO analysis of the overlapping DEGs from AD/PD patients. Functionally grouped network with terms as nodes linked based on their kappa score level (≥ 0.4), where only the label of the most significant term per group is shown. The node size represents the term enrichment significance. Functionally related groups partially overlap

Construction of the PPI network and selection of hub genes

The overlapping DEGs were used to construct a PPI network via the String online and visualized via Cytoscape (Fig. 3A). The top two clusters with the highest clustering scores were obtained via MCODE. The first cluster had 6000 scores, 14 nodes, and 39 edges (Fig. 3B), including SNAP91, DNM3, DNM1, STXBP1, NEFL, AMPH, SYNJ1, GRIA1, GABRG2, ATP2B, CPNE6, CCK, SNCB, and ACTL6B. The second one had 5733 scores, 16 nodes, and 43 edges (Fig. 3C). Then, the PPI network was re-analyzed with the CytoHubba plugin to select the top 10 genes according to the MCC, MNC, and degree topological algorithm, respectively (Table 1). Finally, we obtained five overlapping hub genes after three Algorithms’ Comparison, including synaptosomal-associated protein 25(SNAP25), synapsin 1 (SYN1), synaptotagmin1 (SYT1), growth-associated protein 43 (GAP43) and synaptosomal-associated protein 91 (SNAP91) (Table 1). Collectively, the five hub genes encoded all presynaptic proteins and were both down-regulated in AD and PD patients.

Fig. 3
figure 3

PPI network and modular analysis of the overlapping DEGs. A PPI network was constructed by a total of 174 DEGs via the STRING database, and visualized by Cytoscape. The nodes represent proteins, the edges represent the interaction of the proteins. B, C Top two PPI clusters in MCODE analysis. degree cutoff = 2, node score cutoff = 0.2, k-core = 2 and max. Depth = 100

Table 1 Top 10 genes screened by three algorithms in CytoHubba

Validation of representative hub genes by real-time qPCR in 6-month AD and PD mice

We used RT-qPCR to validate the expression of the five hub genes in 6-month AD and PD mice and their age-matched control mice separately. As shown in Fig. 4A–E, we found that the expression level of all the selected hub genes was more significantly decreased in AD+ mice than in AD− mice (SNAP25, 0.74 ± 0.15, P = 0.0005; SYN1, 0.70 ± 0.23, P = 0.0037; SYT1, 0.67 ± 0.30, P = 0.0129; GAP43, 0.68 ± 0.27, P = 0.0079; SNAP91, 0.80 ± 0.16, P = 0.0064). Similarly, we also found that the expression level of all the hub genes was more significantly decreased in PD+ mice (SNAP25, 0.56 ± 0.25, P = 0.0004; SYN1, 0.85 ± 0.10, P = 0.0010; SYT1, 0.71 ± 0.28, P = 0.0166; GAP43, 0.62 ± 0.25, P = 0.0010; SNAP91, 0.68 ± 0.33, P = 0.0210, Fig. 4F–J). To explore the significance of the hub genes, we used ANOVA to detect hub genes separately. We found that the expression of SNAP25 and SNAP91 were most significant in the HIP in AD mice, while SNAP25 and GAP43 in the SN with PD mice (AD: PSNAP25 = 4.11E−03, PGAP43 = 2.18E−13; PD: PSNAP25 = 0.017, PSNAP91 = 0.031). All above, the results revealed those synapse-relative hub genes had changed in 6-month AD and PD mice. SNAP25, as a common hub gene in AD and PD, has the most significant changes.

Fig. 4
figure 4

Validation of five hub genes in 6-month AD and PD mice samples. The mRNA levels of five hub genes were validated in 6-month AD+/PD+ mice samples and age-matched control mice samples by real-time qPCR. Five hub genes include VAMP2, SNAP25, SYT1, SYN1, NSF, and STX1A. All hub genes significantly declined in both AD+/PD+ mice samples, compared to those in control samples (AD− and PD−). Relative quantification (RQ) values were calculated. Data are the mean ± s.e.m from 8 mice in each group. A T test was used to analyze the difference between the two groups. *0.01 ≤ P < 0.05, **0.001 ≤ P < 0.01, ***P < 0.001

Olfactory deficit in the 6-month AD and PD mice

The olfactory dysfunction of mice was evaluated using odor detection and odor discrimination tests. The procedures are shown in Fig. 5A. Previous studies had demonstrated that AD or PD mice at 6-month-old showed abnormal aggregation of pathological proteins and atypical clinical symptoms [26,27,28]. Thus, the 6-month-old mice in the early stages of AD and PD were used in these olfactory function tests. It is well-known that the mice exhibit habituation for the same odor, while dishabituation for the novel odor. During the odor detection test, AD− mice exhibited normal habituation for the same odor and an increase in exploring time for novel odor (P < 0.05, Fig. 5B). However, the exploring time of AD+ mice showed no obvious change for the same or novel odors (P > 0.05, Fig. 5C). Similar results were found in 6-month PD− and PD+ mice (Fig. 5F, G). During the odor discrimination test, AD− mice showed significant habituation for the same odor, and an obvious dishabituation for novel odor (P < 0.05, Fig. 5D). Conversely, AD+ mice did not display the habituation for the same odor nor dishabituation for similar and novel odor (P > 0.05, Fig. 5E). A similar pattern was also found in PD− and PD+ mice (Fig. 5H, I). In addition, the cross-habituation indexes of AD and PD mice were calculated based on the odor discrimination tests. Compared to AD− mice, the cross-habituation index of AD+ mice robustly declined, as well as that in PD+ mice (all P < 0.01, Fig. 5J, K).

Fig. 5
figure 5

6-month AD and PD mice showed olfactory defects. A Behavioral diagram of odor detection and odor discrimination. B, C Exploring times of AD+ mice were decreased than those of AD- mice in the odor detection task. D, E Exploring time of AD+ and AD− mice in odor discrimination task. F, G Exploring times of AD+ mice were decreased than that of AD− mice in the odor detection task. H, I Exploring time of AD+ and AD− mice in odor discrimination task. J, K Discrimination index of AD and PD mice declined. Data are the mean ± s.e.m from 8 mice in each group. A T test was used to analyze the difference between the two groups. *0.01 ≤ P < 0.05, **0.001 ≤ P < 0.01, ***P < 0.001

Correlation analysis of the SNAP25 mRNA between the OB and OE samples and the HIP or SN samples

To define the changes of SNAP25 in the OB and OE of AD and PD mice, we used real time-qPCR to validate the expression levels of SNAP25 in 6-month AD+ and PD+ mice and age-matched AD− and PD− mice separately. As shown in Fig. 6A, B, we found that the expression of SNAP25 decreased both in the OB and OE of AD+ mice compared to AD− mice (in OB, 0.30 ± 0.21, P < 0.0001; in OE, 0.46 ± 0.28, P = 0.0001). Similarly, we found that the expression of SNAP25 in the OB and OE significantly decreased in PD+ mice compared to PD− mice (in OB, 0.44 ± 0.20, P = 0.0001; in OE, 0.45 ± 0.28, P = 0.0002, Fig. 6C, D). In addition, we correlated the RQ values of SNAP25 in the HIP/SN with those in the OB and OE of AD and PD mice. As shown in Fig. 6E, F, they exhibited a high correlation (HIP vs OB, r = 0.7734, P = 0.0244; HIP vs OE, r = 0.8099, P = 0.0148). Similar results were found in PD mice (SN vs OB, r = 0.7979, P = 0.0176; SN vs OE, r = 0.7945, P = 0.0185, Fig. 6G, H). Our results showed that the expressions of SNAP25 in the OB and OE decreased in the early stages in AD and PD mice and were highly correlated with those in representative brain areas of AD and PD mice.

Fig. 6
figure 6

Correlation of expression of SNAP25 between the OB/OE and HIP/SN samples. A, B Real-time qPCR displayed that the expressions of SNAP25 in the OB and OE were significantly declined in 6-month AD+ mice samples than those in AD- mice samples. C, D Real-time qPCR revealed the expressions of SNAP25 in the OB and OE decreased in 6-month PD+ mice samples than those in PD- mice samples. E, F Pearson correlative analysis revealed that SNAP25’s expressions in the OE and OB have a high correlation with those in the HIP in AD+ mice. G, H Pearson correlative analysis revealed that SNAP25’s expressions in the OE and OB have a high correlation with those in the HIP in PD+ mice. Expression of SNAP25 was normalized by β-actin. Relative quantification (RQ) values were calculated. Data are the mean ± s.e.m from 8 mice in each group. A T test was used to analyze the difference between the two groups. *0.01 ≤ P < 0.05, **0.001 ≤ P < 0.01, ***P < 0.001. Pearson correlative analysis was used to analyze the correlation of SNAP25 expression between the OB/OE and the HIP/SN

Virtual screening of potential compounds targeting SNAP25

SNAP25 is a part of the soluble N-ethylmaleimide-sensitive fusion protein attachment protein receptor (SNARE) complex which is critical for synaptic vesicle fusion in the presynaptic membrane, a key step in the release of neurotransmitters. SNAP25 is composed of 206 amino acid residues. When synaptic vesicles are close enough to the presynaptic membrane, VAMP, syntaxin, and the cytoplasmic region of SNAP25 bind to each other to form a tight complex. In addition, SNAP25 interactions with SYT1, take part in the mechanism of membrane fusion initiated by Ca2+, and are involved in the fusion pore dynamics regulating the final steps of transmitter release [29]. SNARE functional interfaces are usually large and extended, making it difficult for small molecules to modify their conformation significantly. In the present study, we screened small-molecule drugs targeting SNAP25 via virtual screening. Based on the crystal structure of SNAP25, we first predict its potential functional pockets (Fig. 7A). On top of that, ten computationally plausible small-molecule drugs were captured according to their binding energies. The specific information on the residues of drug–protein interactions was listed (Table 2). Among them, Pazopanib with a binding energy of − 9.2 kcal/mol was the best-scoring drug. Moreover, Pazopanib was able to perfectly embed in the predicted activity pocket and interacted with SNAP25 protein residues at ALA18, HIS139, GLU143, ILE195, MET14, GLN15, ALA18, LYS189, ALA191, LEU192, and ARG198 via van der Waals interactions and hydrogen bonding (Fig. 7B, C). These results suggested that SNAP25 could be used as a potential early intervention target for AD and PD.

Fig. 7
figure 7

Virtual screening results between SNAP-25 and FDA-approved drugs. A Structure of the rat neuronal SNARE complex (PDB accession#: 3HD7[22] and binding pocket predicted by PROTEINS PLUS (https://proteins.plus/). B H-bond receptor surface between SNAP-25 and Pazopanib. C Residues of interaction between Pazopanib and SNAP-25. The other corrected information of virtual screening is shown in Table 2

Table 2 Docking score and predicted protein–ligand interaction of the top ten compounds selected in virtual screening

Discussion

In this study, we identified common gene expression profiles in early stage AD and PD via a publicly available GEO database and found five common core molecules associated with synapse in AD and PD data sets, including SNAP25, SYN1, SYT1, GAP43, and SNAP91. All these hub genes were downregulated in AD and PD patient samples. SNAP25 significantly decreased in the HIP of 6-month AD mice and the SN of 6-month PD mice, as well as in the OE and OB samples of both types of mice. Moreover, the decline in the OE and OB displayed a strong association with those in the HIP and SN samples. Through molecular mocking, we further predicted the potential compounds targeting SNAP25. These data suggest that SNAP25 may be a potential target for early intervention in AD and PD.

In the present study, we identified 174 overlapping DEGs (28 up-regulated and 146 down-regulated) as common DEGs in two data sets. In enrichment analysis against GO terms and KEGG pathways, we found that five of the top 10 GO terms were relative to synaptic structure and function (Fig. 2A). Similar results can also be found in ClueGo enrichment analysis (Fig. 2B). Our results suggest that synaptic dysfunction was involved in the early occurrence of AD and PD, consistent with previous studies showing that synaptic dysfunction occurred in AD and PD, even in the asymptomatic phase [30, 31]. Moreover, we obtained 5 common hub genes (SNAP25, SYN1, SYT1, GAP43, SNAP91) in AD and PD databases. All of them expressed presynaptic proteins and played crucial roles in synaptic formation and neurotransmitter release. Like those in AD and PD data sets, the mRNA level of five hub genes significantly declined in HIP samples of early stage AD mice and SN samples of PD mice (Fig. 4A–J). Among them, the change of SNAP25 was the most obvious. Previous studies found that the genetic polymorphism of SNAP25 was correlated with AD and PD progression [32, 33]. In AD and PD patients, the decline of SNAP25 levels in the cortex caused neuronal degeneration [34]. In AD patients with mild cognitive impairment, SNAP25 was elevated in cerebrospinal fluid and had a strong association with t-tau [35, 36].

Olfactory dysfunction is a common symptom of AD and PD [23, 37]. More than 90% of AD and PD patients had early onset olfactory dysfunction. Olfactory dysfunction precedes cognitive decline, which reflects the deterioration in the cognitive function of patients [24]. Several studies have found that pathological proteins accumulated in the OB induced synaptic dysfunction and cognitive impairment [25]. In the present study, we validated the expression of SNAP25 in several samples across multiple brain regions of AD and PD mice. Interestingly, SNAP25 decreased in the OB and OE samples, as well as in the HIP and SN samples (Figs. 3A, F; 5A–D). Moreover, the expressions of SNAP25 in the OB and OE were strongly correlated with those in the HIP of AD mice and the SN of PD mice samples (Fig. 5E–H). Several studies have demonstrated that SNAP25 in cerebrospinal fluid was a potential synaptic biomarker, and was associated with cognitive decline in AD and PD patients [35, 38, 39]. Considering that the OE is localized in the nasal cavity, it is convenient to collect test samples or pursue drug inventions in the early stages of AD and PD. In the present study, we predicted the potential compounds bonded with SNAP25 via virtual screening (Fig. 7 and Table 2). Thus, it suggests that SNAP25 in the OE and OB may be a potential drug target for the prevention of AD and PD progression.

Conclusion

A drawback of this paper is that even though we identified SNAP25 as a potential and available drug invention target for AD and PD, further validation of the effectiveness and sensitivity of SNAP25 as a potential target in the OE is still lacking. Meanwhile, the predicted activation center of SNAP25 was originally from the SNARE complex, due to the absence of a separate crystal structure for SNAP25. Thus, we are aware that SNAP25 as a potential target of AD and PD needs a large number of clinical samples and both in vivo and in vitro experiments for further validation.

In summary, we here highlighted the potential comorbidity mechanism of AD and PD and identified the synaptic protein SNAP25 as a potential common target for AD and PD. Furthermore, we also predicted the potential compounds targeting SNAP25 for possible early intervention in AD and PD.

Materials and methods

Animals

AβPP/PS-1 mice (AD) were purchased from the Beijing HuFuKang bioscience company (Beijing, China). A53T transgenic mice (PD) were purchased from Shanghai Model Organisms Center, Inc. Experiments were performed on 6-month AD+ and PD+ mice, and the corresponding control mice (AD− and PD−). The animals were housed on a 12/12 light/dark cycle with ambient temperature (21 ± 1 °C), humidity (50 ± 5%), and food and water ad libitum. All operations and procedures were approved by the animal committee of Xuzhou Medical University.

Data set selection and preprocessing

The gene expression profiles of GSE5281 and GSE8397 were downloaded from the Gene Expression Omnibus (GEO) database. There are 87 late-onset AD patients’ data and 74 age-matched control subjects’ data in the GSE5281 (GPL570) data set. GSE8397 (GPL96) data sets contain 30 PD patients’ data and 17 age-matched control patients’ data. We used GEO2R (https://www.ncbi.nlm.nih.gov/geo/geo2r/), an interactive online tool, to identify differentially expressed genes (DEGs). With default settings, the adjusted P values were selected to decrease the false positive rate using the Hochberg false discovery rate and the Benjamini method. A P value < 0.05 and an absolute log fold-change (FC) greater than 1 for the DEGs were used as the cutoff criteria. Up/down-regulated overlapping DEGs in the GSE5281 and GSE8397 data sets were identified by Venn Diagrams. The DEGs were visualized by volcano maps using OrginPro2021b.

Functional enrichment analysis of the overlapping DEGs

Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were performed to predict the biological functions of DEGs. In the present study, Metascape (http://metascape.org/) and the plugin ClueGo of Cytoscape were used to perform functional enrichment analysis of the overlapping DEGs with a P value < 0.01. The top 10 items were selected to show the possible functional role of the overlapping DEGs. In the calculated process of ClueGo, a kappa coefficient was calculated to reflect the functional correlations between paths or terms based on gene overlaps between pathways or GO terms. The kappa threshold defaulted to 0.4. Functionally similar entries were displayed in the same color. The threshold for enrichment significance was P < 0.05. The results of the enrichment analysis of the overlapping DEGs were visualized by bubble diagrams.

PPI network construction of the overlapping DEGs and Hub gene selection

The PPI network of the overlapping DEGs was constructed using the STRING online (http://string-db.org/), which predicted experimental interactions of proteins [40], with a confidence score ≥ 0.4 for significant differences. The PPI network was visualized by Cytoscape (URL). Molecular complex detection (MCODE), a plugin of Cytoscape, was used to obtain core modules of the PPI network. The screening threshold was set at 2 as the degree cutoff, 0.2 as the node score cutoff, 2 as the k-core, and 100 as the max depth. The top two clusters were obtained as final core modules. Furthermore, the PPI network was re-analyzed by the CytoHubba plugin of Cytoscape software (version 3.7.2). The top 10 hub genes were selected from three algorithms, including the maximal clique centrality (MCC), maximum neighborhood component (MNC), and degree. Finally, five overlapping hub genes were obtained as final hub genes.

RNA extraction and RT-qPCR

Fresh tissues of the HIP, SN (ventral dorsal midbrain), OB, and OE were obtained, respectively, from 6-month AD and PD mice and age-matched wild-type mice. Total RNA was extracted using a TRIzol reagent (KeyGen, Nanjing, China). Reverse transcription of total RNA and qPCR were performed according to the manufacturer’s recommendations of the reagent kit (MCE, New Jersey, USA). The primers are shown in Additional file 1: Table S1. All primers were used in the HIP and SN, and only the SNAP-25 primer was used in the OB and OE samples. Relative quantification (RQ) values were calculated using the 2−ΔΔCt method [41]. RQ values of the hub genes were calculated based on the ratio between RQ values of hub genes in AD and PD mice and normalization to β-actin in the control group.

Olfactory behavior test

The odor detection and odor discrimination tests were used to test the olfactory function of mice. Isoamyl acetate (diluted 1 × 10–6 in paraffin oil) was used in the odor detection test. Meanwhile, acetaldehyde, pentanol, and heptanone (diluted 1 × 10–3 in paraffin oil) were used in the odor discrimination test. Before the experiment started, the mice were habituated in the testing cage with a cotton applicator for 3 min. During the odor detection experiment, the mice received three successive trials of water followed by one trial of isoamyl acetate. Each trial lasted for 1.5 min, with an inter-trial interval of 30 s. During the odor discrimination experiment, the mice were sequentially exposed to acetaldehyde, pentanol, and heptanone. Each odor of the three odors was given three successive trials. Each trial lasted for 1 min with an inter-trial interval of 30 s. Effective sniffing was defined as the distance between the mouse's nose and Q-tip being less than 1 cm. The exploring time of the mice in each trial can be obtained via the sum of the sniffing time. 8 mice at 6-month age in each group were used in these experiments. The odor discrimination index of mice was calculated by the difference between the sniffing time of mice in the third trial of every odor and the sniffing time in the first trial of subsequent odor [42].

Screening of potential compounds targeting SNAP25

Considering that the crystal structure of the independent SNAP25 has not been completely dissected, we used the SNAER complex containing the SNAP25 structure for the follow-up study. The PDB file containing SNARE protein information was acquired from the RCSB Protein Data Bank (PDB accession#: 3HD7) [43]. The potential active site residues of SNAP25 in the SNAER complex were predicted by the Site Finder module of the Molecular Operating Environment (MOE) software and the DoGSiteScorer module of protein.plus [44]. The structure of SNAP25 was optimized by adding hydrogen atoms and charges and prepared for subsequent virtual screening. Structure information for 1810 FDA-approved drugs was provided by the ZINC database (https://zinc20.docking.org/). Structural optimizations of small-molecule drugs, including 3D coordinates, hydrogen atoms, and Minimize Energy are done by Chemoffice calculations. First, we performed a shape-based virtual screening using Sailvina software to search for conformations in a grid space of size less than or equal to 22.5 Å, centered on x, y, z = 176.037, 48.139, 152.503. To eliminate false positive results, we restricted the cutoff value (total score ≤ − 6) in the docking score function. Next, we used the Induced Fit function in the DOCK module of MOE to re-optimize the virtual screening results obtained above to ensure functional binding of the drug to the target SNAP25. Functional annotation for residue information of drug-protein interactions is mainly analyzed by PLIP (Protein–Ligand Interaction Profiler, https://plip-tool.biotec.tu-dresden.de/plip-web/plip/index). Visualization and analysis of SNARE pocket prediction and ligand–protein interactions were produced using PyMOL (https://pymol.org/2/), Discovery Studio, and MOE software.

Statistical analysis

SPSS 22.0 and OriginPro 2021b software were used for statistical data analyses. For the data of RT-qPCR, the student’s t test was used for comparison between the two groups. Analysis of variance (ANOVA) was used to explore the significance of the hub genes. Bonferroni corrections were applied to correct for multiple comparisons. The data of the olfactory behavior tests were analyzed using the paired samples t test. For correlation analysis of SNAP25 in HIP and OB/OE samples, the Pearson correlation coefficient was calculated. The statistical graph was drawn with GraphPad Prism8. Quantitative data are shown by the mean ± SEM., and P < 0.05 was considered statistically significant.

Availability of data and materials

The databases that support the findings of this study are publicly available in (GEO) at (https://www.ncbi.nlm.nih.gov/geo/). All data analyzed in this work are included in this article and its supplementary information files.

Abbreviations

AD:

Alzheimer’s disease

PD:

Parkinson’s disease

GEO:

Gene Expression Omnibus

DEGs:

Differentially expressed genes

NDDs:

Neurodegenerative diseases

PPI:

Protein–protein interaction

HIP:

Hippocampal

SN:

Substantial nigra

OB:

Olfactory bulb

OE:

Olfactory epithelium

RT-qPCR:

Real-time quantitative polymerase chain reaction

GO:

Gene Ontology

KEGG:

Kyoto Encyclopedia of Genes and Genomes

MOE:

Molecular Operating Environment

SNARE:

Soluble N-ethylmaleimide-sensitive fusion protein attachment protein receptor

References

  1. Jia L, Quan M, Fu Y, Zhao T, Li Y, Wei C, Tang Y, Qin Q, Wang F, Qiao Y, Shi S, Wang YJ, Du Y, Zhang J, Zhang J, Luo B, Qu Q, Zhou C, Gauthier S, Jia J. Dementia in China: epidemiology, clinical management, and research advances. Lancet Neurol. 2020;19(1):81–92.

    Article  PubMed  Google Scholar 

  2. Baldacci F, Mazzucchi S, Della Vecchia A, Giampietri L, Giannini N, Koronyo-Hamaoui M, Ceravolo R, Siciliano G, Bonuccelli U, Elahi FM, Vergallo A, Lista S, Giorgi FS, Hampel H. The path to biomarker-based diagnostic criteria for the spectrum of neurodegenerative diseases. Expert Rev Mol Diagn. 2020;20(4):421–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Dubois B, Hampel H, Feldman HH, Scheltens P, Aisen P, Andrieu S, Bakardjian H, Benali H, Bertram L, Blennow K, Broich K, Cavedo E, Crutch S, Dartigues JF, Duyckaerts C, Epelbaum S, Frisoni GB, Gauthier S, Genthon R, Gouw AA, Habert MO, Holtzman DM, Kivipelto M, Lista S, Molinuevo JL, O’Bryant SE, Rabinovici GD, Rowe C, Salloway S, Schneider LS, Sperling R, Teichmann M, Carrillo MC, Cummings J, Jack CR Jr. Preclinical Alzheimer’s disease: Definition, natural history, and diagnostic criteria. Alzheimer’s Dementia J Alzheimer’s Assoc. 2016;12(3):292–323.

    Article  Google Scholar 

  4. Lashley T, Schott JM, Weston P, Murray CE, Wellington H, Keshavan A, Foti SC, Foiani M, Toombs J, Rohrer JD, Heslegrave A, Zetterberg H. Molecular biomarkers of Alzheimer’s disease: progress and prospects. Dis Models Mech. 2018;11(5):dmm031781.

    Article  Google Scholar 

  5. Kovacs GG, Milenkovic I, Wöhrer A, Höftberger R, Gelpi E, Haberler C, Hönigschnabl S, Reiner-Concin A, Heinzl H, Jungwirth S, Krampla W, Fischer P, Budka H. Non-Alzheimer neurodegenerative pathologies and their combinations are more frequent than commonly believed in the elderly brain: a community-based autopsy series. Acta Neuropathol. 2013;126(3):365–84.

    Article  CAS  PubMed  Google Scholar 

  6. Beach TG, Monsell SE, Phillips LE, Kukull W. Accuracy of the clinical diagnosis of Alzheimer disease at National Institute on Aging Alzheimer Disease Centers, 2005–2010. J Neuropathol Exp Neurol. 2012;71(4):266–73.

    Article  PubMed  Google Scholar 

  7. Moussaud S, Jones DR, Moussaud-Lamodière EL, Delenclos M, Ross OA, McLean PJ. Alpha-synuclein and tau: teammates in neurodegeneration? Mol Neurodegen. 2014;9:43.

    Article  Google Scholar 

  8. Atri A. The Alzheimer’s disease clinical spectrum: diagnosis and management. Med Clin North Am. 2019;103(2):263–93.

    Article  PubMed  Google Scholar 

  9. Elbaz A, Carcaillon L, Kab S, Moisan F. Epidemiology of Parkinson’s disease. Revue Neurol. 2016;172(1):14–26.

    Article  CAS  Google Scholar 

  10. Braak H, Del Tredici K, Rüb U, de Vos RA, Jansen Steur EN, Braak E. Staging of brain pathology related to sporadic Parkinson’s disease. Neurobiol Aging. 2003;24(2):197–211.

    Article  PubMed  Google Scholar 

  11. Gustavsson T, Syvänen S, O’Callaghan P, Sehlin D. SPECT imaging of distribution and retention of a brain-penetrating bispecific amyloid-β antibody in a mouse model of Alzheimer’s disease. Transl Neurodegen. 2020;9(1):37.

    Article  CAS  Google Scholar 

  12. Hawkes CH, Del Tredici K, Braak H. A timeline for Parkinson’s disease. Parkinsonism Related Disord. 2010;16(2):79–84.

    Article  Google Scholar 

  13. Fearnley JM, Lees AJ. Ageing and Parkinson’s disease: substantia nigra regional selectivity. Brain J Neurol. 1991;114(Pt 5):2283–301.

    Article  Google Scholar 

  14. Kovács T, Cairns NJ, Lantos PL. Olfactory centres in Alzheimer’s disease: olfactory bulb is involved in early Braak’s stages. NeuroReport. 2001;12(2):285–8.

    Article  PubMed  Google Scholar 

  15. Wu N, Rao X, Gao Y, Wang J, Xu F. Amyloid-β deposition and olfactory dysfunction in an Alzheimer’s disease model. J Alzheimer’s Dis JAD. 2013;37(4):699–712.

    Article  CAS  PubMed  Google Scholar 

  16. Dintica CS, Marseglia A, Rizzuto D, Wang R, Seubert J, Arfanakis K, Bennett DA, Xu W. Impaired olfaction is associated with cognitive decline and neurodegeneration in the brain. Neurology. 2019;92(7):e700–9.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Suzuki H, Sugiura S, Nakashima T, Teranishi M, Shimono M, Murotani K, Sakurai T, Uchida Y, Saji N. Cognitive impairment is correlated with olfactory identification deficits in older Japanese adults: a cross-sectional study using objective and subjective olfactory measures. Geriatr Gerontol Int. 2022;22(11):924–9.

    Article  PubMed  Google Scholar 

  18. Iannilli E, Stephan L, Hummel T, Reichmann H, Haehner A. Olfactory impairment in Parkinson’s disease is a consequence of central nervous system decline. J Neurol. 2017;264(6):1236–46.

    Article  PubMed  Google Scholar 

  19. Rey NL, Wesson DW, Brundin P. The olfactory bulb as the entry site for prion-like propagation in neurodegenerative diseases. Neurobiol Dis. 2018;109(Pt B):226–48.

    Article  CAS  PubMed  Google Scholar 

  20. Greicius MD, Srivastava G, Reiss AL, Menon V. Default-mode network activity distinguishes Alzheimer’s disease from healthy aging: evidence from functional MRI. Proc Natl Acad Sci USA. 2004;101(13):4637–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Wang K, Liang M, Wang L, Tian L, Zhang X, Li K, Jiang T. Altered functional connectivity in early Alzheimer’s disease: a resting-state fMRI study. Hum Brain Map. 2007;28(10):967–78.

    Article  Google Scholar 

  22. Song Y, Wu H, Chen S, Ge H, Yan Z, Xue C, Qi W, Yuan Q, Liang X, Lin X, Chen J. Differential abnormality in functional connectivity density in preclinical and early-stage Alzheimer’s disease. Front Aging Neurosci. 2022;14: 879836.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Barresi M, Ciurleo R, Giacoppo S, Foti Cuzzola V, Celi D, Bramanti P, Marino S. Evaluation of olfactory dysfunction in neurodegenerative diseases. J Neurol Sci. 2012;323(1–2):16–24.

    Article  PubMed  Google Scholar 

  24. Fullard ME, Tran B, Xie SX, Toledo JB, Scordia C, Linder C, Purri R, Weintraub D, Duda JE, Chahine LM, Morley JF. Olfactory impairment predicts cognitive decline in early Parkinson’s disease. Parkinsonism Related Disord. 2016;25:45–51.

    Article  Google Scholar 

  25. Cao L, Schrank BR, Rodriguez S, Benz EG, Moulia TW, Rickenbacher GT, Gomez AC, Levites Y, Edwards SR, Golde TE, Hyman BT, Barnea G, Albers MW. Aβ alters the connectivity of olfactory neurons in the absence of amyloid plaques in vivo. Nat Commun. 2012;3:1009.

    Article  PubMed  Google Scholar 

  26. Benseny-Cases N, Álvarez-Marimon E, Aso E, Carmona M, Klementieva O, Appelhans D, Ferrer I, Cladera J. In situ identification and G4-PPI-His-Mal-dendrimer-induced reduction of early-stage amyloid aggregates in Alzheimer’s disease transgenic mice using synchrotron-based infrared imaging. Sci Rep. 2021;11(1):18368.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. He F, Qi G, Zhang Q, Cai H, Li T, Li M, Zhang Q, Chen J, Ming J, Tian B, Zhang P. Quantitative phosphoproteomic analysis in alpha-synuclein transgenic mice reveals the involvement of aberrant p25/Cdk5 signaling in early-stage Parkinson’s disease. Cell Mol Neurobiol. 2020;40(6):897–909.

    Article  CAS  PubMed  Google Scholar 

  28. Takahashi K, Ohsawa I, Shirasawa T, Takahashi M. Early-onset motor impairment and increased accumulation of phosphorylated alpha-synuclein in the motor cortex of normal aging mice are ameliorated by coenzyme Q. Exp Gerontol. 2016;81:65–75.

    Article  CAS  PubMed  Google Scholar 

  29. Toft-Bertelsen TL, Ziomkiewicz I, Houy S, Pinheiro PS, Sørensen JB. Regulation of Ca2+ channels by SNAP-25 via recruitment of syntaxin-1 from plasma membrane clusters. Mol Biol Cell. 2016;27(21):3329–41.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Colom-Cadena M, Spires-Jones T, Zetterberg H, Blennow K, Caggiano A, DeKosky ST, Fillit H, Harrison JE, Schneider LS, Scheltens P, de Haan W, Grundman M, van Dyck CH, Izzo NJ, Catalano SM. The clinical promise of biomarkers of synapse damage or loss in Alzheimer’s disease. Alzheimer’s Res Ther. 2020;12(1):21.

    Article  Google Scholar 

  31. Lista S, Hampel H. Synaptic degeneration and neurogranin in the pathophysiology of Alzheimer’s disease. Expert Rev Neurother. 2017;17(1):47–57.

    Article  CAS  PubMed  Google Scholar 

  32. Guerini FR, Farina E, Costa AS, Baglio F, Saibene FL, Margaritella N, Calabrese E, Zanzottera M, Bolognesi E, Nemni R, Clerici M. ApoE and SNAP-25 polymorphisms predict the outcome of multidimensional stimulation therapy rehabilitation in Alzheimer’s disease. Neurorehabil Neural Repair. 2016;30(9):883–93.

    Article  PubMed  Google Scholar 

  33. Agliardi C, Guerini FR, Zanzottera M, Riboldazzi G, Zangaglia R, Sturchio A, Casali C, Di Lorenzo C, Minafra B, Nemni R, Clerici M. SNAP25 gene polymorphisms protect against Parkinson’s disease and modulate disease severity in patients. Mol Neurobiol. 2019;56(6):4455–63.

    Article  CAS  PubMed  Google Scholar 

  34. Hoerder-Suabedissen A, Korrell KV, Hayashi S, Jeans A, Ramirez DMO, Grant E, Christian HC, Kavalali ET, Wilson MC, Molnár Z. Cell-specific loss of SNAP25 from cortical projection neurons allows normal development but causes subsequent neurodegeneration. Cereb Cortex (New York NY: 1991). 2019;29(5):2148–59.

    Google Scholar 

  35. Halbgebauer S, Steinacker P, Hengge S, Oeckl P, Abu Rumeileh S, Anderl-Straub S, Lombardi J, Von Arnim CAF, Giese A, Ludolph AC, Otto M. CSF levels of SNAP-25 are increased early in Creutzfeldt-Jakob and Alzheimer’s disease. J Neurol Neurosurg Psychiatry. 2022;93:1059–65.

    Article  Google Scholar 

  36. Kivisäkk P, Carlyle BC, Sweeney T, Quinn JP, Ramirez CE, Trombetta BA, Mendes M, Brock M, Rubel C, Czerkowicz J, Graham D, Arnold SE. Increased levels of the synaptic proteins PSD-95, SNAP-25, and neurogranin in the cerebrospinal fluid of patients with Alzheimer’s disease. Alzheimer’s Res Ther. 2022;14(1):58.

    Article  Google Scholar 

  37. Park H, Kim H, Kim S, Cha H. The association between olfactory function and cognitive impairment in older persons with cognitive impairments: a cross-sectional study. Healthcare (Basel, Switzerland). 2021;9(4):399.

    PubMed  Google Scholar 

  38. Öhrfelt A, Brinkmalm A, Dumurgier J, Zetterberg H, Bouaziz-Amar E, Hugon J, Paquet C, Blennow K. A novel ELISA for the measurement of cerebrospinal fluid SNAP-25 in patients with Alzheimer’s disease. Neuroscience. 2019;420:136–44.

    Article  PubMed  Google Scholar 

  39. Brinkmalm A, Brinkmalm G, Honer WG, Frölich L, Hausner L, Minthon L, Hansson O, Wallin A, Zetterberg H, Blennow K, Öhrfelt A. SNAP-25 is a promising novel cerebrospinal fluid biomarker for synapse degeneration in Alzheimer’s disease. Mol Neurodegen. 2014;9:53.

    Article  Google Scholar 

  40. Shen S, Kong J, Qiu Y, Yang X, Wang W, Yan L. Identification of core genes and outcomes in hepatocellular carcinoma by bioinformatics analysis. J Cell Biochem. 2019;120(6):10069–81.

    Article  CAS  PubMed  Google Scholar 

  41. Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative C(T) method. Nat Protoc. 2008;3(6):1101–8.

    Article  CAS  PubMed  Google Scholar 

  42. Rankin CH, Abrams T, Barry RJ, Bhatnagar S, Clayton DF, Colombo J, Coppola G, Geyer MA, Glanzman DL, Marsland S, McSweeney FK, Wilson DA, Wu CF, Thompson RF. Habituation revisited: an updated and revised description of the behavioral characteristics of habituation. Neurobiol Learn Memory. 2009;92(2):135–8.

    Article  Google Scholar 

  43. Khvotchev M, Soloviev M. SNARE Modulators and SNARE Mimetic Peptides. Biomolecules. 2022;12(12):1779.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Volkamer A, Griewel A, Grombacher T, Rarey M. Analyzing the topology of active sites: on the prediction of pockets and subpockets. J Chem Inf Model. 2010;50(11):2041–52.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by the Faculty Startup Fund from Xuzhou Medical University, Xuzhou Science and Technology Innovation Project (KC20090, Jiangsu Province, China), Xuzhou Medical University Open Fund Project (XYKF202114), Key Projects of Jiangsu University Student Innovation and Entrepreneurship Training Program (202110313021Z), and National Demonstration Center for Experimental Basic Medical Science Education (Xuzhou Medical University).

Author information

Authors and Affiliations

Authors

Contributions

YZ and ST designed the study; YZ, ST, and QW wrote the paper; XY and QS analyzed data; QW and CX performed RT-qPCR; LX, HD, FW, and YG bred AD and PD mice and performed the behavioral tests; JL, XX, and KL worked for the genotyping. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Qi Shen, Xiaobo Yu or Yingwei Zheng.

Ethics declarations

Ethics approval and consent to participate

All procedures were approved by the Animal Care and Use Committee of Xuzhou Medical University (Protocol Number: 202007A136).

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

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

Supplementary Information

Additional file 1: Table S1.

Primer list in the RT-qPCR.

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wang, Q., Tao, S., Xing, L. et al. SNAP25 is a potential target for early stage Alzheimer’s disease and Parkinson’s disease. Eur J Med Res 28, 570 (2023). https://doi.org/10.1186/s40001-023-01360-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40001-023-01360-8

Keywords