Gene expression profiling in the stress control brain region hypothalamic paraventricular nucleus reveals a novel gene network including Amyloid beta Precursor Protein
© Tsolakidou et al; licensee BioMed Central Ltd. 2010
Received: 23 March 2010
Accepted: 8 October 2010
Published: 8 October 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
Differentially expressed genes in the PVN between mouse strains (examples, +: higher, -: lower expressed in DBA/2J vs C57BL/6J)
DNA segment, Chr 10, Brigham & Women's Genetics 1379 expressed
LUC7-like 2 (S. cerevisiae)
Mitochondrial ribosomal protein S18C
Mitogen activated protein kinase kinase kinase 12
Pituitary tumor-transforming 1
Histocompatibility 2, Q region locus 1
Glutamate-ammonia ligase (glutamine synthase)
Sparc/osteonectin, cwcv and kazal-like domains proteoglycan 3
Gamma-aminobutyric acid (GABA-A) receptor, subunit alpha 1
Amyloid beta (A4) precursor protein
Guanine nucleotide binding protein, alpha o
Potassium inwardly-rectifying channel, subfamily J, member 9
RIKEN cDNA 1700112M01 gene
ACTH and gene expression profiles after forced swimming
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 genes that are up- or down-regulated 4 h after stress in C57BL/6J mice
Alpha 1 microglobulin/bikunin
NADH dehydrogenase (ubiquinone) 1 alpha subcomplex, 8
P21 (CDKN1A)-activated kinase 2
Ca < 2+ > dependent activator protein for secretion
Mitogen-activated protein kinase kinase kinase kinase 3
Calcium binding protein 39
Polymerase (DNA directed), delta 1, catalytic subunit
Calcium channel, voltage-dependent, P/Q type, alpha 1A subunit
BCL2-associated athanogene 4
Opioid growth factor receptor-like 1
Examples of Genes up- or down-regulated 8 h after stress in C57BL/6J mice
A disintegrin and metalloprotease domain 10
Insulin induced gene 1
Brain-specific angiogenesis inhibitor 1-associated protein 2
LUC7-like 2 (S. cerevisiae)
Heat shock protein 1, alpha
Serine/arginine-rich protein specific kinase 2
Protein tyrosine phosphatase 4a2
Gap junction membrane channel protein alpha 12
Sprouty homolog 2 (Drosophila)
Proteolipid protein (myelin) 1
Kinesin family member 3B
Expressed in non-metastatic cells 1, protein
Guanine nucleotide binding protein, alpha o
Isocitrate dehydrogenase 3 (NAD+) beta
Sodium channel, type IV, beta
Regulator of G-protein signaling 2
ATP-binding cassette, sub-family F (GCN20), member 2
Examples of Genes up- or down-regulated 4 h after stress in DBA/2J mice
Hematopoietically expressed homeobox
Protein phosphatase 3, catalytic subunit, beta isoform
Pre-B-cell colony-enhancing factor 1
Nuclear receptor subfamily 3, group C, member 1
Ribosomal protein S6
Bone morphogenetic protein 6
Protein kinase C, zeta
Thyrotropin releasing hormone
Internexin neuronal intermediate filament protein, alpha
Nuclear factor of activated T-cells, cytoplasmic, calcineurin-dependent 1
Guanine nucleotide binding protein, alpha inhibiting 2
Mitogen activated protein kinase 3
Vesicle-associated membrane protein 2
Microtubule-associated protein 1 light chain 3 alpha
Examples of Genes up- or down-regulated 8 h after stress in DBA/2J mice
Homer homolog 3 (Drosophila)
Zinc finger, A20 domain containing 2
Chromodomain helicase DNA binding protein 4
Solute carrier family 6 (neurotransmitter transporter, betaine/GABA), member 12
Protein regulator of cytokinesis 1
Toll-like receptor 7
Transcription factor 7-like 2, T-cell specific, HMG-box
Heat shock protein 1A
Calcium channel, voltage-dependent, gamma subunit 1
Nuclear factor I/A
Telomerase binding protein, p23
Amyloid beta (A4) precursor protein
Cyclin-dependent kinase inhibitor 1B (P27)
A disintegrin and metalloprotease domain 10
A disintegrin and metalloprotease domain 10
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
Ontologically sorted genes regulated 4 h and 8 h after stress
C57BL/6 Stress-4 h
C57BL/6 Stress-8 h
DBA/2J Stress-4 h
DBA/2J Stress-8 h
mitochondrial, biosynthetic and metabolic proteins
vesicular transport and components of cytoskeleton
receptors and signal transduction molecules
transcription and mRNA molecules
ion channels and ion transport
molecules with other function
Number of genes regulated 4 h and 8 h after stress in C57BL/6J or DBA/2J
Stressed - 4 h
common S4 h - S8 h
Stressed - 8 h
common between C57BL/6 and DBA/2J
Possible signalling pathways elicited after forced swimming in the PVN
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
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 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.
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).
- 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.PubMedView ArticleGoogle Scholar
- 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.Google Scholar
- 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.PubMedView ArticleGoogle Scholar
- Holsboer F: The corticosteroid receptor hypothesis of depression 1204. Neuropsychopharmacology. 2000, 23: 477-501. 10.1016/S0893-133X(00)00159-7.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Heuser I, Lammers CH: Stress and the brain. Neurobiology of Aging. 2003, 24: S69-S76. 10.1016/S0197-4580(03)00048-4.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Gillies G, Lowry PJ: Corticotropin-Releasing Hormone and Its Vasopressin Component. Frontiers in Neuroendocrinology. 1982, 7: 45-75.Google Scholar
- 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.PubMedView ArticleGoogle Scholar
- Landgraf R: The involvement of vasopressin system in stress-related disorders. CNS Neurol Disord Drug Targets. 2006, 5: 167-179. 10.2174/187152706776359664.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Ising M, Holsboer F: Genetics of stress response and stress-related disorders. Dialogues Clin Neurosci. 2006, 8 (4): 433-444. Ref Type: GenericPubMed CentralPubMedGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Porsolt RD, Lepichon M, Jalfre M: Depression - New Animal-Model Sensitive to Antidepressant Treatments. Nature. 1977, 266: 730-732. 10.1038/266730a0.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Belzung C, Philippot P: Anxiety from a phylogenetic perspective: is there a qualitative difference between human and animal anxiety?. Neural Plast. 2007, 59676-Google Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Neves SR, Ram PT, Iyengar R: G protein pathways. Science. 2002, 296: 1636-1639. 10.1126/science.1071550.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Neves SR, Ram PT, Iyengar R: G protein pathways. Science. 2002, 296: 1636-1639. 10.1126/science.1071550.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Panegyres PK: The functions of the amyloid precursor protein gene. Reviews in the Neurosciences. 2001, 12: 1-39.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- Mucke L, Abraham CR, Masliah E: Neurotrophic and neuroprotective effects of hAPP in transgenic mice. Neurobiology of Alzheimer'S Disease. 1996, 777: 82-88.Google Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Churchill GA: Fundamentals of experimental design for cDNA microarrays. Nature Genetics. 2002, 32: 490-495. 10.1038/ng1031.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar