Evolution of sociality by natural selection on variances in reproductive fitness: evidence from a social bee
© Stevens et al; licensee BioMed Central Ltd. 2007
Received: 29 November 2006
Accepted: 29 August 2007
Published: 29 August 2007
The Central Limit Theorem (CLT) is a statistical principle that states that as the number of repeated samples from any population increase, the variance among sample means will decrease and means will become more normally distributed. It has been conjectured that the CLT has the potential to provide benefits for group living in some animals via greater predictability in food acquisition, if the number of foraging bouts increases with group size. The potential existence of benefits for group living derived from a purely statistical principle is highly intriguing and it has implications for the origins of sociality.
Here we show that in a social allodapine bee the relationship between cumulative food acquisition (measured as total brood weight) and colony size accords with the CLT. We show that deviations from expected food income decrease with group size, and that brood weights become more normally distributed both over time and with increasing colony size, as predicted by the CLT. Larger colonies are better able to match egg production to expected food intake, and better able to avoid costs associated with producing more brood than can be reared while reducing the risk of under-exploiting the food resources that may be available.
These benefits to group living derive from a purely statistical principle, rather than from ecological, ergonomic or genetic factors, and could apply to a wide variety of species. This in turn suggests that the CLT may provide benefits at the early evolutionary stages of sociality and that evolution of group size could result from selection on variances in reproductive fitness. In addition, they may help explain why sociality has evolved in some groups and not others.
The Central Limit Theorem (CLT) states that variance among sample means is inversely related to sample size. Furthermore, as sample sizes increase, sample means become more normally distributed, even if the underlying probability density function is asymmetrical or deviates from normality in other ways. Wenzel and Pickering  suggested that the CLT entails benefits for sociality in that larger colonies are able to better predict future food acquisition. For example, consider a species where food brought back to the group by foragers is shared among immature offspring. Each foraging trip can be regarded as a stochastic sample of the environment, with the number of foraging trips per unit time corresponding to sample size . Larger group sizes could enable the number of foraging trips per unit time to increase, either because the number of foragers increases, or because group living allows dedicated foragers that are not required to spend time and effort on other nest related tasks (such as cleaning, brood tending, and nest defence) to make more trips. In such cases, a greater number of foraging trips should lead to lower variation among colonies in foraging yield . Consequently, larger groups should produce clutch sizes that match future food acquisition and therefore run a lower risk of either under- or over-producing eggs relative to future resources. If future food acquisition is unpredictable, colonies may either produce more eggs than can be fed and then having to abort some (a 'no wasted food' strategy), or they may under-produce eggs and be unable to fully exploit all the resources that may become available (a 'no wasted brood' strategy) .
In addition to the arguments put forward by Wenzel and Pickering , Gillespie  argued that variance in offspring number can be as important as the mean number of offspring produced when calculating the fitness of alternative genotypes and this holds when variance is due to both environmental and developmental sources of stochasticity. Gillespie showed that selection will favour lower variance in offspring number and the strength of selection will be equal to that on mean numbers in the case of environmental fluctuations. This suggests a very large effect of group size if increasing the group size lowers the variation in offspring production.
The idea that group living may entail benefits because of a purely statistical principle is intriguing as it adds a completely new light to the evolution of sociality. Traditionally, studies looking at the advantages of sociality have focused on benefits deriving from ecological, ergonomic or genetic factors [3–9]. Wenzel and Pickering  found a lower percentage of brood abortion events in larger colonies, but did not directly address the variance in brood numbers in relation to the number of females present. Since then there have only been two studies based on a social wasp (Ropalidia marginata) that have explicitly examined the CLT in relationship to group living [10, 11]. Although Shakarad and Gadagkar  established a decreasing variance in brood production with increasing colony size, they were not able to distinguish between benefits derived as a result of the CLT or from a better prevention of brood loss through protection from predation. Naug and Wenzel  used behavioural observations from five colonies of the same species and showed that an increase in predictability of foraging success (food supply to the nest) was consistent with the CLT predictions. It appears from these studies that, at least for R. marginata, assessing predictability in reproductive success (directly or indirectly) has been problematic.
Our approach is two-pronged. Firstly, we investigate the CLT-derived predictions that total brood weight will show less variation and become more normally distributed as colony size increases. These patterns are not predicted by ecological or ergonomic models, which instead predict that brood rearing becomes more efficient or that brood predation is decreased as colonies become larger. Secondly, we investigate whether initial clutch sizes correspond better with final brood sizes for larger colonies than for smaller colonies.
Shapiro-Wilk tests for deviation from normality in the distribution of total brood weights for different colony sizes in the late spring (6 November) and early summer (4 December) samples
Colony size (females/nest)
Sample size (number of colonies)
Re-sample comparison (Sig.)
Given the widely varying sample sizes for different-sized colonies, we used a re-sampling procedure to ask whether colony-size variation in Shapiro-Wilk values (Table 1) may have been attributable to differences in sample sizes, rather than any underlying differences in normality. The proportion of re-sampled 1-female colony brood weight data with a higher Shapiro-Wilk value than the corresponding multi-female sample approximates the likelihood that the multi-female data is more normally distributed than the 1-female data, and is given in Table 1 for each colony size and sample date. The only cases where these proportions were not less than 5% (i.e. P < 0.05) involved 4-female colonies. However, these involved samples sizes of only 10 and 5, for late spring and early summer respectively, where the power of our re-sampling technique is likely to be low. Consequently, our analyses indicate that variation in measures of normality in brood weights is not simply due to variation in sample sizes.
In our second approach, we used ANCOVA with standard deviation as the dependent variable, sample date as the treatment, and mean brood weight and colony size as covariates. Variances were not heteroscedastic (Levene's test, P = 0.388), and there was a significant effect of colony size (F1,8 = 6.891, P = 0.03) and mean brood weight (F1,8 = 34.586, P < 0.001), but no effect of sample date (F1,8 = 0.113, P = 0.746). We then pooled the two samples and regressed the standard deviation of brood weight for each colony size-sample date combination onto mean brood weight (F1,10 = 25.387, P = 0.001, β = 0.847) and saved the residuals. These residuals were then regressed onto colony size with results indicating a significant negative relationship (F1,10 = 5.575, P = 0.040, β = -0.598). The magnitude of β here suggests a substantial negative effect of colony size on the standard deviation once the effects of mean brood weight are taken into account, while the positive β value for the initial regression indicates that the standard deviation increases with mean brood weight, as expected.
Consequently, we further assessed this statistically by comparing per capita brood numbers for each colony size category between spring and summer samples using Mann-Whitney U-tests, restricting analyses to colonies with four or fewer females because of small sample sizes for larger nests. Because we used a non-parametric test for comparisons between spring and summer colonies with the same number of adults, these tests also correspond to comparing total brood size. Brood numbers significantly declined between spring and summer for 1-female colonies (Mann-Whitney U test Z = -2.338, P = 0.019, asymptotic 2-tailed significance) and 2-female colonies (Z = -2.561, P = 0.010), but not for 3-female (Z = -0.223, P = 0.823) or 4-female colonies (Z = -1.467, P = 0.149). Thus, in small colonies the number of offspring being reared by summer was significantly smaller than the initial clutch size, indicating that brood had been reduced, while for large colonies final offspring numbers corresponded with initial clutch size.
Our study provides several lines of evidence for the existence of benefits that derive from the CLT, rather than from ecological, ergonomic, or genetic factors. Firstly, we have shown that colonies with an increasing number of females showed an increasingly normal distribution of brood weight. For colonies with the same number of females, the brood weight became more normally distributed as the season progressed (Figure 2; Table 1). Secondly, we have shown a clear decrease in variation in total brood weight with increasing colony sizes (Figure 3). Interestingly, the decrease in CV with increasing colony size appears to be lower for the early summer sample (Figure 4 and Fligner-Killeen tests). This is not surprising, since by early summer all colony size categories would have experienced a greater number of foraging trips than in late spring, and under the CLT this should lead to greater normality in the distribution of total food income.
The mean per capita number of brood differed with colony size early in the season, and was higher for small colonies than for larger colonies. By early summer, when brood was more mature, this was no longer the case (Figure 5). This suggests that some small colonies may be initially over-producing brood, a 'no wasted food' strategy. In contrast, larger colonies appear better able to predict future food income and final clutch sizes and therefore better able to avoid/lessen the costs of aborting excess brood in the event of low food income, a 'no wasted brood' strategy.
The benefits that the CLT might provide to larger colonies of E. nigrescens cannot be quantified from our data because we do not know the costs involved in abortion of brood, but these are likely to be substantial. Brood are reared progressively, so that resources are sequestered in the developing brood over time, and if not all brood can be reared then already-invested resources will be lost. Oophagy and larval cannibalism in E. nigrescens are extremely rare  so it is unlikely that investment in aborted eggs and larvae is recovered to any large extent. Importantly, egg size in allodapines is very large (egg length is approximately 30% of total body length, see Figure 1C) compared to most other social insects , so that over-production of eggs is likely to be costly, and abortion of larvae, where even further resources have been invested, will be more so. Therefore, colonies that can better match brood numbers to food acquisition should have a selective advantage. Such benefits derived from the CLT clearly refute Alexander's  assertion that there are three automatic and universal costs to group living but no universal benefit.
Even if colonies are not able to benefit from better matching egg production to expected food income, the lower variation in mean brood weight should provide benefits via the arguments suggested by Gillespie  for selection on lower variance in offspring numbers. The arguments put forward by Wenzel and Pickering , and supported by our results, suggest that lower variation in brood weight is due to the effect of group size on variation in food acquisition, ultimately deriving from stochasticity in the environment. Gillespie's  models suggest that the selective effects of this may be equivalent to changes in mean brood numbers.
Group living, especially where it involves behavioural specialisation, is likely to provide benefits from many sources, such as enhanced task efficiency [3, 4, 9, 18–20] and protection against predators [21–23]. However, the CLT has the potential to enable benefits for group living without involving any behavioural specialisation and requires only that larger groups have higher numbers of foraging trips and that the incoming food is shared by the group [1, 11]. Indeed, Naug and Wenzel  demonstrated that an increase in predictability of foraging success (food supply to the nest) in the social wasp R. marginata was consistent with the CLT predictions. Such basic requirements suggest that the CLT may provide some benefits at the early evolutionary stages of sociality, and that selection can act on variances in reproductive fitness [see ].
Given the potential for the CLT to provide benefits to group living, it is surprising that more studies have not assessed it as an explicit factor shaping group performance as a function of group size. However, increased reliability and predictability in offspring production has recently been suggested as a factor influencing colonial nesting in birds , and may be important in shaping the group size of communally breeding birds. Reed and Walters  investigated the effect of the number of helpers on variance in reproductive success in Red-cockaded Woodpeckers (Picoides borealis) and, using the literature available, in other communally breeding birds. They found that the presence of helpers was often associated with higher rather than lower variances, but pointed out that this result could have been confounded by effects of habitat quality, a factor that was ruled out in our experimental set up . In contrast, data presented by Woxvold  suggest a lower variance with increasing group size for apostle birds (Struthidea cinerea), emphasizing that breeding systems, among other factors, need to be considered.
There are likely to be a large number of factors that influence when and how group living evolves [e.g. [5, 20, 27]], but the non-random distribution of social origins with respect to life-history traits  suggest that some factors are more important than others. The ability to pool group resources when provisioning offspring means that foraging bouts, which represent samples from the environment, increase with group size. The importance of large sample sizes are well known to biologists – they increase our ability to more accurately estimate population parameters, and ultimately this ability derives from the CLT. It would be surprising if this statistical property has been unexploited during the tinkerings of natural selection. The importance of CLT needs to be further explored in other social groups to determine its wider potential facilitating role in social evolution.
Location and data collection
The study was carried out in Cobboboonee State Forest, Victoria, Australia (38°77'S, 141°35'E). Occupied nests (Xanthorrhea minor flower scapes, Figure 1B) were collected from an area approximately 400 m by 2000 m on 31 July 1996 , when colonies were overwintering and before egg production had begun. We placed equal numbers of nests (randomly allocated) in four nearby sites . Over three collecting trips (21 September (early spring), 6 November (late spring), 4 December (early summer) 1996) an equal number of E. nigrescens nests were collected from the four sites . Nests were stored at 10°C and each collection processed within two days. For each colony we obtained a weight (mg) for the total number of eggs and for the total number of larvae. Because the number of adult females also relates to the possible information gathered about future food acquisition, irrespective of who is reproductive at any one time , we define brood per capita as 'the number of brood relative to the number of adult females', which was calculated by dividing the total brood (eggs, small larvae, medium larvae, large larvae, pre-pupae, and pupae) in the nest by the number of adult females in the nest. For each site and field collection, the mean and standard deviations in the number of brood per capita and brood weight were calculated for 1-female, 2-female, 3-female, and 4 (or more)-female colonies.
Data were analysed with SPSS version 11 for Macintosh. Where assumptions of homoscedasticity could not be met and where any variable was ordinal or categorical, non-parametric tests were used. Previous studies of Exoneura show that per capita brood number increases with colony size [22, 29, 30] and this is largely due to avoidance of total brood loss to enemies-at-the-nest  – allodapine brood are not enclosed within cells (see Figure 1C–E) and are therefore highly vulnerable to predation . To avoid this confounding ecological benefit of group living, we restricted our analyses to nests with brood. The nest contents did not differ between sites for any of the sampling dates (Kruskal-Wallis 1-way ANOVA, all P > 0.05)  and were thus pooled into one collection per sampling date.
We assessed normality of total brood weights for 1, 2, 3 and 4-female colonies using quantile-quantile (Q-Q) plots and Shapiro-Wilk tests. However, the power of such tests may vary with sample sizes, and in our case these varied widely with colony size. In order to examine whether sample size was sufficient to explain the outcome of the Shapiro-Wilk tests, we used a re-sampling procedure, which compared measures of normality for total brood weight in 1-female colonies with 2, 3 and 4-female colonies, larger colonies were not used because they were too few for meaningful comparison. We re-sampled, with replacement, total brood weights from 1-female colonies using the same sample sizes that we had for colonies with 2, 3 or 4 females (Table 1). Colonies with 2 females were more common than 1-female colonies (Table 1), so for these we randomly selected a subset of the 2-female colonies to produce the same sample size as for 1-female colonies and calculated the Shapiro-Wilk values for these subsets. Re-sampling with replacement asks whether an underlying distribution of brood weights for 1-female colonies are likely to yield more-normal distributions for 2, 3 or 4-female colonies for each of the respective sample sizes. For each protocol, we used 100 random re-samples and determined the proportion of cases that had higher values of the Shapiro-Wilk statistic than was calculated for the multi-female colonies. This proportion represents a one-tailed significance test because in each case we explicitly asked whether multi-female colonies had more normally distributed brood weights than 1-female colonies.
Variation around sample means can be estimated using the standard deviation but interpretation of the SD can be problematic if variation increases with the mean, a common phenomenon . In such cases the coefficient of variation (CV) is a more useful tool [10, 11, 32]. The CV is independent of the unit of measure, and is the standard deviations (SD) expressed as a proportion of the means (y), so CV = SD/y, often converted to be a percentage of the mean; but here we use it as a proportion. We corrected for small sample sizes by multiplying the CV by (1 + 1/4n) . However, calculating the CV without assessing if the standard deviations are affected by the mean may result in an over-compensation and produce false trends . Therefore, the relationship between the mean and SD of brood weight was analysed using an ANCOVA.
Comparing CV across different groups is not straightforward, and most existing methods have high Type I error rates and are not robust to differences in underlying distributions, such as differences in kurtosis or skewness . The modified Fligner-Killeen test  suffers least from these drawbacks but is designed for only pairwise comparisons and requires that the smallest sample size for any comparison be ≥7 . We used the Fligner-Killeen test to compare CV for 1-female colonies with colonies of 2 and 3-females. Small sample sizes precluded comparisons between colonies in the other size categories.
As a further check on results from analysis of the CV, we examined the SD itself. For this procedure we calculated the SD of total brood weight for each colony size for each sample date, and then regressed the SD onto mean brood weight, taking sample dates into account. Residuals were saved and then regressed onto colony size to determine if colony size contributed to the SD independently of mean brood weight.
We apply several biologically realistic assumptions here, namely that (1) brood weight is a proxy for foraging effort; (2) the existence of differing male and female metabolic rates do not affect our data interpretation; and (3) different metabolic rates do not influence normality of brood weights (increasing normality in larger colonies and over time). These assumptions hold for several reasons [see [33–37]]. Most importantly, if sex ratios are constant for different sized colonies, this would be equivalent to changing the total weight of each colony by a set amount and this would not affect our analyses since we used the CV (see above), which removes the effect of mean colony brood weights. Moreover, because sex ratios vary strongly with colony size in E. nigrescens [15, 35, 36] gender differences in metabolic rates will result in a greater variance in brood weights for larger colonies where there is the possibility for greater variation in sex ratios. Because larger colonies have a much greater incidence of male immatures, then there would be a greater effect of differential metabolic rates, leading to a greater range in individual larval weights and hence a more platykurtic distribution of colony brood weights. However, because 1-female colonies have in most cases only females , the distribution of brood weights should tend towards normality.
We are grateful to Nathalie Loussert, Nick Bull, Remko Leys, Adam Cronin, and Trevor Lehmeyer for help with field and/or lab work, and to David Penny, Michael Hendy, Philipp Langer, Mark Blows, Klaus Schliep, and three anonymous reviewers for helpful discussions and/or comments on previous drafts. This research was supported by Flinders University and Australian Research Council grants to MPS.
- Wenzel JW, Pickering J: Cooperative foraging, productivity, and the central limit theorem. Proc Natl Acad Sci USA. 1991, 88: 37-38. 10.1073/pnas.88.1.36.View ArticleGoogle Scholar
- Gillespie JH: Natural selection for variances in offspring numbers: a new evolutionary principle. Am Nat. 1977, 111: 1010-1014. 10.1086/283230.View ArticleGoogle Scholar
- Hamilton WD: The genetic evolution of social behaviour, I and II. J Theoret Biol. 1964, 7: 1-52. 10.1016/0022-5193(64)90038-4.View ArticleGoogle Scholar
- Michener CD: Reproductive efficiency in relation to colony size in Hymenopterous societies. Insect Soc. 1964, 11: 317-342. 10.1007/BF02227433.View ArticleGoogle Scholar
- Wilson EO: The insect societies. 1971, The Belknap Press of Harvard University Press, Cambridge, MassachusettsGoogle Scholar
- Anderson C, Ratnieks FLW: Task partitioning in insect societies. I. Effect of colony size on queueing delay and colony ergonomic efficiency. Am Nat. 1999, 154: 521-535. 10.1086/303255.View ArticlePubMedGoogle Scholar
- Langer P, Hogendoorn K, Keller L: Tug-of-war over reproduction in a social bee. Nature. 2004, 428: 844-847. 10.1038/nature02431.View ArticlePubMedGoogle Scholar
- Wilson EO, Hölldobler B: Eusociality: origin and consequences. Proc Natl Acad Sci USA. 2005, 102: 13367-13371. 10.1073/pnas.0505858102.PubMed CentralView ArticlePubMedGoogle Scholar
- Fjerdingstad EJ, Crozier RH: Worker caste diversity in social insects. Am Nat. 2006, 167: W390-400. 10.1086/499545.View ArticleGoogle Scholar
- Shakarad M, Gadagkar R: Colony founding in the primitively eusocial wasp, Ropalidia marginata (Hymenoptera: Vespidae). Ecol Entomol. 1995, 20: 273-282.View ArticleGoogle Scholar
- Naug D, Wenzel J: Constraints on foraging success due to resource ecology limit colony productivity in social insects. Behav Ecol Sociobiol. 2006, 60: 62-68. 10.1007/s00265-005-0141-5.View ArticleGoogle Scholar
- Schwarz MP: Local resource enhancement and sex ratios in a primitively social bee. Nature. 1988, 331: 346-347. 10.1038/331346a0.View ArticleGoogle Scholar
- Schwarz MP, Silberbauer LX, Hurst PS: Intrinsic and extrinsic factors associated with social evolution in allodapine bees. The evolution of social behavior in insects and arachnids. Edited by: Choe JC, Crespi BJ. 1997, Cambridge: Cambridge University Press, 333-346.View ArticleGoogle Scholar
- Donnelly SM, Kramer A: Testing for multiple species in fossil samples: an evaluation and comparison of tests for equal relative variation. Am J Phys Anthropol. 1999, 108: 507-529. 10.1002/(SICI)1096-8644(199904)108:4<507::AID-AJPA8>3.0.CO;2-0.View ArticlePubMedGoogle Scholar
- Stevens MI: Egg stockpiling strategies and benefits to group living in Exoneura bicolor Smith (Apidae, Xylocopinae). Thesis. 1997, Flinders University, Adelaide, AustraliaGoogle Scholar
- Michener CD: Size and form of eggs of allodapine bees. J Entomol Soc South Africa. 1973, 36: 281-285.Google Scholar
- Alexander RD: The evolution of social behavior. Annu Rev Ecol Syst. 1974, 5: 325-383. 10.1146/annurev.es.05.110174.001545.View ArticleGoogle Scholar
- Frank SA: Policing and group cohesion when resources vary. Anim Behav. 1996, 52: 1163-1169. 10.1006/anbe.1996.0263.View ArticleGoogle Scholar
- Karsai I, Wenzel JW: Productivity, individual-level and colony-level flexibility, and organization of work as consequences of colony size. Proc Natl Acad Sci USA. 1998, 95: 8665-8669. 10.1073/pnas.95.15.8665.PubMed CentralView ArticlePubMedGoogle Scholar
- Amdam GV, Csondes A, Fondrk MK, Page RE: Complex social behaviour derived from maternal reproductive traits. Nature. 2006, 439: 76-78. 10.1038/nature04340.PubMed CentralView ArticlePubMedGoogle Scholar
- Hogendoorn K: Socio-economics of brood destruction during supersedure in the carpenter bee Xylocopa pubescens. J Evol Biol. 1996, 9: 931-952. 10.1046/j.1420-9101.1996.9060931.x.View ArticleGoogle Scholar
- Schwarz MP, Bull NJ, Hogendoorn K: Evolution of sociality in the allodapine bees: a review of sex allocation, ecology and evolution. Insect Soc. 1998, 45: 349-368. 10.1007/s000400050095.View ArticleGoogle Scholar
- Soucy SL, Giray T, Roubik DW: Solitary and group nesting in the orchid bee Euglossa hyacinthina (Hymenoptera, Apidae). Insect Soc. 2003, 50: 248-255. 10.1007/s00040-003-0670-8.View ArticleGoogle Scholar
- Kosciuch KL, Langerhans RB: Evolution of coloniality via commodity selection: What about variance?. The Auk. 2004, 121: 257-262. 10.1642/0004-8038(2004)121[0257:EOCVCS]2.0.CO;2.View ArticleGoogle Scholar
- Reed MJ, Walters JR: Helper effects on variance components of fitness in the cooperatively breeding red-cockaded woodpecker. The Auk. 1996, 113: 608-616.View ArticleGoogle Scholar
- Woxvold IA: Breeding ecology and group dynamics of the apostlebird. Aust J Zool. 2004, 52: 561-581. 10.1071/ZO04031.View ArticleGoogle Scholar
- Queller DC: Extended parental care and the origin of eusociality. Proc Roy Soc Lond, B. 1994, 256: 105-111. 10.1098/rspb.1994.0056.View ArticleGoogle Scholar
- Choe J, Crespi B: The Evolution of Social Behaviour in Insects and Arachnids. 1997, Cambridge: Cambridge University PressView ArticleGoogle Scholar
- Tierney SM, Cronin AL, Loussert N, Schwarz MP: The biology of Brevineura froggatti and phylogenetic conservatism in Australian allodapine bees (Apidae, Allodapini). Insect Soc. 2000, 47: 96-97. 10.1007/s000400050016.View ArticleGoogle Scholar
- Tierney SM, Schwarz MP, Neville T, Schwarz PM: Social behaviour in the African bee genus Macrogalea and implications for the origin of sociality in the Allodapini. Biol J Linn Soc. 2002, 76: 211-224. 10.1046/j.1095-8312.2002.00059.x.View ArticleGoogle Scholar
- Hogendoorn K, Zammit J: Benefits of cooperative breeding through increased colony survival in an allodapine bee. Insect Soc. 2001, 48: 392-397. 10.1007/PL00001796.View ArticleGoogle Scholar
- Sokal RR, Rohlf FJ: Biometry. 1995, San Francisco: FreemanGoogle Scholar
- Klostermeyer EC, Mech SJ, Rasmussen WB: Sex and weight of Megachile rotundata (Hymenoptera: Megachilidae) progeny associated with provision weights. J Kans Entomol. 1973, 46: 536-548.Google Scholar
- Danforth BN: Provisioning behavior and the estimation of investment ratios in a solitary bee, Calliopsis (Hypomacrotera) persimilis (Cockerell) (Hymenoptera: Andrenidae). Behav Ecol Sociobiol. 1990, 27: 159-168. 10.1007/BF00180299.View ArticleGoogle Scholar
- Schwarz MP: Female-biased sex ratios in a facultatively social bee and their implications for social evolution. Evolution. 1994, 48: 1684-1697. 10.2307/2410257.View ArticleGoogle Scholar
- Bull NJ, Schwarz MP: Brood insurance via protogyny: a source of female biased sex allocation. Proc Roy Soc Lond, B. 2001, 268: 1869-1874. 10.1098/rspb.2001.1687.View ArticleGoogle Scholar
- Bosch J, Vicens N: Body size as an estimator of production costs in a solitary bee. Ecol Entomol. 2002, 27: 129-137. 10.1046/j.1365-2311.2002.00406.x.View ArticleGoogle Scholar
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.