Skip to main content

Mixed messages: wild female bonobos show high variability in the timing of ovulation in relation to sexual swelling patterns

Abstract

Background

The evolution of primate sexual swellings and their influence on mating strategies have captivated the interest of biologists for over a century. Across the primate order, variability in the timing of ovulation with respect to females’ sexual swelling patterns differs greatly. Since sexual swellings typically function as signals of female fecundity, the temporal relation between ovulation and sexual swellings can impact the ability of males to pinpoint ovulation and thereby affect male mating strategies. Here, we used endocrine parameters to detect ovulation and examined the temporal relation between the maximum swelling phase (MSP) and ovulation in wild female bonobos (Pan paniscus). Data were collected at the Luikotale field site, Democratic Republic of Congo, spanning 36 months. Observational data from 13 females were used to characterise female swelling cycles (N = 70). Furthermore, we measured urinary oestrone and pregnanediol using liquid chromatography–tandem mass spectrometry, and used pregnanediol to determine the timing of ovulation in 34 cycles (N = 9 females).

Results

We found that the duration of females’ MSP was highly variable, ranging from 1 to 31 days. Timing of ovulation varied considerably in relation to the onset of the MSP, resulting in a very low day-specific probability of ovulation and fecundity across female cycles. Ovulation occurred during the MSP in only 52.9 % of the analysed swelling cycles, and females showed regular sexual swelling patterns in N = 8 swelling cycles where ovulation did not occur. These findings reveal that sexual swellings of bonobos are less reliable indicators of ovulation compared to other species of primates.

Conclusions

Female bonobos show unusual variability in the duration of the MSP and in the timing of ovulation relative to the sexual swelling signal. These data are important for understanding the evolution of sexual signals, how they influence male and female mating strategies, and how decoupling visual signals of fecundity from the periovulatory period may affect intersexual conflict. By prolonging the period during which males would need to mate guard females to ascertain paternity, the temporal variability of this signal may constrain mate-guarding efforts by male bonobos.

Background

In social animals, interactions between males and females are often driven by reproductive interests [14]. In species where males compete for access to females, male behaviour may be influenced by female secondary sexual characteristics (e.g., [57]). Males who are sensitive to female sexual signals often change their behaviour in ways that maximise their mating opportunities and reproductive success [8]. Females may incur high costs from displaying sexual signals if the signals incite aggressive or coercive behaviour from males [9]. Per contra, females may derive benefits from sexual signals if they buffer aggression from males [10], facilitate recruiting agonistic support in conflicts, enhance access to food sources [11], elicit parental care [12], or minimise the risk of infanticide by enabling females to mate polyandrously [13]. Depending on the mating strategies of males, females may exploit sexual signals in order to maximise the benefits that they derive [14]. Females may do so either by displaying prolonged sexual signals which temporally exceed the phase of elevated fecundity [15, 16], or by obscuring cyclic changes in fecundity [17]. Empirical evidence suggests that female sexual signals offer a particularly effective strategy to manipulate the behaviour of males [1820].

One type of sexual signal, which is absent in most taxa but common in nonhuman primates, is female sexual skin swellings (hereafter sexual swellings). In primate species that possess this morphological trait, the skin surrounding the female genitalia changes in size, shape, turgidity, and colour during the follicular phase, and usually culminates in maximum size and turgidity around the periovulatory period [21, 22]. Among the catarrhine primates, some species have particularly conspicuous or exaggerated sexual swellings, e.g., chimpanzees [16, 23], macaques [24], baboons [25], and bonobos [26]. The majority of species that possess exaggerated sexual swellings live in multimale-multifemale groups [27]. This lends support to the theory that these morphological signals play a role in intersexual communication and function to influence mating patterns.

Several hypotheses pertaining to the evolution and function of exaggerated sexual swellings have been reviewed in the primate literature [2830] with the breadth of hypotheses reflecting the variance in the reliability of this trait. The reliable indicator hypothesis [31] proposes that sexual swellings are honest signals of female quality or condition. For example, sexual swellings seem to reflect aspects of long-term reproductive value in female chacma baboons [32]. The graded-signal hypothesis [28] proposes that variation in the size of exaggerated sexual swellings conveys information about the probability of ovulation, allowing males to adjust their mating strategies according to their individual monopolisation potential. Across primate species, the precision with which sexual swellings signal ovulation varies considerably along a continuum. Near one end of the continuum, sexual swellings are relatively precise indicators of ovulation through which females signal their reproductive status and fecundity to conspecifics (e.g., [3336]). Thereby, females attract mating partners and potentially incite male-male competition [36]. Further along the continuum, sexual swellings indicate ovulation with less precision and reliability [37], although the highest probability of ovulation still occurs when swellings are at their peak size during the maximum swelling phase (MSP) [16, 38, 39]. In these species, males often adjust their mate-guarding efforts in accordance with the signal, and females are able to bias or confuse paternity depending on the amount of variability in the timing of ovulation and in the duration of the MSP. The more a female’s MSP duration exceeds the period of fecundity and the greater the variability in the timing of ovulation within the MSP, the less precisely swellings indicate the timing of ovulation. In turn, it becomes increasingly difficult for males to time their mating efforts accurately within the period of high fecundity [30].

Early invasive studies of captive nonhuman primates revealed that cycles of sexual swelling are modulated by hormonal events of the menstrual cycle, with tumescence of the sexual swelling caused by an increase in oestrogens, while progesterone can act synergistically with oestrogens or inhibit their effect and thereby result in swelling detumescence [23, 4043]. However, a female’s social environment also can affect ovarian function and swelling size [4446]. Flexibility in the expression of sexual swellings has been noted in studies where different social cues or events caused females to display situation-dependent swellings [47, 48] or deceptive swellings in the absence of ovulation [49]. Variability in the flexibility and accuracy of this signal therefore merits further investigation.

The same hormones underlying female sexual swellings, i.e., oestrogens and progesterone, are also used in many studies of human and nonhuman primates to determine the occurrence of ovulation and to pinpoint its timing [16, 50]. One species that is notable in terms of variability in timing of ovulation in relation to sexual swelling patterns and female reproductive strategies is the bonobo (Pan paniscus). To date, only a few studies have measured ovarian hormones in female bonobos, and all were conducted exclusively in captivity [5154]. These studies indicated that the variability in the timing of ovulation in relation to sexual swelling patterns in bonobos is considerably higher than in other species, with ovulation sometimes occurring outside the MSP [53]. However, because energy intake can alter hormone levels and reproductive functioning, and is generally higher and more regular in captive settings, field studies are critical for furthering our understanding of the reproductive endocrinology and sexual signalling of wild bonobos.

Bonobos live in multimale-multifemale societies in which females exhibit extended sexual receptivity and mating is polygynous and polyandrous. During the ovarian cycle, female bonobos show conspicuous changes in the size, shape, and firmness of their sexual swellings and the skin surrounding their perinea [55, 56]. Bonobos are a particularly interesting species in which to study signal accuracy and honesty in tandem with reproductive endocrinology, since females’ sexual swellings often display an unusually lengthy MSP both in captivity [57, 58] and in the wild [59]. Furthermore, wild female bonobos have a shorter period of lactational amenorrhoea compared to chimpanzees [6062] and often resume having sexual swelling cycles within 1–2 years following parturition. Finally, bonobos are an interesting species in which to study sexual signalling given that females are often dominant over males [63], males do not engage in sexual coercion, and mating behaviour is used in a variety of contexts [64, 65].

In this study, we investigated the relationship between ovarian hormones and sexual swellings in wild female bonobos. Our specific objectives were to: (1) quantify the duration of sexual swelling cycles of wild female bonobos; (2) assess the intraindividual and interindividual variability in the duration of the MSP; (3) quantify reproductive hormone levels by measuring urinary oestrone and pregnanediol to generate hormone profiles across swelling cycles; (4) examine the temporal relation between sexual swellings and ovulation; and (5) investigate factors which may influence the duration of the MSP. Our final aim was to assess the predictability of ovulation in wild bonobos by examining how accurately sexual swellings signal ovulation and fecundity in female bonobos.

Methods

Study site and subjects

Data were collected from the Bompusa community of wild bonobos at the Luikotale field site, located near the southern sector of Salonga National Park, Democratic Republic of the Congo [66]. Data collection spanned a three-year period, from December 2010 to December 2013. During this study, the Bompusa community consisted of 13–14 mature females (estimated age > 10 yrs; Table 1), 6–7 mature males (estimated age > 10 yrs), and 11–19 dependent offspring. We based age estimates on physical features such as body size, dentition, physical signs of aging, and sexual swellings. Females were considered to be nulliparous if they immigrated to the community as young females, were not observed to have given birth, and if the appearance of their nipples suggested that they had not lactated previously. All bonobos were fully habituated and individually recognisable, and all females were long-term, permanent residents, with the exception of one female (Djulie) who immigrated to the community in 2012. Eight females gave birth during the study period.

Table 1 Identity and characteristics of mature females in the Bompusa community

Assessment of sexual swellings and the maximum swelling phase (MSP)

Following a previously established methodology for bonobos [57, 67], sexual swellings of all mature females observed each day were visually inspected and assigned a score (1 to 4) based on relative tumescence, firmness, labial occlusion, and lustre. Swelling stage 4 was termed the maximum swelling phase (MSP), and was characterised by fully taut and swollen tissue of the sexual swelling, a sturdy appearance during locomotion, labial occlusion, and lustre [67].

The duration of the MSP was calculated from the first to the last day of swelling stage 4. In several cycles, female swelling scores temporarily dropped to stage 3 during the MSP. If the swelling returned to stage 4 following the temporary drop to stage 3 and did not detumesce beyond stage 3, then these temporary periods (range: 1–4 days) of swelling stage 3 were included in the overall duration of the MSP. Only cycles that did not have a sample gap greater than one day at the onset or end of the MSP were used to analyse variation in the duration of the MSP. A total of 70 swelling cycles from 13 females fulfilled this criterion (Table 1). When there was a one-day sample gap, we used a more conservative estimate of the MSP and only counted the number of days that we observed a female with swelling stage 4.

Urine sample collection

Urine sample collection targeted a subset of nine focal females (see Table 1), representing females of different ages and reproductive states, e.g., cycling, pregnant, and lactating. During swelling stages 1–3, samples were collected from focal females at regular intervals, e.g., every two to three days, to monitor hormone levels and to enable detection of any ovulations outside the MSP. We attempted to collect urine samples every day at the onset of stage 4 (the MSP) and for at least five days after the onset of detumescence. Urine was collected from females using a previously described technique [68]. In brief, the underside of large leaves was used to capture urine from arboreal females. If this was not possible, urine was collected by aspiration from ground vegetation using disposable plastic pipettes. Only urine not contaminated with faeces was collected. Urine was pipetted into 2-ml polypropylene tubes and stored in liquid nitrogen within 12 h of collection. Samples were transferred from liquid nitrogen to dry ice for transport to the Max Planck Institute for Evolutionary Anthropology (MPI-EVA) in Leipzig, Germany, where they were stored at –20 °C until analysis.

Hormone extraction and measurement

We selected a subset of 710 urine samples for analysis in the Endocrine Laboratory of the MPI-EVA. Samples were selected with the objective to equally represent each of the nine focal females with several cycles collected from different points in time. We followed the methods described by Hauser et al. [69] to extract oestrone (E1)–an oestrogen metabolite–and pregnanediol (Pd)–a urinary metabolite of progesterone–amongst other steroid hormones. In brief, 100 μl of each urine sample was spiked pre-extraction with a 50 μl solution containing deuterated internal standards. Progesterone-d9 (Dr. Ehrenstorfer, Augsburg, Germany) and Estrone-d4 (Sigma Chemical Co., St. Louis, MO, USA) were chosen to quantify and correct for matrix loss in Pd and E1 respectively. To de-conjugate the steroid glucuronides, we performed an enzymatic hydrolysis using β-glucuronidase from E. coli (Sigma Chemical Co., St. Louis, MO, USA). This was followed by solvolysis to cleave the steroid sulphate conjugates.

We measured levels of urinary E1 and Pd following the liquid chromatography–tandem mass spectrometry (LC-MS/MS) method [69] to generate profiles of hormone levels across female swelling cycles. This was done using a high performance Waters Alliance 2695 (Waters, Milford, MA) for chromatographic separation. Subsequently, mass spectrometry was performed on a Quattro Premier XE (Micromass, Manchester, UK) with a Z-spray ESI interface. For more details see Additional file 1. We measured creatinine in each sample [70, 71] and indexed hormone values with creatinine to control for variations in the volume and concentration of the voided urine. All hormone values are expressed in ng/mg creatinine (Cr).

Interpretation of hormone profiles and definitions

The timing of ovulation and corresponding onset of the luteal phase was calculated from urinary Pd excretion profiles. A sustained rise in Pd levels above a defined threshold value of two standard deviations above the mean of the preceding three to five baseline values [72] was used to non-invasively detect ovulation, in adherence with other studies [16, 17, 73]. Based on the finding that the Pd rise in urine occurs one to three days after the serum luteinising hormone peak [50, 74] and to account for the time lag in the excretion of hormones into urine, the day preceding the Pd rise was defined as the day of ovulation and designated as day zero. However, there may be an error of one day in this determination.

The fecund phase, or periovulatory period, was defined as the day of ovulation plus the three preceding days [16]. Since sperm remain viable in the female tract for several days, copulations that occur two to three days prior to ovulation can result in conception. This window of fecundity, when conception is most likely to occur, was delineated based on findings from studies of captive nonhuman primates [75] and humans [76]. Cycles were termed anovulatory if a sustained rise in Pd was not detected in the urine. The interovulatory interval (IOI) was calculated from the day of ovulation to the day of ovulation in the next cycle (N = 8 IOI). The interswelling interval (ISI) was used to infer duration of swelling cycles, and was calculated from the onset of the MSP to the onset of the MSP in the subsequent cycle (N = 37 ISI). We used the IOI and the ISI as proxies for cycle duration. We did not use the intermenstrual interval, because menstruation could not be detected reliably in all females. The duration of the luteal phase was calculated from the day of ovulation to the first day of observed menstruation, in cycles where menstruation was detected. Hormone profiles were plotted against female swelling scores to examine the temporal relation between ovulation and swelling cycles. In cases where ovulation occurred outside of the MSP, we associated ovulation with the MSP that was closest in terms of the number of days.

Statistical analysis

We estimated the day-specific probability of ovulation by dividing the number of times we observed ovulation on a particular cycle day by the total number of cycles examined. This was calculated in accordance with Deschner et al. [16] using the equation:

$$ P\left(T = t\right)=\frac{n_t}{n},\kern0.75em t = 1,\ 2,\ 3\dots, $$

where t represents a specific cycle day (relative to the start of the MSP), n t is the number of cycles in which ovulation occurred on day t, and n is the total number of cycles.

Likewise, the day-specific probability of fecundity was estimated following Deschner et al. [16] using the equation:

$$ P\left(X(f)=1\right)={\displaystyle \sum_{t=f}^{f+3}}P\left(\mathrm{T}=\mathrm{t}\right), $$

where (X(f) = 1) represents a day on which a female could conceive, and P(T = t) is as stated above. The day-specific probability of fecundity, also referred to as the probability of conception [77, 78], is a measure of the probability that copulation could lead to conception on any given day.

Models and test predictors

We ran six analyses using linear mixed models (LMMs) and Generalised Linear Mixed Models (GLMMs) [79, 80]. All models were fitted in R version 3.2.4 [81] using the functions lmer or glmer of the package lme4 [82]. We assessed collinearity amongst predictors by deriving Variance Inflation Factors (VIFs) [83, 84], using the function “vif” of the package “car” [85] based on standard linear models lacking the random effects. For each model, we first assessed the significance of the fixed effects as a whole [86], by comparing the fit of the full model to a null model using a likelihood ratio test [87]. The null models lacked the fixed effects. We then determined the significance of the individual fixed effects using likelihood ratio tests [88], comparing the full model with reduced models, dropping the fixed effects one at a time. For each model, we obtained model stability by comparing estimates obtained from the full model with estimates from models with the levels of the random effects excluded one at a time. Since the estimates did not vary greatly [89], all model results were robust.

Female dominance rank and social status can influence ovarian hormone levels [90], the duration of the swelling phase [91], and the duration of cycles and interbirth intervals [91, 92]. Therefore, we included female rank as a fixed effect in all models. Social dominance was assessed and ranks were generated (see Table 1) using the ADAGIO method, version 1.1 [93]. Dominance ranks ranged from one (highest rank) to nine (lowest rank). Female ranks were z-transformed to a mean of zero and a standard deviation of one, prior to fitting each model [94].

MSP duration model

Previous studies of nonhuman primates have proposed that female parity and reproductive state may influence the duration of the MSP (e.g., [53, 95, 96]). Based on these findings, we fitted a LMM to investigate to what extent these factors influenced the duration of female bonobos’ MSPs at Luikotale. As fixed effects, we included female parity as a factor with two levels (“multiparous” and “primiparous”), female reproductive state as a factor with two levels (“cycling”, i.e., experiencing ovulatory cycles, and “not cycling”, e.g., pregnant), number of days since parturition, and female dominance rank as a quantitative predictor. Because the number of days since parturition was skewed and we wanted to avoid outliers that would bias the results, we square root transformed this variable. To control for potential seasonal variation we also included the sine and cosine of the Julian date (after multiplying it by 2 × π and then dividing by 365.25, to convert date into a circular variable). Such a representation of season allowed us to model the response showing a sinusoidal periodicity with a period duration of one year; that is, the response peaking once per year (for more details see [97]). As a random effect, we included female identity (ID). To keep type one error rate at the nominal level of 0.05, random slopes [88, 94] of days since parturition as well as sine and cosine of date within female ID were included in the model. Random slopes of the other fixed effects could not be included, because they varied either rarely within females (e.g., reproductive state) or not at all (female parity).

The sample size for this model was 53 MSPs from 11 females. Since MSP duration was rather skewed, we square root transformed it before fitting the model. This resulted in residuals fulfilling the assumptions of normality and homogeneity (verified by visual inspection of a QQ-plot and residuals plotted against fitted values). Collinearity, assessed by VIFs, appeared to be a minor issue between parity and female rank (maximum VIF: 3.5). Therefore, we fitted two additional LMMs: one excluding the test predictor parity, and a second excluding female rank. These models were fitted and checked in the same way as the main model. Collinearity was not an issue in these additional models (maximum VIF: 1.2).

We tested for absence of influential cases by excluding females one at a time from the data and comparing the estimates derived with those obtained for the full data set, which revealed the model to be stable. To test the overall effect of the fixed effects [86], we compared the full model with a null model that comprised only the effects of season and the random effects, using a likelihood ratio test [87]. Furthermore, to test for significant interindividual variation above and beyond the four fixed effects, we compared the full model to a reduced model lacking only the random intercept term of female ID. The sample size for this reduced model was the same as the full model.

ISI duration model

We fitted a GLMM with poisson error distribution and log link function to investigate variation in the ISI duration. The sample size for this model was 37 ISIs from 13 females. As fixed effects, we included female parity as a factor with three levels (“multiparous”, “nulliparous”, and “primiparous”), female reproductive state as a factor with two levels (“cycling” and “early lactation”), and female rank as a quantitative predictor. Following O’Malley et al. [98], we defined early lactation as 0–24 months following parturition, based on evidence that lactation, and the energetic burden of lactation, are most intense in chimpanzees during the first 24 months following parturition [99, 100]. We chose not to include days since parturition as a fixed effect so that the nulliparous females could be included in this model. As a random effect, we included female ID. We could not control for seasonal variation or days since parturition in this model because including them resulted in model stability problems and the model being too complex given the sample size [84]. As an overall test of the effect of the three fixed effects, we compared the full model with a null model that lacked the fixed effects and comprised only the random effect of female ID. Collinearity was not an issue (maximum VIF: 2.1). The model was overdispersed (dispersion parameter = 1.74) which makes the model anticonservative. However, given the model results, this did not represent a problem (see Results section).

Female rank and occurrence of ovulation model

We fitted a GLMM with binomial error structure and a logit link function to investigate the occurrence of ovulation. Specifically, we tested the influence of female rank and reproductive state on whether or not a cycle was ovulatory (yes/no). The sample size for this model was 34 cycles from nine females. As fixed effects, we included female reproductive state as a factor with two levels (“cycling” and “early lactation”) and female rank as a quantitative predictor. As a random effect, we included female ID. Collinearity was not an issue (maximum VIF: 1.0). To determine the significance of the fixed effects, we compared the full model with a null model that lacked the fixed effects and comprised only the random effect of female ID.

Female rank and timing of ovulation model

