Introduction

A basic assumption of socioecological theory is that females compete for access to nutritional resources and males compete for access to females (Kappeler and van Schaik 2002). Predation risk and feeding competition are thought to be among the most important ecological factors that contribute to the evolution of female social relationships (van Schaik 1989; van Schaik and van Hooff 1983). Predation risk is argued to lead diurnal female primates to form groups for safety (van Schaik 1983), which in turn leads to intragroup competition for essential resources (Terborgh and Janson 1986; van Schaik 1989). Relationships among females are thought to be influenced by the type and strength of the feeding competition they face within and between their social groups, and this is affected by the distribution of resources (Chapman et al. 1995; Isbell 1991; Koenig 2002; Sterck et al. 1997).

Food intake can have important effects on a female’s survival and lifetime reproductive success and so the ability to dominate and exclude competitors from food resources is advantageous under certain conditions (Sterck et al. 1997; Trivers 1972). When resources are evenly distributed or plentiful, indirect (scramble) competition is expected (Kappeler and van Schaik 2002). Display of aggressive behavior and dominance over others is not expected to be advantageous under this type of competition and so female dominance relations are expected to be ill-defined or not identifiable (Janson 1985; Janson and van Schaik 1988), and hierarchies (if detectable) are expected to be indistinct and non-linear (Sterck et al. 1997). However, when resources are clumped and monopolizable then direct (contest) competition is expected to occur (Janson and van Schaik 1988). This type of competition is likely to result in agonistic interactions between females and the formation of hierarchical social relationships (Sterck et al. 1997). Females with despotic relationships have established dominance relations with one another and linear dominance hierarchies (Sterck et al. 1997).

Egalitarian relationships are expected between females in folivorous species as the primary food resource (leaves) is thought to be evenly distributed, abundant, and therefore unmonopolizable (Isbell 1991; cf. Koenig et al. 1998). Frugivory is typically thought to favour despotic relationships because ripe fruits tend to be high-quality, monopolizable and, therefore, contestable. Some studies have shown that agonism occurs disproportionately during feeding in primates, particularly when feeding on fruits (Cords 2000; Klass and Cords 2015). Studies have purported that dietary categories (e.g., degree of frugivory) can provide information as to the contestability and distribution of food and therefore give an indication of the expected level of feeding competition within a social group (McKenna 1979; Snaith and Chapman 2007). However, in a cross-taxa analysis, the expected association between high rates of agonism and frugivory in female primates was not generally supported (Wheeler et al. 2013) and some folivorous female primates show relatively low rates of agonism but still form decided dominance relationships (Gorilla beringei: Robbins et al. 2005; Trachypithecus phayrei: Koenig et al. 2004).

In addition to the distribution and quality of food, ecology can also indirectly affect female social structure by its effect on group size (Terborgh and Janson 1986) and there is notable variation in group size across the primate Order (Campbell et al. 2011; Kappeler and van Schaik 2002). The intensity of intragroup competition levels may increase with group size (Klass and Cords 2015; Koenig and Borries 2006; Terborgh and Janson 1986) and studies have found that primates living in larger groups are characterized by higher rates of agonism (Cowl and Shultz 2017; Wheeler et al. 2013). Furthermore, the demographic factor, female philopatry, is thought to be involved in the formation of female relationships (Isbell 2004; Kappeler and van Schaik 2002; Sterck et al. 1997). In species in which females disperse, female relationships are predicted to be egalitarian. In species in which females remain in their natal group, relationships vary depending on the levels of intra- and intergroup contest competition (Sterck et al. 1997).

Although the broad-scale descriptive classifications of female dominance relationships used by Sterck et al. (1997) are useful, there is considerable variation in how these relationships manifest themselves between (and sometimes even within) species (Klass and Cords 2015; Wikberg et al. 2013). This observed variation is beyond the descriptive scope proposed by Sterck et al. (1997) and species often do not conform to some predictions of the socioecological models described above (Klass and Cords 2015; Robbins et al. 2005; Wheeler et al. 2013; Wikberg et al. 2013), leading researchers to position species into the categories which they most closely match. For example, in a study on a population of predominantly folivorous Colobus vellerosus (a species with facultative female dispersal: Wikberg et al. 2014), occasional foraging on temporally available, clumped foods lead females to form individualistic dominance hierarchies of intermediate strength and high directional consistency (Wikberg et al. 2013).

The shortcomings of current socioecological models highlight the importance of considering quantitative measures of female dominance behavior (e.g., hierarchical steepness and linearity, directional consistency in dominance interactions and rates of agonism), as these may better represent the variation we see across primate species (Klass and Cords 2015). Unlike hierarchical linearity and steepness (Klass and Cords 2011), the directional consistency of female dominance relationships (measured via the directional consistency index [DCI]), is not sensitive to unknown relationships and is therefore considered the most accurate measure of despotism (Koenig et al. 2013; van Hooff and Wensing 1987). Additionally, because variation in ecology alone is not able to fully explain the observed differences in agonism and dominance patterns among female primates, it is necessary to also consider the effects of phylogeny when examining the effects of ecology on behavior (Klass and Cords 2015; Thierry 2008).

