Figures
Abstract
Animal vocal signals may provide information about senders and mediate important social interactions like sexual competition, territory maintenance and mate selection. Hence, it is important to understand whether vocal signals provide accurate information about animal attributes or status. Gibbons are non-human primates that produce loud, distinctive and melodious vocalizations resembling more those of birds than of other non-human primates. Wild gibbons are characterized by flexibility in social organization (i.e., pairs and multimale units) as well as in mating system (i.e., monogamy and polyandry). Such features make them a suitable model to investigate whether the physiology (hormonal status) and socio-demographic features find their correspondence in the structure of their songs. By combining male solo song recordings, endocrine outputs using non-invasive fecal androgen measures and behavioral observations, we studied 14 groups (10 pair-living, 4 multimale) of wild white-handed gibbons (Hylobates lar) residing at Khao Yai National Park, Thailand. We collected a total of 322 fecal samples and recorded 48 songs from 18 adult animals. Our results confirmed inter-individuality in male gibbon songs, and showed a clear correlation between androgen levels and song structures. Gibbons with higher androgen levels produced calls having higher pitch, and similarly adult individuals produced longer calls than senior males. Thus, it is plausible that gibbon vocalizations provide receivers with information about singers' attributes.
Citation: Barelli C, Mundry R, Heistermann M, Hammerschmidt K (2013) Cues to Androgens and Quality in Male Gibbon Songs. PLoS ONE 8(12): e82748. https://doi.org/10.1371/journal.pone.0082748
Editor: David Reby, University of Sussex, United Kingdom
Received: September 12, 2013; Accepted: October 28, 2013; Published: December 18, 2013
Copyright: © 2013 Barelli et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This research was supported by grants from Alexander von Humboldt Foundation (post-doctoral fellowship to CB) and the German Research Foundation (grant HE2699/6-1 to MH). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Vocalizations are assumed to provide information about senders, including their identity [1]–[3], rank [4], [5], age, sex or size [6], [7]. Given the importance of vocal signals in mediating crucial social interactions (i.e., sexual competition, territorial maintenance, partner or parent/young recognition [8]) in birds, anurans and mammals vocal signals (e.g., duration and fundamental frequency) should provide accurate information about animal attributes or status. They may also convey honest information because only individuals in better condition should be more capable to afford any associated costs with signal production.
In birds, hormones play a central role in singing behavior. Studies on birds of the temperate-zone have shown that incidence and complexity of songs are closely related to changes in androgen (i.e. testosterone) levels [9], [10]. These behavioral changes find their correlates in neuronal structural changes of birds' brains, which is one of the most spectacular examples of neuroplasticity [11]. Apart from the endocrine component, age, social status and reproductive success have also been reported to be correlated with vocal performance [12]–[14].
Similar to birds, mammals' vocal displays are also often mediated by androgens [15], [16]. These hormones affect vocalizations through actions on motivational centers and vocal motor pathways in the central nervous system [16] or via modulation of peripheral structures involved in signal production. For example, changes in androgens may influence the pitch of vocalizations in giant pandas (Ailuropoda melaneuca: [17]) or affect the trill rate in Neotropical singing mice (genus Scotinomys: [18]). The influence of androgens on the structure of mammalian vocalizations and pronunciation [19], [20] is relatively well documented for humans [20]–[22], but in non-human primates no direct evidence has yet been found (baboons: [23], [24]). Age, possibly due to age-related vocal fold growth, as well as social status may have an influence on fundamental frequency (F0; e.g., humans [25], [26]; baboons [4], [6], [27]; chimpanzees [28]). However, there is evidence for a high level of independence between F0 and age (e.g., humans [29]; non-human primates [30]).
Based on the hypothesis that vocal fold morphology is not exclusively relevant to vocal differences [31], [32], vocal folds' mechanical properties may shape specific vocal signals such as F0 [32], [33]. Moreover, androgen levels, as well as other features (e.g., age, social status), may directly act on laryngeal muscles and connective tissue of vocal folds by constantly remodeling them [34], [35]. Changes in hormonal levels have an obvious physiological explanation during adolescence when androgen variations have a permanent impact in the length or tension of the vocal folds [32]. However, after maturation, androgens may still vary and fluctuate in a reversible manner. Androgen changes may thus have also a psychological component affecting the vocal production which an individual may use as part of a social interaction strategy (e.g., humans [21]).
All in all, it is plausible that inter-individual variation of vocal signals and their related components in anatomy, physiology and behavior can lead to vocal differences between individuals and identify individuals, which in turn may convey information to different types of receivers, for example, to potential competitors (humans [26]; red deer [36]; domestic dogs [37]) or to potential mates (humans [38]–[40]; other mammals [41]).
Among non-human primates, gibbons are of particular interest since they produce distinctive vocalizations (‘song’) which are species-specific [42], [43]. Male songs (‘solos’), which occur in addition to gibbons' well known female-male duets, are loud enough to be heard up to a kilometer away [44] and can last for up to 4 hours [45]. Gibbon solo songs may function in one or more of several ways, including home range defense against neighboring males and solitary conspecific strangers as well as communicating with candidate mates and strengthening pair bonds [44], [46]–[50]. As a consequence, gibbon solo songs have developed clear adaptations to improve long-distance transmission [51]–[54] making them very different from other non-human primate vocalizations, but resembling those of rainforest birds.
Despite the existence of some basic information about male gibbon vocalizations, to date it is completely unknown whether individual attributes of the caller are encoded in the acoustic structure of their songs. Thus, we investigated the wild white-handed gibbon (Hylobates lar) population residing at Khao Yai National Park, Thailand, which is characterized by flexibility in both social organization (i.e., single-male/single-female groups and groups with one female and more than one adult male [55], [56]) and mating system (monogamy and polyandry occur, as well as extra-pair copulations and conceptions [55], [57], [58]). Moreover, a recent study has revealed a close association between androgen levels and social organization, with higher androgen concentrations being found in males living in pairs rather than those living in groups with more than one male [59]. All these features make this population particularly suitable to examine whether male solo songs' structure and its specific acoustic parameters are related to physiological status (i.e., androgen levels), as well as to socio-demographic features (i.e., social status or age). We first assess (i) inter-individual differences in song structure between gibbon males. Later, we combine male solo song recordings and endocrine outputs using non-invasive fecal androgen measures and behavioral observations. We describe and identify the relationship between song structure (e.g., F0 and duration) and (ii) androgen levels, (iii) social status, and (iv) age. Specifically, considering the spectral domain, we expect to find some sort of association between androgens and pitch, whereas the direction is difficult to predict. Regarding the temporal domain, a positive relationship between call duration and androgens would be expected. Considering the different social units (i.e., pairs and unifemale/multimale units) and previous findings that males living in pairs showed higher androgen levels than those living in groups with more than one male, we also predict that pitch varies between different types of males (i.e., the only male in pairs, primary or secondary males in unifemale/multimale units).
Materials and Methods
Ethics Statement
All fieldwork was performed without direct contact or interaction with the individuals and under research/collecting permit number 2.3/2210 to CB, issued by the National Research Council of Thailand (NRCT) and the National Park, Wildlife and Plant Conservation Department (DNP) of Thailand. Fecal material was imported into Germany under the permit V3/19j 06.23 issued from the Hessian Ministry for Environment, Energy, Agriculture and Consumerism. A CITES permit was not required as animal excretory products (urine, feces) do not fall under CITES regulations.
Gibbon population and study site
At time of data collection, the study population residing in the Mo Singto-Klong E-Tau area, located in the central part of the Khao Yai National Park, Thailand (2,168 km2; 14°26′ N, 101°22′ E; ∼150 km NE of Bangkok), hosted a total of 24 males living in 14 well-known habituated groups. Khao Yai National Park consists of a sandstone plateau ranging between 600 m and 1000 m above sea level and covered by seasonal evergreen forest [60]. The study area is a continuous forest on a hilly terrain of 8.5 km2 which hosts approximately 15.9 individuals of white-handed gibbons per km2 [57]. Group size ranges between two and six individuals, usually consisting of a mated pair with their putative offspring, or multimale units with a single breeding female and two or more sexually mature males [57], [61]. Group home ranges tend to be of approximately 400–500 meters in diameter. Out of the 24 males available, we were successful in recording solo songs from 18 animals (three males performed only duets, while for 3 others acoustic data were not available; Table 1). Among those, 14 males were considered adults, both primary and secondary males (for definition see below), and four were classified as subadults, which were fully grown males but still residing in their natal group (Table 1).
Social organization, social status and age class
Due to the flexible social organization of the Khao Yai gibbons [55], [56], during the study period 10 of the 14 focal groups consisted of a single adult female, a single adult male and up to three offspring, reflecting a pair-living social organization (cf. [62]). The remaining four groups were considered unifemale/multimale units, with a single adult female and at least two adult males unrelated to the resident female (Table 1; [55], [58], [59]).
We classified all study males as either 'primary' or 'secondary' [54], [57]. Primary males were the only males in the pair-living groups or those males in unifemale/multimale groups [55] which predominantly engaged in singing duets with the adult female and performed the majority of copulations with them (N = 9). All remaining males in unifemale/multimale units, who rarely sang or copulated with the group's adult female [55] were considered as secondary (N = 5).
To test predictions about the relationship between age and acoustic features, we classified males into two age categories: adult (8–25 years of age) and senior (exceeding 25 years of age [59]). Since a few males were mature in age (i.e., falling in the adult category) but still residing in their natal group (i.e., falling in the subadult category), we decided to exclude them from analysis and present only descriptive information.
Fecal sampling and androgen analysis
To examine the correlation between male androgen status and acoustic parameters, we determined fecal androgen metabolite levels for each adult male sampled and related them to their vocal structure. While recording vocalizations, three teams of field assistants collected fecal samples from target animals (all individuals are well known since decades and recognized by individual natural markers [57]) for endocrine analysis of 18 males (including males classified as subadults). Each team followed four to five of the 14 study groups from dawn to dusk (mean observation time: 8 h/day) for three months each. All samples were collected between 6:30 and 14:00 directly following defecation (for details see [59]). A total of 322 samples uncontaminated by urine (15 samples per male; range 5–29) were collected, kept on ice until arrival at the field station and then stored frozen until transport to the endocrinology laboratory of the German Primate Center (DPZ) for measurement of androgen content. Specifically, samples were analyzed for immunoreactive epiandrosterone (EA), a major metabolite of testosterone in primate feces [63]–[65], using an enzyme immunoassay (EIA) recently validated for monitoring androgen output in the white-handed gibbon [59]. Fecal extractions and assay procedures were carried out as previously described [59], [63]. Intra-assay coefficients of variation (CV) for a high- and low-concentrated quality control were 7.1% and 8.5%, respectively. Corresponding figures for inter-assay CVs were 12.7% and 15.9%.
Acoustic recordings and analysis
We collected most of the acoustic data within a period of six months, between October 2008 and March 2009, with a second set of recordings being undertaken between April and May 2010. We recorded gibbon vocalizations ad libitum using a Sennheiser directional microphone (K6 power module and ME66 recording head with MZW66 pro windscreen) and two Marantz solid state recorders (PMD 660 and 670). We daily followed focal groups from dawn till dusk (average 8 hrs/day) and, whenever a male started singing, we recorded his vocalization within a distance of 5–20 meters. Information regarding subject identity and context was always spoken onto the tape or noted down into spreadsheets. Sounds were recorded in mono format with 16-bit resolution and 44.1-kHz sampling rate.
Vocalizations were characterized by a number of structural and temporal parameters. We included temporal measurements because changes in androgen levels could also lead to motivational changes which likely influence the temporal structure of primate vocalization. We defined as 'element' the single note uttered by a singing individual, while a sequence of undefined number of elements, separated by a short interval of time between each other, was classified as 'call'. Combinations of call sequences identified male 'song' for each individual gibbon (Fig. 1). To obtain an adequate frequency resolution, we down-sampled files from 44.1 kHz to 8 kHz. By using SASLab Pro 5.1 (Avisoft Bioacoustics, Berlin, Germany), we estimated several parameters describing the frequency modulation of F0 which in gibbons is the frequency with the highest amplitude [66], [67]. We used the automatic parameter measurement tool to extract acoustic parameters from spectrograms (FFT length = 256, frequency resolution = 31 Hz, temporal resolution = 16 ms (overlap = 50%), window type = Hamming). For each element we measured: (i) the initial peak of fundamental frequency (defined as 'start F0'), (ii) the final peak of fundamental frequency (end F0) and (iii) the maximum peak of fundamental frequency (max F0). In addition, we calculated three temporal measures: (iv) duration (in seconds) of each element from the initial to the final F0, (v) duration (in seconds) between consecutive elements, and finally (vi) the temporal location (in seconds) of max F0 divided by the element duration (Fig. 1). Depending on the background noise we used a flexible threshold (ranging between −5 and −20 dB, mean value: 12.8) to distinguish between noise and signal. We combined the frequency measurements per call element to characterize changes at the call level. Beside mean values per element, we also included maximum of a call and variation within a call to account for variability between call elements. Together with call duration we had 22 acoustic parameters to characterize the gibbon calls in frequency and temporal domain (Table 2). For the 14 animals included into the acoustic analysis, we recorded a total of 48 songs, 784 calls and 3,993 elements.
Statistical analysis
Factor analysis.
To remove redundancy between the acoustic parameters we first ran a Factor Analysis (FA) on parameters derived from calls. This approach was justified as indicated by large correlations between the acoustic parameters, Bartlet's test of sphericity (χ2 = 30707, df = 231; P<0.001) and the Kaiser-Meier-Olkin measure of sampling adequacy (0.711 [68]). Before running the FA, and in order to achieve an approximately symmetrical distribution, we checked the distribution of each of the 22 acoustic parameters and transformed variables when required (Table 2). The FA was run with varimax rotation, and we used the regression method to obtain scores for each of the factors derived. In the subsequent analyses we used the derived factor scores as measures of the acoustic properties of the songs recorded.
Inter-individuality.
We used a Discriminant Function Analysis (DFA) to test for differences between calls of different individuals, and a permuted DFA to account for non-independence of calls recorded at the same day (pDFA [69]). Moreover, since vocal recordings of the same individual were sometimes collected during the same days, we also permuted calls day-wise between subjects. The DFA included a total of 10 individuals for which at least 23 calls were recorded (i.e., at least one more than the number of the acoustic parameters: [N = 22]). Prior to the DFA, we transformed variables as described above (Table 2). The total number of calls included in this analysis was 647.
To derive the discriminant functions and to balance the individual contribution, we used 23 randomly selected calls per individual, while to reduce the impact of any random selection we ran 100 random selections and averaged the results. P-values of the DFA were based on cross-validated calls and determined using 1,000 permutations into which the original data were included as one permutation. The pDFA was conducted in R using a script written by one of us (RM). To estimate the contribution of the individual acoustic parameters to the discriminability between males, we ran a DFA in SPSS (version 15) including all calls.
Relationship between acoustic parameters, socio-demographic features and fecal androgens.
To investigate whether the acoustic structure of the call parameters (i.e., the factor scores derived from the FA) varied according to socio demographic features (i.e., social organization, social status and age) and male androgen levels, we used General Linear Mixed Models (GLMM [70]). We ran six separate models, each with one of the six factor scores describing the acoustic features of the calls as the response. We included into these models male status (pair-living male; primary male in uni-female/multimale units; secondary male in uni-female/multimale units), age (adult; senior) and fecal androgen level (see below) as fixed effects while group, subject identity, date (nested with subject) and song were included as random effects. To test whether variation in call parameters due to varying androgen levels was happening at the level of between subjects variation (i.e., effects of androgen being largely a function of overall differences between subjects with regard to their average androgen levels) and/or within subjects (i.e., call parameters varying as a functions of short term variations of androgen levels within subjects) we used within subjects centering [71]. More precisely, male androgen levels were represented in the models by two terms: one being the average androgen level per subject and one being the actual androgen values, centered to a mean of zero per subject (by subtracting from each value the mean androgen level of the respective subject). To control the possibility that the effect of varying androgen levels on the acoustic structure of the calls varied between males we also included a random slope component of (within subjects centerd) androgen level within males [72].
Due to the presumed time lag for metabolites excretion into feces in gibbons [57], [73], we considered day 3 after the vocal recording to be the day at which the fecal sample best reflects the androgen level at the day of recording. Thus, as a measure of androgen level we used values from fecal sample collected closest to this optimal day, whereby we considered only fecal samples which were collected between the day of recording and 7 days later. When several samples fulfilled this criterion we averaged the values. In total, we had 295 samples collected at the optimal day, 62 and 64 collected at the day of recording or the next day, respectively, and 3, 5, 22, 52 at days 4 to 7 after the recordings. Prior to the analysis, we z-transformed values of fecal androgen levels to a mean of zero and a standard deviation of one. We included only one type of song, male solos, and only adult males (we excluded subadults from the model), leaving to a total of 503 calls from 34 songs recorded on 24 days from 13 adult animals out of 10 groups.
The models were fitted using Gaussian error function and identity link. We checked for the assumptions of normally distributed and homogenous residuals by visually inspecting histograms and qq-plots of the residuals as well as residuals plotted against fitted values. None of these indicated severe deviations from these assumptions (assumptions were checked only after the autocorrelation term had been included; see below). Estimating the significance of fixed effects in mixed models is controversial [74]. Here we estimated P-values using Markov-chain Monte-Carlo (MCMC) analysis, presumably the most reliable method currently available [70].
Song structure was likely to show temporal autocorrelation (i.e., calls recorded closer to one another in time being more similar to one another than those recorded more distantly). Such temporal autocorrelation may lead to non-independent residuals and could potentially greatly devalue the validity of the model. Hence, we incorporated autocorrelation into the model by first running the full model (as described above) and retrieving the residuals from it. Subsequently, and separately for each data point, we calculated a weighted average of the residuals of all other data points from the same male with the weight of the residuals being proportionate to their time lag to the specific data point. The resulting variable, 'autocorrelation term', was then included as an additional fixed effect. The weight function had the shape of a Gaussian function with a mean at a time lag equal to zero. Its standard deviation was determined such that the likelihood of the full model with the autocorrelation term included was maximized. Time was measured on a continuous scale considering the actual time each particular call began. To achieve an easy and interpretable estimate we z-transformed the autocorrelation term to a mean of zero and standard deviation of one before including it in the model.
We tested for model stability, by excluding subjects one by one and comparing the estimates. This showed instability issues in some of them, but given the non significance of most of the estimates, this result could be expected. In fact, the only term that remained significant after the correction for multiple testing showed little variations in the revealed P-values.
Since testing the impact of fixed effects (fecal androgen level, social organization, social status and age) on the six factors required accounting for multiple testing. We hence applied Simes' method [75] for this purpose.
The FA and a single DFA were run in SPSS (15.0). All other analyses were conducted in R (2.15.2 [76]). GLMMs were calculated using the function lmer of the R package lme4 [77], MCMC P-values were derived using the function pvals.fnc of the R package languageR [78], the pDFA was based on the function lda of the R-package MASS [79], and the autocorrelation term was derived using a self-written function.
Results
Factor analysis
The FA revealed six factors with Eigenvalues >1 together explaining 82.4% of the total variance (Table 2). Based on the loadings of acoustic parameters on the rotated components we were able to characterize and label these factors.
All acoustic parameters describing the element pitch showed high loadings on Factor 1 (pitch). Factor 2 (element duration) showed high loadings of parameters describing mean and maximum duration of elements. Acoustic parameters which showed high loadings on Factor 3 (pitch variation) described the variation in F0 of elements. The three high loadings on Factor 4 (location of max F0) described mean, maximum and variation of the location of maximum F0. Total call duration loaded exclusively on Factor 5 (call duration) which, in addition, showed stronger loadings of acoustic parameters describing interval duration between elements. Factor 6 (start F0) showed only high loadings of parameters describing the start F0.
Inter-individual differences
The pDFA revealed that calls differed between individuals (average percentage of correctly assigned cross classified calls: 74.6%, chance level = 10%, P = 0.001; Fig. 2). Running a DFA on all calls revealed two discriminant functions with Eigenvalues >1. Variables with high absolute loadings (≥0.5) on any of the first two discriminant functions were (i) maximum element duration, (ii) variation in element duration and (iii) maximum of mean F0.
Full names of gibbons recorded are reported at the x-axis as well as total number of calls used in the analysis on top of the graph. Males are represented in order of increasing mean androgen levels. Indicated are medians (horizontal lines) and quartiles (vertical lines).
Androgens, social status and age on gibbon calls
We found a clear link between fecal androgen levels and Factor 1 'element F0' (Table 3). Gibbons with higher average androgen levels produced calls having call elements with significantly higher pitch (Factor 1; Table 3; see also Appendix, Table S1). None of the other acoustics properties tested co-varied with androgen levels.
We also found that among adult males those of senior age had lower call duration (Factor 5; Table 3; Appendix, Table I). No obvious relation among any of the remaining call parameters considered was found between males belonging to different social status (Table 3).
Although only qualitative data were available, subadults (males already mature but still residing in their natal groups) presented interesting similarities to senior males (i.e., number of elements per call, number of call per song, start and maximum F0; Table 4). Indeed subadults differed from anybody else in call duration, duration of intervals between elements and element duration (Table 4).
Discussion
Our study aimed to investigate wild white-handed male gibbon solo songs with respect to individuality, hormonal underpinning and relationship to socio-demographic features such as social status and age. First, we confirm that male gibbon songs exhibit significant differences among individuals and such variation is expressed in terms of song characteristics. Individual differences are common in a variety of other primate vocalizations, including gibbons (male songs [54]; female songs [80]), baboons [81], [82], and chimpanzees [83], as well as in other animal species [84]–[86]. Many playback studies have shown that listeners can use this information to distinguish between group members [4], [87], [88].
Our results show a highly significant correlation between androgen levels and vocal pitch. Males with higher androgen levels produced elements having higher pitch. The relationship found between male androgen levels and their vocal parameters shows the opposite of what is known about the influence of testosterone on the male voice in humans [21], [89]. It is known that during puberty, elevated testosterone may act through androgen receptors on vocal folds and causes their growth [90], [91]. By lengthening and thickening the vocal folds, fundamental frequency becomes lower [92] and, as a consequence, men vocalize with lower pitch, approximately half of the F0 in women [93]. However, in studies of non-human primate vocalizations, data from two studies on chacma baboons indirectly suggest a positive relationship between male testosterone levels and vocal pitch [23], [24]. While running around or leaping through the trees, male chacma baboons give so-called 'wahoos' which are loud and highly costly calls. Fischer and colleagues [4] showed that high ranking males produce 'wahoos' with higher F0 and longer 'hoo' elements than low ranking males. When males fall in rank the 'hoo' syllables become shorter and F0 declines. In a subsequent study, measuring testosterone levels of chacma baboon males, it has been shown that high ranks are strongly correlated with higher testosterone levels and that with a decline in rank also testosterone levels decrease [23], [24]. Thus, although indirectly, these two studies strongly suggest that in chacma baboons an increase in testosterone is likely associated with an increase in F0, as observed in male gibbons of this study. Differently from chacma baboons, white-handed gibbons do not show an increase in element duration, but a significant increase in call pitch. Another study on male loud call characteristics of Thomas langurs (Presbytis thomasi) revealed an increase in tonal units and duration with increasing testosterone levels [94]. To our knowledge, those are the only studies of non-human primates which focused on a direct relation between androgen levels and vocal structure. Although our data and those of others are correlational, one possible interpretation is that androgen levels lead to a higher motivation to sing which may trigger structural changes in vocalizations, as it has been found in complex hyrax (Procavia capensis) songs [95]. Independent of the proximate mechanism, it is likely that group members use such structural differences to assess motivation and status of the caller.
Androgen receptors have been shown in laryngeal muscles and connective tissue of vocal folds as well as in the brain [34], [35]. Thus, it could be that higher androgen levels affect vocal folds and change the vibration characteristics although the excitation frequency remains the same. Chances in laryngeal muscles and connective tissue via androgen receptors are reversible to some extent and therefore in accordance with the observed reversed pattern for the relation between F0, androgen and age. In addition, the relations between vocal fold morphology and mechanical properties studied in several species can support the hypothesis that differences in mechanical properties can explain our observed vocal differences [96], [97]. However, it is also possible that elevated androgen levels could lead to a higher motivation to call. This higher motivation could lead to calls having higher amplitude and thus an increase in call pitch [98], [99]. All in all, these studies on male call structure in relation to androgen levels suggest that male vocal signals in gibbons are consistent with an effect of androgens and hence could function as a reliable signal of male competitive abilities.
Socio-demographic features did not seem to have any impact on song parameters, except for age. Although previous findings have revealed differences in androgen levels among males living in pairs than those living in unifemale/multimale units [59], no obvious differences were found in vocal signals between males living under different social conditions. Subadult and juveniles do more likely show high variation in androgens levels compared to adults; however, we did not consider them in the analysis. Thus, age appeared to have some sort of effect on vocal production also in a rather restricted narrow androgen concentration range (among adults of different age) when vocal anatomy maturation was certainly already completed. Assuming that singing is a costly signal difficult to produce and which may reflect relative quality of individuals, only physically fit males should be able to perform song for a longer period [100]. A comparison of call duration of senior and adult males showed that adult males sang longer than senior ones. Such finding, combined with the result that males having higher androgens sang with higher pitch voice, suggests that solo songs may potentially provide honest signals of male quality important in mate choice [101].
As in other animal species, male vocalizations could also help to estimate male fighting ability without engaging in direct contest for both mate partner and territory [102]. By identifying individuals' age, males should attempt to assess asymmetries in fighting ability before engaging in escalated and potentially costly fights [103]–[105]. However, since the type of song examined were not given during the context of group encounters, but in the early morning when no encounters took place, we can speculate that they could likely function in assessing more male quality rather than territorial defense.
Supporting Information
Table S1.
Results of the GLMMs with factor scores as the responses, and androgen, age, social status as predictors. The models account for androgen levels using two fixed effects, one accounting for varying androgen levels between subjects (average androgen levels per male) and one for the within subjects variation of androgen levels (androgen levels centered to a mean of zero per subject). Subadults were excluded from the data analyzed. Note that the P-values are not corrected for multiple testing.
https://doi.org/10.1371/journal.pone.0082748.s001
(DOC)
Acknowledgments
We are grateful to the National Research Council of Thailand (NRCT), the National Park, Wildlife and Plant Conservation Department (DNP), the superintendent of Khao Yai National Park of Thailand (Research Permit No: 2.3/2210) for permissions to conduct this research, U.H. Reichard for agreeing to conduct our project at the study site and proving data on age class, and K. Hodges in supporting the project and its development. We would like to thank the Thai field assistants, T. Wolf and several international volunteers for assistance in the field, as well as A. Heistermann and P. Kiesel for assistance in hormone analysis. Many thanks go to the editor D. Reby, as well as to C. Crockford, C. Neumann and B. Wheeler for commenting on previous versions of this manuscript and T. Riede for insightful comments that greatly improved the manuscript.
Author Contributions
Conceived and designed the experiments: CB KH. Analyzed the data: RM KH. Wrote the paper: CB KH. Performed the hormonal analysis: MH. Obtained permission for conducting research in Thailand: CB.
References
- 1. Semple S (2001) Individuality and male discrimination of female copulation calls in the yellow baboon. Anim Behav 61: 1023–1028.
- 2. Owren MJ, Rendall D (2003) Salience of caller identity in rhesus monkey (Macaca mulatta) coos and screams: perceptual experiments with human (Homo sapiens) listeners. J Comp Psychol 117: 380–390.
- 3. Price T, Arnold K, Zuberbühler K, Semple S (2009) Pyow but not hack calls of the male putty-nosed monkey (Cercopithcus nictitans) convey information about caller identity. Behaviour 146: 871–888.
- 4. Fischer J, Kitchen D, Seyfarth R, Cheney D (2004) Baboon loud calls advertise male quality: acoustic features and relation to rank, age, and exhaustion. Behav Ecol Sociobiol 56: 140–148.
- 5. Vannoni E, McElligott AG (2008) Low frequency groans indicate larger and more dominant fallow deer (Dama dama) males. PLoS ONE 3: e3113.
- 6. Pfefferle D, Fischer J (2006) Sounds and size: identification of acoustic variables that reflect body size in hamadryas baboons, Papio hamadryas. Anim Behav 72: 43–51.
- 7. Ey E, Hammerschmidt K, Seyfarth RM, Fischer J (2007) Age- and sex-related variations in clear calls of Papio ursinus. Int J Primatol 28: 947–960.
- 8.
Owings D, Morton E (1998) Animal vocal communication: a new approach. Cambridge: Cambridge University Press.
- 9.
Catchpole CK, Slater PJB (2008) Bird song. Cambridge, UK: Cambridge University Press.
- 10. Voigt C, Leitner S (2013) Testosterone-dependency of male solo song in a duetting songbird - Evidence from females. Horm Behav 63: 122–127.
- 11. Ball GF, Riters LV, Balthazart J (2002) Neuroendocrinology of song behavior and avian brain plasticity: multiple sites of action of sex steroid hormones. Front Neuroend 23: 137–178.
- 12. Ballentine B (2009) The ability to perform physically challenging songs predicts age and size in male swamp sparrows, Melospiza georgiana. An Behav 77: 973–978.
- 13. Botero CA, Rossman RJ, Caro LM, Stenzler LM, Lovette IJ, et al. (2009) Syllable type consistency is related to age, social status and reproductive success in the tropical mockingbird. An Behav 77: 701–706.
- 14. de Kort SR, Eldermire ERB, Valderrama S, Botero CA, Vehrencamp SL (2009) Trill consistency is an age-related assessment signal in banded wrens. Proc R Soc B 276: 2315–2321.
- 15. Moore FL, Boyd SK, Kelley DB (2005) Historical perspective: hormonal regulation of behaviors in amphibians. Horm Behav 48: 373–383.
- 16. Bass AH, Remage-Healey L (2008) Central pattern generators for social vocalization: Androgen-dependent neurophysiological mechanisms Horm Behav. 53: 659–672.
- 17. Charlton BD, Keating JL, Kersey D, Rengui L. Huang Y, et al. (2011) Vocal cues to male androgen levels in giant pandas. Biol Lett 23: 71–74.
- 18. Pasch BAS, George P, Campbell P, Phelps SM (2011) Androgen-dependent vocal performances influences female preference in Neotropical singing mice. Anim Behav 82: 177–183.
- 19. Doupe AJ, Kuhl PK (1999) Birdsong and human speech: common themes and mechanisms. Annu Rev Neurosci 22: 567–631.
- 20. Evans S, Neave N, Wakelin D, Hamilton C (2008) The relationship between testosterone and vocal frequencies in human males. Physiol Behav 93: 783–788.
- 21. Dabbs JM, Mallinger A (1999) High testosterone levels predict low voice pitch among men. Pers Indiv Diff 27: 801–804.
- 22. King A, Ashby J, Nelson C (2001) Effects of testosterone replacement on a male professional singer. J Voice 15: 553–557.
- 23. Kitchen DM, Cheney DL, Seyfart RM, Beehner J (2009) A link between fecal testosterone and an honest signal – the loud ‘wahoo’ vocalization of chacma baboons. Am J Phys Anthropol 48: 166.
- 24. Gesquiere LR, Learn NH, Simao MC, Onyango PO, Alberts SC, et al. (2011) Life at the top: rank and stress in wild male baboons. Science 333: 357–360.
- 25. Puts DA, Gaulin SJC, Verdolini K (2006) Dominance and the evolution of sexual dimorphism in human voice pitch. Evol Hum Behav 27: 283–296.
- 26. Puts D, Hodges C, Cardenas R, Gaulin S (2007) Men's voices as dominance signals: vocal fundamental and formant frequencies influence dominance attributions among men. Evol Hum Behav 28: 340–344.
- 27. Fischer J, Hammerschmidt K, Cheney D, Seyfarth R (2002) Acoustic features of male baboon loud calls: influences of context, age and individuality. J Acoust Soc Am 111: 1465–1474.
- 28. Riede T, Arcadi AC, Owren MJ (2007) Nonlinear acoustics in pant hoots and screams of common chimpanzees (Pan troglodytes): Vocalizing at the edge. J Acoust Soc Am 121: 1758–1767.
- 29. Collins SA (2000) Men's voices and women's choices. Anim Behav 60: 773–780.
- 30. Ey E, Pfefferle D, Fischer J (2007) Do age- and sex-related variations reliably reflect body size in non-human primate vocalizations? A review. Primates 48: 253–267.
- 31. Riede T, Titze IR (2008) Vocal fold elasticity of the Rocky Mountain elk (Cervus elaphus nelsoni) – producing high fundamental frequency vocalization with a very long vocal fold. J Exp Biol 211: 2144–2154.
- 32. Titze IR (2011) Vocal fold mass is not a useful quantity for describing F0 in vocalization. J Speech Lang Hear Res 54: 520–522.
- 33. Riede T, Brown C (2013) Body size, vocal fold length and fundamental frequency - Implications for mammal vocal communication. Nova Acta Leopoldina NF 111 380: 1–20.
- 34. Newman SR, Butler J, Hammond EH, Gray SD (2000) Preliminary report on hormone receptors in the human vocal fold. J Voice 14: 72–81.
- 35. Voelter C, Kleinsasser D, Joa P, Nowack I, Martínez R, et al. (2008) Detection of hormone receptors in the human vocal fold. Eur Archiv Oto-Rhino-Laryngol 265: 1239–1244.
- 36. Reby D, McComb K, Cargnelutti B, Darwin C, Fitch WT, et al. (2005) Red deer stags use formants as assessment cues during intrasexual agonistic interactions. Proc R Soc Lon B 272: 941–947.
- 37. Taylor AM, Reby D, McComb K (2010) Size communication in domestic dogs (Canis familiaris) growls. Anim Behav 79: 205–210.
- 38. Feinberg DR, Jones CC, Little AC, Burt DM, Perrett DI (2005) Manipulations of foundamental and formant frequencies influence the attractiveness of human male voices. Anim Behav 69: 561–568.
- 39. Apicella CL, Feinberg DR, Marlowe FW (2007) Voice pitch predicts reproductive success in male huntergatherers. Biol Lett 3: 682–684.
- 40. Apicella CL, Feinberg DR (2009) Voice pitch alters mate-choice-relevant perception in hunter-gatherers. Proc R Soc Lon B 276: 1077–1082.
- 41.
Charlton BD (2008) Female mate choice in nonhuman animals. In: Weber EA, Krause LH, editors. Animal Behaviour: new research. New York: Nova Science Publishers Inc. pp. 35–56.
- 42. Marshall J, Marshall E (1976) Gibbons and their territorial songs. Science 193: 235–237.
- 43.
Geissmann T (2000) Gibbon song and human music from an evolutionary perspective. In: Wallin NL, Merker B, Brown S editors. The Origins of Music. Cambridge MA: The MIT Press. pp. 103–123.
- 44. Mitani J (1985) Responses of gibbons (Hylobates muelleri) to self, neighbor, and stranger song duets. Int J Primatol 6: 193–200.
- 45. Raemaekaers PM, Raemaekaers JJ (1985) Field playback of loud calls to gibbons (Hylobates lar): territorial, sex-specific and species-specific responses. Anim Behav 33: 481–493.
- 46.
Preuschoft H, Chivers D, Brockelman WY, Creel N (1984) The Lesser Apes: Evolutionary and Behavioural Biology. Edinburgh: Edinburgh University Press.
- 47. Mitani JC (1987) Territoriality and monogamy among agile gibbons (Hylobates agilis). Behav Ecol Sociobiol 20: 265–269.
- 48. Geissmann T (1999) Duet songs of the siamang, Hylobates syndactylus: II. Testing the pair-bonding hypothesis during a partner exchange. Behaviour 136: 1005–1039.
- 49. Geissmann T, Orgeldinger M (2000) The relationship between duet songs and pair bonds in siamangs, Hylobates syndactylus. Anim Behav 60: 805–809.
- 50.
Lappan S, Whittaker DJ (2009) The Gibbons: New Perspectives on Small Ape Socioecology and Population Biology. New York: Springer-Verlag.
- 51. Padgham M (2004) Reverberation and frequency attenuation in forests-implications for acoustic communication in animals. J Acoust Soc Am 115: 402–410.
- 52. Nemeth E, Dabelsteen T, Pedersen SB, Winkler H (2006) Rainforests as concert halls for birds: Are reverberations improving sound transmission of long song elements? J Acoust Soc Am 119: 620–626.
- 53. Schneider C, Hodges K, Fischer J, Hammerschmidt K (2008) Acoustic niches of Siberut primates. Int J Primatol 29: 601–613.
- 54. Sun G-Z, Huang B, Guan Z-H, Geissmann T, Jiang X-L (2011) Individuality in male song of wild black crested gibbons (Nomascus concolor). Am J Primatol 73: 431–438.
- 55. Barelli C, Heistermann M, Boesch C, Reichard UH (2008) Mating patterns and sexual swellings in pair-living and multimale groups of wild white-handed gibbons (Hylobates lar). Anim Behav 75: 991–1001.
- 56.
Reichard UH (2009) Social organization and mating system of Khao Yai white handed gibbons, 1992–2006. In: Lappan SM, Whittaker D, editors. The Gibbons: New Perspectives on Small Ape Socioecolog and Population Ecology. Berlin: Springer. pp 347–384.
- 57. Barelli C, Heistermann M, Boesch C, Reichard UH (2007) Sexual swellings in wild white-handed gibbon females (Hylobates lar) indicate the probability of ovulation. Horm Behav 51: 221–230.
- 58.
Barelli C, Mastudaira K, Roos C, Wolf T, Heistermann M, et al.. (2013) Extra-pair paternity confirmed in wild white-handed gibbons. Am J Primatol. DOI: 10.1002/ajp.22180.
- 59. Barelli C, Heistermann M (2012) Sociodemographic correlates of fecal androgen levels in wild male white-handed gibbons (Hylobates lar). Int J Primatol 33: 784–798.
- 60.
Lynam AJ, Round PD, Brockelman WY (2006) Status of birds and large mammals in Thailand's Dong Phayayen – Khao Yai Forest Complex. Biodiversity Research and Training (BRT) Program and Wildlife Conservation Society, Bangkok.
- 61. Barelli C, Reichard UH, Mundry R (2011) Is grooming used as a commodity in wild white-handed gibbons, Hylobates lar? Anim Behav 82: 801–809.
- 62. Kappeler PM, van Schaik CP (2002) Evolution of primate social systems. Int J Primatol 23: 707–740.
- 63. Möhle U, Heistermann M, Palme R, Hodges JK (2002) Characterization of urinary and fecal metabolites of testosterone and their measurement for assessing gonadal endocrine function in male nonhuman primates. Gen Comp Endocrinol 129: 135–45.
- 64. Girard-Buttoz C, Heistermann M, Krummel S, Engelhardt A (2009) Seasonal and social influences on fecal androgen and glucocorticoid excretion in wild male long-tailed macaques (Macaca fascicularis). Physiol Behav 98: 168–175.
- 65. Weingrill T, Willems EP, Zmmermann N, Steinmetz H, Heistermann M (2011) Species-specific patterns in fecal glucocorticoid and androgen levels in zoo-living orangutans (Pongo spp.). Gen Comp Endocrinol 172: 446–457.
- 66. Thinh VN, Hallam C, Roos C, Hammerschmidt K (2011) Concordance between vocal and genetic diversity in crested gibbons. BMC Evol Biol 11: 36.
- 67. Koda H, Nishimura T, Tokuda IT, Oyakawa C, Nihonmatsu T, et al. (2012) Soprano singing in gibbons. Am J Phys Anthropol 149: 347–355.
- 68.
McGregor PK (1992) Quantifying responses to playback: one, many, or composite multivariate measures? In: McGregor PK editor. Playback and Studies of Animal Communication. New York, London: Plenum Press.
- 69. Mundry R, Sommer C (2007) Discriminant function analysis with nonindependent data: consequences and an alternative. Anim Behav 74: 965–976.
- 70.
Baayen RH (2008) Analyzing Linguistic Data. Cambridge: Cambridge University Press. 368 p.
- 71. van de Pol M, Wright J (2009) A simple method for distinguishing within- versus between-subject effects using mixed models. Anim Behav 77: 753–758.
- 72. Schielzeth H, Forstmeier W (2009) Conclusions beyond support: overconfident estimates in mixed models. Behav Ecol 20: 416–420.
- 73. Pirovino M, Heistermann M, Zimmermann N, Zingg R, Clauss M, et al. (2011) Fecal glucocorticoid measurements and their relation to rearing, behavior, and environmental factors in the population of pileated gibbons (Hylobates pileatus) held in European zoos. Int J Primatol 32: 1161–1178.
- 74. Bolker BM, Brooks ME, Clark CJ, Geange SW, Poulsen JR, et al. (2008) Generalized linear mixed models: a practical guide for ecology and evolution. Trends Ecol Evol 24: 127–135.
- 75. Sarkar SK, Chang C-K (1997) The Simes method for multiple hypothesis testing with positively dependent test statistics. J Am Stat Assoc 92: 1601–1608.
- 76.
R Core Team (2012) R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria.
- 77.
Bates D, Maechler M, Bolker B (2013) lme4: Linear mixed-effects models using S4 classes. R package version 0.999999-2.
- 78.
Baayen RH (2011) languageR: Data sets and functions with “Analyzing Linguistic Data: A practical introduction to statistics”. R package version 1.4.
- 79.
Venables WN, Ripley BD (2002) Modern Applied Statistics with S. Fourth Edition. New York: Springer.
- 80. Oyakawa C, Koda H, Sugiura H (2007) Acoustic features contributing to the individuality of wild agile gibbon (Hylobates agilis agilis) songs. Am J Primatol 69: 777–790.
- 81. Owren MJ, Seyfarth RM, Cheney DL (1997) The acoustic features of vowel-like grunt calls in chacma baboons (Papio cyncephalus ursinus): Implications for production processes and functions. J Acoust Soc Am 101: 2951–2963.
- 82. Fischer J, Hammerschmidt K, Cheney DL, Seyfarth RM (2001) Acoustic features of female chacma baboon barks. Ethology 107: 33–54.
- 83. Notman H, Rendall D (2005) Contextual variation in chimpanzee pant hoots and its implication for referential communication. Anim Behav 70: 177–190.
- 84. McCowan B, Hooper SL (2002) Individual acoustic variation in Belding's ground squirrel alarm chirps in the High Sierra Nevada. J Acoust Soc Am 111: 1157–1160.
- 85. McComb K, Reby D, Baker L, Moss C, Sayialel S (2003) Long-distance communication of acoustic cues to social identity in African elephants. Anim Behav 65: 317–329.
- 86. Yin S, McCowan B (2004) Barking in domestic dogs: context specificity and individual identification. Anim Behav 68: 343–355.
- 87. Cheney DL, Seyfarth RM (1980) Vocal recognition in free-ranging vervet monkeys. Anim Behav 28: 362–364.
- 88. Hammerschmidt K, Fischer J (1998) The vocal repertoire of Barbary macaques: A quantitative analysis of a graded signal system. Ethology 104: 203–216.
- 89. Harries M, Hawkins S, Hacking J, Hughes I (1998) Changes in the male voice at puberty: vocal fold length and its relationship to the fundamental frequency of the voice. J Laryngol Otol 112: 451–4.
- 90. Harries MLL, Walker JM, Williams DM, Hawkins S, Hughes IA (1997) Changes in the male voice at puberty. Arch Dis Chil 77: 445–447.
- 91. Fitch WT, Giedd J (1999) Morphology and development of the human vocal tract: a study using magnetic resonance imaging. J Acoust Soc Am 106: 1511–1522.
- 92. Jenkins JS (1998) The voice of the castrato. Lancet 351: 1877–1880.
- 93.
Titze IR (2000) Principles of Voice Production. Iowa City: National Center for Voice and Speech.
- 94. Wich SA, van der Post DJ, Heistermann M, Möhle U, van Hooff JARAM, et al. (2003) Life-phase related changes in male loud call characteristics and testosterone levels in wild Thomas langurs. Int J Primatol 24: 1251–1265.
- 95. Koren L, Geffen E (2009) Complex call in male rock hyrax (Procavia capensis): a multi-information distributing channel. Behav Ecol Sociobiol 63: 581–590.
- 96. Chan RW, Fu M, Young L, Tirunagari N (2007) Relative contributions of collagen and elastin to elasticity of the vocal fold under tension. Ann Biomed Eng 35: 1471–1483.
- 97. Riede T (2010) Elasticity and stress relaxation of rhesus monkey (Macaca mulatta) vocal folds. J Exp Biol 213: 2924–2932.
- 98. Morton ES (1977) On the occurrence and significance of motivation structural rules in some bird and mammal sounds. Am Nat 111: 855–869.
- 99. Fichtel C, Hammerschmidt K (2003) Responses of squirrel monkeys to their experimentally modified mobbing calls. J Acoust Soc Am 113: 2927–2932.
- 100. Christie PJ, Mennill DJ, Ratcliffe LM (2004) Pitch shifts and song structure indicate male quality in the dawn chorus of black-capped chickadees. Behav Ecol Sociobiol 55: 341–348.
- 101. Reby D, Charlton BD, Locatelli Y, McComb K (2010) Oestrous red deer hinds prefer male roars with higher fundamental frequencies. Proc R Soc Lon B 277: 2747–2753.
- 102. Ten Cate C, Slabbekoorn H, Ballintijn MR (2002) Birdsong and male–male competition: Causes and consequences of vocal variability in the collared dove (Streptopelia decaocto). Ad Stud Behav 31: 31–75.
- 103. Reby D, McComb K (2003) Anatomical constraints generate honesty: acoustic cues to age and weight in the roars of red deer stags. Anim Behav 65: 519–530.
- 104. Enquist M, Leimar O, Ljungberg T, Mallner Y, Segerdahl N (1990) A test of the sequential assessment game: fighting in the cichlid fish Nannacara anomala. Anim Behav 40: 1–14.
- 105. Maynard Smith J, Harper D (1995) Animal signals: models and terminology. J Theor Biol 177: 305–311.