About the Authors:
Mirjam Bunck
Contributed equally to this work with: Mirjam Bunck, Ludwig Czibere
Affiliation: Max Planck Institute of Psychiatry, Munich, Germany
Ludwig Czibere
Contributed equally to this work with: Mirjam Bunck, Ludwig Czibere
* E-mail: [email protected]
Affiliation: Max Planck Institute of Psychiatry, Munich, Germany
Charlotte Horvath
Affiliation: Max Planck Institute of Psychiatry, Munich, Germany
Cornelia Graf
Affiliation: Max Planck Institute of Psychiatry, Munich, Germany
Elisabeth Frank
Current address: University of Wollongong, Wollongong, New South Wales, Australia
Affiliation: Max Planck Institute of Psychiatry, Munich, Germany
Melanie S. Keßler
Affiliation: Max Planck Institute of Psychiatry, Munich, Germany
Chris Murgatroyd
Affiliation: Max Planck Institute of Psychiatry, Munich, Germany
Bertram Müller-Myhsok
Affiliation: Max Planck Institute of Psychiatry, Munich, Germany
Mariya Gonik
Affiliation: Max Planck Institute of Psychiatry, Munich, Germany
Peter Weber
Affiliation: Max Planck Institute of Psychiatry, Munich, Germany
Benno Pütz
Affiliation: Max Planck Institute of Psychiatry, Munich, Germany
Patrik Muigg
Affiliation: University of Innsbruck, Innsbruck, Austria
Markus Panhuysen
Affiliation: Affectis Pharmaceuticals, Martinsried, Germany
Nicolas Singewald
Affiliation: University of Innsbruck, Innsbruck, Austria
Thomas Bettecken
Affiliation: Max Planck Institute of Psychiatry, Munich, Germany
Jan M. Deussing
Affiliations Max Planck Institute of Psychiatry, Munich, Germany, Helmholtz Zentrum, Institute of Developmental Genetics, Munich/Neuherberg, Germany
Florian Holsboer
Affiliation: Max Planck Institute of Psychiatry, Munich, Germany
Dietmar Spengler
Affiliation: Max Planck Institute of Psychiatry, Munich, Germany
Rainer Landgraf
Affiliation: Max Planck Institute of Psychiatry, Munich, Germany
Introduction
Despite decades of intensive research, anxiety and depression disorders remain two of the most common and debilitating psychiatric illnesses. Their causes and appropriate medication are still elusive [1], and the number of affected people continues to increase [2], [3]. Depression is a highly heritable disorder; epidemiologic studies have shown that about 40–50% of the risk for depression is genetically determined [4], [5] and that depression is in 60% of all cases comorbid with anxiety disorders [6], [7]. These findings are strongly supported by familial and twin studies [8]. However, the considerable evidence for heritability has not been accompanied by satisfactory progress towards identifying genetic risk factors underlying anxiety and depression. An improved understanding of the genetic contributions would have significant implications for the diagnosis and treatment of neuropsychiatric disorders.
In order to identify shared genotype/phenotype relationships in animals and humans, the same gene should affect analogous phenotypes in both species. Thus, a useful approach should depend on conserved gene function, the presence of functional polymorphisms in this gene in both species and the choice of an appropriate analogous phenotype in both the model species and in humans [9]. The gene for the nonapeptide vasopressin (AVP) fulfills all these requirements. It is evolutionary conserved as both antidiuretic hormone after peripheral secretion and anxiogenic stress-related neuropeptide after central release [10]–[12]. Moreover, the presence of functional polymorphisms and subsequent deficits in expression/processing in rodents and patients leads to signs of both central diabetes insipidus and emotional alterations, e.g. decreased anxiety-related behavior [13]–[16]. Finally, trait anxiety with the genetic predisposition towards either high (HAB) or low (LAB) anxiety-related behavior is a heuristically and clinically relevant phenotype, with a high degree of face, construct and predictive validity [17].
Preclinical studies, including findings from the HAB/LAB rat model [18]–[21] and voles [22], as well as clinical observations [23] support a direct involvement of centrally released AVP in anxiety/depression-like behaviors and disorders. To make use of the advantages of the full spectrum of molecular-genetic approaches, CD1 mice were selectively inbred for either low or high anxiety-related behavior in analogy to the HAB/LAB rat model, while maintaining a high degree of similarity in non-selected traits. Along with their divergence in anxiety-related behavior, which is reflected in a variety of distinct behavioral paradigms, these mouse lines also exhibit differences in tests of other metrics, including depression-like and explorative behaviors [24]. The final goal of our selective breeding strategy is to identify candidate genes of anxiety and to prove their causal involvement in the anxiety phenotype, while simultaneously reducing the possible influence of genetic drift and the concomitant risk of indicating false associations between the gene of interest and a given phenotype, a risk particularly high in inbred lines or strains. Therefore, a series of measures have been employed in order to minimize the risk related to genetic drift: (i) in each selectively bred line, several independent sublines are run in parallel; (ii) identified genetic polymorphisms are tested in a freely-segregating F2 panel and are accepted to be causal if they co-segregate with the phenotype; (iii) in addition to the divergent lines, normal anxiety (NAB) mice are used to represent intermediate anxiety-related behavior; and (iv) replication is important in terms of the impact of founder effects and genetic drift and should be run concurrently with the main breeding protocol in the same or another species, thus also facilitating complementary inter-species genetics [9].
Therefore, if Avp gene expression is causally related to the HAB vs. LAB phenotype, genetic polymorphisms contributing to this phenomenon should be detectable not only in selectively bred rats but also in mice. Indeed, LAB mice display a deficit in hypothalamic AVP [25]. In light of the functional impact of polymorphisms in the Avp regulatory regions for both gene expression and anxiety-related behavior in rats [18], we here addressed the question whether polymorphisms in the Avp promoter of LAB mice may be identified and contribute to the neuroendocrine and behavioral divergence demonstrated in HAB vs. LAB mice. This will provide new insights into the genetic predisposition that shapes emotionality in both inter-individual variation and psychopathology.
Materials and Methods
Animals
All animals tested were bred in the animal facilities of the Max Planck Institute of Psychiatry in Munich (Germany) as described previously [24]. Briefly, >250 animals from >25 litters of outbred Swiss CD1 mice purchased from Charles River (Sulzfeld, Germany) were used as a starting point for selective and bidirectional breeding for anxiety-related behavior on the elevated plus-maze (EPM) at the age of seven weeks. Males and females that spent either the least or most time on the open arms of the EPM were mated to establish the HAB and LAB lines, respectively. This intra-strain approach was chosen to ensure that the HAB/LAB lines show a maximum divergence in the selected trait, while maintaining a high degree of similarity in non-selected traits. The animals were routinely tested at the age of seven weeks, with HAB mice spending <10% and LAB >50% of the total test time on the open arms of the EPM. NAB mice were bred for intermediate anxiety-related behavior. As >80% of CD1 mice spend between 25% and 35% of their time on the open arms of the EPM, this range was chosen for the selection of NAB mice without any overlapping either with LAB or with HAB animals. While CD1 mice of the parental generation were used as NAB controls by Krömer et al. [24] and Kessler et al. [25], we then decided to start inbreeding them in parallel with HAB and LAB mice, aiming at further reducing all variables that are unrelated to anxiety, such as slight differences in body weight between outbred vs. inbred animals.
All experiments were carried out on adult male or female HAB, LAB, NAB, HABxLAB F1 intercross and F2 mice. In the present study, HAB and LAB mice of generations 12–23 and NAB mice of generation 1 were used. Importantly, in a wide variety of tests and parameters, the scores of NAB and F1 intercross controls were found to be similar, if not identical, and between the scores of HAB and LAB mice ([24]; unpublished results).
The animal studies were both approved by local authorities and conducted according to current regulations for animal experimentation in Germany and the European Union (European Communities Council Directive 86/609/EEC).
EPM test
The EPM consisted of a plus-shaped platform elevated 37 cm above the floor, with two open (30×5 cm), two closed (30×5×15 cm) arms and a connecting central zone (5×5 cm) made of grey PVC. The open arms were lit by white light of 300 lux intensity, the neutral zone of 60 lux intensity and the closed arms of 5 lux intensity. The maze was cleaned with water containing a detergent and dried before each trial.
Animals were directly transferred from their homecage to the test apparatus, starting with the animal placed on the central part facing one closed arm. During the 5-min exposure, various parameters were recorded by means of a video/computer system. Analysis was performed post hoc by an observer blind to the line using Plus-maze V2.0 (Ernst Fricke, Munich, Germany) [24], [26].
Tail-suspension test
In the tail-suspension test, animals were suspended at the end of their tail to a bar at a height of 37 cm above the floor with an adhesive tape [27]. Four animals were tested in parallel during each test. The behavior was videotaped and the duration of total immobility (immobility of all four extremities) scored by an observer blind to the line using Eventlog 1.0 (EMCO Software, Reykjavik, Iceland).
Forced swim test
In the forced swim test each mouse was placed for 6 min into a glass cylinder (height: 23 cm; diameter: 11 cm) containing 15 cm of water at 22–23°C. For each videotaped trial, immobility time (floating) was analyzed by an observer blind to the line using Eventlog 1.0 (EMCO Software). A mouse was judged immobile when it stopped any slight movements, except those that are necessary to keep its head above water [28].
Open field
The round open field, having a diameter of 60 cm, was made of black and grey PVC and was optically divided into an inner (diameter: 30 cm) and an outer area to allow quantification of the time spent in each compartment. To devise the test situation as less aversive as possible and to make the test suitable to measure locomotor activity independent of anxiety, the arena was illuminated with dim light of 60 lux. Before each trial, the chamber was cleaned with water containing a detergent. The animals' behavior was videotaped during the 5-min test time and the parameter path length was analyzed with the tracking software ANY-maze version 4.30 (Stoelting, Wood Dale, IL, USA).
Gene and peptide expression profiling
If not indicated otherwise, all animals used for gene and peptide expression profiling or receptor radioautography were sacrificed by decapitation under isoflurane anesthesia (Abbott, Wiesbaden, Germany), between 9a.m. and 12p.m. The brains were removed, snap-frozen in N-methylbutane (Carl Roth, Karlsruhe, Germany) and stored at −20°C or −80°C.
Tissue sampling for microarray and quantitative real-time PCR.
The tissue micropunching technique [29] was applied for acquiring tissue samples from specific regions with micropunchers of 0.5 mm and 1 mm diameter (Fine Science Tools GmbH, Heidelberg, Germany). Punches were collected from 200 µm sections between −0.56 mm and −0.96 mm from bregma, medially 0.8 mm above the ventral tissue edge around the dorsal end of the 3rd ventricle to obtain the PVN (Ø = 1 mm) and bilaterally from the optic tract (Ø = 0.5 mm) in order to acquire tissue from the SON. All coordinates were based on the mouse brain atlas [30]. Tissue punches were kept at −80°C until further processing.
RNA extraction and reverse transcription from tissue punches.
Tissue punches were homogenized with a pipette in 300 µl Trizol (Tri Reagent, Sigma-Aldrich). Then 30 µl bidistilled water, 1 µl linear acrylamide (5 mg/ml, Ambion, Austin, TX) and 60 µl chloroform were added before continuing RNA extraction with a standard protocol, using n-propanol for RNA precipitation overnight. The yield of total RNA was 0.3–1.5 µg.
All isolated RNA was further processed for microarray analysis using the Ambion Amino Allyl Message Amp aRNA kit (Ambion) in two rounds of amplification. For qPCR, not more than 1 µg of total RNA was reverse transcribed with Superscript II (Invitrogen, Karlsruhe, Germany) after DNAse treatment. All steps required for reverse transcription were performed according to the manufacturer's protocol. For quality control, a small aliquot of cDNA was analyzed on an agarose gel.
Microarray analysis.
Whole-genome gene expression data were acquired using a dual-color experimental design (cyanine 3 and cyanine 5) on an in-house microarray platform (MPI24K) with pooled samples (six LAB vs. six HAB male mouse brain pools) using ten technical replicates including a dye-swap of the hybridized samples on five microarray slides. Pretreatment of microarray probes, assay hybridization, scanning, data evaluation and statistical analysis were performed as described previously [31], [32]. Briefly, microarray slides were scanned in a PerkinElmer ScanArray 4000 (PerkinElmer Life and Analytical Sciences, Shelton, CT) and quantified using QuantArray software (GSI Lumonics, Billerica, MA) with a fixed-circle quantification protocol. Data normalization was performed in multiple steps based on the raw fluorescence intensities ensuring the elimination of artifacts resulting from the unequal distribution of probes or from unbalanced fluorescence intensities within one microarray slide.
Quantitative PCR analysis.
cDNA of male or female, LAB, HAB or NAB mice was analyzed by quantitative real-time PCR (qPCR), using the LightCycler® FastStart DNA MasterPLUS SYBR Green I reagent (Roche Diagnostics, Mannheim, Germany) or the QuantiFast SYBR Green PCR Kit (Qiagen, Hilden, Germany) according to the manufacturer's instructions and the respective oligonucleotide primers for Avp, Avpr1b, Oxtr. For reaction conditions and primer sequences see Table 1. Experiments were performed in duplicates with the Lightcycler®2.0 instrument (Roche Diagnostics) under the following PCR conditions: initial denaturation at 95°C for 10 min, then 40 cycles of denaturation (95°C for 10 s), annealing (56–65°C for 4–5 s) and elongation (72°C, 7–13 s) for the Roche-kit. Similar conditions were applied for the Qiagen-kit, only the annealing and extension phase were standardly substituted by a combined annealing/extension phase at 60°C for 30 s. At the end of every run, a melting curve (50–95°C with 0.1°C/s) was generated in order to ensure the homogenity of the PCR product.
[Figure omitted. See PDF.]
Table 1. Forward and reverse primers used for quantitative PCR and their specific annealing temperatures.
https://doi.org/10.1371/journal.pone.0005129.t001
Crossing points were calculated by the LightCycler®Software 4.0 (Roche Diagnostics) using the absolute quantification fit points method. Threshold and noise band were set to the same level in all compared runs.
Relative gene expression was determined by the 2−ΔΔCT method [33] using the real PCR efficiency calculated from an external standard curve. The crossing point was normalized to the housekeeping genes Gapdh, Hprt and Atp2b1 (Table 1), respectively and values were calculated relative to the expression mean of LAB mice.
Avp and Oxytocin mRNA in situ hybridization.
Brains were sectioned into 14 µm slices with a cryostat (Microm HM 500, Microm, Walldorf, Germany). Three sets of sections at PVN level were taken and used for Avp or Oxytocin (Oxt) in situ hybridization (ISH) according to established protocols [20], [34].
Briefly, sections were dehydrated in increasing ethanol concentrations, delipidated with chloroform, rinsed in ethanol and air-dried afterwards. A highly specific 48-bp oligonucleotide sequence directed against the last 16 amino acid coding region of the glycoprotein that AVP does not share with OXT (5′ GGG CTT GGC AGA ATC CAC GGA CTC TTG TGT CCC AGC CAG CTG TAC CAG 3′; [35], [36]) was used for hybridization with Avp. Detection of Oxt mRNA was performed by hybridization of a similarly specific 48-bp oligonucleotide (5′ CTC GGA GAA GGC AGA CTC AGG GTC GCA GGC GGG GTC GGT GCG GCA GCC 3′; [35]). The oligonucleotides were labeled by using [35S]ATP (NEN DuPont, Bad Homburg, Germany) and terminal transferase (Tdt, Boehringer, Mannheim, Germany) and purified by tRNA (Sigma, Germany) precipitation. Tissue sections (5 sections per slide) were saturated with 100 µl of hybridization buffer (deionized formamide, 20×SSC, dextran sulfate, 0.2 M Na-phosphate buffer, Dehnhard's solution, 20% sarcosyl, salmon sperm DNA, 5 M DTT (DL-Dithiothreitol)) containing 106 cpm/100 µl for the Avp-35S and 0.5×106 cpm/100 µl for the Oxt-35S oligoprobe. Coverslipped sections were incubated in a moist chamber for 18–22 h at 45°C. After several washes in 1×SSC, the sections were dehydrated and air-dried [20]. For all following steps see section of data analysis below.
AVP immunohistochemistry in the PVN.
Immunohistochemistry was performed as described previously [37]. In brief, animals were deeply anaethesized by an overdose of sodium pentobarbital (200 mg/kg) by i.p. injection prior to cardiac perfusion. Brains of these unstressed mice were removed and postfixed at 4°C in 4% paraformaldehyd over night and then stored in 0.2 M phosphate buffer until sectioning. Immunofluorescence staining was used to visualize AVP. Vibratome (Series 1000, Ted Pella Inc., Redding, CA) sections (50 µm) were preincubated in immunobuffer containing 5% normal donkey serum (Jackson ImmunoResearch, West Grove, PA) and then incubated for 48 h at room temperature with a polyclonal rabbit antibody to AVP (dilution 1∶10000, Peninsula Laboratories, Belmont, CA). Sections were exposed for 2.5 h to a Cy2-conjugated donkey anti-rabbit IgG (dilution 1∶100, Jackson ImmunoResearch). Negative controls were performed by omission of the primary antibodies. No staining was observed in any of the control sections. Images were recorded on a fluorescence microscope (Olympus BX-51, Vienna, Austria), equipped with a U-M41001 filter system (excitation filter, 455–495 nm; dichromatic mirror, 505 nm; emission filter, 510–555 nm) for Cy2 (excitation maximum 492 nm; emission maximum 510 nm), using a digital camera (Olympus DP50) and analySIS® image processing software (Soft Imaging Systems, Münster, Germany).
AVP receptor 1a autoradiography.
For AVP receptor 1a (V1a) autoradiography, frozen cryostat brain sections (14 µm) of HAB and LAB mice were utilized and processed as described elsewhere [20]. Before incubation, the slides were fixed using 0.1% PFA and washed in Tris buffer. For competitive displacement of endogenous AVP and for receptor binding, slides were incubated (60 min, RT) with the AVP V1a-R antagonist 125I-lin AVP (NEX310, NEN DuPont, Boston, MA; 200 cpm/µl). After several washes, the slides were air-dried.
Data analyses.
For all ISHs and receptor binding analyses, the sections were exposed to a Kodak BioMax film (Amersham, Braunschweig, Germany) for 1–3 days (increased exposure time for Oxt due to the lower amount of labelled oligoprobe). Autoradiographs were quantified using an image analysis optical software (Optimas 5.2, Optimas Corporation, Seattle, WA). Autoradiograms from ISH and receptor binding studies were analyzed using the optical density readings (grey intensity). Three to six sections were analyzed and the highest expression (measured as hybridization signal minus background) per animal was taken for calculating mean hybridization per nucleus and line [20]. Due to the high expression of Avp and Oxt mRNAs in the hypothalamic PVN and SON, we focused on these nuclei in LAB, HAB, NAB or F1 animals. For V1a receptor autoradiography, the amygdala, bed nucleus of the stria terminalis, hippocampus, lateral septum, PVN, SON, suprachiasmatic nucleus and ventral pallidum were analyzed.
Avp sequencing
Genomic DNA was purified from tail tissue clips of HAB and LAB mice from the 20th generation, using the NucleoSpin Tissue kit (Macherey-Nagel, Düren, Germany). The Avp coding sequence [25] as well as 3 kbp up- and downstream the coding region were amplified using 12 primer pairs (MWG Biotech, Ebersberg, Germany; Table 2). All PCR amplifications were carried out by using Taq polymerase (Fermentas, St. Leon-Rot, Germany) in 25-µl reactions under the following conditions: initial denaturation at 94°C for 4 min, 40 cycles of denaturation (94°C for 1 min), annealing (52–56°C for 1 min) and extension (72°C for 1 min), concluded by a final extension of 10 min at 72°C. Purification of the PCR products was performed with the NucleoFast kit (Macherey-Nagel). Then, the PCR products were ligated into the pGEM-T vector (Promega Corp., Madison, WI), according to the manufacturer's instructions and transformed into competent E. coli DH5α cells via electroporation. Transformed bacteria were collected in 1 ml SOB medium, 150 µl of the suspension were incubated overnight at 37°C on lysogeny broth/ampicillin agar plates inoculated with 5-bromo-4-chloro-3-indolyl-β-D-galactoside (X-Gal: 25 µl; 10%) and isopropyl-β-D-thiogalactopyranoside (IPTG: 30 µl; 0.1 M) solved in dimethyl sulfoxide. Colonies were picked and their plasmid inserts were amplified using T7 and SP6 primers in a 35- cycles PCR with an annealing temperature of 49°C and otherwise identical conditions, as described above. 5 µl of each PCR product were analyzed on a 1.5% agarose gel, another 2 µl were purified with the NucleoFast 96 PCR kit (Macherey-Nagel) and used for the sequencing reaction (BigDye Terminator Kit, Applied Biosystems, Foster City, CA) that was performed according to the manufacturer's instructions. Sequencing reaction products were purified with the Sequencing Reaction Cleanup Kit (Millipore, Billerica, MA) and resolved by electrophoresis on a 3730 DNA Analyzer (Applied Biosystems). DNA sequences around the identified polymorphic sites were analyzed via the Transcription Element Search System (TESS, http://www.cbil.upenn.edu/cgi-bin/tess/tess). Putative binding sites of transcription factors were additionally checked for their function and occurrence within the brain according to the NCBI database (http://www.ncbi.nlm.nih.gov/).
[Figure omitted. See PDF.]
Table 2. Forward and reverse primers used for Avp sequencing and their specific annealing temperatures.
https://doi.org/10.1371/journal.pone.0005129.t002
Allele-specific transcription analysis
PVN and SON were dissected from the hypothalamus of four female F1 mice that were heterozygous for the recently described C(40)T (rs50049109) single-nucleotide polymorphism (SNP) [25]. Total RNA was extracted using a Trizol (Invitrogen) chloroform extraction protocol and reverse transcribed with Omniscript (Qiagen GmbH, Hilden, Germany). 208 bp of the Avp transcript containing the SNP at position C(40)T were amplified by PCR in 40 identical cycles of denaturation at 94°C for 1 min, annealing at 56°C for 1 min and elongation at 72°C for 1 min using Taq polymerase (Fermentas) and 5′ GAG CAG AGC CTG AGC TGC ACA CAG T 3′ as forward and 5′ AGC AGA TGC TTG GTC CGA AGC ACG 3′ as reverse primers (MWG Biotech). Amplified fragments were ligated into a vector and transfected into competent E. coli cells that were grown overnight as described above.
From each agar plate, 25 colonies containing insert DNA (according to white/blue selection) were picked and sequenced. The specific alleles from each transcript were determined by the C(40)T polymorphism.
Avp promoter deletion and signal peptide polymorphism in a freely-segregating F2 panel
Biallelic DNA sequence polymorphisms were assessed using a custom-designed oligonucleotide pool and the GoldenGate assay on an Illumina Sentrix Array Matrix platform (Illumina, San Diego, CA). In the parental generation, the two polymorphisms located in the coding sequence or the promoter of Avp display opposite homozygosity in HAB (+/+ and CC) vs. LAB (−/− and TT) animals. A total of 508 F2 offspring of a HABxLAB F1 intercross progeny were phenotyped and genotyped at these two loci based on DNA isolated from tail tips as described recently [25]. Genotyping was performed using the manufacturer's protocol and genotypes were determined using the Illumina BeadStudio software 3.1.0.0 and the genotyping module 3.1.12. Twenty four samples were genotyped twice, assuring data quality and reproducibility.
Statistical analyses
The presented data are given as means+SEM and were statistically analyzed using SPSS 12.0. Two independent groups were compared using the Mann-Whitney U-test. For three independent groups, the nonparametric Kruskal-Wallis H-test was used with post hoc Mann-Whitney U-test comparisons, including sequential Bonferroni correction of the respective p-values, if applicable. For correlation analysis, a Spearman-Rho correlation analysis was performed. For allele-specific transcription analysis, deviation from the expected value was estimated using the chi-square test.
For the F2 panel, association analysis was performed using a permutation approach (WG-Permer, www.mpipsykl.mpg.de/wg-permer). Trait values were rank-transformed to protect against artifacts. Analysis was done using a genotypic model, i.e. the three genotypic classes possible for each of the phenotypes were treated as a separate class each and global test on equality of the three means was performed.
Results
Behavioral phenotype and Avp expression in LAB, HAB and NAB mice
EPM.
In accordance with our breeding strategy, LAB, HAB and NAB male mice showed significant differences in the percent of total test time spent on the open arms of the EPM (Fig. 1A), with LABs spending most time (about 60%) on the open arms, HABs least time (about 5%) and NAB animals displaying an intermediate phenotype (about 30%; p<0.01 for all group comparisons).
[Figure omitted. See PDF.]
Figure 1. Alterations in behavior and vasopressin (Avp) expression.
Phenotypic measures of HAB/NAB/LAB mice reflecting behavior (A) on the elevated plus-maze (EPM), in the (B) tail-suspension test (TST) and (C) forced swim test (FST) and corresponding vasopressin (Avp) expression patterns as measured by (D) in situ hybridization (ISH) in the paraventricular nucleus (PVN) in male mice and by (E) quantitative PCR (qPCR) in the PVN of male and female mice. Data are shown as means+SEM; * p<0.05; ** p<0.01; *** p<0.001; Mann-Whitney-U test.
https://doi.org/10.1371/journal.pone.0005129.g001
Tail suspension and forced swim test.
Male LAB animals exhibited a dramatically decreased immobility in the tail suspension test (Fig. 1B) (p<0.01 vs. HAB and NAB). Similarly, LAB mice showed decreased floating time in the forced swim test relative to HAB and NAB mice (p<0.01 each; Fig. 1C).
Gene and protein expression analysis and validation
Analysis of genome wide expression revealed a significantly decreased Avp mRNA level in the PVN (p<0.05) and showed a tendency towards the same direction in the SON (p<0.1) of LAB compared to HAB mice. Importantly, no differences were observed for Avpr1a, Avpr1b, Oxt or Oxtr (p>0.1) in the respective brain regions investigated (Table 3).
[Figure omitted. See PDF.]
Table 3. Regulation of gene products of the vasopressin and oxytocin systems from the MPI24K platform-based microarray experiment of HAB vs. LAB mice.
https://doi.org/10.1371/journal.pone.0005129.t003
Confirming microarray data, ISH revealed a decreased expression of Avp in both the PVN (Fig. 1D) and SON (Fig. 2A) in male LAB compared to both HAB and NAB mice (PVN: p<0.01 vs. HAB and NAB; SON: p<0.05 vs. HAB and NAB). Similarly, qPCR analysis showed a reduced Avp expression in the PVN (Fig. 1E) and SON (Fig. 2B) of both male and female LABs compared to HAB or NAB animals (PVN of male mice: p<0.05 vs. HAB; SON of male mice: p<0.05 vs. HAB; PVN of female mice: p<0.01 vs. HAB; SON of female mice: p<0.05 vs. HAB and NAB). To address the question whether the described differences in AVP transcript levels actually translate into biologically relevant AVP peptide contents, we carried out immunohistochemistry with an antibody directed against AVP. Confirming ISH data, these experiments showed a decreased amount of AVP peptide in LAB compared to HAB mice in the PVN (Fig. 3).
[Figure omitted. See PDF.]
Figure 2. Gene expression of HAB, NAB, LAB, and HABxLAB F1 intercross mice.
Vasopressin (Avp) mRNA expression in the supraoptic nucleus (SON) measured by (A) in situ hybridization (ISH) and (B) quantitative PCR (qPCR) under basal conditions in male and female mice. Data are shown as means+SEM; T p<0.1; * p<0.05; ** p<0.01; Mann-Whitney-U test.
https://doi.org/10.1371/journal.pone.0005129.g002
[Figure omitted. See PDF.]
Figure 3. Vasopressin (AVP) immunohistochemistry analysis.
AVP fluorescence antibody-staining of the paraventricular nucleus (PVN) flanking the 3rd ventricle (3V) in HAB and LAB mice. For corresponding in situ hybridization, see Fig. 7.
https://doi.org/10.1371/journal.pone.0005129.g003
In contrast to AVP, the structurally related hypothalamic neuropeptide OXT failed to reveal any significant differences among the lines (Fig. 4) as detected by ISH. Also, no differences in the expression of Avpr1b or Oxtr (p>0.1) could be confirmed using qPCR (Table 3). Likewise, V1a autoradiography failed to reveal any difference between the HAB and LAB lines in a variety of brain regions critically involved in the regulation of anxiety- and depression-related behaviors, thus further strengthening the concept that events upstream of the receptor contribute to the behavioral phenotype of HAB and LAB mice (Table 3).
[Figure omitted. See PDF.]
Figure 4. Oxytocin (Oxt) mRNA expression.
Oxt in the paraventricular nucleus (PVN) measured by in situ hybridization under basal conditions in male HAB, NAB and LAB mice. Data are shown as means+SEM.
https://doi.org/10.1371/journal.pone.0005129.g004
Polymorphic AVP alleles
Next we investigated whether differences in AVP expression originate from differences in the sequence of HAB and LAB alleles that could underpin differences in the recruitment of proteins driving transcription. Sequencing of the Avp gene resulted in the identification of nine polymorphic sites, of which eight correspond to SNPs and one comprises a 12 bp deletion (Fig. 5).
[Figure omitted. See PDF.]
Figure 5. Vasopressin (Avp) gene sequence of HAB vs. LAB mice.
Polymorphic sites are indicated with positions from transcription start (−1 to −2600 bp) in (A) the promoter region (two SNPs and deletion in LABs); (B) within the Avp coding sequence from transcription start (1 to 1960 bp; three SNPs); (C) in the intergenic region between Avp and oxytocin (Oxt) from the end of the last exon (+1 to +2600 bp; three SNPs). Exons and untranslated regions (UTRs) are indicated by boxes (exons shaded, UTRs completely filled black or white).
https://doi.org/10.1371/journal.pone.0005129.g005
The upstream promoter region contains two SNPs: T(−2521)C with LABs carrying cytosine and HABs thymine and C(−1422)T with LABs carrying thymine and HABs cytosine. Furthermore, LAB mice miss a 12 bp segment mapping to Δ(−2180–2191). The gene-coding locus contains three SNPs: C(40)T with LABs carrying thymine and HABs cytosine, A(1431)G with LABs carrying guanine and HABs adenine and T(1527)C with LABs carrying cytosine and HABs thymine (see also [25]). The downstream enhancer region, also referred to as the intergenic region between Avp and Oxt, contains three SNPs: A(+399)G with LABs carrying guanine and HABs adenine, T(+476)G with LABs carrying guanine and HABs thymine and C(+2444)A with LABs carrying adenine and HABs cytosine. A detailed in silico search for potential transcription factor binding sites at the polymorphic loci led to the identification of a number of candidates that are summarized in Table 4.
[Figure omitted. See PDF.]
Table 4. Polymorphisms in the promoter and downstream enhancer region between the HAB and LAB specific sequence with probable binding factors.
https://doi.org/10.1371/journal.pone.0005129.t004
In all our studies, involving LAB, HAB, NAB and F2 mice, the C(40)T SNP and 12 bp deletion were in complete linkage, i.e. the homozygous LAB-specific combination of alleles (TT) at the 40th bp of the transcript-coding sequence was strictly associated with the absence of the 12 bps of the promoter. In analogy, if the transcript-coding alleles were CT or CC, the −2180–2191 locus bore the +/− and +/+ allelic constellation, respectively. The most frequent genotypic combination in NAB mice was the HAB-specific one (Table 5).
[Figure omitted. See PDF.]
Table 5. Allele frequency in a NAB population (n = 165) for the SNP in the Avp signal peptide and the strictly linked promoter deletion.
https://doi.org/10.1371/journal.pone.0005129.t005
Allele-specific Avp transcription in heterozygous HABxLAB F1 intercross mice
To obtain evidence for a functional role of these polymorphisms in contributing to different levels of AVP expression in male and female LAB vs. HAB mice, we studied the activity of each allele in an F1 intercross. Under this condition, both alleles are contained in the same cellular background, thus eliminating differential synaptic input as a confounding factor. In support of an allele-dependent effect, heterozygous F1 mice revealed intermediate levels of Avp mRNA expression in the SON in comparison to LABs and HABs (p<0.1 F1 vs. HAB and LAB, p<0.05 LAB vs. HAB; Fig. 2A). In line with these findings, the LAB-specific allele displayed a significantly decreased expression relative to the HAB-specific one in both the PVN and SON (p<0.001; Fig. 6), consistent with Avp mRNA expression data of HAB vs. LAB mice (see Figs. 1D, E; 2A, B).
[Figure omitted. See PDF.]
Figure 6. Allele-specific vasopressin (Avp) expression.
Proportion of HAB allele-specific vs. LAB allele-specific Avp transcripts from heterozygous F1 (HABxLAB intercross) mice in the paraventricular nucleus (PVN; χ2 = 14.4 p = 1.4e-4) and the supraoptic nucleus (SON; χ2 = 15.2 p = 4.8e-5). Data are shown as means+SEM; *** p<0.001.
https://doi.org/10.1371/journal.pone.0005129.g006
Correlation of Avp mRNA in the PVN and anxiety/depression-like behaviors of HAB, F1 and LAB mice
In light of the above data pointing to a dual role of hypothalamic AVP in contributing to HAB and LAB phenotypes, respectively, we performed correlation analyses with the aim to obtain a more detailed picture of AVP's role in anxiety- and depression-like behaviors. Avp mRNA expression in the PVN of male LAB, F1 and HAB mice under basal conditions is significantly correlated with anxiety-related behavior (r = 0.619; p<0.001; Fig. 7A) as well as depression-like behavior (r = 0.655; p<0.001; Fig. 7B).
[Figure omitted. See PDF.]
Figure 7. Correlation between gene expression and behavior.
Correlation of vasopressin (Avp) mRNA expression in the paraventricular nucleus (PVN) with (A) anxiety-related behavior and (B) depression-like behavior of HAB, F1 and LAB male mice under basal conditions. For corresponding immunohistochemistry, see Fig. 3.
https://doi.org/10.1371/journal.pone.0005129.g007
Association between the phenotype and polymorphisms at the AVP locus
To further test their functional involvement, the Avp-related polymorphisms were analyzed in a freely-segregating F2 panel (n = 508). In a total of 258 male mice with HAB grandmothers, 49 animals were identified as carrying the LAB-typical polymorphisms, i.e. the 12 bp deletion and the C(40)T in the combination −/− and TT, 67 animals the HAB-typical +/+ and CC and 142 animals the heterozygous +/− and CT genotypes. In all F2 animals, the allele distribution roughly reflected a 1∶2∶1 ratio (Table 6). If these polymorphisms have an impact on anxiety-related behavior, then homozygous −/− and TT F2 animals should behave less anxious than their +/+ and CC counterparts. Indeed, a co-segregation could be detected with −/− and TT mice being significantly less anxious on the EPM (p<0.05; Fig. 8A). Importantly, this association of the genetic polymorphisms with anxiety levels was independent of locomotor activity, measured in the open field (Fig. 8B). In the remaining 250 F2 male mice with LAB grandmothers, no association was detectable (data not shown).
[Figure omitted. See PDF.]
Figure 8. Association between anxiety-related behavior and the Avp-related genotypes in the freely-segregating F2 panel.
Phenotypic indices of male F2 mice with HAB grandmothers and Avp-related genotypes (+/+) and CC (HAB-typical homozygotes, n = 67), (+/−) and CT (heterozygotes, n = 142), or (−/−) and TT (LAB-typical homozygotes, n = 49). (A) % time spent on the open arms of the elevated plus-maze (EPM) and (B) total distance traveled in the open field (OF). LAB-typical homozygotes were less anxious, independent of locomotor activity. Data are shown as means+SEM; * p<0.05.
https://doi.org/10.1371/journal.pone.0005129.g008
[Figure omitted. See PDF.]
Table 6. Allele frequency in a freely-segregating F2 panel (HABxLAB intercross offspring, n = 508) for the SNP in the Avp signal peptide and the strictly linked promoter deletion.
https://doi.org/10.1371/journal.pone.0005129.t006
Discussion
In this study, we identified polymorphisms in the Avp gene underlying a decreased expression of Avp mRNA in the hypothalamus of LAB animals under basal conditions compared to HAB, NAB and F1 mice, irrespective of gender. Both correlative and associative evidence indicates that the hypomorphic Avp allele causally contributes to the locomotion-independent reduction of anxiety-related behavior typical of LAB mice.
Importantly, differences in Avp expression were not accompanied by changes in V1a and V1b receptors that could confound the analysis of the behavioral involvement of AVP in these lines. Several groups independently reported an association of the V1a receptor with anxiety-related behavior. V1a receptor knockout mice, although at times providing inconsistent results, showed impaired social interaction, social recognition and reduced anxiety-related behavior [38]–[40]. Furthermore, V1a antisense RNA targeting the septum made rats less anxious [21]. Similarly, V1b antagonism was found to reduce anxiety- and depression-like behaviors [41], [42]. In the present study, expression profiling in combination with qPCR and receptor autoradiography revealed no V1a/b-related difference between HAB and LAB mice in a variety of anxiety-associated brain regions, suggesting that line-specific divergences in behavior are primarily due to events upstream of the receptor, i.e. differential expression and release of AVP.
Despite the high sequence homology of the evolutionary related and tail-to-tail oriented Avp and Oxt genes that are separated by 4 kbp, the so-called intergenic region [43], and despite the anxiolytic potency of the central OXT system [44], [45], differences in gene expression were only restricted to the former. This is consistent with the idea that our bidirectional breeding protocol selects for genetic variation in specific anxiety-related signaling pathways rather than disrupting hypothalamic functions in general.
To identify the molecular mechanisms underlying the reduced Avp expression, we sequenced the Avp-coding region and the surrounding structures. Nine polymorphic loci were identified reliably differing between HAB and LAB mice; all structural parts, the upstream promoter, the gene locus and the downstream enhancer region, contained three polymorphisms each. How may the polymorphisms in the promoter region translate into different expression profiles and phenotypes? While the C(−1422)T SNP is not located in a known transcription factor binding site, the T(−2521)C SNP and the Δ(−2180–2191) deletion are located directly in the center of a binding site for nuclear factor 1 (NF-1). This well-described transcription factor in the brain has the potential to enhance or impair gene expression, depending on the concomitance of a C/EBP binding site [46]–[49]. At the Δ(−2180–2191) deletion site, the binding of NF-1 would be impaired. The T(−2521)C SNP could enhance binding, but has no active concomitant C/EBP site, thus ultimately contributing to reduced gene expression. Thus, both polymorphic sites are likely to lead to an impaired Avp expression in mice bearing the LAB-specific DNA sequence.
While there were no polymorphisms identified in the non-coding (untranslated region and intronic) sequences, all three SNPs in the gene locus were located in the coding sequence. Two of these are silent mutations (A(1431)G and T(1527)C), and the third one (C(40)T) causes a substitution of alanine to valine in the third amino acid of the AVP signal peptide. As already reported by Keßler et al. [25], this genetic marker co-segregated in an F2 panel with symptoms of central diabetes insipidus in LAB mice and, partially, with anxiety-related behavior, further confirming the strict linkage between the promoter deletion Δ(−2180–2191) and C(40)T.
Also SNPs in the downstream enhancer region can have an impact on expression. There is a repetition of motifs from the 178-bp region between +370 bp and +480 bp [43], where two further SNPs have been identified. Analysis of this region resulted in the identification of a binding site for c-Ets-2 near a binding site of C/EBPbeta [48]. In the center of the c-Ets-2 binding site, the A(+399)G polymorphism would allow for an enhanced transcription rate in the HAB-specific but not the LAB-specific DNA sequence. Finally, even polymorphisms in the coding region are likely to contribute to lower Avp mRNA content in LAB mice by negatively influencing mRNA secondary structure and stability [50].
Still, we had to make sure that Avp under-expression in LAB mice is definitely due to the described genetic polymorphisms rather than to differences in other variables, including epistatic gene-gene interactions or synaptic input to the hypothalamic nuclei, likely to differ between the lines. Therefore, to eliminate the enormous complexity of the intact organism, we decided to use allele-specific transcription analysis, which had already been successfully applied in HAB/LAB rats [18]. By cross-mating HAB and LAB mice, we produced heterozygous F1 animals that host both HAB and LAB line-specific alleles in each cell. The analysis revealed a strongly decreased expression of Avp but not Oxt by the LAB-specific allele compared to the HAB-specific one, thus further confirming that line-specific polymorphisms are causally related to line-specific expression profiles of Avp and, consequently, the phenotype of trait anxiety.
However, while AVP expression is clearly under genetic control, the regulation of anxiety-related behavior is more complex. Compounding the difficulties associated with detecting multiple small or moderate genetic effects, the latter operate via subtle interactions with each other and the environment, thus further masking the functional significance of given polymorphisms as described here. Therefore, we used different approaches to examine the possible contribution of the hypomorphic allele-driven AVP deficit to trait anxiety, including correlational analyses in HAB/F1/LAB mice and association studies in F2 animals. Importantly, experimental crosses such as an F2 intercross are not subject to allelic association between genes on different chromosomes, thus eliminating spurious associations between markers and phenotypes [51]. Indeed, in addition to correlative evidence, a significant, though moderate association between the promoter deletion and the C(40)T SNP and locomotion-independent anxiety-related behavior was demonstrated in a freely-segregating F2 panel, indicative of a causal contribution to the phenotype. It is of note in this context that these F2 animals were assigned according to their Avp-related genotype, while all the other (anxiety-relevant) genes freely segregated, explaining why differences in anxiety-related behavior are rather modest (Fig. 8A) compared to those seen in Fig. 1A.
As mentioned before, the individual phenotype is known to be shaped by phenomena that are or are not mediated by sequence variations in DNA, i.e. inheritance is both genetic and environmental, with the latter primarily relying on the quality of mother-infant interaction [52]. To reduce possible influences of different epigenetic effects, deriving from the fact that half of the F2 had LAB and the other half HAB grandmothers, split calculations according to the grandmaternal background were performed for associations. Indeed, a significant association at the Avp locus was found for F2 mice with HAB grandmothers. Thus, while the polymorphisms identified in the Avp gene are unlikely to generally cause hypo-anxiety, their variation may contribute to the severity of the phenotype, further confirming the molecular basis of many small or moderate genetic effects giving rise to trait anxiety. In other words, while the link between genes, the brain and behavior remains weak, variation in the Avp locus may be one contributor, among many others, to trait anxiety. Several dozen quantitative trait loci, containing one gene or several genes, may be involved in anxiety regulation, and most of them are yet to be identified [53], [54].
Further studies will focus on the importance of the (grand)parental background for anxiety-related behavior, as genes regulated by (grand)parent-specific epigenetic modifications can lead to monoallelic gene expression [55], [56] or are simply X- or Y-linked. As there is a higher prevalence of anxiety and depression in women than in men [57], the (grand)parental influence on anxiety-related behavior, as revealed in our F2 panel, might also be due to a modulating effect of a sex-specific locus. HAB/LAB mice might thus represent a model to further investigate the complex pattern of genetic vs. epigenetic inheritance.
As LAB mice exhibit a deficit not only in AVP precursor processing and neuropeptide release from both dendrites and axon terminals [25], but also in Avp expression, additive effects at multiple levels are likely to produce a more or less “global” deficit in bioactive AVP. Remarkably, recent studies in Brattleboro rats and humans confirmed that an AVP deficit may be accompanied by symptoms of central diabetes insipidus, reduced anxiety-related or attenuated depression-like behavior [58] and signs of diminished agoraphobia [15], respectively, suggesting complementary inter-species genetics.
Determining the most common Avp-related genotypes in the NAB line revealed that the HAB-specific polymorphisms represent the most common genetic variant. More than 70% of the NAB mice were found to be homozygous for the HAB-specific sequence, whereas less than 1% carried the LAB-specific allele homozygously. Nevertheless, the respective distribution is in Hardy-Weinberg equilibrium, although with the LAB-specific allele at a decreased frequency. The fact that the NAB line's most common genotype for the Avp gene and surrounding sequence fully corresponds to the HAB-specific sequence helps to explain the often observed discrepancy between NAB and F1 mice, both of which are often used as controls. In contrast to strictly intermediate F1, NAB mice displayed in about 70% of all animals a HAB-like Avp expression and depression-like behavior in the tail-suspension test (Fig. 1B, D, E and [24]), which has also been demonstrated recently in an unselected CD1 population [59], thus further emphasizing the functional implications of centrally released AVP.
The application of similar breeding strategies in two different species (Rattus norvegicus and Mus musculus) resulted in similar findings concerning the identification of polymorphisms in and around the Avp gene and both the expression and release patterns of AVP. This was shown to be likely involved in mediating anxiety-related behavior. While these results cannot entirely exclude an effect of genetic drift, they strongly point to an evolutionary conserved dynamic role of AVP along a continuum, with hypomorphic Avp alleles finally contributing to low anxiety- and possibly depression-like behaviors and hypermorphic alleles contributing to hyper-anxiety. Together with human studies [23], [60], they do not only ensure the construct validity of the HAB/LAB model, but additionally provide a strong impetus for Avp expression as a potential biomarker and therapeutic target in psychopathology.
Acknowledgments
We would like to thank Alexandra Wigger, Wolfgang Jacob, Andrea Steiner, Claudia Kühne, Sabine Damast and Markus Nußbaumer for expert technical assistance and their contributions to this study.
Author Contributions
Conceived and designed the experiments: BM LC EF CM BMM TB DS RL. Performed the experiments: BM LC CH CG EF MSK CM PM. Analyzed the data: BM LC CH CG EF MSK MG PW BP TB. Contributed reagents/materials/analysis tools: EF MSK CM BMM MG PW BP PM MP NS TB JD FH DS RL. Wrote the paper: BM LC RL.
Citation: Bunck M, Czibere L, Horvath C, Graf C, Frank E, Keßler MS, et al. (2009) A Hypomorphic Vasopressin Allele Prevents Anxiety-Related Behavior. PLoS ONE4(4): e5129. https://doi.org/10.1371/journal.pone.0005129
1. Nestler EJ, Barrot M, DiLeone RJ, Eisch AJ, Gold SJ, et al. (2002) Neurobiology of depression. Neuron 34(1): 13–25.EJ NestlerM. BarrotRJ DiLeoneAJ EischSJ Gold2002Neurobiology of depression.Neuron34(1)1325
2. Ustun TB, Ayuso-Mateos JL, Chatterji S, Mathers C, Murray CJ (2004) Global burden of depressive disorders in the year 2000. Br J Psychiatry 184: 386–92.TB UstunJL Ayuso-MateosS. ChatterjiC. MathersCJ Murray2004Global burden of depressive disorders in the year 2000.Br J Psychiatry18438692
3. Kessler RC, McGonagle KA, Zhao S, Nelson CB, Hughes M, et al. (1994) Lifetime and 12-month prevalence of DSM-III-R psychiatric disorders in the United States. Results from the National Comorbidity Survey. Arch Gen Psychiatry 51(1): 8–19.RC KesslerKA McGonagleS. ZhaoCB NelsonM. Hughes1994Lifetime and 12-month prevalence of DSM-III-R psychiatric disorders in the United States. Results from the National Comorbidity Survey.Arch Gen Psychiatry51(1)819
4. Fava M, Kendler KS (2000) Major depressive disorder. Neuron 28(2): 335–41.M. FavaKS Kendler2000Major depressive disorder.Neuron28(2)33541
5. Sanders AR, Detera-Wadleigh SD, Gershon ES (1999) Molecular genetics of mood disorders. In: Charney DS, Nestler EJ, Bunney BS, editors. Neurobiology of mental illness. 1st ed. New York: Oxford University Press. pp. 299–316.AR SandersSD Detera-WadleighES Gershon1999Molecular genetics of mood disorders.DS CharneyEJ NestlerBS BunneyNeurobiology of mental illness. 1st edNew YorkOxford University Press299316
6. Ressler KJ, Nemeroff CB (2000) Role of serotonergic and noradrenergic systems in the pathophysiology of depression and anxiety disorders. Depress Anxiety 12: Suppl 12–19.KJ ResslerCB Nemeroff2000Role of serotonergic and noradrenergic systems in the pathophysiology of depression and anxiety disorders.Depress Anxiety12Suppl 1219
7. Landgraf R (2001) Neuropeptides and anxiety-related behavior. Endocr J 48(5): 517–33.R. Landgraf2001Neuropeptides and anxiety-related behavior.Endocr J48(5)51733
8. Gillespie NA, Kirk KM, Evans DM, Heath AC, Hickie IB, et al. (2004) Do the genetic or environmental determinants of anxiety and depression change with age? A longitudinal study of Australian twins. Twin Res 7(1): 39–53.NA GillespieKM KirkDM EvansAC HeathIB Hickie2004Do the genetic or environmental determinants of anxiety and depression change with age? A longitudinal study of Australian twins.Twin Res7(1)3953
9. Kas MJ, Fernandes C, Schalkwyk LC, Collier DA (2007) Genetics of behavioural domains across the neuropsychiatric spectrum; of mice and men. Mol Psychiatry 12(4): 324–30.MJ KasC. FernandesLC SchalkwykDA Collier2007Genetics of behavioural domains across the neuropsychiatric spectrum; of mice and men.Mol Psychiatry12(4)32430
10. Ring RH (2005) The central vasopressinergic system: examining the opportunities for psychiatric drug development. Curr Pharm Des 11(2): 205–25.RH Ring2005The central vasopressinergic system: examining the opportunities for psychiatric drug development.Curr Pharm Des11(2)20525
11. Landgraf R (2006) The involvement of the vasopressin system in stress-related disorders. CNS Neurol Disord Drug Targets 5(2): 167–79.R. Landgraf2006The involvement of the vasopressin system in stress-related disorders.CNS Neurol Disord Drug Targets5(2)16779
12. Frank E, Landgraf R (2008) The vasopressin system - from antidiuresis to psychopathology. Eur J Pharmacol 583(2–3): 226–42.E. FrankR. Landgraf2008The vasopressin system - from antidiuresis to psychopathology.Eur J Pharmacol583(2–3)22642
13. Birkett SD, Pickering BT (1993) The aberrant vasopressin precursor in rats with central diabetes insipidus (Brattleboro rats). Response to stimuli. Ann N Y Acad Sci 689: 83–90.SD BirkettBT Pickering1993The aberrant vasopressin precursor in rats with central diabetes insipidus (Brattleboro rats). Response to stimuli.Ann N Y Acad Sci6898390
14. Bohus B, de Wied D (1998) The vasopressin deficient Brattleboro rats: a natural knockout model used in the search for CNS effects of vasopressin. Prog Brain Res 119: 555–73.B. BohusD. de Wied1998The vasopressin deficient Brattleboro rats: a natural knockout model used in the search for CNS effects of vasopressin.Prog Brain Res11955573
15. Bruins J, Kovacs GL, Abbes AP, Burbach JP, van den Akker EL, et al. (2006) Minor disturbances in central nervous system function in familial neurohypophysial diabetes insipidus. Psychoneuroendocrinology 31(1): 80–91.J. BruinsGL KovacsAP AbbesJP BurbachEL van den Akker2006Minor disturbances in central nervous system function in familial neurohypophysial diabetes insipidus.Psychoneuroendocrinology31(1)8091
16. Verbalis JG (2003) Diabetes insipidus. Rev Endocr Metab Disord 4(2): 177–85.JG Verbalis2003Diabetes insipidus.Rev Endocr Metab Disord4(2)17785
17. Landgraf R, Kessler MS, Bunck M, Murgatroyd C, Spengler D, et al. (2007) Candidate genes of anxiety-related behavior in HAB/LAB rats and mice: focus on vasopressin and glyoxalase-I. Neurosci Biobehav Rev 31(1): 89–102.R. LandgrafMS KesslerM. BunckC. MurgatroydD. Spengler2007Candidate genes of anxiety-related behavior in HAB/LAB rats and mice: focus on vasopressin and glyoxalase-I.Neurosci Biobehav Rev31(1)89102
18. Murgatroyd C, Wigger A, Frank E, Singewald N, Bunck M, et al. (2004) Impaired repression at a vasopressin promoter polymorphism underlies overexpression of vasopressin in a rat model of trait anxiety. J Neurosci 24(35): 7762–70.C. MurgatroydA. WiggerE. FrankN. SingewaldM. Bunck2004Impaired repression at a vasopressin promoter polymorphism underlies overexpression of vasopressin in a rat model of trait anxiety.J Neurosci24(35)776270
19. Keck ME, Welt T, Muller MB, Uhr M, Ohl F, et al. (2003) Reduction of hypothalamic vasopressinergic hyperdrive contributes to clinically relevant behavioral and neuroendocrine effects of chronic paroxetine treatment in a psychopathological rat model. Neuropsychopharmacology 28(2): 235–43.ME KeckT. WeltMB MullerM. UhrF. Ohl2003Reduction of hypothalamic vasopressinergic hyperdrive contributes to clinically relevant behavioral and neuroendocrine effects of chronic paroxetine treatment in a psychopathological rat model.Neuropsychopharmacology28(2)23543
20. Wigger A, Sanchez MM, Mathys KC, Ebner K, Frank E, et al. (2004) Alterations in central neuropeptide expression, release, and receptor binding in rats bred for high anxiety: critical role of vasopressin. Neuropsychopharmacology 29(1): 1–14.A. WiggerMM SanchezKC MathysK. EbnerE. Frank2004Alterations in central neuropeptide expression, release, and receptor binding in rats bred for high anxiety: critical role of vasopressin.Neuropsychopharmacology29(1)114
21. Landgraf R, Gerstberger R, Montkowski A, Probst JC, Wotjak CT, et al. (1995) V1 vasopressin receptor antisense oligodeoxynucleotide into septum reduces vasopressin binding, social discrimination abilities, and anxiety-related behavior in rats. J Neurosci 15(6): 4250–8.R. LandgrafR. GerstbergerA. MontkowskiJC ProbstCT Wotjak1995V1 vasopressin receptor antisense oligodeoxynucleotide into septum reduces vasopressin binding, social discrimination abilities, and anxiety-related behavior in rats.J Neurosci15(6)42508
22. Young LJ (2001) Oxytocin and vasopressin as candidate genes for psychiatric disorders: lessons from animal models. Am J Med Genet 105(1): 53–4.LJ Young2001Oxytocin and vasopressin as candidate genes for psychiatric disorders: lessons from animal models.Am J Med Genet105(1)534
23. Scott LV, Dinan TG (2002) Vasopressin as a target for antidepressant development: an assessment of the available evidence. J Affect Disord 72(2): 113–24.LV ScottTG Dinan2002Vasopressin as a target for antidepressant development: an assessment of the available evidence.J Affect Disord72(2)11324
24. Kromer SA, Kessler MS, Milfay D, Birg IN, Bunck M, et al. (2005) Identification of glyoxalase-I as a protein marker in a mouse model of extremes in trait anxiety. J Neurosci 25(17): 4375–84.SA KromerMS KesslerD. MilfayIN BirgM. Bunck2005Identification of glyoxalase-I as a protein marker in a mouse model of extremes in trait anxiety.J Neurosci25(17)437584
25. Kessler MS, Murgatroyd C, Bunck M, Czibere L, Frank E, et al. (2007) Diabetes insipidus and, partially, low anxiety-related behaviour are linked to a SNP-associated vasopressin deficit in LAB mice. Eur J Neurosci 26(10): 2857–64.MS KesslerC. MurgatroydM. BunckL. CzibereE. Frank2007Diabetes insipidus and, partially, low anxiety-related behaviour are linked to a SNP-associated vasopressin deficit in LAB mice.Eur J Neurosci26(10)285764
26. Lister RG (1987) The use of a plus-maze to measure anxiety in the mouse. Psychopharmacology (Berl) 92(2): 180–5.RG Lister1987The use of a plus-maze to measure anxiety in the mouse.Psychopharmacology (Berl)92(2)1805
27. Steru L, Chermat R, Thierry B, Simon P (1985) The tail suspension test: a new method for screening antidepressants in mice. Psychopharmacology (Berl) 85(3): 367–70.L. SteruR. ChermatB. ThierryP. Simon1985The tail suspension test: a new method for screening antidepressants in mice.Psychopharmacology (Berl)85(3)36770
28. Porsolt RD, Bertin A, Jalfre M (1978) “Behavioural despair” in rats and mice: strain differences and the effects of imipramine. Eur J Pharmacol 51(3): 291–4.RD PorsoltA. BertinM. Jalfre1978“Behavioural despair” in rats and mice: strain differences and the effects of imipramine.Eur J Pharmacol51(3)2914
29. Palkovits M (1973) Isolated removal of hypothalamic or other brain nuclei of the rat. Brain Res 59: 449–50.M. Palkovits1973Isolated removal of hypothalamic or other brain nuclei of the rat.Brain Res5944950
30. Paxinos G, Franklin K (2001) The mouse brain in stereotaxic coordinates. 2nd Edition ed. San Diego, San Francisco, New York, Boston, London, Sydney, Tokyo: Academic Press. G. PaxinosK. Franklin2001The mouse brain in stereotaxic coordinates. 2nd Edition edSan Diego, San Francisco, New York, Boston, London, Sydney, TokyoAcademic Press
31. Deussing JM, Kuhne C, Putz B, Panhuysen M, Breu J, et al. (2007) Expression profiling identifies the CRH/CRH-R1 system as a modulator of neurovascular gene activity. J Cereb Blood Flow Metab 27(8): 1476–95.JM DeussingC. KuhneB. PutzM. PanhuysenJ. Breu2007Expression profiling identifies the CRH/CRH-R1 system as a modulator of neurovascular gene activity.J Cereb Blood Flow Metab27(8)147695
32. Tsolakidou A, Trumbach D, Panhuysen M, Putz B, Deussing J, et al. (2008) Acute stress regulation of neuroplasticity genes in mouse hippocampus CA3 area - Possible novel signalling pathways. Mol Cell Neurosci 38(3): 444–52.A. TsolakidouD. TrumbachM. PanhuysenB. PutzJ. Deussing2008Acute stress regulation of neuroplasticity genes in mouse hippocampus CA3 area - Possible novel signalling pathways.Mol Cell Neurosci38(3)44452
33. Livak KJ, Schmittgen TD (2001) Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) Method. Methods 25(4): 402–8.KJ LivakTD Schmittgen2001Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) Method.Methods25(4)4028
34. Muller MB, Keck ME, Binder EB, Kresse AE, Hagemeyer TP, et al. (2003) ABCB1 (MDR1)-type P-glycoproteins at the blood-brain barrier modulate the activity of the hypothalamic-pituitary-adrenocortical system: implications for affective disorder. Neuropsychopharmacology 28(11): 1991–9.MB MullerME KeckEB BinderAE KresseTP Hagemeyer2003ABCB1 (MDR1)-type P-glycoproteins at the blood-brain barrier modulate the activity of the hypothalamic-pituitary-adrenocortical system: implications for affective disorder.Neuropsychopharmacology28(11)19919
35. Ivell R, Richter D (1984) Structure and comparison of the oxytocin and vasopressin genes from rat. Proc Natl Acad Sci U S A 81(7): 2006–10.R. IvellD. Richter1984Structure and comparison of the oxytocin and vasopressin genes from rat.Proc Natl Acad Sci U S A81(7)200610
36. Villar MJ, Ceccatelli S, Ronnqvist M, Hokfelt T (1994) Nitric oxide synthase increases in hypothalamic magnocellular neurons after salt loading in the rat. An immunohistochemical and in situ hybridization study. Brain Res 644(2): 273–81.MJ VillarS. CeccatelliM. RonnqvistT. Hokfelt1994Nitric oxide synthase increases in hypothalamic magnocellular neurons after salt loading in the rat. An immunohistochemical and in situ hybridization study.Brain Res644(2)27381
37. Hetzenauer A, Sinnegger-Brauns MJ, Striessnig J, Singewald N (2006) Brain activation pattern induced by stimulation of L-type Ca2+-channels: contribution of Ca(V)1.3 and Ca(V)1.2 isoforms. Neuroscience 139(3): 1005–15.A. HetzenauerMJ Sinnegger-BraunsJ. StriessnigN. Singewald2006Brain activation pattern induced by stimulation of L-type Ca2+-channels: contribution of Ca(V)1.3 and Ca(V)1.2 isoforms.Neuroscience139(3)100515
38. Bielsky IF, Hu SB, Ren X, Terwilliger EF, Young LJ (2005) The V1a vasopressin receptor is necessary and sufficient for normal social recognition: a gene replacement study. Neuron 47(4): 503–13.IF BielskySB HuX. RenEF TerwilligerLJ Young2005The V1a vasopressin receptor is necessary and sufficient for normal social recognition: a gene replacement study.Neuron47(4)50313
39. Bielsky IF, Hu SB, Szegda KL, Westphal H, Young LJ (2004) Profound impairment in social recognition and reduction in anxiety-like behavior in vasopressin V1a receptor knockout mice. Neuropsychopharmacology 29(3): 483–93.IF BielskySB HuKL SzegdaH. WestphalLJ Young2004Profound impairment in social recognition and reduction in anxiety-like behavior in vasopressin V1a receptor knockout mice.Neuropsychopharmacology29(3)48393
40. Egashira N, Tanoue A, Matsuda T, Koushi E, Harada S, et al. (2007) Impaired social interaction and reduced anxiety-related behavior in vasopressin V1a receptor knockout mice. Behav Brain Res 178(1): 123–7.N. EgashiraA. TanoueT. MatsudaE. KoushiS. Harada2007Impaired social interaction and reduced anxiety-related behavior in vasopressin V1a receptor knockout mice.Behav Brain Res178(1)1237
41. Griebel G, Simiand J, Serradeil-Le Gal C, Wagnon J, Pascal M, et al. (2002) Anxiolytic- and antidepressant-like effects of the non-peptide vasopressin V1b receptor antagonist, SSR149415, suggest an innovative approach for the treatment of stress-related disorders. Proc Natl Acad Sci U S A 99(9): 6370–5.G. GriebelJ. SimiandC. Serradeil-Le GalJ. WagnonM. Pascal2002Anxiolytic- and antidepressant-like effects of the non-peptide vasopressin V1b receptor antagonist, SSR149415, suggest an innovative approach for the treatment of stress-related disorders.Proc Natl Acad Sci U S A99(9)63705
42. Engin E, Treit D (2008) Dissociation of the anxiolytic-like effects of Avpr1a and Avpr1b receptor antagonists in the dorsal and ventral hippocampus. Neuropeptides 42(4): 411–21.E. EnginD. Treit2008Dissociation of the anxiolytic-like effects of Avpr1a and Avpr1b receptor antagonists in the dorsal and ventral hippocampus.Neuropeptides42(4)41121
43. Fields RL, House SB, Gainer H (2003) Regulatory domains in the intergenic region of the oxytocin and vasopressin genes that control their hypothalamus-specific expression in vitro. J Neurosci 23(21): 7801–9.RL FieldsSB HouseH. Gainer2003Regulatory domains in the intergenic region of the oxytocin and vasopressin genes that control their hypothalamus-specific expression in vitro.J Neurosci23(21)78019
44. Veenema AH, Neumann ID (2008) Central vasopressin and oxytocin release: regulation of complex social behaviours. Prog Brain Res 170: 261–76.AH VeenemaID Neumann2008Central vasopressin and oxytocin release: regulation of complex social behaviours.Prog Brain Res17026176
45. Waldherr M, Neumann ID (2007) Centrally released oxytocin mediates mating-induced anxiolysis in male rats. Proc Natl Acad Sci U S A 104(42): 16681–4.M. WaldherrID Neumann2007Centrally released oxytocin mediates mating-induced anxiolysis in male rats.Proc Natl Acad Sci U S A104(42)166814
46. Ji C, Chen Y, Centrella M, McCarthy TL (1999) Activation of the insulin-like growth factor-binding protein-5 promoter in osteoblasts by cooperative E box, CCAAT enhancer-binding protein, and nuclear factor-1 deoxyribonucleic acid-binding sequences. Endocrinology 140(10): 4564–72.C. JiY. ChenM. CentrellaTL McCarthy1999Activation of the insulin-like growth factor-binding protein-5 promoter in osteoblasts by cooperative E box, CCAAT enhancer-binding protein, and nuclear factor-1 deoxyribonucleic acid-binding sequences.Endocrinology140(10)456472
47. Alikhani-Koupaei R, Fouladkou F, Fustier P, Cenni B, Sharma AM, et al. (2007) Identification of polymorphisms in the human 11beta-hydroxysteroid dehydrogenase type 2 gene promoter: functional characterization and relevance for salt sensitivity. Faseb J 21(13): 3618–28.R. Alikhani-KoupaeiF. FouladkouP. FustierB. CenniAM Sharma2007Identification of polymorphisms in the human 11beta-hydroxysteroid dehydrogenase type 2 gene promoter: functional characterization and relevance for salt sensitivity.Faseb J21(13)361828
48. Chakrabarty A, Roberts MR (2007) Ets-2 and C/EBP-beta are important mediators of ovine trophoblast Kunitz domain protein-1 gene expression in trophoblast. BMC Mol Biol 8: 14.A. ChakrabartyMR Roberts2007Ets-2 and C/EBP-beta are important mediators of ovine trophoblast Kunitz domain protein-1 gene expression in trophoblast.BMC Mol Biol814
49. Kraus RJ, Shadley L, Mertz JE (2001) Nuclear factor 1 family members mediate repression of the BK virus late promoter. Virology 287(1): 89–104.RJ KrausL. ShadleyJE Mertz2001Nuclear factor 1 family members mediate repression of the BK virus late promoter.Virology287(1)89104
50. Shabalina SA, Ogurtsov AY, Spiridonov NA (2006) A periodic pattern of mRNA secondary structure created by the genetic code. Nucleic Acids Res 34(8): 2428–37.SA ShabalinaAY OgurtsovNA Spiridonov2006A periodic pattern of mRNA secondary structure created by the genetic code.Nucleic Acids Res34(8)242837
51. Cervino AC, Li G, Edwards S, Zhu J, Laurie C, et al. (2005) Integrating QTL and high-density SNP analyses in mice to identify Insig2 as a susceptibility gene for plasma cholesterol levels. Genomics 86(5): 505–17.AC CervinoG. LiS. EdwardsJ. ZhuC. Laurie2005Integrating QTL and high-density SNP analyses in mice to identify Insig2 as a susceptibility gene for plasma cholesterol levels.Genomics86(5)50517
52. Champagne FA (2008) Epigenetic mechanisms and the transgenerational effects of maternal care. Front Neuroendocrinol 29(3): 386–97.FA Champagne2008Epigenetic mechanisms and the transgenerational effects of maternal care.Front Neuroendocrinol29(3)38697
53. Turri MG, DeFries JC, Henderson ND, Flint J (2004) Multivariate analysis of quantitative trait loci influencing variation in anxiety-related behavior in laboratory mice. Mamm Genome 15(2): 69–76.MG TurriJC DeFriesND HendersonJ. Flint2004Multivariate analysis of quantitative trait loci influencing variation in anxiety-related behavior in laboratory mice.Mamm Genome15(2)6976
54. Willis-Owen SA, Flint J (2006) The genetic basis of emotional behaviour in mice. Eur J Hum Genet 14(6): 721–8.SA Willis-OwenJ. Flint2006The genetic basis of emotional behaviour in mice.Eur J Hum Genet14(6)7218
55. Delaval K, Feil R (2004) Epigenetic regulation of mammalian genomic imprinting. Curr Opin Genet Dev 14(2): 188–95.K. DelavalR. Feil2004Epigenetic regulation of mammalian genomic imprinting.Curr Opin Genet Dev14(2)18895
56. Wilkins JF, Haig D (2003) What good is genomic imprinting: the function of parent-specific gene expression. Nat Rev Genet 4(5): 359–68.JF WilkinsD. Haig2003What good is genomic imprinting: the function of parent-specific gene expression.Nat Rev Genet4(5)35968
57. Jacobi F, Wittchen HU, Holting C, Hofler M, Pfister H, et al. (2004) Prevalence, co-morbidity and correlates of mental disorders in the general population: results from the German Health Interview and Examination Survey (GHS). Psychol Med 34(4): 597–611.F. JacobiHU WittchenC. HoltingM. HoflerH. Pfister2004Prevalence, co-morbidity and correlates of mental disorders in the general population: results from the German Health Interview and Examination Survey (GHS).Psychol Med34(4)597611
58. Mlynarik M, Zelena D, Bagdy G, Makara GB, Jezova D (2007) Signs of attenuated depression-like behavior in vasopressin deficient Brattleboro rats. Horm Behav 51(3): 395–405.M. MlynarikD. ZelenaG. BagdyGB MakaraD. Jezova2007Signs of attenuated depression-like behavior in vasopressin deficient Brattleboro rats.Horm Behav51(3)395405
59. Touma C, Bunck M, Glasl L, Nussbaumer M, Palme R, et al. (2008) Mice selected for high versus low stress reactivity: A new animal model for affective disorders. Psychoneuroendocrinology 33(6): 839–62.C. ToumaM. BunckL. GlaslM. NussbaumerR. Palme2008Mice selected for high versus low stress reactivity: A new animal model for affective disorders.Psychoneuroendocrinology33(6)83962
60. Neumann ID, Landgraf R (2008) Advances in vasopressin and oxytocin - from genes to behaviour to disease. Prog. Brain Res. 170, 1st ed. Oxford, Amsterdam: Elsevier. ID NeumannR. Landgraf2008Advances in vasopressin and oxytocin - from genes to behaviour to disease. Prog. Brain Res. 170, 1st edOxford, AmsterdamElsevier
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2009 Mirjam et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License (the “License”), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
Background
To investigate neurobiological correlates of trait anxiety, CD1 mice were selectively bred for extremes in anxiety-related behavior, with high (HAB) and low (LAB) anxiety-related behavior mice additionally differing in behavioral tests reflecting depression-like behavior.
Methodology/ Principal Findings
In this study, microarray analysis, in situ hybridization, quantitative real-time PCR and immunohistochemistry revealed decreased expression of the vasopressin gene (Avp) in the hypothalamic paraventricular (PVN) and supraoptic (SON) nuclei of adult LAB mice compared to HAB, NAB (normal anxiety-related behavior) and HABxLAB F1 intercross controls, without detecting differences in receptor expression or density. By sequencing the regions 2.5 kbp up- and downstream of the Avp gene locus, we could identify several polymorphic loci, differing between the HAB and LAB lines. In the gene promoter, a deletion of twelve bp Δ(−2180–2191) is particularly likely to contribute to the reduced Avp expression detected in LAB animals under basal conditions. Indeed, allele-specific transcription analysis of F1 animals revealed a hypomorphic LAB-specific Avp allele with a reduced transcription rate by 75% compared to the HAB-specific allele, thus explaining line-specific Avp expression profiles and phenotypic features. Accordingly, intra-PVN Avp mRNA levels were found to correlate with anxiety-related and depression-like behaviors. In addition to this correlative evidence, a significant, though moderate, genotype/phenotype association was demonstrated in 258 male mice of a freely-segregating F2 panel, suggesting a causal contribution of the Avp promoter deletion to anxiety-related behavior.
Discussion
Thus, the identification of polymorphisms in the Avp gene promoter explains gene expression differences in association with the observed phenotype, thus further strengthening the concept of the critical involvement of centrally released AVP in trait anxiety.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer