The endocrine stress response is linked to one specific locus on chromosome 3 in a mouse model based on extremes in trait anxiety
© Gonik et al.; licensee BioMed Central Ltd. 2012
Received: 19 August 2011
Accepted: 29 October 2012
Published: 31 October 2012
The hypothalamic-pituitary-adrenal (HPA) axis is essential to control physiological stress responses in mammals. Its dysfunction is related to several mental disorders, including anxiety and depression. The aim of this study was to identify genetic loci underlying the endocrine regulation of the HPA axis.
High (HAB) and low (LAB) anxiety-related behaviour mice were established by selective inbreeding of outbred CD-1 mice to model extremes in trait anxiety. Additionally, HAB vs. LAB mice exhibit comorbid characteristics including a differential corticosterone response upon stress exposure. We crossbred HAB and LAB lines to create F1 and F2 offspring. To identify the contribution of the endocrine phenotypes to the total phenotypic variance, we examined multiple behavioural paradigms together with corticosterone secretion-based phenotypes in F2 mice by principal component analysis. Further, to pinpoint the genomic loci of the quantitative trait of the HPA axis stress response, we conducted genome-wide multipoint oligogenic linkage analyses based on Bayesian Markov chain Monte Carlo approach as well as parametric linkage in three-generation pedigrees, followed by a two-dimensional scan for epistasis and association analysis in freely segregating F2 mice using 267 single-nucleotide polymorphisms (SNPs), which were identified to consistently differ between HAB and LAB mice as genetic markers.
HPA axis reactivity measurements and behavioural phenotypes were represented by independent principal components and demonstrated no correlation. Based on this finding, we identified one single quantitative trait locus (QTL) on chromosome 3 showing a very strong evidence for linkage (2ln (L-score) > 10, LOD > 23) and significant association (lowest Bonferroni adjusted p < 10-28) to the neuroendocrine stress response. The location of the linkage peak was estimated at 42.3 cM (95% confidence interval: 41.3 - 43.3 cM) and was shown to be in epistasis (p-adjusted < 0.004) with the locus at 35.3 cM on the same chromosome. The QTL harbours genes involved in steroid synthesis and cardiovascular effects.
The very prominent effect on stress-induced corticosterone secretion of the genomic locus on chromosome 3 and its involvement in epistasis highlights the critical role of this specific locus in the regulation of the HPA axis.
To warrant for adequate stress reactions, the hypothalamic-pituitary-adrenal (HPA) axis is activated upon exposure to stressors. As a consequence of such a reaction, corticosterone (CORT) and cortisol, respectively, are secreted in mammals. This is the first phase of the physiological stress response, which already initiates the depletion of further stress hormone secretion via negative feedback mechanisms[1–3]. As demonstrated by various preclinical and clinical studies, exaggerated responses of the HPA axis can lead to continuously elevated stress, considered as one underlying cause in the aetiology of anxiety and depression disorders in humans[3, 4]. Anxiety is also reflected at a cognitive level by the individuals' behaviour and strictly connected to physiological responses of central nervous, neuroendocrine and cardiovascular systems.
To shed light on the genetic underpinnings of the polygenic, multifactorial trait of anxiety, including its neuroendocrine correlates, multiple approaches have been applied up to date. Genetic engineering has helped to create knock-out and knock-in mice. These models have contributed a considerable amount of knowledge to describe the effects of neurotransmitter and neuromodulator systems, e.g., the corticotropin releasing hormone (CRH), vasopressin (AVP), tachykinin-class neuropeptides or serotonin (5-HT)[6–9]. However, the high number of regulatory mechanisms that are required to maintain the stress response and its return to homeostasis suggests that each molecular player contributes slightly to each phenotype, making it difficult to highlight all the mechanisms involved[10, 11].
Compared to genetic engineering, selective breeding strategies have the advantage of keeping the integrity of the genome and simulating the interplay between multiple systems, thus being closer to the clinical situation. In our laboratory, mice of the outbred CD-1 strain were selectively bred for > 20 generations, based on their anxiety-related behaviour on the elevated plus-maze (EPM), resulting in high (HAB) vs. low (LAB) anxiety-related behaviour lines. In addition to differences in anxiety-related behaviour, these mice also exhibit comorbid phenotypic differences in other behavioural tests, such as the dark/light box, open field (OF), tail-suspension (TST) or forced swim tests (FST), the latter two being indicative of depression-like behaviour. Hence, using the HAB/LAB mouse model furthered the identification of candidate genes of trait anxiety, including Avp, Enoph1, Glo1 and Tmem132d[13–16].
We here generated a HAB x LAB F2 intercross and focussed on endocrine system responses as a trait of interest. To reveal the link of genomic loci with the modulation of CORT responses upon stress exposure, we performed linkage analysis followed by an association study to specify the localisation of the QTL. Implementation of multipoint oligogenic segregation and linkage analysis by use of Bayesian Markov chain Monte Carlo (MCMC) method in the presence of missing values increases the power of the study by including the information from a total three-generation pedigree dataset and by sophisticated IBD sharing probability estimations. In the present study, the focus is on the genetic analysis of endocrine phenotypes, while other findings of the linkage of anxiety-related and depression-like behaviour phenotypes (Czibere, unpublished) are largely neglected at this stage. The findings of this study have translational potential for human studies dealing with the dysregulation of the HPA axis, which is regarded as one of the key phenotypes in major depression.
Animals were kept under standard housing conditions in groups of two to four per cage (clear plastic, Makrolon type II, 23 × 16 × 14 cm) with a 12-h light/dark cycle (lights on at 07:00 am). Breeding pairs and newborn litters were housed in Makrolon type III (38 × 22 × 15 cm) cages until weaning of the pups at four weeks of age.
All behavioural tests were carried out between 08:00 am and 01:00 pm to minimise circadian variations of HPA axis activity. To investigate different traits including anxiety-related behaviour, exploratory drive, locomotor activity, CORT secretion profile and depression-like behaviour, various test paradigms were used. Starting at the age of seven weeks, all animals (F0, F1 and F2) underwent a series of tests including EPM (5 min, 300 Lux on the open arms), TST (6 min), OF (5 min, 60 Lux in the central compartment), elevated platform test (EPF; 5 min), HPA axis reactivity test (HPA-RT; also called stress reactivity test), and FST (6 min, water at 23°C), as already described previously[12, 19–21]. Animals were tested once in each paradigm in the order described above with an intertest interval of 48 h. Throughout the test battery, we tried to avoid exposure of animals to two similar tests on two consecutive days. Therefore, as a test for depression-like behaviour, TST was performed 48 h after the first anxiety-related behaviour test (EPM) on day one. To examine long-lasting effects of the relatively stressful TST on the OF test, we exposed 9 CD-1 mice to TST stress two days prior to OF and compared them to 10 unexposed controls. An independent sample of 12 HAB and 11 LAB mice from generation 18 was used to reassure the phenotypic divergence of HAB vs. LAB mice in the HPA-RT.
Described just briefly, the procedure for the HPA-RT was as follows. Mice were taken out of their home cages and a first blood sample was taken from their tail veins using glass capillaries. The procedure was limited to 2 min to allow acquisition of basal (i.e., unstressed) values for CORT. The animals were then immobilised for 15 min in 50 ml plastic tubes (with holes at the front and end), and the tubes were covered by a lid to ensure darkness for the restraint stress period. After this 15-min stressor, a second blood sample was taken from a fresh cut less caudally than the first one. Blood samples were processed to obtain at least 10 μl of blood plasma, which was then quantified for CORT using a radioactive immunosorbent assay (DRG Diagnostics, Marburg, Germany) according to the manufacturer’s protocol. All of the samples were assayed in the same batch at the same time. According to the manufacturer’s recommendation, intra-assay variability was reduced to the least possible extent by measuring duplicates. Single measurement values were considered valid, if they were within the range of the standards supplied by the manufacturer.
All behavioural analyses were performed by an experimenter blind to the background of the animals and videotaped for evaluation. Body weight was measured directly after EPM testing. At the age of 14 weeks, animals were single-housed for four days, and their 24-h water consumption was measured as described previously. Then, the animals were decapitated under isoflurane anaesthesia and tail tips were collected for further analysis. Videotapes were scored by an observer blind to the background of the animals using “plus-maze” software for the EPM (E. Fricke, Munich, Germany), Anymaze (Stoelting Co, Wood Dale, IL) for the OF, and Eventlog 1.0 (EMCO software, Reykjavik, Iceland) for TST, EPF and FST, respectively. In the EPM test, an open arm entry is an entry with the forepaws only, corresponding to an entry of > 30% of the full body size, whereas a full open arm entry is an entry with the fore- and hind paws, thus corresponding to > 95% of the full body size.
The results of this study will focus on the phenotypic analyses and on the genetic linkage of CORT-related measurements from the HPA-RT, comprising the CORT concentrations before (basal CORT) and after (stress CORT) the 15-min restraint and the increase in CORT levels (stress CORT minus basal CORT) due to stressor exposure. As the increase in CORT is calculated as mentioned above, it is likely to be correlated, though not identical, with stress CORT, i.e., the two values will not be truly independent.
Selection of suitable single-nucleotide polymorphisms (SNPs) for HAB vs. LAB mice
In an unbiased whole-genome assay, SNPs that consistently differ between HAB and LAB animals (i.e., each line carries the opposite alleles homozygously) were identified, thus allowing for genotyping of F2 mice. To assess a high number of SNPs, the Mouse Medium Density Linkage Panel (Illumina, San Diego, CA) was chosen to determine 1,449 genotypes simultaneously, giving an overview of SNPs available as specific markers for HAB vs. LAB animals. Providing the basis for subsequent genotype-phenotype associations in F2 mice, all male and female F0 (N = 16 each) and some F1 mice (N = 18) were genotyped, with the latter serving as controls. For SNPs, where HAB and LAB mice would show the opposite genotype homozygously, F1 mice should exclusively show heterozygous genotypes.
DNA samples were isolated from tail tips using the NucleoSpin Tissue kit (Macherey-Nagel, Düren, Germany). Samples were processed according to the manufacturer’s instructions (Illumina Golden Gate Assay for Sentrix Array Matrix workflow sheets). Fluorescence signals of hybridised samples were captured in the Illumina BeadStation, and genotypes were called from fluorescence intensity clusters using the Illumina BeadStudio software Ver. 126.96.36.199 with the genotyping module Ver. 3.1.12. All intensity clusters were inspected individually and adjusted manually, if necessary. Analysis of genotypes focussed on the detection of valid genotypic SNP markers to distinguish HAB from LAB mice.
Genotyping of F2 mice
Based on the previously described experiment using the Medium Density Linkage Panel, a custom designed oligo pool (384 SNPs, Golden Gate Assay, Illumina) was set up for genotyping of all 521 F2 mice and, in addition, the 16 HAB and 16 LAB F0 mice, respectively. Therefore, 267 SNPs were chosen from the HAB vs. LAB mouse SNP identification experiment and another 116 SNPs were selected from the MGI database to fill up unmapped gaps and to screen selected genes from previously published or unpublished own studies (Additional file1: Table S1). One SNP (rs3659551) was genotyped twice serving as internal control. Evaluation was performed as described for the Medium Density Linkage Panel. All genomic locations mentioned refer to genome assembly NCBI37/Mm9.
Two groups comparison (Mann-Whitney-U test) was conducted for the parental F0 HAB and LAB mice to examine differences in the endocrine stress measures.
A principal component (PC) analysis was applied to examine the phenotypic structure of the F2 mice. PC analysis was used to reveal the phenotypes giving the highest loadings to the total phenotypic variance and to check for correlations between the HPA axis reactivity and the behavioural phenotypes. PCs selected by the stricter broken-stick criterion were kept for analysis. PC analysis was followed by an orthogonal varimax rotation, which resulted in each retained factor bearing a small number of large loadings and a large number of zero (or small) loadings, i.e., the variance of the loadings was maximised. Loadings with absolute values > 0.7 were treated as major loadings. PC analysis and varimax rotation were performed using the open-source statistical software package R (http://www.r-project.org/). PC analysis was carried out using the normalised phenotypic data of male mice only.
The Loki package[23, 24] for multipoint IBD estimations and oligogenic model for segregation and linkage analysis were used to benefit from the entire three-generation pedigree and to implement reversible jump MCMC techniques for likelihood estimations and Bayesian approaches. A quantitative trait was modelled as being genetically controlled by diallelic QTL and affected by the sex covariate. The description by Cox et al. served as the basis for conversion of physical to genomic distances on the autosomal chromosomes.
The Oligogenic model suggests the presence of multiple putative QTLs. The advantage of the Bayesian strategy, in contrast to other linkage methods, is that the number of putative QTLs is not fixed a priori. In contrast, this number was treated as an unknown model parameter. The posterior distribution for the number of underlying QTLs was calculated by using an iterative MCMC sampling technique.
The Bayes factor, or L-score, was computed to determine linkage as follows. For each sampling iteration (100,000 iterations were used for each analysis), the prior probability of finding a QTL linked to a 1 cM bin was 1 / t, where t was the total map length of the genome (approximately 1400 cM). If, for a particular iteration, there were n QTLs in the model, the prior probability, P, of at least one QTL located in the bin was 1 − (1 – 1 / t) n. The posterior probability, Q, is 1 or 0 depending on whether at least one QTL is located in the 1 cM bin. The score for each bin was estimated by averaging Q / P over an interval of iterations, which is supposed to perform a relatively stable outcome (75,000 – 100,000 iteration interval was used). Genomic regions that have a high probability of containing a QTL were expected to have considerably elevated L-scores compared to surrounding regions. The Bayes factors (L-scores) cannot be quantified in terms of a LOD score, but guidelines for estimating the importance of the Bayes factor are given by Kass and Raftery. Twice the natural logarithm of the Bayes factor is on the same scale, as are the familiar device and likelihood ratio test statistics. 2 log (L-score) scores of a range of 2 to 6 were considered positive signals, 6 to 10 were considered strong, and over 10 very strong.
To quantify the linkage signal in terms of more common LOD scores, multipoint parametric linkage analysis of quantitative traits was performed by MQScore_SNP software package. A priori, a diallelic autosomal QTL is supposed to be placed near each marker. Parameters of the model are estimated during linkage analysis by maximisation of the LOD score. This score compares the null hypothesis that SNP and QTL segregate independently with the alternative hypothesis, i.e., close linkage between QTL and each marker.
To estimate the location of the identified QTL in terms of common confidence interval and to test for epistatic interactions, QTL mapping in the F2 intercross was performed using the open-source statistical software package R/qtl. Since there is no standard procedure available to determine confidence intervals in multipoint linkage analysis, the single-QTL scan was performed by standard interval mapping (maximum likelihood algorithm). The confidence interval of QTL location was estimated by 95% Bayes credible intervals, which are supposed to depend less on marker density. The QTL-by-QTL epistatic interactions were characterized by the p-values for the LOD scores of corresponding epistatic models. The interactions were accepted to be significant if they did not exceed the adjusted threshold of 0.01. To assess the adjusted p-values, a traditional permutation test (N = 1000) was performed for the genomic region of interest.
To localise the most prominent markers more exactly, we additionally implemented the WG-Permer software (http://www.wg-permer.org) for rapid association analysis on the autosome chromosomes and R (http://www.r-project.org) for association analysis on X chromosome in male F2 mice. The Bonferroni correction for multiple testing was applied. The level of significance for all association tests was set to 10-4. An allelic model was used to make the results comparable with those of the linkage analysis mentioned above.
The gene functional classification tool available at the “Database for Annotation, Visualization and Integrated Discovery” (DAVID;http://david.abcc.ncifcrf.gov/)[32, 33] was used to retrieve the functional gene categories overrepresented among the genes situated in the region of the identified QTL. Analysis of the enriched functional-related gene groups was carried out with the DAVID-based Fisher's exact test, and results were Bonferroni corrected.
Characteristic phenotypes of HAB vs. LAB mice (F0 mice, N = 8 per gender and line); values are represented as mean ± standard error of means (SEM); corrected p-values refer to Bonferroni correction and are sorted by p-values of comparison in male mice
HAB males mean±SEM
LAB males mean±SEM
corrected p (males)
HAB females mean±SEM
LAB females mean±SEM
EPM - percent time on the open arm
EPM - percentage of open arm entries
TST - total immobility time [s]
EPM - closed arm entries
EPM - full open arm entries
EPM - open arm entries
TST - latency to 1st immobility [s]
PHYS - total fluid intake [ml]
OF - total distance travelled [m]
PHYS - daily water intake [ml/g]
PHYS - urine osmolality [osmol/kg]
EPF - frequency of head dips
HPA-RT - increase in CORT [μg/l]
EPM - latency to 1st open arm entry [s]
EPF - latency to 1st head dip [s]
HPA-RT - stress CORT [μg/l]
EPF - total grooming time [s]
FST - total swimming time [s]
FST - total floating time [s]
OF - total time in the inner zone [s]
EPF - frequency of grooming
FST - latency to 1st floating [s]
HPA-RT - basal CORT [μg/l]
OF - entries to the inner zone
EPM - all arm entries
FST - total struggling time [s]
PHYS - body weight [g]
EPM - sum of total arm entries
Principal component (PC) analysis of behavioural phenotypes in HAB x LAB derived F2 mice sorted by behavioural test paradigms
HPA-RT – stress CORT [μg/l]
HPA-RT – increase in CORT [μg/l]
EPM – entries to the middle open arm
EPM – entries to the distal open arm
EPM – time on the distal open arm [s]
EPM – time on the middle open arm [s]
EPM – full open arm entries
EPM – time spent on the proximal open arm [s]
EPM – percentage of open arm entries
EPM – percent time on the open arm
EPM – total time on the open arm [s]
EPM – all open arm entries
FST – total swimming time [s]
FST – total floating time [s]
OF – total distance travelled [m]
OF – mean travelling speed [m/s]
OF – total distance in the outer zone [m]
OF – entries to the outer zone
OF – total time in the outer zone [s]
OF – total time in the inner zone [s]
OF – total distance in the inner zone [m]
OF – latency to 1st entry of the outer zone [s]
PHYS – daily water intake per g bodyweight [ml/g*d]
PHYS – total daily water intake [ml/d]
Importantly, even if the preceding test is stressful, it is unlikely to have any impact on behaviour two days later. Indeed, CD-1 mice tested in the TST and CD-1 controls without prior testing showed similar behavioural indices in the OF (Additional file2: Figure S1).
The comorbidity of anxiety-related and depression-like phenotypes observed in HAB vs. LAB mouse comparisons was not reflected by the PC analysis in our F2 animals, leading to high loadings of anxiety-related behaviour as measured in the EPM in PC2, in the OF in PC3 and of depression-like behaviour in PC4 (Table2). In a similar manner, phenotypes reflecting locomotor activity or exploratory drive in the EPM and OF showed their highest loadings on two separate anxiety-independent PCs (PC5 and PC1, Table2).
The PC analysis clearly indicated that two of three endocrine parameters, which both reflect the CORT response upon restraint stress (stress CORT and increase in CORT), showed a high correlation (adjusted R2 = 0.99). These two endocrine parameters demonstrated the highest loadings (−0.97 and −0.96, respectively) on PC6. Hardly any other PC bore comparably high loadings except for the inner/outer zone time spent in the OF (PC3) and entries to the middle open arm in the EPM (PC1). Thus, PC6, accounting for 4.6% of the total phenotypic variance, is representative of CORT responses to a stressor as measured in the HPA-RT with basal CORT measurements, largely not affecting the overall outcome in F2 mice, leading to similar stress CORT and increase in CORT values. Based on these findings, any of these two endocrine measurements (stress CORT or increase in CORT) might be used as a compound score for the CORT response.
Genotyping of the F2 animals at 384 SNPs (average call rate 93%) resulted in 27 unidentifiable loci, 89 showed no informative value concerning the F2 panel, and one SNP was genotyped twice (showing identical results), leaving 267 SNPs for further testing with informative value for genotype-phenotype associations. From these markers, 233 SNPs were taken from the HAB vs. LAB mouse SNP identification experiment and 35 SNPs from those that were additionally screened on the 384 Sentrix Array Matrix (more than 25% of the newly added SNPs were informative for our mouse lines). Eleven markers were only specific for one LAB-subline, thereof two (rs29341895, rs29354500) from the list of SNPs added after the initial screening experiment. These 11 SNPs were not included in the phenotype-genotype linkage due to a shift of allele frequencies, leaving a total of 256 SNPs. This resulted in an average SNP marker density of 1 SNP per 11 Mbp. There were one genomic region of 53 Mbp on chromosome 3 and further five regions between 40 and 45 Mbp in size on chromosomes 2, 4, 10 and 6 for which we did not have any SNP markers. According to their genotypes, all F0 animals could clearly be clustered to their respective HAB or LAB origin.
Association analysis of markers on the X chromosome for hypothalamic pituitary adrenal axis activity-related phenotypes (p-values for increase in CORT, stress CORT and basal CORT)
increase in CORT
Importantly, to exclude parent-of-origin effects, calculations made separately for all F2 mice with a HAB F0 mother (N = 273) and with a LAB F0- mother (N = 261) did not show any differences in the effects shown on chromosome 3 (nominal p-value for offspring of HAB F0 mother rs13477268: p < 10-28; for offspring of LAB F0 mother: p < 10-25).
Panther functional enrichment analysis of the genomic region strongly linked to the CORT response phenotype
Enrichment score: 4.76
Bonferroni corrected p
androgen/ oestrogen/ progesterone biosynthesis
Our results highlight a specific genomic region on chromosome 3 as a major contributor to the variation of CORT response derived from the HAB/LAB mouse model. The strong contribution of genomic factors of this region seems to be essential for regulating stress-related reactivity of the HPA axis.
Here, we described first of all a difference in CORT secretion in HAB vs. LAB mice upon a 15-min restraint stressor, with LAB mice secreting more CORT compared to HAB mice. Thus, HAB mice appear to have a rather blunted stress response. Although, in general, high levels of CORT response are described to be associated with high anxiety states[34, 35], the opposite is true for HAB vs. LAB mice. This finding was replicated in CD-1 mice bred for differences in CORT responses to a stressor, with lower HPA axis reactivity mice showing more passive exploratory behaviour. In addition, a stronger increase in CORT has been demonstrated after social defeat in LAB than in HAB rats, probably due to the more active coping style of the former in response to any kind of stressor.
We could not detect any major difference in locomotion of HAB vs. LAB mice in the EPM, indicating that the animals’ anxiety-related phenotype is mainly independent of effects based on differing locomotion. Though, we cannot completely reject a confounding influence of locomotion on behaviour, especially as HAB vs. LAB male mice exhibited nominally significant differences in locomotion in the OF. On the other hand, we could also detect nominally significant differences in exploratory behaviour on the EPF for HAB vs. LAB males, a test much less affected by differences in locomotion.
Linkage studies were applied in rats and mice to a wide variety of complex phenotypes, including alcohol preference, anxiety-related, depression-like and exploratory behaviours or HPA axis function[37–40]. Based on the study by Williams et al., where a comparable breeding strategy led to the identification of about 300 genotypic markers useful for complex trait analysis, our mouse population seems to fulfil the requirements for linkage analysis. In addition, loading of the phenotypes of different behavioural tests in F2 on independent PCs underlines the good segregation of characteristic traits from the original HAB and LAB mouse lines.
Anxiety phenotypes of EPM and OF tests mainly loaded on two separate PCs showing limited comparability and predictability of measurements from different anxiety tests in freely segregating F2 animals. This is consistent with findings in other studies[42, 43] supporting the idea that different paradigms reflect – at least partially – different facets of trait anxiety. Similarly, in the F2 generation, phenotypes representing anxiety, depression and endocrine stress response loaded on separate PCs, reflecting strongly reduced comorbidity of these traits compared to F0 HAB and LAB mice. This resembles findings of another F2 study, where a connection between depression-like behaviour and stress reactivity was rejected. However, as solely based on low loadings of other behavioural phenotypes in the stress-related phenotypic cluster, a connection between these phenotypes cannot be completely denied, being in line with the idea of anxiety and stress reactivity representing multifactorial and polygenic traits with minor contributions of many genes. It is, therefore, likely that associated phenotypes and analogous comorbidities in the clinical situation such as those between anxiety, depression-like behaviours and HPA axis/stress reactivity share contributing genetic factors that are hardly detectable. This is also underlined by a recent study, proposing a network approach to analyse overlapping symptoms in the concept of comorbidity.
As stress CORT and increase in CORT are highly correlated (due to the low levels of basal CORT, resulting in almost near-to-identical values for the two parameters) and placed on the independent PC6 representing the only prominent loadings, both might be considered equally representative of the neuroendocrine system response. Stress CORT is exemplarily shown (Figure4) for genetic analysis. Identified by sequential linkage analysis of trait cosegregation in a full three-generation pedigree and F2 intercross mapping, a CORT-related QTL was detected on mouse chromosome 3 between 41.3 cM (closest marker: rs13477268 at 40.5 cM) and 43.3 cM (closest marker: rs4138887 at 45.3 cM) with a linkage peak at 42.3 cM. Two-QTL scans showed an evidence for epistasis between the region of the peak location at 42.3 cM and the locus at 35.3 cM (closest marker: rs6376008 at 37.9 cM). The association analysis of freely segregating F2 mice derived from a HAB x LAB cross identified a prominent trait co-occurrence with rs13477268 and two highly correlated nearby markers, rs4138887 and rs6376008. The high correlation of nearby SNPs (rs13477268, rs4138887 and rs6376008) points to the presence of a broad linkage disequilibrium (LD) block, indicating that the identified group of significant markers represents the effect of the true locus underlying the neuroendocrine stress response.
Taken together, our findings highlight the region from 37.9 cM to 45.2 cM of mouse chromosome 3 as an endocrine stress response-specific locus. The additionally identified regions on the X chromosome, only showing nominally significant association with the traits stress CORT and increase in CORT, provide an impetus for the involvement of sex-specific mechanisms in influencing the respective trait, which is supported by other studies.
Comparative genome analyses of mouse chromosome 3 showed rat and human homologues on chromosomes 2 (rat) and 1 and 4 (human), respectively. The markers that were linked to stress CORT and increase in CORT are not located near to any neurotransmitter receptor gene or other candidate genes for HPA axis regulation mapped to date in mice (according to NCBI, Ensembl and UCSC databases).
A particularly interesting candidate gene for the trait of neuroendocrine stress response is Pde4dip, a locus characterised by the maximum likelihood estimate of QTL location and involved in epistasis. The encoded protein, myomegalin, is responsible for cardiac contractility during adrenergic stimulation. The rat homologues of our identified QTL (Figure4A) are highly overlapping with previously identified QTLs for cardiovascular disease-related phenotypes. These overlaps include rat QTLs for blood pressure and heart rate and bradycardia. As these rat QTLs also include the genomic locus of Pde5a, which is capable of modulating acute and chronic cardiac stress responses[50, 51], represented by rs13477379 in our study, it makes this gene an interesting additional candidate, although not included in the region of strongest linkage. Particularly as acute CORT responses show a pronounced effect on the cardiovascular system, the previously described cardiovascular QTLs in this genomic region are in line with our findings.
The identified QTL is located in a region that also overlaps with a QTL for ethanol-induced CORT responses in mice and is homologous to a QTL for fasting cortisol levels in humans. Both ethanol administration and fasting have been shown to provoke an increase in CORT secretion and glucocorticoid production, similar to acute stress responses[55, 56]. Finally, the 7-cM region, identified to be linked to stress CORT in the present study, is completely homologous with one of the QTLs (Srcrt-1) Solberg et al. described to influence stress CORT in their rat model, with the difference that the linkage in our study can be limited to a single, smaller region, which is strongly linked to the stress CORT phenotype. In case of Solberg et al., Srcrt-1 shows the second strongest effect on stress CORT, behind a marker on rat chromosome 6 (homologous to mouse chromosome 12) of five linked genomic loci altogether.
The strongly linked 7-cM genomic region on chromosome 3 harbours at least five genes involved in steroid, most prominently glucocorticoid, synthesis (Table4), emphasising a strong effect on the stress CORT phenotype. Although these genes are primarily focussed on androgen, oestrogen and progesterone synthesis (Table4), their enzyme products are known to play a vital role in the synthesis of all glucocorticoids, among them CORT (KEGG pathway mmu00140;http://www.genome.jp/kegg/). This strongly supports a conserved genomic and, thereby, heritable effect on stress reactivity. Interestingly, a recent study indicates that corticosteroids can be synthesised endogenously in the brain. Therefore, the genomic region identified in our study might also affect corticosteroids directly in the brain. The second functionally enriched group (Table4), albeit showing only nominal significance, is related to cholesterol biosynthesis, which is the first metabolite required for glucocorticoid synthesis. Another genomic site containing genes coding for glutamate receptors and other transmembrane proteins (NCBI) is related to the locus at 35.3 cM, which is in epistatic interaction with the identified locus at 42.3 cM.
While many mechanisms of genomic and non-genomic regulation remain largely unclear in terms of activation vs. inhibition of the HPA axis, our results strongly point to a genetic influence of the 7-cM region. The additionally identified regions on the X chromosome, only showing nominally significant linkage to the stress CORT trait, provide an impetus for the involvement of sex-specific mechanisms in influencing the respective trait, which is supported by other studies.
Taken together, based on a HAB x LAB mouse intercross panel, we have identified one major QTL on chromosome 3 linked to both stress-induced CORT level and the level of increase in CORT, providing a strong basis for the heritability of neuroendocrine stress reactivity and emphasising a critical role of the genes relevant for steroid synthesis. In addition to the neuroendocrine effect, Pde4dip covering the position of the linkage peak and Pde5a placed in a nearby region overlapping with known cardiovascular QTLs, suggest a comorbid cardiovascular effect driven by the identified locus. Further supported by epistatic effects, our findings highlight a novel aspect towards the understanding of the well-orchestrated mechanisms of stress response regulation. They will contribute to further facilitate the search for candidate genes that are important for the regulation of neuroendocrine systems under stress conditions and, simultaneously, provide deeper insights into their phenotypic outcome, behaviour.
The authors would like to thank Markus Nußbaumer, Marina Zimbelmann, Sabine Damast, and Susann Sauer for perfect technical assistance and Ramona Zeh, Hendrik Stein, Julia Baier, Jan-Michael Heinzmann and Kathrin Halfter for their help with phenotyping. Last but not least, the authors would like to thank the anonymous reviewers for their suggestions to improve this manuscript.
All authors and the study were funded by the Max Planck Society. The funder had neither a role in study design, collection, analysis and interpretation of data, in writing the manuscript, nor in the decision to submit the manuscript for publication.
- Herman JP, Cullinan WE: Neurocircuitry of stress: central control of the hypothalamo-pituitary-adrenocortical axis. Trends Neurosci. 1997, 20: 78-84. 10.1016/S0166-2236(96)10069-2.View ArticlePubMedGoogle Scholar
- Ardayfio P, Kim KS: Anxiogenic-like effect of chronic corticosterone in the light–dark emergence task in mice. Behav Neurosci. 2006, 120: 249-256.View ArticlePubMedGoogle Scholar
- Kudielka BM, Wust S: Human models in acute and chronic stress: assessing determinants of individual hypothalamus-pituitary-adrenal axis activity and reactivity. Stress. 2010, 13: 1-14. 10.3109/10253890902874913.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle Scholar
- Rivest S: Interactions between the immune and neuroendocrine systems. Prog Brain Res. 2010, 181: 43-53.View ArticlePubMedGoogle Scholar
- Severini C, Improta G, Falconieri-Erspamer G, Salvadori S, Erspamer V: The tachykinin peptide family. Pharmacol Rev. 2002, 54: 285-322. 10.1124/pr.54.2.285.View ArticlePubMedGoogle Scholar
- Anisman H, Merali Z, Stead JD: Experiential and genetic contributions to depressive- and anxiety-like disorders: clinical and experimental studies. Neurosci Biobehav Rev. 2008, 32: 1185-1206. 10.1016/j.neubiorev.2008.03.001.View ArticlePubMedGoogle Scholar
- Keszthelyi D, Troost FJ, Masclee AA: Understanding the role of tryptophan and serotonin metabolism in gastrointestinal function. Neurogastroenterol Motil. 2009, 21: 1239-1249. 10.1111/j.1365-2982.2009.01370.x.View ArticlePubMedGoogle Scholar
- Hodges MR, Richerson GB: The role of medullary serotonin (5-HT) neurons in respiratory control: contributions to eupneic ventilation, CO2 chemoreception, and thermoregulation. J Appl Physiol. 2010, 108: 1425-1432. 10.1152/japplphysiol.01270.2009.PubMed CentralView ArticlePubMedGoogle Scholar
- Flint J, Valdar W, Shifman S, Mott R: Strategies for mapping and cloning quantitative trait genes in rodents. Nat Rev Genet. 2005, 6: 271-286.View ArticlePubMedGoogle Scholar
- Belzung C, Philippot P: Anxiety from a phylogenetic perspective: is there a qualitative difference between human and animal anxiety?. Neural Plast. 2007, 2007: 59676-PubMed CentralPubMedGoogle Scholar
- Kromer SA, Kessler MS, Milfay D, Birg IN, Bunck M, Czibere L, Panhuysen M, Putz B, Deussing JM, Holsboer F, Landgraf R, Turck CW: Identification of glyoxalase-I as a protein marker in a mouse model of extremes in trait anxiety. J Neurosci. 2005, 25: 4375-4384. 10.1523/JNEUROSCI.0115-05.2005.View ArticlePubMedGoogle Scholar
- Bunck M, Czibere L, Horvath C, Graf C, Frank E, Kessler MS, Murgatroyd C, Muller-Myhsok B, Gonik M, Weber P, Putz B, Muigg P, Panhuysen M, Singewald N, Bettecken T, Deussing JM, Holsboer F, Spengler D, Landgraf R: A hypomorphic vasopressin allele prevents anxiety-related behavior. PLoS One. 2009, 4: e5129-10.1371/journal.pone.0005129.PubMed CentralView ArticlePubMedGoogle Scholar
- Ditzen C, Varadarajulu J, Czibere L, Gonik M, Targosz BS, Hambsch B, Bettecken T, Kessler MS, Frank E, Bunck M, Teplytska L, Erhardt A, Holsboer F, Muller-Myhsok B, Landgraf R, Turck CW: Proteomic-based genotyping in a mouse model of trait anxiety exposes disease-relevant pathways. Mol Psychiatry. 2010, 15: 702-711. 10.1038/mp.2008.146.View ArticlePubMedGoogle Scholar
- Erhardt A, Czibere L, Roeske D, Lucae S, Unschuld PG, Ripke S, Specht M, Kohli MA, Kloiber S, Ising M, Heck A, Pfister H, Zimmermann P, Lieb R, Putz B, Uhr M, Weber P, Deussing JM, Gonik M, Bunck M, Kessler MS, Frank E, Hohoff C, Domschke K, Krakowitzky P, Maier W, Bandelow B, Jacob C, Deckert J, Schreiber S, et al, et al: TMEM132D, a new candidate for anxiety phenotypes: evidence from human and mouse studies. Mol Psychiatry. 2011, 16: 647-663. 10.1038/mp.2010.41.View ArticlePubMedGoogle Scholar
- Kessler MS, Murgatroyd C, Bunck M, Czibere L, Frank E, Jacob W, Horvath C, Muigg P, Holsboer F, Singewald N, Spengler D, Landgraf R: Diabetes insipidus and, partially, low anxiety-related behaviour are linked to a SNP-associated vasopressin deficit in LAB mice. Eur J Neurosci. 2007, 26: 2857-2864. 10.1111/j.1460-9568.2007.05917.x.View ArticlePubMedGoogle Scholar
- Almasy L, Blangero J: Human QTL linkage mapping. Genetica. 2009, 136: 333-340. 10.1007/s10709-008-9305-3.PubMed CentralView ArticlePubMedGoogle Scholar
- Cryan JF, Holmes A: The ascent of mouse: advances in modelling human depression and anxiety. Nat Rev Drug Discov. 2005, 4: 775-790. 10.1038/nrd1825.View ArticlePubMedGoogle Scholar
- Czibere L, Baur LA, Wittmann A, Gemmeke K, Steiner A, Weber P, Pütz B, Ahmad N, Bunck M, Graf C, Widner R, Kühne C, Panhuysen M, Hambsch B, Rieder G, Reinheckel T, Peters C, Holsboer F, Landgraf R, Deussing JM: Profiling trait anxiety: Transcriptome analysis reveals cathepsin B (Ctsb) as a novel candidate gene for emotionality in mice. PLoS One. 2011, 6: e23604-10.1371/journal.pone.0023604.PubMed CentralView ArticlePubMedGoogle Scholar
- Touma C, Bunck M, Glasl L, Nussbaumer M, Palme R, Stein H, Wolferstatter M, Zeh R, Zimbelmann M, Holsboer F, Landgraf R: Mice selected for high versus low stress reactivity: a new animal model for affective disorders. Psychoneuroendocrinology. 2008, 33: 839-862. 10.1016/j.psyneuen.2008.03.013.View ArticlePubMedGoogle Scholar
- Frank E, Kessler MS, Filiou MD, Zhang Y, Maccarrone G, Reckow S, Bunck M, Heumann H, Turck CW, Landgraf R, Hambsch B: Stable isotope metabolic labeling with a novel N-enriched bacteria diet for improved proteomic analyses of mouse models for psychopathologies. PLoS One. 2009, 4: e7821-10.1371/journal.pone.0007821.PubMed CentralView ArticlePubMedGoogle Scholar
- Jackson DA: Stopping rules in principal components analysis: a comparison of heuristical and statistical approaches. Ecology. 1993, 74: 2204-2214. 10.2307/1939574.View ArticleGoogle Scholar
- Heath SC: Markov chain Monte Carlo segregation and linkage analysis for oligogenic models. Am J Hum Genet. 1997, 61: 748-760. 10.1086/515506.PubMed CentralView ArticlePubMedGoogle Scholar
- Heath SC, Snow GL, Thompson EA, Tseng C, Wijsman EM: MCMC segregation and linkage analysis. Genet Epidemiol. 1997, 14: 1011-1016. 10.1002/(SICI)1098-2272(1997)14:6<1011::AID-GEPI75>3.0.CO;2-L.View ArticlePubMedGoogle Scholar
- Cox A, Ackert-Bicknell CL, Dumont BL, Ding Y, Bell JT, Brockmann GA, Wergedal JE, Bult C, Paigen B, Flint J, Tsaih SW, Churchill GA, Broman KW: A new standard genetic map for the laboratory mouse. Genetics. 2009, 182: 1335-1344. 10.1534/genetics.109.105486.PubMed CentralView ArticlePubMedGoogle Scholar
- Green P: Reversible jump Markov chain Monte Carlo. Biometrika. 1995, 82: 711-732. 10.1093/biomet/82.4.711.View ArticleGoogle Scholar
- Kass RE, Raftery AE: Bayes factors. J Am Stat Assoc. 1995, 90: 773-795. 10.1080/01621459.1995.10476572.View ArticleGoogle Scholar
- Dietrich W, Katz H, Lincoln SE, Shin HS, Friedman J, Dracopoli NC, Lander ES: A genetic map of the mouse suitable for typing intraspecific crosses. Genetics. 1992, 131 (2): 423-447.PubMed CentralPubMedGoogle Scholar
- Axenovich TI, Aulchenko YS: MQScore_SNP software for multipoint parametric linkage analysis of quantitative traits in large pedigrees. Ann Hum Genet. 2010, 74: 286-289. 10.1111/j.1469-1809.2010.00576.x.View ArticlePubMedGoogle Scholar
- Broman KW, Wu H, Sen S, Churchill GA: R/qtl: QTL mapping in experimental crosses. Bioinformatics. 2003, 19: 889-890. 10.1093/bioinformatics/btg112.View ArticlePubMedGoogle Scholar
- Broman KW, Sen S: A Guide to QTL Mapping in R/QTL. 2009, New York: SpringerView ArticleGoogle Scholar
- Huang DW, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID Bioinformatics Resources. Nature Protoc. 2009, 4: 44-57.View ArticleGoogle Scholar
- Huang DW, Sherman BT, Lempicki RA: Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009, 37: 1-13. 10.1093/nar/gkn923.PubMed CentralView ArticleGoogle Scholar
- Maccari S, Darnaudery M, Morley-Fletcher S, Zuena AR, Cinque C, van Reeth O: Prenatal stress and long-term consequences: implications of glucocorticoid hormones. Neurosci Biobehav Rev. 2003, 27: 119-127. 10.1016/S0149-7634(03)00014-9.View ArticlePubMedGoogle Scholar
- Wigger A, Sanchez MM, Mathys KC, Ebner K, Frank E, Liu D, Kresse A, Neumann ID, Holsboer F, Plotsky PM, Landgraf R: Alterations in central neuropeptide expression, release, and receptor binding in rats bred for high anxiety: critical role of vasopressin. Neuropsychopharmacology. 2004, 29: 1-14. 10.1038/sj.npp.1300290.View ArticlePubMedGoogle Scholar
- Frank E, Salchner P, Aldag JM, Salome N, Singewald N, Landgraf R, Wigger A: Genetic predisposition to anxiety-related behavior determines coping style, neuroendocrine responses, and neuronal activation during social defeat. Behav Neurosci. 2006, 120: 60-71.View ArticlePubMedGoogle Scholar
- Bice P, Valdar W, Zhang L, Liu L, Lai D, Grahame N, Flint J, Li TK, Lumeng L, Foroud T: Genomewide SNP screen to detect quantitative trait loci for alcohol preference in the high alcohol preferring and low alcohol preferring mice. Alcohol Clin Exp Res. 2009, 33: 531-537. 10.1111/j.1530-0277.2008.00866.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Salimov RM, Markina NV, Perepelkina OV, Poletaeva II: Exploratory behavior of F2 crosses of mouse lines selected for different brain weight: a multivariate analysis. Prog Neuropsychopharmacol Biol Psychiatry. 2004, 28: 583-589. 10.1016/j.pnpbp.2004.01.016.View ArticlePubMedGoogle Scholar
- Solberg LC, Ahmadiyeh N, Baum AE, Vitaterna MH, Takahashi JS, Turek FW, Redei EE: Depressive-like behavior and stress reactivity are independent traits in a Wistar Kyoto x Fisher 344 cross. Mol Psychiatry. 2003, 8: 423-433. 10.1038/sj.mp.4001255.View ArticlePubMedGoogle Scholar
- Solberg LC, Baum AE, Ahmadiyeh N, Shimomura K, Li R, Turek FW, Takahashi JS, Churchill GA, Redei EE: Genetic analysis of the stress-responsive adrenocortical axis. Physiol Genomics. 2006, 27: 362-369. 10.1152/physiolgenomics.00052.2006.View ArticlePubMedGoogle Scholar
- Williams RW, Bennett B, Lu L, Gu J, DeFries JC, Carosone-Link PJ, Rikke BA, Belknap JK, Johnson TE: Genetic structure of the LXS panel of recombinant inbred mouse strains: a powerful resource for complex trait analysis. Mamm Genome. 2004, 15: 637-647. 10.1007/s00335-004-2380-6.View ArticlePubMedGoogle Scholar
- van der Staay FJ, Schuurman T, van Reenen CG, Korte SM: Emotional reactivity and cognitive performance in aversively motivated tasks: a comparison between four rat strains. Behav Brain Funct. 2009, 5: 50-10.1186/1744-9081-5-50.PubMed CentralView ArticlePubMedGoogle Scholar
- Post AM, Weyers P, Holzer P, Painsipp E, Pauli P, Wultsch T, Reif A, Lesch KP: Gene-environment interaction influences anxiety-like behavior in ethologically based mouse models. Behav Brain Res. 2011, 218: 99-105. 10.1016/j.bbr.2010.11.031.View ArticlePubMedGoogle Scholar
- Cramer AO, Waldorp LJ, van der Maas HL, Borsboom D: Comorbidity: a network perspective. Behav Brain Sci. 2010, 33: 137-150. 10.1017/S0140525X09991567.View ArticlePubMedGoogle Scholar
- Walker JJ, Terry JR, Tsaneva-Atanasova K, Armstrong SP, McArdle CA, Lightman SL: Encoding and decoding mechanisms of pulsatile hormone secretion. J Neuroendocrinol. 2010, 22: 1226-1238. 10.1111/j.1365-2826.2010.02087.x.View ArticlePubMedGoogle Scholar
- Brudno M, Poliakov A, Salamov A, Cooper GM, Sidow A, Rubin EM, Solovyev V, Batzoglou S, Dubchak I: Automated whole-genome multiple alignment of rat, mouse, and human. Genome Res. 2004, 14: 685-692. 10.1101/gr.2067704.PubMed CentralView ArticlePubMedGoogle Scholar
- Uys G, Ramburan A, Loos B, Kinnear C, Korkie L, Mouton J, Riedemann J, Moolman-Smook J: Myomegalin is a novel A-kinase anchoring protein involved in the phosphorylation of cardiac myosin binding protein C. BMC Cell Biol. 2011, 12: 18-10.1186/1471-2121-12-18.PubMed CentralView ArticlePubMedGoogle Scholar
- Rapp JP: Genetic analysis of inherited hypertension in the rat. Physiol Rev. 2000, 80: 135-172.PubMedGoogle Scholar
- Alemayehu A, Breen L, Krenova D, Printz MP: Reciprocal rat chromosome 2 congenic strains reveal contrasting blood pressure and heart rate QTL. Physiol Genomics. 2002, 10: 199-210.View ArticlePubMedGoogle Scholar
- Zhang M, Kass DA: Phosphodiesterases and cardiac cGMP: evolving roles and controversies. Trends Pharmacol Sci. 2011, 32: 360-365. 10.1016/j.tips.2011.02.019.PubMed CentralView ArticlePubMedGoogle Scholar
- Takimoto E, Belardi D, Tocchetti CG, Vahebi S, Cormaci G, Ketner EA, Moens AL, Champion HC, Kass DA: Compartmentalization of cardiac beta-adrenergic inotropy modulation by phosphodiesterase type 5. Circulation. 2007, 115: 2159-2167. 10.1161/CIRCULATIONAHA.106.643536.View ArticlePubMedGoogle Scholar
- Bechtold AG, Scheuer DA: Glucocorticoids act in the dorsal hindbrain to modulate baroreflex control of heart rate. Am J Physiol Regul Integr Comp Physiol. 2006, 290: R1003-R1011.View ArticlePubMedGoogle Scholar
- Roberts AJ, Phillips TJ, Belknap JK, Finn DA, Keith LD: Genetic analysis of the corticosterone response to ethanol in BXD recombinant inbred mice. Behav Neurosci. 1995, 109: 1199-1208.View ArticlePubMedGoogle Scholar
- Ober C, Abney M, McPeek MS: The genetic dissection of complex traits in a founder population. Am J Hum Genet. 2001, 69: 1068-1079. 10.1086/324025.PubMed CentralView ArticlePubMedGoogle Scholar
- Glover M, Pruett SB: Role of corticosterone in immunosuppressive effects of acute ethanol exposure on Toll-like receptor mediated cytokine production. J Neuroimmune Pharmacol. 2006, 1: 435-442. 10.1007/s11481-006-9037-z.View ArticlePubMedGoogle Scholar
- Rohner-Jeanrenaud F: Hormonal regulation of energy partitioning. Int J Obes Relat Metab Disord. 2000, 24 (Suppl 2): S4-S7.View ArticlePubMedGoogle Scholar
- Clarke TR, Bain PA, Burmeister M, Payne AH: Isolation and characterization of several members of the murine Hsd3b gene family. DNA Cell Biol. 1996, 15: 387-399. 10.1089/dna.1996.15.387.View ArticlePubMedGoogle Scholar
- Higo S, Hojo Y, Ishii H, Komatsuzaki Y, Ooishi Y, Murakami G, Mukai H, Yamazaki T, Nakahara D, Barron A, Kimoto T, Kawato S: Endogenous synthesis of corticosteroids in the hippocampus. PLoS One. 2011, 6: e21631-10.1371/journal.pone.0021631.PubMed CentralView ArticlePubMedGoogle Scholar
- Miller WL, Auchus RJ: The molecular biology, biochemistry, and physiology of human steroidogenesis and its disorders. Endocr Rev. 2011, 32: 81-151. 10.1210/er.2010-0013.PubMed CentralView ArticlePubMedGoogle Scholar
- Chrousos GP: Stress and sex versus immunity and inflammation. Sci Signal. 2010, 3: pe36-10.1126/scisignal.3143pe36.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.