Email updates

Keep up to date with the latest news and content from BMC Biology and BioMed Central.

Journal App

google play app store
Open Access Research article

Quality prevails over identity in the sexually selected vocalisations of an ageing mammal

Elodie Briefer1, Elisabetta Vannoni23 and Alan G McElligott1*

Author Affiliations

1 Queen Mary University of London, School of Biological and Chemical Sciences, Mile End Road, London E1 4NS, UK

2 Zoologisches Institut, Universität Zürich, Switzerland

3 Current address: Anatomisches Institut, Universität Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland

For all author emails, please log on.

BMC Biology 2010, 8:35  doi:10.1186/1741-7007-8-35


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1741-7007/8/35


Received:13 October 2009
Accepted:9 April 2010
Published:9 April 2010

© 2010 Briefer et al; licensee 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.

Abstract

Background

Male sexually selected vocalisations generally contain both individuality and quality cues that are crucial in intra- as well as inter-sexual communication. As individuality is a fixed feature whereas male phenotypic quality changes with age, individuality and quality cues may be subjected to different selection pressures over time. Individuality (for example, morphology of the vocal apparatus) and quality (for example, body size and dominance status) can both affect the vocal production mechanism, inducing the same components of vocalisations to convey both kinds of information. In this case, do quality-related changes to the acoustic structure of calls induce a modification of vocal cues to identity from year to year? We investigated this question in fallow deer (Dama dama), in which some acoustic parameters of vocalisations (groans) code for both individuality and quality.

Results

We carried out a longitudinal analysis of groan individuality, examining the effects of age and dominance rank on the acoustic structure of groans of the same males recorded during consecutive years. We found both age- and rank-related changes to groans; the minimum values of the highest formant frequencies and the fundamental frequency increased with the age of males and they decreased when males became more dominant. Both age- and rank-related acoustic parameters contributed to individuality. Male quality changed with age, inducing a change in quality-related parameters and thus, a modification of vocal cues to male individuality between years.

Conclusions

The encoding of individuality and quality information in the same components of vocalisations induces a tradeoff between these two kinds of signals over time. Fallow deer vocalisations are honest signals of quality that are not fixed over time but are modified dynamically according to male quality. As they are more reliable cues to quality than to individuality, they may not be used by conspecifics to recognize a given male from one year to another, but potentially used by both sexes to assess male quality during each breeding season.

Background

Male vocalisations are a key component of sexual selection in many species. They can potentially provide receivers (both competitors and potential mates), with important characteristics about senders [1]. Two biologically significant characteristics of the caller that are likely to be encoded in sexually selected calls are individual identity and aspects of male quality [2-5].

Vocal cues to individuality generally result from inter-individual differences in the vocal production anatomy or in the way it is operated by individuals [6]. These cues are distinctive characteristics that allow senders to be recognised by their conspecifics [7,8]. Recognition plays a crucial role in a variety of social contexts. During the breeding season, males may benefit from signalling identity in the context of intrasexual communication through decreased aggressive competition over status in dominance hierarchies, and in the context of intersexual communication when being recognizable improves mating success [6,9,10].

Quality cues in vocalisations are usually a consequence of physical or physiological constraints that affect sound production [11]. Quality signals inform receivers about the overall phenotypic and genetic constitution of the sender [7]. Therefore, such signals can be used by males to broadcast their competitive abilities to rivals and females [1]. Mammal vocalisations, including human voice, have been shown to encode various indicators of quality such as body size, weight, age, dominance rank and fatigue [3,4,12-15].

For a given individual, identity is a fixed feature, whereas some of its phenotypic quality traits can change between years and even within years according to external or internal factors (life-history evolution) [16]. Therefore signals of individuality may be subjected to different selection pressures over time and have different properties compared to signals of quality [2,10]. Tradeoffs may exist between signals of identity that should be stable over time to allow individual recognition both within a year and between years, and signals of quality that should change as a function of the phenotypic quality of the individual in order to convey reliable information (honest signals). These tradeoffs could be partially resolved by using different features that are free to vary independently of one another to signal individuality and quality (segregation of information) [17].

According to the source-filter theory of voice production [18], mammal vocalisations are generated by vibrations of the vocal folds (source) and are subsequently filtered by the supralaryngeal vocal tract (filter). The source determines the fundamental frequency of the call (hereafter F0), whereas the filter shapes the source signal by selectively amplifying certain frequencies and dampening out others. This filtering mechanism produces peaks called formants [18,19]. The source and filter are usually independent processes as they can vary independently of each other [20]. F0 and formant frequencies can thus contain independent information and could be good candidates to encode identity and quality cues separately. However, acoustic cues to quality and individuality seem in several mammalian species to be linked to both source and filter characteristics [21,22].

When identity and quality cues are not conveyed in different and independent components of the vocalisations, quality-related changes to the acoustic structure of calls may induce modifications of the vocal cues to identity. We investigated this hypothesis in fallow deer, in which some acoustic parameters of groans code for both individuality and quality (body size and dominance status) [14,21].

The fallow deer (Dama dama) is a highly polygynous, size-dimorphic and long-lived species in which males only vocalise for a short discrete period each year. As they age, their phenotypic quality changes. They usually reach prime-age (in terms of mating) at six to seven years old, with survival and reproduction probabilities generally declining thereafter [23]. Groans are sexually selected calls mainly produced by mature males (≥ four years old). They are directed at other males, and are also important for mate attraction [24,25]. Fallow bucks groan only during the breeding season (late September to early November in the northern hemisphere), and are silent for the rest of the year. The breeding season corresponds to an intense period of competition amongst males for access to females and of vocal activity during which males may call at rates of more than 3,000 groans per hour. During this period, body condition of socially mature males declines due to increased energy expenditure and decreased foraging time [26]. Strong differences in phenotypic quality and mating success exist between males; those that invest more in display have higher mating and survival rates without losing more body condition compared to other males or suffering any long-term consequences (in terms of decreased survival or reduction in fecundity) [23,27].

We investigated if the vocal cues to male identity changed between years along with changes in dominance status and age. We first tested whether the phenotypic quality of males changed between years and the impact of these quality changes on mating success. We then examined age- and rank-related changes to the acoustic parameters of groans. Finally, we determined if these changes were related to modifications in the vocal cues to individuality. We predicted that if male quality changes between years, honest vocal quality traits should change accordingly. Such modifications should have important implications for signaling identity and quality in sexual selection.

Methods

Study site and population

The study was carried out on a herd of European fallow deer in Phoenix Park (53° 22' N, 6° 21' W), Dublin, Ireland. The age and identity of males used in this study were known because they had been tagged as fawns, as part of the management of the deer by the park authorities [28].

Behavioural observations, mating success and dominance relationships

We conducted behavioural observations during the breeding season in 2000, 2002, 2003 and 2004. All-event recordings of agonistic interactions and matings were carried out every day from dawn to dusk (circa 11 hours per day). There were usually 7 to 13 observers in the field at all times and the coverage of animals was maximised. The measure of mating success for each male (years 2000, 2002 and 2003) was based on the number of observed copulations during the breeding season, which provides a very good estimator of reproductive success in fallow deer [29].

The outcomes of the agonistic interactions were used to calculate the dominance rank of each male by applying the David's score method [30]. This method is appropriate when interactions are recorded over a long period of time, because it takes into account repeated interactions between dyad members that may determine win/loss asymmetries [31]. Dominance ranks were calculated for males that interacted with at least 10% of other mature males. The rank value obtained for each male was then divided by the number of males present each year (55.5 ± 1.5, N = four years) to standardize this variable (Rank, one measure per male per year). The breeding season can be divided into two main periods according to the availability of mating opportunities: the prerut, which begins when all males have cleaned the velvet from their antlers and ends on the day before the first mating, and the rut, which is the period between the days on which the first and last matings occur [26]. Rank values available for males in 2002 and 2004 were prerut values whereas those available for males in 2000 and 2003 were rut values. Dominance ranks of fallow bucks can be measured soon after they regrow their antlers each year and before they become vocal. Male rank is thus well established before the rut so that prerut and rut rank values are highly correlated [26,32].

Recordings and selection of groans

Recordings were made using a directional microphone (Sennheiser MKH-70; Old Lyme, CT, USA)connected to a digital audio tape recorder (Sony TCD-D100; Tokyo, Japan). Groans were recorded between dawn and sunset at a distance of 10 to 50 m from the vocalising animal. Vocalisations were imported into a computer using Avisoft-SASLab Pro v.4.38 (Avisoft Bioacoustics, Berlin, Germany) [33] at a sampling rate of 22.05 KHz and saved in WAV format, and at 16-bit amplitude resolution [34]. The recordings that did not contain energy above 8 KHz (as visualised on a spectrogram) were downsampled to 16 KHz for a better frequency resolution. Narrowband spectrograms of groans (FFT method, window length = 0.03 s, time step = 1,000, frequency step = 250, frequency resolution = 20 Hz, Gaussian window shape, dynamic range = 50 dB) were edited using Praat v.5.0.47 DSP Package (P. Boersma and D. Weenink, University of Amsterdam, The Netherlands) [35]. Vocalisations with high levels of background noise were not considered for analyses.

During the rut, fallow bucks feed very little, lose approximately 26% of their body weight, and the acoustic structure of groans is affected by fatigue [15,27]. We therefore analysed recordings taken between 9 October and 23 October when only a small proportion (15% or less) of the total number of matings had usually occurred, and when the majority of agonistic interactions among males were non-contact displacements [15,36]. For each male, selected groans had been recorded on the same day or on different days. Vocalisations produced within the same bout are more likely to be acoustically homogeneous than vocalisations emitted in different bouts [37]. To minimize this problem, the majority (84%) of the groans were randomly selected from different bouts and the remaining groans (16%) belonged to long vocalisation bouts but were never consecutive. We included in our analysis males between five and eight years of age because they had reached their maximum size and mass, and because this range of ages includes the males that are socially mature and account for the vast majority of vocalisations and matings [23,28,32].

Acoustic analysis

We analysed the vocalisations of 14 males recorded during either two (11 males) or three (3 males) breeding seasons (2000, 2002, 2003 and 2004). Groans are short and low-pitched vocalisations characterized by a pulse-train structure, which is visible on both the envelope and the spectrogram of the signal (Figure 1). The pulses represent the vibrations of the vocal folds and determine the fundamental frequency (F0) of the call. From the envelope of the call (amplitude vs. time), we measured mean, minimum and maximum F0 (F0mean, F0min and F0max), the number of pulses (Pulses), the duration of the groan (Duration) and a measure of cycle-to-cycle frequency variation (F0 perturbation), known as Jitter [19]. F0 is equivalent to the inverse of interpulse interval and this can be measured as the distance between consecutive pulse onsets [21]. Distances between pulses were measured using pitch analysis (Sound: To manipulation command) in Praat. We then calculated the values of F0 along the groan and averaged these values to obtain F0mean. Because all groans showed at least a modest frequency inflection, F0min and F0max were also included in the analyses. We also measured Pulses, and calculated Duration as the distance between the onset of the first pulse and the end of the last one. Finally, we quantified the variation of F0 along the call by calculating Jitter. Jitter has been used as a measure of voice quality in mammals [38], and in a similar manner to other source-related features, this parameter may contribute to vocal distinctiveness [21,39]. In our analysis, Jitter was calculated by dividing the average absolute difference between consecutive frequencies by the mean F0 per groan (peak-picking method) [40,41]. From the spectrum of the call (frequency vs. amplitude), we measured the frequency of maximum amplitude (Frmax) using peak detection with the package Seewave v.1.5.4 in R v.2.9.0 (R Development Core Team, Vienna, Austria) [42]. Frmax has been shown to signal individual identity in mammals [43-45].

thumbnailFigure 1. Age- and rank-related changes in groans. Spectrograms (above) and oscillograms (below) of three groans produced by the same male at six, seven and eight years old. Pulses are visible both on the spectrograms and oscillograms and are indicated as black vertical lines. Formants are visible as horizontal frequency bands on the spectrograms. Age-related changes to F0min, and rank-related changes to F6min are shown: F0min was measured for each groan as the inverse of the maximum interpulse interval (grey rectangle) on the oscillogram and increases with age (grey rectangle decreases); F6min is indicated with a black arrow and increases when the male becomes lower-ranking.

Finally, the minimum frequencies of the first six formants (F1min-F6min) were estimated using Linear Predictive Coding analysis (LPC; Sound: To Formant (burg) command) in Praat. Formants are evident as horizontal frequency bands on the spectrogram of the groan (frequency vs. time and amplitude indicated as a colour scale; Figure 1). We compared the outputs of the LPC analysis with visual inspections of spectrograms to check if Praat was accurately tracking the formants. By performing a single LPC analysis on each groan, higher formants (F4 to F6) were better detected and therefore more accurately measured than lower formants (F1 to F3). We therefore conducted a double or a triple LPC analysis on each groan in order to get the best estimations of all formants. We first carried out an LPC analysis (time step = 0.01 s, maximum number of formants = three to five, maximum formant = 500 to 800 Hz, window length 0.04 to 0.07 s) to measure the frequencies of the first three formants (F1 to F3). Then we performed a second LPC analysis (time step = 0.01 s, maximum number of formants = four to seven maximum formant = 1,400 to 2,300 Hz, window length 0.04 to 0.07 s) to estimate the next three formant frequencies (F4 to F6). When the sixth formant was not detected by the second LPC analysis, we conducted a third LPC analysis (time step = 0.01 s, maximum number of formants = four to seven maximum formant = 1,900 to 2,300 Hz, window length 0.04 to 0.07 s). The decrease in formant frequencies along the groan reflects the elongation of the vocal tract occurring during vocalisations [46]. We calculated F1min to F6min from each groan by averaging the values over the last part of the call when formants become flat, after deleting spurious values and correcting for octave jumps. Minimum formant frequencies occur when the larynx is pulled down to the maximum extent, and therefore provide reliable information about the caller's body size [15,46]. We estimated the minimum spacing of the formants, known as minimum formant dispersion (Dfmin) [12] using the method described by [4].

Statistical analyses

Data on dominance ranks and mating success were not available for all males whose vocalisations had been analysed. As a result, sample sizes (number of groans and number of males) varied among analyses including Age, Rank and Mating Success.

Age and dominance rank effects on mating success and on acoustic parameters of groans