Reproductive competition can be a strong selective force in females (Rosvall 2011) but most studies investigating sexual competition focus on males (Clutton-Brock 2009). Studies in primates have revealed associations between the development of female ornamentation and weaponry and the form of reproductive competition between females (Clutton-Brock and Harvey 1976; Clutton-Brock and Huchard 2013). A study in wild female chacma baboons revealed that aggression between females was associated with access to mates (Huchard and Cowlishaw 2011). Females received aggression from other female group members at higher rates when they were reproductively active (displaying a sexual swelling or being mate-guarded by a male). Where relationships between females were traditionally considered to be the result of competition for access to resources, this suggests that sexual competition also has some involvement in shaping female social relationships (Huchard and Cowlishaw 2011). When investigating aspects of intrasexual competition between females, it is thus important to consider the impacts of a species’ mating system, alongside ecological factors relating to social structure and dominance relationships (e.g., Huchard and Cowlishaw 2011). Indeed, sexual competition, in addition to competition for access to food resources, may shape female social relationships, even in large polygynandrous primate groups in which females are the high investing sex and males are the principal competitors (Huchard and Cowlishaw 2011).

Aggression between females is typically less frequent and intense than between males (Clutton-Brock and Huchard 2013) but this general pattern is often not the case for cooperatively breeding primate species where a single dominant female is able to monopolize reproduction (Garber 1994; Saltzman et al. 2009). Among primates, cooperative breeding is largely limited to the small-bodied callitrichine monkeys characterized by monogamous or polyandrous mating systems. Variance in reproductive success may be higher among females than males in cooperative breeders and so reproductive competition is likely to be a major selective factor for monogamous and polyandrous females (Garber 1994). Indeed, in monogamous species both sexes often play a role in territory and mate defense and are therefore mostly intolerant of same-sex intruders (Garber 1994), whereas in polyandrous species, it is typically the breeding female that is intolerant of same-sex competitors and males have a more relaxed approach to each other (Garber 1994; van Hoof and van Schaik 1994). These higher stakes for females could explain why females of polyandrous and monogamous species tend to be highly territorial and intolerant of same-sex intruders (Garber 1994).

Another social factor possibly influenced by a species’ ecology is the relationships between males and females. In mammalian species, male dominance over females is the norm, and this is true of most primates (Dunham 2008; Kappeler 1993; Pereira et al. 1990). In some primate species, however, the relationships between the sexes are egalitarian in nature or there is no clear dominance of one sex over the other and the sexes are described as codominant (e.g., callitrichines and gibbons; Carpenter 1940; Koba et al. 2012; Reichard and Barelli 2008; Savage et al. 1988). Such dominance appears to correlate with monogamous and polyandrous mating systems, with some exceptions (Eulemur collaris: Balestri et al. 2014; Brachyteles arachnoides: Strier 1990). Female dominance, while rare, is most prevalent among Lemuriformes (Eichmueller et al. 2013; Kappeler 1993; Pereira et al. 1990) and occurs in some haplorhines such as Pan paniscus (Parish 1994) and some species of Saimiri (Mitchell 1994), and this pattern does not appear to be associated with particular mating systems (Kappeler 1991, 1993; Petty and Drea 2015).

While the hypothesized ecological factors contributing to aspects of social relationships among females and between females and males in primates have been discussed extensively, the role that underlying proximate mechanisms may play in regulating behaviors promoting particular dominance relationships is poorly understood. In many of the hypotheses that attempt to explain the link between social variation and ecological factors, agonistic behavior plays an integral part (Klass and Cords 2015; Sterck et al. 1997). Animal studies implicate prenatal androgen effects (PAEs) as having some influence over the shaping of an individual’s tendency towards future aggressive behavior (Mazur and Booth 1998), and higher PAEs have been known to increase offspring competitive abilities (Burley and Foster 2004). In species in which competition for resources is high, the ability to out-compete challengers is highly advantageous. Therefore, high PAEs may enhance competitive behaviors in female primates that live in large groups (especially those with despotic dominance relationships), and if degree of frugivory is indeed a good proxy for the intensity of direct competition experienced by females, then one might expect greater selection for PAEs in females of highly frugivorous species.

