Age related decline in female lar gibbon great call performance suggests that call features correlate with physical condition
© Terleph et al. 2015
Received: 18 August 2015
Accepted: 21 December 2015
Published: 5 January 2016
White-handed gibbons (Hylobates lar) are small Asian apes known for living in stable territories and producing loud, elaborate vocalizations (songs), often in well-coordinated male/female duets. The female great call, the most conspicuous phrase of the repertoire, has been hypothesized to function in intra-sexual territorial defense. We therefore predicted that characteristics of the great call would correlate with a caller’s physical condition, and thus might honestly reflect resource holding potential (RHP). Because measurement of RHP is virtually impossible for wild animals, we used age as a proxy, hypothesizing that great call climaxes are difficult to produce and maintain over time, and that older adults will therefore perform lower quality great calls than young adults. To test this we analyzed the great call climaxes of 15 wild lar gibbon females at Khao Yai National Park, Thailand and 2 captive females at Leo Conservation Center, Greenwich, CT.
Findings show that call climaxes correlate with female age, as young animals (n = 8, mean age: 12.9 years) produced climaxes with a higher frequency range (delta F0), maximum F0 frequency and duty cycle than old animals (n = 9, mean age: 29.6 years). A permuted discriminant function analysis also correctly classified calls by age group. During long song bouts the maximum F0 frequency of great call climaxes’ also decreased. Additional data support the hypothesis that short high notes, associated with rapid inhalation as an individual catches its breath, reflect increased caller effort. Older females produced more high notes than younger females, but the difference only approached statistical significance, suggesting that calling effort may be similar across different ages. Finally, for the first time in this species, we measured peak intensity of calls in captive females. They were capable of producing climaxes in excess of 100 dB at close range (2.7 m).
Age and within-bout differences in the lar gibbon great call climax suggest that call features correlate with physical condition and thus the call may have evolved as an honest signal in the context of intra-sexual territorial defense and possibly also in male mate choice via sexual selection, although further testing of these hypotheses is necessary.
Territorial vocalizations may honestly indicate an individual’s resource holding potential (RHP) if call qualities are determined by physiological constraints [1, 2]. Characteristics of vocalizations, including both spectral and temporal parameters, can vary with body size, age, rank, hormonal state or fatigue level [3–13], all of which are measures of physical condition and can thus be assumed to also be correlates of RHP. Associations between physical condition and vocal qualities have been studied most often in male animals (exceptions: [14, 15]), a sampling bias that may be due to the fact that males show greater variation in RHP, because they are usually thought to be under greater sexual selection pressure than females. More direct physiological measures of RHP, such as energetic output during contests, including ventilation rate or overall O2 consumption (review: ) would be extremely useful, as they could be directly related to field-based data, such as long term territory tenure or reproductive success.
Because direct measures of RHP are difficult to obtain from wild animals under natural conditions , field researchers may have to rely on indirect measures, like body size, to assess physical condition. Unfortunately, body size can likewise be difficult to measure directly in wild animals, and because of this impediment, in primate studies age has been used as a proxy for body condition [6, 18]. In fact, age may be just as reliable an indicator of RHP as size, because regardless of one’s size, physical decline that influences RHP comes with advanced age. Primate respiratory function, for example, is associated with overall health  and progressively decreases with age [20–22]. An individual’s age should affect not only RHP but also the individual’s ability to produce physiologically demanding vocal signals, because changes in respiratory functioning associated with age include a reduction in the lungs’ vital capacity due to loss of elasticity and distensibility, and a reduction in abdominal muscle mass .
Measurable characteristics of song quality that are potentially influenced by lung capacity include call intensity, call rate, and note lengths . These vocal qualities are predicted to change with age and can therefore serve as indirect measures of RHP. While measuring call intensity from wild animals is difficult, measuring its spectral correlates is not. Two vocal parameters associated with loud calls that may signal RHP are the maximum fundamental frequency (F0) and the degree of change in F0 within a call note (delta F0). This is because F0 is dependent on subglottal pressure [25, 26]; so frequency and intensity are not controlled independently in the loud vocalizations of mammals , and receivers from a variety of taxa have been shown to respond to high intensity sounds as a fitness cue [27–30]. Furthermore, in the red deer (Cervus elaphus), for example, females prefer male calls with a higher F0 , and in some primates high-ranking individuals produce loud calls with a higher F0 than lower ranking animals [10, 32]. High F0 can also signal hormonal condition: male lar gibbons (Hylobates lar) show a positive association between androgen levels and F0 , and in other primates, e.g. male chacma baboons (Papio ursinus), androgen levels strongly correlate with social rank, and decline following a lowering in rank .
In addition to loud, high F0 signals potentially indicating that an animal has large or powerful lungs, in animals with large larynxes, like most primates, reaching a high F0 also requires extreme muscular strength to stiffen the vocal folds. Thus receivers may additionally interpret loud, high pitched calls as a fitness signal associated with larynx strength, which may correspond with overall body strength , although to our knowledge this hypothesis has so far not been tested in primates.
Temporal aspects of loud calls can also be associated with caller condition or state. Examples include young adults producing longer calls than more senior individuals , and high-ranking males generating not only longer notes , but also higher call rates, longer bouts and greater participation in competitive calling bouts . Note duration or rate may also correlate with a signaler’s immediate energy levels, potentially providing a more useful cue to an animal’s current physical condition than static body size cues [4, 36].
Gibbons (family Hylobatidae) are small Asian apes that produce loud, complex vocalizations in song bouts, and females produce the great call, often in the context of well-coordinated duets with a mate [37–39]. Great calls likely evolved for distance communication, as they are audible up to 1 km or further through the forest [40–42], although the intensity of these calls has previously been measured in only one gibbon species, the siamang (Symphalangus syndactylus) . The great call is thought to serve primarily as a spacing mechanism between females in inter-group communication, excluding females from each other’s territories [40, 44–47], and as such has had the potential to evolve into an honest signal of RHP.
The great calls of all gibbon species have a conspicuous climax near the end, characterized by notes produced at the highest speed and/or pitch, and a peak in amplitude . For lar gibbons (Hylobates lar) the notes of the call climax contribute substantially to inter-individual variability and span a wide frequency range . The closely spaced notes of the call climax are also the loudest notes produced by females, suggesting that the climax may be difficult to produce. If this is the case, then the great call climax potentially serves as an honest signal of a caller’s physical condition and thus may have evolved in the context of resource competition between females, the social context in which great calls are often given , or in the context of mate choice by males.
Although age has not been verified experimentally to correlate with RHP in white-handed gibbons, in order to provide support for the hypothesis that the lar gibbon’s great call climax is an honest index of caller condition, we used age as the most appropriate and available proxy for RHP, given that physical strength generally declines with age, including a reduction in respiratory function [20–23]. We predicted that, because F0 correlates with call intensity, older animals would show (i) a lower maximum F0 frequency, (ii) a lower F0 frequency range in the highest pitched climax note (delta F0), and (iii) a smaller ratio of note to inter-note interval durations (duty cycle). If these climax parameters are constrained by a caller’s condition and thus RHP, it should become increasingly more difficult for an aging individual to maximize them, relative to a younger individual. We also predict that, if older animals have more difficulty maintaining robust call climaxes, (iv) this will be revealed by an increase in brief, high-pitched inspirations, hereafter referred to (in italics) as ‘high notes’. Our previous research has shown that high notes occur between the longer, regular expiration notes of lar gibbon great calls, in particular following the climax notes that rise in frequency or that are maintained at high frequency and intensity . They result from brief inhalations as the animal “catches its breath” between longer notes, with an average of two high notes occurring per great call (range: 0–5) . We tested the hypothesis that high notes correspond with an increased vocal effort by analyzing the length and F0 frequency of great call notes that immediately precede high notes, versus great call notes that do not, predicting that (v) high notes will follow longer, higher pitched great call climax notes that span a wider frequency range.
To further test the hypothesis that our measures of the great call climax are potential indices of RHP, we analyzed how vocal parameters change within individual singing bouts, predicting that (vi) as a bout of great calls progresses, the maximum F0, delta F0, and climax duty cycle will decrease, and the number of high notes will increase, possibly because an individual’s vocal ability becomes taxed over the course of singing.
Finally, we also report, for the first time in this species, the peak intensity of great call climaxes, using measurements from captive animals. By quantifying the high intensity of great call climaxes, we lend support to the hypothesis that production of great call climaxes pushes females to their physical limits.
Subjects and recordings
We made recordings from a large group of adult female lar gibbons, spanning a wide age range (9–45.2 years). 15 of the subjects were wild animals, recorded at the Mo Singto—Klong E-Tau long-term research site at Khao Yai National Park, Thailand (101°22’E, 14°26’N), 130 km NE of Bangkok in June-July 2013 and July-August 2014, under a research permit issued by the National Research Council of Thailand (NRCT) and the Department of National Parks, Wildlife and Plant Conservation (DNP) of Thailand. All wild lar gibbons have been habituated to the presence of researchers for decades, in some cases for their entire lifetime. We also obtained recordings from two captive female lar gibbons at the LEO Zoological Conservation Center in Greenwich, Connecticut, USA. Each captive animal was housed with an adult male in one of two large outdoor enclosures. We made all recordings with a digital solid-state recorder (Marantz PMD661, Kanagawa, Japan) and a Sennheiser microphone: K6 power module and ME67 recording head and windscreen (Sennheiser Electronic, Wedemark, Germany). All recordings were 24 bit, with a 48 kHz sample rate.
All research was approved by the Institutional Animal Care Committee of Sacred Heart University, and adhered to the American Society of Primatologists Principles for the Ethical Treatment of Non-Human Primates. We avoided direct interaction with animals in all cases, but to minimize recording distances we usually approached wild gibbons close enough to record from directly underneath the tree in which an individual was vocalizing (<50 m), and in direct line of sight. We recorded all calls between 6:00 and 11:00 AM, as the females engaged in morning duets with males, often in response to the calls of nearby animals. We identified individual wild animals by territory location and individual pelage color and markings based on long-term photo-records and long-term knowledge of individuals .
Ages of all females studied
Age at first vocal recording (years/months)
Window of accuracy (days)
Females with known birth date ( N = 8)
Females with estimated birth date ( N = 9)
Our method of calculating female age based on first parturition potentially underestimates the age of these 4 animals, if individuals had a birth prior to their natal dispersal and/or before their emigration into the study population. However, secondary female dispersal has not yet been observed in the Khao Yai white-handed gibbon population and our approach to calculate female ages was conservative as it resulted in minimum ages. Moreover, a margin of error of only a few years is unlikely to be consequential, particularly for the purpose of grouping animals into the broad categories of young and old, given the long lifespan of the species and the broad range of ages sampled. The maximum lifespan of captive Hylobates has been reported to be 35.6 years , but this is an underestimate, as two of the wild subjects in our study exceeded that age, at 38.2 and 45.2 years, and a captive gibbon (H. muelleri) has been reported to have lived 60 years .
In addition, we combined the measures described above to create an index of call quality for each animal, to obtain a single measure that could be compared across all ages. We created the index by ranking animals, from 1 to 17 (n = 17 animals), in each of the first three measurement parameters (delta F0, maximum F0 and duty cycle). The lowest rank for a given measure, a rank of 1, was thus assigned to the animal that had the smallest value of that measure, and a rank of 17 (or the highest rank) to the animal with the highest value. Consistent with our predictions, high note production was ranked in the reverse order, so that the lowest ranked animal produced the greatest average number of high notes per call. All four ranking positions were then averaged for each animal, to produce its individual ranking index.
Two age groups were compared, using separate between-subject t-tests (two tails), to test the hypothesis that broadly differing age classes are distinguishable, even if more precise ages cannot always be predicted based on call characteristics. Too few subjects were available to divide the sample further. Thus, all animals below the median age of 17.6 years were categorized as young adults (n = 8, mean age: 12.9 years, range: 9–16.7 years), and all others as old adults (n = 9, mean age: 29.6 years, range: 17.6–45.2 years) (Table 1). For each of the measures we first calculated within-animal means, and then compared these means across all 17 animals. We corrected for multiple t-tests using the method described by Benjamini and Hochberg , which adjusted the P level of significance from 0.05 to 0.0375. In order to explore the overall relationships between age and the call parameters across the entire population, we also made separate regressions for each measurement parameter, as a function of caller age.
In order to determine how often great calls can be accurately assigned to either the young adult or old adult age groups, we also used the call parameters in a nested, permuted discriminant function analysis (pDFA) . We conducted the pDFA using a function written in R  by R. Mundry, and based on the function Ida of the R package MASS . Delta F0 was omitted as a classification parameter for this analysis, as it was highly correlated with the maximum F0 frequency measure. Three animals (D, E and G) were also excluded from the analysis because their call sample sizes were too small (2–3 calls). The remaining 14 animals had a mean of 15.9 ± 1.8 great calls per animal (range: 6–27).
In order to test the hypothesis that high notes correspond with an increase in vocal effort, we used separate paired T-tests to compare the delta F0, maximum F0, and the note duration of climax notes, depending on whether or not they were followed by a high note. We predicted that when a female exaggerates a climax note by increasing these parameters, it will require her to produce a rapid inspiration (high note), and thus high notes will follow longer notes that span a wider frequency range and reach a higher pitch.
The results of our comparison of song parameters by age allowed us to predict that these parameters would change over the course of song bouts, within individuals. We were able to record bouts containing ≥ 8 great calls from seven females, all from the wild population (mean: 9.9 calls per bout; range: 8–12, mean bout length 32.8 min, range 12–54 min). Two animals produced two bouts each; for each of these animals, we averaged the note parameters from each of their bouts and used those averages to represent that animal in subsequent comparisons. We calculated the mean delta F0, maximum F0, climax duty cycle and number of high notes per great call from each animal’s bout. Since bouts differed in the number of great calls produced, in order to generate regressions that compared measurement parameters over the course of bouts, we plotted each parameter over the last 7 great calls of a given bout, as a percentage of the average value obtained from the calls that preceded these 7 calls, at the start of the bout (n = 1–5 great calls, depending on the number of calls in a given bout). We also compared parameters obtained from the first four and last four great calls of bouts, predicting that delta F0, maximum F0, and duty cycle would be decreased, and high note production increased from early to late in bouts. We used matched-sample, two-tailed t-tests to compare the early versus late bout parameters.
We recorded the peak intensity (dB SPL, A scale, 125 ms ‘fast’ time weighting) of great call climaxes from the two captive animals (ages 9.0 and 25.9 years) in their outdoor enclosures, each housed with an adult male. We obtained dB measures with an SPL meter (Checkmate CM-150, Galaxy Audio, Wichita, KS) with a built in ½ inch electret condenser microphone and a windscreen (frequency range 31.5Hz-8KHz, measuring level range 30-130 dB, accuracy ±1.5 dB), and only collected data when the animals were oriented towards the microphone. Animals were separated from the microphone by wire mesh, all calls were produced and recorded from heights of over 2 m, and no major sound reflective structures were in the vicinity of the recording site, so the data are essentially free-field. As subjects were free to move within their enclosures, we had to measure from a range of distances, between 2.7 and 6.1 m. In order to control for amplitude differences due to the effect of spherical spreading, the dB of all calls obtained from distances greater than the nearest (2.7 m) call were adjusted to that distance standard by applying the inverse distance law. In order to calculate the mean dB of all peak call notes for a given individual, each dB measure was converted to intensity, these intensity measures were then averaged, and that average was then converted back into an average dB measure. Peak call amplitudes were compared between the two animals with a two-tailed between-subjects t-test.
Age and great call climaxes
The maximum F0 effects described above may be due to a general lowering of voice in older females, rather than solely due to a failure to reach higher pitch in the peak notes. Lowering of voice has been identified in other female primates, as a correlate of increased body size , but it is not known if gibbons increase in size throughout adulthood. In order to investigate the possibility that overall F0 decreased with age, we also measured the lowest F0 frequency of the climax, but found no difference between the age groups. The lowest F0 frequency was 528 ± 7 Hz for the young adult group and 531 ± 10 Hz for the older females (T (15) = 0.20, p = 0.84). Similarly, there was no correlation between this measure and caller age, across the entire population (r = 0.006) (F (1, 15) = 0.001; p = 0.98). Based on this result, it is unlikely that the lower maximum F0 of old females was related to a general lowering of voice with increasing age, but rather to a failure to produce high climaxes.
Great call climaxes from young adult females had a higher duty cycle (mean 0.76 ± 0.01, range 0.73 – 0.81) than those of old females (mean 0.71 ± 0.01, range 0.66 – 0.79) (T (15) = 2.98, p = 0.01), meaning that young females had a longer duration of notes relative to the intervals between them. The overall negative correlation (r = − 0.425) across the entire population, however, barely failed to reach significance (F (1, 15) = 3.32; p = 0.09).
In contrast to the F0 and duty cycle measures, where young females were expected to show higher performance, old females were predicted to produce more high notes. Since high note averages were derived from count data, we confirmed that our high note data were normally distributed before we analyzed them, by performing a Shapiro-Wilk W Test (W = 0.95, p = 0.47). Although a trend for young females to produce fewer high notes than old females existed, it was not significant. The mean rate of high note production by young adult females was 2.17 ± 0.28 and that of the old females was 2.72 ± 0.31 (T (15) = 1.38, p = 0.18). We also did not find a significant correlation between age and the number of high notes produced per great call (r = 0.36) (F (1, 15) = 2.19; p = 0.16).
High note inspirations follow longer climax notes that span a wider frequency
Because high notes are produced by inspiration and more often follow longer expiration climax notes of higher maximum F0 and delta F0, high notes may reflect a greater effort invested in producing the note that precedes them. High notes may be produced involuntarily, as a female rapidly inhales to refill lungs that were depleted by the preceding note. Thus, the production of high notes potentially indicates that a female is having difficulty producing a high quality climax, and that females with waning stamina may be forced to produce more high notes during a song bout.
Climaxes change within song bouts
Maximum F0 decreased over the course of song bouts (r = − 0.78) (F (1, 6) = 9.3; p = 0.023). Mean maximum F0 at the beginning of bouts was 1488 ± 66 Hz, and at the end 1430 ± 67 Hz. (T (6) = 2.79, p = 0.032). This frequency decrease across bouts is consistent with the notion that the production of high frequency call climaxes may be difficult to sustain. A trend of decreasing delta F0 over the course of song bouts barely failed to reach significance, (r = −0.69) (F (1, 6) = 5.5; p = 0.057), and mean delta F0 did not significantly differ between the beginning (640 ± 118 Hz) and the end (566 ± 74 Hz) of bouts: (T (6) = 1.13, p = 0.30). A trend of increasing high note production from the beginning (1.63 ± 0.32) to the end of bouts barely failed to reach statistical significance (1.84 ± 0.3) (T (6) =2.12, p = 0.078), and the overall trend over the course of the entire bouts (r = 0.66) was also not significant (F (1, 6) = 4.7; p = 0.074). We did not find a significant trend in duty cycle change over the course of bouts (r = 0.13) (F (1, 6) = 0.11; p = 0.11). Duty cycle was 0.679 ± 0.04 at the beginning of bouts, and 0.683 ± 0.04 at the end (T (6) = 0.39, p = 0.70). Extended great call bouts thus result in a decrease in maximum F0, but not in robust changes to the other call parameters.
Climaxes can exceed 100 dB
We also measured the intensity of great call climaxes. Each of the two captive animals from whom call intensity was measured was capable of producing climax notes in excess of 100 dB SPL at 2.7 m. The younger animal’s peak call intensity (mean = 105.3 dB, range 101.4-106.9 dB, n = 3) was higher than that of the older female (mean = 98.5 dB, range 80.3-103.5 dB, n = 11): t (12) = 2.2, p = 0.047), demonstrating that maximum call amplitude can differ between individuals. Whether or not such differences are associated with age will require analysis of a significantly larger sample. This finding reveals that the amplitude of lar gibbons’ great calls is comparable to vocalizations that are known for their high intensity from other species, such as the Mantled howling monkey (Alouatta palliata) calls: >90 dB at 5 m , or the roars of lions (Panthera leo): approximately 114 dB at 1 m .
Our study is the first to demonstrate the effect of age on the production of gibbon great calls. The great call is arguably the most significant vocalization in the female gibbon’s vocal repertoire, in part due to its long-range propagation. If the differences in the call parameters that we measured correlate with physical condition, then young females may be capable of signaling high RHP to potential rivals. These hypotheses have not yet been fully tested, but our interpretation is compatible with field observations of female-female replacement, because among the females of known ages that have been observed to date, only young females have replaced old females, and not vice versa (N = 3; ). Here we have analyzed cross-sectional data of vocalizations from multiple individuals, but in order to provide a more reliable estimation of long-term within-individual variation, a longitudinal study involving data collected at different ages from the same individuals would be needed. Cross-sectional studies run the risk of being biased by differential mortality , but that is not likely to be the case for the species we studied, which is very long-lived, with low adult mortality.
Within-bout changes to the climax of the territorial great call and measures of the high intensity of climaxes support the hypothesis that these vocalizations are produced at the limit of callers’ vocal abilities, and thus may indicate a caller’s condition or stamina. We do not yet know if this is the case, and further testing of this hypothesis would be technically challenging. Appropriate tests would require experiments that directly measure individual physical condition and relate it to our great call measures, in addition to field data showing that great calls predict, for example, territory size, quality, length of tenure on a territory, or lifetime reproductive success. Nevertheless, additional characteristics of the great call climax, and the contexts in which it occurs, are in line with the hypothesis. Climaxes have been described as “penetrating and unmusical” , and, unlike the other notes of the species’ repertoire, often contain multiple harmonics and rapid jumps in pitch, unstable and noisy qualities that correspond with the chaotic vocal-fold vibration that typically occurs when a mammal’s vocal folds are taxed by reaching their amplitude and/or frequency limit [64, 65]. Great calls also often fail to reach full pitch on the climax , climax notes are sometimes ‘left out’ while the preceding rising and following falling notes are still produced (pers. observ.), and when calls are occasionally stalled or entirely aborted, it is typically just before the climax , further suggesting that the climax is the most difficult part of the call to produce. It is also noteworthy that, during development, young adolescent and subadult lar females may practice great calls in temporal synchrony with their mother’s calling, yet usually do not produce the climax portion of the call until they are fully mature or, more often, until they sing independently of their mother [66, 67].
The great call differences by age that we have documented could reflect differences in caller motivation rather than physical ability, what game theoretical models describe as ‘aggressiveness signaling’ . For example, older adult females may produce softer, lower-pitched climaxes to avoid conflicts with neighbors, and not because they are incapable of producing more robust calls. Alternatively, younger adults may be more motivated to call because they have more recently acquired a territory and thus may have more to lose than an older individual who has been a long-time resident. Territory ownership is critical for female reproduction, because over the approximately three decades that the population analyzed here has been monitored, no female has successfully reproduced who was not a resident on a territory. Our findings, however, suggest that it may not be the case that the age differences in call climaxes that we have found are related to motivation. High note inhalations likely correlate with caller effort because they follow climax notes, the loudest notes that female lar gibbons produce , and as we have shown here, when an individual produces a high note, it follows a climax note of higher average frequency, higher delta F0 and longer duration. Thus, since both young and old animals produced a similar number of high notes in their great call climaxes, they likely exert a similar effort to produce them. If the effort exerted in call production remains high throughout life, then the differences we have measured in other call parameters are more likely to be attributable to age-related physiological decline than to motivation. Similarly, high notes could potentially maintain a signal’s honesty throughout life and even within bouts, if attempts to exaggerate a climax result in the production of more of these notes, thus revealing to listeners the difficulties an individual of comparatively lower quality may have producing a climax. Playback studies would be an effective method of testing the functional significance of high notes.
If the great call climax contains cues that make it an honest index of a caller’s physical condition, this does not guarantee that other animals attend to such cues, or that they are even capable of detecting the differences that we describe. It is theoretically possible that evolutionary constraints have acted to enforce honest call climaxes, but in the absence of selection specifically for honesty . However, as loud territorial calls that correlate with RHP have been shown to influence conspecific behavior in multiple taxa [8–10, 28–30], and because there is strong competition between gibbon females for possession of a territory , which was likely the ultimate reason for hylobatids living in pairs [69, 70], it is reasonable to expect that an evolutionary pressure existed for female calls to signal RHP. The current study did not test whether listeners can detect and differentially respond to the great call features that we have quantified, for which the most appropriate test would be to measure responses to great calls in playback experiments. Field observations do, however, suggest that other animals attend to the climax. The female’s mate usually ceases calling until just after the great call climax , and neighboring animals also often pause their own duet singing until after the conclusion of another female’s great call climax (pers. observ.). The fact that males stop their calling to listen to neighboring females’ call climaxes further suggests that one possible function of the great call may be to advertise a female’s condition to prospective mates outside of her current pair bond. This is in line with mounting evidence that lar gibbons employ a mixed reproductive strategy that includes females mating polyandrously, both within multi-male groups, with neighboring males from outside the group, and with floaters arriving from more distant locations [33, 52, 53].
Psychophysical data on primate auditory perception suggest that the maximum F0 and delta F0 differences that we quantified here (the smallest being the change in maximum F0 frequency within bouts, at approximately 4 %) are large enough to be detectable by conspecific listeners. In general, auditory thresholds are similar across many primates , but still await measurement for gibbons. Macaques (Macaca fuscata and M. mulatta) can detect frequency changes of approximately 2–3 % for pure-tone stimuli below 2 kHz [73, 74]. Frequency discrimination is similar for baboons , and even greater for humans . Hence, if lar gibbons show frequency sensitivity like that of other primates, the variation in maximum F0 and delta F0 that we report here between age groups and over the course of bouts should be detectable. Similarly, detecting the production of high notes is likely to be trivial for immediate neighbors, as human observers easily hear high notes from distances that are typical of encounters between neighboring animals . But high notes should not be discriminable at greater distances, as they are of high frequency and low amplitude and therefore attenuate rapidly in a forest environment. If high note detection is only possible between immediate neighbors, then those neighbors may be able to discriminate calls more accurately than distant animals, who generally pose more of a threat to one’s territory than an immediate neighbor . However, if high notes are involuntarily produced from rapid inhalation when a female “catches her breath” after producing a robust climax note, then natural selection would likely act to minimize this potential cue of low RHP. More tests are needed to explore the significance of great call high notes, especially because we did not find age differences in their production.
Our measurements of peak call intensity were all obtained from the highest pitched expirations of the call climax, demonstrating that intensity and frequency are correlated in the great call, and likely interdependent. Similarly, a positive correlation between amplitude and F0 occurs in loud vocalizations from multiple taxa [10, 77, 78], and high F0 has been positively associated with social rank as well . As amplitude and F0 are often correlated, receivers could theoretically attend to either one, or both parameters. For lar gibbon great calls, however, maximum F0 is likely to be more reliably transmitted to receivers than intensity. This is because call intensity should vary much more than frequency, both between and within calls, as the caller is often changing her directional orientation. These changes are in fact typical during climaxes, when a female may perform acrobatic locomotor displays that include leaping between tree branches . Additional studies are needed to analyze the relationship between intensity and the spectral and temporal features of great calls. Of the up to 19 currently recognized gibbon species , call intensity has only been measured in one other species, the siamang (Symphalangus syndactylus) . The maximum intensity reported for siamang calls was between 95 to 113 dB SPL, but measurements were taken at varying distances that were not reported for each call, so a direct comparison to the lar female’s peak call intensity, which we report here, is not possible.
Although many studies that analyze the temporal features of calls have measured note duration [10, 13, 81] or note rates [82, 83], we chose a more complex parameter, the duty cycle of call climaxes. We chose this as a potential correlate of caller age because selection does not seem to have favored either long note durations or enough notes to emphasize a rapid delivery rate in the lar gibbon call climax. In regard to note durations, in all gibbon species whose great calls have been analyzed, the climaxes are characterized by multiple notes of short duration, and this appears to be the ancestral condition for hylobatids [48, 79]. Furthermore, lar gibbon climax notes are high-pitched and loud, features associated with a rapid emptying of the lungs and therefore not expected to correlate with a maximization of note duration, especially because the climax is spread over multiple notes. We also did not expect note rate to differ by age or across bouts because, relative to other gibbon species [48, 79], the lar gibbon’s call climax has only a few notes that are produced at a relatively low rate. In contrast, even if selection has not specifically favored long notes or a rapid delivery rate in lar great call climaxes, a high ratio of note to inter-note interval duration (duty cycle) is potentially an index of caller condition because it provides a single metric that summarizes the duration in which multiple loud, high pitched call notes can be sustained, relative to the time needed to refill the lungs between them. It is reasonable to assume that potential indices of caller condition like this, and others , can include integrals of multiple features, especially in the context of a complex, multi-component vocalization such as the great call. Furthermore, if gibbons attend to multiple features of their neighbor’s calls, they may be able to assess age and thus potentially RHP with greater accuracy, as the overall call quality indices that we report here suggest.
Our within-bout results revealed a decrease in maximum F0, similar to the negative age-effect. But unlike the effect of age, no change in duty cycle was observed over the course of bouts. Overall, the within-bout changes that we observed support the notion that the production of great call climaxes with a high F0 cannot be sustained over prolonged bouts. Unfortunately, our within-bout sample contained too few callers to allow for age comparisons across extended bouts.
It is not known what specific proximate mechanisms are responsible for the age-related changes to vocal output that we report. Reproductive hormonal condition is one likely candidate. Although there is an association between androgen level and call qualities in male lar gibbons , the effect of hormones on female great calls has not yet been studied, and is technically challenging as hormone levels fluctuate considerably across the menstrual cycle . Changes in female vocal quality are known to occur throughout life in humans, including during menstruation, pregnancy and menopause [85, 86], and estrogen deprivation associated with menopause changes the mucous membrane lining of the vocal tract which can lead to voice changes . Hormonal state is also correlated with changes to the voice in yellow baboons (Papio cynocephalus cynocephalus) , and Barbary macaques (M. sylvanus) .
Changes to the female voice, especially those mediated by sex steroid hormones, may influence mate choice by males. Barbary macaque males, for example, respond more strongly to female calls that signal an increased probability of conception , and in humans the female voice is rated as more attractive by men as the probability of conception increases . Similar hormonally-influenced vocal cues may be of particular importance in gibbon signaling, as they have the potential to provide information to neighboring males when visual cues, such as sexual swellings, are difficult to assess. Lar gibbons produce only moderate sexual swellings , and inhabit dense forest environments that limit the effectiveness of visual signals. In contrast to visual signals, the great call climax is an effective long-range signal. Although white-handed gibbons are primarily pair-living, long-term demographic data from Khao Yai show that females change a social partner several times over the reproductive life-span . Thus, selection may have promoted a long-range signal in lar gibbon females that could not only signal her RHP, but also might increase her attractiveness to males in the area, thereby increasing her own reproductive success.
Here we have shown that a vocal signal, the great call climax of female lar gibbons, often exceeds 100 dB SPL at close range and shows declines in frequency with age and within bouts, as well as age-related declines in temporal characteristics. If the measures that we used are predictive of an animal’s RHP, which is not yet known, then these findings suggest that the great call might serve as an index of caller condition in the context of female-female competition in territorial defense, and potentially also for attracting a mate(s). Further studies are needed to explore these possibilities.
Thanks to the National Research Council of Thailand for research permits, and the Department of the National Parks, Wildlife, and Plant Conservation for park permits. We also thank K. Homsud and K. Srinoppawan, Superintendents of Khao Yai National Park, Marcella Leone, Director of the LEO Zoological Conservation Center, and R. Mundry, for his assistance with the pDFA test.
Goldberg Research Grant, Nacey Maggioncalda Foundation and a URCG grant, Sacred Heart University to T. Terleph.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Vehrencamp SL. Handicap, index, and conventional signal elements of bird song. In: Espmark Y, Amundsen T, Rosenqvist G, editors. Animal Signals: Signaling and Signal Design in Animal Communication. Trondheim: Tapir Publishers; 2000.Google Scholar
- Maynard Smith J, Harper D. Animal Signals. Oxford: Oxford U Press; 2003.Google Scholar
- Davies NB, Halliday TR. Deep croaks and fighting assessment in toads Bufo bufo. Nature. 1978;274:683–5.View ArticleGoogle Scholar
- Clutton-Brock TH, Albon SD. The roaring of red deer and the evolution of honest advertising. Behaviour. 1979;69:145–70.View ArticleGoogle Scholar
- Ryan MJ, Brenowitz EA. The role of body size, phylogeny, and ambient noise in the evolution of bird song. Am Nat. 1985;126:87–100.View ArticleGoogle Scholar
- Inoue M. Age gradations in vocalization and body weight in Japanese monkeys (Macaca fuscata). Folia Primatol. 1988;51:76–86.PubMedView ArticleGoogle Scholar
- Bee MA, Perrill SA, Owen PC. Size assessment in simulated territorial encounters between male green frogs (Rana clamitans). Behav Ecol Sociobiol. 1999;45:177–84.View ArticleGoogle Scholar
- Sanvito S, Galimberti F, Miller EH. Vocal signalling of male southern elephant seals is honest but imprecise Anim Behav. 2007;73:287–299. Google Scholar
- 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–30.View ArticleGoogle Scholar
- Fischer J, Kitchen DM, Seyfarth RM, Cheney DL. Baboon loud calls advertise male quality: acoustic features and their relation to rank, age, and exhaustion. Behav Ecol Sociobiol. 2004;56:140–8.View ArticleGoogle Scholar
- Ey E, Hammerschmidt K, Seyfarth RM, Fischer J. Age-and sex-related variations in clear calls of Papio ursinus. Int J Primatol. 2007;28:947–60.View ArticleGoogle Scholar
- Vannoni E, McElligott AG. Low frequency groans indicate larger and more dominant fallow deer (Dama dama) males. PLoS One. 2008;3:e3113.PubMedPubMed CentralView ArticleGoogle Scholar
- Erb WM, Hodges JK, Hammerschmidt K. Individual, contextual, and age-related acoustic variation in simakobu (Simias concolor) loud calls. PLoS One. 2013;8(12):e83131.PubMedPubMed CentralView ArticleGoogle Scholar
- Semple S, McComb K. Perception of female reproductive state from vocal cues in a mammal species. P Roy Soc Lond B Bio. 2000;267:707–12.View ArticleGoogle Scholar
- Langmore NE, Davies NB, Hatchwell BJ, Hartley IR. Female song attracts males in the alpine accentor Prunella collaris. P Roy Soc Lond B Bio. 1996;263:141–6.View ArticleGoogle Scholar
- Briffa M, Sneddon LU. Physiological constraints on contest behavior. Funct Ecol. 2007;21:627–37.View ArticleGoogle Scholar
- Ulijaszek SJ. Human energetic methods in biological anthropology. Yearb Phys Anthropol. 1992;35:215–42.View ArticleGoogle Scholar
- Hammerschmidt K, Freundstein T, Jurgens U. Vocal development in squirrel monkeys. Behaviour. 2001;138:1179–204.View ArticleGoogle Scholar
- Wahba WM. Influence of aging on lung function - clinical significance of changes from age twenty. Anesth Analg. 1983;62:764–76.PubMedView ArticleGoogle Scholar
- Awan SN. The aging female voice: acoustic and respiratory data. Clin Linguist Phonet. 2006;20:171–80.View ArticleGoogle Scholar
- Quanjer PH, Stanojevic S, Cole TJ, Baur X, Hall GL, Culver BH, et al. Multi-ethnic reference values for spirometry for the 3–95-yr age range: the global lung function 2012 equations. Eur Respir J. 2012;40:1324–43.PubMedPubMed CentralView ArticleGoogle Scholar
- Herring MJ, Avdalovic MV, Quesenberry CL, Putney LF, Tyler NK, Ventimiglia FF, et al. Accelerated structural decrements in the aging female rhesus macaque lung compared with males. Am J Physiol Lung Cell Mol Physiol. 2013;304:125–34.View ArticleGoogle Scholar
- Sataloff RT, Rosen DC, Hawkshaw M, Spiegel JR. The three ages of voice, the aging adult voice. J Voice. 1997;2:156–60.View ArticleGoogle Scholar
- Fitch WT, Hauser MD. Unpacking “honesty”: vertebrate vocal production and the evolution of acoustic signals. In: Simmons A, Fay RR, Popper AN, editors. Acoustic Communications. New York: Springer; 2002. p. 65–137.Google Scholar
- Titze IR, Riede TA. cervid vocal fold model suggests greater glottal efficiency in calling at high frequencies. PLoS Comput Biol. 2010;6:e1000897.PubMedPubMed CentralView ArticleGoogle Scholar
- Titze I. On the relation between subglottal pressure and fundamental frequency in phonation. J Acoust Soc Am. 1989;85:901–6.PubMedView ArticleGoogle Scholar
- Ritschard M, Riebel K, Brumm H. Female zebra finches prefer high amplitude song. Anim Behav. 2010;79:877–83.View ArticleGoogle Scholar
- Latimer W, Sippel M. Acoustic cues for female choice and male competition in Tettigonia cantans. Anim Behav. 1987;35:887–900.View ArticleGoogle Scholar
- Arak A. Female mate selection in the natterjack toad: active choice or passive attraction? Behav Ecol Sociobiol. 1988;22:317–27.Google Scholar
- Castellano S, Rosso A, Laoretti F, Doglio S, Giacoma C. Call intensity and female preferences in the European green toad. Ethology. 2000;106:1129–41.View ArticleGoogle Scholar
- Reby D, Charlton BD, Locatelli Y, McComb K. Oestrous red deer hinds prefer male roars with higher fundamental frequencies. P Roy Soc Lond B Bio. 2010;277:2747–53.View ArticleGoogle Scholar
- Riede T, Arcadi AC, Owren MJ. Nonlinear acoustics in the pant hoots of common chimpanzees (Pan troglodytes): vocalizing at the edge. J Acoust Soc Am. 2007;121:1758–67.PubMedView ArticleGoogle Scholar
- Barelli C, Mundry R, Heistermann M, Hammerschmidt K. Cues to androgens and quality in male gibbon songs. PLoS One. 2013;8:e82748.PubMedPubMed CentralView ArticleGoogle Scholar
- Gesquiere LR, Learn NH, Simao MC, Onyango PO, Alberts SC, Altmann J. Life at the top: rank and stress in wild male baboons. Science. 2011;333:357–60.PubMedPubMed CentralView ArticleGoogle Scholar
- Kitchen DM, Seyfarth RM, Fischer J, Cheney DL. Loud calls as indicators of dominance in male baboons (Papio cynocephalus ursinus). Behav Ecol Sociobiol. 2003;53:374–84.Google Scholar
- Fitch WT, Hauser MD. Vocal production in nonhuman primates: Acoustics, physiology, and functional constraints on "honest" advertisement. Am J Primatol. 1995;37:191–219.View ArticleGoogle Scholar
- Marshall JT, Marshall ER. Gibbons and their territorial songs. Science. 1976;193:235–7.PubMedView ArticleGoogle Scholar
- Haimoff EH. Acoustic and organizational features of gibbon songs. In: Preuschoft H, Chivers DJ, Brockelman WY, Creel N, editors. The lesser apes: evolutionary and behavioural biology. Edinburgh: Edinburgh University Press; 1984.Google Scholar
- Preuschoft H, Chivers DJ, Brockelman WY, Creel N, editors. The lesser apes: evolutionary and behavioural biology. Edinburgh: Edinburgh University Press; 1984.Google Scholar
- Raemaekers JJ, Raemaekers PM, Haimoff EH. Loud calls of the gibbon (Hylobates lar) - repertoire, organization and context. Behaviour. 1984;91:146–89.View ArticleGoogle Scholar
- Reichard U, Sommer V. Group encounters in wild gibbons (Hylobates lar): agonism, affiliation, and the concept of infanticide. Behaviour. 1997;134:1135–74.View ArticleGoogle Scholar
- Asensio N, Brockelman WY, Malaivijitnond S, Reichard UH. White-handed gibbon (Hylobates lar) core area use over a short-time scale. Biotropica. 2014;46:461–9.View ArticleGoogle Scholar
- Todd NPM, Merker B. Siamang gibbons exceed the saccular threshold: intensity of the song of Hylobates syndactylus. J Acoust Soc Am. 2004;115:3077–80.PubMedView ArticleGoogle Scholar
- Mitani JC. Gibbon song duets and intergroup spacing. Behaviour. 1985;92:59–96.View ArticleGoogle Scholar
- Raemaekers JJ, Raemaekers PM. Field playback of loud calls to gibbons (Hylobates lar): territorial, sex-specific and species-specific responses. Anim Behav. 1985;33:481–93.View ArticleGoogle Scholar
- Cowlishaw G. Song function in gibbons. Behaviour. 1992;121:131–53.View ArticleGoogle Scholar
- Cheyne SM, Chivers DJ, Sugardjito J. Covariation in the great calls of rehabilitant and wild gibbons Hylobates agilis albibarbis. Raffles Bull Zool. 2007;55:201–7.Google Scholar
- Geissmann T. Evolution of communication in gibbons (Hylobatidae). Faculty II: Ph.D: Zürich University Anthropological Institute, Philosoph; 1993.Google Scholar
- Terleph TA, Malaivijitnond S, Reichard UH. Lar gibbon (Hylobates lar) great call reveals individual caller identity. Am J Primatol. 2015;77:811–21.PubMedView ArticleGoogle Scholar
- Reichard UH. Social organization and mating system of Khao Yai white-handed gibbons, 1992–2006. In: Lappan S, Whittacker DJ, editors. The gibbons: new perspectives on small ape socioecology and population biology. NY: Springer; 2009.Google Scholar
- Brockelman WY, Reichard U, Treesucon U, Raemaekers JJ. Dispersal, pair formation and social structure in gibbons (Hylobates lar). Behav Ecol Sociobiol. 1998;42:329–39.View ArticleGoogle Scholar
- Barelli C, Boesch C, Heistermann M, Reichard UH. Female white-handed gibbons (Hylobates lar) lead group movements and have priority of access to food resources. Behaviour. 2008;145:965–81.View ArticleGoogle Scholar
- Reichard UH, Ganpanakngan M, Barelli C. White handed gibbons of Khao Yai: social flexibility, reproductive strategies, and a slow life history. In: Kappeler PM, Watts D, editors. Long-term field studies of primates. Berlin: Springer; 2012.Google Scholar
- Judge DS, Carey JR. Postreproductive life predicted by primate patterns. J Gerontol Biol Sci. 2000;55:201–9.View ArticleGoogle Scholar
- Geissmann T, Geschke K, Blanchard BJ. Longevity in gibbons (Hylobatidae). Gibbon J. 2009;5:81–92.Google Scholar
- Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Series B Stat Methodol. 1995;57:289–300.Google Scholar
- Mundry R, Sommer C. Discriminant function analysis with nonindependent data: consequences and an alternative. Anim Behav. 2007;74:965–76.View ArticleGoogle Scholar
- R Core Team. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2015.Google Scholar
- Venables WN, Ripley BD. Modern Applied Statistics with S. 4th ed. New York: Springer; 2002.View ArticleGoogle Scholar
- Pfefferle D, Fischer J. Sounds and size: identification of acoustic variables that reflect body size in hamadryas baboons. Papio hamadryas Anim Behav. 2006;72:43–51.View ArticleGoogle Scholar
- Whitehead JM. Long-distance Vocalizations and Spacing in Mantled Howling Monkeys, Alouatta palliata. Chapel Hill: PH D: University of North Carolina; 1985.Google Scholar
- Larom D, Garstang M, Payne K, Raspet R, Lindeque M. The influence of surface atmospheric conditions on the range and area reached by animal vocalizations. J Exp Biol. 1997;200:421–31.PubMedGoogle Scholar
- Nussey DH, Coulson T, Festa-Bianchet M, Gaillard JM. The evolutionary ecology of senescence: measuring senescence in wild animal populations, towards a longitudinal approach. Funct Ecol. 2008;22:393–406.View ArticleGoogle Scholar
- Wilden I, Herzel H, Peters G, Tembrock G. Subharmonics, biphonation, and deterministic chaos in mammal vocalization. Bioacoustics. 1998;9:171–96.View ArticleGoogle Scholar
- Fitch WT, Neubauer J, Herzel H. Calls out of chaos: The adaptive significance of nonlinear phenomena in mammalian vocal production. Anim Behav. 2002;63:407–18.View ArticleGoogle Scholar
- Brockelman WY, Schilling D. Inheritance of stereotyped gibbon calls. Nature. 1984;312:634–6.PubMedView ArticleGoogle Scholar
- Koda H, Lemasson A, Oyakawa C, Rizaldi PJ, Pamungkas J, Masataka N. Possible role of mother-daughter vocal interactions on the development of species-specific song in gibbons. PLoS One. 2013;8:e71432.PubMedPubMed CentralView ArticleGoogle Scholar
- Hurd PL. Resource holding potential, subjective resource value, and game theoretical models of aggressiveness signaling. J Theoret Biol. 2006;241:639–48.View ArticleGoogle Scholar
- Reichard UH. Social monogamy in gibbons: the male perspective. In: Reichard UH, Boesch C, editors. Monogamy mating strategies and partnerships in birds, humans and other mammals. Cambridge: Cambridge U Press; 2003. p. 190–213.View ArticleGoogle Scholar
- Lukas D, Clutton-Brock TH. The evolution of social monogamy in mammals. Science. 2013;341:526–30.PubMedView ArticleGoogle Scholar
- Geissmann T. Duet-splitting and the evolution of gibbon songs. Biol Rev. 2002;77:57–76.PubMedView ArticleGoogle Scholar
- Fay RR. Comparative psychoacoustics. Hearing Res. 1988;34:295–306.View ArticleGoogle Scholar
- Sinnott JM, Owren MJ, Petersen MR. Auditory frequency discrimination in primates: species differences (Cercopithecus, Macaca, Homo). J Comp Psychol. 1987;101:126–31.View ArticleGoogle Scholar
- Pfingst BE. Comparison of spectral and nonspectral frequency difference limens for human and nonhuman primates. J Acoust Soc Am. 1993;93:2124–29.PubMedView ArticleGoogle Scholar
- Hienz RD, Jones AD, Weerts EM. The discrimination of baboon grunt calls and human vowel sounds by baboons. J Acoust Soc Am. 2004;116:1692–7.PubMedView ArticleGoogle Scholar
- Temeles EJ. The role of neighbors in territorial systems – when are they dear enemies. Anim Behav. 1994;47:339–50.View ArticleGoogle Scholar
- Beckers GJL, Suthers RA, ten Cate C. Mechanisms of frequency and amplitude modulation in ring dove song. J Exp Biol. 2003;206:1833–43.PubMedView ArticleGoogle Scholar
- Lopez PT, Narins PM, Lewis ER, Moore SW. Acoustically induced call modification in the white-lipped frog. Leptodactylus albilabris Anim Behav. 1988;36:1295–308.View ArticleGoogle Scholar
- Geissmann T. Gibbon Songs and Human Music from an Evolutionary Perspective. In: Wallin NL, Merker B, Brown S, editors. The Origins of Music. Cambridge: MIT Press; 2000. p. 103–23.Google Scholar
- Anandam MV, Groves CP, Molur S, Rawson BM, Richardson MC, Roos C, et al. Species accounts of Hylobatidae. In: Mittermeier RA, Rylands AB, Wilson DE, editors. Handbook of the mammals of the world, 3. Primates. Barcelona: Lynx Edicions; 2013. p. 778–91.Google Scholar
- Wich SA, van der Post DJ, Heistermann M, Mohle U, van Hooff JARAM, Sterck EHM. Life-phase related changes in male loud call characteristics and testosterone levels in wild Thomas langurs. Int J Primatol. 2003;24:1251–65.View ArticleGoogle Scholar
- Boinski S, Campbell AF. Use of trill vocalizations to coordinate troop movement among white-faced capuchins—a 2nd field test. Behaviour. 1995;11–12:875–901.View ArticleGoogle Scholar
- Koda H, Shimooka Y, Sugiura H. Effects of caller activity and habitat visibility on contact call rate of wild Japanese macaques (Macaca fuscata). Am J Primatol. 2008;70:1055–63.PubMedView ArticleGoogle Scholar
- Barelli C, Heistermann M, Boesch C, Reichard UH. Sexual swellings in wild white-handed gibbon females (Hylobates lar) indicate the probability of ovulation. Horm Behav. 2008;51:221–30.View ArticleGoogle Scholar
- Bryant GA, Haselton MG. Vocal cues of ovulation in human females. Biol Letters. 2009;5:12–5.View ArticleGoogle Scholar
- Caruso S, Roccasalva L, Sapienza G, Zappala M, Nuciforo G, Biondi S. Laryngeal cytological aspects in women with surgically induced menopause who were treated with transdermal estrogen replacement therapy. Fertil Steril. 2000;74:1073–9.PubMedView ArticleGoogle Scholar
- Semple S, McComb K, Alberts S, Altmann J. Information content of female copulation calls in yellow baboons. Am J Primatol. 2002;56:43–56.PubMedView ArticleGoogle Scholar
- Pipitone RN, Gallup GG. Women's voice attractiveness varies across the menstrual cycle. Evol Hum Behav. 2008;29(4):268–74.View ArticleGoogle Scholar