# Sex and grooming as exchange commodities in female bonobos' daily biological market

Simone Anzà, Elisa Demuru, Elisabetta Palagi

► To cite this version:

Simone Anzà, Elisa Demuru, Elisabetta Palagi. Sex and grooming as exchange commodities in female bonobos' daily biological market. Scientific Reports, 2021, 11, pp.19344. 10.1038/s41598-021-00136-6. hal-04680591

HAL Id: hal-04680591

https://hal.science/hal-04680591v1

Submitted on 28 Aug 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for updates

OPEN

# Sex and grooming as exchange commodities in female bonobos' daily biological market

Simone Anzà<fnref n="1" />,⁶✉, Elisa Demuru²,³,⁶✉ & Elisabetta Palagi⁴,⁵

The Biological Market Theory (BMT) posits that cooperation between non-human animals can be seen as a mutually beneficial exchange of commodities similarly to what observed in human economic markets. Positive social interactions are commodities in non-human animals, and mutual exchanges fulfilling the criteria of the BMT have been shown in several species. However, the study of biological markets suffers from methodological limitations that are mainly linked to the difficulty of clearly identifying the currencies and their exchanges in the short-term. Here, we test whether bonobo females are more attractive during their maximum swelling phase, whether they exchange grooming and Genito-Genital Rubbing (GGR) on a daily level of analysis, and whether these daily exchanges fulfil the BMT criteria. Females engaged more in GGR when their sexual swelling was in the maximum phase. Moreover, they exchanged grooming and sex according to the daily "market fluctuations" associated with swelling status. Females in the minimum phase (low-value) increased their probability to engage in GGR with females in the maximum phase (high-value) by grooming them preferentially. In line with the supply/demand law, the female grooming strategy varied depending on the daily number of swollen females present: the higher the number of swollen females, the lower the individual grooming preference. As a whole, our study confirms BMT as a valid model to explain daily commodity exchanges as a function of the temporary value of traders, and underlines the importance of a day-by-day approach to unveil the presence of a biological market when the value of traders frequently changes.

The Biological Market Theory (BMT) posits that cooperation between non-human animals can be seen as a mutually beneficial exchange of commodities (i.e. goods or services) similarly to what is observed in human economic markets<fnref n="1" />. Within this context, a commodity must be monopolizable and cannot be obtained by force, so that individuals should develop an alternative strategy to conflict to get the commodity². Within the same class of traders, individuals are chosen as partners depending on the commodities available on the market, therefore the best strategy is either to increase the value of the commodity offered, or to demand less in return in order to outcompete the other traders³. As it occurs in human markets⁴, the exchange rate of commodities in a Biological Market depends on the law of supply and demand: when a certain commodity is rare, competition over it is high, as its value is perceived as being high as well. Under these circumstances, the other individuals are inclined to "pay more" to get the rare high-value commodity³.

Social interactions are a valuable commodity in non-human animals, and mutual exchanges fulfilling the criteria of the BMT have been shown in several species (cleaners-client fishes⁵; red bishops, *Euplectes orix*⁶; wood mice, *Apodemus sylvaticus*⁷; meerkats, *Suricata suricatta*⁸).

The BMT has also been used to explain several types of short- and long-term exchanges in non-human primates (hereafter primates), with most studies focusing on allo-grooming as a commodity to exchange (e.g., for sex in sifaka, *Propithecus verreauxi*⁹; for sex in gibbons, *Hylobates lar*¹⁰; for sex in chimpanzees, *Pan troglodytes*¹¹; for grooming in bonobos, *Pan paniscus*¹²). Allo-grooming is a frequent and valuable behaviour among primates where it can serve the functions of reinforcing social bonds¹³,¹⁴, fostering coalitions¹⁵, and reducing tension¹⁶. As highlighted by Sánchez-Amaro and Amici¹⁷, the study of biological markets in primates suffers

<footnote n="1">¹Behavioral Ecology Department, University of Goettingen, Goettingen, Germany. ²Laboratoire Dynamique Du Langage, CNRS UMR 5596, University of Lyon 2, Lyon, France. ³Equipe de Neuro-Ethologie Sensorielle, ENES/CRNL, CNRS UMR 5292, Inserm UMR S1028, University of Lyon/Saint-Etienne, Saint-Etienne, France. ⁴Unit of Ethology, Department of Biology, University of Pisa, Pisa, Italy. ⁵Natural History Museum, University of Pisa, Pisa, Italy. ⁶These authors contributed equally: Simone Anzà and Elisa Demuru. ⁷email: simone.anza@uni-goettingen.de; elisa.demuru@cnrs.fr</footnote>
from theoretical and methodological limitations that are mainly linked to the difficulty of identifying clear exchanges in cognitively and socially complex species. Indeed, primates are capable of keeping track of both short- and long-time exchanges¹⁸ and live in markets where social commodities are not only highly variable but also strongly interconnected¹⁹. Concerning the temporal scale, there is a total lack of studies focussing on the exchanges on a daily level of analysis.

Here, we focus on the bonobo (*Pan paniscus*) to fill this gap by exploring the possible existence of a biological market regulating the exchanges of social interactions in females living in two stable social groups. To reach this goal, we focussed our analyses on the daily fluctuation of grooming and socio-sexual interactions as a function of the swelling status of actor and receiver. As far as we know, it has never been demonstrated that the Biological Market Theory explains the exchange of commodities in bonobos, although it has been hypothesized²⁰.

Together with chimpanzees (*Pan troglodytes*), bonobos are human closest living relatives. The two *Pan* species live in multi-male/multi-female social groups characterized by fission-fusion dynamics and show male philopatry and female dispersal²¹. This dispersion pattern causes a higher level of relatedness among males than among females, nevertheless bonobo females are highly cohesive and share strong affiliative bonds²²⁻²⁴. Female cohesiveness and affiliation allow bonobo females to be dominant over males—or co-dominant according to some authors²⁵—which is a notable exception in a male-philopatric species²⁶. Female bonds are established and maintained through a wide array of affiliative behaviour such as play²⁷, consolation²⁸, agonistic support²³, allo-grooming²⁶, face engagement²⁹, and socio-sexuality³⁰. Bonobos are well known for making use of sexual behaviour for social purposes and, particularly, to decrease social tension during competitive contexts, such as food distribution or aggressive interactions³¹. Bonobo females frequently engage in Genito-Genital Rubbing (GGR), a socio-sexual behaviour during which females embrace each-other, frequently face to face, and rub their genitals by moving their hips side to side³². Recently, Moscovice et al.³⁰ showed that GGR plays an essential role in enhancing female cooperation as this behaviour stimulates the release of oxytocin. Moreover, this face-to-face interaction favours the perception and replication of facial expressions through the phenomenon of Rapid Facial Mimicry (RFM³³). RFM is an unconscious phenomenon that represents a basic layer of empathy³⁴ and is a powerful mechanism that helps establishing an emotional bridge between two subjects.

Socio-sexual interactions involving females are enhanced by the presence of the maximum sexual swelling, a conspicuous and temporary signal that in bonobos is highly attractive for both males and females³⁵,³⁶. The ano-genital region of female bonobos goes through changes in size and turgidity along the menstrual cycle, and the attractiveness changes consequently, reaching its peak during the maximum swelling phase³⁵,³⁶. Compared to other species, bonobos' maximum sexual swelling is not a reliable signal of ovulation³⁷ and can even be present during anovulatory periods, such as pregnancy and lactation²⁶. It has been proposed that bonobos' maximum swelling features—together with a posture possibly functioning as amplifier signal³⁸—might have been selected to increase female-female sexual interactions³⁵. Through socio-sexuality, bonobo females get integrated in a group, establish affiliative bonds and increase social cohesiveness which, in turn, allow them to form coalitions that increase their social status²⁶,²⁷.

In the present study, we aimed at testing whether daily grooming and GGR rates across adult females are affected by the different value that females acquire based on their swelling phase, and whether such commodities are exchanged according to the Biological Market Theory (BMT).

**Prediction 1—basic assumption.** For a biological market to be established, it is necessary to first demonstrate that a desirable resource is present. If the maximum swelling phase increases the attractiveness of females and the motivation to engage in GGR³⁵,³⁶, we predict GGR rates to be higher when at least one of the females in the dyad is in the maximum swelling phase.

**Prediction 2.** If grooming is a commodity that is exchanged on a daily basis, we predict that grooming rates show daily fluctuations following the supply/demand law. Specifically, we expect female grooming rates to decrease as the number of available females in maximum swelling phase increases (*Prediction 2a*). Moreover, we expect dyadic grooming preference to vary on a daily basis as a function of actor’ and receiver’s swelling state (*Prediction 2b*).

**Prediction 3.** If grooming and GGR are commodities that are exchanged according to the BMT, we expect the daily frequency of GGR to fluctuate as a function of the number of females showing maximum swelling (supply/demand law) (*Prediction 3a*). Moreover, we expect GGR to be influenced by both the daily grooming rate (commodity) and the swelling phase (temporary value). Specifically, we predict the positive effect of grooming on the frequency of GGR to be maximum when the actor is in the minimum and the receiver is in the maximum swelling phase (*Prediction 3b*).

## Methods

**Ethic statement.** This study was approved by Universities of Pisa and Parma (Animal Care and Use Board). Since the study was purely observational with no human-animal interaction and/or manipulation whatsoever, the committee waived the need for an ethical permit.

**Study sites and subjects.** The study was part of a long-term research project on bonobo social behavior. Data were collected at the Apenheul Primate Park (Apeldoorn, The Netherlands) and at La Vallée des Singes (Romagne, France) where bonobos lived in inter-connected indoor enclosures with free access to an outdoor naturalistic forested island (Supplementary Table S1, Supplementary Information: Housing and living condition). In Apenheul, observations were carried out between August and November 2009 by two observers (one
being ED) for 6 days a week, seven hours per day. The group was composed by 5 adult females, 2 adult males, and 3 immature subjects (see Supplementary Information).

In La Vallée des Singes, observations were carried out between June and August 2012 by two observers (one being ED), 6 days a week from 9:00 to 18:00. The group was composed by 5 adult females, 4 adult male, and 1 immature subject (see Supplementary Information).

**Data collection.** In Apenheul, behavioural data were collected by using a voice-recorder and were then transcribed on database sheets. The techniques used were Focal Animal Sampling (30-min/Focal for a total of 25 h per individual) and All Occurrences Sampling (a total of 502 h)³⁹.

In La Vallée des Singes, behavioural data were collected by using camcoders and by applying the Focal Group Sampling method³⁹ for a total of 166 h of videos. The Focal Group Sampling method consists in randomly selecting a sub-group of subjects and recording their behaviours. In our case, we applied this methodology to collect the videos. Videos were then analysed and transcribed on database sheets. The techniques used to analyse the videos were: Focal Animal Sampling (20 h/individual) and All Occurrences Sampling (200 h)³⁹.

For both groups, allo-grooming sessions were collected through the Focal Animal Sampling method and were used to establish grooming preference. Sexual and agonistic interactions were collected through the All-Occurrences Sampling method and the latter were used to determine group hierarchy (see below for hierarchy estimation). Kinship through the maternal line was known.

Data were collected and analysed according to a blind coding protocol in which the observers were not aware of the hypotheses and predictions that would have been tested. Before systematic data collection and video analyses, the observers and coders underwent a training phase during which they followed the same focal animal simultaneously and then compared data. The training phase lasted until the Cohen's Kappa values measured for all the behaviours defined in Supplementary Table 2 reached 0.85⁴⁰.

**Definition of the sexual swelling phases.** Daily observations of perineal swelling and menstruation were carried out by the bonobo keepers. Keepers from Apenheul Primate Park recorded daily appearance of sexual swelling stage, according to the following rating system: 1 (maximum detumescence for noncycling females); 2 (minimum swelling phase for cycling females); 3 and 4 (intermediate swelling phase for cycling females); and 5 (maximum swelling phase for cycling females). Using the same parameters, keepers from La Vallée des Singes recorded daily appearance of sexual swellings stage according to the following rating system: 1 (maximum detumescence for noncycling females); 2 (minimum swelling phase for cycling females); 3 (intermediate swelling phase for cycling females); and 4 (maximum swelling phase for cycling females). For data unification, we standardized sexual swelling phases in two categories: minimum (phase 1+2) and maximum swelling phase (Apenheul: phase 5; La Vallée des Singes: phase 4). For each adult female showing monthly variations in the swelling cycle (Apenheul group: N = 3; La Vallée group: N = 4—See Supplementary Table 1), the sexual swelling status was recorded by the zoo keepers every day, therefore we analyzed the dyadic grooming preference on a daily basis.

