Skip to main content

Age-associated bidirectional modulation of gene expression in single identified R15 neuron of Aplysia



Despite the advances in our understanding of aging-associated behavioral decline, relatively little is known about how aging affects neural circuits that regulate specific behaviors, particularly the expression of genes in specific neural circuits during aging. We have addressed this by exploring a peptidergic neuron R15, an identified neuron of the marine snail Aplysia californica. R15 is implicated in reproduction and osmoregulation and responds to neurotransmitters such as acetylcholine, serotonin and glutamate and is characterized by its action potential bursts.


We examined changes in gene expression in R15 neurons during aging by microarray analyses of RNAs from two different age groups, mature and old animals. Specifically we find that 1083 ESTs are differentially regulated in mature and old R15 neurons. Bioinformatics analyses of these genes have identified specific biological pathways that are up or downregulated in mature and old neurons. Comparison with human signaling networks using pathway analyses have identified three major networks [(1) cell signaling, cell morphology, and skeletal muscular system development (2) cell death and survival, cellular function maintenance and embryonic development and (3) neurological diseases, developmental and hereditary disorders] altered in old R15 neurons. Furthermore, qPCR analysis of single R15 neurons to quantify expression levels of candidate regulators involved in transcription (CREB1) and translation (S6K) showed that aging is associated with a decrease in expression of these regulators, and similar analysis in three other neurons (L7, L11 and R2) showed that gene expression change during aging could be bidirectional.


We find that aging is associated with bidirectional changes in gene expression. Detailed bioinformatics analyses and human homolog searches have identified specific biological processes and human-relevant signaling pathways in R15 that are affected during aging. Evaluation of gene expression changes in different neurons suggests specific transcriptomic signature of single neurons during aging.


Aging is a ubiquitous process that involves the progressive deterioration of physiological functions. Many of the behavioral changes associated with aging in humans are also observed in lower organisms. These changes include, but are not limited to, a decline in learning and memory, locomotor ability, olfactory sensitivity, and circadian rhythmicity [1, 2]. Since behavior is largely affected by the functional state of the nervous system, age-related behavioral declines could very likely indicate the perturbations in the brain. As an organism ages, most tissues undergo multilayered changes and are subject to cellular damage that accumulate with age. It has been well established that there are changes in gene expression in both model organisms and humans as they age. DNA microarray studies on whole worms and flies have also been used to profile transcriptional changes of aging [35]. Large-scale transcriptional changes in mice and primates identified a number of functional genes that change expression levels during aging [6, 7].

Understanding the gene expression changes that occur during aging in the human brain is of particular interest due to its relation to both normal and pathological neurodegeneration and is particularly challenging due to the inherent complexity of the brain at the anatomical and cellular level. Like other biological systems, the complexities of nervous systems at different levels have been studied using a variety of animal models [813]. These animal models have been instrumental with the identification of key neuronal functions in humans and have demonstrated conserved signaling pathways that are critical in aging and neurodegenerative diseases [1418].

Fundamental brain functions are performed by subsets of neurons – neuronal circuits – that are classically defined as distinct networks of hundreds of neurons. The development of neural circuits is initiated with synaptic connections and refined by systematic neuronal activity. Due to the close relationship between neural connectivity and neural activity throughout the brain, it is essential to consider how neural circuitry function at the single neuronal level. Recent studies suggest that experience can dictate the number of neurons or circuits regulating a specific behavior [19, 20]. These studies emphasize the significance of probing gene functions at the circuit and single neuron level.

To understand aging at the single neuron level, we explored the marine snail Aplysia californica that are characterized by a simple central nervous system (CNS). Aplysia CNS has been extensively studied to understand molecular and physiological changes during memory storage and neural circuitries underlying specific behaviors. For example, Kandel and colleagues identified neural circuits underlying the gill withdrawal reflex response and have identified several key determinants of learning and memory storage [21, 22]. Studies by Peretz and coworkers have described the age dependent sensitivity of the gill withdrawal reflex and osmoregulation [2326]. A key advantage of using Aplysia is that specific behavior is mediated by a simple circuitry [21, 22, 2729]. Moreover the changes associated with a specific behavior can be mapped back to the level of single neurons, thus providing access to molecular and biochemical mechanisms of basic neuronal functions at single-neuron resolution. This possibility offer one of the biggest advantages of employing Aplysia for studying molecular and cellular basis of behavior when compared to other model organisms such as worms, flies and mice. Exploring these advantages, age-associated changes in gene expression of identified neurons R2 and LPl5 from abdominal ganglia of young (3–4 months) and old (7–10 months) animals were recently described [30]. Significant changes in gene expression in old R2 neurons, and a unique set of genes differentially expressed during aging were identified [31].

Senescence in Aplysia is characterized by declining weight, reproduction (egg laying), reflex response and finally death in 11–13 months [32]. Impairment of the long-term retention of habituation, prevention of acquisition of sensitization in the gill withdrawal reflex (GWR) and reduction of the response to food stimuli [3335] are associated with aging in Aplysia. At the macromolecular level aging associated changes include transcriptional and translational changes [31] in response to circadian rhythm [31], and the change in electrophysiological properties in identified neurons [34].

Despite these studies the details of how aging alters gene expression and specific signaling pathways in single neurons and circuitries remains to be determined. To address this, we have focused on age-associated gene expression changes in Aplysia neuron R15. The R15 neuron is one of the best-studied neurons used to understand single neuron burst dynamics [3639], neuronal circuitry functions and cellular basis of behavior [40]. R15 has been a subject of extensive electrophysiological studies due to its role in egg laying [4143], neurochemical modulation [44, 45].

Using custom designed microarrays [36] we compared the gene-expression profiles of single R15 neurons from mature (6 months; reproductively mature) and old (11–12 months; close to senescence) animals and identified several genes that are bidirectionally regulated. Our selection of the specific ages (6 month group and 11–12 months group) for this study is based on the fact that previous studies in Aplysia demonstrated robust changes in behavioral readouts and electrophysiological properties occurring between 4–6 months when the animal is attaining sexual maturity [2, 24, 35], and the behavioral and physiological expressions start to get weakened in animals 250 days (8.5 months) of age and older [37]. Therefore we assumed that by selecting ages 6 month and 11–12 months, we would be able to capture the significant differences in gene expression and signaling pathways. Next, we mapped these differentially expressed genes to human functional pathways and identified three networks involved in cell death, neurological diseases and cell morphology. Finally, using single neuron qPCR analyses, we found that aging alters expression of key transcription regulators (CREB1, CREB2) and translation modulator (S6 kinase) in identified neurons R15, L7, L11 and R2 in unique ways.


Identification of differentially expressed genes during aging of single identified neuron R15

Earlier studies have shown that Aplysia exhibit behavioral changes associated with aging [33]. We assumed that age-dependent changes in specific behaviors in Aplysia are associated with specific molecular changes in single neurons. To investigate age associated molecular changes in single neurons, we studied R15 neurons isolated from two age groups that are six months apart, sexually mature (6 months old from hatching) and old (11–12 months from hatching). Aplysia reach senescence in laboratory conditions by 12–14 months [38, 39]. As described in the “Background” section, R15 has been subjected to extensive electrophysiological characterization to gain insights into the biophysical mechanisms underlying neuron endogenous bursting [4547], and is implicated in reproduction, and osmoregulation [4143, 48].

The overall strategy employed to identify differentially expressed genes is shown in Figure 1. To identify age-associated changes, we followed an unbiased gene expression analysis using custom designed microarray that contain a 44,000-oligonucleotide array (60-mer oligonucleotide sequences; Agilent Technologies) that were designed from each non-redundant sequence in the Aplysia EST database [36, 49].

Figure 1
figure 1

