Original Research ARTICLE
VSNL1 co-expression networks in aging include calcium signaling, synaptic plasticity, and Alzheimer’s disease pathways
- 1Department of Biostatistics, University of Pittsburgh, Pittsburgh, PA, USA
- 2Department of Psychiatry, University of Pittsburgh, Pittsburgh, PA, USA
- 3Campbell Family Mental Health Research Institute, Centre for Addiction and Mental Health (CAMH), Departments of Psychiatry and Pharmacology & Toxicology, University of Toronto, Toronto, ON, Canada
- 4Department of Neurology, University of Pittsburgh, Pittsburgh, PA, USA
- 5VISN 4 Mental Illness Research, Education and Clinical Center (MIRECC), VA Pittsburgh Healthcare System, Pittsburgh, PA, USA
The visinin-like 1 (VSNL1) gene encodes visinin-like protein 1, a peripheral biomarker for Alzheimer disease (AD). Little is known, however, about normal VSNL1 expression in brain and the biologic networks in which it participates. Frontal cortex gray matter obtained from 209 subjects without neurodegenerative or psychiatric illness, ranging in age from 16 to 91, was processed on Affymetrix GeneChip 1.1 ST and Human SNP Array 6.0. VSNL1 expression was unaffected by age and sex, and not significantly associated with SNPs in cis or trans. VSNL1 was significantly co-expressed with genes in pathways for calcium signaling, AD, long-term potentiation, long-term depression, and trafficking of AMPA receptors. The association with AD was driven, in part, by correlation with amyloid precursor protein (APP) expression. These findings provide an unbiased link between VSNL1 and molecular mechanisms of AD, including pathways implicated in synaptic pathology in AD. Whether APP may drive increased VSNL1 expression, VSNL1 drives increased APP expression, or both are downstream of common pathogenic regulators will need to be evaluated in model systems.
Alzheimer disease (AD) is the most prevalent form of dementia in the United States. It is characterized clinically by declining memory, progressive loss of cognitive ability, and behavioral changes. The incidence of AD increases rapidly with increasing age (1). This suggests a role of brain aging in risk for AD, although the basis of this age-dependence is not established. Neuropathologically, the hallmarks of AD are deposition of extracellular amyloid plaques that are predominantly composed of amyloid-β (Aβ) peptide and intracellular neurofibrillary tangles (NFTs) comprised of hyperphosphorylated microtubule-associated protein tau (2, 3). Other pathologic changes include synapse and neuron loss, and reactive gliosis (4).
Among these pathologies, the strongest correlate of cognitive impairment in individuals with AD is loss of synapses across neocortical regions (5, 6), with excitatory synapses onto dendritic spines particularly affected (7, 8). Substantial evidence now indicates that aggregation of Aβ into soluble oligomers is a primary source of synaptotoxicity in AD (9–14). Although studies continue to elucidate how Aβ acts to eliminate dendritic spines, there is evidence that soluble Aβ inhibits mechanisms of long-term potentiation and/or engages mechanisms of long-term depression, including reducing NMDA receptor dependent Ca2+ influx, mGlutamate receptor (mGluR) activation, and low level caspase-3 activation (10, 12). The final common mechanism for these pathways converge on altered endocytotic recycling of GluR, resulting in reduced synaptic expression of GluR1 and GluR2 containing AMPA receptors and synaptic NMDA receptor (15).
Recently, several studies have identified visinin-like protein 1 (Vilip1), a protein encoded by the visinin-like 1 (VSNL1) gene, as a biomarker of AD. Vilip1 concentrations in cerebrospinal fluid and plasma are elevated in AD subjects relative to normal controls (16, 17) and to non-AD dementia subjects (17). Higher levels of cerebrospinal fluid Vilip1 also predicted a faster rate of cognitive decline (18). Other biomarkers for AD, such as cerebrospinal fluid measures of Aβ, tau, and phospho-tau (19), are strongly implicated in the pathogenesis of AD by genetic, post-mortem, animal model, and in vitro studies. To date, the evidence for Vilip1 is much more limited. We reported that genetic variations in VSNL1 were associated with risk for psychosis in AD (20), a phenotype characterized by more rapid cognitive deterioration than seen in AD subjects without psychosis (21, 22). Qualitative studies have reported that Vilip1 can be detected in association with neuritic plaques and NFTs in neocortex of AD subjects (23), and may contribute to phosphorylation of tau and Ca2+-mediated cell death (24).
Vilip1 is a highly brain expressed member of the visinin-like protein subfamily of neuronal calcium sensors (25). Vilip1, like other subfamily members appears to modify receptor recycling (26). For example, the closely related subfamily member, hippocalcin, is necessary for NMDA receptor dependent long-term depression via GluR endocytosis (27, 28). Whether Vilip1 has effects on synaptic plasticity processes implicated in synapse loss in AD, such as GluR recycling, long-term potentiation, and long-term depression, is not known. However, Vilip1 has a higher affinity for Ca2+ than calmodulin, suggesting it may respond to the lower Ca2+ levels, which induce long-term depression (25).
The above findings are consistent with the hypothesis that VSNL1/Vilip1 may contribute to the risk for AD, possibly via age-dependent alterations in expression or by affecting processes that contribute to synapse or neuronal loss. To date, however, very little is known about normal VSNL1 expression in brain, whether it is modulated by genetic variation, and the brain-related biologic co-expression networks in which VSNL1 participates. To begin to address these questions, we assessed VSNL1 expression in two regions of frontal cortex obtained from 209 subjects spanning the adult age range, and without evidence of psychiatric or neurodegenerative illness. We found that VSNL1 expression was present throughout the adult life span and was unaffected by age, sex, and common genetic variants in cis and trans. VSNL1 co-expression networks included KEGG pathways for calcium signaling, AD, and pathways implicated in synaptic pathology in AD.
Materials and Methods
All of the brain specimens were collected during autopsies conducted at the Allegheny County Office of the Medical Examiner with permission obtained from the subjects’ next-of-kin. The protocol used to obtain consent was approved by the University of Pittsburgh Institutional Review Board (IRB) and Committee for Oversight of Research Involving the Dead. An independent committee of experienced clinicians made consensus DSM-IV diagnoses for each subject, using information obtained from clinical records and structured interviews with surviving relatives. These procedures were IRB approved. Samples from a total of 212 subjects without any DSM-IV diagnosis (i.e., including no diagnosis of a cognitive disorder) were obtained for use in this study.
Upon brain collection, ~2 cm coronal blocks from the right hemisphere were cut through the rostro-caudal extent of the brain and stored at −80°C. The RNA integrity (RIN) of each brain was assessed by chromatography (Agilent Bioanalyzer, Santa Clara, CA, USA). Samples were obtained from two prefrontal cortex (PFC) regions: Brodmann areas (BA) 11 and 47. These areas were selected based on prior findings showing robust age-related changes in gene expression that were highly correlated with other PFC regions (e.g., BA9) (29). Gray matter samples containing all six layers and excluding white matter were harvested from three to four consecutive 20 μm sections and stored in Trizol reagent.
Total RNA was extracted from frozen BA11 and BA47 samples stored in TRIZOL and were processed for microarray analysis using GeneChip Human Gene 1.1 ST from Affymetrix according to manufacturer’s protocol (http://www.affymetrix.com). Gene expression data were extracted using Expression Console build 22.214.171.124. The normalization method is based on quantile normalization to eliminate batch effects. Data from arrays were processed by RMA method. Gene expression probes were processed at gene-level and taken in log2 scale for further analysis. After normalization, 33,297 gene-level probes remained. Three samples were removed from study, one due to poor array quality, one due to outlier effect, and one due to XXY genotype, resulting in a final sample size of 209 subjects (Table 1).
DNA samples for 204 of the subjects were available, and processed on Affymetrix SNP 6.0 arrays, assessing genotype at 906,600 SNPs. Genotype calls were generated using Affymetrix Genotyping Console version 4.1.3. For intensity quality control (QC), we used Contrast QC, which is the per sample QC test in the Affymetrix SNP 6.0 intensity array; two samples were removed after QC, leaving a final sample of 202 subjects for genotype measures. Birdseed v2 algorithm was used for genotyping, using the EM algorithm to drive a maximum likelihood fit of a two dimensional Gaussian mixture model.
Linear regression models were fitted to assess the association between VSNL1 expression level and demographic factors (age, sex, race). These models included as covariates the technical factors brain pH and RIN, as both were significantly associated with VSNL1 expression: p = 7.21e-09 (BA11) and 1.94e-04 (BA47) for pH association and p = 3.66e-12 (BA11) and 1.16e-08 (BA47) for RIN association. There was no significant association of VSNL1 expression with post-mortem interval (PMI, p = 0.26 and 0.36 in BA11 and BA47, respectively) and thus PMI was not included in the models.
eQTL analysis was performed in the 169 Caucasian subjects for whom genotype data were available. PCA analysis did not indicate population substructure within this cohort. All subjects had call rates >98%. SNPs were filtered out using the following criteria: (i) sample missing rate >10%, (ii) Minor allele frequency (MAF) <5%, and (iii) p value of Hardy–Weinberg equilibrium (HWE) test <10−3. SNPs were defined as in cis if they were located within 50 kbp of VSNL1. All other SNPs were defined as in trans.
The eQTL model adjusted for age, pH, and RIN values (each was significantly associated with gene expression in both BA11 and BA47) since the effects of those covariates may confound eQTL findings. The eQTL model with three covariates for a given genotype was:
in which: Gij: gene expression of gene i of subject j; αi: intercept term of gene i; γi: effect of the selected genotype to gene i based on additive model; Xj: genotype of subject i, 0 (homozygous major alleles), 1 (heterozygous calls), and 2 (homozygous minor alleles); βik: the effect of covariates k; k = 1 (age), 2 (pH), and 3 (RIN) in gene i; Sjk: the value of covariate k of subject j; k = 1 (age), 2 (pH), and 3 (RIN);
We used the “Matrix eQTL” R package, a recent computationally efficient package for eQTL analysis (30), to detect the desired trans-eQTLs. We applied the adaptive weighted (AW) Fisher’s method (31) for meta-analysis to combine eQTL p-values from two brain regions (BA11 and BA47). The AW Fisher’s method has the advantage to distinguish study homogeneity and heterogeneity by assigning 0 or 1 study (brain region) weights for each eQTL. Detected eQTLs from the meta-analysis can have three possible resulting weights: (1,1) meaning detected eQTL in both brain regions; (1,0) showing detected eQTL in BA11 but not in BA47 and vice versa for (0,1) weight. In order to avoid heterogeneity of eQTL findings across brain regions, we selected only eQTLs that were identified in both brain regions with (1,1) weights in the AW Fisher’s method.
Co-expression pathway analyses
To reveal the co-expression structure of VSNL1, we selected the top 400 positively and negatively correlated genes in both brain regions. The minimum observed correlation among the selected genes was 0.64, and the associated p-values suggest the corresponding false discovery rate (FDR) is <1e-04. Pathway analysis was performed via over representation analysis based on Fisher’s exact test using pathway information retrieved from KEGG, BIOCARTA, and REACTOME databases. The Benjamini Hochberg procedure was used to control FDR for pathways (32).
There were no significant associations of VSNL1 expression in either BA11 or BA47 with age, race, or gender (Figure 1).
Figure 1. VSNL1 expression as a function of age, race, and sex. The correlation of VSNL1 expression levels with age are shown in the first panels for each area. The second and third panels show the density of expression levels as a function of race and sex. No associations were significant. BA, Brodmann area.
VSNL1 eQTL Analysis
There was no significant association of VSNL1 expression with 40 genotyped SNPs located in the cis-regulatory region (Table S1 in Supplementary Material). Examination of the association of VSNL1 expression with SNPs in trans revealed 27 SNPs with suggestive evidence of association (p < 10−6, Figure 2; Table S2 in Supplementary Material). However, no SNPs reached the threshold for genome wide significance (p < 10−8).
Figure 2. Manhattan plot of SNP associations with VSNL1 expression. Associations are shown for each region, and for the weighted average of the two regions in the meta-analysis. BA, Brodmann area.
VSNL1 Co-Expression Networks
We next evaluated VSNL1 co-expression. Although VSNL1 expression was itself not age dependent, many of its potential co-expression partners show age-dependent alterations in expression, including genes related to neurodegeneration (29, 33). Therefore, we separately evaluated co-expression in subjects under and over 50 years of age. The 400 genes with the greatest positive and negative correlations with VSNL1 expression in each age group are shown in Table S3 in Supplementary Material.
The KEGG, BIOCARTA, and REACTOME pathways showing the most significant loading for genes positively and negatively correlated with VSNL1 expression are shown in Tables 2 and 3, respectively (all pathways are shown in Table S4 in Supplementary Material). In subjects under 50, VSNL1 was positively correlated with genes in the KEGG pathway for AD and the BIOCARTA pathway for P35 signaling in AD (P35Alzheimers). The association with the KEGG pathway was driven by correlations with amyloid precursor protein (APP), ATP2A2, CALM1, CDK5, CHP, GSK3B, ITPR1, NDUFAB1, NDUFV2, PPP3CA, PPP3CB, and UQCRH. The association with the P35Alzheimers pathway was driven by three of the same genes (APP, CDK5, and GSK3B) as well as PPP2CA. In subjects over 50, the co-expression of VSNL1 with these AD pathways was weakened, with CDK5, CHP, and NDUFAB1 no longer among the top co-expressed genes.
Other neuronal calcium sensor family members have been shown to mediate processes dependent on glutamate receptor trafficking, such as long-term potentiation and long-term depression (27, 28, 34). Although the protein product of VSNL1, Vilip1, is known to alter availability of α4,β2-nicotinic acetylcholine receptors at the cell membrane via affects on endocytic trafficking (35), a similar effect on glutamate receptors has not been shown. It is of some interest, therefore, that among the top pathways identified by genes positively correlated with VSNL1 expression are the KEGG pathway for long-term potentiation and REACTOME pathways for trafficking of AMPA receptors. These associations include positive correlations with several AMPA receptor subunits: GRIA1, GRIA2, and GRIA3. Although not among the top pathways identified, there was also a significant positive correlation of VSNL1 expression with the KEGG pathway for long-term depression in both age groups (Table S4 in Supplementary Material).
We evaluated the hypothesis that VSNL1 may contribute to the development of AD by assessing whether VSNL1 demonstrated age-dependent expression and by determining the brain-related biologic co-expression networks in which VSNL1 participates. VSNL1 expression was present throughout the adult life span, but was not correlated with subject age. VSNL1 co-expression networks included AD pathways and pathways implicated in synaptic pathology in AD, including long-term potentiation, long-term depression, and trafficking of AMPA receptors. These latter findings provide an unbiased link in support of the hypothesis that VSNL1/Vilip1 may participate in molecular mechanisms of AD.
We found that VSNL1 is expressed throughout the adult lifespan in human frontal cortex and is independent of subject age. We, and others, have previously reported that there is a substantial overlap between genes demonstrating age-related changes in expression and genes involved in the pathogenesis of neurodegenerative illnesses, including AD (36). The lack of changes in VSNL1 expression with age in our subjects suggests that any contribution of VSNL1 to the development of AD is not via this age-dependent mechanism. However, this does not preclude that VSNL1 expression could vary with age in other brain regions relevant to AD pathogenesis, such as the hippocampus. Also, it remains possible that post-transcriptional processing of VSNL1 into Vilip1 protein may vary with age. Evaluating that alternative will require future study of brain Vilip1 levels in an aging cohort.
We had previously identified an association between genetic variation in VSNL1 and the psychotic phenotype of AD (20). Because one important mechanism by which genetic variants may affect risk for neurodegenerative disease is to alter the transcription of their gene products in brain (37), we evaluated whether VSNL1 expression in frontal cortex is associated within or nearby the VSNL1 locus. We did not detect an association of VSNL1 expression with any of the tested SNPs. However, the possibility that nearby SNPs, other than the ones tested, affect VSNL1 expression cannot be excluded. To assess this likelihood, we used data from the 1000 Genomes project to estimate what proportion of all common variants (minor allele frequency ≥5%) in VSNL1 were correlated with the SNPs tested in our analysis. Only 70% of the common variants were tagged by one or more of our SNPs with an r2 ≥ 0.8. The 30% of poorly tagged common variants includes the most strongly-associated SNP in our GWAS, rs4038131. It is also possible that genetic variants, including those evaluated in the current study, may alter VSNL1 expression only in the presence of neurodegenerative pathology. However, it should be noted that recent tests of SNPs in the AD risk genes ABCA7, BIN1, CD2AP, CD33, CLU, CR1, EPHA1, MS4A6A, MS4A6E, and PICALM for associations with the expression of their respective genes in brain tissue from AD and healthy control subjects have been largely negative (38, 39).
Although the focus of our study was on risk for AD, genetic variants in VSNL1, altered brain VSNL1 mRNA expression, and altered brain levels of Vilip1 protein have been reported in schizophrenia (40–42). Of particular relevance to the current report was the finding that microRNA miR-181b, which is elevated in schizophrenia, can downregulate VSNL1 expression in model systems and correlates with reduced VSNL1 expression in the superior temporal gyrus in subjects with schizophrenia (40). We were able to replicate the negative correlation between VSNL1 expression and expression of pre-miR-181b genes in our normal aging cohort in both BA11 and BA47 (Figure S1 in Supplementary Material). However, how this normative regulation of VSNL1 by miR-181b might be changed in the presence of AD pathology and/or genetic variation in VSNL1, and in particular whether they may interact to underlie the association of VSNL1 genetic variation with psychosis in AD, is not known.
VSNL1 was co-expressed with genes in several distinct pathways. This included co-expression with genes identified with AD pathways and genes involved in mediating synaptic plasticity mechanisms such as long-term potentiation, long-term depression, and trafficking of AMPA receptors. These findings support a role for Vilip1 in these processes, as co-expressed genes tend to be related functionally (43). However, the specific nature of any functional relationship cannot be readily inferred, as co-expression can arise from many sources (44). For example, gray matter homogenates are a mixture of many cell types. Multiple genes that are specifically enriched in a single cell type, e.g. microglia, may be identified as co-expressed (43, 45). Alternatively, genes may share a mechanism regulating their transcription such as a transcription factor in common (46), proximity within the linear sequence of DNA leading to synchronous transcription (47), or colocalization within the spatial configuration of chromosomes (48). Similarly, epigenetic control of transcription via histone acetylation and methylation can lead to co-expression of neighboring genes (49). Finally, mRNA degradation via microRNA binding may also lead multiple targets of a given microRNA to demonstrate correlated expression (50).
In summary, VSNL1 is co-expressed with functional groups and gene transcripts in AD pathways, including APP itself and pathways implicated in synaptic pathology in AD. These findings provide an unbiased link between VSNL1 and molecular mechanisms of AD. However, whether increased expression of APP may drive increased VSNL1 expression, increased VSNL1 expression drives increased APP expression, or both are concurrently altered downstream of another factor will need to be evaluated in model systems. For example, evaluation of VSNL1 expression in transgenic mouse models of AD, or evaluation of APP expression after knockdown of VSNL1 would be indicated. Similarly, whether altering VSNL1 expression modifies synaptic pathology in AD would benefit from testing within in vitro and/or genetic mouse models.
Conflict of Interest Statement
Dr. Robert A. Sweet served as a consultant for Lilly, USA. Mr. Chien-Wei Lin, Mr. Lun-Ching Chang, Dr. George C. Tseng, Ms. Caitlin M. Kirkwood, and Dr. Etienne L. Sibille have no conflicts of interest or financial interests to declare. The Review Editor Konasale Prasad declares that, despite being affiliated to the same institution as authors Chien-Wei Lin, Lun-Ching Chang, George C.Tseng, Caitlin M. Kirkwood, Etienne L. Sibille and Robert A. Sweet, the review process was handled objectively and no conflict of interest exists.
This work was supported by Veterans Health Administration Grant BX000452 and NIH Grants MH071533, AG05133, AG044070, AG027224, MH084060, MH093723-01, and MH094862. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health, the Department of Veterans Affairs, or the United States Government.
The Supplementary Material for this article can be found online at http://www.frontiersin.org/Journal/10.3389/fpsyt.2015.00030/abstract
4. Ingelsson M, Fukumoto H, Newell KL, Growdon JH, Hedley-Whyte ET, Frosch MP, et al. Early Abeta accumulation and progressive synaptic loss, gliosis, and tangle formation in AD brain. Neurology (2004) 62:925–31. doi:10.1212/01.WNL.0000115115.98960.37
5. Terry RD, Masliah E, Salmon DP, Butters N, DeTeresa R, Hill R, et al. Physical basis of cognitive alterations in Alzheimer’s disease: synapse loss is the major correlate of cognitive impairment. Ann Neurol (1991) 30:572–80. doi:10.1002/ana.410300410
7. Baloyannis SJ, Costa V, Mauroudis I, Psaroulis D, Manolides SL, Manolides LS. Dendritic and spinal pathology in the acoustic cortex in Alzheimer’s disease: morphological and morphometric estimation by Golgi technique and electron microscopy. Acta Otolaryngol (2007) 127:351–4. doi:10.1080/00016480601126986
8. Grutzendler J, Helmin K, Tsai J, Gan WB. Various dendritic abnormalities are associated with fibrillar amyloid deposits in Alzheimer’s disease. Ann N Y Acad Sci (2007) 1097:30–9. doi:10.1196/annals.1379.003
13. Lue LF, Kuo YM, Roher AE, Brachova L, Shen Y, Sue L, et al. Soluble amyloid beta peptide concentration as a predictor of synaptic change in Alzheimer’s disease. Am J Pathol (1999) 155:853–62. doi:10.1016/S0002-9440(10)65184-X
14. Näslund J, Haroutunian V, Mohs R, Davis KL, Davies P, Greengard P, et al. Correlation between elevated levels of amyloid beta-peptide in the brain and cognitive decline. J Am Med Assoc (2000) 283:1571–7. doi:10.1001/jama.283.12.1571
15. Hsieh H, Boehm J, Sato C, Iwatsubo T, Tomita T, Sisodia S, et al. AMPAR removal underlies Abeta-induced synaptic depression and dendritic spine loss. Neuron (2006) 52:831–43. doi:10.1016/j.neuron.2006.10.035
16. Lee JM, Blennow K, Andreasen N, Laterza O, Modur V, Olander J, et al. The brain injury biomarker VLP-1 is increased in the cerebrospinal fluid of Alzheimer disease patients. Clin Chem (2008) 54:1617–23. doi:10.1373/clinchem.2008.104497
17. Tarawneh R, D’Angelo G, Macy E, Xiong C, Carter D, Cairns NJ, et al. Visinin-like protein-1: diagnostic and prognostic biomarker in Alzheimer disease. Ann Neurol (2011) 70:274–85. doi:10.1002/ana.22448
19. Jack CR Jr, Knopman DS, Jagust WJ, Shaw LM, Aisen PS, Weiner MW, et al. Hypothetical model of dynamic biomarkers of the Alzheimer’s pathological cascade. Lancet Neurol (2010) 9:119–28. doi:10.1016/S1474-4422(09)70299-6
20. Hollingworth P, Sweet R, Sims R, Harold D, Russo G, Abraham R, et al. Genome-wide association study of Alzheimer’s disease with psychotic symptoms. Mol Psychiatry (2012) 17:1316–27. doi:10.1038/mp.2011.125
21. Ropacki SA, Jeste DV. Epidemiology of and risk factors for psychosis of Alzheimer’s disease: a review of 55 studies published from 1990 to 2003. Am J Psychiatry (2005) 162:2022–30. doi:10.1176/appi.ajp.162.11.2022
23. Braunewell K, Riederer P, Spilker C, Gundelfinger ED, Bogerts B, Bernstein HG. Abnormal localization of two neuronal calcium sensor proteins, visinin-like proteins (vilips)-1 and -3, in neocortical brain areas of Alzheimer disease patients. Dement Geriatr Cogn Disord (2001) 12:110–6. doi:10.1159/000051244
24. Schnurra I, Bernstein HG, Riederer P, Braunewell KH. The neuronal calcium sensor protein VILIP-1 is associated with amyloid plaques and extracellular tangles in Alzheimer’s disease and promotes cell death and tau phosphorylation in vitro: a link between calcium sensors and Alzheimer’s disease? Neurobiol Dis (2001) 8:900–9. doi:10.1006/nbdi.2001.0432
26. Braunewell KH, Klein-Szanto AJ. Visinin-like proteins (VSNLs): interaction partners and emerging functions in signal transduction of a subfamily of neuronal Ca2+ -sensor proteins. Cell Tissue Res (2009) 335:301–16. doi:10.1007/s00441-008-0716-3
29. Erraji-Benchekroun L, Underwood MD, Arango V, Galfalvy H, Pavlidis P, Smyrniotopoulos P, et al. Molecular aging in human prefrontal cortex is selective and continuous throughout adult life. Biol Psychiatry (2005) 57:549–58. doi:10.1016/j.biopsych.2004.10.034
33. Glorioso C, Oh S, Douillard GG, Sibille E. Brain molecular aging, promotion of neurological disease and modulation by Sirtuin5 longevity gene polymorphism. Neurobiol Dis (2011) 41:279–90. doi:10.1016/j.nbd.2010.09.016
34. Saab BJ, Georgiou J, Nath A, Lee FJ, Wang M, Michalon A, et al. NCS-1 in the dentate gyrus promotes exploration, synaptic plasticity, and rapid acquisition of spatial memory. Neuron (2009) 63:643–56. doi:10.1016/j.neuron.2009.08.014
35. Lin L, Jeanclos EM, Treuil M, Braunewell KH, Gundelfinger ED, Anand R. The calcium sensor protein visinin-like protein-1 modulates the surface expression and agonist sensitivity of the alpha 4beta 2 nicotinic acetylcholine receptor. J Biol Chem (2002) 277:41872–8. doi:10.1074/jbc.M206857200
37. Zou F, Chai HS, Younkin CS, Allen M, Crook J, Pankratz VS, et al. Brain expression genome-wide association study (eGWAS) identifies human disease-associated variants. PLoS Genet (2012) 8:e1002707. doi:10.1371/journal.pgen.1002707
38. Karch CM, Jeng AT, Nowotny P, Cady J, Cruchaga C, Goate AM. Expression of novel Alzheimer’s disease risk genes in control and Alzheimer’s disease brains. PLoS One (2012) 7:e50976. doi:10.1371/journal.pone.0050976
39. Allen M, Zou F, Chai HS, Younkin CS, Crook J, Pankratz VS, et al. Novel late-onset Alzheimer disease loci variants associate with brain gene expression. Neurology (2012) 79:221–8. doi:10.1212/WNL.0b013e3182605801
40. Beveridge NJ, Tooney PA, Carroll AP, Gardiner E, Bowden N, Scott RJ, et al. Dysregulation of miRNA 181b in the temporal cortex in schizophrenia. Hum Mol Genet (2008) 17:1156–68. doi:10.1093/hmg/ddn005
41. Martins-de-Souza D, Gattaz WF, Schmitt A, Rewerts C, Marangoni S, Novello JC, et al. Alterations in oligodendrocyte proteins, calcium homeostasis and new potential markers in schizophrenia anterior temporal lobe are revealed by shotgun proteome analysis. J Neural Transm (2009) 116:275–89. doi:10.1007/s00702-008-0156-y
42. Braunewell KH, Dwary AD, Richter F, Trappe K, Zhao C, Giegling I, et al. Association of VSNL1 with schizophrenia, frontal cortical function, and biological significance for its gene product as a modulator of cAMP levels and neuronal morphology. Transl Psychiatry (2011) 1:e22. doi:10.1038/tp.2011.20
44. Gaiteri C, Ding Y, French B, Tseng GC, Sibille E. Beyond modules and hubs: the potential of gene coexpression networks for investigating molecular mechanisms of complex brain disorders. Genes Brain Behav (2014) 13:13–24. doi:10.1111/gbb.12106
45. Hawrylycz MJ, Lein ES, Guillozet-Bongaarts AL, Shen EH, Ng L, Miller JA, et al. An anatomically comprehensive atlas of the adult human brain transcriptome. Nature (2012) 489:391–9. doi:10.1038/nature11405
48. Homouz D, Kudlicki AS. The 3D organization of the yeast genome correlates with co-expression and reflects functional relations between genes. PLoS One (2013) 8:e54699. doi:10.1371/journal.pone.0054699
49. Horvath S, Zhang Y, Langfelder P, Kahn RS, Boks MP, van Eijk K, et al. Aging effects on DNA methylation modules in human brain and blood tissue. Genome Biol (2012) 13:R97. doi:10.1186/gb-2012-13-10-r97
50. Gennarino VA, D’Angelo G, Dharmalingam G, Fernandez S, Russolillo G, Sanges R, et al. Identification of microRNA-regulated gene networks by expression analysis of target genes. Genome Res (2012) 22:1163–72. doi:10.1101/gr.130435.111
Keywords: visinin-like 1, visinin-like protein 1, Alzheimer disease, co-expression networks, calcium signaling, synaptic plasticity
Citation: Lin C-W, Chang L-C, Tseng GC, Kirkwood CM, Sibille EL and Sweet RA (2015) VSNL1 co-expression networks in aging include calcium signaling, synaptic plasticity, and Alzheimer’s disease pathways. Front. Psychiatry 6:30. doi: 10.3389/fpsyt.2015.00030
Received: 14 November 2014; Accepted: 12 February 2015;
Published online: 09 March 2015.
Edited by:Margit Burmeister, University of Michigan, USA
Reviewed by:Ju Wang, University of Virginia, USA
Konasale Prasad, University of Pittsburgh School of Medicine, USA
Thomas G. Schulze, University Medical Center Göttingen, Germany
Copyright: © 2015 Lin, Chang, Tseng, Kirkwood, Sibille and Sweet. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Robert A. Sweet, Departments of Psychiatry and Neurology, University of Pittsburgh, Biomedical Science Tower, Room W-1645, 3811 O’Hara Street, Pittsburgh, PA 15213-2593, USA e-mail: firstname.lastname@example.org
†Chien-Wei Lin and Lun-Ching Chang have contributed equally to this work.