**Calculation of daily grooming preference.** The daily grooming preference was calculated by the daily grooming duration (A grooms B) divided by the sum of the daily directional grooming duration performed by the same individual towards all other subjects including B [A grooms (B+C+D)]. The daily grooming preference therefore expresses the proportion of time that one individual spends grooming another one, weighted on the total amount of grooming provided by the same subject on a daily basis. Our study focused on individuals offering grooming, so polyadic grooming sessions were split into dyadic grooming sessions when the focal individual was receiving grooming from multiple partners. It is nevertheless worth noting that polyadic grooming sessions are relatively infrequent in bonobos and very rare among bonobo females⁴¹. Moreover, there is a positive correlation between the dyads involved in dyadic and in polyadic session⁴¹.

**Hierarchy.** We evaluated hierarchical relationships on the basis of dyadic and decided conflicts that occurred over the whole period of observation by using a winner/loser matrix to assess the rank according to the Normalized David's Scores (NDS⁴²). NDS were calculated on the basis of a dyadic dominance index (Dij) in which the observed proportion of wins (Pij) is corrected for the chance occurrence of the observed outcome. The chance occurrence of the observed outcome was calculated on the basis of a binomial distribution with each animal having an equal chance of winning or losing in every agonistic encounter⁴². The correction is necessary when the interaction numbers greatly differed across dyads.

**Statistics.** To compare GGR frequencies of dyads involving females in the minimum swelling phase to dyads involving at least one female in the maximum swelling phase, we ran a paired randomization test at the dyadic level with a number of 10,000 shuffles (allowing us to achieve an accuracy of 0.001 of the probability values) to avoid errors due to non-independence of the data⁴³. The test allowed us to understand whether the maximum swelling condition is attractive for other females (*Prediction 1—Basic assumption*). We estimated the individual frequency of GGR per day for each swelling phase (e.g., number of GGR<sub>MaxSwelling</sub>/number of days in maximum swelling). We used the software Resampling Procedures 1.3 by David C. Howell (freeware), that provides a t-value (difference between the means of the samples, standardized by the standard error) when comparing two dependent groups, with the probability to obtain such values under the null hypothesis.

We ran all Generalized Linear Mixed Models (GLMM) models in R (version 4.0.1)⁴⁴ using the RStudio interface (version 1.3.959⁴⁵). Prior to fitting the model, all the covariates were z-transformed to a mean of 0 and a standard deviation of 1 to achieve more easily interpretable estimates and facilitate model convergence⁴⁶. To
rule out collinearity we determined Variance Inflation Factors (VIFs) using the function *vif* of the R package *car* (version 3.0-3^47^) applied to a standard linear model excluding the random effects and the interactions. Both Model 1 and Model 2 presented below have been tested with the inclusion of female parity as factor shaping grooming preference and GGR. However, due to very high collinearity with actor and receiver's age (VIF *parity*.act = 64.18; VIF *parity.rec* = 76.84) we excluded them from all the analyses. After removal of parity, all the VIF values were considered acceptable (maximum VIFs across models = 3.31). Model stability was assessed using a function that excluded the levels of the random effects one at a time^48^. For all models, we performed a full vs control model comparison, with the control model containing only random effects, random slopes and control predictors to test whether the inclusion of explanatory predictors improved model fit. Tests of the individual predictors were derived using likelihood ratio tests^49^ (R function *drop1* with argument 'test' set to "Chisq"). We conducted post-hoc comparisons using Tukey tests, and report only significant tests.

We ran Model 1 to test *Prediction 2a* (effect of the daily number of available females showing maximum swelling on grooming preference) and *Prediction 2b* (effect of actor' and receiver's swelling state on daily dyadic grooming preference). Model 1 investigated whether daily grooming preferences (as dependent variable) varied depending on the following predictors: the daily number of females in the maximum swelling condition, the sexual swelling condition (binomial) of actor and receiver as well as their two-way interaction, and the group (binomial). *Parity* was not included due to high collinearity with *age* (see above). Considering the response range 0–1, we modelled it with a Beta error distribution and logit link function after data transformation which excluded exact values of 0 and 1^50^. The number of females with the maximum sexual swelling included all the studied females that a given actor could access a given day. We included as control predictors the two-way interaction of actor and receiver's age (numeric) and the two-way interaction of actor and receiver's rank (continuous variable). We included the two-way interaction of actor and receiver's rank to control for a possible influence of rank similarity/dissimilarity on the daily grooming preference. We included the two-way interaction of actor and receiver's age to control for a possible influence of age similarity/dissimilarity. We removed non-significant two-way interactions and tested their main predictors. Identity of actor, receiver and dyad were included as random effects. We excluded random slopes due to convergence issues. Model 1 included 94 cases of daily grooming preferences involving 16 different dyads. In a given dyad AB, the grooming preference of subject A towards B is not equal to the preference of B towards A, allowing our model to estimate the effect of AB controlling for BA, and vice versa.

We ran Model 2 to test *Prediction 3a* (effect of the daily number of available females showing maximum swelling on daily GGR frequency) and *Prediction 3b* (effect of daily grooming preference on the daily frequency of GGR depending on actor' and receiver's swelling phase). Model 2 investigated the effect of the number of swollen females, daily grooming preference (continuous variable) and sexual swelling (binomial) on the daily occurrence of GGR (numeric), by considering the directionality of the interaction. We modelled the response with a Gamma error distribution and log link function. To calculate the occurrence of GGR, we analyzed the successful sexual invitations (i.e., those leading to a GGR) performed by the actor towards the receiver. Gestures, body postures and facial expressions used as sexual invitations are listed and described in Supplementary Table S2. All the GGR events not preceded by any clear sexual invitation were halved and assigned to both the actor and the receiver. We included as test predictors the following variables: the daily number of females in the maximum swelling condition, the three-way interaction of daily grooming preference and sexual swelling condition of both the actor and receiver, and the group. *Parity* was not included due to high collinearity with *age* (see above). We included as control predictors the two-way interaction of actor and receiver's age and the two-way interaction of actor and receiver's rank. We removed the not significant two-way interaction *age<sub>act</sub>* *age<sub>rec</sub>* and repeated the model to assess the effect of main age predictors. Identity of actor, receiver and dyad were included as random effects and the *number of swollen females* and *grooming preference* as random slopes. Because Model 2 did not converge, we excluded from the random slopes *grooming preference* and re-ran Model 2 successfully. Model 2 included 74 cases involving 13 dyads.

## Results

**Genito-genital rubbing.** Females engaged significantly less in GGR when they were both in the minimum swelling phase (mean ± SD: 0.32 ± 0.60), compared to when at least one female of the dyad was in the maximum swelling phase (mean ± SD: 0.61 ± 0.69) (Randomization paired t<sub>9</sub> = -2.275, *P* = 0.027) (*Prediction 1 supported*).

**Model 1—daily grooming preference.** Overall, the predictors *number of swollen females*, the two-way interaction *swelling<sub>act</sub>* \* *swelling<sub>rec</sub>* and *group*, had an effect on the daily grooming preference, after controlling for rank and age effect (full vs control model comparison Model 1: $\chi^2$ = 11.92, df = 5, *p* = 0.035, pseudo-R<sup>2</sup> = 0.32). The model did not include the non-significant interactions of *age<sub>act</sub>* \* *age<sub>rec</sub>* (*p* = 0.746) and *rank<sub>act</sub>* \* *rank<sub>rec</sub>* (*p* = 0.726) and assesses the effect of age and rank as main control predictors. The daily grooming preference was significantly affected by the number of swollen females available the same day: the higher the number of swollen females, the lower the daily grooming preference (*Prediction 2a supported*) (Table 1, Fig. 1). An increase of 1 in the number of females in the maximum swelling phase decreased the probability of grooming from 88 to 80%. Moreover, the two-way interaction *swelling<sub>act</sub>* \* *swelling<sub>rec</sub>* had a significant effect on the daily grooming preference and increased the likelihood of grooming of 2.27 times (+69%, Table 1). After removal of non-significant interactions, all main control predictors were not significant. We ran post-hoc comparisons to estimate the differences in the effect of actor' and receiver's swelling condition on daily grooming preferences (Fig. 2). Specifically, when both the actor and the receiver were in the maximum swelling condition the dyadic grooming preference was higher compared to when the actor was in the minimum and the receiver was in the maximum swelling condition (*p* = 0.009 Tukey adjusted, Fig. 2) (*Prediction 2b supported*).
<table><caption>**Table 1.** Variables in Model 1 explaining the daily grooming preference while controlling for rank and age of both actor and receiver. LTR = likelihood-ratio test. Control variables are in italics and explanatory variables with significant *P* values are in bold and marked with *<sup>a</sup>Swelling condition "min" as reference category.*

<sup>b</sup>Group1 as reference category. <sup>c</sup>z-transformed to a mean of 0 and a standard deviation of 1.</caption><thead><tr><th>Fixed variables</th><th>Estimates</th><th>SE</th><th>Z-value</th><th>2.5% C.I</th><th>97.5% C.I</th><th>LRT</th><th>df</th><th>P</th></tr></thead><tbody><tr><td>Intercept</td><td>1.95</td><td>0.50</td><td>3.93</td><td>0.98</td><td>2.92</td><td>-</td><td>-</td><td>-</td></tr><tr><td>N. of swollen females</td><td>-0.52</td><td>0.26</td><td>-2.03</td><td>-1.02</td><td>-0.02</td><td>4.41</td><td>1</td><td>0.036*</td></tr><tr><td><sup>a</sup>Swelling act<sub>(max)</sub> * Swelling rec<sub>(max)</sub></td><td>1.13</td><td>0.57</td><td>1.98</td><td>0.01</td><td>2.25</td><td>3.94</td><td>3</td><td>0.047*</td></tr><tr><td><sup>a</sup>Swelling act<sub>(max)</sub></td><td>0.22</td><td>0.52</td><td>0.43</td><td>-0.80</td><td>1.24</td><td>-</td><td>-</td><td>-</td></tr><tr><td><sup>a</sup>Swelling rec<sub>(max)</sub></td><td>-0.56</td><td>0.60</td><td>-0.94</td><td>-1.74</td><td>0.62</td><td>-</td><td>-</td><td>-</td></tr><tr><td><sup>b</sup>Group2</td><td>-0.13</td><td>0.51</td><td>-0.25</td><td>-1.13</td><td>0.88</td><td>0.06</td><td>1</td><td>0.805</td></tr><tr><td><sup>c</sup>Rank act</td><td>0.06</td><td>0.22</td><td>0.29</td><td>-0.36</td><td>0.49</td><td>0.08</td><td>1</td><td>0.773</td></tr><tr><td><sup>c</sup>Rank rec</td><td>0.47</td><td>0.28</td><td>1.69</td><td>-0.07</td><td>1.02</td><td>3.11</td><td>1</td><td>0.078</td></tr><tr><td><sup>c</sup>Age act</td><td>-0.11</td><td>0.22</td><td>-0.51</td><td>-0.54</td><td>0.31</td><td>0.24</td><td>1</td><td>0.622</td></tr><tr><td><sup>c</sup>Age rec</td><td>-0.23</td><td>0.24</td><td>-0.97</td><td>-0.70</td><td>0.23</td><td>0.98</td><td>1</td><td>0.322</td></tr></tbody></table>

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Effect of the number of females in maximum swelling phase on daily grooming preference (red solid line) at the average of all other predictors. Gray band indicates 95% Confidence Intervals. Data points are plotted using jitter function at 0.15 (width) to avoid complete overlapping and improve graph interpretability.</figcaption></figure>

**Model 2 – Sex, swelling and grooming.** Overall, the number of swollen females and the three-way interaction *Grooming*\* *Swelling*<sub>act</sub>\* *Swelling*<sub>rec</sub> had an effect on daily GGR distribution after controlling for rank and age (full vs control model comparison of Model 2: $\chi^2 = 22.75$, df = 10, *p* = 0.011, conditional R<sup>2</sup> = 0.40). The model did not include the non-significant interaction of *age*<sub>act</sub>\* *age*<sub>rec</sub> (*p* = 0.707) and assesses the effect of age as main predictors. The daily GGR occurrence was significantly affected by the number of females in the maximum swelling phase: the higher the number of swollen females, the lower the GGR occurrence (–29.83% probability of GGR with an increase of 1 female in maximum swelling phase; *p* = 0.002, Table 2) (*Prediction 3a supported*). Moreover, the daily GGR occurrence was significantly affected by the three-way interaction of daily grooming preference and sexual swelling phase of both the actor and the receiver (*p* = 0.011; Table 2, Fig. 3). Actors in maximum swelling phase (MAX) had 1.02 times higher (+2.94%) odds of engage in GGR than actors in minimum swelling phase (MIN), while receivers in MAX had 1.90 times higher (+90.33%) odds than receiver in MIN. Additionally, receivers in MAX had 1.84 times higher (+84.89%) odds of engage in GGR when compared with actors in the same condition, highlighting the role of receivers in accepting sexual contacts. Interestingly, an increase of 1 SD of grooming preference when actor and receiver where both in the minimum swelling phase (MIN-MIN) was associated with an effect size of 0.54 (–45.70%), suggesting that grooming is not exchanged with GGR when both females are in the minimum swelling phase. Consistently with the Biological Market Theory, the grooming preference of actor in minimum and receiver in maximum swelling phase (MIN-MAX) had an effect 1.32 times higher (+32.01%) when compared with MAX-MIN, and 4.43 times higher (+342.61%) when compared with MIN-MIN.

Post-hoc analyses confirmed the differences in the effect of grooming preference on daily GGR occurrence according to different sexual swelling conditions. Specifically, when the actor was in the minimum swelling condition and the receiver was in the maximum swelling condition, grooming preference increased daily GGR
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Predicted effect of the interaction of actor and receiver's swelling state on the dyadic grooming preference. Dots with borders and whiskers indicate estimated mean and standard errors of categorical responses, while colored lines indicate variation in receiver's effect according to actor's swelling. The figure is plotted with all other predictors at their average value. Dots without borders indicate raw data plotted with *jitter* function at 0.1 (width). Parentheses indicate significant post-hoc comparison.</figcaption></figure>

<table><caption>**Table 2.** Variables in Model 2 explaining the occurrence of daily GGR while controlling for rank and age of both actor and receiver. *LTR* likelihood-ratio test. Control variables are in italics, variables with significant *P* values are in bold and marked with *. *Grooming* = daily dyadic grooming preference. *Swelling*, *act,rec* = sexual swelling condition of actor or receiver. <sup>a</sup>Dummy coded with “min” swelling condition as reference category.</caption><thead><tr><th>Fixed variables</th><th>Estimates</th><th>SE</th><th>Z-value</th><th>2.5% C.I</th><th>97.5% C.I</th><th>LRT</th><th>df</th><th>P</th></tr></thead><tbody><tr><td>Intercept</td><td>1.51</td><td>0.07</td><td>21.38</td><td>1.37</td><td>1.65</td><td>-</td><td>-</td><td>-</td></tr><tr><td><strong>N. of swollen females</strong></td><td><strong>-0.35</strong></td><td><strong>0.04</strong></td><td><strong>-9.22</strong></td><td><strong>-0.43</strong></td><td><strong>-0.28</strong></td><td><strong>9.20</strong></td><td><strong>1</strong></td><td><strong>0.002</strong></td></tr><tr><td><sup>a,b</sup>Grooming * Swelling<sub>act(max)</sub> * Swelling<sub>rec(max)</sub></td><td><strong>-0.91</strong></td><td><strong>0.07</strong></td><td><strong>-13.78</strong></td><td><strong>-1.04</strong></td><td><strong>-0.78</strong></td><td><strong>6.41</strong></td><td><strong>1</strong></td><td><strong>0.011</strong></td></tr><tr><td><sup>a,b</sup>Grooming * Swelling<sub>act(max)</sub></td><td>0.60</td><td>0.06</td><td>9.36</td><td>0.47</td><td>0.72</td><td>-</td><td>-</td><td>-</td></tr><tr><td><sup>a,b</sup>Grooming * Swelling<sub>rec(max)</sub></td><td>0.88</td><td>0.06</td><td>13.94</td><td>0.75</td><td>1.00</td><td>-</td><td>-</td><td>-</td></tr><tr><td><sup>a,b</sup>Swelling<sub>act(max)</sub> * Swelling<sub>rec(max)</sub></td><td>-0.40</td><td>0.07</td><td>-5.77</td><td>-0.53</td><td>-0.26</td><td>-</td><td>-</td><td>-</td></tr><tr><td><sup>b</sup>Grooming</td><td>-0.61</td><td>0.06</td><td>-10.48</td><td>-0.72</td><td>-0.50</td><td>-</td><td>-</td><td>-</td></tr><tr><td><sup>a</sup>Swelling<sub>act(max)</sub></td><td>0.03</td><td>0.07</td><td>0.43</td><td>-0.10</td><td>0.16</td><td>-</td><td>-</td><td>-</td></tr><tr><td><sup>a</sup>Swelling<sub>rec(max)</sub></td><td>0.64</td><td>0.07</td><td>9.34</td><td>0.51</td><td>0.78</td><td>-</td><td>-</td><td>-</td></tr><tr><td><sup>b</sup>Group2</td><td>0.20</td><td>0.07</td><td>2.79</td><td>0.06</td><td>0.34</td><td>0.33</td><td>1</td><td>0.563</td></tr><tr><td><sup>c</sup>Rank act * rank rec</td><td>0.37</td><td>0.05</td><td>6.84</td><td>0.26</td><td>0.47</td><td><strong>7.03</strong></td><td><strong>1</strong></td><td><strong>0.008</strong></td></tr><tr><td><sup>c</sup>Rank act</td><td>-0.37</td><td>0.05</td><td>-6.99</td><td>-0.47</td><td>-0.26</td><td>-</td><td>-</td><td>-</td></tr><tr><td><sup>c</sup>Rank rec</td><td>-0.51</td><td>0.06</td><td>-9.00</td><td>-0.62</td><td>-0.40</td><td>-</td><td>-</td><td>-</td></tr><tr><td><sup>c</sup>Age act</td><td><strong>0.48</strong></td><td><strong>0.05</strong></td><td><strong>8.88</strong></td><td><strong>0.38</strong></td><td><strong>0.59</strong></td><td><strong>10.77</strong></td><td><strong>1</strong></td><td><strong>0.001</strong></td></tr><tr><td><sup>c</sup>Age rec</td><td>0.21</td><td>0.05</td><td>3.94</td><td>0.11</td><td>0.32</td><td>2.15</td><td>1</td><td>0.142</td></tr></tbody></table>

