- Research article
- Open Access
Does life history shape sexual size dimorphism in anurans? A comparative analysis
BMC Evolutionary Biology volume 13, Article number: 27 (2013)
The evolution of sexual size dimorphism (SSD) is likely constrained by life history. Using phylogenetic comparative methods, we examined correlations between SSD among anurans and their life history traits, including egg size, clutch size, mating combat, and parental care behaviour. We used sexual dimorphism index (SDI = Body-sizefemale /Body-sizemale –1) as the measurement for SSD. Body size, life history and phylogenetic data were collected from published literature. Data were analysed at two levels: all anuran species and within individual families.
Female-biased SSD is the predominant form in anurans. SSD decreases along with the body size increase, following the prediction of Rensch’s rule, but the magnitude of decrease is very small. More importantly, female body size is positively correlated with both fecundity related traits, egg size and clutch size, and SDI is also positively correlated with clutch size, suggesting fecundity advantage may have driven the evolution of female body size and consequently leads to the evolution of female-biased SSD. Furthermore, the presence of parental care, male parental care in particular, is negatively correlated with SDI, indicating that species with parental care tend to have a smaller SDI. A negative correlation between clutch size and parental care further suggests that parental care likely reduces the fecundity selection pressure on female body size. On the other hand, there is a general lack of significant correlation between SDI and the presence of male combat behaviour, which is surprising and contradictory to previous studies.
We find clear evidence to support the ‘fecundity advantage hypothesis’ and the ‘parental care hypothesis’ in shaping SSD in anurans. Nevertheless, the relationships of both parental care and combat behaviour to the evolution of SSD are complex in anurans and the extreme diversity of life history traits may have masked some potential interesting relationships. Our study represents the most comprehensive study of SSD in anurans to date.
Sexual size dimorphism (SSD), where males and females differ in body size, is the most conspicuous difference between the two sexes in many animals. Females are larger than males in most animals, while the opposite pattern is predominant in birds and mammals [1, 2]. SSD has been hypothesized as adaptation of males and females to their disparate reproductive roles and associated differences in ecology and life history strategies [2, 3]. A widely accepted explanation of female-biased SSD is that selection for fecundity favors larger females, because they can produce more offspring or can have larger energy storage for reproduction (“the fecundity advantage hypothesis” ). In contrast, a male-biased SSD is expected when natural selection for female fecundity is less intense than sexual selection on male body size. In this scenario, large size gives advantages to males in contests for females or for female choice of larger males . Selective advantages may bias SSD to one way or another; and the extent of SSD is likely constrained by the shared genomes of the two sexes . A widespread SSD pattern in the animal kingdom is that the magnitude of SSD decreases with body size in groups where females are the larger sex, but increases with body size in groups that males are the larger sex (Rensch’s rule) .
Selection on many aspects of life history undoubtedly plays a major role in determining direction and magnitude of SSD, and studies on many animal groups have attempted to explain SSD by their life history traits. Fecundity-related life history traits, such as egg size and clutch size, are probably the most evident [1, 5], but others are often controversial . For example, parental care is likely associated with the body size of the care-providing parent, and it may lead to dimorphism between two sexes. In some cases, selection may favor larger parents because they provide more resources to offspring (e.g.). In other cases, parental care activities may result in the small body size of the care-providing parents because the energy expenditure on breeding activities may constrain growth (e.g.).
We know very little about the impact of life history traits on the evolution of SSD in anuran amphibians, despite the fact that they have a tremendous diversity of life history traits [8–10]. Shine (1979) presented the most comprehensive study to date that tested adaptive hypotheses of SSD in relation to mating-related life history traits . Assessing 589 anuran species, Shine concluded that female-biased SSD is the common pattern in anurans with only approximately 10% species showing male-biased SSD. Furthermore, male-biased SSD was strongly correlated with the presence of male combat behaviour and the appearance of male weaponry, such as spines and tusks . Several recent case studies supported this “male combat hypothesis”. For example, in two distantly related species, the tusked frog (Adelotus brevis) and the fanged frog (Limnonectes kuhlii), males are larger than conspecific females and have larger paired projections in their low jaw [12, 13]. Males defend their calling sites by attacking rival males using their fangs or tusks; hence larger body size and larger fangs or tusks may provide males with a better chance of winning battles and enhancing their reproductive success . Zheng et al.  also found a positive correlation between the presence of weapon-like keratinized maxillary spines and male-biased SSD in megophryid frogs of the genus Leptobrachium. Hudson et al. confirmed that the Emei mustache toads (L. boringii) indeed use the spines (“moustache”) as weaponry to fight their rivals. Other life history traits, such as parental care, length of breeding/growth season, developmental rate/time and ages of sexual maturity have also been explored in anurans [16–19].
Several phylogenetic comparative methods (PCMs) have been developed [20, 21] since Shine’s study in 1979 , which provide better tools to analyse correlated evolution.Closely related species share phenotypic similarities inherited from a common ancestor, so direct correlation studies that treat each species as an independent data point tend to overestimate the degrees of freedom for statistical tests. Differing from traditional cross-species analysis, PCMs correct this statistical non-independence within a phylogenetic framework. The advent of DNA sequence analysis and the development of phylogenetic methods in the last few decades have also led to an explosion of new phylogenies of many major organism groups. Anurans are no exception; several large-scale phylogenies have been reconstructed (e.g.[22, 23]). Additionally, tremendous amounts of data on body size and life history have been accumulated for many anuran species since Shine’s study. Together, these new tools and new data provide an excellent opportunity to utilize phylogenetically corrected analysis to re-examine the role of life history in shaping SSD in anurans.
We attempt a comprehensive evaluation of SSD in all anurans with two objectives. First, we assess the overall patterns of SSD and test for Rensch’s rule in anurans. Second, we investigate correlated evolution between SSD and life history traits, particularly traits that are related to fecundity, combat, and parental care behaviour. We used size dimorphism index (SDI; Body-sizefemale/Body-sizemale -1)  to represent SSD, and all tests were controlled for evolutionary relatedness among species using PCMs. We made several predictions based on current understanding of SSD evolution and observations in anurans, assuming that inter-specific differences reflect a history of selection at an intra-specific level: 1) Egg size and clutch size are positively correlated with female body size and SDI (fecundity advantage hypothesis). 2) The presence of male mating combat behaviour is positively correlated with male body size and negatively correlated with SDI (male combat hypothesis). 3) Rensch’s rule predicts a negative correlation between the ratio of female to male body size and the mean body size of a species for a group with female-biased SSD as the predominate form. How parental care is associated with SSD is uncertain; it may be positively or negatively correlated with the body size of the care providing sex and hence SSD.
Body size, life history, and phylogenetic data were collected from published literature. For behavioural data, only species with sufficiently detailed descriptions were used. Thus, many species in Shine’s study  were not included. In total, body size and life history data were collected for 688 anuran species, including 534 species with body size data and at least one life history trait. For the remaining 154 species, we did not have body size information. Among them, 119 species had egg size and clutch size data and were used for correlation analysis; 36 species had only mating combat and/or parental care information and they were not used for correlation analysis but used for detecting phylogenetic signals of these traits. In total, the 688 species represented approximately 10% of the described extant anuran species . The majority of the data came from six groups: the superfamily Dendrobatoidea (family Aromobatidae and family Dendrobatidae; 56 species), the families Bufonidae (70 species), Dicroglossidae (37 species), Hylidae (191 species), Megophryidae (36 species), and Ranidae (115 species). The remaining 183 species were scattered across 35 other families. Species names and family names followed the “Amphibian Species of the World version 5.5” website as of March 2012 . All body size and life history data used in this study are presented in Additional file 1, and their sources are provided in Additional file 2.
Body size and life history traits
Snout-vent length (SVL) was used to represent anuran body size. Different publications provided different levels of detail for body size measurement, and consequently, the body size data were separated into two categories: mean body size (‘mean’) and body size range (‘range’). ‘Mean’ referred to body size data given as a mean SVL, and it required measurements with a minimum sample size of three for each sex per species. ‘Range’ included body size provided as maximum and minimum values, and also included body size measured on only one or two individuals in each sex. If more than one study provided ‘mean’ data, an unweighted average across studies was calculated for each sex. This criterion was also applied for other life history traits that were reported in multiple sources. Body size of a species was represented by the mean of both sexes.
All analyses pertaining to body size and sexual size dimorphism were conducted twice: once with only the high quality data, the ‘mean’ data, and once with pooled ‘mean’ and ‘range’ data. For ‘range’ data, the midpoint of the range was used to represent the body size in each sex. Although the ‘range’ data are not as accurate as the ‘mean’ data, they represent an additional 154 species (Table 1), and sample size is important for statistical tests (e.g. AIC scores).
Sexual size dimorphism
SSD was treated as a continuous trait using size dimorphism index (SDI) . SDI was estimated by taking the ratio of female to male body size and subtracting one (Body-sizefemale/Body-sizemale -1). The index has been widely used in many previous studies because of its intuitive appeal for both direction and degree of SSD as well as other advantages .
Egg size and clutch size
Egg size was represented by the diameter of the ovum. Clutch size referred to the number of eggs a female laid in a clutch. When females produced multiple clutches in a year, we only used the reported size of a single clutch for that species. The mean values of egg size and clutch size were preferentially used, but the median of a range was accepted when a mean value was not available. Egg size and clutch size are potentially correlated with female body size and SDI, and therefore, we considered both as explanatory variables in a linear multiple regression analysis. In the full model of the regression, we considered the effects of both traits and their interaction.
The presence of mating combat was defined as observing aggressive physical contests between individuals of the same sex to gain access to mates or attractive territories . The absence of a record of mating combat may result from lack of observation. Therefore, only species for which other mating behaviours had been observed with no mention of combat were considered as exhibiting no form of mating combat. Such mating information includes descriptions of nest building, mate searching, advertisement call, male satellite behaviour, and breeding aggregation. In contrast to Shine , tusks and spines on male’s body were not viewed as signals of combat potential, as they might be primarily developed and generally used for other purposes, such as aids in amplexus, anti-predation or foraging adaptations . Both male and female anurans are known to partake in combat , therefore, each sex was analysed separately. Species, in which both males and females exhibited combat behaviour, were included in both male and female combat analyses.
Parental care, defined as any form of post-ovipositional parental investment that increases the survival of offspring , was reviewed based on descriptions in the literature of parental attendance of clutches, building of foam nests when laying eggs, carrying eggs, transporting tadpoles, larval development in or on parents’ bodies, and female deposition of unfertilized eggs to feed tadpoles . The absence of parental care was defined as no evidence for any of these behaviours when other breeding information was available. This breeding information includes the descriptions of breeding habitat, egg development habitat, breeding behaviour during egg-laying, egg developmental time, egg size, and clutch size. Parental care in the two sexes was analysed separately. Species with both parents exhibiting parental care were included in both analyses.
Phylogenetic comparative analysis
We used two PCMs to remove phylogenetic autocorrelation between species: multiple regression with a phylogenetic generalized least squares (PGLS) model [28, 29] and phylogenetic independent contrasts (PIC) [20, 30]. All analyses were conducted using the ‘Comparative analysis of phylogenetics and evolution package in R (caper)’ . PCMs require a phylogenetic tree, and we used Pyron & Wiens’ (2011) tree of amphibians , which represents the most recent and comprehensive phylogeny of anurans. Anuran classification is undergoing major revisions (e.g.), and some formerly recognized species have been split into several new species. In these cases, we matched the source populations of life history data with the species in the phylogeny.
Anuran families differ dramatically in mating systems and reproductive strategies, and so the correlated evolution of certain life history traits may be specific to certain families. Accordingly, most analyses were conducted at two levels: across all anuran species (all-anuran analysis) and within each of the six families or superfamily (within-family analysis). Species in the superfamily Dendrobatoidea, including the families Dendrobatidae and Aromobatidae, have many life history similarities [8, 32], therefore, they were analysed together. It was noteworthy that, for each correlation analysis, only species with data for both traits were used; species with missing data were automatically pruned from the analyses in ‘caper’. Consequently, the sample size in each correlation analysis on given traits was always smaller than the number of species in a group. Since each comparison had a different composition of species, a Bonferroni correction could not be applied .
Testing phylogenetic signals of traits
We used the D statistic  and a phylogenetic parameter λ as the measures of phylogenetic signals of life history traits. D is applicable only for binary traits, and was carried out with the ‘phylo.d’ function in ‘caper’. D typically varies between 0 and 1. A D of 0 indicates that a trait evolves on a tree following the Brownian model (strong phylogenetic signal), and a D of 1 indicates that a trait evolves following a random model (no phylogenetic signal). D can be negative, which means that a trait evolves in a conserved way: more conserved than predicted by the Brownian model. Additionally, we conducted a simulation (1000 permutations) to test whether an estimated D was significantly different from the predictions of a random or a Brownian style evolution.
The parameter λ is applicable for both continuous and binary traits, which reveal the dependence among species for a given phylogeny and a given trait. Its optimum value and confidence limits were estimated using the ‘pgls’ function in ‘caper’ with a maximum likelihood method when performing a correlation analysis. λ shows the strength of the phylogenetic signal of a phylogenetic generalized linear model, which ranges between 0 and 1. When λ=0, the variation of a trait is modeled as a function of an independent evolution along the branches leading to the tips. When λ=1, the trait shows variation expected under the Brownian model. The estimated λ was used not only for measuring strength of phylogenetic signal, but also for transformation of internal branch lengths for phylogenetic linear model analysis.
Multiple regression analysis using a phylogenetic generalized least squares model
PGLS includes a phylogenetic tree as a covariance matrix in a linear model . It is capable of evaluating multiple causative variables simultaneously and incorporating polytomies .
We first used a model selection approach to test the relative importance of all life history traits and their interactions. The SDI was set as the response variable, and initially egg size, clutch size, male combat, and parental care were included as explanatory variables in the model. Other measured variables either overlap with male combat (i.e. male scramble competition and male territory defence) or with parental care (i.e. male and female parental care), so they are not independent predictors and cannot be used in the same model. To further test the combat hypothesis, female combat and the two different forms of male combat were analysed separately by replacing male combat with each of them. Similarly, female parental care and male parental care were analysed separately by replacing parental care with each of them to further explore the parental care hypothesis. From the saturated full model, including four explanatory variables and their interactions, we excluded the non-significant interactions and variables step by step. The reduced model was then compared to the unreduced model using the ‘anova’ function, and the difference in AIC scores was used to determine the significance of the exclusion. Model comparison requires complete data for each species. Consequently, approximately 68% of the species were excluded from this analysis because of missing data in at least one of the variables. This dataset with no missing values is denoted as the ‘reduced dataset’. The analysis was conducted for all-anurans and for each of the six (super)families, and was carried out using the ‘pgls’ function in ‘caper’.
Branch lengths were transformed by three scaling parameters, λ, κ, and σ in order to make the data fit the Brownian model of evolution. Although all of them could be optimized by ‘caper’, altering three parameters simultaneously would make the biological interpretation difficult . Consequently, we only estimated the parameter λ, because it was significantly different from its upper and lower bound values (0 and 1) for all traits. The other two parameters were kept constant with their default values. The λs for ‘range + mean’ data and ‘mean’ data were estimated separately because of different accuracy in body size measurements. The λs used in the family-level analysis were the same λs estimated for all-anurans, because they were based on more data and thus were more likely to be accurate.
We further tested the correlation between SDI and each of the explanatory variables independently (egg size and clutch size were analysed in combination), as well as between several pairs of life history traits. Compared with the model selection approach, this simple linear regression analysis only excluded species with missing values of the focal traits, and therefore, utilized the most available data and likely had higher statistical power to detect a correlation. The full dataset used in this set of analyses is denoted as the ‘extended dataset’. It was noteworthy that since several traits (e.g. SDI) were repeatedly used in several linear regressions, the Type I errors of this set of analyses were inevitably inflated.
Comparative analysis with phylogenetic independent contrasts
PIC only allows pairwise correlation analysis. We used the ‘crunch’ function in ‘caper’ to calculate standardized contrasts of all continuous traits. This function then used a simple linear model to test for a correlation between contrasts of two continuous traits. We used the ‘brunch’ function to test for correlation between a continuous trait and a categorical trait. ‘Brunch’ estimates independent contrast values of a continuous trait on the nodes where the state of the categorical variable changes . A one-sample t-test was then used to determine if the mean of the independent contrasts of the continuous response variable was significantly different from zero. After running ‘crunch’ and ‘brunch’, the ‘caic.diagnostics’ function in ‘caper’ carried out diagnostic tests for the robustness of contrasts and identified the outlier contrasts that had absolute studentised residuals greater than three. The outlier contrasts might bias the correlation analysis, and so were filtered out using the ‘caic.robust’ function before the final analyses . Ultimately, all diagnostic tests of our correlation analyses were insignificant, indicating our data met the assumption of evolution under a Brownian model.
Testing Rensch’s rule
Rensch’s rule predicts that the magnitude of SSD decreases along with the increase of body size in anurans where females are generally the larger sex . This pattern was detected as a slope significantly smaller than one when male body size was on the x-axis in a major-axis regression between the independent contrasts of logarithm transformed male and female body size . A major-axis regression through the origin was conducted in R .
General pattern of SSD in anurans
This study presents the most comprehensive study of SSD in anurans to date. Of the 688 species that we could gather data for, body size data were available for 534 species, within which 380 had ‘mean’ data, and 154 had ‘range’ data. The sample size for each trait was different; all traits and their sample sizes are presented in Table 1. The distributions of all continuous traits (SDI) or their logarithm-transformed forms (body size, egg size, and clutch size) were close to normal. Figure 1 presents the distribution of SDI and log body size.
Body size (SVL) had an average of 47.92 mm for the ‘range + mean’ data and 49.06 mm for the ‘mean’ data (Figure 1), and ranged from 13.35 mm in the blue-bellied poison frog (Ranitomeya minutus) to 260.00 mm in the goliath frog (Conraua goliath) (Additional file 1). The majority of species (97%) had body sizes smaller than 100 mm.
Sexual dimorphism in anurans was predominantly biased towards females (Figure 1). The distribution of SDI was close to normal, with an average of 0.163 for the ‘range + mean’ data and 0.156 for the ‘mean’ data. Male size was found to significantly exceed the size of conspecific females in only nine species, although 58 species (10.9%) with ‘range + mean’ data and 48 species (12.6%) with ‘mean’ data had negative SDIs. The magnitude of male-biased SSD was similar to the magnitude of female-biased SSD. Females could be almost twice the size of conspecific males (e.g. 188% in the odorous frog Odorrana schmackeri; Additional file 1), while the African bullfrog (Pyxicephalus adspersus), which demonstrated an extreme male-biased anuran SSD, had male sizes 176% that of conspecific females .
Species with negative SDI (male-biased SSD) were present in 19 anuran families (Aromobatidae, Bufonidae, Bombinatoridae, Centrolenidae, Cryptobatrachidae, Dendrobatidae, Dicroglossidae, Hylidae, Leiuperidae, Leptodactylidae, Limnodynastidae, Mantellidae, Megophryidae, Microhylidae, Myobatrachidae, Pelobatidae, Petropedetidae, Pipidae, Pyxicephalidae, Ranidae, Rhacophoridae, and Scaphiopodidae) and distributed across six continents (Asia, Africa, Australia, Europe, North America, and South America) in both temperate and tropical zones. Some of these species were not included in the final analysis because of a lack of other life history information.
The impact of phylogeny on the evolution of life history traits
The D statistics of all but one of the binary traits fell between 0 and 1 (Additional file 3). The evolution of all traits was significantly different from random evolution, and the evolution of three binary traits (female combat, parental care, and male parental care) did not differ significantly from the Brownian model (p > 0.05). D of ‘male parental care’ was < 0, indicating the evolution of this trait was more conservative than predicted by the Brownian evolution model. Lambdas (λ) were between 0 and 1 for all models, and were significantly different from the two bounds, indicating that the evolution of all the traits depended on their phylogenetic history (not random), but might not strictly follow the Brownian model.
Clearly, the evolution of all traits was phylogenetically constrained at various degrees. Therefore, including the phylogenetic perspective in the comparative analysis was essential. In the meantime, the evolutionary history of these traits did not strictly follow the Brownian evolution, on which all our analyses were based, and therefore, branch length transformation was necessary to best-fit the evolution of these traits to the Brownian model.
Correlations between SSD, body size, and life history traits
Comparative analysis using PGLS
Based on the ‘reduced dataset’, the model selection analysis for all-anurans found that the best model was SDI to be negatively correlated with parental care, and with male parental care in particular (Table 2). Parental care was the only significant, and the most influential, variable that might explain SDI (p < 0.01). Both ‘range + mean’ data and ‘mean’ data supported this correlation. Family level analysis revealed several correlative relationships within specific families (Table 2): a positive correlation between SDI and egg size in the families Bufonidae and Ranidae, a negative correlation between SDI and parental care in the family Hylidae, a positive correlation between SDI and clutch size and negative correlations between SDI and female parental care in the superfamily Dendrobatoidea, and a negative correlation between SDI and parental care in the family Hylidae.
Simple linear regression analysis based on the ‘extended dataset’ of all-anurans again showed that SDI was negatively correlated with parental care (p < 0.01), particularly with male parental care (p < 0.01; Table 3). In addition, SDI was negatively correlated with egg size (p = 0.05). Furthermore, female body size was positively correlated with both egg size and clutch size (p < 0.01), and parental care, particularly female parental care, was negatively correlated with clutch size (p < 0.01). Nevertheless, none of the combat behaviour was correlated with SDI or body size, except for a marginal correlation between SDI and male territory defence (p = 0.09; Table 3). More detailed results are provided in Additional file 4.
Overall, the presence of parental care behaviour was clearly associated with SDI and fecundity, but its correlation with body size was less evident. The fecundity traits were correlated with SDI and female body size. On the other hand, considering previous findings in anurans (e.g. ) and patterns of SSD in other animal groups (e.g. mammals and birds ), the lack of significant correlation between male combat behaviour and SSD was surprising.
Comparative analysis using PIC
PIC analysis produced similar results to that of PGLS. For example, it detected a negative correlation between SDI and parental care, male parental care in particular, and a positive correlation between female body size and both egg size and clutch size in all-anuran analysis. The two analyses also agreed on several correlations at the family level, e.g. clutch size was positively correlated with SDI in the superfamily Dendrobatoidea, and parental care was negatively correlated with SDI in the family Hylidae.
PIC analysis revealed several additional significant correlative relationships within individual families (Additional file 4). In Dendrobatoidea, male body size and male combat behaviour were positively correlated (p = 0.04), and female body size and female parental care were negative correlated (p < 0.05). In Hylidae, SDI was negatively correlated with female parental care (p < 0.05). Additionally, PIC had several minor disagreements with PGLS on the importance of different fecundity traits. For example, PGLS suggested that SDI was positively correlated with egg size in Bufonidae, but PIC suggested a positive correlation with clutch size. Complete results of PIC analysis for both all-anurans and individual families are presented in Additional file 4.
The independent contrasts of log female body size and log male body size were significantly correlated with each other in the all-anuran analysis and in every within-family analysis (p < 0.01, Figure 2). The slopes were significantly lower than one in the all-anuran analysis and in the family Megophryidae, which was consistent with the prediction of Rensch’s rule (Figure 2). The slope and its upper confidence limit in the all-anuran analysis, however, were close to one (slope = 0.96, upper confidence limit = 0.99), and therefore, the negative correlation between body size and SSD was statistically significant but biologically weak in anurans. With an increase of body size, male body size increased only slightly faster than female body size. Slope of one was within the confidence intervals of all other families.
Sexual size dimorphism in anurans
Female-biased SSD is the predominant pattern in anurans (Figure 1). This is congruent with Shine’s early assessment . Deviation or reversal from this common pattern has independently evolved many times and its occurrence is not restricted to certain lineages or in certain geographic regions. SSD is relatively small in magnitude, and in the species with the largest SSD, females do not exceed twice the size (snout-vent length) of conspecific males. Male-biased SSD in anurans, however, are among the most extreme cases. Male African bullfrogs exceed 1.7 times the size of conspecific females, which is comparable to the most extreme known case in animals, in which the males are twice the size of females (cichlid fish Lamprologus callipterus) .
The SSD in anurans shows a weak trend that is consistent with Rensch’s rule. The trend is statistically significant but biologically weak in anurans as a whole, and it is apparent only in the family Megophryidae (Figure 2). The weak, or lack of, support for Rensch’s rule in anurans is not particularly surprising. Supports for the rule are mostly seen in groups that have a male-biased SSD , and exceptions to the rule are commonplace . Anurans have a high diversity of life history traits, which may also contribute to the weak support for Rensch’s rule, as many of these traits may alter the general pattern of body size dimorphism.
Correlated life history traits
The negative correlation between SDI and parental care provides clear evidence for the ‘parental care hypothesis’ that parental care directly contributes to the patterns of SSD (Tables 2 and 3). The idea that male parental care may drive the evolution of male-biased SSD in anurans has recently been suggested several times, but has never been formally tested (e.g.[16, 42]). Similar to fish , but unlike most other vertebrate groups, male anurans are the primary care providers (Table 1). The correlation between SDI and parental care may be mostly contributed by the species with male care in anurans, as a significant correlation between SDI and male parental care was found in all-anurans. Nevertheless, significant negative correlation between SDI and female parental care was also found in the (super) families Dendrobatoidea (PGLS) and Hylidae (PIC).
The mechanisms linking parental care and SSD remain to be tested. Several hypotheses have been proposed. First, parental care may improve the survival of offspring and hence ‘relax’ fecundity selection on female body size, which may reduce the magnitude of female-biased SSD or reverse its direction. Parental care and clutch size were negatively correlated in anurans (Table 3 & Additional file 4), which is consistent with the prediction of this ‘relaxed’ fecundity selection scenario. Summers et al. also reported a significant positive correlation between parental care and egg size in anurans, although we only detected a marginal correlation between male parental care and egg size (Table 3). This correlation has also been documented in many other vertebrate groups, e.g. fishes . Furthermore, male-biased SSD is often present in lineages where females do not produce large numbers of eggs, such as those in the families Bombinatoridae, Megophryidae, the superfamily Dendrobatoidea, and some clades in the family Hylidae (Additional file 1). Second, some forms of parental care, such as protecting offspring from predation, may demand strength or size of the parents [8, 10]. Males in some species with male-biased SSD, such as the African bullfrog (Pyxicephalus adspersus), were observed to defend tadpoles against vertebrate predators . Third, the care-giving parents may benefit from large amount of fat storage in their body, and hence large body size. Looking after young may require additional energy consumption and/or a reduction in energetic intake . Case studies also demonstrate significant body weight loss during the breeding period in Emei mustache toads (Leptobrachium boringii) . We did not find support for the last two hypotheses because of the lack of any positive correlation between parental care and body size in both genders. Furthermore, the negative correlations between female body size and female parental care (PIC; Additional file 4) in Dendrobatoidea, and between SDI and female parental care in both Dendrobatoidea (PGLS; Table 2) and Hylidae (PIC; Additional file 4) are consistent with the first hypothesis, but not with the last two. Therefore, our data support the relaxed fecundity hypothesis, although it does not necessarily reject the other two alternatives.
Fecundity advantage appears to select for large females in anurans, and considering that nearly 90% of all anurans have female-biased SSD, it certainly is a leading selecting force in shaping SSD in this group. The increase of female body size is strongly linked to the increase of egg size and/or clutch size in all-anurans as well as in individual families (Table 3 & Additional file 4). This pattern is congruent with observations from the majority of animal groups , and is consistent with the fecundity advantage hypothesis  (but see ). Furthermore, significant correlation between SDI and the fecundity-related traits (egg size and/or clutch size) was found in all-anurans and several families (Additional file 4). The positive correlations between these fecundity-related traits with SDI and female body size provide support for fecundity driven SSD evolution. Female body size may have driven the evolution of female-biased SSD through selection on egg size and clutch size. On the other hand, the relaxation of selection on female fecundity and accordingly on female body size, may result in male-biased SSD.
Surprisingly, our data do not support the ‘male combat hypothesis’. We did not detect any significant correlations between SDI and any male combat behavioural traits. Significant correlation between male body size and male combat behaviour was detected only in the superfamily Dentrobatoidea (PIC, Additional file 4). This weak, or lack of, correlation is contradictory to the well-established concept that male-male competition (sexual selection) drives to increase male body size and leads to male-biased SSD. Within anurans, Shine’s early evaluation and several case studies also associated the combat behaviour with male-biased SSD [11–13, 15]. This apparent contradiction is unlikely the result of small sample sizes, and hence low statistical power in our analysis. We have 599 species with combat (or lack of combat) information (Table 1). In the all-anuran analysis with PIC, the number of contrasts between SDI and male combat is 73 (‘range + mean’ data) and 51 (‘mean’ data only), respectively (Additional file 4). We are convinced that the lack of support for the ‘male combat hypothesis’ in anurans is real and there are several potential causes. First, male body size may influence mating success in some anuran species, but may have no impact in others . Second, the common occurrence of alternative mating tactics in anurans, such as satellite male behaviour , may reduce the selection pressure on male body size. Small males may gain fitness by using alternative mating tactics. For example, while large male bullfrogs (Rana catesbeiana) are territory defenders, small males are often satellites, and satellite males may account for as much as 20% of the total mating in a population . Third, the extreme diversity of mating systems, life history traits and their plasticity in anurans may mask the potential contribution of male-male competition to body size evolution. The few case studies that clearly link them together likely represent extreme cases, as a consequence of researchers sampling bias. For example, in all three aforementioned case studies [12, 13, 15], males have weaponry structure and are significantly larger than conspecific females. These characters make them more attractive to be study objects than others.
It is worthwhile to note that the three alternative hypotheses, parental care, fecundity advantage, and male combat, are not mutually exclusive. In nature, different selection forces may complement or oppose each other, and their interactions may vary among taxa. Therefore, the sparse or lack of significant correlation in many cases may reflect the true complexity of various selection forces that are involved in the evolution of body size and SSD.
Sexual size dimorphism in anurans is predominantly female-biased and fecundity advantage may play a significant role in determining this pattern. Paternal care in anurans may lead to deviations from this common pattern and contribute to the evolution of male-biased SSD. This discovery provides a new insight in understanding the evolution of SSD: natural selection, as well as sexual selection, can drive male-biased SSD. More research in fish, amphibians and invertebrates are needed to understand the general patterns of SSD in animals.
A large number of detailed case studies may be necessary to understand the alternative functional explanations of SSD. Correlations are not causations, but well-established correlations provide guidance for future experimental work, and only the latter can unambiguously establish causation relationships. Species in the superfamily Dendrobatoidea, and the families Dicroglossidae and Megophryidae exhibit both male-biased and female-biased SSD in a large number of clades (e.g.). They may provide a fertile ground for both family level comparative studies and detailed case studies.
Availability of supporting data
The data supporting the results of this article are included within the article and its additional files.
XH is currently a PH.D candidate at the Queen’s University and her primary interests are the evolutionary consequences of sperm senescence in Drosophila melanogaster. JF is a herpetologists and is interested in a broad range of evolutionary questions.
Akaike information criterion
Phylogenetic comparative method
Phylogenetic generalized least squares model
Phylogenetic independent contrasts
Sexual dimorphism index
Sexual size dimorphism
Andersson M: Sexual Selection. 1994, Princeton University Press, Princeton
Fairbairn DJ, Blanckenhorn WU, Szekely T: Sex, Size and Gender Roles: Evolutionary Studies of Sexual Size Dimorphism. 2007, Oxford University Press, New York
Shine R: Ecological causes for the evolution of sexual dimorphism: a review of the evidence. Q Rev Biol. 1989, 64: 419-461.
Rensch B: Die Abhangigkeit der relativen Sexualdifferenz von der Körpergroße. Bonner Zoologische Beiträge. 1950, 1: 58-69.
Andersson M, Iwasa Y: Sexual selection. Trends Ecol Evol. 1996, 11: 53-58.
Ralls K: Mammals in which females are larger than males. Q Rev Biol. 1976, 51: 245-276.
Ng MY, Wilbur HM: The cost of brooding in Plethodon cinereus. Herpetologica. 1995, 51: 1-8.
Duellman WE, Trueb L: Biology of Amphibians. 1994, The Johns Hopkins University Press, Baltimore
Kupfer A: Sexual size dimorphism in amphibians: an overview. Sex, size and gender roles: evolutionary studies of sexual size dimorphism. Edited by: Fairbairn DJ, Blanckenhorn WU, Szekely T. 2007, Oxford University Press, New York, 50-59.
Wells KD: The Ecology and Behaviour of Amphibians. 2007, The University of Chicago Press, Chicago
Shine R: Sexual selection and sexual dimorphism in the Amphibia. Copeia. 1979, 297-306.
Katsikaros K, Shine R: Sexual dimorphism in the tusked frog, Adelotus brevis (Anura: Myobatrachidae): the roles of natural and sexual selection. Biol J Linn Soc. 1997, 60: 39-51.
Tsuji H, Matsui M: Male-male combat and head morphology in a fanged frog (Rana kuhlii) from Taiwan. J Herpetol. 2002, 36: 520-526.
Zheng Y, Li S, Fu J: A phylogenetic analysis of the frog genera Vibrissaphora and Leptobrachium and the correlated evolution of nuptial spine and reversed sexual size dimorphism. Mol Phylogenet Evol. 2008, 46: 695-707.
Hudson C, He X, Fu J: Keratinized nuptial spines are used for male combat in the Emei moustache toad (Leptobrachium boringii). Asian Herpetological Res. 2011, 2: 142-148.
Zheng Y, Deng D, Li S, Fu J: Aspects of the breeding biology of the Omei mustache toad (Leptobrachium boringii): polygamy and paternal care. Amphibia-Reptilia. 2010, 31: 183-194.
Woolbright LL: Sexual selection and size dimorphism in anuran Amphibia. Am Nat. 1983, 121: 110-119.
Roff DA: Life History Evolution. 2002, Sinauer Associates Incorporation, Sunderland
Monnet J-M, Cherry MI: Sexual size dimorphism in anurans. Proc Royal Soc London Series B. 2002, 269: 2301-2307.
Felsenstein J: Phylogenies and the comparative method. Am Nat. 1985, 125: 1-15.
Harvey PH, Pagel MD: The Comparative Method in Evolutionary Biology. 1991, Oxford University Press, Oxford
Frost DR, Grant T, Faivovich J, Bain RH, Haas A, Haddad CFB, De Sa RO, Channing A, Wilkinson M, Donnellan SC, Raxworthy CJ, Campbell JA, Blotto BL, Moler P, Drewes RC, Nussbaum RA, Lynch JD, Green DM, Wheeler WC: The amphibian tree of life. Bull Am Mus Nat Hist. 2006, 297: 1-370.
Pyron RA, Wiens JJ: A large-scale phylogeny of Amphibia including over 2800 species, and a revised classification of extant frogs, salamanders, and caecilians. Mol Phylogenet Evol. 2011, 61: 543-583.
Lovich JE, Gibbons JW: A review of techniques for quantifying sexual size dimorphism. Growth Dev Aging. 1992, 56: 269-281.
Frost DR: Amphibian Species of the World: An Online Reference. Version 5.5. 2011, American Museum of Natural History, New York, USA, http://research.amnh.org/herpetology/amphibia/ (Accessed in March 2012)
Futuyma DJ: Evolution. 1993, Sinauer Associates Incorporation, Sunderland
Wells KD: The social behaviour of anuran amphibians. Anim Behav. 1977, 25: 666-693.
Freckleton RP, Harvey PH, Pagel M: Phylogenetic analysis and comparative data: a test and review of evidence. Am Nat. 2002, 160: 712-726.
Pagel M: Inferring the historical patterns of biological evolution. Nature. 1999, 401: 877-884.
Garland TJ, Harvey PH, Ives AR: Procedures for the analysis of comparative data using phylogenetically independent contrasts. Syst Biol. 1992, 41: 18-32.
Orme D, Freckleton R, Thomas G, Petzoldt T, Fritz S, Isaac N: Comparative Analyses of Phylogenetics and Evolution in R, version 0.5. 2012, http://caper.r-forge.r-project.org/,
Grant T, Frost DR, Caldwell JP, Gagliardo R, Haddad CFB, Kok PJR, Means DB, Noonan BP, Schargel WE, Wheeler WC: Phylogenetic systematics of dart-poison frogs and their relatives (Amphibia: Athesphatanura: Dendrobatidae). Bull Am Mus Nat Hist. 2006, 299: 1-262.
Zar JH: Biostatistical Analysis. 1996, Prentice Hall, Upper Saddle River, 3
Fritz SA, Purvis A: Selectivity in mammalian extinction risk and threat types: a new measure of phylogenetic signal strength in binary traits. Conserv Biol. 2010, 24: 1042-1051.
Pagel M: A method for the analysis of comparative data. J Theor Biol. 1992, 156: 431-442.
Burt A: Comparative methods using phylogenetically independent contrasts. Oxf Surv Evol Biol. 1989, 6: 33-53.
Abouheif E, Fairbairn DJ: A comparative analysis of allometry for sexual size dimorphism: assessing Rensch’s rule. Am Nat. 1997, 149: 540-562.
R Development Core Team: R: A language and environment for statistical computing. 2011, R Foundation for Statistical Computing, Vienna, Austria, http://www.R-project.org/,
Ye CY, Fei L, Hu SQ: Rare and Economic Amphibians of China. 1993, Sichuan Science and Technology publishing House, Chengdu
Cook CL, Ferguson JWH, Telford SR: Adaptive male parental care in the giant bullfrog. Pyxicephalus adspersus. J Herpetol. 2001, 35: 310-315.
Schütz D, Taborsky M: Giant males or dwarf females: what determines the extreme sexual size dimorphism in Lamprologus callipterus?. J Fish Biol. 2000, 57: 1254-1265.
Ho TC, Lathrop A, Murphy RW, Orlov NL: A redescription of Vibrissaphora ailaonica with a new record in Vietnam. Russ J Herpetol. 1999, 6: 48-54.
Breder CM, Rosen DE: Modes of reproduction in fishes. 1966, Natural History Press, New York
Summers K, McKeon CS, Heying H: The evolution of parental care and egg size: a comparative analysis in frogs. Proc Royal Soc London Series B. 2006, 273: 687-692.
Kolm N, Ahnesjö I: Do egg size and parental care coevolve in fishes?. J Fish Biol. 2005, 66: 1499-1515.
Balinsky BI, Balinsky JB: On the breeding habitat of the South African bullfrog, Pyxicephalus adspersus. South Afr J Anim Sci. 1954, 51: 55-58.
Shine R: The evolution of large body size in females: a critique of Darwin’s “fecundity advantage” model. Am Nat. 1988, 131: 124-131.
Gerhardt HC, Daniel RE, Perrill SA, Schramm S: Mating behaviour and male mating success in the green treefrog. Anim Behav. 1987, 35: 1490-1503.
Howard RD: The evolution of mating strategies in bullfrog, Rana catesbeiana. Evolution. 1978, 32: 850-871.
Fu J, Weadick CJ, Bi K: A phylogeny of the high-elevation Tibetan megophryid frogs and evidence for the multiple origins of reversed sexual size dimorphism. J Zool. 2007, 273: 315-325.
The authors would like to thank J Bogart, T Nudds, B Robinson and C Hudson for insightful comments and suggestions on earlier versions of this manuscript. This project is supported by a NSERC (Canada) discovery grant to JF.
The authors declare that they have no competing interests.
XH collected most of the data, conducted most analyses and drafted the manuscript. JF conceived the project, helped with the data collection and analysis, and finalized the manuscript. Both authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Mean body size in each sex, mean egg size, mean clutch size, mating combat, and parental care behavior in 688 anuran species. SDI = sexual dimorphism index [(female body size/male body size) – 1]. * R = range body size data, M = mean body size data. $ The presence of a trait is coded as “1”, and the absence of a trait is coded as “0”; missing data are represented by “-”. (DOCX 103 KB)
Additional file 2: References for body size, egg size, clutch size, mating combat, and parental care in 688 anuran species. The general mating or breeding descriptions are used to define the absence of mating combat or parental care behaviour, respectively. (DOCX 106 KB)
Additional file 3: The D statistic for all binary traits. The non-significant p values are in bold, which means the traits are under Brownian evolution. (DOCX 52 KB)
Additional file 4: Results from simple linear regression analyses on the ‘extended dataset’ using phylogenetic generalized least squares (PGLS) model and phylogenetic independent contrasts (PIC). PIC analyses are conducted at both all-anuran and family levels; outliers, which deviated from the majority of the data points by more than three interquartile ranges from the quartiles, were excluded. All species with relevant data are included in each pairwise correlation analysis. Only correlations with degrees of freedom ≥ 3 are presented. P < 0.10 are in bold. * Dicro = Dicroglossidae; the family has ‘mean’ data only. (DOCX 136 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Han, X., Fu, J. Does life history shape sexual size dimorphism in anurans? A comparative analysis. BMC Evol Biol 13, 27 (2013). https://doi.org/10.1186/1471-2148-13-27
- Sexual size dimorphism
- Phylogenetic comparative analysis
- Life history
- Male combat
- Parental care
- Fecundity advantage
- Rensch’s rule