It is not yet known whether PAEs play a role in intersexual dominance patterns in primates. In the female dominant Hapalemur alaotrensis, females are more aggressive than males (Waeber and Hemelrijk 2003) and a study on six captive Eulemur species found that in female dominant species, females directed more dominance and aggressive behavior toward their male partner than they received whereas codominant females showed less dominance behavior than and equal rates of aggression as their male counterparts (Petty and Drea 2015). Species characterized by female dominance also had more masculine androgen profiles than females of codominant species and the authors suggest a hormonal mechanism may underlie the evolution of female dominance (Petty and Drea 2015). Therefore, in comparison to male dominated or codominant species, females belonging to species in which female dominance is the rule could be expected to be more behaviorally masculinized, as selection should favor increased expression of aggressive and competitive behavior.

The ratio between the lengths of the second and fourth digits of the hands (2D:4D ratio) is a proposed negative correlate of PAEs and positive correlate of prenatal estrogen effects (POEs) (Zheng and Cohn 2011). Artificially increasing PAEs in female Macaca mulatta resulted in masculinized (lower) right hand 2D:4D ratios, and the degree of this masculinization correlated strongly with prenatal androgen-dependent trait anogenital distance (Abbott et al. 2012; Thankamony et al. 2016), although males have been observed with higher 2D:4D ratios than females in this species (Baxter et al. 2018). In 37 species of anthropoid primate, greater sexual dimorphism in 2D:4D ratio was observed in non-pair-bonded species in comparison with pair-bonded species and all species displayed significant sex differences in 2D:4D ratio (Nelson and Shultz 2010). This pattern was not mirrored in a study investigating sexual dimorphism in the second and fourth metacarpal ratio in six species of nonhuman primate, but it is likely that the growth of the metacarpals and phalanges respond differently to PAEs (Hart 2018). Low 2D:4D ratios (inferred high PAEs) are associated with higher social rank and behaviors linked to dominance and aggression in cercopithecine female nonhuman primates (Howlett et al. 2012, 2015; Hurd et al. 2011; Nelson et al. 2010) and humans (Benderlioglu and Nelson 2004). Using classifications of female social relationships (Sterck et al. 1997), female 2D:4D ratio was found to decrease with increasing intrasexual competition levels (Nelson and Shultz 2010) and these findings are concurrent with higher PAEs being of adaptive value for female primates experiencing high levels of contest competition (Nelson and Shultz 2010).

Studies suggest there is a link between high 2D:4D ratios (low inferred PAEs) and the increased postnatal expression of behaviors linked with oestrogen and its associated neuropeptides oxytocin (OT) and vasopressin (VA) (Fink et al. 2007; Fisher et al. 2010; Williams et al. 2003). Many of the behaviors characteristic of a monogamous mating system such as pair bonding, partner preference, and parental care are all facilitated through the actions of estrogen, OT, and VA (French et al. 2018; Lee et al. 2009; Neumann 2008), potentially resulting in increased selection for reduced PAEs and/or increased POEs in monogamous and polyandrous female primates. Consistent with this, pair-bonded anthropoid primates have generally higher 2D:4D ratios (indicate lower PAEs) than non-pair-bonded species (Nelson and Shultz 2010). However, in describing species only in terms of pair-bonded or non-pair-bonded, much of the variation in mating systems is lost. Primates display a diverse array of mating systems and distinguishing between these (and the different forms they take) may elucidate how variation in female sexual competition levels affects the strength of selection for PAEs. Because the classifications of female social relationships used by Sterck et al. (1997) have been criticized due to inconsistencies between model predictions and observed behavior (Klass and Cords 2015; Thierry 2008; Wikberg et al. 2013) and the 2D:4D ratio measurement methods used by Nelson and Shultz (2010) could be improved upon in terms of accuracy and reliability (Allaway et al. 2009; Kemper and Schwerdtfeger 2009; Ranson et al. 2013; Voracek et al. 2007), a reexamination of the relationship between female social structure and PAEs, preferably using quantitative measures of female dominance behavior and more reliable 2D:4D ratio data, is desirable.

This study investigates variation in 2D:4D ratio in both strepsirrhine and haplorhine female primates. We investigated variation in 2D:4D ratio in female primates using improved 2D:4D measurement methods and quantitative measures of female social structure. We aimed to investigate the relationship between variation in PAEs (inferred from 2D:4D ratios) and female intrasexual competition and intersexual dominance relationships across primates. We tested five hypotheses, based on the assumption that 2D:4D ratio correlates negatively with PAEs. First, if food contestability affects female competitive relationships, females in species that depend on foods thought to elicit contest competition should face selection for higher PAEs than females in species that do not depend on such foods. Thus, we predict that highly frugivorous species have lower 2D:4D ratios than less frugivorous species. Second, if direct competition increases with group size, we predict that species that live in larger groups will display lower 2D:4D ratios than those that live in smaller groups. Third, if behavioral traits associated with higher PAEs are favored in females in species characterized by a high degree of despotism and high rates of female–female agonism, then we predict that the directional consistency in dominance interactions and rate of agonism will both be negatively related to 2D:4D ratio. Fourth, if species in which females experience higher competition for mates are likely to benefit from higher PAEs, we predict that monogamous females will have the highest 2D:4D ratios, followed by polyandrous females, with polygynous and polygynandrous females having the lowest 2D:4D ratios. Finally, PAEs play a vital role in the masculinization of behavioral tendencies (Goy et al. 1988; Pomerantz et al. 1985, 1986; Wallen 2005; Wallen and Hassett 2009), so we hypothesize that higher PAEs may be important for the maintenance of intersexual dominance relationships among primates (particularly female dominance) and predict that female 2D:4D ratios will be lowest in species in which females are the dominant sex (Table I).

Table I The socioecological variables investigated and their predicted relationships with PAEs and the 2D:4D ratio in female primates

Methods

Study Subjects

We collected hand images between March 2016 and November 2017 on 488 individual female primates across 71 species housed in 29 zoos, wildlife/safari parks, and primate research centres in the UK, Germany, Italy, Netherlands, Spain, and one sanctuary in South Africa. We also included published 2D:4D ratio data from 25 wild female chacma baboons (Papio ursinus) (Howlett et al. 2015). Animals included in the study were of juvenile (53 individuals), subadult (21 individuals), and adult (439 individuals) ages. The use of animals that are not fully mature in this study is justified as the 2D:4D ratio is fixed early in prenatal development (Galis et al. 2010) and does not change appreciably during puberty (Králik et al. 2014; Manning 2002; Manning et al. 2003, 2004).

2D:4D Ratio Measurements

We used the digital photographic and computer-assisted image analysis software method as described in Howlett et al. (2015) to obtain 2D:4D ratio measurements. We collected images of primate hands using the “free photo” method as described in Howlett et al. (2015) and photographs of the ventral surface of the animals’ hands were taken against a clear Perspex® sheet on occasions when animals were being handled for veterinary treatment (Fig. 1). For each individual, we identified three photographs per hand in which digits were in the optimal positions and for purposes of data reliability (Allaway et al. 2009), the first author carried out all measurements of 2D:4D ratios from these photographs using the software program ImageJ.

Fig. 1
figure 1

Image showing an example of computer-assisted measurement of hands using ImageJ software. White lines indicate the path of measurement for each digit from the basal crease to the tip of the extended digit. The images above present a Perspex® photograph of the left hand of an anesthetized male coppery titi monkey (Callicebus cupreus).

For each individual, we measured the second and fourth digits five times per photograph, giving a total of 15 measurements for each digit, and we used the mean of these 15 measurements as the measurement for that digit. We calculated the 2D:4D ratio of each hand by dividing the length of the second digit by the length of the fourth digit for each individual. We calculated the mean 2D:4D ratio (M2D:4D) by averaging the right 2D:4D ratio (R2D:4D) and left 2D:4D ratio (L2D:4D). We then took the mean of individual within species 2D:4D ratio measurements to obtain R2D:4D, L2D:4D, and M2D:4D data for each species. These species means are the 2D:4D ratio measures we used in analyses and the final data set comprised the R2D:4D, L2D:4D, and M2D:4D of 513 individual female primates across 71 species (Table II).

Table II Data set with female right, left and mean 2D:4D with standard deviation (SD) for each primate species, substrate use, mean female body mass (g), mean percentage (%) of fruit in diet, mean group size, female rate of agonism (per hour of observation time), female directional consistency index (DCI), and intersexual dominance patterns

Anatomical Considerations

Variation in primate hand morphology is associated with substrate use (Jouffroy et al. 1993; Lemelin and Schmitt 1998; Richmond 2007) and so we collated substrate use data for each species from the literature. Primates were classified as arboreal (68–100% arboreal), arboreal/terrestrial (34–67% arboreal), or terrestrial (0–33% arboreal). Where substrate use data were not available as percentages, we used categorical classifications given by the authors of the publications (Electronic Supplementary Material [ESM] Table SI).

In humans, it has been suggested that men have lower 2D:4D ratios because they have longer digits than women (Kratochvíl and Flegr 2009; Lolli et al. 2017; cf. Manning 2010), implying that a link exists between 2D:4D ratio and body size. The 2D:4D ratio is fixed early in prenatal development and is generally stable throughout postnatal growth (Galis et al. 2010; Manning and Fink 2018), making it unlikely that allometry has an extensive influence over its development. However, to account for possible effects of body size on 2D:4D ratio, we collected mean adult female body mass (in grams) data for each species and controlled for this variable in all analyses (Table II; ESM Table SII).

Diet and Group Size