<sup>b</sup>Dummy coded with Group1 as reference category. <sup>c</sup>-transformed to a mean of 0 and a standard deviation of 1.
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Effect of the three-way interaction *grooming*\* *swelling*<sub>act</sub> \* *swelling*<sub>rec</sub>. (a) Indicates actor in the minimum swelling phase when interacting with a receiver in the minimum swelling phase (in blue) or with a receiver in the maximum swelling phase (in red). (b) Indicates actor in the maximum swelling phase when interacting with a receiver in the minimum swelling phase (in yellow) or with a receiver in the maximum swelling phase (in green). Shaded bands indicate 95% Confidence Intervals. Darker points represent overlapping of several data points.</figcaption></figure>

occurrence when compared with both actor and receiver in the minimum swelling condition (MIN-MAX vs MIN-MIN: *p* < 0.001, Fig. 3). After removing the non-significant two-way interaction *age*<sub>act</sub> \* *age*<sub>rec</sub>, only the main control predictor *age*<sub>act</sub> was significant (*p* = 0.001, Table 2, Supplementary Figure 1). Moreover, the two-way interaction of the control predictors *rank*<sub>act</sub> \* *rank*<sub>rec</sub> was significant (*p* = 0.008, Table 2, Supplementary Figure S2), with a positive effect on the daily occurrence of GGR.

Discussion

Taken together, our results confirm the key role of sexual swelling in bonobos and provide new insights on the exchange of grooming with sex on a daily basis. More particularly, our study shows that grooming and sex represent valuable exchange commodities for bonobo females and that their exchange fluctuates depending on females' sexual swelling status. These daily fluctuations can be explained by the supply/demand law¹. In line with other studies on sexual swelling and socio-sexual behaviour in bonobos³⁵,³⁶, our results confirmed that females engaged significantly more in Genito-Genital Rubbing (GGR) when at least one female of the dyad was in the maximum swelling phase, compared to when both of them were in the minimum swelling phase (*Prediction 1*). When in the maximum swelling phase, females are more attractive for both males and females and tend to be invited more frequently to engage in sexual interactions than females in the minimum swelling phase²⁶. This is in line with the hypothesis proposing that the peculiar features of bonobo females' sexual swelling, and in particular its extremely long duration and presence during anovulatory periods, might have been selected to increase female-female sexual interactions³⁵. Socio-sexuality in bonobo females help establishing and maintaining strong social bonds which, in turn, allow them to form coalitions that increase their social status and centrality³⁰.

According to the Biological Market Theory derived from the Market Theory applied in economics, the exchange rate of commodities depends on the supply/demand law, and when a certain commodity becomes common, its perceived value decreases. Consistent with the supply/demand law, the daily grooming preference was significantly affected by the number of females showing the sexual swelling at the maximum phase during the same given day: the higher the number of swollen females, the lower the daily grooming preference (*Prediction 2a*, Fig. 1). Interestingly, when both the actor and the receiver were in the minimum swelling phase, and there were no females available in the maximum swelling phase, the grooming preference was very high. Via Model 1 we considered the directional grooming preference and we included the preference of both dyads AB and BA. Therefore, our results suggest that females exchanged grooming for grooming when both in the minimum swelling phase. Moreover, the significant effect of the interaction *swelling*<sub>act</sub> \* *swelling*<sub>rec</sub> showed higher grooming preference when both the actor and the receiver were in the maximum swelling phase, compared to when the actor was in the minimum and the receiver was in the maximum swelling phase. These results confirm the higher attractiveness of females showing the maximum swelling phase (*Prediction 1*, *2b*), suggest that females tend to exchange grooming for grooming when their value is “low” (MIN) and that high-value females (MAX) prefer “high-value” grooming partners as well.
Our results on the effect of the number of swollen females on grooming preference are in line with a study on male-female grooming in captive chimpanzees¹¹ showing that females in the maximum swelling phase receive more grooming. Although our investigation focussed on female-female sexual interactions, our results are in line with those of chimpanzees given that the highest daily preference was achieved when both actor and receiver were in the maximum swelling phase. The results coming from the two *Pan* species are comparable considering the high attractiveness of the maximum sexual swelling in both chimpanzees and bonobos although in chimpanzees the sexual swelling is attractive for males only, while in bonobos this signal is attractive for both sexes.

In many primate groups, higher-ranking subjects typically receive more grooming than lower-ranking ones⁵¹⁻⁵³ and most of the grooming occurs between subjects with similar rank⁵⁴. Among bonobos, the distribution of grooming is related to several variables (e.g., age, rank, and sex) and there are some indications that high-ranking and older females are preferred grooming partners⁵⁵. Subjects within the same group may therefore compete for grooming access to high-ranking groupmates, since those can be the best coalition partners⁵²,⁵⁴. Alternatively, a grooming distribution related to rank can be explained by females' preference in grooming subjects of similar social class (e.g., rank and age) as an adaptive strategy derived from greater compatibility between subjects⁵⁶. These mechanisms are beyond the aim of our investigation but can still drive the variation of grooming preferences and mask the effect of the daily market. For this reason, we included in our models rank and age of both actor and receiver as control predictors. However, our results on the daily grooming preference showed no effect of actor and receiver's rank and age, suggesting that the factors shaping grooming preference among bonobo females on a daily basis may not coincide with those affecting such preference on the long term. Moreover, it has been showed that the steepness of the hierarchy varies between different bonobo groups in captivity²⁰,²⁵. Particularly, these studies demonstrated that grooming was preferentially directed towards dominant individuals in those groups where the hierarchy was very steep, whereas grooming was more reciprocal in those groups showing a mild hierarchy. It is worth noting that the steepness of our study groups was extremely low (Apenheul: 0.346; La Vallée: 0.305 unpublished data).

The daily occurrence of GGR was significantly affected by the number of females in the maximum swelling phase: the higher the number of swollen females, the lower the daily occurrence of GGR (*Prediction 3a*). Post-hoc analyses revealed no significant effect of grooming preference on the daily occurrence of GGR when the actor was in the maximum swelling condition and the receiver was in the maximum or minimum swelling condition (Fig. 3b). On the other hand, when an actor in the minimum swelling phase preferentially groomed a receiver in the maximum swelling phase, it significantly increased the probability of engaging in GGR compared to when the receiver was in the minimum swelling phase (*Prediction 3b*, Fig. 3a). This might indicate that: (1) when both partners are in the minimum swelling phase, none of them has the “economic power” to influence the exchange of goods and grooming is therefore not exchanged with sex (Fig. 3a; MIN-MIN), and that (2) lower-value females (MIN) need to provide more grooming to get the access to higher-value partners (Fig. 3a; MIN-MAX). However, this strategy seems to work only when is the receiver (of grooming/sexual invitations) who benefits from the mismatch, especially if we consider that an increase of grooming performed by a high-value actor does not increase the chances of engaging in sexual interactions: a high-value actor has a constant chance of engaging in sexual interactions despite receiver's value (Fig. 3b; MAX-MIN, MAX-MAX). The results from Model 1 and Model 2 are consistently stressing that grooming can be performed to balance the mismatch in females' value, and when the partners have equal value, grooming can be exchanged for grooming as well.