We fitted a GLMM with binomial error structure and a logit link function to investigate the timing of ovulation. Specifically, we tested the influence of female rank and reproductive state on whether or not ovulation occurred during the MSP (yes/no). The sample size for this model was 26 cycles from nine females. As fixed effects, we included female reproductive state as a factor with two levels (“cycling” and “early lactation”) and female rank as a quantitative predictor. As a random effect, we included female ID. Collinearity was not an issue (maximum VIF: 1.0). To determine the significance of the fixed effects, we compared the full model with a null model which comprised only the random effect of female ID.

Results

Female sexual swelling cycles, interswelling interval (ISI), and interovulatory interval (IOI)

All females had cyclical fluctuations in the relative degree of tumescence, firmness, lustre, and labial occlusion of their sexual swellings. Although there was variation in female swelling characteristics, e.g., the absolute size of sexual swellings, the four different swelling stages could be reliably distinguished within each female. The mean cycle duration, inferred from the ISI, was (\( \overline{X} \) ± SD) 41.2 ± 13.8 days (range: 20–74 days; N = 37 cycles; N = 13 females). Mean cycle duration calculated from the IOI was very similar, (\( \overline{X} \) ± SD) 40.8 ± 6.8 days (range: 30–51 days; N = 8 cycles; N = 6 females). Mean duration of the luteal phase was (\( \overline{X} \) ± SD) 9.5 ± 1.2 days (range: 8–11 days; N = 6 cycles; N = 4 females). Overall, there was no significant effect of female parity, rank, and reproductive state on the ISI duration (full-null model comparison: χ 2 = 0.90, df = 4, p = 0.925; see Additional file 2: Table S1).

In four out of six pregnancies, females had sexual swelling cycles including a MSP during at least the first three months of pregnancy (range: 1–6 months; N = 6 females). Females resumed showing sexual swelling cycles as early as three months following parturition (\( \overline{X} \) ± SD = 7.0 ± 4.8 months; range: 3–14 months post-parturition; N = 5 females).

Duration of the maximum swelling phase (MSP)

Duration of the MSP ranged from 1–31 days (\( \overline{X} \) ± SD = 10.6 ± 6.8 days; N = 70 cycles; N = 13 females; Fig. 1). MSP duration was highly variable both within and between females (Fig. 2). The full-null model comparison did not reveal a significant effect of female parity, number of days since parturition, female rank, and reproductive state on the MSP duration (χ 2 = 7.01, df = 4, p = 0.136; Additional file 2: Table S2). The model excluding parity was also nonsignificant (full-null model comparison: χ 2 = 5.66, df = 3, p = 0.129; Additional file 2: Table S3). Only the model which excluded female rank had a marginally nonsignificant effect on the MSP duration (full-null model comparison: χ 2 = 6.36, df = 3, p = 0.095; Additional file 2: Table S4). More specifically, we found a tendency for the MSP duration to increase as the number of days since parturition increased. Although this predictor had the strongest effect on the duration of the MPS, it was not statistically significant (see Additional file 2: Table S4). Furthermore, MSP duration did not vary significantly between individuals (χ 2 = < 0.001, df = 1, p = 1.00). There was little variation between mean values of MSP duration from different reproductive states (Table 2). The mean MSP duration was shorter for females during early lactation, i.e., within 24 months of parturition, than during other reproductive states. This corroborates the finding of a tendency for MSP duration to increase as days since parturition increases, as suggested by the results of the MSP model that excluded female rank.

Fig. 1
figure 1

Frequency distribution of the duration of females’ MSPs (N = 70 cycles; 13 females)

Fig. 2
figure 2

Intra-individual variation in MSP duration. Shown are medians (thick horizontal lines), quartiles (boxes), percentiles (2.5 and 97.5 %, vertical lines), as well as the minimum and maximum of each female’s MSP (circles). Females are arranged from highest-ranking (left) to lowest-ranking (right)

Table 2 Mean duration of the maximum swelling phase (MSP) for each reproductive state

Hormone profiles

Representative profiles of urinary E1 and Pd in relation to the pattern of sexual swelling during a nonconception cycle in an individual female are depicted in Fig. 3a. All presumably ovulatory profiles showed a well-defined pattern of Pd excretion, characterised by consistently low levels during the follicular phase followed by markedly elevated levels during the luteal phase. Although profiles of E1 were more variable between and within cycles, levels of E1 gradually increased and reached a discernible peak one to three days before the defined postovulatory Pd rise in 96.2 % of ovulatory cycles.

Fig. 3
figure 3

Profiles of urinary oestrone (open circles) and pregnanediol (closed circles) in relation to changes in sexual swelling size (open bars) in two individual cycles of two female bonobos: (a) a presumably ovulatory cycle from Susi; and (b) an anovulatory swelling cycle from Zoe. The day relative to the first day of the MSP (day 0) is indicated on the x-axis

In contrast to the presumably ovulatory cycles, in eight swelling cycles (collected from N = 5 non-pregnant females) there were only small fluctuations in levels of E1, and increases in E1 were not followed by a significant or sustained rise in Pd (Fig. 3b). The corresponding profiles of urinary Pd also showed a very different pattern in these eight cycles compared to the ovulatory cycles, with no significant rise that would be indicative of ovulation and the start of the luteal phase. In some of these cycles, the rise in E1 may be suggestive of follicular development comparable to the follicular phase of an ovulatory cycle. However, given the consistently low levels of Pd in these cycles, it is highly unlikely that ovulation and a subsequent luteal phase occurred. These eight swelling cycles were therefore considered to be anovulatory. Anovulatory cycles occurred at various points in time throughout a female’s interbirth interval, ranging from 10.8 to 42.3 months after parturition (see Additional file 3: Table S8). The full-null model comparison did not reveal a significant effect of female rank and reproductive state on whether a cycle was ovulatory or anovulatory (χ 2 = 0.45, df = 2, p = 0.799; Additional file 2: Table S5).

Composite profiles

We created composite profiles of urinary E1 and Pd, representing the 26 ovulatory cycles of the nine focal females (Fig. 4). Concentrations of E1 and Pd varied between females and between cycles of individual females, with the largest variation occurring during periods when hormones were at peak values. Mean E1 levels fluctuated between 9.6 ± 2.0 and 17.8 ± 15.4 ng/mg Cr during days –14 to –9 (early follicular phase), then increased to a peak mean level of 63.5 ± 40.1 ng/mg Cr two days before the Pd rise. Mean levels of Pd were consistently low on days during the early follicular phase (range of \( \overline{X} \) ± SD: 83.7 ± 38.9 to 182.8 ± 95.3 ng/mg Cr), then increased rapidly to a mean peak concentration of 3303.3 ± 2597.6 ng/mg Cr during the luteal phase.

Fig. 4
figure 4

Composite profiles of urinary E1 (open circles) and Pd (closed circles), derived from 26 ovulatory cycles from 9 females. Values are presented as mean ± SD and aligned to the day of ovulation (day 0), one day prior to the significant postovulatory rise in Pd, indicating the onset of the luteal phase

Timing of ovulation relative to the maximum swelling phase (MSP)

The timing of ovulation relative to the onset and end of the MSP varied considerably, both between females and within females (N = 26 ovulatory cycles; see Fig. 5). In 18 cycles (69.2 % of ovulatory cycles), ovulation occurred during the MSP, between 2 and 18 days after the onset of the MSP (which corresponded to 1 to 9 days before the end of the MSP). In the remaining eight ovulatory cycles (30.8 %), ovulation occurred outside of the MSP, either during swelling stage 3 before the onset of the MSP (N = 2) or during swelling stage 3 after the onset of detumescence (N = 6). Overall, there was no significant effect of female rank and reproductive state on whether or not ovulation occurred during the MSP (full-null model comparison: χ 2 = 4.46, df = 2, p = 0.107; Additional file 2: Table S6). The MSP duration in each of the eight anovulatory swelling cycles (Fig. 6) was within the range of the MSP duration of the ovulatory cycles.

Fig. 5
figure 5

Duration and timing of the MSP relative to the day of ovulation (day 0) in ovulatory cycles. Asterisks indicate conception cycles; the shaded area indicates the fecund phase. Thick black bars indicate the MSP; thin black bars indicate swelling stage 3; white bars indicate sampling gaps. Female names and cycle numbers are shown on the y-axis. (1)Indicates that this cycle occurred during early lactation. All other cycles occurred > 24 months postparturition

Fig. 6
figure 6

Duration of the MSP in anovulatory swelling cycles. Thick black bars indicate the MSP; thin black bars indicate swelling stage 3; white bars indicate sampling gaps. Female names and cycle numbers are shown on the y-axis. (1)Indicates that this swelling cycle occurred during early lactation. All other swelling cycles occurred > 24 months postparturition

Day-specific probability of ovulation and fecundity

Due to the high variability in the timing of ovulation relative to the MSP, the day-specific probability of ovulation was very low on any given day of a female’s cycle (Fig. 7a). Days three and seven of the MSP revealed the highest probability of ovulation (0.08). The probability of fecundity (i.e., of copulation potentially resulting in conception) was highest on days three and four of the MSP (0.24), and steadily decreased to probabilities below 0.10 by day fourteen. The high variability in the timing of ovulation combined with the fact that some females had long MSPs resulted in very low day-specific probabilities of fecundity. When all cycles from our study were pooled together, high fecundity was not concentrated around specific days of the MSP. This resulted in a high proportion of cycle days when the probability of fecundity was low, but above zero, for the females.

Fig. 7
figure 7

Probability of ovulation (black lines, black circles) and fecundity (red lines, white circles) relative to the first day of the MSP of: (a) bonobos (N = 34 cycles); and (b) chimpanzees at Taï National Park (N = 33 cycles) (calculations based on data from [16])

Discussion

This study provides detailed information about reproductive hormone levels–urinary E1 and Pd–the timing of ovulation, and associated swelling patterns in wild female bonobos. Assessment of sexual swelling patterns revealed that the duration of females’ MSPs was highly variable. Furthermore, there was high variability in the timing of ovulation relative to the onset of the MSP, with ovulation occurring outside of the MSP in almost one third of the ovulatory cycles. Since ovulation did not occur in an additional eight swelling cycles, ovulation occurred during the MSP in only 52.9 % of the analysed cycles. These findings resulted in a very low probability of ovulation occurring on any particular day of a female’s MSP, suggesting that the sexual swellings of wild bonobos are a poor indicator of the timing of ovulation within cycles, and do not always signal fecundity across cycles.

Female dominance rank was not found to have a significant effect on the ISI duration, the MSP duration, the occurrence of ovulation, or the timing of ovulation. Some studies suggested that lower-ranking females may experience suppressed reproductive functioning due to harassment from higher-ranking individuals [92]. If social subordination inhibits the oestrogen-induced surge in luteinising hormone, this can result in anovulatory cycles or premature termination of menstrual cycles [101]. However, aggressive interactions between females are relatively infrequent in bonobos [63, 102]. Dominant females show tolerance toward lower-ranking females, and conflicts are often resolved in non-agonistic ways without overt aggression. Additionally, since several females shared the same dominance rank in our study (Table 1), there was not a steep dominance hierarchy among the females. Access to food resources and the nutritional status of female primates can influence their reproductive hormone levels [103, 104], fecundity (e.g., [105]), and the size and duration of sexual swellings [30]. However, differences in energetic condition and access to resources may be mitigated by the tolerant nature of bonobos [106], and by reduced levels of feeding competition compared to other species [107] (but see [108]). The relationship between feeding ecology and female rank in bonobos remains an area for future investigation.

Duration of the interswelling interval (ISI) and the interovulatory interval (IOI)

The mean duration of ovarian cycles for females in this community, determined from the ISI and the IOI, is within the range previously reported for the intermenstrual interval of bonobos in captivity (\( \overline{X} \) = 33.8–49 days) [52, 57, 58, 109111] and in the wild (\( \overline{X} \) = 42 days) [59]. Differences in cycle duration are likely due to variability in follicular phase duration, as has been found in other studies [52, 112], and as can be inferred from the relatively constant duration of the luteal phase found in our study. Several factors including lactation [52, 113], dominance rank [91], and age [112, 114] have been found to influence duration of the follicular phase. More specifically, mean cycle duration was significantly longer for low-ranking gelada baboons [92] and for lactating captive bonobos [53]. Although there was considerable interindividual and intraindividual variation in our study, dominance rank and reproductive state did not have a significant effect on the ISI duration. The observed variation in the ISI duration merits further investigation.

There was less variability in the IOI compared to the ISI. This suggests that ovulation may be relatively stable temporally, while the onset of the MSP, and swelling patterns in general, vary according to external influences (e.g., social or ecological factors) or due to physiological factors. For example, the consumption of plant steroids can influence endogenous hormone levels [115] and resulted in the suppression of sexual swellings in wild baboons [116]. In other studies, elevated cortisol levels were linked to extended cycle lengths [117] and reproductive suppression [118, 119]. If the IOI remains constant across multiple cycles, variability in the timing of ovulation relative to the MSP might be caused by variability in the onset and duration of the MSP. Accurate assessment of this phenomenon would require the collection of multiple consecutive cycles from the same female.

High variability in the duration of the maximum swelling phase (MSP)

The mean MSP duration (\( \overline{X} \) ± SD = 10.6 ± 6.8 days) was shorter in the Luikotale community compared to wild bonobos at Wamba (\( \overline{X} \) = 14.6 days [59]) and bonobos in captivity (\( \overline{X} \) = 16.0 ± 6.8 days [53]; \( \overline{X} \) = 11.5 days [52]). This could be due in part to our conservative estimate of the MSP duration when there was a sample gap of one day at the onset or end of the MSP. Additionally, studies show that captive primates who are well-nourished and have improved nutritional condition, compared to their wild counterparts, have enhanced reproductive and hormonal functioning [104, 120]. A shorter mean duration of the MSP could be caused by fluctuations in food availability and energetic costs. These fluctuations are likely higher for bonobos in the wild, compared to bonobos in captivity or at research sites during periods when they were provisioned with food [59].

Although there was considerable variation in the MSP duration of the Luikotale bonobos (range: 1–31 days), this finding corroborates other studies of captive bonobos, MSP range: 4–26 days [52] and 3–30 days [53], as well as a study of wild bonobos which reported a MSP range of 3–22 days [59]. The variability in MSP duration for bonobos exceeds the variability found in populations of wild chimpanzees at Taï National Park, Côte d’Ivoire (\( \overline{X} \) = 10.9 days, range: 6–18 days, [16]), at Mahale (\( \overline{X} \) = 11.3 days, range: 5.0–14.0 days, [121]), and at Gombe (\( \overline{X} \) = 9.6 days, range: 7–17 days, [122]).

Consistent with Reichert et al. [53], neither female parity nor reproductive state had a significant effect on the MSP duration. This finding differs from several other species of nonhuman primates, including some populations of chimpanzees, where nulliparous females often have longer MSPs compared to parous females [95, 123]. At Luikotale, females’ MSPs generally increased in duration as time since parturition increased; however, the number of days since parturition did not have a statistically significant effect on MSP duration. Shorter MSPs following parturition could be caused by the energetic demands of lactation [124] and an associated decrease in ovarian steroid production [99]. However, some of our females had swelling cycles with prolonged MSPs while lactating. Contrary to similar studies on other species of nonhuman primates [16, 38], we found that MSP duration did not vary significantly between individuals. This may be due to the large intraindividual variability in MSP duration across different cycles of some females. High intraindividual and interindividual variability in MSP duration contribute to bonobo sexual swellings being less precise signals of ovulation, and render it difficult for male bonobos to predict the timing of ovulation.

Anovulatory swelling cycles in wild bonobos

There are only a few cases of anovulatory cycles in nonhuman primates reported in the literature [52, 125, 126]. During our study, we observed eight swelling cycles with no indication that ovulation occurred, yet females displayed sexual swellings with MSPs that resembled patterns of normal, ovulatory cycles (see Fig. 6 and Additional file 3: Table S8). Generally, changes in the size of female sexual swellings are regulated by oestradiol and progesterone [41, 43]. In the absence of a pronounced pre-ovulatory rise in oestrogen, swelling tumescence may result from changes in oestrogen and progestin receptor concentration in the sexual skin [127], in combination with small fluctuations in oestrogen and progesterone levels. In comparison to other species which do not display sexual swellings decoupled from ovulatory cycles, the receptor sensitivity and density of bonobo sexual swellings may be different, or may fluctuate in a different way in relation to hormone excretion. Furthermore, it is possible that oestrogen is metabolised at or near oestrogen receptors in target tissues, e.g., in female sexual swellings, but that this metabolism is not reflected in urinary measurements of oestrogen metabolites [128].

In addition to bonobo sexual swellings being a relatively weak intracycle signal of the fecund phase, the presence of tumescent sexual swellings when females are not ovulating exemplifies the low intercycle reliability of this signal across the interbirth interval. As reported in other studies [62, 65], females displayed maximally tumescent sexual swellings during periods of gestation. They also displayed MSPs as early as three months (this study) to eight months [62] following parturition, during which time ovulation is unlikely. These findings parallel studies that found that female primates displayed situation-dependent sexual swellings during certain events, e.g., encounters with strange males or group takeovers by a new male [47, 129]. Other researchers have proposed that sexual swellings may function as a social passport during intergroup transfer and immigration to enhance social integration of females [130132]. Since young, nulliparous females often do not give birth for several years after immigration, it is possible that they are displaying sexual swelling cycles during this time without ovulating, or that the quality of these cycles is not sufficient for conception. Females may display maximally tumescent swellings which are decoupled from ovulation during these situations to appear receptive and attractive to males [133]. If tumescent sexual swellings are perceived to signal female fecundity and the ability to conceive, they could facilitate social interactions and integration with both males and females, e.g., mothers of potential male mating partners. Furthermore, sexual swellings decoupled from ovulation may be used in a strategic way by females with young, dependent offspring, enabling females to appear sexually attractive and receptive to males without incurring the risk of conceiving. Given these potential social functions of sexual swellings, it is possible that sexual swellings during periods of low or zero fecundity may not have been selected against, as they might facilitate female immigration and social interactions.

High variability in the timing of ovulation results in low predictability of ovulation and fecundity

Previous studies spanning several species of primates (e.g., Macaca nigra, Papio cynocephalus: reviewed in [28]; Pan troglodytes spp: [16, 134, 135]) reported some variation in the temporal relation between ovulation and sexual swellings; however, ovulation almost always occurred during the second half of the MSP. A study on captive bonobos found greater variability in the timing of ovulation relative to patterns of sexual swelling; however, the variability was limited to the second half of the MSP and post-detumescence [53]. Our results from wild bonobos show even more variability in the timing of ovulation, with ovulation occurring before, after, or on almost any day of the MSP. Given that ovulation occurred during the MSP in only 52.9 % of the analysed swelling cycles, female bonobos appear to be an extreme example of variability in the timing of ovulation relative to the sexual swelling signal. If we conceptualise the signal reliability of sexual swellings as a continuum, species with swellings that reliably or accurately signal the timing of ovulation would be distributed at one end of the continuum. Our findings suggest that wild bonobos occupy a position towards the opposite end of the continuum, where sexual swellings indicate ovulation with much less reliability and accuracy than in other species.

Consequently, the day-specific probabilities of ovulation and fecundity for female bonobos were very low, especially when compared to the same probabilities in female Western chimpanzees at Taï [16] (Fig. 7b). At its highest calculated value (0.24), the probability of a female bonobo being fecund, i.e., able to conceive, was two and a half times lower than in other species of primates, e.g., Pan troglodytes verus: 0.64 [16] and Hylobates lar: 0.73 [38]. The low predictability of ovulation in wild bonobos may hinder male mate guarding of females, especially when several females display MSPs simultaneously [133, 136]. Overlap in females’ MSPs often occurs in bonobos due to the lengthy duration of the MSP within a cycle and the high number of swelling cycles within interbirth intervals of females. Reproductive synchrony and temporal overlap in female receptivity or oestrous have been found to inhibit male monopolisation potential in other species [137, 138], and may affect male mating strategies in bonobos as well.

Broader Implications

In species where the timing of ovulation within a cycle can be more accurately predicted, males may be able to mate guard and monopolise fecund females during days when they are able to conceive. In a number of primate species, mate guarding is often used as a form of indirect sexual coercion by males [139] to constrain with whom a female can mate, and thereby ensure that the mate-guarding male sires a female’s offspring. In communities of chimpanzees where ovulation usually occurs near the end of the MSP and is thus relatively predictable, some studies found high frequencies of male-male mate competition and corresponding high levels of testosterone in high-ranking males when females exhibited MSPs (e.g., [140] but see [141]).

Given the low predictability of ovulation and fecundity in wild bonobos based on sexual swelling patterns alone, males may have to attend to other cues and signals to correctly pinpoint a female’s fecund phase and time their mating efforts effectively. Chemosignals [142144], behavioural cues [125], and vocal cues [145, 146] may play ancillary roles in signalling female fecundity in bonobos, as has been found in other species. If males are using other signals and focussing their mating efforts accordingly, this may offer one explanation as to why a high proportion of copulations are observed outside the MSP in bonobos, compared to other species [58, 60, 147]. Our results show that ovulation occurred outside the MSP in over 30 % of the ovulatory cycles we analysed. If males are able to detect ovulation by behavioural, olfactory, or other signals, this could cause males to solicit females outside the MSP and result in a weaker correlation between male solicitations and the MSP.

When ovulation is not tightly linked to the MSP, and if males are not able to discern and predict the window of fecundity in a female’s cycle using other signals, then it is more difficult and costly for a male to monopolise a female over a long period of time [148]. The costs associated with mate guarding a female throughout an extended MSP may outweigh the benefits [30]. Consequently, the high degree of uncertainty in bonobo sexual swellings may constrain or eliminate mate guarding by males. This might have resulted in male bonobos adopting alternative mating strategies to increase their mating success, e.g., deferring to females in feeding contexts, grooming females, or being perpetually willing to copulate with females. Accordingly, alternative mating strategies may lead to males investing more into affiliative relationships with females [149] rather than competing with other males for mating opportunities. In support of this theory, a recent study in this community found that high-ranking male bonobos with the highest mating success did not have elevated levels of testosterone during periods of mate competition [68]. In contrast to many communities of chimpanzees, these results suggest that male bonobos are not using intrasexual aggression to compete over potentially fecund females.

By displaying tumescent sexual swellings during extended periods of time, female bonobos may prolong the period during which they are attractive to males, and thereby increase their ability to confuse paternity [150]. Furthermore, if high costs associated with mate guarding and aggressing females lead to a decrease in these behaviours by males, females may be less constrained in expressing female mate choice compared to other species [151]. If female bonobos can manipulate male mating strategies or benefit from having sexual swellings that are less reliable and of longer duration, then one might ask why female chimpanzees do not have sexual swellings similar to those of female bonobos? The cost-of-sexual-attraction hypothesis may offer one explanation [107]. This hypothesis posits that the intensity of female sexual attractiveness, exhibited via sexual swellings, is primarily driven by levels of within-group scramble competition. If ranging in large groups and maintaining sexual swellings is costly for females due to low food availability, insufficient food between fruit trees, e.g., terrestrial herbaceous vegetation, or anthropogenic disturbances, then females might intensify their sexual attractiveness during shorter periods of time [152]. This, in turn, could lead to elevated levels of sexual coercion by males [153]. On the other hand, a low cost of grouping could lead to a lower intensity of sexual attractiveness, i.e., longer MSPs, and could increase the number of cycles per conception, resulting in a lower level of sexual coercion and more social benefits for females. Although this hypothesis seems plausible, ecological data to support it are still lacking. Empirical data pertaining to the quality and quantity of available food resources and the cost of grouping in wild bonobos will contribute to understanding how sexual swellings affect ranging and behaviour patterns in bonobos versus chimpanzees.

The extent to which male bonobos attend to female sexual swellings and use these signals to time their mating efforts and strategies remain to be investigated. Further study will assess whether there is evidence of female mate choice in wild bonobos, whether females modify their mating strategies across the ovarian cycle, and to what extent sexual swellings serve other social functions in this species.

Conclusions

We found a weak temporal relationship between sexual swellings and the timing of ovulation in wild female bonobos, which resulted in a very low day-specific probability of ovulation and fecundity during a female’s MSP. Bonobo sexual swellings appear to send mixed messages to males, as they do not always signal fecundity or imminent ovulation. Hence they are only probabilistic signals, rather than reliable indicators of ovulation. High variability in the relation between this sexual signal and the timing of ovulation may make it difficult for males to accurately time their mating efforts, if they use sexual swellings alone to assess female fecundity. By inference, it is likely that other sources of information, e.g., behavioural cues from females, trigger reproductive investment of males. If the temporal inflation and variability of sexual swellings in relation to ovulation constrains mate-guarding efforts by male bonobos, this could enable females to express mate choice without being constrained by males. Such a scenario would support models that emphasise differences in patterns of sexual conflict between bonobos and chimpanzees [1, 154].

Abbreviations

Cr, creatinine; E1, Oestrone; GLMM, generalised linear mixed model; IOI, interovulatory interval; ISI, interswelling interval; LMM, linear mixed model; MSP, maximum swelling phase; Pd, pregnanediol.