To test predictions related to frugivory, we collated data from the literature on the percentage of fruit in each species’ diet (Table II; ESM Table SIII). Percentage of fruit in diet data describes the percentage of time spent feeding/foraging on fruit based on observational data. Nuts, seeds, and seed pods are typically high value, patchily distributed and contestable in the same way as fruits and so were also included in this category. We obtained some data from sources using different but equivalent methods and calculated the percentage time spent feeding/foraging on fruit from these data, e.g., where data were given as relative frequencies. Equivalent methods are those based on scan sampling, where the interval between samples is short relative to the mean duration of the behavior. This method provides estimates of time spent feeding on food items (Martin and Bateson 2007), as do methods in which dietary intake is determined through the proportion of all scan samples in which a food item was recorded as being fed upon (Davies et al. 1999).

We obtained data on mean species group size (includes both sexes and all age classes) from the literature. If more than one value was available for a species (e.g., representing different populations), we took the mean of these as the mean group size value for the species as a whole. Where data described the size of mixed-species groups, we included only the number of individuals of the species in question. For some species that live in large multilevel groups, mean group size data describes smaller grouping levels, e.g., band level of social structure in Papio hamadryas (Table II; ESM Table SIV).

Female Social Variables

We collated data on DCI and rate of agonism from the literature. DCI describes the proportion of interactions in the more common direction in each dyad (i.e., when the outcome of an agonistic encounter is in favor of the individual in the dyad who most often wins) among females in a social group. It is calculated by subtracting the number of interactions in the rarer direction from those in the more common direction, divided by the total number of interactions. DCI varies on a scale from 0 to 1 (low to high) with a DCI value of 1, meaning there is complete unidirectionality in the outcome of dominant/agonistic interactions across all dyads (van Hooff and Wensing 1987). We also obtained data on rates of agonism per hour of observation time (rate of agonism). This was based on data collected through continuous focal sampling of known females. If more than one value was available for a species (representing different groups/populations), we used the mean of these as the DCI and rate of agonism values for the species as a whole. Data on DCI and rates of agonism were available for 13 and 8 species in the 2D:4D ratio data set respectively (Table II; ESM Table SV).

Mating System Variables

We collated data from various sources on mating system and mating system subcategory for each species (ESM Table SVI). Categories for mating system are monogamy, polyandry, polygyny, and polygynandry. Polygynous and polygynandrous mating systems were then broken down into further subcategories depending on the form in which these systems present in each species (ESM Table SVII; Howlett and Wheeler unpubl. data). As primate mating systems and group size are linked, such that monogamous species often live in smaller groups while the largest groups tend to be associated with polygynandry, we also controlled for group size in the mating system analyses to ensure that any relationships between mating variables and 2D:4D ratios are not a consequence of group size effects.

Intersexual Dominance Patterns

We obtained data on intersexual dominance patterns for all 71 species in the data set from the literature (Table II; ESM Table SVIII). We categorized species as female dominant, male dominant, or codominant. Female dominant species are those in which females are consistently able to dominate and/or receive submissive signals from their male counterparts without reversals in these roles. Male dominant describes species that display the opposite condition. Codominant describes those species in which there is no clear dominance pattern between the sexes or where an alpha male and female pair are codominant over other group members of both sexes.

Statistical Methods

We used Shapiro–Wilk tests throughout when assessing normality of the data. Data for female R2D:4D, L2D:4D, and M2D:4D were normally distributed, as were percentage fruit in diet and rate of agonism. Female body mass and group size were not normally distributed and so these variables were log-transformed. DCI data were also not normally distributed and were arcsine transformed. All other variables are categorical. We tested for multicollinearity of the predictor variables in each analysis using the variance inflation factor (VIF). Variance inflation factors (VIF) of between 1 and 1.344 were obtained indicating there is no evidence for multicollinearity among the independent variables.

2D:4D Ratio Measurement Reliability

We used the single measures values of the intraclass correlation coefficient (ICC) set to the ‘absolute agreement’ definition to test intraobserver reliability. ICC values showed that 2D:4D ratio measurements were highly repeatable for both hands (for R2D:4DICC = 0.992, F512, 7168 = 1891.487, P = <0.001; for L2D:4D ICC = 0.992, F512, 7168 = 1774.794, P = <0.001). We investigated differences in female R2D:4D and L2D:4D using a paired t-test (two-tailed). R2D:4D and L2D:4D in females were not significantly different (t70 = -0.754, P = 0.435) and were tightly correlated (Pearson correlation: r = 0.954, df = 71, P = <0.001); therefore we also used M2D:4D in subsequent analyses.

Phylogenetic Signal

We tested for phylogenetic signal using Pagel’s λ and the packages Devtools and Models of Trait Macroevolution on Trees (motmot) in the statistical software program R version 3.4.1 “Single candle” and 10kTrees (version 3) phylogeny with associated taxonomy from GenBank for the phylogenetic trees (Arnold et al. 2010). We found statistically significant Pagel’s λ values for all variables other than percentage of fruit in diet, DCI, and rate of agonism (Table III). Since 2D:4D ratio measures had significant amounts of phylogenetic signal, the use of phylogenetically controlled methods is appropriate.