Females in the minimum swelling phase exhibit a short-term strategy of exchanging commodities shaped by the supply/demand law and receiver's value. Given that females in the maximum swelling phase are more sexually (and therefore socially) attractive, females in the minimum swelling phase must augment the commodity they provide to increase their possibility of engaging in sexual interactions with higher-value females. Furthermore, by strategically offering more grooming, lower-value females can even outcompete higher-value ones to gather “sexual favours” from the most attractive partners.

In conclusion, our study confirms the BMT in explaining the daily exchange of commodities and its fluctuations as a function of the temporary value of traders, and underlines the importance of a day-by-day approach to unveil the presence of a biological market when the value of traders rapidly changes.

## Data availability

The raw data supporting the conclusion of this article are provided as supporting material.

Received: 11 June 2021; Accepted: 15 September 2021

Published online: 29 September 2021

## References

1. Noë, R. & Hammerstein, P. Biological markets: supply and demand determine the effect of partner choice in cooperation, mutualism and mating. *Behav. Ecol. Sociobiol.* **35**, 1–11 (1994).

2. Lewis, R. J. Beyond dominance: the importance of leverage. *Q. R. Biol.* **77**, 149–164 (2002).

3. Hammerstein, P. & Noë, R. Biological trade and markets. *Philos. Trans. R. Soc. B.* **371**, 20150101 (2016).

4. Barclay, P. Strategies for cooperation in biological markets, especially for humans. *Evol. Hum. Behav.* **34**, 164–175 (2013).

5. Bshary, R. The cleaner fish market. In *Economics in Nature: Social Dilemmas, Mate Choice and Biological Markets* (eds Noë, R. et al.) 146–172 (Cambridge University Press, 2001).

6. Metz, M., Klump, G. M. & Friedl, T. W. Temporal changes in demand for and supply of nests in red bishops (*Euplectes orix*): dynamics of a biological market. *Behav. Ecol. Sociobiol.* **61**, 1369–1381 (2007).

7. Stopka, P. & Macdonald, D. W. The market effect in the wood mouse, *Apodemus sylvaticus*: selling information on reproductive status. *Ethology* **105**, 969–982 (1999).

8. Kutsukake, N. & Clutton-Brock, T. The number of subordinates moderates intrasexual competition among males in cooperatively breeding meerkats. *Proc. R. Soc. B* **275**, 209–216 (2008).
9. Norscia, I., Antonacci, D. & Palagi, E. Mating first, mating more: biological market fluctuation in a wild prosimian. *PLoS ONE* 4, e4679 (2009).

10. Barelli, C., Reichard, U. & Mundry, R. Is grooming used as a commodity in wild white-handed gibbons, *Hylobates lar?*. *Anim. Behav.* **82**, 801–809 (2011).

11. Koyama, N. F., Caws, C. & Aureli, F. Supply and demand predict male grooming of swollen females in captive chimpanzees, *Pan troglodytes*. *Anim. Behav.* **84**, 1419–1425 (2012).

12. Surbeck, M. & Hohmann, G. Social preferences influence the short-term exchange of social grooming among male bonobos. *Anim. Cogn.* **18**, 573–579 (2014).

13. Dunbar, R. I. M. *Reproductive Decisions: An Economic Analysis of Gelada Baboon Social Strategies* (Princeton University Press, 1984).

14. Schino, G. Grooming and agonistic support: A meta-analysis of primate reciprocal altruism. *Behav. Ecol.* **18**, 115–120 (2007).

15. Seyfarth, R. M. & Cheney, D. L. Grooming, alliances and reciprocal altruism in vervet monkeys. *Nature* **308**, 541–543 (1984).

16. Aureli, F., Preston, S. D. & de Waal, F. B. M. Heart rate responses to social interactions in free-moving rhesus macaques (*Macaca mulatta*): A pilot study. *J. Comp. Psychol.* **113**, 59–65 (1999).

17. Sanchez-Amaro, A. & Amici, F. Are primates out of the market?. *Anim. Behav.* **110**, 51–60 (2015).

18. Tomasello, M. & Call, J. *Primate Cognition* (Oxford University Press, 1997).

19. Gilby, I. C. Cooperation among non-kin: reciprocity, markets and mutualism. In *The Evolution of Primate Societies* (eds Mitani, J. C. *et al.*) 514–530 (The University of Chicago Press, 2013).

20. Stevens, J. M., Vervaecke, H., De Vries, H. & van Elsacker, L. The influence of the steepness of dominance hierarchies on reciprocity and interchange in captive groups of bonobos (*Pan paniscus*). *Behaviour* **142**, 941–960 (2005).

21. Kano, T. *The last ape* (Pygmy chimpanzee behavior and ecology. Stanford University Press, 1992).

22. Parish, A. R. Female relationships in bonobos (*Pan paniscus*) Evidence for bonding, cooperation, and female dominance in a male-philopatric species. *Hum. Nat.* **7**, 61–96 (1996).

23. Tokuyama, N. & Furuichi, T. Do friends help each other? Patterns of female coalition formation in wild bonobos at Wamba. *Anim. Behav.* **119**, 27–35 (2016).

24. Moscovice, L. R. *et al*. Stable and fluctuating social preferences and implications for cooperation among female bonobos at LuiKotale, Salonga National Park, DRC. *Am. J. Phys. Anthropol.* **163**, 158–172 (2017).

25. Stevens, J. M., Vervaecke, H., De Vries, H. & van Elsacker, L. Sex differences in the steepness of dominance hierarchies in captive bonobo groups. *Int. J. Primatol.* **28**(6), 1417–1430 (2007).

26. Furuichi, T. Female contributions to the peaceful nature of bonobo society. *Evol. Anthropol.* **20**, 131–142 (2011).

27. Palagi, E., Paoli, T. & Tarli, S. B. Short-term benefits of play behavior and conflict orevention in *Pan paniscus*. *Int. J. Primatol.* **27**, 1257–1270 (2006).