We investigated the effect of Age on Rank and on Mating Success, the effect of Rank on Mating Success, and Age and Rank effects on acoustic parameters of groans, using Generalized Linear Mixed Models (GLMMs) fitted with Restricted Estimate Maximum Likelihood (REML, lme function in R) [47]. Based on scatterplots showing the relationship between the parameters used in the different GLMMs, we fitted as fixed effects linear, quadratic or log terms. To fit a quadratic term as a fixed effect, we z-transformed the variable and included all terms (linear and quadratic) in the model [48]. When more than one model was applied using different terms (linear, quadratic or log) to test a given relationship and both were significant, we assessed the best model based on the Akaike Information Criterion adjusted for small sample size (AICc) [14,47]. For a given model, the value of AICc is a measure of the information lost when the model is used to explain a particular pattern. Therefore, the model with the smallest AICc value is estimated to best fit the data set relative to other models considered [49]. In this case, only best models (lowest AICc) are presented in the results. Q-Q plots and scatterplots of the residuals of the dependent variables were inspected visually to ensure their normal distribution. The calculation of a coefficient of determination R2 for GLMM is not obvious because of the presence of random effects. We thus estimated R2 following [50] to describe the way models fitted the observed data, as follows: R2 = 1 - exp (-2/n (logLM - logL0)), where n is the number of observations (groans), logLM is the standard log-likelihood of the model (which include fixed and random effects) and logL0 is the standard log-likelihood of the intercept-only model.

We carried out a first series of models to measure the effect of Age on Rank and on Mating Success, and the effect of Rank on Mating Success. In these models, male individual identity was fitted as a random term to control for repeated measurements of the same individual between years. In the models assessing the effect of Age on Rank and on Mating Success, linear (Age) and quadratic (Age2) terms were fitted as fixed effects (N = 12 individuals, 25 Rank values). In the model assessing the effect of Rank on Mating Success, log term (log(Rank)) was fitted as a fixed effect (N = 10 individuals, 16 Mating Success values). As the relationship between Age and Rank was significant, instead of including the original Mating Success values in the model testing the effect of Age on Mating Success, we included the residuals extracted from the model testing the effect of log(Rank) on Mating Success, which are thus independent of Rank. In the same way, in the model testing the effect of Rank on Mating Success, we included the residuals extracted from the model testing the effect of Age on Mating Success, which are thus independent of Age.

We then carried out a second series of models to examine Age- and Rank-related changes to the acoustic parameters of groans. In these GLMMs, year of recording nested within individual identity was fitted as a random term to control for repeated measurements of the same individual each year and across years. The date of recording (9 to 23 October) was included as a covariate (linear (Date) only or both linear (Date) and quadratic (Date2) terms) following [15] to control for individual acoustic variation over the breeding season. In the models assessing the effect of age on acoustic parameters, we fitted as fixed effects a linear term (Age) for Frmax, Pulses, Duration, and parameters related to F0 (F0mean, F0min, F0max and Jitter), and both linear (Age) and quadratic (Age2) terms for parameters related to formants (F1min-F6min and Dfmin; N = 305 groans, 9.84 ± 0.09 per individual per year, 14 individuals). In the models assessing the effect of Rank on acoustic parameters, a log term (log(Rank)) for Frmax, Pulses, Duration, and parameters related to F0 (F0mean, F0min, F0max and Jitter), and both linear (Rank) and quadratic (Rank2) terms for parameters related to formants were fitted as fixed effects (F1min-F6min and Dfmin; N = 245 groans, 9.80 ± 0.12 per individual per year, 12 individuals).

Between-years modification of the vocal cues to individuality