Preparation of R15 single neurons from Aplysia abdominal ganglia for genomic analysis. Aplysia central nervous system is organized into nine different ganglia. Identified neurons of abdominal ganglia are shown. Single R15 neuron was isolated from anesthetized Aplysia and immediately stored in Trizol at -80°C. Total RNA was isolated using standard Trizol RNA isolation method and amplified using MessageAmp™ II aRNA Amplification Kit. Quality control and quantification of freshly isolated CNS was performed using a Bioanalyzer and by spectrophotometry.

Single R15 neurons were isolated from abdominal ganglia of mature and old animals, each group in biological triplicates. We carried out dual color microarray analyses following linear amplification of total RNAs from individual neurons. Scatter plots and correlation coefficients for each single neuron were calculated (Additional file 1: Figure S2A). The correlation between the biological replicates in each group was calculated using principle component analysis (PCA) mapping as shown in Figure 2A. Figure 2A shows the first three principal components of microarray analysis data (PC1, PC2 and PC3) in X, Y and Z respectively, and demonstrate correlation between the expression profile of the three single neurons per group (Mature Red # 1–3, Old blue# 1–3). Similarly the PCA was calculated for the technical replicates for each neuron and the first three principal components (PC1, PC2 and PC3) in X, Y and Z respectively showed a positive correlation (Additional file 1: Figure S2B). Our analysis of microarray data have identified 1083 probe sets that are differentially expressed based on the criteria of a fold change of ≥2 and adjusted p < 0.05 between mature versus old neurons (Figure 2B).

Figure 2
figure 2

Heat map analysis of microarray data showing differentially expressed genes between mature and old single R15 neurons. (A) Principle Component Analysis (PCA) mapped scatter plot: The global gene expression profiles of the R15 neurons from both mature and old Aplysia were analyzed by PCA. The figure represents the first two principal components of microarray analysis data (PC1, PC2 and PC3) in X, Y and Z respectively, and demonstrated the expression profile of the three single neurons per group (Mature Red # 1–3, Old blue# 1–3). (B) Heat map analysis of microarray data showing hierarchical clustering of 1083 differentially expressed probes between R15 mature and old single neurons. Each group has three single neurons tested. Red or green colors indicate differentially up or downregulated genes, respectively. The mean signals were background corrected and transformed to the log2 scale. Genes with at least 2-fold changes with p < 0.05 at the 95% confidence level were considered as significant.

The 1083 differentially expressed probes (Additional file 2: Table S1) were annotated using a new Aplysia sequence resource database generated by the Institute for Genome Sciences (IGS), at the University of Maryland. The sequences were generated by The Broad Institute and were distributed to the community ( through the Institute of Genome Sciences (IGS), University of Maryland. After considering a 2-fold change cut-off, we identified 572 differentially expressed transcripts, 247 were enriched in R15 mature and 325 were enriched in R15 old neurons. This corresponded to ~56% upregulated transcripts and ~43% downregulated transcripts during aging. The top 9 transcripts that are enriched in mature or old are listed in Table 1. The complete list of differentially expressed transcripts in mature and old are listed in Additional file 2: Table S2. The differentially expressed genes include, kinases (MAP kinase), phosphatases (phosphoglycolate phosphatase), ion channels (C9 serotonin receptor), neuropeptides (L11), transcription factors (CREB2) protein synthesis regulators (S6K), and those involved in regulating metabolism (adenylate cyclase).

Table 1 Age associated differentially expressed genes in mature and old R15 single neurons

GO analysis of differentially expressed genes in R15

To understand the biological significance of the microarray data the annotated Aplysia sequences were mapped to gene ontology (GO) biological process category. We identified several GO categories specific among genes upregulated or downregulated with age (Figure 3A and B). We found 178 genes upregulated in the mature neurons that were mapped into 22 biological process categories at level three of the GO analysis. However, there are only 18 GO categories for 181 genes that are downregulated in the mature neurons (Additional file 2: Table S3).

Figure 3
figure 3

Functional annotation of differentially expressed genes based on Gene Ontology (GO) classification. Pie chart representing numbers of upregulated (A) from a total of 178 gene IDs and downregulated (B) from a total of 181 gene IDs in corresponding Gene Ontology Biological Process layer 2 categories in the mature R15 single neuron. Biological processes that are uniquely regulated are labeled using a red rectangle. Numbers show the unique hits corresponding to a biological process category.

Four GO categories (cell-cell signaling (GO:0008219), macromolecule localization (GO:0043170), response to endogenous stimulus (GO:0051234), and response to external stimulus (GO:0051716) are unique to the upregulated gene data set of mature neurons (Figure 3A). With the exception of a few processes, the number of sequences upregulated or downregulated in each GO category is roughly similar, indicating a bidirectional differential expression during aging.

Further analysis of the key functions of specific genes belonging to a particular GO category and the directionality of the expression suggest their possible role in altered biological functions associated with aging. For example, in the category “catabolic process” (GO:0007049) enriched sequences are significantly higher than the downregulated sequences. One of the sequences, kinesin 1b, is downregulated with aging. Kinesin is a motor protein involved in the transport of cargos such as synaptic vesicles and mitochondria from cellbody to synapses and is a critical mediator of long-term memory storage. On the other hand, eukaryotic translation initiation factor 2-alpha kinase 3 that is involved in the repression of global protein synthesis, is upregulated with age as its expression is higher in the neurons of old Aplysia. Key regulators of signaling cascades such as MAP kinase 14 and MAP kinase-activated protein kinase 2, and transcription factor 4 were consistently downregulated in R15 with age. Intriguingly, adenylate cyclase and Churchill were upregulated in R15 with age. Adenylate cyclase is involved in cyclic AMP metabolism whereas Churchill is involved in neural development.

Identification of human homologs of Aplysia genes that are altered during aging

Next, we searched for human homologs of genes that are bidirectionally regulated in R15 during aging. Of the 572 gene IDs from the differentially expressed gene set between mature and old R15 neurons, 190 IDs were successfully mapped into different human pathways using Ingenuity Pathway Analysis (IPA). The top 12 functional categories that met stringent statistical criteria and significant representation of genes are illustrated in Figure 4A. The mapped datasets contain either upregulated or downregulated genes associated with specific functional pathways (Additional file 2: Table S4).

Figure 4
figure 4

Mapping of differentially expressed genes into human pathways. Differentially expressed transcripts showing a +/- 2-fold change or greater were selected from the microarray data were analyzed using IPA software to identify potential homologous human pathways. From the 572 gene IDs of the differentially expressed gene set between R15 mature and old single neurons, 190 human homologs were identified. (A) Top 12 pathways from IPA analysis (B) Network 1, identified from IPA analysis, has 25 genes differentially expressed that are associated with cell signaling, cell morphology, and skeletal muscular system development. Node (gene) and edge (gene relationship) symbols are shown along with network diagram. Red shows genes that are upregulated and green indicates down regulation of the genes. Uncolored notes represent genes that were not identified in our data set as differentially expressed in our analysis, but relevant to this network.

We next focused on the top three signaling networks from the IPA analysis. Network 1 (Figure 4B) has 25 genes differentially expressed and associated with cell signaling, cell morphology, and skeletal muscular system development. FOXO1, (a transcription factor of the Fox family) one of the upregulated genes in the network, has been implicated in aging process [50, 51]. Network 2 has 25 genes differentially expressed that are involved in cell death and survival, cellular function maintenance and embryonic development (Additional file 1: Figure S3A). Two genes in Network 2 ATR (Ataxia Telangiectasia and Rad3 related) and MLH1 (MutL Homolog 1) are upregulated in the mature neuron and earlier studies showed their role in aging and longevity [52, 53]. Network 3 has 15 genes differentially expressed and is associated with neurological diseases, developmental and hereditary disorders (Additional file 1: Figure S3B). Thus our mapping of Aplysia single neuron gene expression signatures into the human database has identified several evolutionarily conserved genes that are important for human aging.

Aging of R15 is associated with specific changes in expression of CREB and S6 kinase genes

Transcription and translation are two molecular processes that undergo dynamic changes during aging process [5458]. To quantitate changes in transcription and translation in the R15 neuron during aging, we carried out quantitative PCR (qPCR) analysis of expression of CREB1 and CREB 2 (transcription factors) and S6 kinase (S6K, a regulator of translation) in mature and old abdominal ganglia (Additional file 1: Figure S1) and individual neurons. Importantly, these genes are identified as downregulated genes in old R15 neuron from our microarray analysis. We find from qPCR analysis that CREB1 in mature and old ganglia did not show any significant differences in expression (Additional file 1: Figure S1) whereas the single neuron analysis showed significant changes in expression during aging. Specifically, CREB1 showed significant change in the expression levels with 10 fold decrease in the old neuron, where as CREB 2 is ~ 3 fold and S6K is ~1.3 fold downregulated in the old neuron (Figure 5B) confirming our microarray data. Student t test was used to evaluate the statistical significance while calculating the fold change between mature and old single neurons and a p-value < 0.05 is considered statistically significant.

Figure 5
figure 5

Aging related changes in the expression of CREB1, CREB2 and S6K in identified single neurons of the abdominal ganglia. Single neuron qPCR analysis of expression of CREB1, CREB2 and S6K mRNAs in old neurons when compared to that of mature neurons are shown in bar graphs. (A) Identified neurons in the abdominal ganglia of Aplysia. Single neurons analyzed in this study are color labeled; (B) In the case of old R15, all three genes are significantly downregulated; (C) In old R2, CREB1 and S6K are significantly upregulated. CREB2 did not change with aging in R2; (D) In old L11, CREB1 and CREB2 do not show any significant down regulation, whereas S6K is downregulated significantly; (E) In old L7, all the three genes are downregulated significantly. Figures 5F-H represent the percent change in gene expression of individual genes in different identified neurons of abdominal ganglia during aging. (F) Significant changes in expression of CREB1 in R15 (90 ± 5%), R2 (80 ± 5%) and L7 (50 ± 5%). The change in L11 was not significant; (G) CREB2 showed significant changes in expression levels in R15 (70 ± 5%) and L7 (70 ± 5%) and but no significant changes in R2 and L11; (H) S6K expression changes significantly in all four neurons, R15 (30 ± 5%), L11 (50 ± 5%), L11 (70 ± 5%), and R2 (400 ± 15%). Ct values of target genes were normalized to 18S rRNA internal standard. Each measurement consists of 5 biological repeats, 4 technical repeats and the error bars indicate standard error of the mean (SEM). Student’s t test determined statistical significance where **P-value < 0.01, ***P-value < 0.001.

Next we examined whether CREB and S6K are downregulated during aging of other neurons in the abdominal ganglia. To address this question we studied gene expression in three identified neurons, R2, L7 and L11 (Figure 5A). L7 and L11 are two motor neurons and R2 is a cholinergic neuron. L7 is well characterized for its role in learning and memory of gill withdrawal reflex [16, 59, 60]. Using qPCR we quantified the CREB1, CREB2 and S6K transcript levels in mature and old neurons. In contrast to R15, mRNAs for CREB1 (~2 fold) and S6K (~5 fold) were upregulated significantly whereas CREB2 mRNA was not changed (Figure 5C) in R2. Interestingly, we found that in old L11 neuron, CREB1 and CREB2 mRNAs were not changing significantly, however S6K mRNA was ~1.5 fold downregulated (Figure 5D) whereas in L7 neurons, CREB1, CREB2 and S6K mRNAs were downregulated significantly (~2 fold, ~3 fold and ~2 fold respectively; Student’s t test; p < 0.05) (Figure 5E).

Aging induced changes in gene expression are neuron specific

To further understand the neuron-specific changes during aging, from single neuron qPCR measurements we calculated the percentage changes in the expression of each transcript in different neurons. We observed 80 ± 10% changes in the expression of CREB1 mRNA in R15 (downregulated) and R2 (upregulated) neurons whereas in L7, CREB1 expression was 50 ± 5% downregulated. Unlike significant changes in CREB1 mRNA levels in R15, R2 and L7 the expression of CREB1 mRNA did not significantly change in L11 single neuron (Figure 5F) during aging. Similarly, we measured the percentage change in CREB2 mRNA levels and found that in R2 and L11, the CREB2 mRNA levels were not affected by aging, however there were significant changes in R15 (60 ± 10%) and L7 (65 ± 10%) (Figure 5G). Interestingly, S6K mRNA expression was affected by aging in all neurons. However the magnitude of change varied significantly among these four neurons. In the R15 and L11 neurons S6K mRNA expression changed 35 ± 10% (downregulated in both) whereas in L7 neuron the percentage change was 60 ± 10% (downregulated). Change in R2 neuron was very different from the rest of the neurons with 400 ± 15% upregulation (Figure 5H).


Aging is a fundamental biological process that has distinct phenotypical, physiological and molecular level characteristics. The gene expression level change hypothesis of aging is one of the principal theories of age-associated changes in cellular functions in many tissues, including the brain [30, 56, 61]. Accordingly, the functional genomics approach is one of the useful ways for describing some of the molecular characteristics of neuronal changes underlying normal aging and neurodegeneration [30, 6265]. Considering the cellular heterogeneity of the brain and that behavior is regulated by specific neural circuitries, single circuit or single neuron expression analysis could be a better way to look into the specific molecular changes associated with behavioral decline during aging. Such studies might provide novel therapeutic targets and accelerate future drug discovery efforts. Our comparative analysis of CREB1 expression during aging in the intact ganglia and individual neurons highlight the importance of single neuron study of aging. While we did not find any statistically significant differences in the expression of CREB1 in the abdominal ganglia (Figure 1) during aging, we observed significant differences in CREB1 expression in our single neuron measurements. In this study, we have explored the age-associated molecular changes in four different identified neurons (R15, R2, L7 and L11) in the abdominal ganglia of Aplysia.

Bidirectional changes in gene expression during aging

Our microarray analyses of gene expression in single R15 neurons have identified a subset of genes, 572 of which are bi-directionally regulated during aging. We studied two age groups, sexually mature (6 months) and old animals that are close to senescence (11–12 months). Of the 572 genes we found that 247 genes were upregulated and 325 genes were downregulated in mature R15 neurons when compared to old R15 neurons. Based on the conservative estimate that Aplysia contain approximately 20,000 genes and that our custom oligonucleotide array correspond to approximately 11,000 genes, we estimate that this array contain about 50-60% expected transcriptome. Therefore the change in the gene expression of 1083 transcripts suggests that approximately 10% of the genes are bi-directionally regulated (≥ 2 fold change, p-value ≤ 0.05) in the R15 neuron during aging. However this change might not reflect the net changes in transcriptome during aging. The changes only reflect mRNA regulation in the two age groups we studied. It is quite possible that noncoding transcriptome might be regulated differently.

In our study, differentially expressed genes from both age groups fall into several functional categories that include transcription, translation, protein synthesis, degradation and metabolism. This alteration of gene expression identified in the R15 neuron is consistent with the molecular aging data from other organisms ranging from fly to human where a considerable number of genes that belong to different functional categories undergo dynamic gene expression during aging [66]. Our finding that specific sets of genes that are upregulated or downregulated in single old neurons, indicate that age-associated changes are bidirectional at the molecular level. The age-linked decline in many physiological processes are in fact associated with an increase in the expression of many genes that could alter specific signaling pathways by generating an imbalance in the functional networks, hence affecting key biological functions. Consistent with this idea, the GO category analysis of microarray data in Figure 3 shows that most of the differentially expressed genes share common GO. GO terms corresponding to cell-cell communication, response to internal stimulus, response to external stimulus, and macromolecule localization was absent in old neurons. These biological processes are vital for the maintenance of neuronal health and we find that during aging these processes become defective. Thus our studies suggest significant signaling imbalance associated with aging.

Regulation of human relevant signaling pathways in R15 during aging

We next analyzed our data in the context of human aging by mapping on to human pathways using Ingenuity Pathway Analysis (IPA). Differentially expressed genes were further mapped to identify signaling pathways homologous to human. 190 human genes were mapped in the IPA gene list marking ~35% of annotated genes differentially expressed in Aplysia. The top three pathways are associated with cell death and survival, cellular function and maintenance, cell signaling and neurological disease. Further supporting our notion that aging is associated with an imbalance in signaling, we find that both upregulated and downregulated genes were mapped to the same network. Thus the up or down regulation of nodes of the networks will affect the output of the network, affecting various biological processes. At the single cell level, few key biological processes that are affected by aging include intra and inter-cellular transport, cell-cell communications, transcriptional and translational regulation and epigenetic regulations, all of which significantly change the functioning profile of a cell. Therefore these results are likely to be relevant in the context of human aging.

Transcription and translation are differentially regulated in single neurons during aging

Transcription and translation are fundamental to cell maintenance and survival. CREB has been studied extensively in the context of brain specific functions such as synaptic plasticity [6769], long-term memory [70] and aging [71, 72]. Recently CREB has been shown to be involved in human aging [73], cognitive and neurodegenerative disorders [74]. In mice, CREB regulate aging related genes such as the A-T-mutated gene (ATM) [72]. Consistent with these observations, we find that the expression of CREB is altered at the single neuron level during aging.

Ribosomal S6 kinase (S6K), a translation regulator, is differentially expressed in R15 during aging. S6K is a downstream target of rapamycin (mTOR) pathway, a key modulator of aging [75]. S6K is involved in several signaling pathways [76, 77] that contribute to various pathological states, including aging-related pathology [78, 79] and extending the life span of mammalian cells [79]. In agreement with these results obtained for mammalian aging, we find from our single neuron microarray and qPCR analysis that S6K is downregulated during aging of Aplysia R15 neuron.

Differential aging of neural circuitries

Moroz and Kohn (2010) have shown, using gene expression analysis, that two cholinergic neurons, LPI1 and R2, age differentially [30]. We studied CREB1, CREB2 and S6K mRNA expression in four neurons (R15, L7, L11 and R2). In support of the differential aging idea, we find that age-dependent changes in expression of these genes are neuron specific. We made two conclusions from the gene expression analysis of these four identified neurons. First, we find that during the aging process expression of genes are differentially affected in neurons. In R15 and L7 neurons, the CREB1, CREB2 and S6K showed significant decrease in expression, whereas in L11 only expression of S6K was affected. In R2, unlike the other three neurons during aging, CREB1 and S6K were upregulated. Second, the extent to which expression of each gene is altered depends on the neuron. For example, CREB1 showed ~85% decrease in expression in R15 whereas L7 showed a 50% decrease. On the other hand, the percentage decrease in CREB2 in R15 and L7 were not significantly different. Interestingly, the CREB2 levels did not change in old L11 and R2. Since these neurons are part of circuitries involved in different behaviors, our results suggest that as in the case of single neurons, different neural circuitries could also age differently and may directly contribute to the specific neuronal correlates of aging. Consistent to this idea, it was shown that different regions of mammalian brain age differently and that changes in specific circuitries are associated with behavioral decline during normal aging [8082].


Our study demonstrates the importance of analysis of single neurons to understand transcriptome changes associated with aging. We conclude from our single neuron analysis that genes undergo changes in expression levels during aging. These changes could be bidirectional- sets of genes upregulated and sets of genes downregulated during aging. Importantly, several molecular changes during aging of mammalian cells are identified in aging Aplysia neurons suggesting that molecular changes associated with aging involve several conserved mechanisms. Furthermore, all neurons do not necessarily undergo similar changes in gene expression. For example, in the case of the R2 neuron, transcriptional and translational regulators that we studied were upregulated during aging whereas expressions of these regulators were downregulated in the other three neurons. Also the magnitude of change in expressions of genes is specific to individual neurons. This now poses a question - how and why aging affects neurons differently? Whole transcriptome and epigenome analyses of identified neurons might address this challenge.


The Institutional Biosafety Committee of The Scripps Research Institute (TSRI) approved all the experimental protocols (IBC Protocol 2010-019R1) described in this study. There are no ethical approvals required for the research using invertebrate animals, such as Aplysia, which are not included in the list of regulated animals. However, we have discussed the experiments with the Institutional Animal Care and Use Committee of TSRI and every effort was made to minimize suffering of Aplysia.

Single neuron isolation and RNA analysis

Mature (6 months) and old (11–12 months) Aplysia were obtained from the NIH/University of Miami National Resource for Aplysia and maintained at 16°C under 12:12 light–dark conditions in the Instant Ocean artificial sea water (ASW). Methodology to isolate single neurons for RNA analysis was recently described by our laboratory (Akhmedov et al., 2013, Journal of Visualized Experiments, Accepted). Briefly, prior to dissection, animals were anesthetized by injection of isotonic MgCl2 (337 mM) at a volume of 50%–60% of their body weight. After a 5–10 minute incubation, abdominal ganglion was removed by surgical operation and subjected to 0.1% protease treatment for 30 minutes at 34.5°C followed by ASW perfusion (flow rate: 150 μl/min) at room temperature. Ganglion was desheathed under a binocular stereomicroscope followed by a wash with 100% ethanol.

Selected identified neurons (R15, L7, L11 and R2) were recognized based on their characteristic position in abdominal ganglia and isolated using glass microelectrodes, and placed in Trizol. RNA was isolated according to the standard Trizol (Invitrogen) protocol and RNA concentration was determined spectrophotometrically on a nanodrop at 230, 260 and 280 nm. The total RNA from single neurons were subjected to two rounds of linear T7 RNA Polymerase–driven transcription using MessageAmp™ II aRNA Amplification Kit from Ambion.

Microarray and bioinformatics analysis

A custom designed microarray [36] was used in gene expression experiments. Microarray analysis was carried out in collaboration with the genomics core at Sanford-Burnham Medical Research Institute, Orlando, Florida. RNA from three animals was used in each group (mature and old) for the analysis. The purity and integrity of the total RNA were analyzed on RNA Nano chip (Agilent Technologies) using Eukaryote Total RNA Nano series protocol. The total RNA was subjected to two rounds of linear IVT-amplification and labeled with Cy3-labeled CTP using Amino Allyl MessageAmp II Amplification kits (Ambion). The resulting Cy3 dye incorporated aRNA was quantified using ND-1000 spectrophotometer (Nano Dro Technologies) and 1.65 ug of labeled aRNA was hybridized on to Agilent’s 4 × 44 k format arrays (Agilent Technologies). After hybridization, the arrays were washed following the manufacturer’s protocol using Gene Expression Wash Pack (Agilent Technologies) and scanned using the Agilent C Scanner. The intensities of the scanned fluorescence images were extracted with Agilent Feature Extraction software version The data reported in this paper have been deposited in the Gene Expression Omnibus (GEO) database ( with the accession number GEO: GSE46618.

The mean signals were background corrected and transformed to the log2 scale. The data was then normalized between arrays by quantile approach. The empirical Bayes moderated t-statistics, which is implemented in the limma Bioconductor package [83], were used for differential expression detection. Genes with at least 2-fold changes and p-value 0.05 were considered as significant. The hierarchical clustering and other statistical analyses were performed using R/Bioconductor [84]. Principal component analysis (PCA) was performed with Partek Genomics Suite (Partek Inc. St. Louis, MO).

To ensure maximum accuracy for annotation of the signature list (1083 microarray probes), the probes were blasted against the assembled RNA transcripts ( First, a BLAST database was prepared. The resulting database was verified to contain the expected 1253635 sequences. Then, a BLASTN was performed against this database using the flat file containing the microarray probes. We chose the "tabular" output format to ensure easy manipulation programmatically.

The method described above produced 737 transcripts that matched the microarray probes (E-value < 1). Out of 737 transcripts, 572 were successfully blasted (tblastx) against nr database (E-value < 10) using the BLAST2GO application ( 247 transcripts were upregulated and 325 were downregulated in R15 mature versus R15 old neurons and these transcripts were studied by IPA analysis. To annotate the above 737 transcripts based on Gene Ontologies, the transcripts were blasted (blastx) against nr database (E-value < 1.0E-3). Mapping function in Blast2GO links different protein IDs in the BLAST hits to Gene Ontology database (GO consortium <>). The GO database contains several million functionally annotated genes for hundreds of species from public resources provided by the NCBI, PIR and GO. To ensure high quality annotation, only ontologies obtained from hits with E-value < 1.0E-6 are used for annotation and shown in the GO graphs. Annotation cutoff was set to 55 and GO weight was set to 5. Filter GO by taxonomy parameter was set to none.

Mapping and annotation was performed based on default parameters in the BLATS2GO application. Additionally, InterProScan function was utilized to improve the gene ontology analysis results. GO analysis for upregulated and downregulated hits shown in the Additional file 1 was performed using GO-Slim and combined graph function in BLAST2GO. The sequence filter was set to 5.

Ingenuity Pathway Analysis (IPA)

The data mining software, Ingenuity Pathway Analysis (IPA), was used for further analysis of the gene list obtained from the BLATS2GO analysis. The top blast hits accessions exported from Blast2GO were converted to Entrez Gene IDs using NCBI Batch Entrez tool, and then imported into IPA. IPA defaults to human however does have mouse and rat genes within it. Only if a human gene ortholog is not available, the mouse gene will be mapped. In our analysis, 190 human genes were mapped to the above gene list in IPA. This list was used for functional and network analysis using the Core Analysis feature in IPA. Red indicates upregulation and green indicates downregulation in R15 mature neurons relative to R15 old neurons.

Quantitative real time PCR (qPCR)

The genes identified by microarray analysis were validated by qPCR. Total RNA was isolated from single neurons as described earlier and was subjected to two rounds of linear IVT-amplification. 1 μg of RNA was used in 40 ul of reaction performed with qScript™ cDNA SuperMix from Quanta BioSciences according to the manufacturer's protocol. All amplifications were primed by pairs of chemically synthesized 18- to 24- mer oligonucleotides designed using freely available primer design software (Primer-3,) to generate target amplicons of 70–110 bp. Sequence information of the forward and reverse primers used are listed in Additional file 2: Table S5. The qPCR mixture had a total volume of 10 μl containing 2 μl of H2O, 2 μl of cDNA, 5 μl of 2X Master Mix, 1.0 μl of 10 μM (each) forward and reverse primer. The reaction was carried out in a 7900HT Fast Real-Time PCR System (Applied Biosystems Carlsbad, CA) under the following conditions: 95°C for 10 minutes, followed by 40 cycles of 95°C for 15 seconds, 60°C for 1 minute. There were five biological replicates and four technical replicates for each biological replicate. Quantification of the target transcripts was normalized to the Aplysia18S reference gene using the Pfaffl method [85]. Student t test was used to detect genes with statistically significant expression levels between mature and old single neurons where *P-value < 0.05, **P-value < 0.01, ***P-value < 0.001.

Availability of supporting data

All the microarray data is deposited at Gene Expression Omnibus (GEO) database ( with the accession number GEO: GSE46618. All the other supporting data (Additional file 1 and Additional file 2) are deposited at LabArchives ( in the folder named Puthanveettil-BMC Genomics).



Artificial Sea Water


Basic Local Alignment Search Tool


Central nervous system


cAMP Response Element-Binding protein


Gene Expression Omnibus


Gene ontology


In Vitro Transcription


Principal component analysis


Quantitative real time PCR


S6 kinase.


  1. Bordner KA, Kitchen RR, Carlyle B, George ED, Mahajan MC, Mane SM, Taylor JR, Simen AA: Parallel declines in cognition, motivation, and locomotion in aging mice: association with immune gene upregulation in the medial prefrontal cortex. Exp Gerontol. 2011, 46 (8): 643-659.

    PubMed Central  PubMed  Google Scholar 

  2. Nakamura TJ, Nakamura W, Yamazaki S, Kudo T, Cutler T, Colwell CS, Block GD: Age-Related Decline in Circadian Output. J Neurosci. 2011, 31 (28): 10201-10205. 10.1523/JNEUROSCI.0451-11.2011.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. Landis GN, Abdueva D, Skvortsov D, Yang J, Rabin BE, Carrick J, Tavaré S, Tower J: Similar gene expression patterns characterize aging and oxidative stress in Drosophila melanogaster. Proc Natl Acad Sci U S A. 2004, 101 (20): 7663-7668. 10.1073/pnas.0307605101.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  4. Lund J, Tedesco P, Duke K, Wang J, Kim SK, Johnson TE: Transcriptional profile of aging in C-elegans. Curr Biol. 2002, 12 (18): 1566-1573. 10.1016/S0960-9822(02)01146-6.

    Article  CAS  PubMed  Google Scholar 

  5. Pletcher SD, Macdonald SJ, Marguerie R, Certa U, Stearns SC, Goldstein DB, Partridge L: Genome-wide transcript profiles in aging and calorically restricted Drosophila melanogaster. Curr Biol. 2002, 12 (9): 712-723. 10.1016/S0960-9822(02)00808-4.

    Article  CAS  PubMed  Google Scholar 

  6. Southworth LK, Owen AB, Kim SK: Aging mice show a decreasing correlation of gene expression within genetic modules. Plos Genet. 2009, 5 (12): e1000776-10.1371/journal.pgen.1000776.

    Article  PubMed Central  PubMed  Google Scholar 

  7. Fraser HB, Khaitovich P, Plotkin JB, Paabo S, Eisen MB: Aging and gene expression in the primate brain. Plos Biol. 2005, 3 (9): 1653-1661.

    Article  CAS  Google Scholar 

  8. Sherwood CC, Gordon AD, Allen JS, Phillips KA, Erwin JM, Hof PR, Hopkins WD: Aging of the cerebral cortex differs between humans and chimpanzees. P Natl Acad Sci USA. 2011, 108 (32): 13029-13034. 10.1073/pnas.1016709108.

    Article  CAS  Google Scholar 

  9. Savvateeva EV, Popov AV, Kamyshev NG, Iliadi KG, Bragina JV, Heisenberg M, Kornhuber J, Riederer P: Age-dependent changes in memory and mushroom bodies in the Drosophila mutant vermilion deficient in the kynurenine pathway of tryptophan metabolism. Ross Fiziol Zh Im I M Sechenova. 1999, 85 (1): 167-183.

    CAS  PubMed  Google Scholar 

  10. Hertweck M, Hoppe T, Baumeister R: C-elegans, a model for aging with high-throughput capacity. Exp Gerontol. 2003, 38 (3): 345-346. 10.1016/S0531-5565(02)00208-5.

    Article  PubMed  Google Scholar 

  11. Papka R, Peretz B, Tudor J, Becker J: Age-dependent anatomical changes in an identified neuron in the CNS of Aplysia californica. J Neurobiol. 1981, 12 (5): 455-68. 10.1002/neu.480120505.

    Article  CAS  PubMed  Google Scholar 

  12. Jucker M, Ingram DK: Murine models of brain aging and age-related neurodegenerative diseases. Behav Brain Res. 1997, 85 (1): 1-25. 10.1016/S0166-4328(96)02243-7.

    Article  CAS  PubMed  Google Scholar 

  13. Cummings BJ, Head E, Ruehl W, Milgram NW, Cotman CW: The canine as an animal model of human aging and dementia. Neurobiol Aging. 1996, 17 (2): 259-268. 10.1016/0197-4580(95)02060-8.

    Article  CAS  PubMed  Google Scholar 

  14. Yu DH, Akalal DBG, Davis RL: Drosophila alpha/beta mushroom body neurons form a branch-specific, long-term cellular memory trace after spaced olfactory conditioning. Neuron. 2006, 52 (5): 845-855. 10.1016/j.neuron.2006.10.030.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Tomchik SM, Davis RL: Dynamics of learning-related cAMP signaling and stimulus integration in the Drosophila olfactory Pathway. Neuron. 2009, 64 (4): 510-521. 10.1016/j.neuron.2009.09.029.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Kandel ER: Neuroscience - The molecular biology of memory storage: a dialogue between genes and synapses. Science. 2001, 294 (5544): 1030-1038. 10.1126/science.1067020.

    Article  CAS  PubMed  Google Scholar 

  17. Barco A, Bailey CH, Kandel ER: Common molecular mechanisms in explicit and implicit memory. J Neurochem. 2006, 97 (6): 1520-1533. 10.1111/j.1471-4159.2006.03870.x.

    Article  CAS  PubMed  Google Scholar 

  18. Giles AC, Rankin CH: Behavioral and genetic characterization of habituation using Caenorhabditis elegans. Neurobiol Learn Mem. 2009, 92 (2): 139-146. 10.1016/j.nlm.2008.08.004.

    Article  PubMed  Google Scholar 

  19. Hill ES, Vasireddi SK, Bruno AM, Wang J, Frost WN: Variable neuronal participation in stereotypic motor programs. Plos One. 2012, 7 (7): e40579-10.1371/journal.pone.0040579.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  20. Popescu IR, Frost WN: Highly dissimilar behaviors mediated by a multifunctional network in the marine mollusk Tritonia diomedea. J Neurosci. 2002, 22 (5): 1985-1993.

    CAS  PubMed  Google Scholar 

  21. Castellucci V, Pinsker H, Kupfermann I, Kandel ER: Neuronal mechanisms of habituation and dishabituation of the gill-withdrawal reflex in Aplysia. Science. 1970, 167 (3926): 1745-1748. 10.1126/science.167.3926.1745.

    Article  CAS  PubMed  Google Scholar 

  22. Castellucci VF, Carew TJ, Kandel ER: Cellular analysis of long-term habituation of the gill-withdrawal reflex of Aplysia californica. Science. 1978, 202 (4374): 1306-1308. 10.1126/science.214854.

    Article  CAS  PubMed  Google Scholar 

  23. Peretz B, Romanenko A, Markesbery W: Functional history of two motor neurons and the morphometry of their neuromuscular junctions in the gill of Aplysia: evidence for differential aging. Proc Natl Acad Sci. 1984, 81 (13): 4232-4236. 10.1073/pnas.81.13.4232.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Peretz B, Ringham G, Wilson R: Age-diminished motor neuronal function of central neuron L7 in Aplysia. J Neurobiol. 1982, 13 (2): 141-151. 10.1002/neu.480130206.

    Article  CAS  PubMed  Google Scholar 

  25. Peretz B: Long-term change of viability of neuron functioning and its possible behavioral consequences in the adult Aplysia. J Physiol-Paris. 1988, 83 (3): 217-223.

    PubMed  Google Scholar 

  26. Skinner TL, Peretz B: Age sensitivity of Osmoregulation and of its neural correlates in Aplysia. Am J Physiol. 1989, 256 (4): R989-R996.

    CAS  PubMed  Google Scholar 

  27. McManus JM, Lu H, Chiel HJ: An in vitro preparation for eliciting and recording feeding motor programs with physiological movements in Aplysia californica. J Vis Exp. 2012, 70: e4320-

    PubMed  Google Scholar 

  28. Perrins R, Weiss KR: A cerebral central pattern generator in Aplysia and its connections with buccal feeding circuitry. J Neurosci. 1996, 16 (21): 7030-7045.

    CAS  PubMed  Google Scholar 

  29. Xin Y, Weiss KR, Kupfermann I: An identified interneuron contributes to aspects of six different behaviors in Aplysia. J Neurosci. 1996, 16 (16): 5266-5279.

    CAS  PubMed  Google Scholar 

  30. Moroz LL, Kohn AB: Do different neurons age differently? Direct genome-wide analysis of aging in single identified cholinergic neurons. Front Aging Neurosci. 2010, 2:

    Google Scholar 

  31. Sloan MA, Levenson J, Tran Q, Kerbeshian M, Block GD, Eskin A: Aging affects the ocular circadian pacemaker of Aplysia californica. J Biol Rhythms. 1999, 14 (2): 151-159. 10.1177/074873099129000542.

    Article  CAS  PubMed  Google Scholar 

  32. Fieber LA, Carlson SL, Capo TR, Schmale MC: Changes in D-aspartate ion currents in the Aplysia nervous system with aging. Brain Res. 2010, 1343: 28-36.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. Bailey CH, Castellucci VF, Koester J, Chen M: Behavioral changes in aging Aplysia: a model system for studying the cellular basis of age-impaired learning, memory, and arousal. Behav Neural Biol. 1983, 38 (1): 70-81. 10.1016/S0163-1047(83)90399-0.

    Article  CAS  PubMed  Google Scholar 

  34. Rattan KS, Peretz B: Age-dependent behavioral changes and physiological changes in identified neurons in Aplysia californica. J Neurobiol. 1981, 12 (5): 469-478. 10.1002/neu.480120506.

    Article  CAS  PubMed  Google Scholar 

  35. Hallahan BJ, Peretz B, Skinner TL: Increased age affects properties characterizing behavioral plasticity in freely behaving Aplysia. Neurobiol Aging. 1992, 13 (2): 217-225. 10.1016/0197-4580(92)90033-T.

    Article  CAS  PubMed  Google Scholar 

  36. Puthanveettil SV, Antonov I, Kalachikov S, Rajasethupathy P, Choi YB, Kohn AB, Citarella M, Yu F, Karl KA, Kinet M, et al: A strategy to capture and characterize the synaptic transcriptome. Proc Natl Acad Sci U S A. 2013, 110 (18): 7464-7469. 10.1073/pnas.1304422110.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  37. Peretz B, Srivatsan M: Differences in aging in two neural pathways: proposed explanations from the nervous system of Aplysia. Exp Gerontol. 1992, 27 (1): 83-97. 10.1016/0531-5565(92)90031-T.

    Article  CAS  PubMed  Google Scholar 

  38. Gerdes R, Fieber LA: Life history and aging of captive-reared California sea hares (Aplysia californica). J Am Assoc Lab Anim Sci. 2006, 45 (1): 40-47.

    CAS  PubMed  Google Scholar 

  39. Stommes D, Fieber LA, Beno C, Gerdes R, Capo TR: Temperature effects on growth, maturation, and lifespan of the california sea hare (Aplysia californica). Contemp Top Lab Anim Sci. 2005, 44 (3): 31-35.

    CAS  PubMed  Google Scholar 

  40. Frazier WT, Kandel ER, Kupfermann I, Waziri R, Coggeshall RE: Morphological and functional properties of identified neurons in the abdominal ganglion of Aplysia Californica. J Neurophysiol. 1967, 30 (6): 1288-1351.

    Google Scholar 

  41. Alevizos A, Weiss KR, Koester J: Synaptic actions of identified peptidergic neuron R15 in Aplysia. I. Activation of respiratory pumping. J Neurosci. 1991, 11 (5): 1263-1274.

    CAS  PubMed  Google Scholar 

  42. Alevizos A, Weiss KR, Koester J: Synaptic actions of identified peptidergic neuron R15 in Aplysia. II. Contraction of pleuroabdominal connectives mediated by motoneuron L7. J Neurosci. 1991, 11 (5): 1275-1281.

    CAS  PubMed  Google Scholar 

  43. Alevizos A, Weiss KR, Koester J: Synaptic actions of identified peptidergic neuron R15 in Aplysia. III. Activation of the large hermaphroditic duct. J Neurosci. 1991, 11 (5): 1282-1290.

    CAS  PubMed  Google Scholar 

  44. Benson JA, Levitan IB: Serotonin increases an anomalously rectifying K + current in the Aplysia neuron R15. Proc Natl Acad Sci U S A. 1983, 80 (11): 3522-3525. 10.1073/pnas.80.11.3522.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  45. Gospe SM, Wilson WA: Dopamine inhibits burst-firing of neurosecretory cell R 15 in Aplysia californica: establishment of a dose–response relationship. J Pharmacol Exp Ther. 1980, 214 (1): 112-118.

    PubMed  Google Scholar 

  46. Canavier CC, Clark JW, Byrne JH: Simulation of the bursting activity of neuron R15 in Aplysia: role of ionic currents, calcium balance, and modulatory transmitters. J Neurophysiol. 1991, 66 (6): 2107-2124.

    CAS  PubMed  Google Scholar 

  47. Marcus EA, Carew TJ: Development and modulation of endogenous bursting in identified neuron R15 of juvenile Aplysia. J Neurobiol. 1991, 22 (4): 418-429. 10.1002/neu.480220409.

    Article  CAS  PubMed  Google Scholar 

  48. Brown RO, Mayeri E: Central actions of R15, a putative peptidergic neuron in Aplysia. J Neurobiol. 1987, 18 (1): 3-13. 10.1002/neu.480180103.

    Article  CAS  PubMed  Google Scholar 

  49. Moroz LL, Edwards JR, Puthanveettil SV, Kohn AB, Ha T, Heyland A, Knudsen B, Sahni A, Yu F, Liu L, et al: Neuronal transcriptome of Aplysia: neuronal compartments and circuitry. Cell. 2006, 127 (7): 1453-1467. 10.1016/j.cell.2006.09.052.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  50. Giannakou ME, Goss M, Junger MA, Hafen E, Leevers SJ, Partridge L: Long-lived Drosophila with overexpressed dFOXO in adult fat body. Science. 2004, 305 (5682): 361-10.1126/science.1098219.

    Article  CAS  PubMed  Google Scholar 

  51. Ogg S, Paradis S, Gottlieb S, Patterson GI, Lee L, Tissenbaum HA, Ruvkun G: The Fork head transcription factor DAF-16 transduces insulin-like metabolic and longevity signals in C. elegans. Nature. 1997, 389 (6654): 994-999. 10.1038/40194.

    Article  CAS  PubMed  Google Scholar 

  52. Ruzankina Y, Pinzon-Guzman C, Asare A, Ong T, Pontano L, Cotsarelis G, Zediak VP, Velez M, Bhandoola A, Brown EJ: Deletion of the developmentally essential gene ATR in adult mice leads to age-related phenotypes and stem cell loss. Cell stem cell. 2007, 1 (1): 113-126. 10.1016/j.stem.2007.03.002.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  53. Kim DJ, Yi SM, Lee SY, Kang HS, Choi YH, Song YW, Park SC: Association between the MLH1 gene and longevity. Human genetics. 2006, 119 (3): 353-354. 10.1007/s00439-006-0148-7.

    Article  PubMed  Google Scholar 

  54. Fraser HB, Khaitovich P, Plotkin JB, Pääbo S, Eisen MB: Aging and gene expression in the primate brain. PLoS Biol. 2005, 3 (9): e274-10.1371/journal.pbio.0030274.

    Article  PubMed Central  PubMed  Google Scholar 

  55. Lund PK, Hoyt EC, Bizon J, Smith DR, Haberman R, Helm K, Gallagher M: Transcriptional mechanisms of hippocampal aging. Exp Gerontol. 2004, 39 (11–12): 1613-1622.

    Article  CAS  PubMed  Google Scholar 

  56. Seroude L: Differential gene expression and aging. ScientificWorldJournal. 2002, 2: 618-631.

    Article  CAS  PubMed  Google Scholar 

  57. Schimanski LA, Barnes CA: Neural protein synthesis during aging: effects on plasticity and memory. Front Aging Neurosci. 2010, 2:

    Google Scholar 

  58. Makrides SC: Protein synthesis and degradation during aging and senescence. Biol Rev Camb Philos Soc. 1983, 58 (3): 343-422. 10.1111/j.1469-185X.1983.tb00394.x.

    Article  CAS  PubMed  Google Scholar 

  59. Cohen TE, Kaplan SW, Kandel ER, Hawkins RD: A simplified preparation for relating cellular events to behavior: mechanisms contributing to habituation, dishabituation, and sensitization of the Aplysia gill-withdrawal reflex. J Neurosci. 1997, 17 (8): 2886-2899.

    CAS  PubMed  Google Scholar 

  60. Kupfermann I, Castellucci V, Pinsker H, Kandel E: Neuronal correlates of habituation and dishabituation of the gill-withdrawal reflex in Aplysia. Science. 1970, 167 (3926): 1743-1745. 10.1126/science.167.3926.1743.

    Article  CAS  PubMed  Google Scholar 

  61. Harries LW, Hernandez D, Henley W, Wood AR, Holly AC, Bradley-Smith RM, Yaghootkar H, Dutta A, Murray A, Frayling TM, et al: Human aging is characterized by focused changes in gene expression and deregulation of alternative splicing. Aging cell. 2011, 10 (5): 868-878. 10.1111/j.1474-9726.2011.00726.x.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  62. Vijg J, Suh Y: Functional genomics of ageing. Mech Ageing Dev. 2003, 124 (1): 3-8. 10.1016/S0047-6374(02)00163-X.

    Article  PubMed  Google Scholar 

  63. Kadakkuzha BM, Puthanveettil SV: Genomics and proteomics in solving brain complexity. Mol BioSyst. 2013, 9 (7): 1821-1807.

    Article  Google Scholar 

  64. Weindruch R, Kayo T, Lee CK, Prolla TA: Microarray profiling of gene expression in aging and its alteration by caloric restriction in mice. J Nutr. 2001, 131 (3): 918S-923S.

    CAS  PubMed  Google Scholar 

  65. Weindruch R, Kayo T, Lee CK, Prolla TA: Gene expression profiling of aging using DNA microarrays. Mech Ageing Dev. 2002, 123 (2–3): 177-193.

    Article  CAS  PubMed  Google Scholar 

  66. Yankner BA, Lu T, Loerch P: The aging brain. Annu Rev Pathol. 2008, 3: 41-66. 10.1146/annurev.pathmechdis.2.010506.092044.

    Article  CAS  PubMed  Google Scholar 

  67. Benito E, Barco A: CREB's control of intrinsic and synaptic plasticity: implications for CREB-dependent memory models. Trends Parasitol. 2010, 33 (5): 230-240.

    CAS  Google Scholar 

  68. Barco A, Marie H: Genetic approaches to investigate the role of CREB in neuronal plasticity and memory. Cell Mol Neurobiol. 2011, 44 (3): 330-349. 10.1007/s12035-011-8209-x.

    Article  CAS  Google Scholar 

  69. Barco A, Alarcon JM, Kandel ER: Expression of constitutively active CREB protein facilitates the late phase of long-term potentiation by enhancing synaptic capture. Cell. 2002, 108 (5): 689-703. 10.1016/S0092-8674(02)00657-8.

    Article  CAS  PubMed  Google Scholar 

  70. Suzuki A, Fukushima H, Mukawa T, Toyoda H, Wu LJ, Zhao MG, Xu H, Shang Y, Endoh K, Iwamoto T, et al: Upregulation of CREB-mediated transcription enhances both short- and long-term memory. J Neurosci. 2011, 31 (24): 8786-8802. 10.1523/JNEUROSCI.3257-10.2011.

    Article  CAS  PubMed  Google Scholar 

  71. Morris KA, Gold PE: Age-related impairments in memory and in CREB and pCREB expression in hippocampus and amygdala following inhibitory avoidance training. Mech Ageing Dev. 2012, 133 (5): 291-299. 10.1016/j.mad.2012.03.004.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  72. Shi Y, Venkataraman SL, Dodson GE, Mabb AM, LeBlanc S, Tibbetts RS: Direct regulation of CREB transcriptional activity by ATM in response to genotoxic stress. Proc Natl Acad Sci U S A. 2004, 101 (16): 5898-5903. 10.1073/pnas.0307718101.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  73. Keleshian VL, Modi HR, Rapoport SI, Rao JS: Aging is associated with altered inflammatory, arachidonic acid cascade, and synaptic markers, influenced by epigenetic modifications, in the human frontal cortex. J Neurochem. 2013, 125 (1): 63-73. 10.1111/jnc.12153.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  74. Saura CA, Valero J: The role of CREB signaling in Alzheimer's disease and other cognitive disorders. Rev Neurosci. 2011, 22 (2): 153-169.

    Article  CAS  PubMed  Google Scholar 

  75. Johnson SC, Rabinovitch PS, Kaeberlein M: mTOR is a key modulator of ageing and age-related disease. Nature. 2013, 493 (7432): 338-345. 10.1038/nature11861.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  76. Fenton TR, Gout IT: Functions and regulation of the 70 kDa ribosomal S6 kinases. Int J Biochem Cell Biol. 2011, 43 (1): 47-59. 10.1016/j.biocel.2010.09.018.

    Article  CAS  PubMed  Google Scholar 

  77. Meyuhas O, Dreazen A: Chapter 3 ribosomal protein S6 kinase: from TOP mRNAs to cell size. Prog Mol Biol Transl Sci. 2009, 90: 109-53.

    Article  CAS  PubMed  Google Scholar 

  78. Zoncu R, Efeyan A, Sabatini DM: mTOR: from growth signal integration to cancer, diabetes and ageing. Nat Rev Mol Cell Bio. 2011, 12 (1): 21-35. 10.1038/nrm3025.

    Article  CAS  Google Scholar 

  79. Selman C, Tullet JMA, Wieser D, Irvine E, Lingard SJ, Choudhury AI, Claret M, Al-Qassab H, Carmignac D, Ramadani F, et al: Ribosomal protein S6 kinase 1 signaling regulates mammalian life span. Science. 2009, 326 (5949): 140-144. 10.1126/science.1177221.

    Article  CAS  PubMed  Google Scholar 

  80. Khaitovich P, Muetzel B, She X, Lachmann M, Hellmann I, Dietzsch J, Steigele S, Do HH, Weiss G, Enard W, et al: Regional patterns of gene expression in human and chimpanzee brains. Genome Res. 2004, 14 (8): 1462-1473. 10.1101/gr.2538704.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  81. Evans SJ, Choudary PV, Vawter MP, Li J, Meador-Woodruff JH, Lopez JF, Burke SM, Thompson RC, Myers RM, Jones EG, et al: DNA microarray analysis of functionally discrete human brain regions reveals divergent transcriptional profiles. Neurobiol Dis. 2003, 14 (2): 240-250. 10.1016/S0969-9961(03)00126-8.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  82. Leal SL, Yassa MA: Perturbations of neural circuitry in aging, mild cognitive impairment, and Alzheimer's disease. Ageing Res Rev. 2013, 12 (3): 823-831. 10.1016/j.arr.2013.01.006.

    Article  PubMed Central  PubMed  Google Scholar 

  83. Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: Article 3-

    Google Scholar 

  84. Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, et al: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5 (10): R80-10.1186/gb-2004-5-10-r80.

    Article  PubMed Central  PubMed  Google Scholar 

  85. Pfaffl MW: A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001, 29 (9): e45-10.1093/nar/29.9.e45.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references


We thank Subramaniam Shyamala Govindarajan and Jian-Liang Li from the Analytical Genomics Core Facility at Sanford Burnham Medical Institute (Lake Nona) for performing microarray and analyzing array data. We thank the funding support from Whitehall Foundation, NIH grant 1 R21 MH096258-01A1 and startup funds from The Scripps Research Institute for carrying out this work.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Sathyanarayanan V Puthanveettil.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

SP, BK and KA conceived of the study. TC generated aged animals. BK, MS and SP wrote the article. MS & AC analyzed the microarray data and conducted GO and IPA analysis. All authors have read and approved the manuscript for publication.

Electronic supplementary material


Additional file 1 Figure S1: Relative expression of CREB1 in mature and old abdominal ganglia. CREB1 mRNA levels in the mature and old abdominal ganglia were quantified and normalized to 18S rRNA internal standard. CREB1 do not show any significant change in expression with aging in the whole abdominal ganglia. The error bars indicate standard error of the mean (SEM). Student’s t test determined statistical significance where P-value < 0.05 was considered significant. Figure S2. Statistical analysis of microarray data. (A) Scatter plots and correlation coefficients for each microarray experiment pair R15 mature and old single neurons. Mature #2 Rep1 and Old #3 Rep1 represent the replicates within each slide (internal replicates). (B) Principle Component Analysis (PCA) mapped scatter plot of technical replicates: The global gene expression profiles of the R15 neurons from both mature and old Aplysia were analyzed by PCA. The figure represented the first two principal components of microarray analysis data (PC1, PC2 and PC3) in X, Y and Z respectively. (Mature Red # 1–3, Old blue# 1–3). Figure S3. Identification of biological network functions by Ingenuity Pathway Analysis. A) Network 2: Top functions of the genes were related to cell death and survival, cellular function maintenance and embryonic development. (B) Network 3: This network has 15 genes differentially expressed and the biological functions associated with network 3 are neurological diseases, developmental and hereditary disorders. Node and edge symbols are shown along with network diagram. Red shows genes that are upregulated and green indicates downregulation of the genes. Uncolored notes represent genes that were not identified in our data set as differentially expressed in our analysis, but relevant to this network. (PPTX 3 MB)


Additional file 2 Table S1: Differentially expressed probe list from the microarray analysis of mature and old single R15 neuron RNA. Table S2. List of annotated Aplysia genes from RNA seq analysis. Table S3. Gene Ontology classification of differentially expressed genes from mature R15 neuron. TableS4. Ingenuity Pathway Analysis (IPA) of differentially expressed genes from mature R15 neuron. Table S5. Primer sequences used for qPCR analysis. (XLSX 5 MB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Kadakkuzha, B.M., Akhmedov, K., Capo, T.R. et al. Age-associated bidirectional modulation of gene expression in single identified R15 neuron of Aplysia. BMC Genomics 14, 880 (2013).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Aplysia californica
  • Bi-directional gene regulation
  • Single neuron transcriptome
  • R15 neuron
  • Neuronal circuitry
  • Signaling networks
  • Aging
  • CREB
  • S6 Kinase