- Research article
- Open Access
What the hyena's laugh tells: Sex, age, dominance and individual signature in the giggling call of Crocuta crocuta
BMC Ecologyvolume 10, Article number: 9 (2010)
Among mammals living in social groups, individuals form communication networks where they signal their identity and social status, facilitating social interaction. In spite of its importance for understanding of mammalian societies, the coding of individual-related information in the vocal signals of non-primate mammals has been relatively neglected. The present study focuses on the spotted hyena Crocuta crocuta, a social carnivore known for its complex female-dominated society. We investigate if and how the well-known hyena's laugh, also known as the giggle call, encodes information about the emitter.
By analyzing acoustic structure in both temporal and frequency domains, we show that the hyena's laugh can encode information about age, individual identity and dominant/subordinate status, providing cues to receivers that could enable assessment of the social position of an emitting individual.
The range of messages encoded in the hyena's laugh is likely to play a role during social interactions. This call, together with other vocalizations and other sensory channels, should ensure an array of communication signals that support the complex social system of the spotted hyena. Experimental studies are now needed to decipher precisely the communication network of this species.
The origin and maintenance of social group structure is a topic of central concern in vertebrate biology [1–4]. Whereas one approach is to understand the processes that can account for sociality over an evolutionary scale [5–8], a proximal point of view aims to decipher the mechanisms by which the social structure of a group is maintained -or not- over an individual lifetime scale [9–13]. Since Darwin's book on the expression of emotion  and followed by numerous studies in the field, it is well known that information on social status can help individuals adjust their behaviour, for example by avoiding useless fights and polishing social interactions [15–17]. Chemical, visual and acoustic signals have been shown to encode information about sex, kinship, individual identity, morphological cues, as well as motivational and physiological states of the sender [18–28]. As some of these cues can potentially be correlated to fighting ability and dominance rank, the information helps congeners evaluate the emitter's social position within the group [18, 29–38]. We extended these studies by examining the information that is present in one of the acoustic communication signals of a unique social carnivore, the spotted hyena Crocuta crocuta. Spotted hyenas are nocturnal social carnivores, typically living in multi-male, multi-female, "clans" of 10 - 90 individuals. Spotted hyenas are efficient hunters. A lone hyena is capable of capturing prey as large as a wildebeest. Hyenas will also hunt collaboratively, for example to catch zebras . But both sole and collaborative hunting can generate intense competition as clan mates will converge on the carcass. Spotted hyenas have a matrilineal social system similar to that of many old world primates . Within spotted hyena clans, there are separate male and female dominance hierarchies, but all females, and their sub-adult offspring, totally dominate all adult immigrant males. Such female dominance persists in the captive colony at the University of California, Berkeley, where we conducted this study. This colony was established in order to permit the study of the endocrine substrates of dominance and aggression, as well as the basis of "masculinization" of the external genitalia of female spotted hyena, which occurs in utero.
In nature, the dominance hierarchies described above determine priority of access to food. The hierarchical position of an individual is fully acquired around 18 months of age and is not correlated with size or fighting ability: the social status of females is determined by their mothers' social rank [42–44], whereas males' one depends on their sequence of arrival within the social group [45, 46]. During their entire life, hyenas form coalitions against competitors to defend their social ranks . Within their clan, individuals have to cope with a complicated network of social congeners, and this promotes the development of effective abilities to discriminate and rank congeners. Although, it is known that hyenas communicate through visual, chemical and acoustical modalities, there is still much to learn about the nature of exchanged information and the way it is encoded into communication signals .
One of the most striking hyena communication channels is acoustical . These animals are well known for their vocalizations that dominate the nightly soundscape in the African savannah. The vocal repertoire of the spotted hyena is large, with more than ten different vocalizations, many of them being graded into each other which makes them difficult to be classified [39, 49]. As hyenas are primarily nocturnal , vocal signalling is a privileged channel, used for both long- and short-range communications. For instance, "whoops," with long inter-whoop intervals, are primarily used to signal separated individuals, supporting within- and between-clans acoustic interaction [51, 52]; conversely, the "grunt"  ("soft growl" in ) is uttered during close meeting of clan mates and remains barely audible after a few meters of propagation. Previous studies have focused on the long-distance "whoop" call, showing that it supports information related to sex and individual identity and thus may allow discrimination between clan members and alien individuals [48, 51, 53]. Although they might play an important role in mediating the relationships within the clan, calls, other than the whoop, of adult hyenas have been neglected. Although inferences have been drawn from the contexts in which the various vocalizations are emitted, the information that such vocalizations contain has not yet been investigated .
Among the vocalizations used during interactions between adult clan-mates, we chose to focus on the "giggle call", often referred as the hyena's laugh (Figures 1 (see additional files 1, 2, 3 &4 for sound) 2 (see additional files 5, 6, 7 &8 for sound); ). Giggles are high pitched sounds emitted in bouts, mainly when hyenas are feeding together on a prey [39, 49]. This call has been described by field observers as a submissive vocalization uttered by an individual in front of a dominant [39, 49]. Although giggles are emitted during close-range interactions between two or more individuals, they are loud and can be easily eavesdropped by other clan-mates. In nature, giggles are commonly emitted during competitions between dominant and subordinate animals such as those that would occur for food at a carcass [39, 49], but identification of individual giggles by human observers is difficult because of the simultaneous emission of such giggles by multiple individuals. In the Berkeley colony, it was possible to record giggles emitted by individual hyenas when presented with food. In the present report, we describe the use of acoustic structure analysis, in order to determine whether the giggle of spotted hyenas encodes information about sex, age, dominance and individual identity, and if it could allow congeners to assess an emitter's social status based on its laughing cues.
The animals that participated in the study are members of a captive breeding colony of spotted hyenas, maintained at the Field Station for the Study of Behavior, Ecology, and Reproduction (FSSBER) at the University of California, Berkeley. Currently, this colony houses 26 hyenas: 14 adult females, 10 adult males, and 2 sub-adult males (less than two years of age at the time of testing). Animals are sexed and their life history is known. The reproductive founders of the Berkeley colony were collected in 1985 in the Narok district of Kenya. The present study was conducted in September-October 2008, and all the animals currently at the field station were born in captivity. The FSSBER is fully accredited by the American Association for the Accreditation of Laboratory Animal Care and the animal care also met all institutional guidelines. The behavioural procedures were approved by the Institutional Animal Care and Use Committee (IACUC) of the University of California, Berkeley.
In the field station, individuals are housed in groups of two (n = 7 dyads) or three (n = 1 one triad) in outdoor or semi-outdoor enclosures (enclosure mean area = 330 m2, min-max = 147-929 m2). Table 1 provides detailed for information about individuals whose giggles were analyzed in the present report. The dominance status within each dyad was assessed on a daily basis by caretakers familiar with the animals. There was perfect agreement between two caretakers who independently assigned dominance ranks, based upon observation of submissive (e.g., retreat, tail between legs) and aggressive acts (e.g., push, stand over, chase). For each dyad, one animal was qualified as 'dominant', the other as 'subordinate' (Table 1).
The colony has been supported by the National Institute of Mental Health and, more recently, the National Science Foundation, for studies of the endocrine substrates of genital masculinization and dominance in the female spotted hyena [41, 54]. Accordingly, the majority of subjects observed in the present study were treated with compounds that blocked the actions of androgens or estrogens in utero, or were gonadectomized at varying stages of development (see Table 1). Thus although the hormonal and developmental state of the animals were not a designed factor in our study, the effect of these treatments had to be examined so that they could be discarded (or alternatively further examined if they had showed interesting results). The animals in the study were initially assigned to the following groups: (A) "control" (no treatment or simple vasectomy; 2 males, 3 females), (B) treatment with "anti-androgens" (flutamide and finasteride, or Casodex and finasteride) in utero (3 males, 7 females; for details of drug treatments, see ); or (C) treatment with an aromatase inhibitor (Letrozole, preventing the synthesis of estrogen) in utero (1 male, 1 female; for details of Letrozole treatment see ). In addition, (D)1 male and 1 female were gonadectomized during the first months of life. When compared with control hyenas, gonadectomized hyenas would have had low concentrations of plasma testosterone, or estradiol, at the time of testing . Male hyenas treated with anti-androgens in utero have reduced concentrations of testosterone as adults ; while female spotted hyenas treated with anti-androgens in utero displayed elevated concentrations of estradiol . In Hyenas, males are mature at 2 years old, and females at 3 (Glickman, pers. obs.).
To test the effect of hormonal treatments, we performed statistical analyses by dividing the animals into the four groups mentioned above. Alternatively, we analyzed the data by dividing the experimental animals into two groups (anti-estrogen and anti-androgens) by including the two gonadectomized animals into these groups (the female into the anti-estrogen and the male in the anti-androgen group). Finally, we also tested the control animals against all experimental animals. These various groupings were done to potentially increase the statistical power of the tests given the small sample size in each group. As described below, none of the analyses revealed an effect of hormonal treatment.
We recorded hyena's vocalizations in the morning, prior to feeding (recording chain: Shure Model 16A condenser microphone placed at 3-5 meters from the vocalizing animal, connected to a MARANTZ PMD670 solid state recorder; sampling frequency = 44100 Hz; frequency response of the recording chain: ± 2.5 dB over the 50 to 15,000 Hz frequency range). Giggling was elicited by teasing the focused animal with a bone or a piece of meat presented through the fence of the enclosure. In most cases (and in particular for subordinate animals), the recorded animal was isolated from its cage-mate to avoid fights and potential injuries.
To limit the potential impact of pseudo-replication, each individual was recorded over a minimum of 4 different days (10 females, 7 males; see Table 1 for details). Giggles are emitted in bouts which consist of a rapid succession of very short calls or notes (Figure 1A & 2). In this study, we focused our analysis on the information that could be extracted from both the average acoustical features of single notes and from features describing the range of giggle notes produced by single animals. Although, there is certainly additional information in the structure of a particular sequence of notes within bouts these were not examined in this analysis; analyzing such sequence effects would require the estimation of transition probabilities for particular note types and this initial data set is too small for such characterization. The notes used for the analysis were selected pseudo-randomly from the recordings to avoid noisy or transitional sounds such as whines within the giggles. The raw data consisted of 254 giggle bouts and a total of 1807 notes (Average notes/bout = 7.1, Standard deviation = 4.1). The final data set consisted of 41 ± 19 (std) notes per individual (minimum = 20; maximum = 98), or a total of 695 notes. Balanced statistical analyses were performed by repeatedly choosing 20 random notes for each animal.
To characterize the giggle notes' complex acoustic structure, we performed the measurement of 13 variables, in both temporal and frequency domains (Figure 1B, C). The acoustical analysis was done using custom software written in MATLAB (Ver 7.6, The Mathworks, Cambridge MA). The time-varying fundamental frequency of notes was extracted using a custom algorithm that combined two separate measurements of fundamental frequency to obtain a best guess using a Bayesian approach. The two methods were: 1) a cepstral analysis where the inverse time corresponding to the highest peak in the Fourier transform of the log spectrogram was taken as the fundamental, and 2) a direct analysis of the short-time Fourier Transform (STFT) of the sound waveform where the smallest distance between all major successive peaks including a peak at zero was taken as the fundamental. The spectrogram and the STFT where obtained with Gaussian shaped windows with a time-frequency scale of 3.2 ms/50 Hz measured by the standard deviation parameter of the Gaussian window in time or frequency. The Bayesian approach consisted of assigning a likelihood to each of these two measurements and a prior. The likelihood was obtained by estimating the probability of the possible fundamental frequencies given the relative size of the peaks in the cepstrum or STFT. The prior depended on previous history, more precisely on the local derivative of the fundamental obtained from the prior three measurements (the prior was uniform for fundamental below 1 kHz for the initial measurements). The algorithm returned the best guess of the fundamental and its probability. Guesses with low probabilities were dismissed as sounds lacking clear periodicity. The performance of the algorithm was verified visually for all calls in the database and was deemed to perform very similarly to what an expert acoustician might extract from analyzing the spectrogram. Using this methodology, we extracted the time-varying fundamental frequency of each note in our data set with a (over-sampled) resolution of 1 ms. The following parameters were then calculated from these data: mean of the fundamental frequencies ('Mean F'), maximum of the fundamental frequencies ('Max F'), minimum of the fundamental frequencies ('Min F') and the coefficient of variation of the fundamental frequencies ('CV F'). We also extracted acoustical parameters from the frequency spectrum obtained for entire notes. The frequency spectrum was calculated with the Welch average periodogram method using a 23.2 ms Hanning window . From the frequency spectrum, we measured the frequencies corresponding to the first, second and third quartiles of energy ('Q1', 'Q2', 'Q3', respectively), the mean frequency ('Mean S'), and the standard deviation of the spectrum ('SD S'). We also extruded three additional measures of spectral shape: skewness, kurtosis and entropy. The spectrum skewness ('Skew') was calculated as ; the spectrum kurtosis ('Kurt') was calculated as , where f i = frequencies, p(fi) = the normalized power spectrum () and σ2 is the variance of the spectrum given by . The spectral entropy ('Ent') was calculated as . These values of spectral shape have been recently proposed for acoustical analysis (Seewave R software http://rug.mnhn.fr/seewave). Sound duration ('Dur') of each note was measured visually from the oscillogram and its corresponding spectrogram.
All measures described above are used to describe the acoustical structure of single notes, or when averaged to describe giggles notes produced by an individual, as a prototypical single giggle note. However, it is clear from experiencing giggle sounds (personal experience), or, for the reader of this article by visual inspection of spectrograms (Figure 1 &2) as well as listening to the example sound clips (see Additional Files) that individual hyenas modulate both the fundamental frequency (pitch) and the type of giggle note (timbre) that they produce. The structure present in the succession of giggles notes in a bout is rich, almost musical, and potentially very informative. In this study, we just began the analysis of this structure by calculating the coefficient of variation (CV) of the mean fundamental frequency ('CVMeanF') and of the spectral mean ('CVMeanS') from all the giggle notes obtained from each individual. In other words, for each hyena, we calculated from all its giggle notes, one CVMeanF and one CVMeanS. These measures provide a simple quantitative measurement of how much each individual varies the pitch and the timbre of the giggle note that it produces. Since the CV is obtained by dividing the standard deviation by the mean, we also obtained the mean of 'Mean F' ('GrandMeanF'), the standard deviation of 'Mean F' ('SDMeanF'), the mean of 'Mean S' ('GrandMeanS') and standard deviation of 'Mean S' ('SDMeanS'). We examined these means and standard deviations to assess their respective contribution to the coefficient of variance.
For all multivariate statistical analyses, the raw values of the 13 acoustical parameters of interest were centred and normalized (i.e. transformed into z-scores) to insure correct weighting since our acoustical parameters had different units (e.g., 'Mean F' and 'CV F'). Normality was visually inspected using the "normplot" command in Matlab.
We analyzed the differences between giggles notes from different individuals or groups of individuals using Discriminant Function Analysis (DFA) . The DFA was performed in Matlab (using the functions 'manova', and 'classify') and repeated in R for further validation as explained below. A DFA is composed of two steps: in the first step a set of discriminant functions is obtained from a training data set; in the second step these functions are used to test classification on a validation set. We chose the linear discriminant functions given by the eigenvectors of the ratio of the between and within covariance matrices; this approach is equivalent to a multivariate analysis of variance (MANOVA). The MANOVA assumes a jointly normal distribution for the parameters that describe the variability between and within individuals and returns the number of significant discriminant functions. The statistical significance of the functions is assessed by assuming normal distributions and calculating Wilk's lambda test statistic. In the second step, these discriminant functions are then used to classify giggle notes chosen from a validation data set. This cross-validation step gives a measure of the effect size (the percent correct) and of the statistical significance by comparing the percent correct and its standard error to chance. The measure of standard error is obtained by analyzing the percent correct assignment of 100 random selections of the original data set divided into a fitting and testing set. In all cases the training set consisted of 19 randomly selected giggle notes per animal and the discrimination was tested on a different randomly selected note (Note 19+1 = 20, the smallest number of notes that we obtained for an animal in our data set). As shown in Figure 3 and Results, 100 random selections were more than sufficient to obtain good estimates. By performing this cross-validation step, not only does one obtain a desirable measure of effect size (the percent correct) but also the assumption of normality is relaxed.
From the cross-validation results, we could also extract a complete confusion matrix: the conditional probability of guessing that the test giggle note came from individual i when in fact it was emitted by j:
From the confusion matrix, one can also obtain a measure of the goodness of the classification (effect size) by estimating the mutual information between guesses and actual values :
Here the probability of the actual values, p(j), is uniform because the same number of testing giggles is used for each individual or group. When the DFA is used to classify individuals (n = 17), p(j) = 1/17 = 0.058. p(i, j) and p(i) are obtained from the confusion matrix and Bayes' theorem. It should be noted that the MI will have a positive bias for small number of testing notes and will have a negative bias (or be bounded) when the number of individuals is too small. For example, for 17 individuals, the upper bound for estimate of MI from the confusion matrix is log2(17) = 4.75 bits. As shown in the results, we were far from the regime were the negative bias becomes important and we corrected for positive bias effects. The measure of effect size provided by the MI is then independent of the number of animals or calls and, moreover, can be normalized by the length of the sound (by division) to get an effect size in bits/s. It can therefore be used to compare results on the information bearing content of communication sounds across experiments and species.
The straightforward DFA described above was used to evaluate the information present in giggle notes about individual identity. However, this procedure needs to be further modified when DFA is used to assess differences among groups of individuals classified by their age, social status or treatment effect. For those analyses, there are two additional statistical considerations: pseudo-replication and nested and/or interacting effects. The pseudo-replication comes from the fact that, when analyzing effects other than individual identity, giggle notes from the same animal are not statistically independent measures. In fact, the individual identity is a nested effect: one needs to consider differences in notes from different individuals as a separate effect that might explain some or all of the differences observed between groups. This point was clearly made by Mundry and Sommer  who suggested addressing the nested individual factor in a DFA by comparing the percent correct obtained in the analysis to the distribution of percent correct values obtained by randomly assigning the group identity to each individual. In our analysis, this distribution was obtained from 1000 randomly created data sets where the group identity of each individual is permuted in each set while preserving the number in each group. Note that in the random sets any information that is dependent on identity (such as belonging to a particular dyad) is also lost. This procedure is called permuted DFA (pDFA). The pDFA was performed in Matlab and also calculated with a R routine provided by Roger Mundry . The statistical conclusions from both the Matlab and R routines were identical although we did observe small differences in the percent correct measures. In this study, we report the results obtained from the Matlab routines.
To further examine the potential interaction between the individual differences and the different conditions (age, sex, dominance, treatment), we performed a multiple linear regression (least-square and robust) with age, sex, dominance and treatment as predictors for the average value of the discriminant functions obtained in the DFA for individual differences. By using the average value for each individual in the regression we eliminated the pseudo replication problem that could lead to inflated significance. To minimize the number of statistical tests, this regression was performed for the first discriminant function, then the second, and so forth until we failed to find an effect. When age was a significant factor, we also used an analysis of covariance (ANCOVA) as an exploratory tool to investigate potential interactions between age and the other independent variables; in the ANCOVA age is treated as a regular scale variable in linear regression and the other conditions are examined to determine whether separate lines for each condition are warranted (the interaction effect). The p-values were then adjusted for multiple tests using the Bonferronni correction.
We also performed six linear regression analyses (least-square and robust) with 'GrandMeanF', 'SDMeanF', 'GrandMeanS', 'SDMeanS' "CVMeanF" and "CVMeanS" as dependent variables and age, sex, dominance and treatment as predictor variables. The p-values were corrected for multiple tests using the Benjamini and Hochberg False Discovery Rate . ANCOVA was then used to examine the presence of interaction effects between age and other predictor variables when age was a significant factor in the linear stepwise regression analysis. Finally, considering that the dominance/subordinate status was assessed for each individual in respect to its cage-mate, we performed paired t-tests considering hyena dyads. For the single triad in our data (see table 1), Jambo was assessed as being dominant over his two cage mates, Tembo and Haji; thus, we took the average acoustical parameters of Tembo and Haji to characterize the subordinate call.
Giggle note acoustic structure
All analyzed giggle notes (n = 695) showed the same general structure, i.e., short, broad-band signals, with some modulation in frequency (Figure 1B). Average note duration was 69 ± 18 ms (mean ± standard deviation). The mean fundamental frequency was 547 ± 146 Hz, with a maximum of 741 ± 180 Hz and a minimum of 399 ± 143 Hz. The fundamental frequency was modulated in time with a coefficient of variation of 0.155 ± 0.101. The energy was concentrated on the lower- harmonics as the frequency values at the first (25%), second (50%) and third (75%) quartiles of energy are 635 ± 398 Hz, 848 ± 556 Hz, and 1196 ± 717 Hz respectively. However, this distribution showed great variation among individuals. For example, the frequency at the third quartile ranged from 754 ± 315 Hz to 2279 ± 996 Hz, depending on individuals. The spectrum mean frequency was 979 ± 495 Hz, with a mean standard deviation of 544 ± 229 Hz. The frequency spectrum showed a positive skewness (skew equals 2.56 ± 1.93), underlying that most of the energy was concentrated over the low part of the signal frequency bandwidth. The leptokurticity of the spectrum (kurtosis of 18.5 ± 24.5) and the low values of spectral entropy (0.52 ± 0.12) quantify the fact that the energy was concentrated at a few frequency values (peaks). Differences in giggle notes from different individuals are illustrated in Figure 1A and 2 where giggle bouts from four individuals are shown: as can be seen the two bouts coming from the same individual bared similarities while two bouts coming from different individuals were clearly more distinct.
Acoustic Structure in Succession of Giggle Notes
Giggle notes are produced in bouts as shown in figures 1 and 2. In our experiment, a typical giggle bout lasted less than 1 second and had approximately 7 notes (Avg = 7.1 notes/bout). The longest bout we recorded lasted 4s and had 28 notes. The shortest had two notes and lasted 200 ms.
The sequence of notes in a giggle bout has also a rich acoustic structure, in the sense that both the pitch and the timbre of giggle notes appear to change in an orderly fashion from the first to the last note in a bout. In this study, we coarsely began to quantify this structure by examining the variability in the mean fundamental and the spectrum mean frequency for each individual. The values for those two measurements are shown in the top panels in Figure 1A for the succession of giggle notes for the first example bout for each animal. In these examples, Ursa varied the mean fundamental (F) and the spectrum mean frequency (S) across its giggle notes more than did Kombo. It can also be clearly observed in the spectrogram that within a giggle bout, Ursa used more different "types" of notes than Kombo (see also the corresponding example sound clips). The same difference can be observed between the giggle notes produce by Winnie (more variable) and Kadogo shown in figure 2. Ursa and Winnie were subordinate animals while Kombo and Kadogo were dominant animals. Ursa and Kombo (shown in figure 1) were control animals while Winnie and Kadogo (shown in figure 2) were treated with anti-androgen and anti-estrogens respectively. We will show below that dominance status and age both affect the variability of giggle notes within animals but that hormonal treatment did not appear to affect the acoustic structure of giggle vocalizations as described by the measures we used. Across all animals, the mean CV for the mean fundamental was 0.19+- 0.05 (Range: 0.12-0.34). The mean CV for the spectrum mean frequency was 0.35+-0.1 (Range: 0.22-0.57).
Individual identity and age of the sender
The DFA identified significant acoustic differences between individuals, calculating eight significant linear discriminant functions that allowed maximizing individual separation (Figure 3A). Table 2 shows the variance explained by each of the eight functions and the variable loadings on each function for the three most important factors. The spectrum mean frequency and the mean fundamental were the two main factors that separated individuals on the first function. The second function relied on the spectrum standard deviation, the note duration and the third quartile of energy, while the third function depended on the spectrum standard deviation, the spectral entropy, and the mean fundamental. Thus, while the first function mostly reflected the fundamental frequency (pitch), the second and third function quantified the spectral envelope (one of the factors that affect timbre).
The results of the cross-validation step showed that, although perfect individual identification was not achieved, the rate of success was greatly above chance (32% versus 6%) (Figure 3C). This is not surprising given the high number of significant acoustical parameters. However, there was also extensive variation in the classification success rate across individuals (average = 32%, individual range from 6% to 62%). In other words, in our sample, some animals appeared to produce similar giggle notes (e.g see Kadogo, Rocco and Nakuru in Figure 3A and 3B) and were thus not easily classified while other animals produced relatively unique vocalizations (e.g. Jambo, Kombo and BJ in Figure 3A and 3B).
To obtain another measure of effect size that could be used to compare with other studies, we also calculated the mutual information from the confusion matrix that shows the joint probability of the actual and predicted individuals in cross-validated data (Figure 3B). The mutual information was 1.09 bits (well below the ceiling value or log2(17) = 4.08 bits). Since this information was obtained for single giggle notes of average duration of 69 ms, the resulting baud rate is 15.8 bits/s. A mutual information of ~1 bit/giggle means that a giggle note carries the same amount of information as a variable that would be able to perfectly divide any hyena as belonging into one of two groups. For giggles, this same amount of information actually allows one to associate the hyena as belonging to one of many groups but not with 100% accuracy (for example in our study as 1 of 17 individuals but with 37% accuracy) .
The post-hoc linear regression with age, gender, dominance and treatment as predictor variables and the first discriminant function as a dependent variable was not statistically significant (F(4,12) = 2.4, p = 0.1) but we noted that the coefficient for age was significant both for the least-square and robust regression methods (p = 0.01 and p = 0.03 respectively). Older animals had lower values for the first function than younger animals (Figure 3D, left panel), which, given the loadings of this function, appears to reflect a lowering of the fundamental with age (see Figure 3D right panel and also below). The ANCOVA confirm the statistically significant effect of age (F(1,15) = 10.4, p = 0.005) but showed no statistically significant interaction between age and the other conditions (sex, dominance, treatment). The linear regression did not reveal any statistically significant factors for the second discriminant function and the other discriminant functions were therefore not tested to minimize the number of tests and prevent the risk of false positives. Thus, while the first function might be used to classify individuals based on their age, the other discriminant functions capture idiosyncratic features of each individual's giggles.
In accordance with these previous analyses, the multiple linear regression performed on 'GrandMeanF', 'SDMeanF', and 'GrandMeanS' showed an influence of age on these parameters (see Table 3), and thus confirmed the age-linked lowering of the giggle note mean fundamental frequency. From the age of 2 to 20 years-old, the mean fundamental frequency dropped from 650 Hz to less than 450 Hz (i.e. -30%). The linear regression line for that relationship is shown on the right panel of figure 3D (t(15) = -2.93, p = 0.01).
The three hormonal treatments were found to result in distinctive giggle notes when the data were analyzed with the MANOVA (P < 10-3; DF(1) = 1.135 CV F + 0.688 Ent - 0.527 Max F F; DF(2) = 0.925 Max F + 0.823 SD S - 0.793 CV F) but the pDFA showed that these difference were primarily the result of individual differences. The pDFA showed that 48% of the notes in the cross-validation set could be classified correctly (versus chance at 25%) but similar performance or better could be achieved with four random groups 12% of the time. Similar non-conclusive results were found when experimental animals are classified along two treatment groups (anti-androgen vs anti-estrogen) or lumped into one group (see methods). With two treatment groups the pDFA yielded a 68% of mean correct assignment versus 55% chance but p = 0.59. The step-wise linear regression analysis for the CVMeanF, CVMeanS, GrandMeanF, SDMeanF, GrandMeanS and SDMeanS showed no effect of treatment (Table 3). Therefore, although we cannot exclude the possibility that treatment had an effect on the acoustical quality of the giggles, if this effect exists it is small and we do not have enough power in our experiment to measure it or it is present in acoustical features not characterized here.
Dominance and Sex
The DFA found one significant linear discrimination function that could be used to classify giggle calls into dominant/subordinate groups (P < 10-3; DF(1) = 2.259 Mean S - 1.392 Q2 - 0.956 Q3). All the variables loading on this function dealt with the distribution with energy among the spectrum or, in other words, the spectral envelope. However, the pDFA showed that the successful classification was due to the nested effect of individual differences. Although, the cross-validation data showed that notes could be classified correctly 58% of the time (versus 50% for chance), the random permutations show that the same classification rates would be obtained 82% of the time (i.e. P = 0.82) in random assignments of individuals to two arbitrary groups (see figure 4A). Thus the results from pDFA did not support the idea that dominance status can be extracted from the acoustical characteristics of single giggles but instead showed that there is enough individual information in giggle notes that particular acoustical dimensions can be found to coarsely separate our 17 hyenas into two arbitrary groups.
On the other hand, the analysis of acoustical structure across giggle notes of single individuals showed a robust effect for dominance: the linear regression analysis showed that the coefficient of variation of the mean spectral frequency ('CVMeanS') calculated from all giggle notes for each individual was higher (ie. giggle notes are more variable) in subordinate animals than in dominant animals and that this effect decreased with age (see Figure 4B and Table 3). The ANCOVA post-hoc analysis confirmed that this interaction effect was significant (see Table 3). When the data was analysed as 8 dyads (and without taking age as a factor), a paired t-test also showed a significant effect for 'CVMeanS' (mean CVMeanS difference = 0.082, t(7) = 2.387, p = 0.048) but not for any other acoustical measure.
Sex was found to be distinctive by MANOVA (P < 10-3; DF(1) = 1.262 Ent -1.228 Q2 + 1.119 Mean S) but the pDFA shows that this separation could also be due to individual differences. In the cross-validated data sets, we found that 57% of notes could be classified according to gender but permutations show that this level of performance would be found 87% of the time. The stepwise regression analysis for the CVMeanF, "CVMeanS" 'GrandMeanF', 'SDMeanF', 'GrandMeanS' and 'SDMeanS' showed no effect of sex (Table 3).
Giggles bear information about the sender
According to the analysis presented above, the hyena's laugh potentially encodes information about individual identity, dominant/subordinate status and age, giving receivers cues to assess the social position of an emitting individual. The discriminant functions used to separate individuals support the idea that information about individual identity is primarily encoded by pitch (1st discriminant function) and by the energy distribution among the frequency spectrum of giggle notes (2nd and 3rd discriminant functions); note duration was also involved (2nd discriminant function). As the first discriminant function correlated with age, pitch was age-dependent. Following our results, the pitch aspect of an individual signature is likely to change when the animal gets older. Besides this, age independent identity and other individual characteristics are mainly encoded by the energy distribution among the frequency spectrum. The dominant/subordinate status appears to be encoded by variations of the spectrum mean frequency within giggles of an individual. The fact that giggle notes are emitted in bouts thus represents a crucial aspect for the potential assessment of this social cue. That way, the multiple giggle notes could not only provide redundant information but also additional (synergistic) information from the structure created by particular note sequence.
We did not find reliable effect of sex even after correcting for the age effect. This result is not universal in hyena vocalizations: an analysis of groans by our group (in prep.) showed a very clear difference between males and females with males showing higher pitch than females. Giggles might therefore be less sexually dimorphic than groans and potentially more influenced by other factors. In addition, it is possible that pre-natal hormonal treatments also affected sexually dependent characteristics.
From a methodological point of view, the results from the classical and pDFAs show that while, information coding of age and individual identity within the giggle note are quite robust, any information that we uncovered about gender, dominant/subordinate status or hormonal treatment could be due solely to individual effects. Our data set and analyses therefore illustrate the importance of a careful use of DFA statistics of acoustic features when searching for vocal signatures, a point that was first stressed by .
Individual identity coding
Previous studies, investigating gregarious mammals' and birds' acoustic signals in the context of mate or parent-offspring recognition, have shown that individual identity can be coded by one or several of the following sound features: the fundamental frequency ("pitch"), the energy distribution among the spectrum ("vocal timbre"), the characteristics of the frequency modulation of the fundamental (e.g., penguins: ; gulls: ; shearwaters: ; marmots: ; fur seals: [69, 70]). All the acoustic signals supporting individual identification are usually highly redundant with regards to the coding process and thus resistant to masking effects. For example, non-nested birds use a secured code for long-range propagation of sender's identity that relies on lower harmonics band and frequency modulation of the fundamental. Both cues resist well to propagation-induced modifications such as filtering or scattering by obstacles. In nest-building species where acoustic exchanges occur mostly at short range, frequency modulation is almost absent and individual identity is encoded by pitch and energy distribution among the spectrum, the latter being highly sensitive to frequency-filtering during long-range propagation . The present analysis showed that giggles encode individuality mostly by these two last cues, the frequency modulation being not primarily involved. This contrasts with another hyena vocalization, the whoop, where individual signature should be encoded by the pronounced frequency modulation, as visual inspections of spectrograms showed . Thus and in spite of instabilities in their frequency spectrum , whoops might be more reliable labels of individual identity than giggles over long distances.
The reasons for giggling: a puzzling situation
Up to now, no systematic study has tried to decipher the behavioural context in which giggles are emitted, nor their role in the society network. Field researchers have observed that these vocalizations are often produced during food contests, by animals that are prevented from securing access to the kill by the intervention of higher ranking individuals. More generally, giggling can occur in non-feeding situations, as the result of simple threat from another individual. Giggles have thus been considered as submissive signals [39, 49]. The present study was done in captivity and our behavioural observations were certainly biased. Nevertheless, we observed that hyenas were giggling rather as a result of frustration (we kept the bone or the piece of meat out of their reach) than of harassment or chase.
Although they are emitted during close-range interactions between two or three individuals, giggles are loud and propagate over great distances (the authors pers. obs.). They are thus extremely susceptible to eavesdropping by remote receivers. Giggles may attract other spotted hyenas as well as lions Panthera leo, and even vultures . During his study in the Ngorongoro (Tanzania), Kruuk  noticed that lions often eat prey previously killed by hyenas. Giggles could benefit an individual hyena is different situations. It is known, that a solitary hyena has no chance when confronted to a lion, whereas a hyena group often can "mob" one or two lions, and maintain or gain access to a carcass . Thus a lone hyena encountering a kill dominated by lions could use its giggle call to rally its clan. The attraction of a neighbouring hyena clan or a lion group can also be an issue, as intense competition and giggling may occur over a kill currently mob by the giggling hyena's clan. In this situation, one can hypothesize that a dominant hyena might allow a subordinate hyena access to food to prevent it from giggling and attracting further attention from unwanted competitors. Or potentially giggling by multiple subordinate animals (giggle chorus) could serve as a distraction for the more dominant animals. However, these hypothetical situations could also potentially incur cost for the emitters who risk being completely deprived of food if the competing clan or lion group takes over. Cooperation and competition are everyday components of a spotted hyena's life. When hearing a giggling individual, clan-mates hyenas could get some information about "who" is currently in a competitive situation (in terms of individual identity, age class, dominant/subordinate status) and decide to join the giggler, or conversely to ignore it or move away. On a larger scale, giggles from a hyena group could attract conspecifics, allowing more successful "mobbing" of lions. It is also interesting to note that the loud giggle call is absent from the vocal repertoire of the sympatric but less social Brown hyena, Hyaena brunnea. Field observations and playback experiments are needed to determine the attraction potential of giggles towards allies versus intra- and inter-specific competitors, and to assess the cost-to-benefit balance of this vocalization. The possibility of monitoring the size and the composition of a feeding group from a distance should also be investigated to assess if a clan could get information about its neighbours by simply eavesdropping on their giggles. This information could also be of interest to researchers with conservation purposes.
Giggles and the spotted hyena's repertoire: acoustic tools for a social network
Primarily emitted during common confrontations, such as those occurring over food, giggles may play an important role in spreading of individual-related information among clan-mates. The social organization of a hyena clan involves a matrilineal hierarchy. The majority of adult males are immigrants and they are "queuing" for social status, i.e., the most recently arrived male is at the bottom of the male hierarchy . Since males depart from their natal clan and attempt to join a new clan when mature, the age information embedded in the giggle could be highly informative.
In our study, hyenas were housed in dyads (and one triad) where a dominant/subordinate ranking was quickly established. The giggle was also produced in a conflict with a human which we attribute to have a neutral social rank. Although the bi-modal grouping and our experimental protocol allowed us to clearly divide animals into subordinate and dominant, it is a situation which is far from the hierarchical rank that could be found in wild clans of 10 to 90 individuals. It is possible that animals in the "middle" of the rank produce different giggle bouts when these are directed to higher ranking versus lower ranking animals. This hypothesis would also imply acoustical structure that is dependent not on morphological characteristics but on the recognition of social context and the possibility of vocal plasticity. This hypothesis could be tested both in the field and in our captive colony by housing animals in larger groups.
The vocal repertoire of spotted hyenas is very large and most of the calls should play a role in the regulation of the social network. Status signalling is almost certainly not restricted to giggles. The whoop is used to transfer information to remote congeners and may allow individual identification of the sender as well as its current emotional state over distances up to 5 km [51, 53, 73, 74]. Not only are whoop calls highly idiosyncratic making but it is also known that mothers will respond to the whoop from their own cub more frequently than to the whoop of a non-kin cub, demonstrating that the individual characteristics in the whoop sounds are recognized and utilized . It is likely that close-range vocalizations like the grunt, groan and growl, are also multi-informative [, Page et al. in prep.]. The use of both 'public' and 'private' signalling [75, 76], would enable hyenas to manage their social interactions with great precision. In this context, giggles might play a dual role, addressed to both nearby clan-mates and remote potential allies.
Spotted hyenas demonstrate high cognitive skills, like their ability to recognize third-party relationships , a characteristic that has been found in a restricted number of animals (e.g., primates ; birds ). Their substantial vocal repertoire should play a very important role, by providing eavesdroppers with a number of important cues about the emitters' identity and characteristics. As the present study showed, the giggle is likely to carry a broad range of messages, certainly not all perfectly reliable, but sufficiently informative to play a role during social interactions. Information carried by vocalizations, together with chemical, tactile and visual channels [48, 80], ensure to the spotted hyena an array of communication signals which underlie its complex social system. More research, with a particular emphasis on experimental studies with playback - or involving manipulations -, is needed to specifically determine whether the information bearing structure in giggle sounds described here is actually used by the spotted hyena and, more generally, to decipher the complex communication network of this species. Comparing the spotted hyena's communication system with the one of other, less social, Hyaenids  would also be of great interest and facilitate understanding of the relationships between sociality and animal signals.
Maynard Smith J, Szathmary E: The major transitions in Evolution. 1995, San Francisco, Freeman
Clutton-Brock T: Breeding together: kin selection and mutualism in cooperative vertebrates. Science. 2002, 296: 69-72. 10.1126/science.296.5565.69.
Krause J, Ruxton GD: Living in groups. 2002, New York, Oxford University Press
Frank SA: All of life is social. Cur Biol. 2007, 17: R648-R650. 10.1016/j.cub.2007.06.005.
Hamilton WD: The genetical evolution of social behaviour. I J Theor Biol. 1964, 7: 1-16. 10.1016/0022-5193(64)90038-4.
Dugatkin IA: Cooperation among animals: an evolutionary perspective. 1997, Oxford, Oxford University Press
West SA, Griffin AS, Gardner A: Evolutionary explanations for cooperation. Cur Biol. 2007, 17: R661-R672. 10.1016/j.cub.2007.06.004.
Hochberg ME, Rankin DJ, Taborsky M: The coevolution of cooperation and dispersal in social groups and its implications for the emergence of multicellularity. BMC Evol Biol. 2008, 8: 238-10.1186/1471-2148-8-238.
Waser PM: Individual recognition, intragroup cohesion and intergroup spacing: evidence from sound playback to forest monkeys. Behaviour. 1977, 60: 28-74. 10.1163/156853977X00270.
Jorgenson JW, Novotny M, Carmack M, Copland GB, Wilson SR, Katona S, Whitten WK: Chemical scent constituents in the urine of the red fox (Vulpes vulpes) during the winter season. Science. 1978, 199: 796-798. 10.1126/science.199.4330.796.
Gosling LM, McKay HV: Competitor assessment by scent-matching: an experimental test. Behav Ecol Sociobiol. 1990, 26: 415-420. 10.1007/BF00170899.
Heg D, Bachar Z, Brouwer L, Taborsky M: Predation risk is an ecological constraint for helper dispersal in a cooperatively breeding cichlid. Proc R Soc Lond B. 2004, 271: 2367-2374. 10.1098/rspb.2004.2855.
Pollard KA, Blumstein DT: Time allocation and the evolution of group size. Anim Behav. 2008, 76: 1683-1699. 10.1016/j.anbehav.2008.08.006.
Darwin C: The expression of the emotions in Man and animals. 1998, Oxford, Oxford University Press first published 1872, revised third
Whiting MJ, Nagy KA, Bateman PW: Evolution and maintenance of social status-signaling badges -Experimental manipulations in lizards. Lizard social behavior. Edited by: Fox SF, Baird TA. 2003, The Johns Hopkins University Press, Baltimore and London
Bokony V, Lendvai AZ, Liker A: Multiple cues in status signalling: the role of wingbars in aggressive interactions of male house sparrows. Ethology. 2006, 112: 947-954. 10.1111/j.1439-0310.2006.01246.x.
Martin J, Moreira PL, Lopez P: Status-signalling chemical badges in male Iberian rock lizards. Func Ecol. 2007, 21: 568-576. 10.1111/j.1365-2435.2007.01262.x.
Clutton-Brock TH, Albon SD: The roaring of red deer and evolution of honest advertisement. Behaviour. 1979, 69: 145-169. 10.1163/156853979X00449.
Hammerschmidt K, Todt D: Individual differences in vocalizations of young Barbary macaques (Macaca sylvanus): A multiparametric analysis to identify critical cues in acoustic signalling. Behaviour. 1995, 132: 381-399. 10.1163/156853995X00621.
Rendall D, Rodman PS, Emond RE: Vocal recognition of individuals and kin in free-ranging rhesus monkeys. Anim Behav. 1996, 51: 1007-1015. 10.1006/anbe.1996.0103.
Welch AM, Semlitsch RD, Gerhardt HC: Call duration as an indicator of genetic quality in male gray tree frogs. Science. 1998, 280: 1928-1930. 10.1126/science.280.5371.1928.
Drickamer LC: Urine marking and social dominance in male house mice (Mus musculus domesticus). Behav Process. 2001, 53: 113-120. 10.1016/S0376-6357(00)00152-2.
Fichtel C, Hammerschmidt K, Jurgens U: On the vocal expression of emotion. A multi-parametric analysis of different states of aversion in the squirrel monkey. Behaviour. 2001, 138: 97-116. 10.1163/15685390151067094.
Rostain RR, Ben-David M, Groves P, Rabdall JA: Why do river otters scent-mark? An experimental test of several hypotheses. Anim Behav. 2004, 68: 703-711. 10.1016/j.anbehav.2003.10.027.
Hirano ZMB, Correa IC, DE Oliveira DAG: Contexts of rubbing behavior in Alouatta guariba clamitans: a scent-marking role?. Amer J Primat. 2008, 70: 575-583. 10.1002/ajp.20531.
Shipley C, Hines M, Buchwald JS: Individual differences in threat calls of northern elephant seal bulls. Anim Behav. 1981, 29: 12-19. 10.1016/S0003-3472(81)80147-9.
Brown RE, Macdonald DW: Social Odours in Mammals. 1985, Oxford, Claredon
Insley SJ, Phillips AV, Charrier I: A review of social recognition in pinnipeds. Aqua Mam. 2003, 29: 181-201. 10.1578/016754203101024149.
Bartholomew GA, Collias NE: The role of vocalization in the social behaviour of the northern elephant seal. Anim Behav. 1962, 10: 7-14. 10.1016/0003-3472(62)90124-0.
McComb K: Female choice for high roaring rates in red deed, Cervus elaphus . Anim Behav. 1991, 41: 79-88. 10.1016/S0003-3472(05)80504-4.
Fitch WT: Vocal tract length and formant frequency dispersion correlate with body size in the rhesus macaques. J Acoust Soc Am. 1997, 102: 1213-1222. 10.1121/1.421048.
Bergman TJ, Beehner JC, Cheney DL, Seyfarth RM: Hierarchical classification by rank and kinship in baboons. Science. 2003, 302: 1234-1236. 10.1126/science.1087513.
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. 10.1006/anbe.2003.2078.
Frey R, Gebler A, Fritsch G: Arctic roars - laryngeal anatomy and vocalization of the muskox (Ovibos moschatus Zimmermann, 1780, Bovidae). J Zool. 2006, 268: 433-448. 10.1111/j.1469-7998.2006.00053.x.
Pfefferle D, Fisher J: Sounds and size: identification of acoustic variables that reflect body size in hamadryas baboons, Papio hamadryas . Anim Behav. 2006, 72: 43-51. 10.1016/j.anbehav.2005.08.021.
Koren L, Mokady O, Geffen E: Social status and cortisol levels in singing rock hyraxes. Horm Behav. 2008, 54: 212-216. 10.1016/j.yhbeh.2008.02.020.
McComb K, Reby D, Baker L, Moss C, Sayialel S: Long-distance communication of acoustic cues to social identity in African elephants. Anim Behav. 2003, 65: 317-329. 10.1006/anbe.2003.2047.
Hollister-Smith JA, Alberts SC, Rasmussen LEL: Do male African elephants, Loxodonta africana, signal musth via urine dribbling?. Anim Behav. 2008, 76: 1829-1841. 10.1016/j.anbehav.2008.05.033.
Kruuk H: The Spotted Hyena. A study of predation and social behavior. 1972, Chicago, Univ Chicago Press
Franck LG: Social organization of the spotted hyaena (Crocuta crocuta): II. Dominance and reproduction. Anim Behav. 1986, 35: 1510-1527. 10.1016/S0003-3472(86)80221-4.
Glickman SE, Cunha GR, Drea CM, Conley AJ, Place NJ: Mammalian sexual differentiation: lessons from the spotted hyena. Tr Endocr Metab. 2006, 17: 349-356. 10.1016/j.tem.2006.09.005.
Holekamp KE, Smale L: Rank acquisition during mammalian social development: the 'inheritance' of maternal rank. Amer Zool. 1991, 31: 306-317.
Jenks S, Weldele M, Frank L, Glickman SE: Acquisition of matrilineal rank in captive spotted hyaenas: emergence of a natural social system in peer-reared animals and their offspring. Anim Behav. 1995, 50: 893-904. 10.1016/0003-3472(95)80092-1.
Engh AL, Esch K, Smale L, Holekamp KE: Mechanisms of maternal rank 'inheritance' in the spotted hyaena, Crocuta crocuta . Anim Behav. 2000, 60: 323-332. 10.1006/anbe.2000.1502.
Smale L, Frank LG, Holekamp KE: Ontogeny of dominance in free-living spotted hyenas: juvenile rank relations with adults. Anim Behav. 1993, 46: 467-477. 10.1006/anbe.1993.1215.
East ML, Hofer H: Male spotted hyenas (Crocuta crocuta) queue for status in social groups dominated by females. Behav Ecol. 2001, 12: 558-568. 10.1093/beheco/12.5.558.
Holekamp KE: Spotted hyenas. Cur Biol. 2006, 16: R944-10.1016/j.cub.2006.10.009.
Holekamp K, Sakal ST, Lundrigan BL: Social intelligence in the spotted hyena (Crocuta crocuta). Phil Tr R Soc B. 2007, 362: 523-538. 10.1098/rstb.2006.1993.
Mills MGL: Kalahari Hyaenas: the Behavioural Ecology of Two Species. 1990, London, Unwin Hyman
Kolowski JM, Katan D, Theis KR, Holekamp KE: Daily pattern of activity in the spotted hyena. J Mammal. 2007, 88: 1017-1028. 10.1644/06-MAMM-A-143R.1.
East ML, Hofer H: Crocuta crocuta. 1991, 42: 637-649.
East ML, Hofer H: Loud calling in a female-dominated mammalian society. II. Behavioral contexts and functions of whooping of spotted hyaenas, Crocuta crocuta . Anim Behav. 1991, 42: 651-669. 10.1016/S0003-3472(05)80247-7.
Theis KR, Greene KM, Benson-Amram SR, Holekamp KE: Sources of variation in the long-distance vocalizations of spotted hyenas. Behaviour. 2007, 144: 557-584. 10.1163/156853907780713046.
Glickman SE, Short RV, Renfree MB: Sexual differenciation in three unconventional mammals: Spotted hyenas, elephants and tammar wallabies. Horm Behav. 2005, 48: 403-417. 10.1016/j.yhbeh.2005.07.013.
Drea CM, Weldele ML, Forger NG, Coscia EM, Frank LG, Licht P, Glickman SE: Androgens and masculinization of genitalia in the spotted hyaena (Crocuta crocuta). 2. Effects of prenatal anti-androgens. J Reprod Fertil. 1998, 113: 117-127.
Place NJ, Glickman SE: Masculinization of female mammals: Lessons from nature. Hypospadias and Genital Development. Edited by: Baskin. 2004, Plenum Press, New York
Glickman SE, Frank LG, Pavgi S, Licht P: Hormonal correlates of 'masculinization' in female spotted hyaenas (Crocuta crocuta): I. Infancy through sexual maturity. J Reprod Fertil. 1992, 95: 451-462.
Place NJ, Holekamp KE, Sisk CL, Weldele ML, Coscia EM, Drea CM, Glickman SE: Effects of treatment with anti-androgens on LH secretion and gonadal steroidgenesis in the spotted hyena, Crocuta crocuta . Biol Reprod. 2002, 67: 1405-1413. 10.1095/biolreprod.102.004226.
Welch PD: The Use of Fast Fourier Transform for the Estimation of Power Spectra: A Method Based on Time Averaging Over Short, Modified Periodograms. IEEE Trans Audio Electroacoustics. 1967, 15: 70-73. 10.1109/TAU.1967.1161901.
Sueur J, Aubin T, Simonis C: Seewave, a free modular tool for sound analysis and synthesis. BioAcoustics. 2008, 18: 213-226.
Krzanowski WJ: Principles of multivariate analysis: A user's perspective. 1988, New-York, Oxford University Press
Shannon CE, Weaver W: The mathematical theory of communication. 1963, Chicago, University of Illinois Press
Mundry R, Sommer C: Discriminant function analysis with nonindependent data: consequences and an alternative. Anim Behav. 2007, 74: 965-976. 10.1016/j.anbehav.2006.12.028.
Benjamini Y, Hochberg Y: Controlling the false discovery rate: A practical and powerful approach to multiple testing. J R Stat Soc B. 1995, 57: 289-300.
Aubin T, Jouventin P: How to identify vocally a kin in a crowd? The penguin model. Adv Stud Behav. 2002, 31: 243-277. full_text.
Aubin T, Mathevon N, Staszewski V, Boulinier T: Acoustic communication in the Kittiwake Rissa tridactyla: potential cues for sexual and individual signatures in long calls. Polar Biol. 2007, 30: 1027-1033. 10.1007/s00300-007-0262-6.
Curé C, Aubin T, Mathevon N: Acoustic convergence and divergence in two sympatric burrowing nocturnal seabirds. Biol J Linn Soc. 2009, 96: 115-134. 10.1111/j.1095-8312.2008.01104.x.
Blumstein DT, Munos O: Individual, age and sex-specific information is contained in yellow-bellied marmot alarm calls. Anim Behav. 2005, 69: 353-361. 10.1016/j.anbehav.2004.10.001.
Phillips AV, Stirling : Vocal individuality in mother and pup South American fur seals, Arctocephalus australis. Mar Mam Sc. 2000, 16: 592-616. 10.1111/j.1748-7692.2000.tb00954.x.
Charrier I, Mathevon N, Jouventin P: How does a fur seal mother recognize the voice of her pup? An experimental study of Arctocephalus tropicalis . J Exp Biol. 2002, 205: 603-612.
Peters G, East ML, Hersel H, Henschel JR, Mills MGL, Wilhelm K, Hofer H: Spotted hyaens whoops: frequent incidence of vocal instabilities in a mammalian loud call. BioAcoustics. 2004, 14: 99-109.
Cooper SM: Optimal hunting group size: the need for lions to defend their kills against loss to spotted hyaenas. Afr J Ecol. 1991, 29: 130-136. 10.1111/j.1365-2028.1991.tb00993.x.
Holekamp KE, Boydston EE, Szykman M, Graham I, Nutt KJ, Birch S, Piskiel A, Singh M: Vocal recognition in the spotted hyaena and its possible implications regarding the evolution of intelligence. Anim Behav. 1999, 58: 383-395. 10.1006/anbe.1999.1157.
Holekamp K, Sakal ST, Lundrigan BL: The spotted hyena (Crocuta crocuta) as a model system for a study of the evolution of intelligence. J Mammal. 2007, 88: 545-554. 10.1644/06-MAMM-S-361R1.1.
McGregor P: Animal communication networks. 2005, Cambridge, Cambridge University Press
Mathevon N, Aubin T, Vielliard J, Da Silva ML, Sebe F, Boscolo D: Singing in the rain forest: How a tropical bird song transfers information. PloS ONE. 2008, 3 (2): e1580-10.1371/journal.pone.0001580.
Engh AL, Siebert ER, Greenberg A, Holekamp KE: Patterns of alliance formation and post-conflict aggression indicate spotted hyenas recognize third party relationships. Anim Behav. 2005, 69: 209-217. 10.1016/j.anbehav.2004.04.013.
Tomasello M, Call J: Primate cognition. 1997, Oxford, Oxford University Press
Vignal C, Mathevon N, Mottin S: Audience drives male songbird response to partner's voice. Nature. 2004, 430: 448-451. 10.1038/nature02645.
Drea CM, Vignieri SN, Kim S, Weldele ML, Glickman SE: Responses to olfactory stimuli in spotted hyenas (Crocuta crocuta): II. Discrimination of conspecific scent. J Comp Psychol. 2002, 116: 342-349. 10.1037/0735-7036.116.4.342.
Watts HE, Holekamp KE: Hyena societies. Cur Biol. 2007, 17: R657-660. 10.1016/j.cub.2007.06.002.
We thank Roger Mundry for his pDFA R routine and four anonymous referees for their careful review and their constructive criticisms. We are grateful to the Institut universitaire de France (NM), and the Miller Institute at UC Berkeley for their support. This study has been performed during a sabbatical leave of NM at the University of California, Berkeley (Visiting Miller Professorship funded by the Miller Institute for Basic Research in Science). During the period when this research was conducted, the Berkeley hyena colony was supported by NSF Grant SGER-0809914 and NSF Grant IOS-0920793.
NM and FT conceived and coordinated the study. All authors performed the work. NM and FT drafted the manuscript. All authors read and approved the final manuscript.
Nicolas Mathevon and Frédéric E Theunissen contributed equally to this work.