We quantified the individual distinctiveness of groans by performing a Principal Component Analysis (PCA) followed by a Multivariate Analysis Of Variance (MANOVA) and Discriminant Function Analyses (DFAs). We used a PCA carried out on all 2000, 2002, 2003 and 2004 recordings to eliminate redundancy due to the high intercorrelation of the acoustic variables in our data and to examine clustering among variables (N = 305 groans, 9.84 ± 0.09 groans per individual per year, 14 different individuals). We retained the Principal Components (PCs) of the PCA with eigenvalues greater than 1 (Kaiser's criterion). The scores of the five extracted PCs (PC1 to PC5) were confirmed for normality (Kolmogorov-Smirnov test) and used as input variables in the MANOVA and DFAs.

We used a MANOVA, with individual identity of bucks nested within year of recording fitted as a random term, to confirm within-year statistical differences in PC1 to PC5 scores across individuals. We then used DFAs with one factor (individual) to quantify the extent to which individuals can be classified on the basis of their groans and which group of variables (which PC) account most for this classification [51]. On the basis of the discriminant functions of the DFA, each PC score (corresponding to a groan) was assigned to the appropriate individual (correct assignment) or to another individual (incorrect assignment). We cross-validated our results by performing a leave-one-out classification, which is an appropriate method for small sample sizes [52-54]. We calculated the percentage of correct classification due to chance by applying a randomization procedure. The expected level of correct assignment was averaged from DFAs performed on 1,000 randomized permutations of the data set [54,55].

We first carried out a DFA and calculated percentages of correct classification on each year separately to assess within-year individual acoustic variation (same-year correct classification rate). As only one of the males recorded in 2000 was selected for the analyses, we carried out DFAs on PC1 to PC5 scores corresponding to 2002, 2003 and 2004 recordings (three different DFAs, N = 10.00 ± 1.00 individual per year).

To investigate between-years modification of vocal cues to male individuality, we then calculated percentages of correct classification of PC scores (corresponding to groans) of one year by a DFA carried out on a previous year (cross-year correct classification rate; function predict in R v.2.9.0) [56]. This procedure allowed us to calculate correct classification rates for groans recorded in 2003 and 2004 by projecting these data onto the discriminant functions generated by the 2002 DFA (N = seven and four individuals respectively), and for groans recorded in 2004 by projecting these data onto the discriminant functions generated by the 2003 DFA (N = eight individuals). Same-year correct classification rates (calculated using leave-one-out classification) and cross-year correct classification rates were then compared using a two-tailed exact permutation test (N = 13 individuals, 19 comparisons). A significant difference between same-year and cross-year correct classification rates would indicate a between-years modification of vocal cues to male individuality, according to the parameters measured in this study.

Statistical analyses were carried out using R v.2.9.0 (R Development Core Team, Vienna, Austria). All tests were two tailed and factors were considered to have a statistically significant influence if P < 0.05. All means are given with SEs.

Results

Age and dominance rank effects on mating success and acoustic parameters of groans

There was a tendency for Age to affect Rank, with males being more dominant at six and seven years old than at five and eight years old (quadratic relationship, GLMM: F1,11 = 3.33, P = 0.09, R2 = 0.23; Figure 2). Males tended to become higher-ranking by 0.14 ± 0.05 (for example, 8.12/58 rank positions) between five and six years old (or six and seven years old depending on the males considered) and to become lower-ranking by 0.12 ± 0.03 (for example, 6.96/58 rank positions) between six and seven years old (or seven and eight years old depending on the males considered). The effect of Age on Mating Success residuals (controlled for Rank) was not significant (GLMM: linear, F1,4 = 0.01, P = 0.93, R2 = 0.0007). The log relationship between Rank and Mating Success residuals (controlled for Age) was significant (GLMM: F1,4 = 10.19, P = 0.03,, R2 = 0.41), with males being more successful when high-ranking than when low-ranking. All males having a rank value ≤0.1 (for example, ≤5/58) obtained ≥14 matings. The most important difference was found for a male that obtained 46 matings when his rank was of 0.07 (4/54) and only six matings when his rank was of 0.26 (15/58). Therefore, male age influences their dominance rank, and higher-ranking males (six to seven years old) gain the most matings.

thumbnailFigure 2. Relationship between the age and dominance rank of males. Lower values of dominance rank indicate higher-ranking males. Males tend to be higher-ranking at six and seven than at five and eight years old (quadratic relationship, P = 0.09). Dots represent mean values per year for each male and lines show repeated measures of the same individual across years.

The results of the GLMMs investigating Age- and Rank-related changes to the acoustic parameters are shown in Table 1 (see also Additional files 1, 2 and 3). There was a significant positive linear relationship between Age and all measures of the F0 contour (F0max, F0mean and F0min), with males producing higher-pitched groans when older. This relationship was stronger for F0min (mean F0min: five years old, N = three, 19.39 ± 0.53 Hz; six years old, N = 12, 20.71 ± 0.50 Hz; seven years old, N = 12, 22.05 ± 0.54; eight years old, N = four, 23.26 ± 0.98; Figures 1 and 3). Age was also significantly related to F4min, F6min and Dfmin, with males having higher minimum formant frequencies and higher Dfmin (corresponding to shorter apparent vocal tracts) at five and eight years old than at six and seven years old (quadratic relationship, Figure 4). The effects of Age on Jitter, Duration, Pulses, Frmax, F1min-F3min and F5min were not significant (Table 1). Rank had a significant effect on F0max, with males having lower-pitched groans when higher-ranking. Rank also had a marginally significant effect on the highest formant frequencies (F5min, quadratic relationship, P = 0.07; F6min, linear relationship, P = 0.05), with males having lower minimum formant frequencies when higher-ranking (Figures 1 and 5). The effects of Rank on F0mean, F0min, Jitter, Duration, Pulses, Frmax, F1min-F4min and Dfmin were not significant (Table 1). Thus, the age and dominance ranks of males shape their groans, with both age and dominance rank affecting some measures of the fundamental frequency and the minimum frequencies of the highest formants.

Additional file 1. Age: 6, rank: 2/58. Groans produced by a male when he was six years old and high ranking (2/58). The fundamental frequency and minimum frequencies of the highest formants are low.

Format: WAV Size: 91KB Download fileOpen Data

Additional file 2. Age: 7, rank: 11/58. Groans produced by the same male (as in additional file 1) when he was seven years old and lower ranking (11/58). The fundamental frequency and minimum frequencies of the highest formants are higher than for the groans in Additional file 1.

Format: WAV Size: 113KB Download fileOpen Data

Additional file 3. Age: 8, rank: 24/58. Groans produced by the same male (as in Additional files 1 and 2) when he was eight years old and his rank had declined further (24/58). The fundamental frequency and minimum frequencies of the highest formants are higher than for the groans in Additional files 1 and 2.

Format: WAV Size: 134KB Download fileOpen Data

Table 1. Age- and rank-related changes in the acoustic parameters of groans.

thumbnailFigure 3. Relationship between the age of males and the minimum fundamental frequency of groans. The minimum fundamental frequency increases as males get older (linear relationship). Dots represent mean values per year for each male and lines show repeated measures of the same individual across years.

thumbnailFigure 4. Relationship between the age of males and the minimum formant dispersion of groans. Males produce groans with higher minimum formant dispersion, indicating shorter apparent vocal tract lengths when they are five and eight years old than when they are seven and eight years old (quadratic relationship). Dots represent mean values per year for each male and lines show repeated measures of the same individual across years.

thumbnailFigure 5. Relationship between male rank and the minimum frequency of the sixth formant of groans. Lower values of dominance rank indicate higher-ranking males. Males tend to produce groans with a lower minimum frequency of the sixth formant when higher-ranking (linear relationship, marginally significant, P = 0.05). Dots represent mean values per year for each male and lines show repeated measures of the same individual across years.

Between-years modification of the vocal cues to individuality

The PCA performed on groans (N = 305 groans, N = 14 individuals, Table 2) generated five components (PC1 to PC5) that exceeded Kaiser's criterion (eigenvalues >1). These five components combined accounted for 77.5% of the variation in the original data set. Filter-related parameters were grouped in two different components that represent different aspects of the vocal tract function (PC1 and PC4); PC1 (UpperF) included higher formant frequencies (F3min-F6min and DFmin), which are mainly related to the length of the vocal tract, whereas PC4 (LowerF) consisted of lower formant frequencies (F1min-F2min), which are generally associated with the shape of the vocal tract. Two components reflected aspects of source characteristics (PC2 and PC5); PC2 (F0) represented measures of the F0 contour (F0max, F0mean and F0min), whereas PC5 (F0&Jitter) included F0min and Jitter, a measure of voice quality based on the stability of vocal folds vibrations. Temporal parameters (Pulses and Duration) were highly correlated with PC3 (TempP). Finally, Frmax was not correlated with any of the five PCs. There was a significant within-year difference among males for principal component scores (MANOVA: F30,274 = 5.23, P < 0.0001). Thus, males of a given year produce groans that differ for five sets of correlated parameters, with each set of parameters reflecting a different component of the sound production apparatus.

Table 2. Factor loadings of the measured acoustic parameters on the first five principal components.

Each of the DFAs carried out on PC1 to PC5 scores corresponding to 2002, 2003 and 2004 recordings (N = 10.00 ± 1.00 individual per year) generated five discriminant functions. For each year, the first two discriminant functions, shown in Table 3, together accounted for 77.16 ± 0.03% of the variation (N = three DFAs). Both functions were highly correlated (|r| ≥0.5) with scores from PC1 (higher formants), PC2 (F0 contour) and PC3 (Pulses and Duration). The second discriminant function was also highly correlated with PC4 (lower formants) and PC5 (Jitter and F0min), but only for the 2004 DFA. The remaining three discriminant functions together accounted for 22.84 ± 0.03% of the variation (N = three DFAs). Each year, cross-validated DFA classified 38.66 ± 0.02% of groans to the correct male (same-year correct classification rate, N = three DFAs). This correct classification rate is higher than the percentage of correct classification due to chance (10.13 ± 0.01%; permutation test: 1,000 permutations, P = 0.0009 for each year). Cross-year correct classification rates of groans (PC scores) by DFAs carried out on the preceding year (22.43 ± 0.07%, correct classification rates of 2003 groans (N = seven individuals) and 2004 groans (N = eight individuals) by 2002 and 2003 DFA respectively) were still higher than chance level (9.81 ± 0.01%, permutation test: 1,000 permutations, 2003, P = 0.0009; 2004, P = 0.02). However, the cross-year correct classification rate of 2004 groans by 2002 DFA (15.0%, N = four individuals) was not significantly higher than chance level (11.09 ± 0.70%, permutation test: 1,000 permutations, P = 0.07). Same-year correct classification rates (36.09 ± 0.05%) were significantly higher than cross-year correct classification rates (21.62 ± 0.05%), indicating a between-years modification of vocal cues to male individuality (exact permutation test, P = 0.02). These results suggest that the vocal signals of male individual identity, which are parameters related to both source (larynx) and filter (vocal tract) characteristics, change by 13% (35 to 22% correct classification rates between consecutive years) between consecutive breeding seasons.

Table 3. Contribution of the acoustic parameters to the individuality of groans.

Discussion

We investigated if the vocal cues to individuality of fallow bucks changed between years because of changes in phenotypic quality, using a longitudinal analysis with the same individuals compared at different ages. As males aged, their dominance ranks changed and rank was a good predictor of mating success. We found both age- and rank-related changes in the acoustic parameters of groans. All measures of the fundamental frequency contour (F0max, F0mean and F0min), the minimum frequencies of the highest formants (F4min, F6min) and the minimum formant dispersion (Dfmin) changed according to the age of males, with males producing groans with higher F0 and higher formant frequencies when older. Rank-related changes to the minimum frequency of the highest formants (F5min and F6min) were marginally significant, and F0max was also affected by rank changes, with males having lower F0max when higher-ranking. Within each breeding season, groans were individually distinctive, with higher formant frequencies (F3min-F6min and DFmin), F0 contour (F0max, F0mean and F0min) and temporal parameters (Pulses and Duration) being the most important cues to male individual identity. We found that these signals of individuality substantially changed between breeding seasons. By including repeated measures and thus isolating the effects of differential mortality and age-specific signalling [57-59], our study is the first to investigate how both individuality and phenotypic quality shape vocalisations over time. Our results show that fallow deer vocalisations are honest signals of quality that are not fixed but are modified dynamically according to male phenotypic quality, showing a very robust example of 'truth in quality advertising' in animal communication [11].

Are fallow bucks using calls to encode their individuality or quality? Whereas the individuality of males is a fixed feature, we showed that one aspect of their phenotypic quality (dominance status) changed with age, and this had an impact on mating success (see also [23]). However, a comparison between the effects of age and dominance rank on vocalisations (Table 1) and the vocal cues to identity (Tables 3) shows that individuality and quality signals are partially conveyed in the same components of groans. According to the acoustic parameters measured in this study, calls were individualized within each breeding season, but only 38.7% (mean N per DFA = 10.0 individuals) of groans were assigned to the correct individual (see also [21], 36.6% and 53.6% with a stepwise procedure, N = 16 individuals). In fallow bucks, both source- and filter-related parameters can vary within individuals, F0 by the action of laryngeal muscles and formants by an elongation the vocal tract through retraction of the larynx [20,46]. This dynamic may result in acoustic parameters being unreliable cues to identity [38] and could explain the poor individuality of fallow buck groans that we found. However, further analyses would be needed to know if vocal cues to identity are present in other groan parameters (for example, shimmer, spectral slope, details of formant pattern, relative formant amplitudes) [38,60,61]. Indeed, Reby et al. [62] found a higher individuality of fallow buck groans than we did (87.9%) using a network classification based on the spectrum characteristics of calls. However, this higher correct classification rate could also be due to the lower number of individuals used in that study (N = four individuals) [62]. In the present study, we also found that the structure of groans changed between years according male quality, inducing a modification of vocal cues to identity. All these results suggest that groans are more reliable signals of quality than of individuality. There is some evidence that fallow deer females attend to male quality in choosing mates [63,64]. Because male quality changes from year to year, we suggest that females should attend to quality cues more than individuality cues to select the best quality males.

We found that F0 (especially F0min), the minimum frequencies of highest formants (F4min and F6min) and the minimum formant dispersion (DFmin) changed as bucks got older, and these are therefore reliable cues to age. F0, which increased with age in our study, is negatively related to vocal fold length and positively related to vocal fold tension. The larynx is not constrained by neighbouring skeletal structures, so that vocal fold length is not particularly dependent on body size [12,20]. Moreover, males included in our study had reached their maximum size and mass, and groans selected for the analyses had been recorded early in the rut before males lost condition [15,27,32]. Therefore, the observed age-related changes to F0 cannot be explained by changes in body size.

In non-human mammals, age-related changes to vocalisations have been well studied during ontogeny [65-67], but studies investigating such changes occurring in adults as individuals get old are sparse. In baboons (Papio cynocephalus ursinus), Fischer et al. [3] found that males produced lower mean F0 when older. In red deer (Cervus elaphus), older males produce lower minimum and maximum F0 [4]. This decrease in F0 with age in red deer may be the consequence of the observed lengthening of the vocal folds throughout life, even after stags stop growing [5]. However, as we found an increase and not a decrease in F0 as males get older, this phenomenon probably does not occur in fallow deer.

The observed increase in F0 with age is more likely the result of decreased testosterone levels. In humans and closely related primates, the larynx is assumed to be a steroid receptor target organ [68-70]. Pubertal androgens cause thickening and lengthening of the vocal folds, and voice F0 is negatively related to circulating testosterone levels in adult humans [71-73]. Rolf and Fisher [74] found that testosterone levels in fallow deer peak during the rut for the first time in the fifth year of life. Moreover, the testosterone level for a seven-year-old male was lower than for a six-year-old male, suggesting a decrease in testosterone level after six years old, which could induce the increase in F0 observed in our study. An alternative hypothesis is that F0 may increase as a consequence of age-related anatomical and physiological changes in the larynx. Indeed, in humans, vocal folds atrophy, ossification of the laryngeal skeleton and increased stiffness of vocal fold tissue cause a rise in F0 with age in males [75,76].

We found that the minimum frequencies of the highest formants (F4min and F6min) and minimum formant dispersion (DFmin) were higher when males were five and eight years old, than when they were six and seven years old. This indicates a shortening of the apparent vocal tract of males as they age. Higher formant frequencies and minimum formant dispersion are positively related to vocal tract length. Unlike the larynx, the vocal tract is constrained by the bones of the skull, jaw, and spinal column, so that its length strongly depends on body size. Thus, in mammals in which the larynx occupies a standard high position, formant frequencies are usually accurate cues to body size [12,77]. However, fallow bucks and red deer stags possess a descended larynx, which is also highly mobile. The larynx can be retracted during vocalisation by the sternothyroid and sternohyoid muscles, inducing an increase in vocal tract length and consequently a lowering of formant frequencies (fallow deer [46] and red deer [78]). In those species, only the minimum formant frequencies that an individual can produce when the larynx is fully retracted toward the sternum may be an accurate cue to body size (red deer [4], fallow deer [14]). In our study, because only mature males were included, the observed age-related changes in minimum frequencies of formants cannot be caused by body size changes with age. Furthermore, our results indicate an increase in formant frequencies with age, not a decrease that could result from an increase in body size with age.

The age-related changes in formant frequencies that we found are more likely caused by a combination of body condition differences in males, and senescence of the muscles responsible for laryngeal retraction. Unlike in fallow deer, humans, elephant seals (Mirounga leonina) and baboons produce calls with lower formants and lower formant dispersion when older [3,67,75], indicating vocal tract lengthening with age. In humans, this vocal tract lengthening has been suggested to be caused by lowering of the larynx in the neck as a result of stretching of ligaments and atrophy of the strap muscles of the neck [75].

We found that males tended to produce groans with lower minimum frequencies of the highest formants (F5min and F6min) when they were higher-ranking. Red deer stags have been shown to modulate vocal tract length and thus their formant frequencies as a function of the threat (body size) posed by opponents [79]. We suggest that a similar phenomenon exists in fallow deer, with bucks modulating vocalisations in relation to their dominance each breeding season. Fallow bucks use conditional competitive strategies to decide when to interact and fight, with fights being more likely to occur between opponents with similar dominance rank [36]. Furthermore, mature males increase and immature males (≤ three years old) decrease groaning rates following playbacks of groans from mature males [80]. These results strongly suggest that fallow bucks adjust their behaviour and vocalisations in response to the presence and quality of competitors. The dominance rank of males is established before they become vocal and remains largely stable across most of the breeding season [26]. Rank-related changes in the minimum values of the highest formant frequencies observed in our study may thus result from males adjusting their vocal tract length to their dominance status each year within the range of frequencies that each of them can physically produce. By doing so, males may reduce social costs associated with agonistic interactions, by avoiding contests with potentially stronger males (badge of status) [7,10]. This hypothesis is even more likely as, contrary to red deer that typically fully retract the larynx down to the sternum during roaring [78], fallow deer usually retract the larynx on average 52% during groaning [46], which could allow more variability in formant frequencies. Therefore, the formant frequencies of fallow buck groans may be even more prone to motivation- or hormone-related modifications than red deer formant frequencies.

F0 can be modulated through the action of the vocal muscles, thereby inducing changes in vocal fold tension [20]. The observed effect of rank-related changes to F0max of groans could thus result from the same mechanism as formants, with bucks modulating groan F0 according to their dominance rank, to avoid unmatched fights. Such F0 modulation occurs in humans, with men lowering or raising voice F0 when speaking to another individual who they perceive to be more or less physically dominant, respectively, to themselves [81]. Vannoni and McElligott [14] found a negative relationship between F0min and dominance rank by investigating the relationship between acoustic parameters and dominance rank for different individuals (cross-sectional method). We did not find that result in the present study and this is probably because we investigated this relationship using a more appropriate longitudinal method, in which groans of the same males were compared when they had different dominance rank values. These two kinds of methods may lead to slightly different results due to a potential confounding effect of differential mortality and age-specific signalling in cross-sectional studies [57,58,82].

Do all these observed age and rank-related changes induce a modification of vocal cues to the identity of bucks between successive breeding seasons? We found that both age-related parameters (F0 contour) and rank-related parameters (minimum frequencies of the highest formants) contributed strongly to the individuality of groans. Consequently, although groans were individualized within each breeding season, vocal cues to male individuality changed by on average 13% between consecutive years. These observed changes have a potential impact on individual recognition by conspecifics. Because males are silent outside the breeding season and live apart from females, individual recognition of a given male's groans by conspecifics may have to be re-established every year.

Females could benefit from mating with familiar males; familiarity with a given male being an indication of his investment in reproduction [8,26,83]. In fallow deer, long-term investment in vocal activity is related to mating success, suggesting that the familiarisation of females to male vocal characteristics is important for reproduction [26]. The present study suggests that this familiarisation may be more important within than between breeding seasons. During the breeding season, high-ranking males start groaning about three weeks before the first females are in oestrous [26] and this could be a strategy to familiarise females with the new sounds of their groans each year.

Conclusions

To conclude, this study shows that individuality and quality signals can be encoded in the same components of vocalisations. This induces a tradeoff between these two kinds of signals over time, leading to one of them being more reliable than the other. Our results show that evolution favours phenotypic quality cues over individuality cues in fallow buck groans, with quality cues being more reliable than individuality cues over time.

Abbreviations

AIC: Akaike Information Criterion; DFA: Discriminant Function Analysis; Dfmin: minimum Formant Dispersion; F0: Fundamental frequency; F0mean, F0min and F0max: mean, minimum and maximum Fundamental frequency respectively; F1min-F6min: minimum frequency of the first to sixth Formant; Frmax: Frequency of maximum amplitude; GLMM: Generalized Linear Mixed Model; LPC: Linear Predictive Coding; MANOVA: Multivariate Analysis Of Variance; PCA: Principal Component Analysis; PC1-PC5: Principal Component one to five; R: coefficient of determination; RELM: Restricted Estimate Maximum Likelihood; SE: Standard Error.

Authors' contributions

EB analysed the data and wrote the manuscript. EV and AGM designed the study, collected the data and wrote the manuscript.

Acknowledgements

We are grateful to Maurus Birrer, Marco Torriani and Petra Zajec for field assistance, Rob Knell, Roger Mundry and Kamran Safi for statistical advice, and Ben Charlton, Tecumseh Fitch, Tom Reader, and the referees for helpful comments. E Briefer is funded by a Swiss National Science Foundation fellowship. We acknowledge the financial support of the Forschungskommission der Universität Zürich to AG McElligott, and the Swiss Academy of Sciences to E Vannoni during some of the data collection. We thank the staff of Phoenix Park and University College Dublin for their assistance.

References

  1. Bradbury JV, Vehrencamp SL: Principles of Animal Communication. Sinauer, Sunderland, MA; 1998. OpenURL

  2. Dale J, Lank BD, Reeve HK: Signaling individual identity versus quality: a model and case studies with ruffs, queleas, and house finches.

    Am Nat 2001, 158:75-86. PubMed Abstract | Publisher Full Text OpenURL

  3. Fischer J, Kitchen D, Seyfarth R, Cheney D: Baboon loud calls advertise male quality: acoustic features and relation to rank, age, and exhaustion.

    Behav Ecol Sociobiol 2004, 56:140-148. Publisher Full Text OpenURL

  4. Reby D, McComb K: Anatomical constraints generate honesty: acoustic cues to age and weight in the roars of red deer stags.

    Anim Behav 2003, 65:519-530. Publisher Full Text OpenURL

  5. Reby D, McComb K: Vocal communication and reproduction in deer.

    Adv Stud Behav 2003, 33:231-264. Publisher Full Text OpenURL

  6. McComb K, Reby D: Communication in terrestrial animals. In Encyclopedia of Neuroscience. Volume 2. Edited by Squire LR. Oxford: Academic Press; 2009::1167-1171. Publisher Full Text OpenURL

  7. Tibbetts EA, Curtis TR: Rearing conditions influence quality signals but not individual identity signals in Polistes wasps.

    Behav Ecol 2007, 18:602-607. Publisher Full Text OpenURL

  8. Reby D, Hewison M, Izquierdo M, Pepin D: Red deer (Cervus elaphus) hinds discriminate between the roars of their current harem-holder stag and those of neighbouring stags.

    Ethology 2001, 107:951-959. Publisher Full Text OpenURL

  9. Seyfarth RM, Cheney DL: Signalers and receivers in animal communication.

    Ann Rev Psychol 2003, 54:145-173. Publisher Full Text OpenURL

  10. Tibbetts EA, Dale J: Individual recognition: it is good to be different.

    Trends Ecol Evol 2007, 22:530-537. Publisher Full Text OpenURL

  11. Fitch WT, Hauser MD: Unpacking "honesty": vertebrate vocal production and the evolution of acoustic signals. In Acoustic Communication. Edited by Simmons AM, Fay RR, Popper AN. New York: Springer; 2002:65-137. OpenURL

  12. Fitch WT: Vocal tract length and formant frequency dispersion correlate with body size in rhesus macaques.

    J Acoust Soc Am 1997, 102:1213-1221. PubMed Abstract | Publisher Full Text OpenURL

  13. Evans S, Neave N, Wakelin D: Relationships between vocal characteristics and body size and shape in human males: an evolutionary explanation for a deep male voice.

    Biol Psychol 2006, 72:160-163. PubMed Abstract | Publisher Full Text OpenURL

  14. Vannoni E, McElligott AG: Low frequency groans indicate larger and more dominant fallow deer (Dama dama) males.

    PloS One 2008, 3:e3113. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. Vannoni E, McElligott AG: Fallow bucks get hoarse: vocal fatigue as a possible signal to conspecifics.

    Anim Behav 2009, 78:3-10. Publisher Full Text OpenURL

  16. Stearns SC: The Evolution of Life Histories. New York: Oxford University Press; 1992. OpenURL

  17. Marler P: Bird songs and mate selection.

    In Proceedings of the Symposium Animal Sounds and Communication, AIBS meetings: September 1958; Bloomington Edited by Tavolga WN. 1960, 348-367. OpenURL

  18. Fant G: Acoustic theory of speech production. The Hague, Netherlands: Mouton; 1960. OpenURL

  19. Titze IR: Principles of Voice Production. Englewood Cliffs, NJ: Prentice Hall; 1994. OpenURL

  20. Fitch WT, Hauser MD: Vocal production in nonhuman primates: acoustics, physiology, and functional constraints on "honest" advertisement.

    Am J Prim 1995, 37:191-219. Publisher Full Text OpenURL

  21. Vannoni E, McElligott AG: Individual acoustic variation in fallow deer (Dama dama) common and harsh groans: a source-filter theory perspective.

    Ethology 2007, 113:223-234. Publisher Full Text OpenURL

  22. Taylor AM, Reby D: The contribution of source-filter theory to mammal vocal communication research.

    J Zool 2010, 280:221-236. Publisher Full Text OpenURL

  23. McElligott AG, Altwegg R, Hayden TJ: Age-specific survival and reproductive probabilities: evidence for senescence in male fallow deer (Dama dama).

    Proc R Soc B 2002, 269:1129-1137. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. McElligott AG, Hayden TJ: Context-related vocalization rates of fallow bucks, Dama dama.

    Anim Behav 1999, 58:1095-1104. PubMed Abstract | Publisher Full Text OpenURL

  25. McElligott AG, Hayden TJ: Postcopulatory vocalisations of fallow bucks: who is listening?

    Behav Ecol 2001, 12:41-46. OpenURL

  26. McElligott AG, O'Neill KP, Hayden TJ: Cumulative long-term investment in vocalization and mating success of fallow bucks, Dama dama.

    Anim Behav 1999, 57:1159-1116. PubMed Abstract | Publisher Full Text OpenURL

  27. McElligott AG, Naulty F, Clarke W, Hayden TJ: The somatic cost of reproduction: what determines reproductive effort in prime-aged fallow bucks?

    Evol Ecol Res 2003, 5:1239-1250. OpenURL

  28. McElligott AG, Hayden TJ: Lifetime mating success, sexual selection and life history of fallow bucks (Dama dama).

    Behav Ecol Sociobiol 2000, 48:203-210. Publisher Full Text OpenURL

  29. Say L, Naulty F, Hayden TJ: Genetic and behavioural estimates of reproductive skew in male fallow deer.

    Mol Ecol 2003, 12:2793-2800. PubMed Abstract | Publisher Full Text OpenURL

  30. David HA: Ranking from unbalanced paired-comparison data.

    Biometrika 1987, 74:432-436. Publisher Full Text OpenURL

  31. Gammell MP, De Vries H, Jennings DJ, Carlin CM, Hayden TJ: David's score: a more appropriate dominance ranking method than Clutton-Brock et al.'s index.

    Anim Behav 2003, 66:601-605. Publisher Full Text OpenURL

  32. McElligott AG, Gammell MP, Harty HC, Paini DR, Murphy DT, Walsh JT, Hayden TJ: Sexual size dimorphism in fallow deer (Dama dama): do larger, heavier males gain greater mating success?

    Behav Ecol Sociobiol 2001, 49:266-272. Publisher Full Text OpenURL

  33. Specht R: Avisoft-SASLab Pro v.4.38. Berlin: Avisoft; 2006.

  34. Vannoni E, Torriani MVG, McElligott AG: Acoustic signalling in cervids: a methodological approach for measuring vocal communication in fallow deer.

    Cogn Brain Behav 2005, 9:551-566. OpenURL

  35. Boersma P, Weenink D: Praat: doing phonetics by computer. [http://www.praat.org/] webcite

    2009.

  36. McElligott AG, Mattiangeli V, Mattiello S, Verga M, Reynolds CA, Hayden TJ: Fighting tactics of fallow bucks (Dama dama, Cervidae): reducing the risks of serious conflict.

    Ethology 1998, 104:789-803. OpenURL

  37. Reby D, Cargnelutti B, Hewison AJM: Contexts and possible functions of barking in roe deer.

    Anim Behav 1999, 57:1121-1128. PubMed Abstract | Publisher Full Text OpenURL

  38. Rendall D: Acoustic correlates of caller identity and affect intensity in the vowel-like grunt vocalizations of baboons.

    J Acoust Soc Am 2003, 113:3390-3402. PubMed Abstract | Publisher Full Text OpenURL

  39. Owren MJ, Linker CD: Some analysis methods that may be useful to acoustic primatologists. In Current Topics in Primate Vocal Communication. Edited by Zimmermann E, Newman, JD, Jürgens U. New York: Plenum Press; 1995:1-27. OpenURL

  40. Horii Y: Fundamental frequency perturbation observed in sustained phonation.

    J Speech Hear Res 1979, 22:5-19. PubMed Abstract | Publisher Full Text OpenURL

  41. Titze IR, Horii Y, Scherer RC: Some technical consideration in voice perturbation measurements.

    J Speech Hear Res 1987, 30:252-260. PubMed Abstract | Publisher Full Text OpenURL

  42. Sueur J, Aubin T, Simonis C: Seewave: time wave analysis and graphical representation.

    R package version 2.9.0 2009. OpenURL

  43. Charrier I, Harcourt RG: Individual vocal identity in mother and pup Australian sea lions (Neophoca cinerea).

    J Mammal 2006, 87:929-938. Publisher Full Text OpenURL

  44. Gwilliam J, Charrier I, Harcourt RG: Vocal recognition in male Australian sea lions.

    J Exp Biol 2008, 211:2288-2295. PubMed Abstract | Publisher Full Text OpenURL

  45. Lemasson A, Boutin A, Boivin S, Nlois-Heulin C, Hausberger M: Horse (Equus caballus) whinnies: a source of social information.

    Anim Cogn 2009, 12:693-704. PubMed Abstract | Publisher Full Text OpenURL

  46. McElligott AG, Birrer M, Vannoni E: Retraction of the mobile descended larynx during groaning enables fallow bucks (Dama dama) to lower their formant frequencies.

    J Zool 2006, 270:340-345. Publisher Full Text OpenURL

  47. Venables WN, Ripley BD: Modern applied statistics with S. 4th edition. New York: Springer; 2002. OpenURL

  48. Tabachnick BG, Fidell LS: Using Multivariate Statistics. 4th edition. New York: Allyn and Bacon; 2001. OpenURL

  49. Burnham KP, Anderson DR: Model Selection and Multimodel Inference: a Practical Information-Theoretic Approach. 2nd edition. New York: Springer; 2002. OpenURL

  50. Magee L: R2 measures based on Wald and likelihood ratio joint significance tests.

    Amer Stat 1990, 44:250-253. Publisher Full Text OpenURL

  51. Johnson RA, Wichern DW: Applied Multivariate Statistical Analysis. Englewood Cliffs, NJ: Prentice-Hall; 1992. OpenURL

  52. Lachenbruch PA, Mickey MR: Estimation of error rates in discriminant analysis.

    Technometrics 1968, 10:1-11. Publisher Full Text OpenURL

  53. Goutte C: Note on free lunches and cross-validation.

    Neural Comput 1997, 9:1211-1215. Publisher Full Text OpenURL

  54. McGarigal KS, Cushman S, Stafford S: Multivariate Statistics for Wildlife and Ecology Research. New York: Springer-Verlag; 2000. OpenURL

  55. Solow AR: A randomization test for misclassification probability in discriminant analysis.

    Ecology 1990, 71:2379-2382. Publisher Full Text OpenURL

  56. R Development Core Team: Vienna, Austria: R Foundation for Statistical Computing. [http://www.R-project.org] webcite

    2009.

    ISBN 3-900051-07-0

  57. Kokko H: Evolutionarily stable strategies of age-dependant sexual advertisement.

    Behav Ecol Sociobiol 1997, 41:99-107. Publisher Full Text OpenURL

  58. Gil D, Cobb JLS, Slater PJB: Song characteristics are age dependent in the willow warbler, Phylloscopus trochilus.

    Anim Behav 2001, 62:689-694. Publisher Full Text OpenURL

  59. Proulx SR, Day T, Rowe L: Older males signal more reliably.

    Proc R Soc B 2002, 269:2291-2299. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  60. Owren MJ, Seyfarth RM, Cheney DL: The acoustic features of vowel-like grunt calls in chacma baboons (Papio cyncephalus ursinus): implications for production processes and functions.

    J Acoust Soc Am 1997, 101:2951-1963. PubMed Abstract | Publisher Full Text OpenURL

  61. Riede T, Mitchell BR, Tokuda I, Owren MJ: Characterizing noise in nonhuman vocalizations: acoustic analysis and human perception of barks by coyotes and dogs.

    J Acoust Soc Am 2005, 118:514-522. PubMed Abstract | Publisher Full Text OpenURL

  62. Reby D, Joachim J, Lauga J, Lek S, Aulagnier S: Individuality in the groans of fallow deer (Dama dama) bucks.

    J Zool 1998, 245:79-84. Publisher Full Text OpenURL

  63. Komers PE, Birgersson B, Ekvall K: Timing of estrus in fallow deer is adjusted to the age of available mates.

    Am Nat 1999, 153:431-436. Publisher Full Text OpenURL

  64. Farrell ME: Courtship, multiple mating and reproductive synchrony in female fallow deer (Dama dama). In Ph.D. thesis. University College Dublin; 2001. OpenURL

  65. Charrier I, Mathevon N, Jouventin P: Fur seal mothers memorise subsequent versions of developing pups' calls: Adaptation to long-term recognition or evolutionary by-product?

    Biol J Linn Soc 2003, 80:305-312. Publisher Full Text OpenURL

  66. Ey E, Hammerschmidt K, Seyfarth RM, Fisher J: Age- and sex-related variations in clear calls of Papio ursinus.

    Int J Primat 2007, 28:947-960. Publisher Full Text OpenURL

  67. Sanvito S, Galimberti F, Miller EM: Development of aggressive vocalizations in male southern elephant seals (Mirounga leonina): maturation or learning?

    Behaviour 2008, 145:137-170. Publisher Full Text OpenURL

  68. Saez SJ, Martin PM: Androgen receptors in human pharyngo-laryngeal mucosa and pharyngo-laryngeal epithelium.

    J Ster Bioch 1976, 7:919-921. Publisher Full Text OpenURL

  69. Aufdemorte TB, Sheridan PJ, Holt GR: Autoradiographic evidence of sex steroid receptors in laryngeal tissues of the baboon (Papio cynocephalus).

    Laryngoscope 1993, 93:1607-1611. Publisher Full Text OpenURL

  70. Ferguson BJ, Hudson WR, McCarty KS: Sex steroid receptor distribution in the human larynx and laryngeal carcinoma.

    Arch Otolaryngol Head Neck Surg 1987, 113:1311-1315. PubMed Abstract | Publisher Full Text OpenURL

  71. Butler GE, Walker RF, Walker RV, Teague P, Riad-Fahmy D, Ratcliffe SG: Salivary testosterone levels and the progress of puberty in the normal boy.

    Clin Endocrinol 1989, 30:587-596. Publisher Full Text OpenURL

  72. Harries M, Hawkins S, Hacking J, Hughes I: Changes in the male voice at puberty: vocal fold length and its relationship to the fundamental frequency of the voice.

    J Laryngol Otol 1998, 112:451-454. PubMed Abstract | Publisher Full Text OpenURL

  73. Evans S, Neave N, Wakelin D, Hamilton C: The relationship between testosterone and vocal frequencies in human males.

    Physiol Behav 2008, 93:783-788. PubMed Abstract | Publisher Full Text OpenURL

  74. Rolf HJ, Fisher K: Serum testosterone (T) and 5-alpha-dihydrotestosterone (DHT) in male fallow deer (Dama dama L.): seasonality and age dependence.

    Comp Biochem Physiol 1990, 95:445-452. Publisher Full Text OpenURL

  75. Linville SE: The sound of senescence.

    J Voice 1996, 2:190-200. Publisher Full Text OpenURL

  76. Hakkesteegt MM, Brocaar MP, Wieringa MH, Feenstra L: Influence of age and gender on the Dysphonia Severity Index.

    Folia Phoniatr Logop 2006, 58:264-273. PubMed Abstract | Publisher Full Text OpenURL

  77. Fitch WT: The phonetic potential of nonhuman vocal tracts: comparative cineradiographic observations of vocalizing animals.

    Phonetica 2000, 57:205-218. PubMed Abstract | Publisher Full Text OpenURL

  78. Fitch WT, Reby D: The descended larynx is not uniquely human.

    Proc R Soc B 2001, 268:1669-1675. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  79. Reby D, McComb K, Cargnelutti B, Darwin C, Fitch W, Clutton-Brock T: Red deer stags use formants as assessment cues during intrasexual agonistic interactions.

    Proc R Soc B 2005, 272:941-947. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  80. Komers PE, Pélabon C, Stenström D: Age at first reproduction in male fallow deer: age-specific versus dominance-specific behaviors.

    Behav Ecol 1997, 8:456-462. Publisher Full Text OpenURL

  81. Puts DA, Gaulin SJC, Verdolini K: Dominance and the evolution of sexual dimorphism in human voice pitch.

    Evol Hum Behav 2006, 27:283-296. Publisher Full Text OpenURL

  82. Endler JA: Natural Selection in the Wild. Princeton, New Jersey: Princeton University Press; 1986. OpenURL

  83. Charlton BD, Huang Y, Swaisgood RR: Vocal discrimination of potential mates by female giant pandas (Ailuropoda melanoleuca).

    Biol Lett 2009, 5:597-599. PubMed Abstract | Publisher Full Text OpenURL