Synaptosomal Protein of 25 kDa (Snap25) Polymorphisms Associated with Glycemic Parameters in Type 2 Diabetes Patients
A possible role of Snap25 polymorphisms in type 2 diabetes mellitus (T2DM) was evaluated by analyzing three SNPs within intron 1 in a region known to affect the gene expression in vitro. Genomic DNA from 1019 Saudi individuals (489 confirmed T2DM and 530 controls) was genotyped for SNPs rs363039, rs363043, and rs363050 in Snap25 using the TaqMan Genotyping Assay. Significantly higher levels of fasting glucose and HbA1c were detected in T2DM patients carrying the rs363050 (AG/GG) genotypes compared to the (AA) genotype (, , and and , , and , resp.). In these same patients, insulin levels were significantly decreased compared to the (AA) individuals (, , and ). Significant associations were detected between rs363050 (AG/GG) genotypes and increasing fasting glucose levels ( and OR: 1.05), HbA1c levels (OR: 5.06 and ), and lower insulinemia ( and OR: 0.95) in T2DM patients. The minor Snap25 rs363050 (G) allele, which results in a reduced expression of Snap25, is associated with altered glycemic parameters in T2DM possibly because of reduced functionality in the exocytotic machinery leading to suboptimal release of insulin.
Blood glucose levels are regulated by insulin release from pancreatic β-cells via exocytosis, a fine-tuned mechanism mediated by complex machinery. During the exocytotic process secretory granules (SG) of insulin fuse with the plasma membrane; this process is mediated by a number of key players including the soluble N-ethylmaleimide-sensitive factor attachment protein receptor (SNARE) proteins. The SNARE-complex includes the two t-SNARE proteins, synaptosomal protein of 25 kDa (Snap25) and syntaxin 1A, as well as the v-SNARE protein VAMP2. To allow exocytosis the amino terminal of Snap25 binds to syntaxin 1A and the carboxy-terminal binds to VAMP2, forming the four-helical bundle that brings SG in close contact with the plasma membrane, thus enabling fusion to occur. Snap25 is contained in insulin-secreting cells  and is known to modulate several processes apart from the actual fusion event, including the activity of potassium voltage gated (Kv)2.1 channels . Notably, Kv2.1 is the prevalent Kv channel in neuroendocrine and endocrine cells, and the Kv2.1 current repolarizes β-cell action potentials upon exposure to glucose to limit Ca2+ entry and insulin secretion.
Because of the pivotal role of Snap25 in insulin release [2–6] we evaluated a possible role of three Snap25 polymorphisms (rs363039, rs363043, and rs363050) that are localized within intron 1 in a region spanning about 13.8 kb and are known to affect the gene expression in vitro [7, 8], in type 2 diabetes mellitus (T2DM). To this end, fasting glycemia, insulinemia, and the concentration of glycated haemoglobin (HbA1c) were correlated with Snap25 polymorphisms in 489 patients with a diagnosis of type 2 diabetes mellitus (T2DM) and 530 age- and sex-matched healthy controls (HC).
2.1. Patients and Controls
A total of 1019 unrelated adult Saudi individuals (489 confirmed T2DM (280 males and 209 females) and 530 controls (319 males and 211 females)) were randomly selected from the Biomarker Screening Project in Riyadh (Riyadh Cohort), a capital-wide epidemiological study of over 17,000 consenting Saudis coming from different Primary Health Care Centers (PHCCs) in Kingdom of Saudi Arabia (KSA). Written informed consent was obtained from all participants before inclusion in the study. The research was approved by the Ethics Committee of the College of Science, King Saud University, Riyadh, KSA, and the methods were carried out in accordance with the approved guidelines. Diagnosis of T2DM was based on World Health Organization (WHO) proposed cut-off: fasting plasma glucose > or = 7.0 mmol/L or 126 mg/dL. Subjects with comorbidities that needed medical attention or with medical complications (coronary artery disease, nephropathy, thyroid diseases, and end stage renal or liver disease) were excluded from the study. Consenting participants were then requested to visit the nearest participating PHCC for questionnaire administration, anthropometric measurement, and blood extraction. A questionnaire focusing on demographic information and past and present medical history was given to all subjects.
2.2. Anthropometry, Blood Collection, and Biochemical Analyses
Subjects were requested to visit their respective PHCCs following overnight fasting (>10 hours) for anthropometry and blood withdrawal. Anthropometry included height (rounded off to the nearest 0.5 cm), weight (rounded off to the nearest 0.1 kg), waist and hip circumference (centimeters), and mean systolic and diastolic blood pressure (millimeters of Hg) (average of 2 readings). Body mass index (BMI) was calculated as weight in kilograms divided by height in square meters. Fasting glucose and lipid profile were measured using a chemical analyzer (Konelab, Espoo, Finland). Serum insulin was assessed using an enzyme-linked immunosorbent assay (ELISA) with an intra-assay variability of <10% and interassay variation of <15%. All fasting samples fell within the detection range.
2.3. Snap25 Genotyping
Genomic DNA was isolated from whole blood using the blood genomicPrep mini spin kit (GE Healthcare Life Sciences, Piscataway, NJ, USA). SNPs rs363043, rs363039, and rs363050 located in intron 1 of Snap25 gene were genotyped using the TaqMan SNP Genotyping Assays (Applied Biosystems by Life Technologies, Foster City, CA, USA) on an ABI PRISM 7000 Sequence Detection System.
2.4. Statistical Analysis
The Kolmogorov-Smirnov (K-S) test was applied to verify normal distribution of numerical variable scores. All numerical scores resulted in being normally distributed, and they were shown as mean and standard deviation (SD). Hardy-Weinberg equilibrium was tested using Chi-square test. ANOVA calculation was performed to analyse biochemical parameters in relationship with Snap25 genotypes stratified for T2DM and HC. Moreover multivariate stepwise logistic regression analysis was performed, using Snap25 polymorphisms as responsible variables and age, gender, and BMI as covariates.
3. Results and Discussion
Fasting glycemia, insulinemia, and the concentration of glycated haemoglobin (HbA1c) were correlated with Snap25 polymorphisms in 489 patients with a diagnosis of T2DM and 530 age- and sex-matched HC (Table 1), whereas no differences were seen in the distribution of the 3 Snap25 SNPs between T2DM and controls (data not shown). Significantly higher levels of fasting glucose and HbA1c were detected in T2DM patients carrying the rs363050 (AG/GG) (57% of cases) genotypes compared to those carrying the (AA) genotype (43% of cases) (, , and and , , and , resp.). In these same patients insulin levels were significantly decreased as well compared to the (AA) individuals (, , and ) (Table 2). In contrast with these results, no associations were observed between any of the Snap25 genotypes and glycemic parameters in HC (, , and ).
The correlation between the rs363050 (G) minor allele and altered glycemic parameters in T2DM patients was even stronger when multivariate stepwise logistic regression analysis was performed. Thus, when rs363050 (AG/GG) versus (AA) was calculated as the responsible variable and age, gender, and BMI were calculated as covariates, a statistically significant association was detected between (AG/GG) genotypes and increasing fasting glucose levels (, OR: 1.05), with a higher risk in female versus male patients (, OR: 1.63). A similar association was observed as well between the (AG/GG) genotypes and increasing HbA1c levels (, OR: 5.06) and between these genotypes and lower insulinemia (, OR: 0.95) in T2DM patients after adjustment for BMI, gender, and age covariates, even if in this case no specific risk was detected in female patients. Again, no such associations were observed in healthy controls. To rule out a possible role played by confounding factors, diabetes duration, hypoglycemic drug use, insulin resistance condition, and anthropometric parameters were compared between patients carrying the two genotypes, and no significant differences emerged (data not shown).
Insulin secretion depends on the integrity of the SNARE exocytotic machinery, as shown by the observation that decreased expression of t-SNARE associates with an alteration of this process . Snap25, a protein that is expressed in pancreatic endocrine cells, is an essential component of such machinery . Data herein shows that the Snap25 rs363050 (G) allele correlates with increasing glycemia and HbA1c concentration, as well as with reduced insulinemia in T2DM patients alone; this allele was recently demonstrated to be associated with a significantly reduced gene expression in vitro compared to the rs363050 (A) allele . Because no associations between glycemic parameters and Snap25 polymorphism were detected in healthy controls this polymorphism is not likely to be a genetic risk factor for T2DM by itself but rather to correlate with altered T2DM-related metabolic parameters once T2DM has been diagnosed. This effect is possibly mediated by impairment in the SNARE complex, resulting in a reduced functionality of the exocytotic machinery and, as a consequence, in a suboptimal release of insulin.
Snap25 plays an important role in exocytosis, and its polymorphisms have been implicated in modulating the release of neurotransmitters in a number of neurologic conditions including autism and Alzheimer’s disease (AD) . Importantly, the integrity of the insulin signaling plays a pivotal role for cognitive function and memory. Thus, in neurons, insulin stimulates growth, survival, differentiation, metabolism, protein synthesis, synapse formation, and plasticity . The preliminary results herein indicate that the same Snap25 SNPs are present in AD and in T2DM; these data add further evidence to the importance played by an intact insulin signaling pathway in the maintenance of well-structured cognitive functions and reinforce the idea that Snap25 polymorphisms play a role in metabolic disorders. Notably, polymorphisms in the human Snap25 gene have been associated with weight gain after antipsychotic treatment and with altered levels of serum triglycerides [11, 12] and, importantly, with the severity of metabolic syndrome in T2DM  These data, finally, confirm previous results showing a role of Snap25 in insulin secretion [2–6]. Because the incidence of AD is increased in T2DM  it will be interesting to explore whether Snap25 alleles-mediated exocytosis alterations bridge these two pathologies, possibly predicting those T2DM patients with a higher risk of developing AD.
Conflict of Interests
The authors have no conflict of interests to declare. The authors have no financial and nonfinancial conflict of interests.
The authors thank the Deanship of Scientific Research, Prolific Research Group Program (PRG-1436-15), Vice Rectorate for Graduate Studies, and Scientific Research in King Saud University (KSU), Riyadh, Saudi Arabia, for funding the study.
P. E. MacDonald, G. Wang, S. Tsuk et al., “Synaptosome-associated protein of 25 kilodaltons modulates Kv2.1 Voltage-dependent K+ channels in neuroendocrine islet β-cells through an interaction with the channel N terminus,” Molecular Endocrinology, vol. 16, no. 11, pp. 2452–2461, 2002.View at: Publisher Site | Google Scholar
S. Nagamatsu, Y. Nakamichi, C. Yamamura et al., “Decreased expression of t-SNARE, syntaxin 1, and SNAP-25 in pancreatic β-cells is involved in impaired insulin secretion from diabetic GK rat islets: restoration of decreased t-SNARE proteins improves impaired insulin secretion,” Diabetes, vol. 48, no. 12, pp. 2367–2373, 1999.View at: Publisher Site | Google Scholar
M. F. Gosso, E. J. C. De Geus, T. J. C. Polderman, D. I. Boomsma, P. Heutink, and D. Posthuma, “Common variants underlying cognitive ability: further evidence for association between the SNAP-25 gene and cognition using a family-based study in two independent Dutch cohorts,” Genes, Brain and Behavior, vol. 7, no. 3, pp. 355–364, 2008.View at: Publisher Site | Google Scholar
F. R. Guerini, C. Agliardi, M. Sironi et al., “Possible association between SNAP-25 single nucleotide polymorphisms and alterations of categorical fluency and functional MRI parameters in Alzheimer's disease,” Journal of Alzheimer's Disease, vol. 42, no. 3, pp. 1015–1028, 2014.View at: Publisher Site | Google Scholar