Memory acquisition and retrieval impact different epigenetic processes that regulate gene expression
© Peixoto et al.; licensee BioMed Central Ltd. 2015
Published: 26 May 2015
A fundamental question in neuroscience is how memories are stored and retrieved in the brain. Long-term memory formation requires transcription, translation and epigenetic processes that control gene expression. Thus, characterizing genome-wide the transcriptional changes that occur after memory acquisition and retrieval is of broad interest and importance. Genome-wide technologies are commonly used to interrogate transcriptional changes in discovery-based approaches. Their ability to increase scientific insight beyond traditional candidate gene approaches, however, is usually hindered by batch effects and other sources of unwanted variation, which are particularly hard to control in the study of brain and behavior.
We examined genome-wide gene expression after contextual conditioning in the mouse hippocampus, a brain region essential for learning and memory, at all the time-points in which inhibiting transcription has been shown to impair memory formation. We show that most of the variance in gene expression is not due to conditioning and that by removing unwanted variance through additional normalization we are able provide novel biological insights. In particular, we show that genes downregulated by memory acquisition and retrieval impact different functions: chromatin assembly and RNA processing, respectively. Levels of histone 2A variant H2AB are reduced only following acquisition, a finding we confirmed using quantitative proteomics. On the other hand, splicing factor Rbfox1 and NMDA receptor-dependent microRNA miR-219 are only downregulated after retrieval, accompanied by an increase in protein levels of miR-219 target CAMKIIγ.
We provide a thorough characterization of coding and non-coding gene expression during long-term memory formation. We demonstrate that unwanted variance dominates the signal in transcriptional studies of learning and memory and introduce the removal of unwanted variance through normalization as a necessary step for the analysis of genome-wide transcriptional studies in the context of brain and behavior. We show for the first time that histone variants are downregulated after memory acquisition, and splicing factors and microRNAs after memory retrieval. Our results provide mechanistic insights into the molecular basis of cognition by highlighting the differential involvement of epigenetic mechanisms, such as histone variants and post-transcriptional RNA regulation, after acquisition and retrieval of memory.
Genome-wide differential gene expression analysis is widely used in discovery-based studies in biology and medicine. The question of how variability impacts reproducibility of genome-wide results has been subject to extensive research . It is known that unwanted variation is often a confounding factor. Unwanted variation refers to other factors that influence the observed gene expression levels besides the one of interest. A typical example is a batch effect, which can occur when some samples are processed differently than others. Batch effects are not the only source of unwanted variance. Unwanted variance in microarrays arising from technical aspects of the methodology is removed using normalization methods such as RMA . The amount of unwanted biological variance depends on the question of interest and is influenced by factors such as heterogeneity in cell-types, variability in responsiveness to stimulus between biological replicates and the simultaneous presence of other stimuli other than the one of interest, such as time of day or other environmental variables. All of these factors are present when studying gene expression in the brain in vivo and are often hard to control. Thus, in the context of brain and behavior a major challenge is to normalize unwanted variation to minimize false discoveries, increase resolution and maximize the potential of discovery-based approaches to contribute biological insight.
Several aspects of brain function are linked to transcriptional changes. Long-term memory formation, for example, is known to require transcription, protein synthesis and epigenetic processes that regulate gene expression [3–8]. How memories are stored and retrieved in the brain is a fundamental question in neuroscience. Thus, characterizing genome-wide the transcriptional changes that occur after memory acquisition and retrieval is of broad interest and importance. Research has shown that there are "sensitive periods" after memory acquisition during which inhibiting mRNA or protein synthesis impairs memory formation. Using contextual fear conditioning as a task, these windows occur immediately or 4 hours after acquisition for memory tested 24 hours later [9, 10], or 12 hours after training for memory tested a week later . Processes that follow retrieval of the memory trace (extinction or reconsolidation) also require transcription and protein synthesis [12–15]. Several studies have used genome-wide approaches such as microarrays to describe changes in coding and non-coding gene expression after memory acquisition or synaptic activity [16–20]. These studies have led to the identification of some genes relevant for memory formation, such as c-rel or miR-182 [18, 21]. It remains unclear to what degree the variety of other stimuli experienced by the brain in vivo hinders reproducibility and limits the applicability of genome-wide technologies to the study of the brain and behavior.
Here, we examined genome-wide gene expression after contextual conditioning in the mouse hippocampus, a brain region essential for memory formation, during all the established sensitive periods for transcriptional inhibition. We show that most of the variance in gene expression is not due to conditioning and that by removing unwanted variance through additional normalization we are able provide novel biological insights. We show for the first time that histone variants are downregulated after memory acquisition, and splicing factors and microRNAs after retrieval. Our results provide mechanistic insights into the molecular basis of cognition by highlighting the differential involvement of epigenetic mechanisms, such as histone variants and post-transcriptional RNA regulation, after acquisition and retrieval of memory.
Results and discussion
We examined genome-wide changes of gene expression in adult, male C57BL6/J mice following a contextual fear conditioning paradigm (FC), a form of learning in which an aversive stimulus (e.g a shock, US) is associated with a neutral context (CS). Re-exposure to the context triggers retrieval of the memory for the context-shock association (CS-US), which is quantified as freezing in mice. FC is highly reproducible among individuals, requiring a single exposure to the CS-US pairing to learn. In addition the timeline of sensitivity for transcriptional inhibition is established, making it an ideal learning task for our genomic study. FC is known to require the hippocampus, a brain region essential for long-term memory formation. Hippocampal tissue was collected at the established sensitive periods for transcriptional inhibition during memory consolidation: 30 minutes (FC30'), 4 hours (FC4), 12 hours (FC12) or 24 hours after FC (FC24), as well as 30 minutes after memory retrieval (RT30'). Animals in the retrieval group showed typical learning of the task, with average freezing of 55% (+/- 10%) after re-exposure to the context. Animals that were handled but not trained were dissected at the same time of day to control for circadian variation in gene expression (CC30', CC4 and CC12). RNA from nine animals per group representing nine independent behavioral experiments conducted at the same time of day (72 samples) was hybridized simultaneously to an Affymetrix gene Titan Mouse 1.1 gene-EST microarray. Pairing of the CS and US was necessary to evaluate proper retrieval of the memory trace. We have previously shown that genome-wide gene expression changes in the hippocampus 30 minutes after exposure to the CS alone are not different from those after CS+US pairing [17, 22]. These results suggest that in this brain region the US alone does not produce significant differences in gene expression and that differences in gene expression observed are likely due to the encoding of the spatial memory for the context alone. Therefore, we did not include animals exposed to only the CS or US in the design.
Individual variability and circadian time are the biggest drivers of variance in gene expression in the hippocampus in vivo
To explore the main sources of variance in the data, we first performed a principal component analysis (PCA)  (Additional file 1). None of the first three principal components (PC), which account for over 65% of the variance, capture the response to the treatment. Understanding the dominant sources of variance is critical to accurately assess the effects of learning in gene expression. The first and second principal components (PC1 and PC2) represent unusual variability in gene expression in individual mice (Additional file 1A and 1B). Interestingly, the contribution of PC1 and PC2 to variance in gene expression is correlated for a subset of genes (Additional file 1C) and may represent the same biological process. Functional annotation analysis of the genes with correlated scores between PC1 and PC2 shows that they represent response to olfactory stimuli, specifically pheromones (Additional file 2), suggesting that individual responses to, or environmental variations in olfactory stimuli are the strongest drivers of differences in gene expression.
There is substantial evidence that memory consolidation is affected by circadian time [24–26]. However, the circadian influence on genome-wide gene expression in the mouse hippocampus is poorly understood. Published genome-wide studies of gene expression in response to activity in the brain that include several time points often do not include controls for time of day [16, 18], and thus it is hard to distinguish the effect of neuronal activity from the circadian effect in such studies. In our analysis, the third principal component (PC3) reveals that circadian time has a strong influence on hippocampal gene expression. Additional file 1D shows that the effect of circadian is similar to the effect of learning at the second (4 h) and third (12 h) sensitive periods, but not immediately after (30 minutes) memory acquisition or retrieval. To characterize gene expression changes in the mouse hippocampus due to circadian time, we compared genome-wide gene expression among our three control time-points (CC30', CC4 and CC12). The greatest number of differences in gene expression was detected between CC30 and CC12. These time points correspond to Zeitgeber times 3 and 15 (ZT3 and ZT15), one time-point during the light phase and one during the dark phase. 1067 probe sets, corresponding to 1019 known genes, were differentially expressed at a false discovery rate (fdr) <0.1. (Additional file 3). To determine whether our dataset corresponds to known circadian regulated genes, we compared it to those genes known to oscillate in the mouse liver according to Hughes and colleagues . Of the 1019 genes regulated by time of day in the hippocampus, 198 oscillate with 24 hour periodicity in the mouse liver (Additional file 4A). This overlap is significantly higher than expected by chance (Fisher's exact p-value of 0.004) based on an expected 15% overlap between any 2 mouse tissues as reported in the mouse gene atlas ). Our dataset is the first genome-wide dataset describing the effects of circadian time on gene expression in the hippocampus. Genes differentially expressed between ZT3 and ZT15 in the hippocampus include Per1, Per2 and Per3, which are known circadian oscillators . Interestingly, genes usually thought to be associated with memory formation, such as Arc, Bdnf [29, 30], CBP [31–35] and p300 [36, 37] also show circadian changes in expression (Additional file 4B).
Memory acquisition and retrieval induce similar, but distinct, genome-wide changes in gene expression 30 minutes after exposure
Memory acquisition and retrieval downregulate different epigenetic processes that modify gene expression
An interesting observation in Figure 3A is that the lower-left quadrant is almost empty, showing little overlap between downregulated genes. Accordingly, genes downregulated after acquisition and retrieval show no overlap in function. Chromatin assembly is downregulated after acquisition (Additional file 7A), exemplified by histone 2A isoforms Hist1h2af, Hist2h2ab, Hist1h2ao and Hist2h2aa1. RNA processing is downregulated after retrieval (Additional data file 7B), exemplified by splicing factors Prpf38b and Rbfox1, and spliceosome kinase Srpk2. Downregulation of genes involved in chromatin assembly observed in our microarray after acquisition is driven by downregulation of Histone 2A isoforms. Based on closer inspection of the probe-level data, Hist2h2ab emerged as the H2A gene most likely regulated by during memory consolidation. Greater than 2-fold downregulation of expression of Hist2h2ab, herein referred to as H2AB, was confirmed by qPCR in an independent cohort of animals (Figure 3B).
Memory retrieval had a bigger impact on processes that regulate RNA processing. The downregulation of splicing factor Rbfox1 was also replicated by qPCR in an independent cohort (Figure 3C). To our knowledge, this is the first time that Rbfox1 has been reported to be regulated by behavior. Rbfox1 is an important regulator of both splicing and transcription in brain development  that controls neuronal excitability . Clinically, Rbfox1 deletion is correlated with developmental delays, learning disabilities and autistic-like features . The results of the functional annotation analysis (Figure 3A), the validation of the down-regulation of Rbfox1 by qPCR (Figure 3C) and the differences in alternative splicing of Per1 observed between FC30' and RT30' (Figure 2B), suggest that post-transcriptional regulation may be of particular importance following memory retrieval. Splicing regulators have been previously reported to change expression during memory formation . However, alternative splicing following memory formation at the genome-wide level has not been previously examined. Greater efforts in future experiments will be directed at understanding the regulation of different spliced isoforms after acquisition and retrieval.
Memory consolidation and retrieval differentially regulate non-coding RNA expression
We also evaluated the regulation of 1,229 non-coding transcripts with well-established annotation. A summary of non-coding RNAs regulated either at FC30' and RT30' can be found in Additional file 11. MicroRNA genes miR-212, miR-132 and miR-219 were selected for further validation. MiR-212 and miR-132 are CREB-dependent microRNAs derived from the same precursor that are induced by LTP  and play an important role in neuronal plasticity . MiR-219 expression is dependent on the activity of NMDA receptors , which play an essential role in the acquisition of spatial memories in the hippocampus [62, 63]. In addition, we selected three non-coding RNAs whose role in brain function has not been studied previously for further validation: two variants of a small nucleolar RNA (snoRNAs Snord14d and Snord14e) and miR-410, one of the microRNAs contained within the Mirg imprinted non-coding RNA cluster. Snord14e and Snord14d represent some of the highest fold changes seen in our microarray while Mirg is strongly expressed in the brain during development . Mirg contains at least 13 microRNAs (UCSC genome browser) including miR-410 whose expression is known to be specific to the central nervous system .
Downregulation of miR-219 was only found to be significant after retrieval (Figure 5A), supporting our previous observation that retrieval downregulates genes involved in RNA processing. MiR-219 is known to regulate protein levels of CAMKIIγ . To investigate if the observed reduction of miR-219 had a functional effect on protein levels of CAMKIIγ we performed western-blots in an independent cohort of animals (n = 10 per group) and found that CAMKIIγ protein levels are indeed significantly increased (p < 0.05) after memory retrieval (Figure 5B). It is not clear if the molecular changes we observe after retrieval correspond to reconsolidation or extinction of the memory trace. Based on available literature [71–74], the observed down-regulation of miR-219 and up-regulation of CAMKIIγ is consistent with the hypothesis that a single brief re-exposure to the context may inhibit NMDAR activity while maintaining or even increasing CAMKII signaling and thus promote memory reconsolidation while inhibiting extinction. Although it is likely that CAMKIIγ plays an important role in CAMKII mediated signaling, its specific function remains unknown.
Our study characterizes gene expression genome-wide, both protein coding and non-coding, at several time-points during memory consolidation and following retrieval of memory. We show that training is not the main source of variance in gene expression. We introduce the removal of unwanted variance though normalization to the study of transcriptional changes genome-wide in the context of brain and behavior. Using this approach, we successfully identify novel gene expression changes following contextual conditioning as well as reproduce the majority of the previously reported changes. The largest changes in gene expression related to memory formation are observed 30 minutes after memory acquisition and 30 minutes post retrieval. We identify a novel activity-dependent histone variant, H2AB, and show that is downregulated following memory acquisition. We point to several gene candidates that have not been previously described to be involved in learning and memory, including transcriptional regulators Sox18, Btg2 and Sik1, and splicing factor Rbfox1. Finally, we examine genome-wide non-coding RNA regulation following memory acquisition and retrieval, pointing to a likely important role of microRNAs miR-132, miR-212, miR-410 and snoRNAs Snord14d and Snord14e in posttranscriptional regulation during both processes as well as a specific role for and miR-219 and its target CAMKIIγ after retrieval. Epigenetic mechanisms that regulate gene expression have been shown to be essential to long-term memory formation. Our study underlines the importance of two currently understudied epigenetic processes to memory storage and retrieval: histone variants and post-transcriptional RNA regulation; the study of which will expand our knowledge of the molecular mechanisms by which the brain maintains long-lasting changes induced by experience.
C57BL/6J adult male mice (2 months of age) were obtained from Jackson Laboratories and housed individually for a week on a 12 hr./12 hr. light/dark schedule with lights on at 7 am (Zeitgeber time (ZT) 0). Food and water were available ad libitum throughout the experiment. Each animal was handled daily for 3 days prior to contextual fear conditioning (FC). Handling consisted of manipulation of the animals for 1-2 minutes per mouse in the same room as the experimental setting without exposure to the context. The conditioning protocol entailed a single 2-second, 1.5-mA foot shock, terminating at 2.5 minutes after placement of the mouse in the chamber, starting at 10 am (ZT3) daily. Plexiglass operant chambers (Med-Associates, 31.8 × 25.4 × 26.7 cm) housed in sound-attenuating boxes equipped with individual fans and lights were used for all fear conditioning experiments and controls. The floor consisted of stainless-steel grid rods 3.2 mm in diameter, spaced 0.5 cm apart. Footshock was delivered by a Med-Associates solid-state shock source and grid floor scrambler that delivered a constant current (1.5 mA). Conditioning was quantified by measuring freezing behavior, using automated scoring software (Clever Systems). For microarray experiments, hippocampal dissections were performed immediately following the behavioral treatment, and alternated between FC and control animals. Tissue was collected at 30 minutes (FC30'), 4 hours (FC4), 12 hours (FC12) or 24 hours after FC (FC24) as well as 30 minutes after testing for retrieval of the memory (RT30'). Testing was performed at 24 hours after training over a 5-minute interval, which is sufficient to induce reconsolidation [75, 76]. The average freezing was 55 +/- 10%. Tissue was immersed in RNAlater (Qiagen) and immediately frozen. Animals that were handled but not trained were dissected at the same time of day (within 30 minutes) to control for variations due to circadian rhythms (CC30', CC4 and CC12). The protocol was repeated over the course of 2 weeks to obtain 9 animals (2 hippocampi) per group, so that 9 independent FC experiments were represented in each time point and all animals for each group were dissected at the same time of day. For subsequent qPCR and western blots the same protocol was followed, but tissue was collected only at FC30', CC30' and RT30'. Tissue for protein extraction was not immersed in RNAlater. All experiments were approved by the Institution of Animal Care and Use Committee of the University of Pennsylvania and were carried out in accordance with all National Institutes of Health guidelines.
RNA extraction was performed using Qiagen RNAeasy Microarray Tissue kit. All RNA extractions were performed the same day within a week of tissue collection. Induction of positive controls Arc, Fos and Dusp1 was confirmed by qPCR in the same samples previous to submission to microarrays analysis (data not shown). RNA was submitted to the University of Pennsylvania molecular profiling core for cRNA preparation, hybridization and scanning. Samples were simultaneously hybridized to an Affymetrix Mouse 1.1 Gene-EST 96 sample array plate. Target preparation and hybridization protocols were conducted as described in the Affymetrix GeneChip Expression Analysis Technical Manual. Biotinylated cRNA were prepared from 3 µg total RNA using the Ambion WT Expression Kit. Hybridization, staining and washing was performed using the GeneTitan Hybridization, Wash and Stain Kit for WT Array Plates. Fluorescent signal scanning was performed using the GeneTitan multi-channel instrument. The average signal from two sequential scans was calculated for each microarray feature. RMA normalization was performed using Affymetrix Power tools.
Principal component analysis (PCA) and differential expression analysis was performed using the R programming language and available packages from the R/Bioconductor project . PCA was performed after all expression values were standardized against the mean. The contribution of the first three principal components to the variance is as follows: PC1: 40%, PC2, 17% and PC3, 8%. The data was normalized against the first principal component to remove global systemic artifacts in all subsequent analyses of differential expression. Differential expression was performed against time of day matched controls. Local false discovery rates calculation using empirical null distributions was performed by applying the locfdr package to pair-wise t-statistics obtained after PC1 normalization. Estimating an empirical null hypothesis has been shown to add power to estimates of differential gene expression in large-scale studies . However this approach has never been used before to study genome-wide changes gene expression in the brain. Genome-wide fuzzy hierarchical clustering as originally detailed by Gasch and Eisen  was implemented using Cluster 3.0 using average linkage on average expression values per condition after RMA normalization. Clustering of individual replicates did not produce any clusters.
Cross-platform ID mapping, functional annotation and functional interaction analyses
Mapping of gene IDs across different platforms as well as enrichment of functional annotation was assessed using the Database for Visualization and Integrative Discovery (DAVID) . Functional annotation was limited to the following sources to increase information and limit redundancy: GO Biological process, GO Molecular Function, KEGG pathways, and SwissProt and Protein Information Resource keywords. Enrichment for each term was defined relative to the all mouse probe-sets present in the microarray, and was defined as a p < 0.05 with at least 3 genes per term per dataset. Fuzzy Heuristical clustering was performed using kappa similarity >0.3 and final group membership of at least 4 functional terms. Enriched functional clusters were defined as enrichment score ≥1.3 (p-value geometric mean between all genes within the cluster <0.05). Functional interaction analysis was performed with STRING , using co-expression, experimental, database and PubMed text-mining data. The cut-off interaction score was 0.4 (medium confidence) and the number of additional interactions equal to double the number of initial nodes. Clusters within the interaction network were obtained using K-means , K was evaluated from 2-6, results are displayed for the best fitting K value (k = 4).
Hippocampi were dissected 1 hour after fear conditioning and flash frozen in liquid nitrogen. Frozen nuclei were thawed on ice and homogenized in 1 mL NIB-250 (15 mM Tric-HCl pH 7.5, 60 mM KCl, 15 mM NaCl, 5 mM MgCl2, 1 mM CaCl2, 250 mM Sucrose, protease inhibitor cocktail (Sigma), 1 mM DTT, 10 mM sodium butyrate) + 0.3% NP-40 in a Type B dounce homogenizer. After 5 minutes on ice, nuclei were pelleted at 600 g for 5 minutes at 4°C and washed in 10:1 NIB-250 without detergent. Nuclei were pelleted and dissolved in 250 ul of 0.4N H2SO4. Acid extraction was performed on a nutator at 4°C for 3 hours, spun at 1500 g for 5 min, and the supernatant was set aside. The pellet was dissolved in another 250 ul of 0.4N H2SO4, rocked for 1 hour at 4°C, and spun at 1500 g for 5 minutes. The two supernatants were combined, 125 ul of 100% TCA was added, and precipitation was allowed to proceed overnight at -20°C. Samples were spun at 3400 g for 10 minutes, aspirated, and washed with 1 mL 0.1% HCl-acetone. Pellets were washed an additional 2× in acetone, allowed to dry, and redissolved in 20 ul H2O. Histone pellets were prepared for mass spectrometry experiments as previously described . In brief, histones were reacted with a 3:1 propionic anhydride/isopropanol mixture, digested with trypsin at a 20:1 protein/trypsin ratio and then reacted with the propionylation reagent one more time to cap the newly generated N-termini. Digested histones were then loaded onto and separated by reversed-phase high-performance liquid chromatography (HPLC) on an EasyLC 1000 system (Thermo, San Jose, CA) using a 75-μm-inner diameter fused silica column packed with 10-15 cm of 5-μm C18 (Michrom, Auburn, CA). The HPLC gradient was a 1-30% buffer B in buffer A (buffer A, 0.1 M acetic acid; buffer B, 95% acetonitrile in 0.1 M acetic acid) for 35 min followed by 30-99% buffer B for 30 min was used to elute peptides, which were ionized into an Orbitrap Velos instrument via electrospray ionization. Peptides were analyzed on The Orbitrap Velos mass spectrometer (ThermoFisher Scientific, San Jose, CA) with full scans of m/z = 290-1200 with a resolution of 30,000, followed by 7 MS/MS spectra collected in the ion trap. All data was manually analyzed and quantified as previously reported .
Ortholog mapping and molecular phylogenetic analysis
Orthologs of mice Hist2h2ab were determined using BLAST against the OrthoMCL database  (ortholog group OG5_126570). All mouse and human sequences within the ortholog group plus mouse and human H2A.Z were aligned using T-coffee . Phylogenetic reconstruction was performed using PHYML  with 100 bootstrap support and displayed using iTOL version 2 . Histone variant macroH2A was not included because the presence of the macro domain produces severe long-branch attraction in the phylogeny reconstruction.
Quantitative RT-PCR (qPCR)
RNA extraction was performed using Qiagen RNAeasy lipid tissue kit with modifications to obtain both small RNAs as well as mRNAs. Briefly, RNA precipitation was performed using 100% EtOH and washes were performed using Qiagen's RWT buffer. Concentration and purity was quantified by NanoDrop spectrophotometry (Thermo Fisher Scientific, Wilmington, DE). For quantitative real-time RT-PCR, reactions were prepared in 384-well optical reaction plates (ABI, Foster City, CA) with optical adhesive covers (ABI). Two technical replicates were used. Reactions were carried out in ViiA7 real-time PCR system (Invitrogen). For mRNA qPCR, generation of cDNA was carried out by the RETROscript kit (Ambion) with 1 µg of RNA as template. Taqman gene expression assays for all genes were obtained from ABI (Invitrogen). Data was normalized to Gapdh prior to calculation of differences.
Relative quantification of gene expression was performed according to ABI's User Bulletin #2. Fold change was calculated using the delta delta Ct method. The data presented is the calculated mean for the biological replicates with n = 8 (i.e., the number of mice examined). We used t-tests to compare fold change values for each gene in each comparison of interest. Two-tailed p-values are reported. Fos induction was used as positive control on all qPCR runs. For microRNA and snoRNAs, qPCR was performed using the miScript system from Qiagen. Reverse transcription was performed using miScript II RT kit using Hiflex buffer. MicroRNA and small nucleolar RNA detection by real-time PCR was performed using the miScript SYBR green PCR Kit. miSCRIPT primer assays were obtained from Qiagen, with the exception of primers for Snord14d and Snord14e. See Additional data file 12 for assays IDs and sequences. Data was normalized to Snord68 prior to calculation of differences. Fold changes were calculated as detailed for mRNAs above.
Western blot analysis
Frozen hippocampal tissue was homogenized in RIPA buffer with protease and phosphatase inhibitors. Proteins were separated by 4-20% Tris-Glycine SDS-PAGE and transferred to polyvinylidene difluoride (PVDF) membranes. Membranes were blocked in 5% BSA-TBST and incubated overnight at 4°C in primary antibody for CAMKIIγ (1:1000 abcam). They were washed, and incubated with appropriate horseradish peroxidase-conjugated goat anti-mouse or anti-rabbit IgG (1:5000, Santa Cruz) for 1 hour. Blots were exposed using ImageQuant LAS 4000 digital imaging system and quantified using ImageJ. Blots were stripped and re-probed with anti-β-tubulin antibody (1:20000, Sigma). Density of CAMKIIγ signal was normalized to β-tubulin levels.
Availability of Supporting Data
Microarray data generated in this study is publicly available through GEO (GSE50423). Mass Spec data is available at ChorusProject.org under project name hippocampi histone in mice or with the following links:
This research was supported by NRSA training grants T32NS007413 (to L.P.; M. Robinson, P.I.), T32HL007953 (to M.W.; A.I. Pack, P.I.) and K12GM081259 (to J.C., Paterson Y., P.I.) T.A. acknowledges funding from the Brush Family Professorship and R01 MH087463. B.A.G. acknowledges funding from NIH Innovator grant DP2OD007447 from the Office of the Director. We thank Stylianos Monos, Erik van Tilborg and Emily JoAnne Davis for technical support and Hannah Schoch, Morgan Bridi, Anthony Olalerin, John Hogenesch and Marcos Frank for helpful discussions.
Cost of article publications fees was covered through grant R01 MH087463 from the NIH/NIMH.
This article has been published as part of BMC Genomics Volume 16 Supplement 5, 2015: Proceedings of the 10th International Conference of the Brazilian Association for Bioinformatics and Computational Biology (X-Meeting 2014). The full contents of the supplement are available online at http://www.biomedcentral.com/bmcgenomics/supplements/16/S5.
- Scherer A: Batch effects and noise in microarray experiments: sources and solutions. 2009, John Wiley & SonsView ArticleGoogle Scholar
- Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003, 4: 249-264. 10.1093/biostatistics/4.2.249.View ArticlePubMedGoogle Scholar
- Agranoff BW, Davis RE, Casola L, Lim R: Actinomycin D blocks formation of memory of shock-avoidance in goldfish. Science. 1967, 158: 1600-1601. 10.1126/science.158.3808.1600.View ArticlePubMedGoogle Scholar
- Flood JF, Rosenzweig MR, Bennett EL, Orme AE: The influence of duration of protein synthesis inhibition on memory. Physiology & behavior. 1973, 10: 555-562. 10.1016/0031-9384(73)90221-7.View ArticleGoogle Scholar
- Zovkic IB, Guzman-Karlsson MC, Sweatt JD: Epigenetic regulation of memory formation and maintenance. Learn Mem. 2013, 20: 61-74. 10.1101/lm.026575.112.PubMed CentralView ArticlePubMedGoogle Scholar
- Peixoto L, Abel T: The role of histone acetylation in memory formation and cognitive impairments. Neuropsychopharmacology. 2013, 38: 62-76. 10.1038/npp.2012.86.PubMed CentralView ArticlePubMedGoogle Scholar
- Nelson ED, Monteggia LM: Epigenetics in the mature mammalian brain: effects on behavior and synaptic transmission. Neurobiol Learn Mem. 2011, 96: 53-60. 10.1016/j.nlm.2011.02.015.PubMed CentralView ArticlePubMedGoogle Scholar
- Alberini CM: Transcription factors in long-term memory and synaptic plasticity. Physiol Rev. 2009, 89: 121-145. 10.1152/physrev.00017.2008.View ArticlePubMedGoogle Scholar
- Bourtchouladze R, Abel T, Berman N, Gordon R, Lapidus K, Kandel ER: Different training procedures recruit either one or two critical periods for contextual memory consolidation, each of which requires protein synthesis and PKA. Learn Mem. 1998, 5: 365-374.PubMed CentralPubMedGoogle Scholar
- Igaz LM, Vianna MR, Medina JH, Izquierdo I: Two time periods of hippocampal mRNA synthesis are required for memory consolidation of fear-motivated learning. J Neurosci. 2002, 22: 6781-6789.PubMedGoogle Scholar
- Bekinschtein P, Katche C, Slipczuk L, Gonzalez C, Dorman G, Cammarota M, Izquierdo I, Medina JH: Persistence of Long-Term Memory Storage: New Insights into its Molecular Signatures in the Hippocampus and Related Structures. Neurotox Res.Google Scholar
- Vianna MR, Igaz LM, Coitinho AS, Medina JH, Izquierdo I: Memory extinction requires gene expression in rat hippocampus. Neurobiol Learn Mem. 2003, 79: 199-203. 10.1016/S1074-7427(03)00003-0.View ArticlePubMedGoogle Scholar
- Power AE, Berlau DJ, McGaugh JL, Steward O: Anisomycin infused into the hippocampus fails to block "reconsolidation" but impairs extinction: the role of re-exposure duration. Learn Mem. 2006, 13: 27-34. 10.1101/lm.91206.PubMed CentralView ArticlePubMedGoogle Scholar
- Mamiya N, Fukushima H, Suzuki A, Matsuyama Z, Homma S, Frankland PW, Kida S: Brain region-specific gene expression activation required for reconsolidation and extinction of contextual fear memory. J Neurosci. 2009, 29: 402-413. 10.1523/JNEUROSCI.4639-08.2009.View ArticlePubMedGoogle Scholar
- Motanis H, Maroun M: Differential involvement of protein synthesis and actin rearrangement in the reacquisition of contextual fear conditioning. Hippocampus. 2012, 22: 494-500. 10.1002/hipo.20915.View ArticlePubMedGoogle Scholar
- Hermey G, Mahlke C, Gutzmann JJ, Schreiber J, Bluthgen N, Kuhl D: Genome-wide profiling of the activity-dependent hippocampal transcriptome. PLoS One. 2013, 8: e76903-10.1371/journal.pone.0076903.PubMed CentralView ArticlePubMedGoogle Scholar
- Keeley MB, Wood MA, Isiegas C, Stein J, Hellman K, Hannenhalli S, Abel T: Differential transcriptional response to nonassociative and associative components of classical fear conditioning in the amygdala and hippocampus. Learning & memory. 2006, 13: 135-142. 10.1101/lm.86906.View ArticleGoogle Scholar
- Levenson JM, Choi S, Lee SY, Cao YA, Ahn HJ, Worley KC, Pizzi M, Liou HC, Sweatt JD: A bioinformatics analysis of memory consolidation reveals involvement of the transcription factor c-rel. The Journal of neuroscience : the official journal of the Society for Neuroscience. 2004, 24: 3933-3943. 10.1523/JNEUROSCI.5646-03.2004.View ArticleGoogle Scholar
- Barnes P, Kirtley A, Thomas KL: Quantitatively and qualitatively different cellular processes are engaged in CA1 during the consolidation and reconsolidation of contextual fear memory. Hippocampus. 2012, 22: 149-171. 10.1002/hipo.20879.View ArticlePubMedGoogle Scholar
- Griggs EM, Young EJ, Rumbaugh G, Miller CA: MicroRNA-182 regulates amygdala-dependent memory formation. J Neurosci. 2013, 33: 1734-1740. 10.1523/JNEUROSCI.2873-12.2013.PubMed CentralView ArticlePubMedGoogle Scholar
- Ahn HJ, Hernandez CM, Levenson JM, Lubin FD, Liou HC, Sweatt JD: c-Rel, an NF-kappaB family transcription factor, is required for hippocampal long-term synaptic plasticity and memory formation. Learn Mem. 2008, 15: 539-549. 10.1101/lm.866408.PubMed CentralView ArticlePubMedGoogle Scholar
- Leach PT, Poplawski SG, Kenney JW, Hoffman B, Liebermann DA, Abel T, Gould TJ: Gadd45b knockout mice exhibit selective deficits in hippocampus-dependent long-term memory. Learn Mem. 2012, 19: 319-324. 10.1101/lm.024984.111.PubMed CentralView ArticlePubMedGoogle Scholar
- Abdi H, Williams LJ: Principal component analysis. Wiley Interdisciplinary Reviews: Computational Statistics. 2010, 2: 433-459. 10.1002/wics.101.View ArticleGoogle Scholar
- Devan BD, Goad EH, Petri HL, Antoniadis EA, Hong NS, Ko CH, Leblanc L, Lebovic SS, Lo Q, Ralph MR, McDonald RJ: Circadian phase-shifted rats show normal acquisition but impaired long-term retention of place information in the water task. Neurobiology of learning and memory. 2001, 75: 51-62. 10.1006/nlme.1999.3957.View ArticlePubMedGoogle Scholar
- Tapp WN, Holloway FA: Phase shifting circadian rhythms produces retrograde amnesia. Science. 1981, 211: 1056-1058. 10.1126/science.7193351.View ArticlePubMedGoogle Scholar
- Stephan FK, Kovacevic NS: Multiple retention deficit in passive avoidance in rats is eliminated by suprachiasmatic lesions. Behavioral biology. 1978, 22: 456-462. 10.1016/S0091-6773(78)92565-8.View ArticlePubMedGoogle Scholar
- Hughes ME, DiTacchio L, Hayes KR, Vollmers C, Pulivarthy S, Baggs JE, Panda S, Hogenesch JB: Harmonics of circadian gene transcription in mammals. PLoS genetics. 2009, 5: e1000442-10.1371/journal.pgen.1000442.PubMed CentralView ArticlePubMedGoogle Scholar
- Su AI, Wiltshire T, Batalov S, Lapp H, Ching KA, Block D, Zhang J, Soden R, Hayakawa M, Kreiman G, et al: A gene atlas of the mouse and human protein-encoding transcriptomes. Proceedings of the National Academy of Sciences of the United States of America. 2004, 101: 6062-6067. 10.1073/pnas.0400782101.PubMed CentralView ArticlePubMedGoogle Scholar
- Steward O, Wallace CS, Lyford GL, Worley PF: Synaptic activation causes the mRNA for the IEG Arc to localize selectively near activated postsynaptic sites on dendrites. Neuron. 1998, 21: 741-751. 10.1016/S0896-6273(00)80591-7.View ArticlePubMedGoogle Scholar
- Zheng F, Zhou X, Moon C, Wang H: Regulation of brain-derived neurotrophic factor expression in neurons. International journal of physiology, pathophysiology and pharmacology. 2012, 4: 188-200.PubMed CentralPubMedGoogle Scholar
- Oike Y, Hata A, Mamiya T, Kaname T, Noda Y, Suzuki M, Yasue H, Nabeshima T, Araki K, Yamamura K: Truncated CBP protein leads to classical Rubinstein-Taybi syndrome phenotypes in mice: implications for a dominant-negative mechanism. Human molecular genetics. 1999, 8: 387-396. 10.1093/hmg/8.3.387.View ArticlePubMedGoogle Scholar
- Wood MA, Kaplan MP, Park A, Blanchard EJ, Oliveira AM, Lombardi TL, Abel T: Transgenic mice expressing a truncated form of CREB-binding protein (CBP) exhibit deficits in hippocampal synaptic plasticity and memory storage. Learning & memory. 2005, 12: 111-119. 10.1101/lm.86605.View ArticleGoogle Scholar
- Alarcon JM, Malleret G, Touzani K, Vronskaya S, Ishii S, Kandel ER, Barco A: Chromatin acetylation, memory, and LTP are impaired in CBP+/- mice: a model for the cognitive deficit in Rubinstein-Taybi syndrome and its amelioration. Neuron. 2004, 42: 947-959. 10.1016/j.neuron.2004.05.021.View ArticlePubMedGoogle Scholar
- Korzus E, Rosenfeld MG, Mayford M: CBP histone acetyltransferase activity is a critical component of memory consolidation. Neuron. 2004, 42: 961-972. 10.1016/j.neuron.2004.06.002.View ArticlePubMedGoogle Scholar
- Bourtchouladze R, Lidge R, Catapano R, Stanley J, Gossweiler S, Romashko D, Scott R, Tully T: A mouse model of Rubinstein-Taybi syndrome: defective long-term memory is ameliorated by inhibitors of phosphodiesterase 4. Proceedings of the National Academy of Sciences of the United States of America. 2003, 100: 10518-10522. 10.1073/pnas.1834280100.PubMed CentralView ArticlePubMedGoogle Scholar
- Oliveira AM, Estevez MA, Hawk JD, Grimes S, Brindle PK, Abel T: Subregion-specific p300 conditional knock-out mice exhibit long-term memory impairments. Learning & memory. 2011, 18: 161-169. 10.1101/lm.1939811.View ArticleGoogle Scholar
- Oliveira AM, Wood MA, McDonough CB, Abel T: Transgenic mice expressing an inhibitory truncated form of p300 exhibit long-term memory deficits. Learning & memory. 2007, 14: 564-572. 10.1101/lm.656907.View ArticleGoogle Scholar
- Efron B: Large-scale simultaneous hypothesis testing: The choice of a null hypothesis. Journal of the American Statistical Association. 2004, 99: 96-104. 10.1198/016214504000000089.View ArticleGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the False Discovery Rate - a Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society Series B-Methodological. 1995, 57: 289-300.Google Scholar
- Farioli-Vecchioli S, Saraulli D, Costanzi M, Leonardi L, Cina I, Micheli L, Nutini M, Longone P, Oh SP, Cestari V, Tirone F: Impaired terminal differentiation of hippocampal granule neurons and defective contextual memory in PC3/Tis21 knockout mice. PLoS One. 2009, 4: e8339-10.1371/journal.pone.0008339.PubMed CentralView ArticlePubMedGoogle Scholar
- Li H, Radford JC, Ragusa MJ, Shea KL, McKercher SR, Zaremba JD, Soussou W, Nie Z, Kang YJ, Nakanishi N, et al: Transcription factor MEF2C influences neural stem/progenitor cell differentiation and maturation in vivo. Proc Natl Acad Sci USA. 2008, 105: 9397-9402. 10.1073/pnas.0802876105.PubMed CentralView ArticlePubMedGoogle Scholar
- Cole CJ, Mercaldo V, Restivo L, Yiu AP, Sekeres MJ, Han JH, Vetere G, Pekar T, Ross PJ, Neve RL, et al: MEF2 negatively regulates learning-induced structural plasticity and memory formation. Nat Neurosci. 2012, 15: 1255-1264. 10.1038/nn.3189.View ArticlePubMedGoogle Scholar
- Berdeaux R, Goebel N, Banaszynski L, Takemori H, Wandless T, Shelton GD, Montminy M: SIK1 is a class II HDAC kinase that promotes survival of skeletal myocytes. Nat Med. 2007, 13: 597-603. 10.1038/nm1573.View ArticlePubMedGoogle Scholar
- Finsterwald C, Carrard A, Martin JL: Role of salt-inducible kinase 1 in the activation of MEF2-dependent transcription by BDNF. PLoS One. 2013, 8: e54545-10.1371/journal.pone.0054545.PubMed CentralView ArticlePubMedGoogle Scholar
- Li S, Zhang C, Takemori H, Zhou Y, Xiong ZQ: TORC1 regulates activity-dependent CREB-target gene transcription and dendritic growth of developing cortical neurons. J Neurosci. 2009, 29: 2334-2343. 10.1523/JNEUROSCI.2296-08.2009.View ArticlePubMedGoogle Scholar
- Katoh Y, Takemori H, Lin XZ, Tamura M, Muraoka M, Satoh T, Tsuchiya Y, Min L, Doi J, Miyauchi A, et al: Silencing the constitutive active transcription factor CREB by the LKB1-SIK signaling cascade. FEBS J. 2006, 273: 2730-2748. 10.1111/j.1742-4658.2006.05291.x.View ArticlePubMedGoogle Scholar
- Jilg A, Lesny S, Peruzki N, Schwegler H, Selbach O, Dehghani F, Stehle JH: Temporal dynamics of mouse hippocampal clock gene expression support memory processing. Hippocampus. 2010, 20: 377-388.PubMedGoogle Scholar
- Cheval H, Chagneau C, Levasseur G, Veyrac A, Faucon-Biguet N, Laroche S, Davis S: Distinctive features of Egr transcription factor regulation and DNA binding activity in CA1 of the hippocampus in synaptic plasticity and consolidation and reconsolidation of fear memory. Hippocampus. 2012, 22: 631-642. 10.1002/hipo.20926.View ArticlePubMedGoogle Scholar
- Lubin FD, Roth TL, Sweatt JD: Epigenetic regulation of BDNF gene transcription in the consolidation of fear memory. J Neurosci. 2008, 28: 10576-10586. 10.1523/JNEUROSCI.1786-08.2008.PubMed CentralView ArticlePubMedGoogle Scholar
- Sakamoto K, Karelina K, Obrietan K: CREB: a multifaceted regulator of neuronal plasticity and protection. Journal of neurochemistry. 2011, 116: 1-9. 10.1111/j.1471-4159.2010.07080.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Oikawa K, Odero GL, Platt E, Neuendorff M, Hatherell A, Bernstein MJ, Albensi BC: NF-kappaB p50 subunit knockout impairs late LTP and alters long term memory in the mouse hippocampus. BMC neuroscience. 2012, 13: 45-10.1186/1471-2202-13-45.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen P, Zhao J, Li G: Histone variants in development and diseases. J Genet Genomics. 2013, 40: 355-365. 10.1016/j.jgg.2013.05.001.View ArticlePubMedGoogle Scholar
- Pina B, Suau P: Changes in histones H2A and H3 variant composition in differentiating and mature rat brain cortical neurons. Dev Biol. 1987, 123: 51-58. 10.1016/0012-1606(87)90426-X.View ArticlePubMedGoogle Scholar
- Santoro SW, Dulac C: The activity-dependent histone variant H2BE modulates the life span of olfactory neurons. Elife. 2012, 1: e00070-PubMed CentralView ArticlePubMedGoogle Scholar
- Fogel BL, Wexler E, Wahnich A, Friedrich T, Vijayendran C, Gao F, Parikshak N, Konopka G, Geschwind DH: RBFOX1 regulates both splicing and transcriptional networks in human neuronal development. Hum Mol Genet. 2012, 21: 4171-4186. 10.1093/hmg/dds240.PubMed CentralView ArticlePubMedGoogle Scholar
- Gehman LT, Stoilov P, Maguire J, Damianov A, Lin CH, Shiue L, Ares M, Mody I, Black DL: The splicing regulator Rbfox1 (A2BP1) controls neuronal excitation in the mammalian brain. Nat Genet. 2011, 43: 706-711. 10.1038/ng.841.PubMed CentralView ArticlePubMedGoogle Scholar
- Mikhail FM, Lose EJ, Robin NH, Descartes MD, Rutledge KD, Rutledge SL, Korf BR, Carroll AJ: Clinically relevant single gene or intragenic deletions encompassing critical neurodevelopmental genes in patients with developmental delay, mental retardation, and/or autism spectrum disorders. Am J Med Genet A. 2011, 2386-2396. 155AGoogle Scholar
- Antunes-Martins A, Mizuno K, Irvine EE, Lepicard EM, Giese KP: Sex-dependent up-regulation of two splicing factors, Psf and Srp20, during hippocampal memory formation. Learn Mem. 2007, 14: 693-702. 10.1101/lm.640307.PubMed CentralView ArticlePubMedGoogle Scholar
- Wibrand K, Panja D, Tiron A, Ofte ML, Skaftnesmo KO, Lee CS, Pena JT, Tuschl T, Bramham CR: Differential regulation of mature and precursor microRNA expression by NMDA and metabotropic glutamate receptor activation during LTP in the adult dentate gyrus in vivo. Eur J Neurosci. 2010, 31: 636-645. 10.1111/j.1460-9568.2010.07112.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Tognini P, Pizzorusso T: MicroRNA212/132 family: molecular transducer of neuronal function and plasticity. Int J Biochem Cell Biol. 2012, 44: 6-10. 10.1016/j.biocel.2011.10.015.View ArticlePubMedGoogle Scholar
- Kocerha J, Faghihi MA, Lopez-Toledano MA, Huang J, Ramsey AJ, Caron MG, Sales N, Willoughby D, Elmen J, Hansen HF, et al: MicroRNA-219 modulates NMDA receptor-mediated neurobehavioral dysfunction. Proc Natl Acad Sci USA. 2009, 106: 3507-3512. 10.1073/pnas.0805854106.PubMed CentralView ArticlePubMedGoogle Scholar
- Tsien JZ, Huerta PT, Tonegawa S: The essential role of hippocampal CA1 NMDA receptor-dependent synaptic plasticity in spatial memory. Cell. 1996, 87: 1327-1338. 10.1016/S0092-8674(00)81827-9.View ArticlePubMedGoogle Scholar
- Nakazawa K, McHugh TJ, Wilson MA, Tonegawa S: NMDA receptors, place cells and hippocampal spatial memory. Nat Rev Neurosci. 2004, 5: 361-372. 10.1038/nrn1385.View ArticlePubMedGoogle Scholar
- Han Z, He H, Zhang F, Huang Z, Liu Z, Jiang H, Wu Q: Spatiotemporal expression pattern of Mirg, an imprinted non-coding gene, during mouse embryogenesis. J Mol Histol. 2012, 43: 1-8. 10.1007/s10735-011-9367-x.View ArticlePubMedGoogle Scholar
- Wheeler G, Ntounia-Fousara S, Granda B, Rathjen T, Dalmay T: Identification of new central nervous system specific mouse microRNAs. FEBS Lett. 2006, 580: 2195-2200. 10.1016/j.febslet.2006.03.019.View ArticlePubMedGoogle Scholar
- Hansen KF, Karelina K, Sakamoto K, Wayman GA, Impey S, Obrietan K: miRNA-132: a dynamic regulator of cognitive capacity. Brain Struct Funct. 2013, 218: 817-831. 10.1007/s00429-012-0431-4.PubMed CentralView ArticlePubMedGoogle Scholar
- Hansen KF, Sakamoto K, Wayman GA, Impey S, Obrietan K: Transgenic miR132 alters neuronal spine density and impairs novel object recognition memory. PLoS One. 2010, 5: e15497-10.1371/journal.pone.0015497.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen L, Zhang J, Feng Y, Li R, Sun X, Du W, Piao X, Wang H, Yang D, Sun Y, et al: MiR-410 regulates MET to influence the proliferation and invasion of glioma. Int J Biochem Cell Biol. 2012, 44: 1711-1717. 10.1016/j.biocel.2012.06.027.View ArticlePubMedGoogle Scholar
- Campbell DB, D'Oronzio R, Garbett K, Ebert PJ, Mirnics K, Levitt P, Persico AM: Disruption of cerebral cortex MET signaling in autism spectrum disorder. Ann Neurol. 2007, 62: 243-250. 10.1002/ana.21180.View ArticlePubMedGoogle Scholar
- Rogelj B, Hartmann CE, Yeo CH, Hunt SP, Giese KP: Contextual fear conditioning regulates the expression of brain-specific small nucleolar RNAs in hippocampus. Eur J Neurosci. 2003, 18: 3089-3096. 10.1111/j.1460-9568.2003.03026.x.View ArticlePubMedGoogle Scholar
- Da Silva WC, Cardoso G, Bonini JS, Benetti F, Izquierdo I: Memory reconsolidation and its maintenance depend on L-voltage-dependent calcium channels and CaMKII functions regulating protein turnover in the hippocampus. Proc Natl Acad Sci USA. 2013, 110: 6566-6570. 10.1073/pnas.1302356110.PubMed CentralView ArticlePubMedGoogle Scholar
- Fiorenza NG, Rosa J, Izquierdo I, Myskiw JC: Modulation of the extinction of two different fear-motivated tasks in three distinct brain areas. Behav Brain Res. 2012, 232: 210-216. 10.1016/j.bbr.2012.04.015.View ArticlePubMedGoogle Scholar
- Myers KM, Carlezon WA, Davis M: Glutamate receptors in extinction and extinction-based therapies for psychiatric illness. Neuropsychopharmacology. 2011, 36: 274-293. 10.1038/npp.2010.88.PubMed CentralView ArticlePubMedGoogle Scholar
- Sanhueza M, Lisman J: The CaMKII/NMDAR complex as a molecular memory. Molecular Brain. 2013, 6: 10-10.1186/1756-6606-6-10.PubMed CentralView ArticlePubMedGoogle Scholar
- Suzuki A, Josselyn SA, Frankland PW, Masushige S, Silva AJ, Kida S: Memory reconsolidation and extinction have distinct temporal and biochemical signatures. J Neurosci. 2004, 24: 4787-4795. 10.1523/JNEUROSCI.5491-03.2004.View ArticlePubMedGoogle Scholar
- von Hertzen LS, Giese KP: Memory reconsolidation engages only a subset of immediate-early genes induced during consolidation. J Neurosci. 2005, 25: 1935-1942. 10.1523/JNEUROSCI.4707-04.2005.View ArticlePubMedGoogle 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: R80-10.1186/gb-2004-5-10-r80.PubMed CentralView ArticlePubMedGoogle Scholar
- Gasch AP, Eisen MB: Exploring the conditional coregulation of yeast gene expression through fuzzy k-means clustering. Genome Biology. 2002, 3: RESEARCH0059-PubMed CentralView ArticlePubMedGoogle Scholar
- Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4: P3-10.1186/gb-2003-4-5-p3.View ArticlePubMedGoogle Scholar
- Franceschini A, Szklarczyk D, Frankild S, Kuhn M, Simonovic M, Roth A, Lin J, Minguez P, Bork P, von Mering C, Jensen LJ: STRING v9.1: protein-protein interaction networks, with increased coverage and integration. Nucleic Acids Research. 2013, 41: D808-815. 10.1093/nar/gks1094.PubMed CentralView ArticlePubMedGoogle Scholar
- Hartigan JA, Wong MA: Algorithm AS 136: A K-Means Clustering Algorithm. Journal of the Royal Statistical Society Series C (Applied Statistics). 1979, 28: 100-108.Google Scholar
- Plazas-Mayorca MD, Zee BM, Young NL, Fingerman IM, LeRoy G, Briggs SD, Garcia BA: One-pot shotgun quantitative mass spectrometry characterization of histones. J Proteome Res. 2009, 8: 5367-5374. 10.1021/pr900777e.PubMed CentralView ArticlePubMedGoogle Scholar
- Fischer S, Brunk BP, Chen F, Gao X, Harb OS, Iodice JB, Shanmugam D, Roos DS, Stoeckert CJ: Using OrthoMCL to assign proteins to OrthoMCL-DB groups or to cluster proteomes into new ortholog groups. Curr Protoc Bioinformatics. 2011, Chapter 6 (Unit 6 12): 11-19.Google Scholar
- Taly JF, Magis C, Bussotti G, Chang JM, Di Tommaso P, Erb I, Espinosa-Carrasco J, Kemena C, Notredame C: Using the T-Coffee package to build multiple sequence alignments of protein, RNA, DNA sequences and 3D structures. Nat Protoc. 2011, 6: 1669-1682. 10.1038/nprot.2011.393.View ArticlePubMedGoogle Scholar
- Criscuolo A: morePhyML: improving the phylogenetic tree space exploration with PhyML 3. Mol Phylogenet Evol. 2011, 61: 944-948. 10.1016/j.ympev.2011.08.029.View ArticlePubMedGoogle Scholar
- Letunic I, Bork P: Interactive Tree Of Life v2: online annotation and display of phylogenetic trees made easy. Nucleic Acids Res. 2011, 39: W475-478. 10.1093/nar/gkr201.PubMed CentralView ArticlePubMedGoogle Scholar
- Jiao XL, Sherman BT, Huang DW, Stephens R, Baseler MW, Lane HC, Lempicki RA: DAVID-WS: a stateful web service to facilitate gene/protein list analysis. Bioinformatics. 2012, 28: 1805-1806. 10.1093/bioinformatics/bts251.PubMed CentralView ArticlePubMedGoogle Scholar
- Di Tommaso P, Moretti S, Xenarios I, Orobitg M, Montanyola A, Chang JM, Taly JF, Notredame C: T-Coffee: a web server for the multiple sequence alignment of protein and RNA sequences using structural information and homology extension. Nucleic Acids Res. 2011, 39: W13-17. 10.1093/nar/gkr245.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/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.