Table III Phylogenetic signal results using Pagel’s λ for each variable in female primates

Phylogenetically Controlled Analysis

The 2D:4D ratio is a proposed proxy for PAEs, which in turn are likely modulated by the social and ecological environment. Therefore, in all analyses, we assigned 2D:4D ratio measures as the dependent variables and all other variables as the independent variables. To test for relationships between female 2D:4D ratios and the variables of interest (substrate use, female body mass, percentage of fruit in diet, group size, DCI, rate of agonism, mating system, mating system subcategory, and intersexual dominance pattern), we carried out phylogenetic generalized least squares (PGLS) analysis using the package Comparative Analysis of Phylogenetics and Evolution in R (caper). We ran the analysis across a block of 200 phylogenetic trees for each variable. Molecular data were not available for five species in the 10KTrees phylogeny (Callicebus cupreus, Mico melanurus, Cercocebus chrysogaster, Cercocebus lunulatus, and Propithecus coronatus) and so, using the R package Analysis of Phylogenetics and Evolution (ape), we added these species into the trees based on their relationships to sister taxa that were present in 10Ktrees phylogeny. Callicebus moloch and Callicebus cupreus are sister species that diverged ca. 3.65 mya (Byrne et al. 2016; Perelman et al. 2011). Mico melanurus is in the Mico subgenus and we included it using information on phylogeny from Garbino (2015). We included Cercocebus chrysogaster using Cercocebus torquatus, from which it diverged ca. 3.33 mya. We included Cercocebus lunulatus using its sister species, the Sooty mangabey (C. torquatus atys). C. lunulatus was widely considered a subspecies of C. t. atys (Groves 2001) until recent taxonomic reassessment elevated it to species level (Mittermeier et al. 2013; Oates et al. 2016). Genetic data suggest that P. coronatus and P. deckenii are the same subspecies and so we substituted P. deckenii for P. coronatus in the phylogeny (Pastorini et al. 2001; cf. Thalmann et al. 2002).

Ethical Note

This research was supported by the British and Irish Association of Zoos and Aquariums and all hand images were collected with the permission of the institutions involved. Our method of obtaining digit measurements through digital photographs is noninvasive and Perpex® images were collected only opportunistically when individuals were undergoing veterinary care. The authors declare that they have no conflict of interest.

Data Availability

All data used in this study are available in the Electronic Supplementary Material (Tables SI–SVIII).

Results

Anatomical Considerations

There were no significant relationships between female body mass or substrate use and any of the 2D:4D ratio measures (Table IV). However, due to the variation in hand morphology and allometry that exists between species, female body mass and substrate use were included as factors in all analyses.

Table IV Results of a PGLS regression testing for an effect of substrate use and female body mass on female 2D:4D ratio in primates

Other Control Variables

Both mating variables were significantly related to group size (Table V) and so group size was also included as an additional factor in analysis on mating systems.

Table V Results of the PGLS regression testing for an effect of group size on mating system and mating system subcategory in female primates while controlling for substrate use and female body mass

Diet and Group Size

There were no significant associations between the percentage fruit in diet or group size and any of the 2D:4D ratio measures (Table VI).

Table VI Results of a PGLS regression testing for an effect of percentage (%) of fruit in diet and mean group size on female primate 2D:4D ratio while controlling for substrate use and female body mass

Female Social Variables

There were no significant associations between either DCI (Fig. 2) or rate of agonism and any of the 2D:4D ratio measures (Table VII).

Fig. 2
figure 2

Association between arcsine dominance consistency index (DCI) and female primate 2D:4D ratio measures. Right 2D:4D, gray diamonds, gray solid trendline; left 2D:4D, white triangles, black dotted trendline; mean, 2D:4D, black circles, black solid trendline.

Table VII Results of a PGLS regression testing for an effect of female primate DCI and rate of agonism on female 2D:4D ratio while controlling for substrate use and female body mass

Mating System Variables

There were significant associations between mating system and all female R2D:4D and M2D:4D but not L2D:4D (Table VIII; Fig. 3). There were also significant associations between mating system subcategory and female R2D:4D and M2D:4D but not L2D:4D (Table VIII; Fig. 4).

Table VIII Results of a PGLS regression testing for an effect of mating system variables on female primate 2D:4D ratio while controlling for substrate use, female body mass, and mean group size
Fig. 3
figure 3

Relationship between female primate mean 2D:4D (group mean ± standard deviation) and species mating system. Number of species (S) and number of individuals (I) included in each category are given below the mating system.

Fig. 4
figure 4

Relationship between female primate mean 2D:4D ratio (group mean ± standard deviation) and mating system subcategory. Number of species (S) and number of individuals (I) included in each category are given below the mating system subcategory.

Intersexual Dominance Relationships

