Introduction

Energetic management, i.e. balancing energy intake and expenditure to maintain a positive energy balance, is a key determinant of animal reproductive success and survival (Hobbs 1989; Bosaeus et al. 2002; Stubbs and Tolkamp 2006). Most wild animals are subject to seasonal variation in food availability and reproductive opportunities which lead to variation in their energy balance (Tsuji et al. 2008; Emery Thompson et al. 2009). In non-human primates (hereafter “primates”), individuals face an additional challenge which is the high energy demand of their exceptionally large brain (Charnov and Berrigan 1993; Isler and van Schaik 2009). Thus, group-living primates have evolved specific behavioural and social adaptions to fulfil their energetic needs and adjust them to constraints such as seasonal variation in food abundance and mating competition (Zhang and Wang 2007; Tsuji et al. 2008; Emery Thompson et al. 2009; Wessling et al. 2018a).

Primates exhibit a wide variety of social organization ranging from solitary to large cohesive multi-male multi-female groups (Mitani et al. 2012). To better understand this variety and the social dynamics within groups, the socioecological model, linking food distribution and abundance to primate social dynamics, was developed in the 1980s (Wrangham 1980; van Schaik 1989). This model states that the distribution of food in space and time primarily affects female aggregation and feeding competition among females whereas males are primarily impacted by the spatio-temporal distribution of females and male-male competition for access to mates (Steenbeek and Sterck 1997; Kappeler 2000; Koenig 2002). The model does not claim that mating competition does not impact female sociality or that males do not compete over food but highlights a sex difference in the magnitude of both effects for each sex. In turn, female energy balance should be more sensitive to variation in ecological factors than males whereas male energy balance should be primarily impacted by the degree of mating competition in the group. Socioecological theories further hypothesise that female gregariousness evolved as a strategy to respond to between-group feeding competition and predation pressure (Wrangham 1980; van Schaik 1989). However, females are limited in their energetic acquisition by within-group feeding competition (van Schaik 1989; Barton and Whiten 1993; Saito 1996; Wittig and Boesch 2003; Murray et al. 2006). This effect is expected to be unequal across conspecifics living in the same community since the high-ranking ones can easily monopolize access to food resources (Greengrass 2005; Pusey et al. 2005; Kahlenberg et al. 2008; Nurmi et al. 2018). Accordingly, high-ranking females are expected to acquire more energy and have higher energy intake, producing a more positive energy balance than low-ranking females (Ellis 1995; Sterck et al. 1997). In contrast, the effect of dominance rank on energy intake in males is expected to be lower than in females, and rank effects on male energy balance are expected mainly in periods of high mating competition where high-ranking males compete to access fertile females (Muller and Wrangham 2004).

Nevertheless, there is a lack of studies directly assessing the impact of feeding competition and mating competition on energetic condition in both sexes within the same study and using a single analysis. More generally, and because of the foundational assumption of the socioecological theories, studies analyzing the impact of food availability on energetic condition often focus on females (Emery Thompson et al. 2012), and studies analyzing the impact of mating competition (measured by the presence of fertile females, the rate of male-male aggression and/or the degree of female monopolization) on energetic conditions focus on males (Emery Thompson et al. 2009; Georgiev 2012; Girard-Buttoz et al. 2014a).

To fill this gap, we quantify how socioecological factors drive male and female western chimpanzee’s (Pan troglodytes verus) energetic condition to better understand the energetic management strategies of this subspecies, characterized by a lower occurrence of female coercion by males than eastern chimpanzees (Pan troglodytes schweinfurthii) (Stumpf and Boesch 2005; Boesch et al. 2008). In chimpanzees, studies directly linking variation in food availability and mating competition to the energetic condition of males and females are limited, and these studies led to conflicting results. In western chimpanzees, the presence of oestrus females overruled ecological factors, and food availability did not affect grouping patterns when these females were present (Anderson et al. 2002). In eastern chimpanzees, males in large parties reduced their feeding time in the presence of oestrus females, and there was no effect of mean daily party size on male feeding time (Georgiev et al. 2014).

Furthermore, these studies rely mostly on behavioural measures of feeding time or food intake which can lead to imprecise estimation of energy gain and provide no information on energy balance (Emery Thompson 2017). To circumvent this limitation, researchers have developed and validated the use of physiological markers of wild primate energetic conditions such as c-peptide, a peptide produced on an equimolar ratio with insulin, a hormone involved in compartmentalizing and storing of energy (Steiner and Oyer 1967; Rubenstein et al. 1969; Melani et al. 1970). Urinary c-peptide (UCP) has been validated as a marker of energetic condition in several monkey and ape species, including chimpanzees (Deschner et al. 2008; Emery Thompson and Knott 2008; Emery Thompson et al. 2009; Girard-Buttoz et al. 2011). These studies highlight the validity of c-peptide as a reliable marker of energy balance (Emery Thompson 2017) since UCP levels covary with energy intake (Girard-Buttoz et al. 2011), body condition (i.e. body mass and body fatness, Girard-Buttoz et al. 2011) and energy expenditure since in wild macaques UCP levels tend to be negatively correlated to the height climbed (Girard-Buttoz et al. 2014a). High levels of UCP indicate positive energy balance and low levels indicate negative energy balance. Recently, investigations of chimpanzee socioecology using UCP revealed some contrasting patterns across and within populations. High-ranking males who compete the most intensely for access to and monopolization of fertile females had lower UCP levels than low-ranking males (Emery Thompson et al. 2009; Georgiev 2012). Males occupying a more productive habitat had higher UCP levels than their neighbouring community in a less productive area (Emery Thompson et al. 2009). Finally, female eastern chimpanzees had lower UCP levels when associating with more males than when associating with more females (Emery Thompson et al. 2014). These studies provide some contrast to a dichotomic view of socioecological theory by showing that (1) males can also be strongly impacted by food availability and (2) male-female competition might affect female energetic condition more than female-female competition.

The specific aim of our study was to understand the effects of feeding competition, dominance rank and mating competition on the energy balance of both males and females. We derived hypothesis and predictions from socioecological theories (summarized in Table 1). We hypothesized that females should be more affected than males in their energetic condition by food availability and that social factors (dominance rank and presence of oestrus females) should have different effects on males’ and females’ energetic conditions. We predicted that high-ranking females will have a more positive energy balance (higher UCP levels) than low-ranking ones and that this effect will be independent of the presence or absence of oestrus females in the party. In males, we predict that high-ranking individuals will have in general higher UCP levels than low-ranking males since they have priority of access to high-quality food (Whitten 1983; Isbell et al. 1999). However, in the presence of oestrus females, the higher ranking individuals, which are the ones which invest more in mating competition (Muller and Wrangham 2004), will have lower UCP levels than low-ranking individuals since mating competition and female monopolization is costly for males (Alberts et al. 1996; Georgiev 2012; Girard-Buttoz et al. 2014a, b). We studied one community of chimpanzees over 12 months in Taï National Park using a combination of phenological measures to assess food availability, behavioural observations to assess energy intake and energy spent on activities and UCP levels as a marker of energy balance.

Table 1 Predictions for our study derived from the socioecological model and summary of the outcomes of our analyses

Material and methods

Study site and subjects

We conducted the study from July 2017 to June 2018 at the Taï Chimpanzee Project (TCP) [5° 45′ N, 7° 07′ W (Wittig 2018)] in the Taï National Park Cote d’Ivoire. We observed North group, which consisted of 21 chimpanzees with 3 adult males, 1 subadult male, 7 adult females, 3 juveniles and 7 infants in July 2017. Males of 15 years and females of 13 years and older were considered adults, males between 12 and 15 were considered subadults (similar to Watts and Pusey 2002; Wroblewski et al. 2009).

Behavioural observations

PDV conducted all-day focal animal sampling (Altmann 1974) on all adult and subadult males (N = 4) and females (N = 7) of the community with an average observation time per individual of 149.4 h (range: 130.8–167 h), resulting in 1643 focal observation hours (see supplementary Table S1). It was not possible to record data blind because our study involved focal animals in the field.