References

  1. Stumpf RM, Martinez-Mota R, Milich KM, Righini N, Shattuck MR. Sexual conflict in primates. Evol Anthropol. 2011;20:62–75.

    Article  PubMed  Google Scholar 

  2. Andersson M. Sexual Selection. Princeton: Princeton University Press; 1994.

    Google Scholar 

  3. Emlen ST, Oring LW. Ecology, sexual selection, and the evolution of mating systems. Science. 1977;197:215–23.

    Article  CAS  PubMed  Google Scholar 

  4. Clutton-Brock TH. Mammalian mating systems. Proc R Soc London Ser B Biol Sci. 1989;236:339–72.

    Article  CAS  Google Scholar 

  5. Girolami L, Bielert C. Female perineal swelling and its effects on male sexual arousal: an apparent sexual releaser in the chacma baboon (Papio ursinus). Int J Primatol. 1987;8:651–61.

    Article  Google Scholar 

  6. Funk DH, Tallamy DW. Courtship role reversal and deceptive signals in the long-tailed dance fly, Rhamphomyia longicauda. Anim Behav. 2000;59:411–21.

    Article  PubMed  Google Scholar 

  7. Gill FB. Ornithology. New York: W. H. Freeman; 2007.

    Google Scholar 

  8. Dixson AF. Primate Sexuality: Comparative Studies of the Prosimians, Monkeys, Apes and Human Beings. Oxford: Oxford University Press; 2012.

    Book  Google Scholar 

  9. Muller MN, Kahlenberg SM, Thompson ME, Wrangham RW. Male coercion and the costs of promiscuous mating for female chimpanzees. Proc R Soc B Biol Sci. 2007;274:1009–14.

    Article  CAS  Google Scholar 

  10. Wrangham RW. On the evolution of ape social systems. Soc Sci Inf. 1979;18:335–68.

    Article  Google Scholar 

  11. Cumming JM. Sexual selection and the evolution of dance fly mating systems (Diptera: Empididae; Empidinae). Can Entomol. 1994;126:907–20.

    Article  Google Scholar 

  12. Alberts SC, Fitzpatrick CL. Paternal care and the evolution of exaggerated sexual swellings in primates. Behav Ecol. 2012;23:699–706.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Hrdy SB. Infanticide among animals: a review, classification, and examination of the implications for the reproductive strategies of females. Ethol Sociobiol. 1979;1:13–40.

    Article  Google Scholar 

  14. Nunn CL, van Schaik CP, Zinner D. Do exaggerated sexual swellings function in female mating competition in primates? A comparative test of the reliable indicator hypothesis. Behav Ecol. 2001;12:646–54.

    Article  Google Scholar 

  15. Aujard F, Heistermann M, Thierry B, Hodges JK. Functional significance of behavioral, morphological, and endocrine correlates across the ovarian cycle in semifree ranging female Tonkean macaques. Am J Primatol. 1998;46:285–309.

    Article  CAS  PubMed  Google Scholar 

  16. Deschner T, Heistermann M, Hodges K, Boesch C. Timing and probability of ovulation in relation to sex skin swelling in wild West African chimpanzees, Pan troglodytes verus. Anim Behav. 2003;66:551–60.

    Article  Google Scholar 

  17. Heistermann M, Ziegler T, van Schaik CP, Launhardt K, Winkler P, Hodges JK. Loss of oestrus, concealed ovulation and paternity confusion in free-ranging Hanuman langurs. Proc R Soc London, Ser B. 2001;268:2445–51.

    Article  CAS  Google Scholar 

  18. Hockham LR, Ritchie MG. Female secondary sexual characteristics: appearances might be deceptive. Trends Ecol Evol. 2000;15:436–8.

    Article  PubMed  Google Scholar 

  19. Small MF. Female Choices: Sexual Behavior of Female Primates. Ithaca: Cornell University Press; 1995.

    Google Scholar 

  20. Davies NB, Hartley IR, Hatchwell BJ, Langmore NEE. Female control of copulations to maximize male help: a comparison of polygynandrous alpine accentors, Prunella collaris, and dunnocks, P. modularis. Anim Behav. 1996;51:27–47.

    Article  Google Scholar 

  21. Dixson AF. Observations on the evolution and behavioral significance of “sexual skin” in female primates. Adv Study Behav. 1983;13:63–106.

    Article  Google Scholar 

  22. Rowell TE. Female reproduction cycles and social behavior in primates. Adv Study Behav. 1972;4:69–105.

    Article  Google Scholar 

  23. Graham CE. Chimpanzee endometrium and sexual swelling during menstrual cycle or hormone administration. Folia Primatol. 1973;19:458–68.

    Article  CAS  PubMed  Google Scholar 

  24. Thierry B, Heistermann M, Aujard F, Hodges JK. Long-term data on basic reproductive parameters and evaluation of endocrine, morphological, and behavioral measures for monitoring reproductive status in a group of semifree-ranging Tonkean macaques (Macaca tonkeana). Am J Primatol. 1996;39:47–62.

    Article  Google Scholar 

  25. Fitzpatrick CL, Altmann J, Alberts SC. Sources of variance in a female fertility signal: exaggerated estrous swellings in a natural population of baboons. Behav Ecol Sociobiol. 2014;68:1109–22.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Dahl JF. The external genitalia of female pygmy chimpanzees. Anat Rec. 1985;211:24–8.

    Article  CAS  PubMed  Google Scholar 

  27. Dixson AF. Sexual selection and sexually dimorphic traits. In: Primate Sexuality: Comparative Studies of the Prosimians, Monkeys, Apes, and Humans. Oxford: Oxford University Press; 1998. p. 170–216.

    Google Scholar 

  28. Nunn CL. The evolution of exaggerated sexual swellings in primates and the graded-signal hypothesis. Anim Behav. 1999;58:229–46.

    Article  PubMed  Google Scholar 

  29. Stallmann RR, Froehlich JW. Primate sexual swellings as coevolved signal systems. Primates. 2000;41:1–16.

    Article  Google Scholar 

  30. Zinner D, van Schaik CP, Nunn CL, Kappeler PM. Sexual selection and exaggerated sexual swellings of female primates. In: Kappeler PM, van Schaik CP, editors. Sexual Selection in Primates: New and Comparative Perspectives. Cambridge: Cambridge University Press; 2004. p. 71–89.

    Chapter  Google Scholar 

  31. Pagel M. The evolution of conspicuous oestrous advertisement in Old World monkeys. Anim Behav. 1994;47:1333–41.

    Article  Google Scholar 

  32. Huchard E, Raymond M, Benavides J, Marshall H, Knapp LA, Cowlishaw G. A female signal reflects MHC genotype in a social primate. BMC Evol Biol. 2010;10:96.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  33. Hamilton WJI. Significance of parental investment by primates to the evolution of male-female associations. In: Taub DM, editor. Primate Paternalism. New York: Van Nostrand Reinhold; 1984. p. 309–35.

    Google Scholar 

  34. Gesquiere LR, Wango EO, Alberts SC, Altmann J. Mechanisms of sexual selection: sexual swellings and estrogen concentrations as fertility indicators and cues for male consort decisions in wild baboons. Horm Behav. 2007;51:114–25.

    Article  CAS  PubMed  Google Scholar 

  35. Brauch K, Pfefferle D, Hodges K, Möhle U, Fischer J, Heistermann M. Female sexual behavior and sexual swelling size as potential cues for males to discern the female fertile phase in free-ranging Barbary macaques (Macaca sylvanus) of Gibraltar. Horm Behav. 2007;52:375–83.

    Article  PubMed  Google Scholar 

  36. Higham JP, Heistermann M, Saggau C, Agil M, Perwitasari-Farajallah D, Engelhardt A. Sexual signalling in female crested macaques and the evolution of primate fertility signals. BMC Evol Biol. 2012;12:89.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Searcy WA, Nowicki S. Introduction. In: The Evolution of Animal Communication: Reliability and Deception in Signaling Systems. Princeton: Princeton University Press; 2005. p. 1–23.

    Google Scholar 

  38. 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. 2007;51:221–30.

    Article  CAS  PubMed  Google Scholar 

  39. Higham JP, Heistermann M, Ross C, Semple S, Maclarnon A. The timing of ovulation with respect to sexual swelling detumescence in wild olive baboons. Primates. 2008;49:295–9.

    Article  PubMed  Google Scholar 

  40. Clark G, Birch HG. Observations on the sex skin and sex cycle in the chimpanzee. Endocrinology. 1948;43:218–31.

    Article  CAS  PubMed  Google Scholar 

  41. Ozasa H, Gould KG. Demonstration and characterization of estrogen receptor in chimpanzee sex skin: correlation between nuclear receptor levels and degree of swelling. Endocrinology. 1982;111:125–31.

    Article  CAS  PubMed  Google Scholar 

  42. Graham CE, Collins DC, Robinson H, Preedy JRK. Urinary levels of estrogens and pregnanediol and plasma levels of progesterone during the menstrual cycle of the chimpanzee: relationship to the sexual swelling. Endocrinology. 1972;91:13–24.

    Article  CAS  PubMed  Google Scholar 

  43. Carlisle KS, Brenner RM, Montagna W. Hormonal regulation of sex skin in Macaca nemestrina. Biol Reprod. 1981;25:1053–63.

    Article  CAS  PubMed  Google Scholar 

  44. Rowell TE. Baboon menstrual cycles affected by social environment. J Reprod Fertil. 1970;21:133–41.

    Article  CAS  PubMed  Google Scholar 

  45. DeVore I. Male dominance and mating behaviour in baboons. In: Beach FA, editor. Sex and Behavior. New York: John Wiley and Sons; 1965. p. 266–89.

    Google Scholar 

  46. Machatschke IH, Wallner B, Dittami J. Impact of social environment on female chimpanzee reproductive cycles. Horm Behav. 2006;50:126–31.

    Article  CAS  PubMed  Google Scholar 

  47. Zinner D, Deschner T. Sexual swellings in female hamadryas baboons after male take-overs: “deceptive” swellings as a possible female counter-strategy against infanticide. Am J Primatol. 2000;52:157–68.

    Article  Google Scholar 

  48. Roberts EK, Lu A, Bergman TJ, Beehner JC. A Bruce effect in wild geladas. Science. 2012;335:1222–5.

    Article  CAS  PubMed  Google Scholar 

  49. Barelli C, Heistermann M, Boesch C, Reichard UH. Mating patterns and sexual swellings in pair-living and multimale groups of wild white-handed gibbons, Hylobates lar. Anim Behav. 2008;75:991–1001.

    Article  Google Scholar 

  50. Steinetz BG, Ducrot C, Randolph C, Mahoney CJ. Determination of the time of ovulation in chimpanzees by measurement of LH, estrone sulfate, and pregnanediol 3α-glucuronide in urine: comparison with serum hormone patterns. J Med Primatol. 1992;21:239–45.

    CAS  PubMed  Google Scholar 

  51. Jurke MH, Hagey LR, Jurke S, Czekala NM. Monitoring hormones in urine and feces of captive bonobos (Pan paniscus). Primates. 2000;41:311–9.

    Article  Google Scholar 

  52. Heistermann M, Möhle U, Vervaecke H, van Elsacker L, Hodges JK. Application of urinary and fecal steroid measurements for monitoring ovarian function and pregnancy in the bonobo (Pan paniscus) and evaluation of perineal swelling patterns in relation to endocrine events. Biol Reprod. 1996;55:844–53.

    Article  CAS  PubMed  Google Scholar 

  53. Reichert KE, Heistermann M, Hodges JK, Boesch C, Hohmann G. What females tell males about their reproductive status: are morphological and behavioural cues reliable signals of ovulation in bonobos (Pan paniscus)? Ethology. 2002;108:583–600.

    Article  Google Scholar 

  54. Sannen A, Van Elsacker L, Heistermann M, Eens M. Certain aspects of bonobo female sexual repertoire are related to urinary testosterone metabolite levels. Folia Primatol. 2005;76:21–32.

    Article  PubMed  Google Scholar 

  55. Dahl JF, Nadler RD, Collins DC. Monitoring the ovarian cycles of Pan troglodytes and P. paniscus: a comparative approach. Am J Primatol. 1991;24:195–209.

    Article  Google Scholar 

  56. Hohmann G, Fruth B. Use and function of genital contacts among female bonobos. Anim Behav. 2000;60:107–20.

    Article  PubMed  Google Scholar 

  57. Dahl JF. Cyclic perineal swelling during the intermenstrual intervals of captive female pygmy chimpanzees (Pan paniscus). J Hum Evol. 1986;15:369–85.

    Article  Google Scholar 

  58. Savage-Rumbaugh ES, Wilkerson BJ. Socio-sexual behavior in Pan paniscus and Pan troglodytes: a comparative study. J Hum Evol. 1978;7:327–44.

    Article  Google Scholar 

  59. Furuichi T. Sexual swelling, receptivity, and grouping of wild pygmy chimpanzee females at Wamba, Zaire. Primates. 1987;28:309–18.

    Article  Google Scholar 

  60. Thompson-Handler N, Malenky R, Badrian N. Sexual behavior of Pan paniscus under natural conditions in the Lomako Forest, Equateur, Zaire. In: Susman RL, editor. The Pygmy Chimpanzee: Evolutionary Biology and Behavior. New York: Plenum Press; 1984. p. 347–68.

    Chapter  Google Scholar 

  61. Nadler RD, Graham CE, Collins DC, Kling OR. Postpartum amenorrhea and behavior of apes. In: Graham CE, editor. Reproductive Biology of the Great Apes: Comparative and Biomedical Perspectives. New York: Academic; 1981. p. 69–81.

    Google Scholar 

  62. Furuichi T, Hashimoto C. Why female bonobos have a lower copulation rate during estrus than chimpanzees. In: Boesch C, Hohmann G, Marchant LF, editors. Behavioural Diversity in Chimpanzees and Bonobos. Cambridge: Cambridge University Press; 2002. p. 156–67.

    Chapter  Google Scholar 

  63. Surbeck M, Hohmann G. Intersexual dominance relationships and the influence of leverage on the outcome of conflicts in wild bonobos (Pan paniscus). Behav Ecol Sociobiol. 2013;67:1767–80.

    Article  Google Scholar 

  64. Paoli T. The absence of sexual coercion in bonobos. In: Muller MN, Wrangham RW, editors. Sexual Coercion in Primates and Humans. Cambridge: Harvard University Press; 2009. p. 410–23.

    Google Scholar 

  65. Kano T. The Last Ape: Pygmy Chimpanzee Behavior and Ecology. Stanford: Stanford University Press; 1992.

    Google Scholar 

  66. Hohmann G, Fruth B. Lui Kotal - A new site for field research on bonobos in the Salonga National Park. Pan Africa News. 2003;10:25–7.

    Google Scholar 

  67. Hohmann G, Fruth B. Dynamics in social organization of bonobos (Pan paniscus). In: Boesch C, Hohmann G, Marchant LF, editors. Behavioural Diversity in Chimpanzees and Bonobos. Cambridge: Cambridge University Press; 2002. p. 138–50.

    Chapter  Google Scholar 

  68. Surbeck M, Deschner T, Schubert G, Weltring A, Hohmann G. Mate competition, testosterone and intersexual relationships in bonobos, Pan paniscus. Anim Behav. 2012;83:659–69.

    Article  Google Scholar 

  69. Hauser B, Deschner T, Boesch C. Development of a liquid chromatography–tandem mass spectrometry method for the determination of 23 endogenous steroids in small quantities of primate urine. J Chromatogr B. 2008;862:100–12.

    Article  CAS  Google Scholar 

  70. Bahr NI, Palme R, Möhle U, Hodges JK, Heistermann M. Comparative aspects of the metabolism and excretion of cortisol in three individual nonhuman primates. Gen Comp Endocrinol. 2000;117:427–38.

    Article  CAS  PubMed  Google Scholar 

  71. Jaffé M. Ueber den Niederschlag, welchen Pikrinsäure im normalen harn erzeugt, und über eine neue Reaction des Kreatinins. Zeitschrift für Physiol Chemie. 1886;10:391–400.

    Google Scholar 

  72. Jeffcoate SL. Use of rapid hormone assays in the prediction of ovulation. In: Jeffcoate SL, editor. Ovulation, Methods for Its Prediction and Detection. Chichester: J. Wiley; 1983. p. 67–82.

    Google Scholar 

  73. Carosi M, Heistermann M, Visalberghi E. Display of proceptive behaviors in relation to urinary and fecal progestin levels over the ovarian cycle in female tufted capuchin monkeys. Horm Behav. 1999;36:252–65.

    Article  CAS  PubMed  Google Scholar 

  74. Munro CJ, Stabenfeldt GH, Cragun JR, Addiego LA, Overstreet JW, Lasley BL. Relationship of serum estradiol and progesterone concentrations to the excretion profiles of their major urinary metabolites as measured by enzyme immunoassay and radioimmunoassay. Clin Chem. 1991;37:838–44.

    CAS  PubMed  Google Scholar 

  75. Behboodi E, Katz DF, Samuels SJ, Tell L, Hendrickx AG, Lasley BL. The use of a urinary estrone conjugates assay for detection of optimal mating time in the cynomolgus macaque (Macaca fascicularis). J Med Primatol. 1991;20:229–34.

    CAS  PubMed  Google Scholar 

  76. Wilcox AJ, Weinberg CR, Baird DD. Timing of sexual intercourse in relation to ovulation. N Engl J Med. 1995;333:1517–21.

    Article  CAS  PubMed  Google Scholar 

  77. Ferreira-Poblete A. The probability of conception on different days of the cycle with respect to ovulation: an overview. Adv Contracept. 1997;13:83–95.

    Article  CAS  PubMed  Google Scholar 

  78. Zinner D, Alberts SC, Nunn CL, Altmann J. Significance of primate sexual swellings. Nature. 2002;420:142–3.

    Article  CAS  PubMed  Google Scholar 

  79. Baayen RH. Analyzing Linguistic Data: A Practical Introduction to Statistics Using R. Cambridge: Cambridge University Press; 2008.

    Book  Google Scholar 

  80. Bolker BM, Brooks ME, Clark CJ, Geange SW, Poulsen JR, Stevens MHH, et al. Generalized linear mixed models: a practical guide for ecology and evolution. Trends Ecol Evol. 2008;24:127–35.

  81. R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2016. https://www.R-project.org/. Accessed 6 Apr 2016.

  82. Bates D, Maechler M, Bolker B, Walker S. lme4: Linear mixed-effects models using Eigen and S4. R package version 1.1-7. 2014.

    Google Scholar 

  83. Quinn GP, Keough MJ. Experimental Design and Data Analysis for Biologists. Cambridge: Cambridge University Press; 2002.

    Book  Google Scholar 

  84. Field A. Discovering Statistics Using SPSS. London: Sage Publications; 2005.

    Google Scholar 

  85. Fox J, Weisberg S. An R Companion to Applied Regression. 2nd ed. Thousand Oaks: Sage Publications; 2011.

    Google Scholar 

  86. Forstmeier W, Schielzeth H. Cryptic multiple hypotheses testing in linear models: overestimated effect sizes and the winner’s curse. Behav Ecol Sociobiol. 2011;65:47–55.

    Article  PubMed  Google Scholar 

  87. Dobson AJ. An Introduction to Generalized Linear Models. 2nd ed. London: Chapman and Hall; 2002.

    Google Scholar 

  88. Barr DJ, Levy R, Scheepers C, Tily HJ. Random effects structure for confirmatory hypothesis testing: keep it maximal. J Mem Lang. 2013;68:255–78.

    Article  Google Scholar 

  89. Nieuwenhuis R, te Grotenhuis M, Pelzer B. influence.ME: tools for detecting influential data in mixed effects models. R J. 2012;4:38–47.

    Google Scholar 

  90. Thompson ME, Wrangham RW, Reynolds V. Urinary estrone conjugates and reproductive parameters in Kibale (Kanyawara) and Budongo (Sonso) chimpanzees. In: Newton-Fisher NE, Notman H, Paterson JD, Reynolds V, editors. Primates of Western Uganda. New York: Springer Science & Business Media; 2006. p. 227–45.

    Chapter  Google Scholar 

  91. Smuts B, Nicolson N. Reproduction in wild female olive baboons. Am J Primatol. 1989;19:229–46.

    Article  Google Scholar 

  92. Dunbar RIM. Determinants and evolutionary consequences of dominance among female gelada baboons. Behav Ecol Sociobiol. 1980;7:253–65.

    Article  Google Scholar 

  93. Douglas PH, Ngonga Ngomo AC, Hohmann G. ADAGIO: Approach for Dominance Assessment in Gregarious Species. http://ngonga.github.io/adagio. Accessed 6 Apr 2016.

  94. Schielzeth H, Forstmeier W. Conclusions beyond support: overconfident estimates in mixed models. Behav Ecol. 2009;20:416–20.

    Article  PubMed  Google Scholar 

  95. Anderson CM, Bielert CF. Adolescent exaggeration in female Catarrhine primates. Primates. 1994;35:283–300.

    Article  Google Scholar 

  96. Wallis J. A survey of reproductive parameters in the free-ranging chimpanzees of Gombe National Park. J Reprod Fertil. 1997;109:297–307.

    Article  CAS  PubMed  Google Scholar 

  97. Surbeck M, Deschner T, Behringer V, Hohmann G. Urinary C-peptide levels in male bonobos (Pan paniscus) are related to party size and rank but not to mate competition. Horm Behav. 2015;71:22–30.

    Article  CAS  PubMed  Google Scholar 

  98. O’Malley RC, Stanton MA, Gilby IC, Lonsdorf EV, Pusey A, Markham AC, et al. Reproductive state and rank influence patterns of meat consumption in wild female chimpanzees (Pan troglodytes schweinfurthii). J Hum Evol. 2016;90:16–28.

  99. Thompson ME, Muller MN, Wrangham RW. The energetics of lactation and the return to fecundity in wild chimpanzees. Behav Ecol. 2012;23:1234–41.

    Article  Google Scholar 

  100. Fahy GE, Richards MP, Fuller BT, Deschner T, Hublin JJ, Boesch C. Stable nitrogen isotope analysis of dentine serial sections elucidate sex differences in weaning patterns of wild chimpanzees (Pan troglodytes). Am J Phys Anthropol. 2014;153:635–42.

    Article  PubMed  Google Scholar 

  101. Bowman LA, Dilley SR, Keverne EB. Suppression of oestrogen-induced LH surges by social subordination in talapoin monkeys. Nature. 1978;275:56–8.

    Article  CAS  PubMed  Google Scholar 

  102. Furuichi T. Agonistic interactions and matrifocal dominance rank of wild bonobos (Pan paniscus) at Wamba. Int J Primatol. 1997;18:855–75.

    Article  Google Scholar 

  103. Thompson ME, Wrangham RW. Diet and reproductive function in wild female chimpanzees (Pan troglodytes schweinfurthii) at Kibale National Park, Uganda. Am J Phys Anthropol. 2008;135:171–81.

    Article  PubMed  Google Scholar 

  104. Knott C. Female reproductive ecology of the apes: implications for human evolution. In: Ellison PT, editor. Reproductive Ecology and Human Evolution. New York: Walter de Gruyter; 2001. p. 429–63.

    Google Scholar 

  105. Ziegler TE, Strier KB, Van Belle S. The reproductive ecology of South American primates: ecological adaptations in ovulation and conception. In: Garber PA, Estrada A, Bicca-Marques JC, Heymann EW, Strier KB, editors. South American Primates: Comparative Perspectives in the Study of Behavior, Ecology, and Conservation. New York: Springer; 2009. p. 191–210.

    Chapter  Google Scholar 

  106. Hohmann G, Fruth B. Food sharing and status in unprovisioned bonobos. In: Wiessner P, Schiefenhövel W, editors. Food and the Status Quest. Providence & Oxford: Berghahn Books; 1996. p. 47–67.

    Google Scholar 

  107. Wrangham R. The cost of sexual attraction: is there a trade-off in female Pan between sex appeal and received coercion? In: Boesch C, Hohmann G, Marchant L, editors. Behavioural Diversity in Chimpanzees and Bonobos. Cambridge: Cambridge University Press; 2002. p. 204–15.

    Chapter  Google Scholar 

  108. Hohmann G, Fowler A, Sommer V, Ortmann S. Frugivory and gregariousness of Salonga bonobos and Gashaka chimpanzees. In: Hohmann G, Robbins M, Christophe B, editors. Feeding Ecology in Apes and Other Primates. Cambridge: Cambridge University Press; 2006. p. 123–59.

    Google Scholar 

  109. Vervaecke H, van Elsackar L, Möhle U, Heistermann M, Verheyen RF. Inter-menstrual intervals in captive bonobos (Pan paniscus). Primates. 1999;40:283–9.

    Article  Google Scholar 

  110. Thompson-Handler NE. The pygmy chimpanzee: sociosexual behavior, reproductive biology and life history patterns. PhD Diss. Yale University. 1990

  111. Neugebauer W. The status and management of the Pygmy chimpanzee Pan paniscus in European zoos. Int Zoo Yearb. 1980;20:64–70.

    Article  Google Scholar 

  112. Lenton EA, Landgren B-M, Sexton L, Harper R. Normal variation in the length of the follicular phase of the menstrual cycle: effect of chronological age. Br J Obstet Gynaecol. 1984;91:681–4.

    Article  CAS  PubMed  Google Scholar 

  113. Reichert K. Linking behaviour and physiology of female bonobos (Pan paniscus). PhD Diss. Universität Bayreuth. 2005

  114. Liu Y, Gold EB, Lasley BL, Johnson WO. Factors affecting menstrual cycle characteristics. Am J Epidemiol. 2004;160:131–40.

    Article  CAS  PubMed  Google Scholar 

  115. Wasserman MD, Chapman CA, Milton K, Gogarten JF, Wittwer DJ, Ziegler TE. Estrogenic plant consumption predicts red colobus monkey (Procolobus rufomitratus) hormonal state and behavior. Horm Behav. 2012;62:553–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  116. Higham JP, Ross C, Warren Y, Heistermann M, MacLarnon AM. Reduced reproductive function in wild baboons (Papio hamadryas anubis) related to natural consumption of the African black plum (Vitex doniana). Horm Behav. 2007;52:384–90.

    Article  CAS  PubMed  Google Scholar 

  117. Rowell TE. Social Behaviour of Monkeys. Harmondsworth: Penguin; 1972.

    Google Scholar 

  118. Adams MR, Kaplan JR, Koritnik DR. Psychosocial influences on ovarian endocrine and ovulatory function in Macaca fascicularis. Physiol Behav. 1985;35:935–40.

    Article  CAS  PubMed  Google Scholar 

  119. Beehner JC, Lu A. Reproductive suppression in female primates: A review. Evol Anthropol. 2013;22:226–38.

    Article  PubMed  Google Scholar 

  120. Parga JA, Lessnau RG. Female age-specific reproductive rates, birth seasonality, and infant mortality of ring-tailed lemurs on St. Catherines Island: 17-year reproductive history of a free-ranging colony. Zoo Biol. 2005;24:295–309.

    Article  Google Scholar 

  121. Matsumoto-Oda A. Mahale chimpanzees: grouping patterns and cycling females. Am J Primatol. 1999;47:197–207.

    Article  CAS  PubMed  Google Scholar 

  122. Tutin CEG, McGinnis PR. Chimpanzee reproduction in the wild. In: Graham CE, editor. Reproductive Biology of the Great Apes: Comparative and Biomedical Perspectives. New York: Academic; 1981. p. 239–64.

    Google Scholar 

  123. Butler H. Evolutionary trends in primate sex cycles. Contrib Primatol. 1974;3:2–35.

    CAS  PubMed  Google Scholar 

  124. Shelmidine N, Borries C, Koenig A. Genital swellings in silvered langurs: what do they indicate? Am J Primatol. 2007;69:519–32.

    Article  PubMed  Google Scholar 

  125. Engelhardt A, Hodges JK, Niemitz C, Heistermann M. Female sexual behavior, but not sex skin swelling, reliably indicates the timing of the fertile phase in wild long-tailed macaques (Macaca fascicularis). Horm Behav. 2005;47:195–204.

    Article  CAS  PubMed  Google Scholar 

  126. Dailey RA, Neill JD. Seasonal variation in reproductive hormones of rhesus monkeys: anovulatory and short luteal phase menstrual cycles. Biol Reprod. 1981;25:560–7.

    Article  CAS  PubMed  Google Scholar 

  127. Onouchi T, Kato J. Estrogen receptors and estrogen-inducible progestin receptors in the sexual skin of the monkey. J Steroid Biochem. 1983;18:145–51.

    Article  CAS  PubMed  Google Scholar 

  128. Zhu BT, Conney AH. Functional role of estrogen metabolism in target cells: review and perspectives. Carcinogenesis. 1998;19:1–27.

    Article  PubMed  Google Scholar 

  129. Hrdy SB, Whitten PL. Patterning of sexual activity. In: Smuts BB, Cheney DL, Seyfarth RW, Wrangham RW, Struhsaker TT, editors. Primate Societies. Chicago: University of Chicago Press; 1987. p. 370–84.

    Google Scholar 

  130. Pusey A. Intercommunity transfer of chimpanzees in Gombe National Park. In: Hamburg DA, McCown ER, editors. The Great Apes. Menlo Park: Benjamin Cummings; 1979. p. 465–79.

    Google Scholar 

  131. Goodall J. The Chimpanzees of Gombe: Patterns of Behavior. Cambridge: The Belknap Press of Harvard University Press; 1986.

    Google Scholar 

  132. Deschner T, Boesch C. Can the patterns of sexual swelling cycles in female Taï chimpanzees be explained by the cost-of-sexual-attraction hypothesis? Int J Primatol. 2007;28:389–406.

    Article  Google Scholar 

  133. Surbeck M, Mundry R, Hohmann G. Mothers matter! Maternal support, dominance status and mating success in male bonobos (Pan paniscus). Proc R Soc B Biol Sci. 2011;278:590–8.

    Article  Google Scholar 

  134. Emery MA, Whitten PL. Size of sexual swellings reflects ovarian function in chimpanzees (Pan troglodytes). Behav Ecol Sociobiol. 2003;54:340–51.

    Article  Google Scholar 

  135. Machatschke IH, Dittami J, Wallner B. Morphometric and hormonal changes during the chimpanzee menstrual cycle. J Med Primatol. 2006;35:331–40.

    Article  CAS  PubMed  Google Scholar 

  136. Hohmann G, Fruth B. Intra- and inter-sexual aggression by bonobos in the context of mating. Behaviour. 2003;140:1389–413.

    Article  Google Scholar 

  137. Carnes LM, Nunn CL, Lewis RJ. Effects of the distribution of female primates on the number of males. PLoS ONE. 2011;6:e19853.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  138. Nunn CL. The number of males in primate social groups: a comparative test of the socioecological model. Behav Ecol Sociobiol. 1999;46:1–13.

    Article  Google Scholar 

  139. Muller MN, Thompson ME, Kahlenberg SM, Wrangham RW. Sexual coercion by male chimpanzees shows that female choice may be more apparent than real. Behav Ecol Sociobiol. 2011;65:921–33.

    Article  Google Scholar 

  140. Muller MN, Wrangham RW. Dominance, aggression and testosterone in wild chimpanzees: a test of the “challenge hypothesis”. Anim Behav. 2004;67:113–23.

  141. Klinkova E, Heistermann M, Hodges JK. Social parameters and urinary testosterone level in male chimpanzees (Pan troglodytes). Horm Behav. 2004;46:474–81.

    Article  CAS  PubMed  Google Scholar 

  142. Drea CM. D’scent of man: a comparative survey of primate chemosignaling in relation to sex. Horm Behav. 2015;68:117–33.

    Article  CAS  PubMed  Google Scholar 

  143. Ziegler TE, Epple G, Snowdon CT, Porter TA, Belcher AM, Küderling I. Detection of the chemical signals of ovulation in the cotton-top tamarin, Saguinus oedipus. Anim Behav. 1993;45:313–22.

    Article  Google Scholar 

  144. Michael RP, Keverne EB. Pheromones and the communication of sexual status in primates. Nature. 1968;218:746–9.

    Article  CAS  PubMed  Google Scholar 

  145. O’Connell SM, Cowlishaw G. Infanticide avoidance, sperm competition and mate choice: the function of copulation calls in female baboons. Anim Behav. 1994;48:687–94.

    Article  Google Scholar 

  146. Semple S, McComb K. Perception of female reproductive state from vocal cues in a mammal species. Proc R Soc London B Biol Sci. 2000;267:707–12.

    Article  CAS  Google Scholar 

  147. Furuichi T, Hashimoto C. Sex differences in copulation attempts in wild bonobos at Wamba. Primates. 2004;45:59–62.

    Article  PubMed  Google Scholar 

  148. Girard-Buttoz C, Heistermann M, Rahmi E, Agil M, Fauzan PA, Engelhardt A. Costs of and investment in mate-guarding in wild long-tailed macaques (Macaca fascicularis): influences of female characteristics and male–female social bonds. Int J Primatol. 2014;35:701–24.

    Article  PubMed  PubMed Central  Google Scholar 

  149. Smuts BB. Sex and Friendship in Baboons. New York: Aldine; 1985.

    Google Scholar 

  150. Higham JP, Semple S, MacLarnon A, Heistermann M, Ross C. Female reproductive signaling, and male mating behavior, in the olive baboon. Horm Behav. 2009;55:60–7.

    Article  CAS  PubMed  Google Scholar 

  151. Furuichi T. Female contributions to the peaceful nature of bonobo society. Evol Anthropol. 2011;20:131–42.

    Article  PubMed  Google Scholar 

  152. Milich KM, Bahr JM, Stumpf RM, Chapman CA. Timing is everything: expanding the cost of sexual attraction hypothesis. Anim Behav. 2014;88:219–24.

    Article  Google Scholar 

  153. Watts DP. Coalitionary mate guarding by male chipanzees at Ngogo, Kibale National Park, Uganda. Behav Ecol Sociobiol. 1998;44:43–55.

    Article  Google Scholar 

  154. Thompson ME. Sexual conflict: Nice guys finish last. Curr Biol. 2014;24:R1125–7.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

