Testosterone (T) regulates many traits related to fitness, including aggression. However, individual variation in aggressiveness does not always relate to circulating T, suggesting that behavioural variation may be more closely related to neural sensitivity to steroids, though this issue remains unresolved. To assess the relative importance of circulating T and neural steroid sensitivity in predicting behaviour, we measured aggressiveness during staged intrusions in free-living male and female dark-eyed juncos (Junco hyemalis). We compared aggressiveness to plasma T levels and to the abundance of androgen receptor (AR), aromatase (AROM) and oestrogen receptor alpha (ORα) mRNA in behaviourally relevant brain areas (avian medial amygdala, hypothalamus and song control regions). We also asked whether patterns of covariation among behaviour and endocrine parameters differed in males and females, anticipating that circulating T may be a better predictor of behaviour in males than in females. We found that circulating T related to aggressiveness only in males, but that gene expression for ORα, AR and AROM covaried with individual differences in aggressiveness in both sexes. These findings are among the first to show that individual variation in neural gene expression for three major sex steroid-processing molecules predicts individual variation in aggressiveness in both sexes in nature. The results have broad implications for our understanding of the mechanisms by which aggressive behaviour may evolve.
In vertebrates, many complex phenotypes that affect survival or reproduction are mediated by the steroid hormone testosterone (T) [1–3]. T plays a central role in translating genotype to phenotype in both males and females [4–6], and this hormone is a likely target of selection in the evolution of hormone-mediated behaviours, such as intrasexual aggression [7,8]. However, individual variation in circulating T often does not relate to individual differences in aggression in both sexes [2,9] (see also [10–12]). Because individual variation is the raw material of evolution, these observations call into question the role of T in the evolution of aggressive behaviour. A hypothesis that may resolve this issue suggests that variability in the cellular and molecular properties of target tissues (i.e. measures of sensitivity to T and its metabolites) are critical to explaining functional variation in behaviour (e.g. gene expression and protein abundance for molecules that process hormones in the brain [13,14]).
Androgen receptor (AR), oestrogen receptors (ORs) and aromatase (AROM) are among the key molecules in initiating the genomic effects of T [15–17]: after binding with a steroid, AR and OR transcriptionally regulate other genes that ultimately affect physiology and behaviour, and AROM (the enzyme converting T to oestradiol) enables T to activate OR-dependent pathways. Variation in abundance of these sex-steroid processing molecules is therefore likely to affect the expression of T-mediated behaviours such as aggression. Measurement of the abundance of AR, OR or AROM transcript or protein may therefore estimate critical aspects of local steroid sensitivity, though other co-factors and repressors also play a role . Many different experimental manipulations and ‘group comparisons’ have linked AR, AROM and OR in the brain to aggressive and social behaviour (e.g. by comparing males and females, more aggressive and less aggressive species, or receptor antagonist-treated individuals and controls [18–22]; see  for review). It has been widely hypothesized that these group differences can be extrapolated to individual differences [9,13,23,24]—for example, that more aggressive individuals will have a greater abundance of AR, AROM or OR alpha (ORα) in behaviourally relevant brain areas. However, most studies linking aggressive behaviour to individual variation in sensitivity to T have pooled control and experimentally manipulated animals, and efforts to relate these parameters among unmanipulated or free-living animals are particularly rare [25,26] (but see [27–29]).
To rectify this knowledge gap and better integrate evolutionary biology with neuroendocrinology, we employed an approach rooted in quantitative genetics [30,31]. We first considered the degree to which aggressive behaviour relates to individual variation in hormone signal, to measures of sensitivity to that signal in target tissues, or to some combination of the two. Because same-sex aggressive behaviour is an important predictor of reproductive success in females as well as males [32,33], we also asked whether the mechanisms underlying individual variation in aggression are similar or different in males and females. If selection acts on variation in signal strength, we might expect correlated responses to selection in both sexes, but if selection acts on variation in sensitivity, then the sexes may be more free to evolve independently [1,6,34].
Species comparisons have shown that male and female T levels are often positively correlated, which suggests that selection on T-mediated traits in one sex could lead to a correlated response in the other sex [6,35,36]. Experimentally or naturally elevated T can be costly , and thus we might expect selection to favour mechanisms of aggression that are somewhat independent of circulating T, particularly when circulating T is low or the costs of T are high (e.g. during the non-breeding season or in females [18,38,39]). Because the costs of T are thought to be particularly high for females [40,41], we predicted that mechanisms of aggression that do not depend upon T levels themselves may be particularly important for females, though experimentally elevated T has been shown to affect female aggression in this species .
We measured natural variation in aggressiveness towards a same-sex conspecific in male and female free-living dark-eyed juncos (Junco hyemalis) early in the breeding season. The dark-eyed junco is a North American sparrow that is well studied with respect to hormones, behaviour and sex differences [3,6,34]. We compared individual differences in aggressiveness with circulating levels of T and with neural gene expression for AR, AROM and ORα. We used real-time quantitative PCR to quantify these measures of sensitivity to T in three socially relevant brain regions [43,44]: the hypothalamus (hypo), ventromedial telencephalon (VmT, a dissection that is largely limited to the avian medial amygdala or nucleus taeniae) and the right posterior telencephalon (PTR, which includes song control nuclei). We asked whether individual differences in behaviour within each sex covaried with the relative abundance of AR, AROM or ORα gene expression, with the prediction that more aggressive individuals would express greater transcript abundance for these genes. Further, owing to greater potential costs of circulating T to females, we predicted that female behaviour would be less a function of circulating T compared with male behaviour. Because changes in gene expression are thought to be a major driver of evolutionary change in complex phenotypes, including aggression [45–47], we interpret the results in terms of their applicability to the mechanisms by which aggressive behaviour may evolve.
2. Material and methods
(a) Behavioural quantification and tissue collection
All subjects were captured during the breeding season (1 May–5 June 2010) in the area surrounding Mountain Lake Biological Station in Virginia, USA (37°22′ N, 80°32′ W), immediately after a short (6 min) simulated territorial intrusion (modified from [10,12]). All females were in the incubation stage. All males were in breeding condition as evidenced by date and enlarged gonads. Prior to the day of experimentation, territories were mapped using a combination of conspecific playback and observation of natural singing perches. During each aggression test, a live, caged same-sex decoy was placed in the approximate centre of the territory (for males) or 1 m from the nest (for females). To simulate natural male–male aggressive encounters, male aggression tests were accompanied with conspecific song broadcast at a natural rate and amplitude (six songs per minute; 85–90 dB at 1 m). Before each trial, a mist net was set up and furled. An observer retreated approximately 15 m before uncovering the live decoy and beginning the trial.
We recorded the following behaviours: the number of times the focal bird attacked the decoy by swooping (‘flyovers’), distance to the decoy, number of songs (for males) and amount of time responding to the decoy (for females; i.e. time spent responding to the intruder instead of incubating). These behavioural parameters are measures of aggressive response or the likelihood of attack in this species and other songbirds [10,48,49].
At the conclusion of the 6-min trial, the mist net was unfurled and we quickly captured the bird (4.5 ± 0.5 min, range: 1.5–15 min after trial) and euthanized it immediately with an overdose of isoflurane followed by decapitation (2.2 ± 0.6 min, range: 1.5–3 min after capture). All subjects appeared within 2.25 min of beginning the trial, and we were successful in capturing all birds. Thus, the possibility that the sample is biased towards only the most aggressive individuals or those that were easy to capture is unlikely (latencies to appear = 0.25 ± 0.08 min). Brains were dissected from the skull using RNAse-free tools and frozen in powdered dry ice within 5 min. All tissues were stored at –80°C, until microdissection.
(b) Molecular methods
In the laboratory, brains were microdissected into functional regions, following Soma et al. [50,51]. After removing the optic tecta, optic chiasm and the hindbrain at the level of the mammillary bodies, we isolated the hypo to the depth of the anterior commissure (including the preoptic area and ventromedial hypothalamus). We collected the ventromedial telencelphalon (VmT) by removing ∼1 mm of the ventromedial portion of the caudal telencephalon, based on the position of nucleus taeniae in other songbirds , and we isolated tissues containing the song control nuclei by focusing on the PTR. We focused on these tissues because they include nuclei that are rich sites for steroid-mediated regulation of social, mating and aggressive behaviour [52–54]. We extracted RNA using the Trizol method (Invitrogen, Carlsbad, CA). After treating each RNA sample with DNAse (Promega, Madison, WI), we used Superscript III (Invitrogen) to synthesize cDNA from each brain area for each individual, from 1 µg of total RNA.
We quantified the expression of AR, AROM and ORα relative to glyceraldehyde-3-phosphate dehydrogenase (GAPDH) using the 2−ΔΔCt method  for each brain area, except that we did not quantify ORα PTR for females. All reactions were run in duplicate using SYBR green low ROX on a Stratagene MX3000P Real-Time PCR System (Agilent Technologies, Santa Clara, CA). The relative amount of transcript is reported as the fold difference relative to a pooled standard, normalized via the amount of GAPDH present in each sample. The pooled cDNA standard came from junco neural tissue collected during pilot work. Each qPCR reaction included 2.5 µl of cDNA diluted 1 : 10. Each reaction was optimized for high amplification efficiency (93.3–116.4%). We controlled for slightly unequal efficiencies using post-hoc corrections with MxPro software (v. 4.10, Agilent). Primer concentrations for each gene were 0.3 µM in a total volume of 25 µl. Primers were based on sequences from zebra finch (Taeniopygia guttata; see electronic supplementary material, table S1) and we confirmed high sequence identity in these regions (95–98%) using junco transcriptome sequences . All reactions used the following thermal profile: 10 min at 95°C, followed by 40 cycles of 30 s at 95°C, 1 min at 60°C and 30 s at 70°C, with a final dissociation phase to confirm product specificity (1 min at 95°C, 30 s at 55°C and 30 s at 95°C).
(c) Hormone assays
We quantified circulating levels of T from trunk blood. Each blood sample was treated with 50 µl of water-based heparin sodium salt solution, stored on ice and centrifuged before plasma was drawn off the top (241 ± 16 µl, range: 90–410 µl). Plasma T concentrations were quantified using a commercially available enzyme immunoassay kit (Assay Design 901–065) that has already been validated in this species . Briefly, we added tritiated T to each sample, and extracted plasma twice with diethyl ether. Extraction efficiencies were used to account for incomplete recoveries (efficiencies: 90.6 ± 0.5%). We used a logistic standard curve and curve-fitting program (Microplate Manager, Bio-Rad Laboratories) to obtain T concentrations, and we mathematically corrected for the additional volume from the water-based heparin solution. Samples were distributed over two plates, with intra-assay variability of 2.3 and 6.3 per cent, and inter-plate variability of 19.7 per cent. A plate correction factor was used to correct for inter-plate variation, based upon standards that were distributed over both plates multiple times .
(d) Statistical analyses
All statistical methods were performed using JMP v. 9.0.2 (SAS Institute, Cary, NC). We checked for normality using Shapiro–Wilk tests, and we transformed plasma T levels (natural log), flyovers (square-root) and transcript abundance (log2) to achieve normality. Log2-fold changes in gene expression include both positive and negative values, where negative values indicate that the focal sample expressed less transcript than the pooled standard. No transformation normalized the female response variable ‘time spent responding’, so all analyses with this behavioural measure used non-parametric Spearman correlations.
Because all individuals were rapidly captured following a 6-min intrusion and the time from initial disturbance to capture was short (10.5 ± 0.5 min), we assumed that measures of gene expression reflect individual differences that were present prior to the intrusion as opposed to plastic transcriptomic responses to the behavioural assay. However, to test this assumption, we used Pearson correlations to ask whether latency to sacrifice was related to transcript abundance in any gene or brain area.
Next, we tested for sex differences in behaviour and endocrine parameters (i.e. circulating T and transcript abundance in all brain areas) using unpaired t-tests, and we tested for correlations among the behavioural response variables (flyovers, latency to 5 and 1 m, songs and time spent responding). Other than latencies to 1 and 5 m, which are inherently correlated, we found no significant correlations among behavioural response measures (n = 32 individuals, p > 0.15), so we opted to treat each behavioural measure separately.
We also checked for significant relationships among T levels and measures of gene expression using Pearson correlations. Because many variables were strongly correlated, we used a two-step approach to constructing generalized linear models to measure the relationship between endocrine parameters and behavioural responses. First, we used partial least-squares regression as an exploratory tool to identify the salient variables predicting each behaviour. This approach examines variance explained by each dependent variable separately, and it is useful when multi-collinearity prevents the use of standard least-squares analyses . We used cross-validation to determine the optimal number of latent variables based on the lowest r.m.s.e. We then eliminated variables with low scores for variable importance in the projection (VIP < 0.8) and low model coefficients (absolute value < 0.12) because these variables are unlikely to explain significant variance in the dependent variable . These variables were then used in generalized linear models (GLMs, with normal error distribution and link function) to predict each behavioural output from each endocrine parameter on the narrowed list. For aggressive behaviours that were shared between the sexes (i.e. distance measures and flyovers), we included sex and a sex × neuroendocrine variable interaction to directly test whether mechanisms of aggression differed between the sexes. We then removed the main effect of sex and the interaction with sex if these variables were not significant on their own.
We observed strong and significant relationships between several neuroendocrine variables and three common measures of aggression: territorial song (in males), time responding (in females) and flyovers (in both sexes, measured in terms of flights directly over the live decoy). In support of our assumption that the short behavioural assay did not cause appreciable changes in transcription prior to sacrifice, we did not observe any significant correlations between latency to sacrifice and transcript abundance for any gene in any brain area (all |r| < 0.18, all p > 0.36).
(a) Relationships among endocrine parameters
As expected, males showed significantly higher circulating T than females (t = 7.63, d.f. = 29, p < 0.0001); so we compared T levels with measures of sensitivity separately for each sex. T levels were related to only one measure of neural sensitivity to steroids, which was the amount of AROM mRNA hypo in males only (r = −0.72, n = 16, p = 0.016; for all other correlations, p > 0.16). The sexes did not differ in behavioural measures (t < 1.30, d.f. = 30, p > 0.20), with the exception of behaviours that are specific to one sex (e.g. songs), and the sexes did not differ in any measure of transcript abundance (t < 1.23, d.f.: 26–30, p > 0.22). Thus, we considered the sexes together when correlating transcript abundance within and among brain areas, though results are qualitatively similar if the sexes are treated separately. There were no significant relationships between brain areas in AR, AROM or ORα gene expression (all |r| < 0.25, d.f.: 26–30, p > 0.17). However, the abundance of AR, AROM and ORα were positively correlated within each brain area (AR v. AROM hypo: r = 0.64, p < 0.0001, PTR: r = 0.62, p < 0.0001, VmT: r = 0.31, p = 0.081, n = 32; AR v. ORα hypo: r = 0.69, p < 0.0001, n = 29, PTR: r = 0.54, p < 0.034, n = 15, VmT: r = 0.54, p = 0.0018, n = 31; AROM v. ORα hypo: r = 0.73, p < 0.0001, n = 29, PTR: r = 0.88, p < 0.0001, n = 15, VmT: r = 0.54, p = 0.0016, n = 31).
(b) Flyovers in males and females in relation to endocrine parameters
Partial least-squares analysis pointed to the relative abundance of AR, AROM and ORα transcript in VmT as potentially important predictors of the number of flyovers directed at the intruder; neither plasma T nor any other measures of transcript abundance were identified as explaining meaningful variance in the number of flyovers (see electronic supplementary material, table S2). GLMs revealed that more aggressive individuals expressed greater AR, AROM and ORα mRNA in VmT, and there were no sex differences in these patterns of covariation (AR GLM with sex: χ2 = 9.19, d.f. = 28, p = 0.027, AR: χ2 = 7.70, p = 0.0055, sex: χ2 = 0.30, p = 0.58, sex × AR: χ2 = 0.46, p = 0.46; AROM GLM with sex: χ2 = 5.10, d.f. = 28, p = 0.16, AROM: χ2 = 3.61, p = 0.057, sex: χ2 = 0.98, p = 0.32, sex × AR: χ2 = 0.02, p = 0.88; ORα GLM with sex: χ2 = 7.11, d.f. = 27, p = 0.069, ORα: χ2 = 5.80, p = 0.016, sex: χ2 = 0.99, p = 0.31, sex × ORα: χ2 = 0.01, p = 0.93; figure 1). Abundance of AR, AROM and OR transcript significantly predicted variance in the number of flyovers among individuals after removing the non-significant effect of sex (AR GLM: χ2 = 8.26, d.f. = 30, p = 0.004; AROM GLM: χ2 = 4.09, d.f. = 30, p = 0.043; ORα GLM: χ2 = 6.11, d.f. = 29, p = 0.0135.
(c) Songs in males in relation to endocrine parameters
Partial least-squares analysis revealed several endocrine variables with potential relationships with the number of songs: circulating T, AR, AROM and ORα in hypo, and OR and AROM in PTR (see electronic supplementary material, table S3). GLMs demonstrated that males singing more songs expressed more AROM and ORα mRNA in PTR (AROM: χ2 = 5.09, d.f. = 14, p = 0.024, OR: χ2 = 5.06, d.f. = 13, p = 0.024), less AROM, AR and ORα in hypo (AROM: χ2 = 9.30, d.f. = 14, p = 0.0023, AR: χ2 = 4.98, d.f. = 14, p = 0.026, OR: χ2 = 5.09, d.f. = 13, p = 0.024), and have higher circulating T (χ2 = 4.75, d.f. = 13, p = 0.029, figure 2).
(d) Time spent responding in females, in relation to endocrine parameters
Non-parametric correlations revealed no detectable relationships between the amount of time females spent responding to the intruder instead of incubating and any endocrine parameter, except for AR mRNA in hypo. Females that spent more time responding expressed marginally less AR mRNA (ρ = −0.47, n = 16, p = 0.066; figure 3).
Many fitness-relevant phenotypes are mediated by T [1–3], and it has been widely suggested that individual differences in hormone-mediated behaviour, particularly those that cannot be explained by circulating T, may instead relate to individual differences in sensitivity to hormones in the brain, such as the abundance of AR, AROM and ORα [9,13,14,23,24]. Despite the logical appeal of this explanation, few studies have demonstrated that individual differences in behaviour correlate with individual differences in measures of neural sensitivity to sex steroids [27–29]. We found that free-living songbirds demonstrate natural individual differences in the abundance of AR, AROM and ORα mRNA in behaviourally relevant brain areas, and that transcript abundance covaries with individual differences in aggressiveness in both males and females. While we cannot yet address whether these message-level patterns translate directly to protein abundance in the brain, our results provide strong correlational evidence that natural, individual differences in gene expression at three major neural targets of sex steroids account for natural variation in aggressive behaviour. Thus our results are consistent with the long-held assertion that behavioural differences may be the result of individual differences in sensitivity to hormones . The observation that androgen- and oestrogen-dependent pathways show largely similar patterns of covariation with behaviour provides a unified sex-steroid-mediated link between aggression and hormone sensitivity in the brain.
Circulating T was also related to aggressive behaviour, but only in males, as revealed by the positive covariation between song and T. This finding confirms previous reports from this species showing that variation in hormone signal predicts variation in aggression . Our observation that female T levels did not covary with any measure of female aggressiveness is consistent with the hypothesis that the regulation of aggression via differences in sensitivity to sex steroids, instead of variation in T levels themselves, may be an adaptive route to the expression of female aggressive behaviour while avoiding the costs of systemically high T [6,37]. These findings have important implications for basic mechanisms of behaviour and for behavioural evolution.
(a) Evolutionary implications of functional individual difference in neural sensitivity
Our approach links aggressive behaviour and three measures of neural sensitivity to sex steroids in naturally varying wild animals. This approach is complementary to the large body of research that seeks to identify the neural mechanisms of aggression by comparing groups of animals that differ in aggression or neural sensitivity to sex steroid [18,21,38,61,62], or by directly manipulating AR, AROM or ORα and aggression [17,19,20,63,64]. Because these studies typically did not quantify or assign functional significance to individual variability in the abundance of receptors or mRNA [9,23,65,66], the implications of these studies have not been fully integrated into an evolutionary framework. To draw inferences about potential sources of phenotypic variation and evolution in nature, a different approach, such as the one used here, is needed to estimate evolutionary parameters [30,31] (i.e. whether there is natural covariation between behaviour and physiology on which selection may act).
This study demonstrates behavioural relevance to natural individual variation in gene expression for AR, AROM and ORα in areas of the brain known to mediate social behaviour (song control nuclei, avian medial amygdala and hypothalamic nuclei), a key step in the cohesive synthesis of evolutionary biology and behavioural neuroendocrinology . This observed covariation between behaviour and physiology constitutes a natural performance gradient, linking behaviour and physiology . The next step in the synthesis will be to relate performance to selection, by asking whether there is covariation among neural sensitivity to sex steroids, behaviour and Darwinian fitness. Together, performance and selection gradients may reveal how selection leads to behavioural evolution via underlying physiology. In males, it is well established that greater same-sex aggression can directly determine access to territories or mates, and females also may benefit from increased aggressiveness in the context of female–female competition for mates, territories or other resources [32,33]. Our results suggest that selection could shape the evolution of aggression via changes in the expression of AR, AROM or ORα in both males and females, to some degree independently of circulating levels of T.
Past research comparing groups of animals that vary in aggression would have predicted a positive association between aggression and the abundance of AR, AROM or ORα among individuals (summarized above), which is similar to the patterns we found in tissues containing the song control nuclei and medial amygdala. The positive correlations between male song and gene expression along oestrogenic pathways in telencephalic tissues (figure 2b,c) and the positive correlations between flyovers and sex steroid sensitivity in amygdalar tissues in both sexes (figure 1a–c) are entirely consistent with research linking these nuclei with the expression of aggression in birds and mammals [43,51,67]. However, previous findings would not have predicted the surprising results found in the hypothalamus, where greater song output was associated with less AR, AROM and ORα transcript (figure 2d–f).
The negative correlation between circulating T levels and the abundance of AROM mRNA in the hypothalamus in males suggests the possibility of auto-regulatory processes. The hypothalamus includes many heterogeneous cell populations, which collectively integrate various internal and external stimuli, while mediating many social and reproductive behaviours, homeostatic functions and biological rhythms . High T can have a suppressive effect on AR mRNA in the brain , and these hypothalamus-wide patterns of covariation may reflect a large-scale buffering of this brain area from the effects of sex steroids for reasons that may be correlated with, though not directly related to, aggressive behaviour.
Previous findings also would have not predicted the observation that females with more AR mRNA in the hypothalamus spent more time incubating instead of responding to an intruder (figure 3), but this observation is also consistent with the view that greater sensitivity to T in the hypothalamus may relate to the effect of sex steroids on other aspects of behaviour or physiology, such as parental care or reproduction. While additional experiments at a finer spatial scale are needed to pursue these details, the patterns we observed provide a critical link between behavioural variation and individual differences in gene expression in the brain.
(b) Does variation reflect plasticity or stable individual differences?
If individual differences in gene expression relate to behaviour, a natural question is whether these differences reflect an immediate response to the experimental conditions or whether they instead reflect variation that was present prior to the simulated intrusion. Because we found no relationship between latency to sacrifice and any measures of gene expression, our results suggest that individual differences in transcript abundance for AR, AROM and ORα were not directly attributable to rapid changes in transcription in the minutes preceding sacrifice. While the non-genomic actions of sex steroids can occur quite rapidly (e.g. less than 15 min ), the time frame in which these animals were collected is unlikely to involve socially induced changes in the expression of these genes because appreciable transcription of all but immediate early genes probably requires more time (i.e. after neuronal activation and the assembly of upstream transcription factors) . Our findings are also consistent with a study that captured wild birds after roughly 30 min of aggressive interactions, which also reports no rapid effect on some of these same genes .
It is currently difficult to speculate on the degree to which the individual differences in gene expression reported here relate to heritable individual differences. On the one hand, levels of gene expression may have a heritable component [21,46,72,73]. On the other hand, environmental and social conditions can affect transcription as well . An exciting possibility for future research is to determine whether the genes that are up- or downregulated in response to environmental stimuli may be the same genes that contribute to phenotypic evolution (e.g. genetic assimilation ). Given the integral role that hormones play in translating environmental stimuli into a phenotypic effect, hormone-metabolizing enzymes and hormone receptors may be good candidates for this process.
(c) Sex similarities and differences in mechanisms of aggression
One of the most intriguing patterns revealed by this study is that neural gene expression for AR, AROM and ORα similarly predict aggression in both sexes, but circulating T levels predict aggressiveness only in males. Comparative studies demonstrate that circulating T levels may reflect a compromise between the sexes [6,35], where selection favours higher T in males than in females. Based on genetic correlations and coevolution between the sexes, we thus might expect to see T-dependent aggressive behaviour in both sexes. Sexually biased or hormone-induced gene expression is thought to provide a solution to the problem posed by different selective pressures on the two sexes , and our results suggest a related solution in which aggression in both sexes depends upon gene expression for the molecules that initiate the genomic effects of sex steroids in the brain. Data presented here indicate that some measures of aggression can vary among individuals independently of variation in circulating T in both sexes (though singing in males does correlate with T). Behavioural independence from circulating T may permit adaptive expression of aggression that can side-step the costs of high systemic T to some degree [1,34,37]. Similar patterns appear to apply where high levels of aggression persist in the non-breeding season, despite low circulating T [38,39], and these mechanisms of aggression that relate to neural sensitivity to sex steroids may allow individuals of both sexes to express aggression somewhat independently of circulating T.
Natural variation in aggressive behaviour in males and females mapped onto individual differences in gene expression for three key molecules that facilitate the effects of T in neural target tissues, and these patterns of covariation appear to be largely shared between androgen- and oestrogen-mediated pathways. Furthermore, male aggressive behaviour, as reflected in song rate, also related to individual variation in circulating T levels, demonstrating that individual differences in hormone signal and sensitivity to signal may both account for the individually variable behaviours observed in nature. Hormones and gene expression are thought to lie at the intersection of the proximate and ultimate [76,77], and our results reveal that there is ample variation in hormone signal and in gene expression on which selection may act to affect aggressiveness. These data therefore establish a prerequisite for the evolution of T-mediated phenotypes via changes in localized gene expression for the key molecules that process sex steroids, and they suggest that phenotypic evolution can occur with some degree of independence from circulating T levels [1,34].
All procedures used in this study were approved by the Bloomington Institutional Animal Care and Use Committee. Many thanks to the Mountain Lake Biological Station for access to field sites; to D. G. Reichard for playback audio files; to N. Tonge, A. Mirzatoni, R. Stewart, M. P. Peterson and the Center for the Integrative Study for Animal Behaviour for assistance in the field and/or laboratory; and to members of the Ketterson and Goodson labs for feedback on earlier versions of the manuscript. This research was supported by an NIH NRSA fellowship to K.A.R. (T32HD049336 and F32HD068222), Indiana Academy of Sciences grant to K.A.R. and E.D.K., NSF IOS-0820055 to E.D.K., NIMH 061994 and NSF IBN-0646459 to B.A.S., and NSF DDIG (IOS-0909834) to C.M.B.B.
- Received February 27, 2012.
- Accepted May 16, 2012.
- This journal is © 2012 The Royal Society