There were significant relationships between all 2D:4D ratio measures and intersexual dominance pattern (Table IX). Species in which females are the dominant sex had lower female R2D:4D, L2D:4D and M2D:4D, followed by females in male dominated species with females in codominant species having the highest 2D:4D ratios (Fig. 5).

Table IX Results of a PGLS regression testing for an effect of intersexual dominance relationships on female primate 2D:4D ratio while controlling for substrate use and female body mass
Fig. 5
figure 5

Association between intersexual dominance patterns and female primate mean 2D:4D (group mean ± standard deviation). Number of species (S) and number of individuals (I) included in each category are given below the intersexual dominance pattern.

Discussion

The support for variation in primate socioecology being underpinned by variation in PAEs was mixed. Predictions regarding the relationship between PAEs and female reproductive competition and intersexual dominance patterns were met. However, we did not find the expected relationships between PAEs and degree of frugivory, group size, or the two quantitative measures of female dominance behavior. This study demonstrates that PAEs may act as a proximate mechanism underlying variation in certain aspects of social behavior in female nonhuman primates but may not be associated with the (potential) effects of diet and group size on social structure. Alternatively, these two factors that are often used to formulate the hypotheses on which classic socioecological theory is based may not, in fact, be reliable predictors of female intrasexual competition levels.

As predicted, monogamous females had the highest 2D:4D ratios, likely as a result of selection for many of the key behaviors underpinning monogamous relationships in primates (e.g., pair bonding and parental care) also requiring selection for reduced PAEs and/or increased POEs (French et al. 2018). Polyandrous females had lower 2D:4D ratios than monogamous females but higher 2D:4D ratios than females characterized by polygyny or polygynandry. Females of species characterized by polygynous mating had the lowest 2D:4D ratios, followed by polygynandrous females. It is not clear why 2D:4D ratios would differ in polygynous and polygynandrous female primates but the pattern may stem from differences in female relationships between the subcategories of the mating systems, if indeed these social differences place differential selective pressures in terms of PAEs. Females characterized by scramble competition polygynandry had the lowest 2D:4D ratios but these values are based on two strepsirrhine species only and data from haplorhine females characterized by this mating system are needed before any conclusions can be made as to the relationship between PAEs and this mating system. Female 2D:4D ratios in species characterized by single male/multifemale polygyny and contest competition polygynandry were similar, with single male/multifemale polygyny 2D:4D ratios being very slightly lower. This is perhaps unsurprising, as group formation in species characterized by these mating systems is likely to be very similar, meaning that females have similar competitive relationships with one another. Although small sample sizes precluded formal analysis, among polygynandrous species there was a tendency for female 2D:4D ratios to be highest in species characterized by cooperative defence polygynandry. There is evidence that the formation of platonic social bonds is associated with OT in adult primates (Zeigler and Crockford 2017). The species characterized by cooperative defense polygyny in this study (Pan troglodytes and Alouatta caraya) are likely to experience increased selection for cooperative and affiliative behaviors (moderated by estrogen and its associated neuropeptides) and increased POEs/reduced PAEs, as females in these species disperse from their natal groups at maturity and need to form relationships with other (usually) unrelated females in the groups they migrate into.

There was a significant relationship between female 2D:4D ratios and intersexual dominance patterns. Female 2D:4D ratios were lower in species characterized by female dominance than those characterized by both male dominance and codominance. The low 2D:4D ratios in species characterized by female dominance may be due to selection for higher PAEs that could result in the behavioral masculinization necessary for dominating males in these species (Dunham 2008). The results lend support to prenatal sex hormones acting as a proximate mechanism underlying the evolution of female dominance in nonhuman primates (Petty and Drea 2015). Codominant species have the highest 2D:4D ratios (lower PAEs) and this could be due to their mating systems, as species in which the sexes are codominant (e.g., gibbons, titi monkeys, marmosets, and tamarins) are often characterized by monogamy or polyandry (Carpenter 1940; Kinzey 1981; Koba et al. 2012).

There were no significant relationships between female 2D:4D ratios and either of the measures of female dominance behavior. Species that showed higher rates of female–female agonism did not have significantly lower 2D:4D ratios than those with lower rates. Likewise, the associations between DCI and the 2D:4D ratio measures were not significant, although a weak trend was discernible in the predicted direction (i.e., females with higher DCI tended to have lower 2D:4D ratios [higher inferred PAEs]), although this was limited to the right hand. As DCI is perhaps the most reliable way to quantify the degree of despotism characterizing a given primate group, this trend is consistent with previous studies that noted associations between levels of female intrasexual competition and 2D:4D ratios in anthropoid primates (Nelson and Shultz 2010). Quantitative measures of female dominance are better able to represent the variation we see in female social relationships across primates than the broad classifications of traditional socioecological models. However, sample sizes were small for both of these variables and so the lack of a significant relationship may result from a lack of statistical power. Further research into the relationship between quantitative female dominance measures and 2D:4D ratios with larger sample sizes could prove productive.