We thank the ICCN for permission to conduct research in the forest near Salonga National Park and for granting permission to export urine samples. Special thanks go to the people of the village of Lompole for granting access to the forest of their ancestors. We sincerely thank Lucas Goldstone, Max Kölbl, Timothy Lewis-Bale, Apoluke Mobembo, and Pauline Toni for assistance with data collection; Liza R. Moscovice and Niina Nurmi for contributing behavioural data towards the assessment of female ranks; Silke Atmaca for assistance with data processing; Vera Schmeling for assistance in the lab; and Roger Mundry and Colleen Stephens for assistance with statistical analysis. We acknowledge with gratitude Verena Behringer, Frances D. Burton, Roger Mundry, Martin Surbeck, and two anonymous reviewers for valuable comments on earlier drafts.

Funding

We gratefully acknowledge financial support from the L.S.B. Leakey Foundation (to PHD) and the Max Planck Society.

Availability of data and materials

Supporting data can be made available upon request.

Authors’ contributions

PHD, GH, and TD conceived and designed the study. PHD and RTB collected the data. PHD and RM performed the laboratory work (hormone extraction and measurement, respectively). PHD analysed the data. All authors contributed to writing the manuscript and approved the final version of it.

Competing interests

The authors declare that they have no competing interests.

Ethics approval and consent to participate