28. Palagi, E., Paoli, T. & Borgognini Tarli, S. M. Reconciliation and consolation in captive bonobos (*Pan paniscus*). *Am. J. Primatol.* **62**, 15–30 (2004).

29. Annicchiarico, G., Bertini, M., Cordoni, G. & Palagi, E. Look at me while having sex! Eye-to-eye contact affects homosexual behaviour in bonobo females. *Behaviour* **157**, 949–970 (2020).

30. Moscovice, L. R. *et al*. The cooperative sex: Sexual interactions among female bonobos are linked to increases in oxytocin, proximity and coalitions. *Hormones Behav.* **116**, 104581 (2019).

31. Clay, Z. & de Waal, F. B. M. Sex and strife: Post-conflict sexual contacts in bonobos. *Behaviour* **152**, 313–334 (2015).

32. Kuroda, S. Rocking gesture as communicative behaviour in the wild pygmy chimpanzee in Wamba, Central Zaire. *J. Ethol.* **2**, 127–137 (1984).

33. Palagi, E., Bertini, M., Annicchiarico, G. & Cordoni, G. Mirror replication of sexual facial expressions increases the success of sexual contacts in bonobos. *Sci. Rep.* **10**, 18979 (2020).

34. de Waal, F. B. M. & Preston, S. D. Mammalian empathy: Behavioural manifestations and neural basis. *Nat. Rev. Neurosci.* **18**, 498–509 (2017).

35. Ryu, H., Hill, D. A. & Furuichi, T. Prolonged maximal sexual swelling in wild bonobos facilitates affiliative interactions between females. *Behaviour* **152**, 285–311 (2014).

36. Paoli, T., Palagi, E., Tacconi, G. & Borgognini Tarli, S. Perineal swelling, intermenstrual cycle, and female sexual behavior in bonobos (*Pan paniscus*). *Am. J. Primatol.* **68**, 333–347 (2006).

37. Douglas, P. H., Hohmann, G., Murtagh, R., Thiessen-Bock, R. & Deschner, T. Mixed messages: wild female bonobos show high variability in the timing of ovulation in relation to sexual swelling patterns. *BMC Ecol. Evol.* **16**, 1–7 (2016).

38. Demuru, E., Pellegrino, F., Dediu, D. & Levréro, F. Foraging postures are a potential communicative signal in female bonobos. *Sci. Rep.* **10**, 1–7 (2020).

39. Altmann, J. Observational study of behaviour sampling methods. *Behaviour* **49**, 227–265 (1974).

40. Kaufman, A. B. & Rosenthal, R. Can you believe my eyes? The importance of interobserver reliability statistics in observations of animal behaviour. *Anim. Behav.* **78**, 1487–1491 (2009).

41. Girard-Buttoz, C. *et al*. Variable use of polyadic grooming and its effect on access to social partners in wild chimpanzees and bonobos. *Anim. Behav.* **168**, 211–224 (2020).

42. de Vries, H., Stevens, J. M. G. & Vervaecke, H. Measuring and testing the steepness of dominance hierarchies. *Anim. Behav.* **71**, 585–592 (2006).

43. Manly, B. F. J. *Randomization, Bootstrap and Monte Carlo Methods in Biology* (Chapman & Hall, 1997).

44. R Development Core Team. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing (2020). http://www.R-project.org/.

45. RStudio Team. RStudio: Integrated Development for R. RStudio, PBC, Boston, MA (2020). http://www.rstudio.com/.

46. Schielzeth, H. Simple means to improve the interpretability of regression coefficients. *Methods Ecol. Evol.* **1**, 103–113 (2010).

47. Fox, J. & Weisberg, S. *An R Companion to Applied Regression* 3rd edn. (Sage, 2019).

48. Nieuwenhuis, R., Te-Grotenhuis, H. F. & Pelzer, B. ME: tools for detecting influential data in mixed effects models. *R J.* **4**, 38 (2012).

49. Barr, D. J., Levy, R., Scheepers, C. & Tily, H. J. Random effects structure for confirmatory hypothesis testing: Keep it maximal. *J. Mem. Lang.* **68**, 255–278 (2013).

50. Smithson, M. & Verkuilen, J. A better lemon squeezer? Maximum-likelihood regression with beta-distributed dependent variables. *Psychol. Methods* **11**, 54–71 (2006).

51. Tiddi, B., Aureli, F. & Schino, G. Grooming up the hierarchy: The exchange of grooming and rank-related benefits in a new world primate. *PLoS ONE* 7, e36641 (2012).

52. Seyfarth, R. M. A model of social grooming among adult female monkeys. *J. Theor. Biol.* **65**, 671–698 (1977).

53. Schino, G. Grooming, competition and social rank among female primates: A meta-analysis. *Anim. Behav.* **62**, 265–271 (2001).

54. Seyfarth, R. M. The distribution of grooming and related behaviours among adult female vervet monkeys. *Anim. Behav.* **28**, 798–813 (1980).

55. Van Elsacker, L., Vervaecke, H. & De Vries, H. The pivotal role of rank in grooming and support behavior in a captive group of bonobos (*Pan paniscus*). *Behaviour* **137**, 1463–1485 (2000).

56. de Waal, F. B. M. & Luttrell, L. M. The similarity principle underlying social bonding among female rhesus monkeys. *Folia Primatol.* **46**, 215–234 (1986).
## Acknowledgements

The authors wish to thank the bonobo keepers of Apenheul and La Vallée des Singes, and particularly the curators Frank Rietkerk and Jean-Pascal Guéry who permitted this study. We thank Francesca Coppola, Antonio Santucci and Rhea Arini for field assistance, Roger Mundry for statistical tools, and Federica Amici for her insightful comments and suggestions which improved our manuscript. Finally, we thank the Journal Club group of the Department of Behavioral Ecology of the University of Goettingen for stimulating discussions on the Biological Market Theory. The 2012 data collection was funded by the University of Parma (Italy).

## Author contributions

E.P. conceived the study. E.D. collected data and analyzed video data. All authors conceived the statistical analyses. S.A. performed statistical analyses. All authors provided critical feedback and helped shape the research, analyses and manuscript.

## Funding

Open Access funding enabled and organized by Projekt DEAL.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-021-98894-w.

Correspondence and requests for materials should be addressed to S.A. or E.D.

Reprints and permissions information is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_6.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2021, corrected publication 2021