Age-associated bidirectional modulation of gene expression in single identified R15 neuron of Aplysia
© Kadakkuzha et al.; licensee BioMed Central Ltd. 2013
Received: 17 May 2013
Accepted: 5 December 2013
Published: 14 December 2013
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.
KeywordsAplysia californica Bi-directional gene regulation Single neuron transcriptome R15 neuron Neuronal circuitry Signaling networks Aging CREB S6 Kinase
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 [3–5]. 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 [8–13]. 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 [14–18].
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 [23–26]. A key advantage of using Aplysia is that specific behavior is mediated by a simple circuitry [21, 22, 27–29]. 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 . Significant changes in gene expression in old R2 neurons, and a unique set of genes differentially expressed during aging were identified .
Senescence in Aplysia is characterized by declining weight, reproduction (egg laying), reflex response and finally death in 11–13 months . 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 [33–35] are associated with aging in Aplysia. At the macromolecular level aging associated changes include transcriptional and translational changes  in response to circadian rhythm , and the change in electrophysiological properties in identified neurons .
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 [36–39], neuronal circuitry functions and cellular basis of behavior . R15 has been a subject of extensive electrophysiological studies due to its role in egg laying [41–43], neurochemical modulation [44, 45].
Using custom designed microarrays  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 . 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 . 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 [45–47], and is implicated in reproduction, and osmoregulation [41–43, 48].
Age associated differentially expressed genes in mature and old R15 single neurons
p38 MAP kinase
PKA type II regulatory subunit
GO analysis of differentially expressed genes in R15
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
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
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, 62–65]. 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 . 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 [67–69], long-term memory  and aging [71, 72]. Recently CREB has been shown to be involved in human aging , cognitive and neurodegenerative disorders . In mice, CREB regulate aging related genes such as the A-T-mutated gene (ATM) . 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 . 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 . 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 . 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 [80–82].
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  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 10.7.3.1. The data reported in this paper have been deposited in the Gene Expression Omnibus (GEO) database (http://www.ncbi.nlm.nih.gov/geo) 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 , 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 . 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 (http://www.aplysiagenetools.org). 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 (http://www.blast2go.com). 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 <http://www.geneontology.org>). 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 . 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 (http://www.ncbi.nlm.nih.gov/geo) with the accession number GEO: GSE46618. All the other supporting data (Additional file 1 and Additional file 2) are deposited at LabArchives (http://www.labarchives.com/) 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
In Vitro Transcription
Principal component analysis
Quantitative real time PCR
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.
- 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 CentralPubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- Skinner TL, Peretz B: Age sensitivity of Osmoregulation and of its neural correlates in Aplysia. Am J Physiol. 1989, 256 (4): R989-R996.PubMedGoogle Scholar
- 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-PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Seroude L: Differential gene expression and aging. ScientificWorldJournal. 2002, 2: 618-631.View ArticlePubMedGoogle Scholar
- Schimanski LA, Barnes CA: Neural protein synthesis during aging: effects on plasticity and memory. Front Aging Neurosci. 2010, 2:Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Vijg J, Suh Y: Functional genomics of ageing. Mech Ageing Dev. 2003, 124 (1): 3-8. 10.1016/S0047-6374(02)00163-X.View ArticlePubMedGoogle Scholar
- Kadakkuzha BM, Puthanveettil SV: Genomics and proteomics in solving brain complexity. Mol BioSyst. 2013, 9 (7): 1821-1807.View ArticleGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Yankner BA, Lu T, Loerch P: The aging brain. Annu Rev Pathol. 2008, 3: 41-66. 10.1146/annurev.pathmechdis.2.010506.092044.View ArticlePubMedGoogle Scholar
- 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.Google Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Saura CA, Valero J: The role of CREB signaling in Alzheimer's disease and other cognitive disorders. Rev Neurosci. 2011, 22 (2): 153-169.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.