Contrary to our prediction, females with a higher percentage of fruit in their diet did not have lower 2D:4D ratios. The hypothesized greater contestability of fruit relative to other food types led to the suggestion that feeding competition should be greatest among highly frugivorous species and lowest among highly folivorous species (Clutton-Brock and Janson 2012; Snaith and Chapman 2007), which has been corroborated by some studies (Janson and Goldsmith 1995; Steenbeek and van Schaik 2001), but not others. For example, intragroup contest competition was not elevated by increased feeding on fruits in Macaca assamensis (Heesen et al. 2013) and no association was found between rates of agonism in female primates and the percentage of fruit in their diet (Wheeler et al. 2013). Folivorous species display a preference for young, tender leaves that can be patchily distributed and therefore contestable (Koenig et al. 1998; Yeager and Kool 2000). The fact that folivores may also experience high levels of contest competition is supported by the observation that female dominance relationships in a highly folivorous colobine monkey were substantially affected by occasionally feeding on seasonal high-quality, contestable foods (Wikberg et al. 2013). The results lend support to dietary categories (or degree of frugivory) not being good predictors of female intrasexual competition in primates. Other resources may be more important for some species, such as reproduction (French and Inglett 1989), and territory (Schülke and Kappeler 2003) and competition between females may occur in reference to these rather than just food. Although these findings were not significant in the Order-wide analysis, the relationship between PAEs and variables relating to competition for food may be apparent between individual females at the species level via the positive impact that PAEs have on an individual’s competitive abilities and dominance rank (Hönekopp et al. 2006; Howlett et al. 2012). For example, in certain conditions, female primates with lower 2D:4D ratios (high PAEs) may be better able to outcompete intragroup rivals with higher 2D:4D ratios (lower PAEs) due to the behavioral predispositions brought about by high PAEs. This is supported by the link between dominance rank and resource acquisition in female primates including olive baboons (Papio anubis: Barton and Whiten 1993) and sooty mangabeys (Cercocebus torquatus atys: Range and Noë 2002).

Females of species characterized by larger group sizes did not have lower 2D:4D ratios than females of species that live in smaller groups. This is somewhat surprising as females living in larger groups have a greater number of intragroup competitors and an increased likelihood of experiencing contest competition (Koenig and Borries 2006) that can result in the higher rates of agonism that are observed among females living in larger groups (Wheeler et al. 2013). A possible explanation for the lack of relationship could be that intergroup competition in smaller primate groups, for example over territory in pair-living and solitary species, also generates selection for higher levels of PAEs, particularly in species where females are actively involved in territory defence (Anzenberger 1992; Raemaekers and Raemaekers 1985). Additionally, there is likely to be substantial selection for cooperative and affiliative behaviors to facilitate the coexistence of primates in large groups and this is likely to moderate selection for higher PAEs (French et al. 2018). Further, depending on relationships between females, larger group sizes could result in increased intragroup scramble competition and, unlike direct contest competition, scramble competition is not likely to result in increased selection for behaviors associated with higher PAEs (e.g., dominance), although there is evidence to suggest that intragroup contest competition also increases with increasing group size (Koenig and Borries 2006; Wheeler et al. 2013). This, together with the fact that dietary composition and group size cannot fully predict female intrasexual competition within and between groups, likely explains the lack of a relationship between these variables and PAEs.

The hormones underlying competitive and cooperative behaviors are largely conserved across all mammals (Soares et al. 2010; Trumble et al. 2015) and PAEs may act as a proximate mechanism underlying the expression of certain aspects of behavior in female primates in ways that are adaptive to their social system. The results stress that, to fully understand the relationship between PAEs and competition in female nonhuman primates, it is necessary to consider reproductive competition and intersexual dominance relationships alongside ecological factors and female sociality. The predictions of this study are based on the assumption that the 2D:4D ratio is a negative correlate of PAEs across primate species, such that lower 2D:4D ratios indicate higher PAEs and higher 2D:4D ratios indicate lower PAEs. However, because some studies have observed species in which males have on average higher 2D:4D ratios (lower inferred PAEs) than females (Baxter et al. 2018; Roney et al. 2004) this suggests that the pattern may not be universal across the Order and in some species for example, higher 2D:4D ratios may indicate higher PAEs. Such reversals in the pattern would bias against our predictions, and this, and not a lack of PAEs, may explain nonsignificant results. However, the fact that tests based on the broader categories of mating systems and intersexual dominance patterns were significant suggests that the expected trend of lower 2D:4D ratios reflecting higher PAEs likely holds across most species in the data set. Additional research into species specific patterns in 2D:4D ratio with consistent, reliable, and comparable measurement methods is needed to investigate this further.