During focal follows, PDV recorded the focal activities continuously and social interaction using Cyber Tracker software (https://www.cybertracker.org/). We considered the following activities divided into (1) activities providing energy: feeding and (2) activities leading to energy expenditure: travelling, climbing, aggressing, resting, grooming, drumming and playing (for detailed definitions of the activity states and behaviour, see Table 2). During each feeding bout, PDV recorded the duration, the type (i.e. plants or meat and for plants, the part consumed such as leaves, fruit and flowers) and species name of the consumed food and the number of food items ingested. When the focal animal was on the ground, PDV counted each food item ingested during the entire feeding bout. When the chimpanzees consumed honey, ant, larvae, small leaves and flowers, PDV counted the number of handfuls ingested. When the focal chimpanzee was in a tree or visibility was reduced on the ground, the total number of items ingested was estimated based on extrapolation derived from counts of the number of food items ingested during 2 min every 10 min. PDV made sure to cover all periods of the feeding bouts (beginning, middle and end). For large fruits, PDV estimated visually the weight of the whole fruit and the proportion obtained by the focal individual and for meat, he also visually estimated the weight of the prey and the proportion of meat obtained by the focal. While the focal was cracking nuts, PDV counted, when possible, the number of hits used to open each nut to calculate the energy spent for cracking nuts (see below). PDV collected (from the same individual food source, i.e. the same tree) and weighed food items to later estimate the quantity of food ingested by the chimpanzees. During every focal day, using the self-recording function of a GPS (Garmin GPSMAP 64s), PDV recorded the GPS track-log of the focal individual and determined later the daily travel distance and velocity. We set the GPS to record location at a time interval automatically adjusted to the focal individual locomotion (i.e. points were recorded more frequently when the individual was travelling and less frequently when the individual was resting). Despite this setting, the GPS sometimes recorded multiple positions while the individual was immobile due to GPS position error. In order to smooth the track-log and delete errors, we cleaned and smoothed the tracks by removing consecutive positions recorded continuously at the same location using an algorithm developed by Janmaat et al. (2013). While travelling, PDV recorded if the focal was carrying an infant or not. He also continuously recorded changes in vertical locomotion in trees or other structures to estimate the daily climbing height of focal individuals. We considered 6 vertical strata adapted from Buzzard (2006): 0 = ground, 1 = 1–5 m, 2 = 6–10 m, 3 = 11–20 m; 4 = 21–30 m and 5 = > 30 m.

Table 2 Definition of the activity states and behaviours

Since the presence of oestrus females can influence the dynamics of the community, including the dynamics of feeding competition (Muller et al. 2006; Sobolewski et al. 2013; Georgiev et al. 2014), PDV and experienced field assistants recorded, each day, the identity of the oestrus females in the community. Chimpanzees live in societies with a high degree of fission-fusion dynamics so that the party composition may change throughout the day depending on social and ecological events (Wrangham et al. 1996). Since community size can affect energy intake in chimpanzees through foraging efficiency (Isabirye-Basuta 1988), PDV recorded variation in the party composition continuously throughout the day when changes occurred to control for the potential effect of party size on chimpanzee energy balance.

Estimation of energy intake

We determined energy intake based on the rate of food items ingested by the daily focal individual. In total, PDV recorded 613 feeding focal hours comprising 3635 different feeding bouts. PDV was able to count the number of food item consumed by the focal individual during 2139 out of those 3635 bouts. We calculated an estimated food intake rate for each food species using the counts of the number of food items ingested. We then extrapolated the feeding rate to feeding time of each feeding bout where the target was in sight and feeding on an identified food item but its mouth was not visible such that PDV could not count the eaten items (i.e. during 1496 bouts). Based on a coefficient of dry weight computed in a previous study on the same chimpanzee population (N’Guessan et al. 2009), we assessed the dry weight of each food item from the fresh weight we measured in the field. For each food species, PDV collected and weighed several items and calculated an average weight per item for each species. Then, we calculated the total daily energy intake of the focal individual using the nutritional value of each food item as measured by N’Guessan et al. (2009) (see Supplementary Table S2) and by applying the following formula.

$$ {\mathrm{DE}}_{\mathrm{int}}=\sum \limits_{f=1}^n\left[{N}_{\mathrm{f}}\times {\mathrm{DW}}_{\mathrm{f}}\times \left(4\times \%{\mathrm{CP}}_{\mathrm{f}}+4\times \%{\mathrm{TNC}}_{\mathrm{f}}+9\times \%{\mathrm{L}}_{\mathrm{f}}+1.6\times \%{\mathrm{NDF}}_{\mathrm{f}}\right)\right] $$
DEint:

the total daily energy intake

N f :

the number of fruit f consumed in each feeding bout

DWf:

dry weight of one food item f consumed in each feeding bout

For a fruit f, we have the following: CP = crude protein; TNC = total non-structural carbohydrates; L = lipid; NDF = neutral detergent fibre.

The details of average daily energy intake estimated for males and females and of food availability each month are shown in Table S3.

Estimation of energy spent on activities

In our study, the energy expenditure only focused on estimating the energy spent on activities. We calculated the daily energy spent on activities (DEAexp) considering the energetic cost for resting metabolic rate (RMR) which is the amount of energy used for maintaining basal metabolic rate (BMR) during resting time and we added the energy expenditure of travelling (Et); climbing (Ec); nut cracking (En); eating (Ee) and social activities (Es). Since we do not know the exact weight of our study chimpanzees and since this information cannot be obtained non-invasively, we used a constant body mass for all the males of 46.3 kg which is the average body mass of a male western chimpanzee (Smith and Jungers 1997). Likewise, we assigned a body mass of 41.6 kg to each female which corresponds to the average body mass of a female western chimpanzee (Smith and Jungers 1997). The cost of each activity is calculated as a percentage of BMR (see below). BMR is estimated to be the same across all the study males and across all the study females. However, the equation below captures individual variation in daily energy spent on activities since it relates to the activity budget of each focal individual.

$$ \mathrm{DEAexp}=\mathrm{RMR}+\mathrm{Et}+\mathrm{Ec}+\mathrm{En}+\mathrm{Ee}+\mathrm{Es} $$

The averages daily energy spent on activities each month per sex are presented in the Supplementary material (Table S3).

Basal metabolic rate and resting metabolic rate

We estimated the BMR from a formula developed for chimpanzees by Pontzer et al. (2016). We referred to the relation between the RMR and BMR developed for orangutans (Pontzer et al. 2010, 2016) and we use the formula of chimpanzee BMR (Pontzer et al. 2016) to estimate the RMR of chimpanzees. We estimate BMR and RMR from the following formula where M is the body mass of the individual in kilogram:

$$ {\displaystyle \begin{array}{c}\mathrm{BMR}\ \left(\mathrm{kcal}/\mathrm{day}\right)=100.17\times {M}^{0.65}\\ {}\mathrm{RMR}\ \left(\mathrm{kcal}/\mathrm{day}\right)=\mathrm{BMR}/0.8\end{array}} $$

Energy expenditure for travelling (Et)

We estimated the travel energy expenditure using the formula implemented by Pontzer et al. (2014) which posits that a chimpanzee travelling with a velocity v consumes a volume of oxygen (VO2) where:

$$ \mathrm{V}{\mathrm{O}}_2=0.22\times v+0.06\ \mathrm{ml}\ {\mathrm{kg}}^{-1}\ {\mathrm{s}}^{-1} $$

We converted the volume of oxygen as 1 l of O2 = 4.8 kcal following Schmidt-Nielsen (1997). Then, we obtained the travelling energy expenditure using the following formula:

$$ \mathrm{Et}=\mathrm{V}{\mathrm{O}}_2\times {10}^{-3}\times 4.8\times M\ \mathrm{kcal}\ {\mathrm{s}}^{-1} $$

When females were carrying an infant while travelling, we adjusted the cost of locomotion for that female by taking into account the additional load, following the approach by Pontzer and Wrangham (2004). Carrying an infant increases the percentage of energy spent on travelling proportionally to the mass of the load (Taylor et al. 1980). For instance, if a mother is moving while carrying an infant weighing 20% of the mother’s body mass, it will spend 20% more oxygen than without the infant load (see Taylor et al. (1980) and Altmann and Samuels (1992) for details). Based on the body mass versus age curves obtained in eastern chimpanzees by Pusey (1990) and Pusey et al. (2005), we estimated the body mass of every carried infant in our study community using the following formula: Massinfant (kg) = 2 * infant age (years).

Energy expenditure for climbing (Ec)

To determine the energy spent on climbing, we used the estimation implemented by Hanna et al. (2008) giving the following formula where h is the height climbed and M is the body mass:

$$ \mathrm{Ec}\ \left(\mathrm{kcal}\right)=119.8\times {10}^{\hat{\mkern6mu} -3}\times 0.239\times M\times h $$

Energy expenditure for nut cracking (Ec)

By counting the number of hits used to open a nut during cracking sessions, we calculated a mean number of hits required to open one nut per individual and per nut species. We then extrapolated this mean across the entire cracking session. Günther and Boesch (1993) calculated the energetic cost of nut cracking in chimpanzees and found that chimpanzees use an amount of 0.04 kcal per hit to crack a nut of Coula edulis. Following Günther and Boesch (1993), we evaluated the energy spent on nut cracking by multiplying the number of hits necessary to crack a nut by the energy spent on one hit.

Energy expenditure for eating (Ee) and social activities (Es)

We estimated the energy used for eating and other activities with a formula of Leonard and Robertson (1997) who stated that for an activity performed during a time t, a coefficient D can be used to estimate the energy spent. The coefficient D varies with the activities: Deating = 1.38 and Dsocial = 2.35. Thus:

$$ \mathrm{E}a=\left(\mathrm{D}a\times \mathrm{BMR}\times t\right)\times {24}^{-1}\ \mathrm{kcal} $$

Assessment of food availability

An experienced field assistant collected phenological data on 944 fruit trees belonging to 82 species frequently eaten by the chimpanzees using the same methods with N’Guessan et al. (2009). We thus estimated the food availability index FAm, for each month “m” with the following formula from Anderson et al. (2002):

$$ {\mathrm{FA}}_{\mathrm{m}}=\sum \limits_{k=1}^n{D}_k{B}_k{P}_{km} $$
D k :

density of plant species k within the study area.

B k :

mean basal area (measured by trunk diameter at breast height) of species k within the study area.

P km :

percentage of trees of species k with mature fruits within the study area in a month.

Assessment of dominance rank

We calculated the dominance hierarchy using a modified version of the Elo-rating method (Neumann et al. 2011) developed by Foerster et al. (2016). We used the occurrence of unidirectional submissive pant-grunt vocalization (given by the lower ranking of the two individuals towards the higher ranking, Bygott 1979) to establish chimpanzee dominance ranks. We used all of the long-term data available from the Taï chimpanzee project consistently collected ad libitum by several field assistants, students and researchers since 2016. We calculated male and female Elo-rating scores separately and combined them in a single dominance hierarchy (see details in Mielke et al. (2017)), since all males are dominant over all females in the Taï chimpanzees (Gomes et al. 2009). We standardized all Elo-rating scores between 0 and 1 with 1 being the highest ranking individual and 0 the lowest ranking, on a daily basis. We then extracted the Elo-rating score of each individual on the focal observation day considered.

Urine sample collection and analyses

In addition to the behavioural data, PDV collected urine samples from the focal individual (one in the early morning and one in the late afternoon). We analyzed 175 samples for this study. PDV pipetted the urine directly from leaves or the ground into a 2-ml cryotubes. We made sure that the urine was not contaminated with faeces since contamination with faeces can severely alter urinary c-peptide measures (Higham et al. 2011). After collection, PDV placed the samples in a thermos cooler and froze the samples in liquid nitrogen within 12 h of collection at the research camp. We transported samples frozen on dry ice to the Laboratory of Endocrinology at the Max Planck Institute for Evolutionary Anthropology in Leipzig, Germany. In the laboratory, we thawed frozen urine samples at room temperature. After shaking for 10 s (VX-2500 Multi-tube Vortexer), we then centrifuged the samples for 5 min at 4400g (Multifuge Heraeus). We analyzed urine samples for c-peptide levels using C-Pep-EASIA KAP0401, which is a commercially available solid-phase enzyme-amplified sensitivity immunoassay kit from DIAsource, designed to measure C-peptide in human serum and already used in wild great apes (Surbeck et al. 2015; Wessling et al. 2018b). We added 100 μl of pure urine in duplicates into the wells and followed assay instructions. Intra- and inter-assay coefficients of variation of low- and high-value quality controls were 9.4% and 7.4% (N = 5) and 9.1% and 7.9% (N = 9), respectively. Samples were diluted and re-measured if their measured values did not fall into the linear range of the assay and simply re-measured if the two duplicate measures differed by more than 10%. C-peptide values were corrected according to creatinine content and are expressed in ng/mg creatinine.

Statistical analyses

We investigated the socioecological factors affecting chimpanzee overall energy balance and the different components of energy balance (i.e. energy intake, feeding time and energy spent on activities) using four different linear mixed models (LMMs). The structure of all the LMMs ran for this study is summarized in Table 3. In Model 1, we used UCP levels as the response variable. Each urine sample constituted a data point (N=175). We ran additional models to assess which components of the energy balance were impacted by the socioecological factors that we tested. We used the total daily energy intake (Model 2a), the percentage of observation time spent feeding (Model 2b) and the total daily energy spent on activities (Model 3) as responses in three separate LMMs. For Models 2a, 2b and 3, each data point was an observation day (N=158 days) and the response variable was assessed for each data point each day.

Table 3 Structure of the five LMMs used in this study

In all models, we tested the effect of four predictor variables: monthly food availability, sex, dominance rank and presence of oestrus females in the party (a proxy for the degree of mating competition) on the response. We assessed for each observation day whether at least one maximally swollen female was present in the party with the focal individual. We incorporated males and females within the same models to be able to include interactions between sex and our test predictors, and used the significance of such interactions to assess whether certain predictors have a stronger relationship to the energetic measures in one sex as compared to the other. We added an interaction between food availability and sex since we predicted that the energetics of females should be more influenced by food availability than males (Clutton-Brock and Harvey 1978; Sterck et al. 1997). We also incorporated a three-way interaction between rank, sex and presence of oestrus females into each of the four LMMs because we expected that the relationship between rank, mating competition and energetic conditions would differ for males and females since females are expected to compete more intensely over food than males while males are expected to allocate more energy to mating competition than females (Trivers 1972; Emlen and Oring 1977; Kappeler 2000).

In all models, we also added the average daily party size around the focal as a control variable. In addition, in Model 1, we added the time of sample collection to account for diurnal variation in UCP levels (Georgiev 2012). In Models 2a and 3, we also added an offset term for the observation time (in hours) to practically model energy intake and energy spent on activities per hour of observation (considering that the longer the individuals are observed each day, the more energy intake and energy spent on activities will be recorded).

In addition to these four models, we wanted to test directly the relationship between our measures of energy intake and energy spent on activities and UCP levels (a measure of energy balance) to assess indirectly the potential error pertaining to behavioural assessment of energetic parameters. We ran a fifth LMM (Model 4) (Table 3). We considered all urine samples (N = 175 samples) corresponding to the days (N = 98 days) when energy intake and energy spent on activities were measured. In all five LMMs, we incorporated individual ID as a random factor to avoid pseudoreplication arising from collecting several data points per individual.

We ran all the five LMMs (Baayen 2008) with a Gaussian error structure and identity link (McCullagh and Nelder 1989) in R software (version 3.5.3; R Development Core Team 2016) using the package “lme4” (Bates et al. 2015). For all models, we included the maximal random slope structure between each fixed predictor (test and control) and each random effect (Baayen et al. 2008; Barr et al. 2013). In each model, we tested for the overall significance of the test predictors by comparing the full model to a null model comprising all control predictors, all the random effects and random slopes but without the test predictors. We tested each full model against its corresponding null model using a likelihood ratio test (LRT, Dobson 2002). We then assessed the significance of each predictor variable using an LRT between the full model and a reduced model comprising all the variables except the one to evaluate. This process was repeated across all variables using the drop1 function. If the LRT revealed that one interaction had a P value > 0.05, we reran the model without this interaction and reassessed the significance of all the predictors.

We checked the model assumptions and verified that for all three models the residuals were normally distributed. We confirmed model stability by removing one level of each random effect at a time and recalculating the estimates of the different predictors which revealed comparable estimates to the original models. Finally, we tested for collinearity issues between the predictor variables using the function vif (variance inflation factor) from the package “car” (Fox and Weisberg 2011) which revealed no collinearity issue (all vif < 1.3 in all the models).

Results

Energy balance

In our first model testing the effect of socioecological factors on chimpanzee energy balance (UCP levels), we found that the full model was not significantly different from the null model (LRT: χ2 = 16.81, df = 11, P = 0.052). However, since the P value was very close to the arbitrary threshold of 0.05, we investigated the significance of the single predictors in the model, keeping in mind that the results should be interpreted with caution.

After removing the non-significant interaction between food availability and sex (χ2 = 0.73, P = 0.393) from Model 1, we found that food availability had a significant positive effect on energy balance in both sexes with UCP levels being higher when food was more abundant (Model 1: χ2 = 4.19, P = 0.040, Table 4, Fig. 1). Except for a few months in the year (e.g. October and November 2017 and February 2018), monthly chimpanzee UCP levels and monthly food availability covaried throughout the year (Fig. S1).

Table 4 Results of the four generalized linear mixed models examining the effect of ecological (food availability) and social (dominance rank, presence of oestrus females) on urinary C-peptide level (Model 1), energy intake (Model 2a), feeding time (Model 2b) and energy spent on activities (Model 3), and of a fifth model assessing the effect of energy intake and energy spent during activities on UCP levels (Model 4) in wild western chimpanzees
Fig. 1
figure 1

Effect of food availability on UCP levels in wild western chimpanzees. Each dot represents the average UCP levels (in ng/mg creatinine) for each individual each month. The area of the dot is proportional to the number of samples collected on that individual that month. The blue line depicts the model prediction for the effect of food availability on UCP levels (Model 1)

In Model 1, the interaction between individual dominance rank, sex and presence of oestrus females was significant (χ2 = 5.62, P = 0.018) indicating that the presence of oestrus female mediated the relationship between dominance rank and UCP levels and that this relationship differed between male and female chimpanzees (Table 4, Fig. 2). More specifically, higher ranking females had higher UCP levels than low-ranking females in the absence of oestrus females (Fig. 2a) but slightly lower UCP levels than low-ranking females when oestrus females were present (Fig. 2b). Among males, high-ranking males had lower UCP levels than low-ranking males (Fig. 2c) in the absence of oestrus females and higher UCP levels than low-ranking males in the presence of oestrus females (Fig. 2d). Throughout the study, we observed only one rank reversal between the number 2 and 3 in males and all other ranks remained stable in both sexes.

Fig. 2
figure 2

Relation between UCP levels and social (dominance rank*, presence of oestrus females) factors in males and females western chimpanzees. (a) UCP levels of females in the absence of oestrus females within the party; (b) UCP levels of females in the presence of oestrus females within the party; (c) UCP levels of males in the absence of oestrus females within the party; (d) UCP levels of males in the presence of oestrus females within the party. The dots in the graph represent the Elo-rating scores of each individual standardized within each sex. The black lines in each of the four panels indicate the model lines following the slope and intercept of the relationship between dominance rank and UCP levels predicted by Model 1 in the four conditions. Also note that while the interaction rank*sex* presence of oestrus female was significant in Model 1, the full-null model comparison in Model 1 was only close to significance (P = 0.052) so the relationships depicted here should be interpreted with caution

In addition to our test predictors, in Model 1, we did not find a significant effect of average daily party size on UCP levels (χ2 = 0.07, P = 0.788); in contrast, time of sample collection had a significant positive effect on UCP levels (χ2 = 5.27, P = .022) with UCP levels being lower in the morning samples than in the afternoon samples (Supplementary material, Fig. S2).

Energy intake and feeding time

Our calculation of energy intake based on behavioural observations of focal chimpanzee feeding behaviour was not significantly affected by the social (dominance rank and the presence of oestrus females) and ecological (food availability) variables we investigated (full-null model comparison in Model 2a, LRT: χ2 = 5.22, df = 9, P = 0.815). In contrast, the socioecological factors significantly affected chimpanzee feeding time (full-null model comparison in Model 2b, LRT: χ2 = 22.15, df = 9, P = 0.008). None of the interaction terms was significant in Model 2b (all P > 0.07). After removing the non-significant interactions from the model, we found that chimpanzees from both sexes spent less time feeding when food availability was high (χ2 = 8.13, P = 0.004, Fig. 3, Table 4) and when they were in large parties (χ2 = 5.51, P = 0.023, Fig. S3, Table 4).

Fig. 3
figure 3

Effect of food availability on time spent feeding in wild western chimpanzees. Each dot represents the average time spent feeding (in percentage of observation time) for each individual each month. The area of the dot is proportional to the number of observation days for that individual that month. The blue line depicts the model prediction for the effect of food availability on feeding time (Model 2b)

Energy that chimpanzees spent on activities

The socioecological factors under investigation in our study also significantly affected the energy that a chimpanzee spent on activities (full-null model comparison in Model 3, LRT: χ2 = 31.68, df = 9, P < 0.001). None of the interaction terms was significant in Model 3 (all P > 0.24). After removing these non-significant interactions from the model, we found that chimpanzees of both sexes spent more energy on activities when food was more abundant (χ2 = 18.23, P < 0.001, Fig. 4, Table 4) and high-ranking individuals of both sexes spent less energy on activities than lower ranking one (χ2 = 4.72, P = 0.030, Fig. 5, Table 4).

Fig. 4
figure 4

Effect of food availability on energy spent on activities in wild western chimpanzees. Each dot represents the average daily energy spent on activities (in kcal/day) for each individual each month. The area of the dot is proportional to the number of observation days for that individual that month. The blue line depicts the model prediction for the effect of food availability on energy spent on activities (Model 3)

Fig. 5
figure 5

Effect of dominance rank on energy spent on activities in wild western chimpanzees. Each dot represents the average daily energy spent on activities (in kcal/day) for each individual. The area of the dot is proportional to the number of observation days for that individual. The blue line depicts the model prediction for the effect of dominance rank on energy spent on activities (Model 3)

Relationship between UCP levels and behavioural measures of energy intake and energy spent on activities

We could not detect a significant effect of our behavioural measures of energy intake and energy spent on activities on UCP levels (full-null model comparison in Model 4, LRT: χ2 = 1.86, df = 2, P = 0.394).

Discussion

Our study provides new insights into the association between socioecological factors and energetic conditions in western chimpanzees. We found an effect of food availability on energy balance in both females and males. As predicted, social factors such as dominance rank and the presence of oestrus females impacted males and females differently. As expected from socioecological theories (Sterck et al. 1997; Koenig 2002), dominance rank was positively correlated with energy balance in females. Yet this positive effect was only found in the absence of oestrus females in the party. The effect disappeared in the presence of oestrus females. In contrast, in males, we found a negative effect of dominance rank on energy balance when there was no oestrus female present but a positive effect in the presence of oestrus females.

It is important to note that despite the significance of our test predictor in Model 1 assessing the socioecological factors influencing UCP levels in chimpanzees, the full versus null model comparison was not significant so the pattern described here should be considered with caution. In the rest of the analysis, we could not isolate a clear contribution of our behavioural measure of energy intake and energy spent on activities to explain the combined effect of sex and dominance rank on energy balance measured using UCP levels. In fact, we could not detect a combined effect of rank and sex on our behavioural measures of chimpanzee energy intake and energy spent on activities with individuals from both sexes feeding for shorter durations and spending more energy on activities when food availability was high and feeding longer in smaller parties. We also found that high-ranking male and female chimpanzees spent less energy on activities. When testing directly the relationship between our hormonal and behavioural energetic measures (Model 4), we could not detect a significant link between UCP levels and energy intake or energy spent on activities.

Interestingly, and contrary to the prediction derived from socioecological theories (Emlen and Oring 1977; Clutton-Brock and Harvey 1978; Clutton-Brock et al. 1989; Sterck et al. 1997; Kappeler 2000; Koenig 2002; Tobias et al. 2012), in our study, the positive effect of food availability on feeding time and energy balance was of a similar magnitude for male and female chimpanzees since the interaction between sex and food availability was not significant in Model 1 and Model 2b. N’Guessan et al. (2009) found the same relationship between energy balance and food availability in both sexes in the same chimpanzee population, but based on behavioural assessment of energy balance, although similar results were found in male eastern (Emery Thompson et al. 2009, 2014) and western chimpanzees (Wessling et al. 2018a), assessing the difference between sexes was not possible, since Emery Thompson et al. (2009) used males only and Emery Thompson et al. (2014) and Wessling et al. (2018a) did not consider the interaction between food availability and sex. In contrast, our study suggests a similar magnitude of this effect in both sexes. Yet, this result calls for caution since the full and null model comparison in Model 1 is slightly above the significance threshold with P = 0.052, and our data set consists of a small sample size. Our results, however, show that, under certain conditions, ecological factors can impact male and female energetics equally, and should encourage future studies to replicate this analysis in larger chimpanzee communities and with larger datasets.

The relationship between UCP levels and food availability was significant but not fully linear and, in certain months, UCP levels were high when the food availability was low and vice-versa (Fig. S1). Those months were characterized either by a low food availability associated with an abundance of highly calorific fruit (i.e. Coula edulis) or by frequent hunting which is relatively independent of food availability (Stanford et al. 1994; Boesch and Boesch-Achermann 2000; Gilby and Wrangham 2007) but may lead to high energy gains. Indeed, it was suggested that the nutritional importance of meat and bone marrow consumption in the chimpanzee diet could be one of the main factors which drove the evolution of chimpanzees hunting behaviour (Boesch 1994; Stanford et al. 1994; Mitani and Watts 2001). This indicates that it was not only the availability of food in large quantities (our measure) which increased individual energy balance but also the food quality (i.e. its nutritional value). Therefore, we encourage future studies to develop phenological indexes which incorporate not only food quantity but also a nutritional value index.

Our results on the behavioural measures of energy intake and energy spent on activities highlight the difficulty in using these measures as a proxy for individual energy balance. In fact, we found that chimpanzees spent less time feeding and spent more energy on activities when food availability was high but still had a more positive energy balance during these periods (higher UCP levels). Our results indicate that chimpanzees may have evolved a behavioural adaptation documented in other studies consisting of feeding longer (N’Guessan et al. 2009) and on lower quality food (Doran 1997) during periods of food scarcity while refraining from engaging in energetically costly activities. These results also emphasize the value of non-invasive physiological measures of wild mammal energetic status in addition to behavioural measures.

Besides the ecological effect, we also found an influence of dominance rank on individual energetic conditions which differs between males and females and was conditional on the presence or absence of oestrus females. In the absence of oestrus females in the parties, the energy balance of high-ranking females was higher than that of low-ranking females. This finding is in line with general socioecological theories predicting strong female-female competition for access to food resources (Emlen and Oring 1977; Clutton-Brock and Harvey 1978; Clutton-Brock et al. 1989; Tobias et al. 2012), and, as a result, an effect of dominance on food acquisition as observed in many primate species (Barton and Whiten 1993; Saito 1996) including western female chimpanzees (Wittig and Boesch 2003) and in spotted hyenas (Watts et al. 2009). This proximate effect of dominance rank on female access to food is likely to be the main mechanism leading to enhanced fitness in high-ranking females documented in several mammalian species (Pusey et al. 1997; King and Allainé 2002; von Holst et al. 2002; Majolo et al. 2012), including primates (Majolo et al. 2012). Our study contributes to our understanding of the mechanisms linking high dominance status and fitness benefits in chimpanzees. In particular, we show that the priority of access to high-quality food by high-ranking females (Pusey et al. 1997; Murray et al. 2006) leads to more positive energy balance in these females (our study) which is one potential main mechanism explaining why high-ranking females produce offspring which are more likely to survive (Goodall 1977, 1986). Surprisingly, in our study, this energetic advantage of being a high-ranking female disappeared when oestrus females were present (60% of the observation days) suggesting that the social dynamics associated with male-male mating competition impacts female energetics as well. The presence of males around oestrus females may increase feeding competition and reduce female feeding opportunities, as observed in eastern chimpanzees where females had lower UCP levels when associating with more males than when associating with more females (Emery Thompson et al. 2014). In several animal species in which males can obtain coercive mating, conflicts of interest between the sexes are likely to arise (Trivers 1972; Parker 1979; Clutton-Brock and Parker 1995), and receptive females are more often the victims of aggressions. Male chimpanzees often use coercion to access fertile females in order to increase their mating success, which in turn incurs a high energetic cost for the female victim (Muller et al. 2007). In chimpanzees, high-ranking females are often the oldest (Wittig and Boesch 2003), and therefore, the most experienced mothers whose offspring are more likely to survival to adulthood (Boesch 1997; Pusey et al. 1997). This may be why males prefer high-ranking females as mating partners (Muller et al. 2006). In turn, high-ranking females receive more aggression when fertile compared to low-ranking females and thus bear a stronger energetic cost which may limit the benefits of dominance rank on their energetic condition. However, the aggression rate of males towards oestrus females is unlikely to be the main driver of the observed pattern in our study community because coercion is rare in Taï (Stumpf and Boesch 2005). Even though it is relatively rare, male aggression towards oestrus females is used by less preferred males in Taï, but does not seem to constrain female mating choice (Stumpf and Boesch 2005). In our study, the rank effect did not only disappear for fertile females since the effect of dominance rank on female energy balance was alleviated across all females. This indicates that the presence of maximally tumescent females modifies the social, foraging and ranging dynamic of the community in a way that removes the energetic advantages of high-ranking females. In fact, in Taï, the presence of oestrus females increases the number of males but also the number of females in the party (Anderson et al. 2002) which potentially generates different competitive dynamics. This is mostly stemming from the high attractiveness of receptive females for both males and females leading to males following the ranging patterns of receptive females (Emlen and Oring 1977; Mitani et al. 1996; Kappeler 2000). In turn, this may lead to a shift in “power” by placing the foraging decisions and the priority of access to richer food patches in the hands of the fertile females rather than systematically in the hands of high-ranking females. However, the reduced rank effect on energy balance across all females during periods of mating competition may be due to the influence of a few popular females, given our sample size. Future studies should investigate this in more detail, focusing on the effect of male-male mating competition on female access to food and grouping patterns using behavioural measures.

In contrast to the female pattern, we observed a negative relationship between dominance and energy balance in males in the absence of maximally tumescent females, with high-ranking males having lower UCP levels than low-ranking males. This pattern could reflect (1) a lower energy intake for high-ranking males which was not captured by our behavioural measures and/or (2) higher energy expenditure for high-ranking males which is unlikely given our finding that high-ranking males spend less energy on activities than low-ranking males (Model 3) and/or (3) a higher metabolic rate in higher ranking than in lower ranking males stemming from higher physiological or psychological stress (Hollstein et al. 2020) as evidenced by higher levels of cortisol reported in high-ranking male eastern chimpanzees (Muller and Wrangham 2004; Gesquiere et al. 2011). Regardless of the precise effects causing more negative energy balance in high-ranking male chimpanzees, our results, in accordance with other studies on eastern chimpanzees (Emery Thompson et al. 2009; Georgiev 2012), emphasize the trade-off between investment into social status (Goodall 1986) and energetic allocation for other functions such as body maintenance (Komdeur 2001; Low 2006; Schubert et al. 2009). Interestingly, the presence of oestrus females reversed the effect of dominance status on males’ energy balance. This specific result differs from what was documented in eastern chimpanzees where high-ranking males had lower energy balance (UCP levels) than low-ranking males also in periods with high male-male mating competition (Emery Thompson et al. 2009; Georgiev 2012). Eastern male chimpanzees, regardless of their rank, appear to suffer an energetic cost when tumescent females were present in the party since they reduced their feeding time (Georgiev et al. 2014). Feeding time reduction might be associated with an increased focus on the sexually receptive female as was documented during mate guarding in macaques and baboons (Bercovitch 1983; Alberts et al. 1996; Girard-Buttoz et al. 2014a, b). In our chimpanzee population in Taï, males coerce females less than in the eastern populations (Boesch et al. 2008) and females have more opportunity to express a choice in their mating partners. Across mammals, including primates, females often prefer to mate with high-ranking males (Cox and Le Boeuf 1977; Andersson 1994; van Schaik 1996; Berglund et al. 1996; Sterck et al. 1997; Knott et al. 2010) which makes it even harder for low-ranking males to access female mating partners in Taï than in eastern chimpanzees where the female choice is reduced. Therefore, low-ranking male chimpanzees in Taï might have to monitor female behaviour (thereby reducing their feeding time), avoid retaliation from high-ranking males via engaging in sneaky copulations (Nurnberg et al. 1994; Setchell 2008) and gain the favour of females to be chosen for these copulations, such as through food sharing (Gomes and Boesch 2009). In turn, this could lead to a higher energetic cost of the presence of oestrus females for low-ranking than for high-ranking males. Future studies should investigate whether low-ranking males are truly constrained in their feeding time during periods of intense male-male competition.

Genetic relatedness is potentially an additional factor contributing to the observed pattern. In fact, our study community was composed by only four males and eight females and the mating competition was possibly stronger between low-ranking males because each of them was genetically related to one or two females in the community. The alpha males copulated with 100% of females while the beta copulated with 87% and the two low-ranking males in turn copulated with 75% of the females showing that the mating competition is not the same for all males. The low-ranking males had kin in the community and did not compete for access to fertile females when those females were kin which may have created periods of lowered mating competition for the high-ranking males. In turn, this may have alleviated some of the energetic costs of mating competition for higher ranking males and explained the observed pattern in our study.

It is important to note that despite its effect on energy balance, we could not find a significant effect of dominance rank on chimpanzee energy intake. We found a significant negative effect of dominance rank on the energy spent on activities by chimpanzees but this effect (unlike the effect on UCP levels) did not differ between males and females and was not depending on the presence of oestrus females. Therefore, it is likely that the lack of fit between the statistical models using behavioural measures of energy intake and energy spent on activities and the results on UCP levels are due to the level of error in estimating energy intake and energy spent on activities using behavioural observations. In fact, several errors can be introduced at several steps in estimating energy intake from digestibility of wild food (Conklin-Brittain et al. 2006), variation in fruit size, portion of the food item eaten, differential nutritional value of the food items in different parts of the canopy (e.g. for fruits (Houle et al. 2007, 2014; Gullo et al. 2014)). With such bias, it is possible that despite the best efforts to collect detailed feeding data, a crude measure of feeding time may more clearly reveal some socioecological phenomenon (linked to the party size), as we observed. Thus, our results on energy spent on activities but even more one energy intake should be taken with caution. It is likely that under other conditions, sex, rank, food availability and presence of oestrus females do not affect energy intake. UCP measures, although prone to some measurement error, appear a more suited integrative measure to reveal the socioecological processes underlying primate energetics, at least in our study.

Conclusion

Altogether, our results indicate that feeding and mating competition affect male and female chimpanzee energetic status. In turn, this suggests that social strategies and behavioural adaptation can prove useful adaptations to sustain high energetic demands due to ecological and social constraints. Our findings emphasize the importance of ecological factors such as food availability for male energetics and we encourage future study to integrate both sexes in socioecological analyses. On the other hand, we also characterized the importance of the degree of mating competition for the energetics not only of males but also of females since the presence of oestrus females alleviated the energetic advantages of being high ranking in female chimpanzees. Our results also emphasize the complexity of the effect of dominance rank on males energy balance and highlight how the community size and degree of genetic relatedness between individuals can be important in shaping this relationship and producing contrasting patterns to those observed in larger communities. To conclude, we would like to encourage future studies on feeding competition and socioecological energetics to incorporate both sexes and compare male and female patterns. This enterprise will provide a more comprehensive and balanced insight into the complex mechanisms impacting wild animal energetic status. Such studies can provide some contrasts and nuances to some tests of general socioecological theories which often focus on the energetic costs of mating competition for males and of feeding competition for females and do not assess these effects in both sexes simultaneously.