All methods used to collect observational data and urine samples were non-invasive, were in compliance with the regulations and guidelines of the Institut Congolais pour la Conservation de la Nature (ICCN), and adhered to the legal requirements of the Democratic Republic of Congo (DRC). Permits for the export of urine samples were issued by the ICCN and approved by the head veterinarian of the DRC. Import permits were obtained from the head veterinarian of Saxonia, Germany.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Pamela Heidi Douglas.

Additional files

Additional file 1:

Further details of hormone extraction and measurement. (PDF 334 kb)

Additional file 2:

Summary of results from the LMMs and GLMMs: Table S1. the ISI duration GLMM; Tables S2–S4. the three MSP duration LMMs; Table S5. the female rank and occurrence of ovulation GLMM; Table S6. the female rank and timing of ovulation GLMM. (PDF 331 kb)

Additional file 3:

Supplemental information: Table S7. Details of ovulatory cycles; Table S8. Details of anovulatory cycles. (PDF 209 kb)

Rights and permissions

Open Access This 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Douglas, P.H., Hohmann, G., Murtagh, R. et al. Mixed messages: wild female bonobos show high variability in the timing of ovulation in relation to sexual swelling patterns. BMC Evol Biol 16, 140 (2016). https://0-doi-org.brum.beds.ac.uk/10.1186/s12862-016-0691-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12862-016-0691-3

Keywords