- Research article
- Open Access
Gene expression profiling in the stress control brain region hypothalamic paraventricular nucleus reveals a novel gene network including Amyloid beta Precursor Protein
BMC Genomics volume 11, Article number: 546 (2010)
The pivotal role of stress in the precipitation of psychiatric diseases such as depression is generally accepted. This study aims at the identification of genes that are directly or indirectly responding to stress. Inbred mouse strains that had been evidenced to differ in their stress response as well as in their response to antidepressant treatment were chosen for RNA profiling after stress exposure. Gene expression and regulation was determined by microarray analyses and further evaluated by bioinformatics tools including pathway and cluster analyses.
Forced swimming as acute stressor was applied to C57BL/6J and DBA/2J mice and resulted in sets of regulated genes in the paraventricular nucleus of the hypothalamus (PVN), 4 h or 8 h after stress. Although the expression changes between the mouse strains were quite different, they unfolded in phases over time in both strains. Our search for connections between the regulated genes resulted in potential novel signalling pathways in stress. In particular, Guanine nucleotide binding protein, alpha inhibiting 2 (GNAi2) and Amyloid β (A4) precursor protein (APP) were detected as stress-regulated genes, and together with other genes, seem to be integrated into stress-responsive pathways and gene networks in the PVN.
This search for stress-regulated genes in the PVN revealed its impact on interesting genes (GNAi2 and APP) and a novel gene network. In particular the expression of APP in the PVN that is governing stress hormone balance, is of great interest. The reported neuroprotective role of this molecule in the CNS supports the idea that a short acute stress can elicit positive adaptational effects in the brain.
Stressful life events are among the most potent factors that can trigger the development of psychiatric disorders such as depression and anxiety disorders [1–3]. Aberrations in the function of the hypothalamus-pituitary-adrenal (HPA) axis, the key control system of the body to balance stress hormones and the response to stress, already exist prior to the onset of clinical symptoms . The functionality of the HPA axis is mainly governed by genetic endowment, but developmental influences and life events, in particular stress experience early in life, can re-program the settings of the HPA axis [5–8].
The hypothalamus, as part of the HPA axis, is the centre of stress response and a region of the brain that integrates different stress signalling neuronal pathways. The hypothalamic paraventricular nucleus (PVN) is the main area of the hypothalamus where the corticotropin-releasing hormone (CRH, also known as corticotropin-releasing factor, CRF), the crucial neuropeptide that activates the secretion of corticotropin (ACTH), is produced and released. This effect, in turn, causes the secretion of glucocorticoids from the adrenal glands [9, 10]. The levels of ACTH and glucocorticoids in the plasma can be used as markers to monitor stress levels . In addition to CRH, other hormonal molecules such as arginine vasopressin (AVP) and oxytocin contribute to the regulation of the HPA axis' activity [12–16].
In major depression a hypothalamic hyper-drive is observed. This is constituted by the elevation of CRH, AVP and oxytocin, which may influence the clinical symptoms. In the PVN of depressed patients the total number of CRH expressing neurons showing co-localisation with AVP and the amount of CRH-mRNA are increased [17–19]. An important hallmark of HPA axis regulation is the negative feedback exerted by glucocorticoids on the production and release of CRH and AVP in the PVN, as well as of ACTH in the pituitary. Substantial evidence has been provided that this attenuation activity, which is an integral part of the stress response, is impaired in patients suffering from major depression .
Recently, we have determined stress-regulated genes in the hippocampus, a higher limbic centre in the brain, of two inbred mouse strains with differential psycho-phenotype, namely C57BL/6J and the DBA/2J by employing microarray analysis . The choice of these mouse strains is based on reports that these strains have differential response to stress [22, 23], differential basal anxiety  and differ in both, their cognitive abilities  and sensitivity to antidepressants [26, 27]. The microarray method, that belongs to the chip-based whole genome technologies, allows for unbiased approaches with the potential to identify new candidate genes and gene networks . Our results had led to the successful identification of stress-regulated genes and suggestion of possible signal transduction pathways involved.
As it is of great interest to study the impact of stress in brain regions directly involved in stress regulation, in this parallel study we have focused on the impact of stress experience on hypothalamic PVN governing the stress response. The PVN area was micropunctured from the brains of C57BL/6J and DBA/2J mice that had been stressed once by forced swimming  and mRNA profiles were determined by microarray analysis. Forced swimming is also being used routinely as test to monitor depression-like behaviour and drug effects . We report that Guanine nucleotide binding protein, alpha inhibiting 2 and Amyloid β (A4) precursor protein (APP) are up-regulated after stress and we suggest a novel gene network involved in stress response in the two mouse strains. This network implies that the expression of APP might be a neuroprotective component of stress adaptation in the PVN.
Basal gene expression differences between C57BL/6J and DBA/2J mice in PVN
To compare expression profiles between the two mouse strains C57BL/6J and DBA/2J and evaluate expression changes at different time points after stress exposure, we used cDNA microarrays. We first evaluated the differences of the two mouse strains in their basal expression profile in the PVN (see Fig. 1A for scheme of all comparisons). The microarray analysis revealed 670 genes with more than 1.4 fold difference in expression.
Among the most pronounced expression differences we identified genes with a variety of function, such as protein kinase activity, genes with extracellular ligand-gated ion channel activity, genes involved in protein homeostasis (e.g. ubiquitin-dependent protein catabolism, protein biosynthesis-, ribosome-, ribonucleoprotein complex-related genes), cell surface (defence response, peptide antigen binding, antigen presentation, endogenous peptide antigen), and chromosomal processes (chromosome segregation, response to DNA damage stimuli, mitosis, DNA repair, DNA metabolism, homologous chromosome segregation, cell cycle, cell division); examples are provided in table 1 (for the whole list see Additional File 1).
ACTH and gene expression profiles after forced swimming
To determine transcriptional profiles after stress, mice had been subjected to forced swimming for 5 min at 8:00 h or 12:00 h and then decapitated at 16:00 h in order to isolate the brains for our microarray studies as previously described  (experimental plan depicted in Fig. 1B). In this previous study the impact of forced swim stress on stimulation of the HPA axis was evaluated by the determination of corticosterone (CORT) levels in the serum. In order to monitor the activation of the pituitary that responds to the CRF release, we also measured the ACTH levels. A set of animals of both strains was decapitated directly after stress exposure at 08:00 h and at 12:00 h (non stressed control animals were also included). The results show a clear increase for ACTH after stress in both mouse strains and at both time points selected in our study (Fig. 2), like previously also for CORT . Thus, the pituitary as well as the adrenals responded immediately to the stressor. When animals were decapitated 4 h or 8 h after the stress, the levels of ACTH were not different from the non stressed controls, indicative of a functional negative feedback mechanism (Fig. 2). The levels in non-stressed and stressed animals were higher at 16:00 than in the morning though, most likely due to the circadian rhythm of the mice.
To reveal gene expression differences, samples from non stressed mice were compared to samples from mice that had been stressed 4 h or 8 h before decapitation. In C57BL/6J mice, 123 genes were ≥1.4 fold regulated 4 h after stress and 88 genes 8 h after stress. In DBA/2J mice, 185 genes and 96 genes were regulated at the respective time points. (Additional File 2, 3, 4 and 5).
Examples of the most interesting regulated genes are highlighted in tables 2, 3, 4, 5: In C57BL/6J mice, among the genes regulated 4 h after stress we distinguish phosphodiesterase 1C (PDE1C), Mitogen-activated protein kinase kinase kinase kinase 3 (MAP4K3) and polymerase (DNA directed) delta 1, catalytic subunit (POLD1) (see Table 2). In the group of genes regulated 8 h after stress we noticed heat shock protein 1, alpha (Hspca), Guanine nucleotide binding protein, alpha o (GNAO1) and the regulator of G-protein signaling 2 (RGS2) (see Table 3).
In DBA/2J mice, among the genes regulated 4 h after stress, mitogen activated protein kinase 3 (MAPK3), Guanine nucleotide binding protein, alpha inhibiting 2 (GNAi2), nuclear factor of activated T-cells, cytoplasmic, calcineurin-dependent 1 (NFATC1) were most striking (Table 4), while among the genes regulated 8 h after stress Amyloid beta (A4) precursor protein (APP), cyclin-dependent kinase inhibitor 1B (P27) (CDKN1B), Transcription factor 7-like 2, T-cell specific, HMG-box (TCF7L2) were most prominent (Table 5).
The reaction to stress on the transcriptome level differs between mouse strains and displays phases
To gain more insight into the potential functions of the regulated genes, each set of regulated genes was sorted according to their ontology. Genes coding for mitochondrial, biosynthetic and metabolic molecules, receptors, signal transduction molecules, transcription and mRNA processing molecules, ion channels and ion transport molecules, vesicular transport molecules and cytoskeleton components are highly represented (Table 6). While the proportion of receptors and signal transduction molecules decreased between 4 h and 8 h after stress, the proportion of mitochondrial, biosynthetic and metabolic proteins increased in both mouse strains.
These observations suggested that the main focus of expression regulation after stress is shifted over time. To further address the question of the durability of the gene expression changes, we compared the genes regulated at 4 h with those regulated at 8 h after stress for each mouse strain. Surprisingly, the set of regulated genes completely changed from 4 h to 8 h after stress in both mouse strains (Table 7). Genes that responded to the stressor after 4 h showed normalized expression levels (i.e. less than 1.4 fold regulation) after 8 h, while other genes showed up as regulated at that time (Table 7). We also observed that the stress-induced changes of the PVN transcriptome were entirely different between C57BL/6J and DBA/2J mice. Interestingly though, a convergence was apparent 8 h after stress (Table 7).
Possible signalling pathways elicited after forced swimming in the PVN
The observation that the percentage of receptors and signalling molecules among the regulated genes decreased from 4 h to 8 h after stress (Table 6), together with the phased reaction of the transcriptome to stress (Table 7), led us to hypothesise that genes regulated at 4 h have pathway connections to genes regulated at 8 h. Employing a pathway building program to test this hypothesis, we identified links between genes responding at 4 h and genes responding at 8 h. For example, GNAi2, found to be up-regulated 4 h after stress in DBA/2J mice, is upstream of APP, which is up-regulated 8 h after stress. Additionally, NFATC1, found to be up-regulated 4 h after stress, is upstream of heat shock protein 70, HSPAIA (HSP70), which increases 8 h after stress (Fig. 3A).
In C57BL/6J mice, p21[CDKN1A]-activated kinase 2 (PAK2), down-regulated 4 h-after stress, inhibits the expression of 'expressed in non-metastatic cells 1 protein' (NME1) via regulation of tumor necrosis factor alpha (TNFα, FIG. 3B). NME1, in turn, is upstream of sprouty homolog 2 (Drosophila; SPRY2), which is regulated 8 h after stress. Another connection between genes that are both regulated 8 h after stress is between GNAO1 and RGS2 (Fig 3B).
Validation of GNAi2 and APP expression and regulation in the PVN
To validate and further analyse the expression changes of the genes GNAi2 and APP that are linked by a pathway in the PVN of DBA/2J mice, we used real-time PCR. RNA samples from the original punctures were amplified (one round of amplification) and subjected to RT-PCR without pooling. The results confirmed the up-regulation of GNAi2 4 h after stress (2.2 fold, Fig. 4A) detected by the microarray (2.1 fold). To test whether this regulation is specific for DBA/2J mice, or may also occur in C57BL/6J mice (but was not detected in our microarray screening, because of the selection filtering), we also tested the respective samples from C57BL/6J mice. The results showed a non-significant increase in this mouse strain, implying that the regulation is rather strain-specific (Fig. 4A).
Similarly to GNAi2, we validated the expression and regulation of APP 8 h after forced swimming by real-time PCR (Fig. 4B), which was found in the microarray analysis.
To visualize the regulation of this expression with spatial resolution in the PVN, in situ hybridization was performed on coronal brain-sections, followed by semi-quantification of the mRNA signal. The results showed a strong signal in the PVN area (Fig. 4C) and the analysis confirmed again the up-regulation to a level of 1.6 fold (data not shown). In situ hybridization showed also specific expression of APP in the PVN (Fig. 4C).
Finally, to test whether stress-regulation of GNAi2 and APP is specific to the PVN, we also analysed the hypothalamic region just anterior and posterior to the PVN. The results show no significant change in the expression levels of GNAi2 and APP for the respective time point (Fig. 4D).
Clustering analysis of GNAi2 - APP
Having corroborated a potential role of the GNAi2 - APP connection in stress response of DBA mice, we used clustering analysis to identify genes that display similar expression changes throughout the conditions we analysed. In setting up the clustering analysis, we considered both up- and down-regulation, because some transcriptional regulators, such as GR, are able to both up- and down-regulate genes. Genes identified in a clustering analysis may be under a common transcriptional control, or influencing each other.
The dendrograms revealed 10 and 15 genes for GNAi2 (two transcripts of this gene were spotted) and 12 for APP (Fig. 5, see also Additional File 6 for the expression profiles). Several genes with known function were among these nearest neighbours.
To test potential, already described connections between these genes, we again used a pathway building program. Interestingly, we found that GNAi2 is located in a signalling cascade, where 'Heart and neural crest derivatives expressed transcript 2' (HAND2) is upstream, APP is downstream and MAP3K2 further downstream. DHDDS and GNAi2 share a common upstream regulator (AGT) as well as a common downstream target (PARK2). GNAi2 and COPS5 have a common target as well. In addition, GNAi2 and APP are found both downstream of another common regulator (RBL2, Fig. 6A). Furthermore, for APP we identified common upstream regulators with PAPOLA (CDC2 and TNF), EN1 as its regulator and HSD17B4 as its target. Interestingly, DGKZ and APP are linked in a feedback loop (Fig. 6B). Thus, the clustering analysis revealed functionally related genes, indeed.
The PVN of the hypothalamus is pivotal in governing physiological stress response. We examined the impact of forced swimming as an acute stressor on gene expression in the PVN of C57BL/6 and DBA/2J mice. These inbred mouse strains have been used as a genetic animal model of depression-like behaviour and are characterized by a different stress-responsiveness, since DBA/2J mice display a stronger behavioural response to stressful conditions [22, 23, 29]. We discovered that the stress-regulated genes code mainly for receptors and signal transduction molecules, as well as numerous biosynthetic molecules. This result is consistent with a previous study of gene expression profiling in the hypothalamus of mice stressed by immobilization, where genes involved in energy and lipid metabolism, apoptosis, signal transduction, DNA repair, protein biosynthesis, and structure integrity-related genes were found. However, this study referred to the entire hypothalamus and the number of genes studied was comparatively small (6016 genes) , while in our study approximately 19900 gene transcripts have been studied. In addition, one should keep in mind that categorically distinct acute stressors elicit distinct transcriptional profiles in the PVN . Nevertheless, a few genes were found in both studies, for example ATP synthase, H+ transporting mitochondrial F1 complex and ribosomal protein L30 .
In principle, all genes that we found differentially expressed or differentially regulated between the two mouse strains are candidates for the explanation of the differential response to some stressors, reflecting the previously proposed differences in corticosteroid signalling . In addition, we checked genes with known functions that could contribute to the strain differences such as POMC1 (Pro-opiomelanocortin-alpha), GR (Glucocorticoid receptor), CRHR1 and CRHR2 (Corticotrophin Releasing Hormone Receptor 1 and 2). Only GR showed differential regulation, namely a down-regulation in DBA/2J mice 4 h after stress. Causality is difficult to prove in feedback systems such as the HPA axis, i.e. when considering the complex connections of compensatory mechanisms that emerged in this evolutionarily old threat response system . Nevertheless, this difference could be indicative of differential stress signalling, as GR is the most important mediator of corticosteroid action.
Since DBA/2J mice also exhibit a higher sensitivity to antidepressants than C57BL/6  we also specifically investigated genes that had been associated with antidepressant response before, such as the immunophilin FKBP5 , which is an efficient regulator of GR , the multidrug resistance protein ABCB1A/B that determines brain tissue penetration of many antidepressant drugs , the serotonin receptor 5 hT2A, and the transporter proteins SLC6A4 and SLC6A15 (soluble carrier family 6 member 4 and 15) . Among those, SLC6A15 was down-regulated 4 h after stress in DBA/2J mice, but not in C57BL/6 mice. Since stress and GR action is intertwined with the action of antidepressants [3, 39, 40], also any of the stress-induced genes could contribute to the action of antidepressants. Like with other screening studies, it is certainly premature to delineate direct candidates for novel antidepressant targets or for diagnostic markers from this study. However, the synopsis of our results together with results from different screening efforts in genetics , proteomics, metabolomics etc., will yield convergence and thus allow selection of the most promising candidates.
Our study provides strong evidence for a time-phased response of the PVN transcriptome to the stressor. We have previously described a phased stress response for the CA3 hippocampal region as well . This suggests that this might be not an area-specific phenomenon, but rather a more general mechanism. Interestingly, a number of genes are regulated by stress in the hippocampal CA3 as well as in the PVN, e.g. DPYSL2, SNAP25, GNAO1 . Like for the stress-regulated genes in CA3, we used also for the regulated genes in the PVN a pathway building program to propose novel signal transduction pathways elicited after stress. This pathway analysis revealed an interesting link between GNAi2 and APP whose stress regulation we validated by real time qPCR.
GNAi2, that displayed increased expression after stress, is a plasma membrane protein and a member of the G(i) proteins that inhibit adenylate cyclase. Many important hormones and neurotransmitters, including acetylcholine, dopamine and serotonin, use the G(i) pathway to evoke physiological responses . In addition to inhibiting adenylate cyclase, G(i) regulates c-Src (and thereby STAT3) and Rap1 pathways, for example. This well studied inhibition may be physiologically relevant, in particular in inhibiting the effects of cAMP to modulate secretion in response to beta-adrenergic stimuli [42, 43]. Despite their established role in modulating cellular processes, very little is known regarding molecular mechanisms underlying the transcriptional regulation of any G-protein. It has been reported previously, that the increase of reactive oxygen species in K562 cells up-regulates Galpha(i2) . Here we report an increase of GNAi2 in DBA/2J mice 4 h after stress, which reveals a more direct impact of stress on GNAi2 expression by psychological stressors not reported so far.
Similarly to GNAi2, we were able to identify an increase of APP mRNA in the PVN of DBA/2J mice, 8 h after exposure to forced swimming. APP is mostly known for being the source of the toxic amyloid-β (Aβ) peptide found in neuritic plaques of Alzheimer's disease (AD) patients. However, it may in addition be a functionally important molecule in its full-length configuration, as well as be the source of numerous fragments with varying effects on neural function . Besides its cellular function (axonal transport, cell adhesion, cholesterol metabolism and gene transcription) throughout the body, it is also reported to participate in a number of important functions in the CNS, i.e. neuronal development, survival and plasticity  and to be a central molecule in many metabolic and regulatory pathways, so that its regulation may impact on a network of genes . Moreover, there is a growing body of evidence that APP may be part of the cellular response to stress, and serve neurotrophic and neuroprotective functions [47, 48], although the neuroprotective role of APP remains controversial.
Our observed up-regulation of APP is in line with the rapid increase in amyloid precursor protein immunoreactivity in the supraoptic and paraventricular nuclei described in the rat hypothalamus after osmotic stress . Furthermore, exposure of wild-type animals to an enriched environment can up-regulate APP-expression , although it remains unclear whether full-length APP or one of its fragments are involved Since APP has been reported to exert beneficial effects for neurons [48, 51–53], the increase after stress might be important as a protective mechanism for the sensitive PVN area.
Linked to the action of APP, we also observed a regulation of the α-secretase ADAM10 8 h after stress. Upon neuronal over-expression in a mouse model for Alzheimer's disease, this protease has been shown to increase the secretion of sAPPα, reduce the formation of Aβ peptides and prevent their deposition in plaques as well as to alleviate impaired long-term potentiation and cognitive deficits . In the stress exposed mice investigated here, there was pronounced strain specificity, i.e. 8 h after stress ADAM10 was strongly down-regulated in C57BL/6J mice (Table 3), but up-regulated in DBA/2J mice (Table 5). Of note, the inhibition of adenyl cyclase, from GNAi2 could also lead to the non amyloidogenic a-secretase pathway, resulting elevated sAPPα by likely shifting to the protein kinase competing signalling pathway. Moreover, this increase of GNAi2 after stress appears also to be strain-specific, because it was found in DBA/J2, but not in C57BL/6J mice [55, 56]. The hypothesis emerging from these observations, i.e. a role of sAPPα in differentially shaping stress response needs to be tested by further experiments, for example by more directly manipulating the GNAi2 signalling pathway, ADAM10 expression and the activity and APP or its metabolite sAPPα, in vitro and possibly also in vivo. These studies could include stress-exposure and antidepressant response in APP transgenic and/or APP knockout mice.
The clustering analysis we have performed and the pathway analysis that followed turned out to be very useful tools and revealed new possible signalling pathways involving GNAi2 and APP. Obviously, mechanisms in addition to gene expression, such as protein phosphorylation, protein-protein binding etc. also operate in regulatory networks. Nevertheless, we believe that the identified network adds significantly towards the understanding of the complicated mechanisms of the response of the PVN to stress. Moreover, we propose that the combination of our results with future results expected from research efforts targeted towards proteins, gene polymorphisms, epigenomes, metabolome etc. will help identifying markers for diagnosis, stratification of subjects, and possibly also novel drug targets.
Given the importance of PVN hypothalamic area for the physiological stress response and the discussed neuroprotective role of APP, the up-regulation of GNAi2 and APP mRNA levels after a mild stress in mice is suggested as it could be an adaptational stress-response in stress-responsive mice. Novel molecular pathways involving stress regulated genes that respond to stress in the PVN area, have been revealed based on clustering and signalling cascade pathway analysis.
Our study was performed with C57BL/6J and DBA/2J male mice with an age of 3-5 months (supplied from Harlan Winkelman, Borchen) and single housed under a 06:00-18:00 light cycle and subjected to forced swimming to induce an emotional stress, as previously described (Tsolakidou et al 2008). Briefly, the animals were placed in a 24 cm high 11 cm diameter cylinder filled with water at 22-25°C for 5 min. Afterwards they were transferred to their own original cage (i.e. single housed again). The animals elicited a mixture of behavioural activity that can be described as climbing, swimming and immobility, with the latter reflecting a passive stress coping behaviour. The experimental procedures were conducted in accordance with the Guide for the Care and Use of Laboratory Animals of the European Union (European Communities Council Directive 86/609/EEC) and approved by the Government of Bavaria, Germany.
The mice were separated into three groups per strain as previously described (Tsolakidou et al 2008) (see Fig. 1B). The first and the second group were subjected to stress at 8:00 or at 12:00, respectively. The third group was not stressed and further on used as reference group for the microarray experiments. All animals were decapitated at 16:00 to avoid possible interference by circadian variations of corticosterone (CORT) levels. Thus, the first and the second group have been actually sacrificed 8 h or 4 h after stress, respectively. Trunk blood was collected for determination of ACTH concentrations (10 animals/condition) and dissected brains from the same animals were frozen on dry ice and stored at -80°C (to be used for expression profiling). To monitor hormone levels acutely after stress, a small set of animals (5/condition) were sacrificed acutely at the respective time points (8:00 and 12:00). Plasma ACTH concentrations were determined in a radioimmune assay (ICN Biomedicals).
Micropuncture and RNA preparation
Micropuncturing of the PVN and adjacent region on coronal tissue sections (200 μm thick) was applied under dry ice-cold conditions. To control for the accuracy of the puncture the sections were stained (cresyl violet) afterwards. Total RNA was extracted from the collected tissue (Trizol protocol, GibcoBRL). Samples from six animals were pooled to minimize the impact of biological variance, which is intrinsic to all organisms and can be substantial even in inbred mice . After two rounds of amplification (Amino Allyl MessageAmp aRNA kit, Ambion) the RNA was labelled with Cy3 or Cy5 dyes (Amersham)
Microarray hybridization and Analysis
Spotted cDNA microarray chips (MPIP 24K microarray platform with sequences from 19900 different genes,  were used, which allow a greater degree of flexibility in the choice of arrayed elements. Internal controls (sequences from other organisms) were used as contribution to ensure the quality of the data . The microarray experiments were performed by competitive hybridisation of two differentially labelled (Cy3 or Cy5) probes of amplified total RNA samples. 10 arrays were used for each comparison, that is 5 technical replicates and a dye-swap with another 5 technical replicates. 20 μg of each Cy3- or Cy5-labeled sample were denatured at 95°C for 3 min in hybridisation buffer (50% formamide, 50 mM sodium phosphate-buffer pH 7.0, 5× Denhard's solution [Sigma, Taufkirchen, Germany], 6× SSC, 0.5% SDS, 0.4 mg/ml murine COT1-DNA [Invitrogen] and 5 μg poly(dA) [Amersham]). The hybridisation was performed in chambers submerged in a water bath at 42°C for 16 h. The arrays were washed for 15 min with 2× SSC/0.2% SDS at 60°C, in 0.5× SSC for 15 min at 60°C, rinsed in 0.2× SSC for 1 min at room temperature, shaken vigorously in 0.05× SSC at room temperature and finally air dried. All slides were scanned immediately afterwards.
Scanning was performed using a ScanArray 4000 laser scanner and ScanArray 3.1 Software (Perkin Elmer, Boston, USA) with a fixed PMT gain of 80%, and 98% (Cy3) or 70% (Cy5) laser power. The QuantArray software 18.104.22.168 (Perkin Elmer, Germany) and the fixed circle-analysis method were used to perform the quantification. Data were imported into a PostgreSQL relational database for further analysis. Raw data were normalized according to the procedure outlined elsewhere  and subjected to a two-sided one sample t-test for significantly differential expression. Candidate genes were screened for using thresholds of |fold regulation| 1.414 and |Z-score| 1.423. In order to exclude any possibility of wrong annotation the majority of the differentially expressed genes, in particular those of major interest were sequenced. All data are deposited at the GEO server (GSE20877).
Ontology Groups Classification and Molecular Pathways Building
To classify and characterize the differentially expressed genes that resulted from the microarray analysis the Gene Ontology data annotation (version 161) (http://www.geneontology.org, http://www.ebi.ac.uk/GOA/) was used. In addition, the text mining program Pathway Studio 5.0 (Ariadne Genomics, Rockville, USA) was used for the identification of pathways linking the regulated genes. Downstream targets of the genes regulated 4 h after stress were looked for in the group of genes regulated 8 h after stress and upstream regulators of the latter in the former group. The main criteria applied for the final selection of the pathways were: 1. the last step of the pathway has to indicate molecular synthesis/expression so as to validate the expression change described on the microarray result, 2. confirmation of the annotation of the selected genes by sequencing and 3. validation of the literature references used by the program.
In situ Hybridisation
18 μm thin cryostat sections prepared from frozen stored brains were thaw-mounted on poly L-Lysine-coated slides, dried and kept at -80°C. [α-35S] UTP-labelled riboprobes for amyloid β (A4) precursor protein (APP) and guanine nucleotide binding protein, alpha inhibiting 2 (GNAi2) genes were prepared by in vitro transcription from corresponding cDNA clones (sense probes were also included for control). The sections were dried, fixed in 4% paraformaldehyde, washed in PBS (3 times) and subjected to acetylation using 0.25% acetic anhydride in 0.1 M triethanolamine-HCL, pH 8.0. After subsequent dehydration in increasing concentrations of ethanol, brain sections were saturated with 100 μl of hybridisation buffer containing ~ 2.5 × 106 cpm 35S-labelled riboprobe. Brain sections were incubated overnight at 62°C or 58°C. Then, the sections were rinsed in 4× SSC (standard saline citrate), treated with RNAse A (20 mg/L) and washed in increasingly stringent SSC solutions at room temperature. Finally sections were washed in 0.1× SSC for 1 h at 67°C or 64°C (respectively) and dehydrated through increasing concentrations of ethanol. Autoradiography was on Biomax MR film (Eastman Kodak Co., NY, USA) for 3-5 days. The autoradiographs were digitised, and relative expression levels were determined by computer-assisted optical densitometry (Image J, Scion Corporation). The average value of 4-6 measurements was calculated from each animal .
A total of 200 ng of amplified RNA from the first amplification round of the microarray analysis was reverse transcribed with Superscript II (Invitrogen, Karlsruhe, Germany) using random hexamer primers (Ambion) according to the manufacturer's protocol. For quality control, a small aliquot of each cDNA was analyzed on an agarose gel.
cDNA of DBA/2J and C57BL/6J mice from the non stressed, 4 h- and 8 h-after stress groups (for DBA/2J: N = 6 in each group; for C57BL/6J: N = 5 for non stressed and 8 h-after stress samples each, N = 6 for the 4 h stress group), was analyzed by qPCR, using the QuantiFast SYBR Green PCR Kit (Qiagen GmbH, Hilden, Germany) according to the manufacturer's instructions. The respective oligonucleotide primers were for GNAi2 (fwd: 5' CACCTCCATCATCCTCTTCC 3', rev: 5' GCACGTGAAGTGCGTGTAGA 3'), APP (fwd: 5' CGGAAGAGATCTCGGAAGTG 3', rev: 5' TGTTCGAACCCACATCTTCA 3') and the house keeping genes glyceraldehydes-3-phosphate dehydrogenase (GAPDH: fwd: 5' CCATCACCATCTTCCAGGAGCGAG 3', rev: 5' GATGGCATGGACTGTGGTCATGAG 3') hypoxanthine guanine phosphoribosyl transferase 1 (HPRT1; fwd: 5' GTCAAGGGCATA TCCAACAACAAAC 3', rev: 5' CCTGCTGGATTACATTAAAGCACTG 3') polymerase (RNA) II (DNA directed) polypeptide B (POLR2B; fwd: 5 'CAAGACAAGGATCATATCT GATGG 3', rev: 5' AGAGTTTAGACGACGCAGGTG 3') or ribosomal protein L13a (RPL13A; fwd: 5' CACTCTGGAGGAGAAACGGAAGG 3', rev: 5' GCAGGCATGAGGC AAACAGTC 3'). Experiments were performed in duplicates on a Lightcycler®2.0 instrument (Roche Diagnostics, Mannheim, Germany) under the following PCR conditions: Initial denaturation at 95°C for 10 min, followed by 40 cycles of denaturation (95°C for 10 s) and a combined annealing and extension phase (60°C for 30 s). At the end of each run, a melting curve (50-95°C with 0.1°C/s) was recorded to ensure the quality of the PCR product. Crossing points (Cp) were calculated by the LightCycler®Software 4.0 (Roche Diagnostics GmbH) using the absolute quantification fit points' method. Threshold and noise band were set to the same level in all runs used in a comparison. Relative gene expression was determined by the 2-ΔΔCT method  using the real PCR efficiency calculated from an external standard curve. Cp were normalized to the housekeeping genes GAPDH, HPRT1, POLR2B or RPL13A). Fold regulation values were calculated relative to the expression mean of basal mice. The calculation Cp (mean (GAPDH, HPRT1, POLR2B or RPL13A)) - Cp (mean (cand. gene)) was used for each animal.
RT-qPCR experiments were analyzed using non-parametric tests (Mann Whitney using SPSS 12.0).
Identifying genes with similar regulation profiles/Clustering Analysis
To identify genes that show similar expression ratios across time points, i.e. genes that might be co-regulated or affecting each other in a common pathway, we used cluster analysis. Cluster analysis allows the grouping of expression profiles with respect to their relative similarity or, in mathematical terms, a "distance". We consider expression profiles (the expression ratios of the different measurements) to be similar--and thus having a small distance--when they fulfil two criteria: 1. show a high absolute (not mean-corrected) correlation, and 2. have either the same or the opposite regulation (up or down) at all corresponding measurements.
Translated into distance between expression profiles this means that expression profiles, that can be scaled (with either a positive or negative factor, but not shifting the zero line) onto each other, have a small distance. This distance measure groups genes with similar regulation patterns as close neighbours in the cluster analysis. For APP and GNAi2 we show the respective neighbourhoods (targeting at least 10 neighbouring genes, but showing all genes in the most distant branch) as depicted by the corresponding dendrograms.
Heim C, Nemeroff CB: The role of childhood trauma in the neurobiology of mood and anxiety disorders: preclinical and clinical studies. Biol Psychiatry. 2001, 49: 1023-1039. 10.1016/S0006-3223(01)01157-X.
Charney D, Manji H: Life Stress, Genes, and Depression: Multiple Pathways Lead to Increased Risk and New Opportunities for Intervention. Science STKE. 2004, re5-10.1126/stke.2252004re5.
de Kloet ER, Joels M, Holsboer F: Stress and the brain: from adaptation to disease. Nat Rev Neurosci. 2005, 6: 463-475. 10.1038/nrn1683.
Holsboer F: The corticosteroid receptor hypothesis of depression 1204. Neuropsychopharmacology. 2000, 23: 477-501. 10.1016/S0893-133X(00)00159-7.
Matthews SG: Early programming of the hypothalamo-pituitary-adrenal axise. Trends in Endocrinology and Metabolism. 2002, 13: 373-380. 10.1016/S1043-2760(02)00690-2.
Darnaudery M, Maccari S: Epigenetic programming of the stress response in male and female rats by prenatal restraint stress. Brain Research Reviews. 2008, 57: 571-585. 10.1016/j.brainresrev.2007.11.004.
Vieau D, Sebaai N, Leeonhardt M, Dutriez-Casteloot I, Molendi-Coste O, Laborie C, Breton C, Deloof S, Lesage J: HPA axis programming by maternal undernutrition in the male rat offspring. Psychoneuroendocrinology. 2007, 32: S16-S20. 10.1016/j.psyneuen.2007.03.014.
Murgatroyd C, Patchev AV, Wu Y, Micale V, Bockmuhl Y, Fischer D, Holsboer F, Wotjak CT, Almeida OF, Spengler D: Dynamic DNA methylation programs persistent adverse effects of early-life stress. Nat Neurosci. 2009, 12: 1559-1566. 10.1038/nn.2436.
De Bellis MD, Baum AS, Birmaher B, Keshavan MS, Eccard CH, Boring AM, Jenkins FJ, Ryan ND: Developmental traumatology Part I: Biological stress systems. Biological Psychiatry. 1999, 45: 1259-1270. 10.1016/S0006-3223(99)00044-X.
Heuser I, Lammers CH: Stress and the brain. Neurobiology of Aging. 2003, 24: S69-S76. 10.1016/S0197-4580(03)00048-4.
Armario A: The Hypothalamic-Pituitary-Adrenal Axis: What can it Tell us About Stressors. CNS Neurol Disord Drug Targets. 2006, 5: 485-501. 10.2174/187152706778559336.
Raadsheer FC, Sluiter AA, Ravid R, Tilders FJH, Swaab DF: Localization of Corticotropin-Releasing Hormone (Crh) Neurons in the Paraventricular Nucleus of the Human Hypothalamus - Age-Dependent Colocalization with Vasopressin. Brain Research. 1993, 615: 50-62. 10.1016/0006-8993(93)91113-7.
Gillies G, Lowry PJ: Corticotropin-Releasing Hormone and Its Vasopressin Component. Frontiers in Neuroendocrinology. 1982, 7: 45-75.
Rivier C, Vale W: Interaction of Corticotropin-Releasing Factor and Arginine Vasopressin on Adrenocorticotropin Secretion Invivo. Endocrinology. 1983, 113: 939-942. 10.1210/endo-113-3-939.
Landgraf R: The involvement of vasopressin system in stress-related disorders. CNS Neurol Disord Drug Targets. 2006, 5: 167-179. 10.2174/187152706776359664.
Engelmann M, Landgraf R, Wotjak CT: The hypothalamic-neurohypophysial system regulates the hypothalamic-pituitary-adrenal axis under stress: An old concept revisited. Frontiers in Neuroendocrinology. 2004, 25: 132-149. 10.1016/j.yfrne.2004.09.001.
Swaab DF, Bao AM, Lucassen PJ: The stress system in the human brain in depression and neurodegeneration. Ageing Research Reviews. 2005, 4: 141-194. 10.1016/j.arr.2005.03.003.
Raadsheer FC, Hoogendijk WJG, Stam FC, Tilders FJH, Swaab DF: Increased Numbers of Corticotropin-Releasing Hormone Expressing Neurons in the Hypothalamic Paraventricular Nucleus of Depressed-Patients. Neuroendocrinology. 1994, 60: 436-444. 10.1159/000126778.
Raadsheer FC, Vanheerikhuize JJ, Lucassen PJ, Hoogendijk WJG, Tilders FJH, Swaab DF: Corticotropin-Releasing Hormone Messenger-RNA Levels in the Paraventricular Nucleus of Patients with Alzheimers-Disease and Depression. American Journal of Psychiatry. 1995, 152: 1372-1376.
Ising M, Holsboer F: Genetics of stress response and stress-related disorders. Dialogues Clin Neurosci. 2006, 8 (4): 433-444. Ref Type: Generic
Tsolakidou A, Trumbach D, Panhuysen M, Putz B, Deussing J, Wurst W, Sillaber I, Holsboer F, Rein T: Acute stress regulation of neuroplasticity genes in mouse hippocampus CA3 area--possible novel signalling pathways. Mol Cell Neurosci. 2008, 38: 444-452. 10.1016/j.mcn.2008.04.005.
Zacharko RM, Lalonde GT, Kasian M, Anisman H: Strain-Specific Effects of Inescapable Shock on Intracranial Self-Stimulation from the Nucleus-Accumbens. Brain Research. 1987, 426: 164-168. 10.1016/0006-8993(87)90436-7.
Zacharko RM, Gilmore W, Macneil G, Kasian M, Anisman H: Stressor Induced Variations of Intracranial Self-Stimulation from the Mesocortex in Several Strains of Mice. Brain Research. 1990, 533: 353-357. 10.1016/0006-8993(90)91363-L.
Ohl F, Roedel A, Binder E, Holsboer F: Impact of high and low anxiety on cognitive performance in a modified hole board test in C57BL/6 and DBA/2 mice. European Journal of Neuroscience. 2003, 17: 128-136. 10.1046/j.1460-9568.2003.02436.x.
Rossi-Arnaud C, Ammassari-Teule M: What do comparative studies of inbred mice add to current investigations on the neural basis of spatial behaviors?. Experimental Brain Research. 1998, 123: 36-44. 10.1007/s002210050542.
Crowley JJ, Lucki I: Opportunities to discover genes regulating depression and antidepressant response from rodent behavioral genetics. Current Pharmaceutical Design. 2005, 11: 157-169. 10.2174/1381612053382278.
Sillaber I, Panhuysen M, Henniger MSH, Ohl F, Kuhne C, Pütz B, Pohl T, Deussing JM, Paez-Pereda M, Holsboer F: Profiling of behavioral changes and hippocampal gene expression in mice chronically treated with the SSRI paroxetine. Psychopharmacology. 2008, 200: 557-572. 10.1007/s00213-008-1232-6.
Porsolt RD, Lepichon M, Jalfre M: Depression - New Animal-Model Sensitive to Antidepressant Treatments. Nature. 1977, 266: 730-732. 10.1038/266730a0.
Giardino L, Puglisi-Allegra S, Ceccatelli S: CRH-R1 mRNA expression in two strains of inbred mice and its regulation after repeated restraint stress. Brain Res Mol Brain Res. 1996, 40: 310-314. 10.1016/0169-328X(96)00092-7.
Lee HC, Chang DE, Yeom M, Kim GH, Choi KD, Shim I, Lee HJ, Hahm DH: Gene expression profiling in hypothalamus of immobilization-stressed mouse using cDNA microarray. Brain Res Mol Brain Res. 2005, 135: 293-300. 10.1016/j.molbrainres.2004.11.016.
Reyes TM, Walker JR, DeCino C, Hogenesch JB, Sawchenko PE: Categorically distinct acute stressors elicit dissimilar transcriptional profiles in the paraventricular nucleus of the hypothalamus. Journal of Neuroscience. 2003, 23: 5607-5616.
Thoeringer CK, Sillaber I, Roedel A, Erhardt A, Mueller MB, Ohl F, Holsboer F, Keck ME: The temporal dynamics of intrahippocampal corticosterone in response to stress-related stimuli with different emotional and physical load: an in vivo microdialysis study in C57BL/6 and DBA/2 inbred mice. Psychoneuroendocrinology. 2007, 32: 746-757. 10.1016/j.psyneuen.2007.05.005.
Belzung C, Philippot P: Anxiety from a phylogenetic perspective: is there a qualitative difference between human and animal anxiety?. Neural Plast. 2007, 59676-
Lucki I, Dalvi A, Mayorga AJ: Sensitivity to the effects of pharmacologically selective antidepressants in different strains of mice. Psychopharmacology. 2001, 155: 315-322. 10.1007/s002130100694.
Binder EB, Salyakina D, Lichtner P, Wochnik GM, Ising M, Putz B, Papiol S, Seaman S, Lucae S, Kohli MA, Nickel T, Kunzel HE, Fuchs B, Majer M, Pfennig A, Kern N, Brunner J, Modell S, Baghai T, Deiml T, Zill P, Bondy B, Rupprecht R, Messer T, Kohnlein O, Dabitz H, Bruckl T, Muller N, Pfister H, Lieb R, Mueller JC, Lohmussaar E, Strom TM, Bettecken T, Meitinger T, Uhr M, Rein T, Holsboer F, Muller-Myhsok B: Polymorphisms in FKBP5 are associated with increased recurrence of depressive episodes and rapid response to antidepressant treatment. Nat Genet. 2004, 36: 1319-1325. 10.1038/ng1479.
Wochnik GM, Ruegg J, Abel GA, Schmidt U, Holsboer F, Rein T: FK506-binding proteins 51 and 52 differentially regulate dynein interaction and nuclear translocation of the glucocorticoid receptor in mammalian cells. J Biol Chem. 2005, 280: 4609-4616. 10.1074/jbc.M407498200.
Uhr M, Tontsch A, Namendorf C, Ripke S, Lucae S, Ising M, Dose T, Ebinger M, Rosenhagen M, Kohli M, Kloiber S, Salyakina D, Bettecken T, Specht M, Putz B, Binder EB, Muller-Myhsok B, Holsboer F: Polymorphisms in the drug transporter gene ABCB1 predict antidepressant treatment response in depression. Neuron. 2008, 57: 203-209. 10.1016/j.neuron.2007.11.017.
Ising M, Lucae S, Binder EB, Bettecken T, Uhr M, Ripke S, Kohli MA, Hennings JM, Horstmann S, Kloiber S, Menke A, Bondy B, Rupprecht R, Domschke K, Baune BT, Arolt V, Rush AJ, Holsboer F, Müller-Myhsok B: A genomewide association study points to multiple loci that predict antidepressant drug treatment outcome in depression. Arch Gen Psychiatry. 2009, 66: 966-975. 10.1001/archgenpsychiatry.2009.95.
Herr AS, Tsolakidou AF, Yassouridis A, Holsboer F, Rein T: Antidepressants differentially influence the transcriptional activity of the glucocorticoid receptor in vitro. Neuroendocrinology. 2003, 78: 12-22. 10.1159/000071701.
Montkowski A, Barden N, Wotjak C, Stec I, Ganster J, Meaney M, Engelmann M, Reul JM, Landgraf R, Holsboer F: Long-term antidepressant treatment reduces behavioural deficits in transgenic mice with impaired glucocorticoid receptor function. J Neuroendocrinol. 1995, 7: 841-845. 10.1111/j.1365-2826.1995.tb00724.x.
Neves SR, Ram PT, Iyengar R: G protein pathways. Science. 2002, 296: 1636-1639. 10.1126/science.1071550.
Khan ZU, Gutierrez A: Distribution of C-terminal splice variant of G alpha(i2) in rat and monkey brain. Neuroscience. 2004, 127: 833-843. 10.1016/j.neuroscience.2004.05.017.
Neves SR, Ram PT, Iyengar R: G protein pathways. Science. 2002, 296: 1636-1639. 10.1126/science.1071550.
Arinze IJ, Kawai Y: Transcriptional activation of the human G alpha(i2) gene promoter through nuclear factor-kappa B and antioxidant response elements. Journal of Biological Chemistry. 2005, 280: 9786-9795. 10.1074/jbc.M414006200.
Turner PR, O'Connor K, Tate WP, Abraham WC: Roles of amyloid precursor protein and its fragments in regulating neural activity, plasticity and memory. Prog Neurobiol. 2003, 70: 1-32. 10.1016/S0301-0082(03)00089-3.
Panegyres PK: The functions of the amyloid precursor protein gene. Reviews in the Neurosciences. 2001, 12: 1-39.
Roch JM, Masliah E, Rochlevecq AC, Sundsmo MP, Otero DAC, Veinbergs I, Saitoh T: Increase of Synaptic Density and Memory Retention by A Peptide Representing the Trophic Domain of the Amyloid Beta/A4 Protein-Precursor. Proceedings of the National Academy of Sciences of the United States of America. 1994, 91: 7450-7454. 10.1073/pnas.91.16.7450.
Mucke L, Abraham CR, Masliah E: Neurotrophic and neuroprotective effects of hAPP in transgenic mice. Neurobiology of Alzheimer'S Disease. 1996, 777: 82-88.
Palacios G, Mengod G, Frey P, Palacios JM: Rapid Increase in Amyloid Precursor Protein Immunoreactivity in the Supraoptic and Paraventricular Nuclei of the Rat Hypothalamus After Osmotic-Stress. Neuroreport. 1995, 6: 265-268. 10.1097/00001756-199501000-00010.
Huber G, Bailly Y, Martin JR, Mariani J, Brugg B: Synaptic beta-amyloid precursor proteins increase with learning capacity in rats. Neuroscience. 1997, 80: 313-320. 10.1016/S0306-4522(97)00120-6.
Gralle M, Ferreira ST: Structure and functions of the human amyloid precursor protein: The whole is more than the sum of its parts. Progress in Neurobiology. 2007, 82: 11-32. 10.1016/j.pneurobio.2007.02.001.
Mattson MP, Duan WZ, Chan SL, Cheng AW, Haughey N, Gary DS, Guo ZH, Lee JW, Furukawa K: Neuroprotective and neurorestorative signal transduction mechanisms in brain aging: modification by genes, diet and behavior. Neurobiology of Aging. 2002, 23: 695-705. 10.1016/S0197-4580(02)00025-8.
Siman R, Card JP, Nelson RB, Davis LG: Expression of Beta-Amyloid Precursor Protein in Reactive Astrocytes Following Neuronal Damage. Neuron. 1989, 3: 275-285. 10.1016/0896-6273(89)90252-3.
Postina R, Schroeder A, Dewachter I, Bohl J, Schmitt U, Kojro E, Prinzen C, Endres K, Hiemke C, Blessing M, Flamez P, Dequenne A, Godaux E, van Leuven F, Fahrenholz F: A disintegrin-metalloproteinase prevents amyloid plaque formation and hippocampal defects in an Alzheimer disease mouse model. Journal of Clinical Investigation. 2004, 113: 1456-1464.
Kirazov L, Loffler T, Schliebs R, Bigl V: Glutamate-stimulated secretion of amyloid precursor protein from cortical rat brain slices. Neurochemistry International. 1997, 30: 557-563. 10.1016/S0197-0186(96)00119-2.
Nitsch RM, Deng A, Wurtman RJ, Growdon JH: Metabotropic glutamate receptor subtype mGluR1 alpha stimulates the secretion of the amyloid beta-protein precursor ectodomain. Journal of Neurochemistry. 1997, 69: 704-712. 10.1046/j.1471-4159.1997.69020704.x.
Churchill GA: Fundamentals of experimental design for cDNA microarrays. Nature Genetics. 2002, 32: 490-495. 10.1038/ng1031.
Deussing JM, Kuehne C, Puetz B, Panhuysen M, Breu J, Stenzel-Poore M, Holsboer F, Wurst W: Expression profiling identifies the CRH/CRH-R1 system as a modulator of neurovascular gene activity. J Cereb Blood Flow Metab. 2007, 27: 1476-1495. 10.1038/sj.jcbfm.9600451.
Hackl H, Cabo FS, Sturn A, Wolkenhauer O, Trajanoski Z: Analysis of DNA microarray data. Current Topics in Medicinal Chemistry. 2004, 4: 1357-1370. 10.2174/1568026043387773.
Yang YH, Dudoit S, Luu P, Lin DM, Peng V, Ngai J, Speed TP: Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. Nucleic Acids Res. 2002, 30: e15-10.1093/nar/30.4.e15.
Schmidt MV, Deussing JM, Oitzl MS, Ohl F, Levine S, Wurst W, Holsboer F, Müller MB, de Kloet ER: Differential disinhibition of the neonatal hypothalamic-pituitary-adrenal axis in brain-specific CRH receptor 1-knockout mice. European Journal of Neuroscience. 2006, 24: 2291-2298. 10.1111/j.1460-9568.2006.05121.x.
Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(T)(-Delta Delta C) method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.
The authors thank Barbara Berning for excellent technical support on the in situ hybridisation. The contribution of the Helmholtz Centre was supported by the Helmholtz Alliance for Mental Health in an Ageing Society (HA-215) and by the federal ministry for education and research (BMBF NGFNplus 01GS08155, 01GS08151).
AT performed the animal experiments (including hormone measurements), micropuncture, RNA preparation, microarray, Ontology Groups Classification, in situ hybridisation, molecular pathways building, interpretation of data, and manuscript writing. LC performed the real time PCR experiments, BP the microarray analysis as well as the clustering analysis, and DT the PathwayStudio analyses. MP contributed to microarray analyses and ontology group classification. JMD supervised microarray experiments, WW PathwayStudio analyses. IS contributed to project design, RL to real time PCR supervision, and FH to shaping the project concept. TR designed the project and wrote the manuscript together with AT.
All the authors have read and approved the final manuscript.
Electronic supplementary material
Authors’ original submitted files for images
About this article
- Forced Swimming
- Inbred Mouse Strain
- Signal Transduction Molecule
- Guanine Nucleotide Binding Protein
- Hypothalamic Paraventricular Nucleus