- Research article
- Open Access
The evolution of activity breaks in the nest cycle of annual eusocial bees: a model of delayed exponential growth
© Mitesser et al; licensee BioMed Central Ltd. 2006
- Received: 25 December 2005
- Accepted: 02 June 2006
- Published: 02 June 2006
Social insects show considerable variability not only in social organisation but also in the temporal pattern of nest cycles. In annual eusocial sweat bees, nest cycles typically consist of a sequence of distinct phases of activity (queen or workers collect food, construct, and provision brood cells) and inactivity (nest is closed). Since the flight season is limited to the time of the year with sufficiently high temperatures and resource availability, every break reduces the potential for foraging and, thus, the productivity of a colony. This apparent waste of time has not gained much attention.
We present a model that explains the evolution of activity breaks by assuming differential mortality during active and inactive phases and a limited rate of development of larvae, both reasonable assumptions. The model predicts a systematic temporal structure of breaks at certain times in the season which increase the fitness of a colony. The predicted pattern of these breaks is in excellent accordance with field data on the nest cycle of the halictid Lasioglossum malachurum.
Activity breaks are a counter-intuitive outcome of varying mortality rates that maximise the reproductive output of primitively eusocial nests.
- Activity Break
- Colony Development
- Eusocial Insect
- Nest Cycle
- Colony Cycle
Some of the most important components of life history decisions refer to the optimal timing of accumulation of resources and their allocation to growth and reproduction. At any time during its life an organism has not only to spend resources on the conflicting requirements for maintenance, somatic growth, and reproduction, but has also to decide on how much and when resources like food or building material should be accumulated in order to maximise reproductive output. Up to now, theoretical studies on life history strategies of eusocial insects have mainly focused on the first aspect: optimal resource allocation [1–4]. However, the obvious and ample variability in seasonal activity patterns within and between species of eusocial insects requires investigating the optimal timing of resource accumulation too.
Seasonal activity patterns vary widely among the species of bees and wasps that have been studied as model organisms for the evolution of sociality in insects [5–7]. Many annual Polistes, Vespa, Xylocopa and Allodape species show continuous colony activity during the whole season [8–11]. This results in a more or less continuous production of offspring as is assumed in the classical model of colony development by Macevicz and Oster (1976). Model predictions have been tested and were met in field data from Polistes and Vespa species .
However, the nest cycle of most halictids (e.g. in the genera Lasioglossum and Halictus) is characterized by several discrete broods that are separated by distinct activity breaks [[6, 12, 13], but see  and ]. During the solitary founding phase, halictid queens construct nests and supply brood cells with pollen and nectar as provisions for their larvae. After a break of a few weeks, during which the nest is closed and no activity outside the nest can be observed, a first worker brood emerges and starts collecting pollen and nectar to provision the eggs that are usually laid by the queen. Subsequent broods are also separated by breaks during which nests are closed and no outside activity can be observed. Activity breaks can last up to three weeks [[12, 16], Weissel et al., submitted]. Usually sexuals emerge in the last brood only, while all other broods consist mainly of workers. There are also species with an intermediate position between continuous growth and discrete broods. In Bombus species, for example, the nest cycle is organized in more or less discrete broods but without activity breaks and nest closure .
Due to temperature-dependence of their activity and resource availability, ectothermic organisms, like insects, have to adjust their life history to the seasonal conditions in temperate latitudes. Reproduction and growth must be completed within a limited time span and the unfavourable period has to be bridged by diapause. Variability in biotic and abiotic conditions during the reproductive period has been assumed to cause changes in activity patterns on a smaller time scale . Usually bees will forage during the day when visibility is good and temperature is high enough for flying and foraging [19–23]. However, the evolutionary transitions to dim-light foraging have occurred repeatedly in bees, and may be associated with the escape from enemies or competitors [24, 25]. The daily activity patterns of the solitary bee Anthophora plumipes has been attributed to variation in the thermal environment as well as quality and quantity of floral resources . The pattern of activity breaks in halictids has accordingly been related to patterns of resource availability and photoperiod [14, 16, 27]. By contrast, Kaitala et al.'s (1990) modelling approach for L. malachurum assumes synchronized nest closure in halictid nest aggregation to be due to an increasing threat of nest usurpation by intraspecific parasites, so called floater queens.
Furthermore one could suppose that activity breaks after the emergence of the first workers just appear when worker mortality is rather high and all workers of a brood have died before emergence of the individuals of a second brood. However, it is clear from field observations that the breaks do not occur simply because all workers of a brood have died. Some workers even survive a complete activity break and continue foraging when the nest is reopened [, Weissel and Strohm unpublished]. On the contrary, breaks occur even though there are still some workers alive in a nest, showing that there has to be some advantage of interrupting foraging activity.
The well-known colony growth model of Macevicz & Oster (1976) for insect colonies identifies the sequential production of workers first and sexuals just before the ending of the flight season (so called bang-bang reproduction) as the optimal investment strategy to maximize colony fitness. Whereas this model assumes instantaneous occurrence of adult progeny the model that we present accounts for a certain development time of the larvae. The results of our model challenge the assumption that only variation in environmental factors governs the emergence of activity breaks. The model explains the evolution of the observed activity patterns rather by an asymmetric interaction between endogenous and exogenous factors of colony development.
A simple model of delayed exponential growth with activity breaks
We use a simple difference equation model to analyse colony development during a season of length L. Two main dependent variables describe the state of a colony: the number of workers (W i ) at time step i and the number of sexuals (S i ) at that time. For simplicity we do not distinguish between male and female sexuals . The colony cycle typically starts in spring with nest founding by inseminated hibernated queens. During the founding phase the queen works alone and performs all the foraging tasks that will be taken over by workers after their emergence later in the season . Thus we start with initial condition W1 = 1 assuming that the founding queen acts like a single worker until the first eggs have developed to adults . The dynamics in the number of nestmates is governed by two mechanisms: mortality and reproduction. Each individual survives from time step i to i + 1 with a probability q i (that might vary with time step i during the season). Resource allocation in each time step (i) is directly proportional to the current worker force (W i ). Each worker can provision c i (worker efficiency) eggs (= brood cells) per time step. We assume that the actual egg laying rate of the queen is only limited by the number of eggs that can be successfully provisioned by the workers of the colony . Adults emerge after a development time T. Halictid colonies suffer from numerous threats during activity periods (see discussion), so nest and especially brood mortality are rather high [Weissel et al., submitted]. As this parameter is not in the main focus of our analysis and field data are not readily available, we use the same survival probability (q i ) for eggs and for adult workers to keep the model simple (Additional numerical calculations have shown that our results differ only quantitatively if we uncouple worker and brood mortality). Additionally we assume that development time (T) does not correlate with either season length (L), caste or onset of development.
The portion u i of resources spent in time step i is allocated to new workers. Consequently, the portion (1 - u i ) is invested in sexuals S i . Thus, the number of workers (Wi+1) at time step i+1 can be calculated as
Wi+1= q i W i + u i ci-TWi-TΠ i i-Tq i (1)
In most halictid bees, life span of adult females is much longer than life span of workers . Female sexuals have to hibernate before nest founding in the following year, while workers live only for several weeks. Thus we neglect mortality of sexuals as has been done by Oster and Macevicz (1976) in most of their analyses and thus we get for the number of sexuals (Si+1) at time step i+1
Si+1= S i + (1 - u i )ci-TWi-TΠ i i-Tq i (2)
These two equations describe the delayed (by development time T) exponential growth of an annual, primitively eusocial bee colony. Fitness of colonies following such nest dynamics can be measured by the final number of sexuals S L . Oster and Wilson (1978) have studied such systems (in time continuous form and without delay) as optimal control problems with control variable u i (allocation in workers) . They found that the (time-dependent) optimal control solution that maximizes S L is switching in u i from 1 to 0 at an optimal point in time (dichotomous bang-bang strategy, SWT = switching time). So the optimal temporal pattern of reproduction consists of two distinct phases: exclusive worker production followed by exclusive sexual production. This result also holds for delay systems [2, 29]. In our simple model the optimal switching point can be found by a simple argument: Switching should take place when an egg just laid can not mature, eclose and contribute to rearing other individuals anymore. From time L - T to L no eggs should be produced at all, because they would not emerge before the season ends. The last contribution of a worker to sexual production can occur at time L - T - 1 and thus the last worker egg should be laid at L - 2T - 1. So we choose
for further analysis.
Analogously additional activity breaks can be inserted into the nest cycle. However, here we focus on a single break as our main results are not changed by the simultaneous consideration of several breaks. To answer the question of whether there are activity breaks that increase fitness when inserted into the nest cycle we analysed the complete B1-B2 parameter space by simple numerical calculations. We interpret B1 and B2 as life history parameters of the queen, who decides when to close and reopen her nest . During a single B1-B2 space simulation all other parameters were kept constant. Computer simulations are conducted with the programming language R .
Estimating model parameters
We calibrated the model with data from the halictid species L. malachurum. In this species a typical season in central Europe lasts for about 100 to 140 days. Since the absolute length of the season did not change our results within this range we choose L = 100. The mean life time of L. malachurum workers is about 24 days . This results in a survival probability of q = 0.95 per time step. Development from egg to adult typically lasts T = 20 time steps . There are no data available about the shelter effect of nest closure, so we studied the effect of within-break survival Q in the range from q (= 0.95) to 1. Worker efficiency c from 0.5 to 1.0 (per time step) results in an output of about 20 to 80 sexuals, similar to typical colony sizes in field observations .
However, increased protection during breaks (fig. 2b–2d) causes a second and a third area of beneficial times for breaks (fig. 2b and 2c) (type I and II). The position of these spots remains constant whereas their area increases with increasing break survival probability rate Q resulting in an extended area of beneficial break timings.
Next we reduced efficiency during breaks but kept survival rate (the gain part in fitness balance) constant. This has a more complex effect (fig. 4b). In the last period of the colony cycle (fig. 4b, S5) brood production has ceased and thus, fitness is not affected by late reduction in efficiency. During the prior period (starting at switching time SWT, fig. 4b, S4) only eggs for sexuals are laid which contribute to fitness in an additive way. Consequently fitness loss (caused by reduced worker efficiency) is directly proportional to number of workers and as worker number increases during this period, fitness decreases. During the solitary phase of colony development (fig. 4b, S1) sensitivity to productivity loss decreases with time (and relative fitness increases). This effect is not changed at all by different mortality rates (not shown in figure 4b). It is solely caused by the decreasing value of eggs developing into workers. The later a worker egg is produced in the season the less it contributes to overall fitness. At the beginning of the intermediate periods (S2 and S3) the first workers emerge and the oscillating pattern of the relative fitness function is governed by the interaction of both processes acting separately in periods S1 and S4.
In our modelling approach the emergence of activity breaks is caused by an increased protection of developing larvae and provisioned brood cells when colonies are closed during activity breaks. The main difference from Macevicz and Oster's (1976) model is the consideration of development time. The essential predictions of our model are (1) that there are activity breaks that increase colony fitness (and might consequently be favoured by selection), (2) a clustered temporal structure of fitness increasing activity breaks similar to that observed in the field and (3) a remarkable stability of the temporal pattern within a wide range of model parameters.
As long as development time is not included in a model of colony dynamics, temporal variation in the nesting cycle can only be explained by variability in external factors: When resource availability is low (e.g. low worker efficiency rate c) and predation risk is high (high mortality rate 1-q), nests should be closed to protect brood and adults of the colony, otherwise the colony should be active. Without such environmental variation the temporal course of worker numbers would always be monotonically increasing (until the optimal switching time is reached) and thus, the beneficial effect of nest closure (decreased mortality) can never outbalance the effect of productivity loss due to the wastage of time. Thus, nest closure could never increase fitness within the simple concept of temporally constant parameters .
However, as soon as a non negligible time for development of the larvae is taken into account [4, 29], growth is delayed and the temporal course of colony size exhibits increasing as well as decreasing phases (fig. 1). With finite development time of larvae, fitness increasing activity breaks evolve as an emergent property and there is no need to assume external factors like environmental variation. Certainly environmental variation may trigger the appearance of activity breaks too, but our model provides a more general explanation that is in excellent accordance with inter- and intraspecific stability of the observed phenomenon even on a geographical scale of a species' range [31–33].
Weissel et. al (submitted, and unpublished) have shown that the temporal pattern of the colony cycle in the halictid bee L. malachurum depends on soil temperature, but not on resource availability or predation intensity. This result is consistent with our model prediction. We found that the potential temporal structure of active and inactive phases is only determined by development time (in relation to season length) as the main time constant of the system (fig. 5). The close relationship between development time and temperature is well known in many insect species in general [35, 36] and in bees . Constant mortality rates and worker efficiency just determine the occurrence of the potential breaks while temporal position of breaks is affected by the effect of soil temperature on development time. Although the influence of frequency dependent selection selection (e.g. the threat of usurpation by floater females) will tend to synchronise the temporal pattern of all externally driven activity breaks , the observed synchronization of colony activity in aggregations of e.g. L. malachurum is readily explained by the shelter a closed nest provides for the developing brood.
Activity breaks will of course reduce any mortality that is experienced by bees during foraging, e.g. by spiders, birds and wasps. In particular some crabronid wasp species of the genus Cerceris are specialised on hunting halictid bees as provisions for their own progeny . We have observed individuals of C. rybyensis nesting within aggregations of L. malachurum so that they could easily find and paralyze workers returning from foraging trips. There are also conopid flies that wait in the nesting area and parasitize adult bees . Other predators that are excluded by closing the nest are the specialized parasite bee Sphecodes monilicornis that violently enters nests and oviposits into brood cells [41, 42] as well as predatory ants that could destroy the whole nest . Notably, vespid wasps whose more or less open nests would not be much better protected by activity breaks do not show such breaks.
Although seasonal activity patterns of annual halictids with discrete broods have been described quite often [15, 34, 43, 44], exact data on the temporal course of active and inactive phases are scarce. Weissel et al. (submitted) provide data on the seasonal activity state of about 1200 colonies of L. malachurum in northern Bavaria observed during a complete season in 2004 (fig. 6). For an appropriate and reasonable choice of survival rates in our model (Weissel et al., submitted) the number and temporal position of the observed activity breaks in the field are in very good accordance with our model predictions (fig. 2c). Discrepancies between field data and model predictions mainly concern the position and length of the first break. This may be due to the fact that the exact time of colony founding or beginning of egg production is difficult to determine in the field. There is also a number of simplifying assumptions in our model that particularly influence the position and length of the first break. First, we have assumed that egg production starts immediately when the colony has been founded. Second, our model does not account for any differences between founding phase and ergonomic (growth) phase of colony development. During the solitary founding phase mortality will probably be higher than later in the season, when the nest is guarded by a worker . Finally, larval development time might take longer when the temperature is lower early in season. As the model assumes temporally constant mortality and development time the predicted timing of the first break can only be an imprecise estimation and a more detailed model would result in a slightly longer first activity break later in the season.
So far we did not analyse the simultaneous optimization of resource allocation and activity strategies. Although the temporal pattern of fitness increasing breaks turned out to be remarkably stable in our model there will be interactions between the optimal switching point (representing the resource allocation strategy, ) and break emergence (representing the resource accumulation strategy). In addition there might be constraints like egg number or egg laying rate limitation of the queens as has been observed and analysed in bumblebees [4, 45–47] that have only a minor influence on activity patterns, but result in completely different optimal switching points and may even cause a different number of broods within the annual colony cycle .
Our model provides an explanation for the switch from continuous colony growth to reproduction with several discrete broods in social insects . In contrast to the continuous growth model of Oster and Wilson  we provide a model for colonies which show a number of discrete broods per season, separated by distinct activity breaks. There is some evidence that the predicted transition between these two types of reproduction might occur in the field. Populations of Halictus ligatus have been observed exhibiting both strategies in different years . The presence of adult workers within closed colonies or even worker survival for more than one brood also gives empirical indication of transitions between continuous and discrete growth .
Although continuous colony growth is usually taken as the standard type of colony dynamics in eusocial insects, there is no reason to assume that this is the primary state and reproduction with discrete brood periods during a season evolved from that primary state. Our model does not make any predictions about the evolutionary sequence of these two types of colony dynamics. When continuous colony growth is assumed to be the primary type of nest dynamics , then reproduction in discrete broods separated by activity breaks might be assumed to be a derived character. However, the contrary order seems more realistic. Hunt & Amdam (2005) analysed discrete broods as an advanced state of multivoltine reproduction of solitary species. According to their analysis social species evolved most probably from multivoltine solitary species with discrete brood events in the course of the season, as can be found in many solitary bees and wasps. Thus, the first social species most likely had discrete broods rather than continuous colony growth. On this account discrete brood reproduction can be interpreted as an evolutionary link between non-social and advanced eusocial insects like many Apidae . The analysis of multivoltinism within a social context [43, 49] illuminates particular aspects of the first potential transition between non-social reproduction to discrete broods while our modelling approach allows one to determine the necessary conditions for the evolution of continuous reproduction. A more detailed analysis of model parameters (season length, development time, efficiency and mortality) would be useful to determine the optimal reproductive pattern within the whole parameter space.
Activity breaks are not necessarily caused by extrinsic influences. Nests are closed, whenever the resulting loss of productivity is outweighed by the benefit from increased protection of the accumulated brood. For all species and environmental conditions where the decrease in brood mortality caused by nest closure surpasses a critical level it is not necessary to assume any exogenous variability in predation or resource availability as trigger for the alternating phases of activity and inactivity. Nest cycle dynamics itself is sufficient to predict fitness increasing breaks. As productivity loss and survival increase during breaks greatly differ between species, considerable variation in activity patterns is found in annual eusocial insects ranging from species with continuous activity during the whole season to those that show extended periods of inactivity summing up to more than half the season.
Nest dynamics was modelled as a set of coupled difference equations. The model was implemented in the programming language R version 1.7 . The implementation was conducted straight forward as an iteration of the difference equations (1 and 2) with time-dependent parameters (eqns. 3, 4 and 5). Field data collection and analysis is described in the legend of figure 6.
This study was supported by the Deutsche Forschungsgemeinschaft (SFB 554, TP C6). We would like to thank Jürgen Liebig for stimulating discussions on the topics and two anonymous referees for helpful advice and improving the manuscript.
- Macevicz S, Oster G: Modeling Social Insect Populations II: Optimal Reproductive Strategies in Annual Eusocial Insect Colonies. Behavioral Ecology and Sociobiology. 1976, 1: 265-282. 10.1007/BF00300068.View ArticleGoogle Scholar
- Mirmirani M, Oster G: Competition, kin selection, and evolutionary stable strategies. Theoretical Population Biology. 1978, 13: 304-339. 10.1016/0040-5809(78)90049-7.View ArticlePubMedGoogle Scholar
- Brian M, Clarke R, Jones R: A numerical model of an ant society. Journal of Animal Ecology. 1981, 50 (2): 387-405. 10.2307/4062.View ArticleGoogle Scholar
- Beekman M, Lingeman R, Kleijne F, Sabelis M: Optimal timing of the production of sexuals in bumblebee colonies. Entomologia Experimentalis et Applicata. 1998, 88: 147-154. 10.1023/A:1003401628843.View ArticleGoogle Scholar
- Ross K, Mathews R, Eds: The social behavior of wasps. 1991, New York: Comstock, IthacaGoogle Scholar
- Michener C: Social insects: an evolutionary approach to castes and reproduction, Heidelberg: Springer Verlag. 1990, 77-122. chap. Reproduction and caste in social halictine beesGoogle Scholar
- Michener C: The bees of the world. 2000, Baltimore and London: The Johns Hopkins University PressGoogle Scholar
- Miyano S: Life table of colonies and workers in a paper wasp, Polistes chinensis antennalis, in central Japan. Researches on Population Ecology. 1980, 22: 69-88.View ArticleGoogle Scholar
- Archer M: The life history and colonial characteristics of the hornet, Vespa crabro L. (Hym., Vespinae). Entomologist's Monthly Magazine. 1993, 129 (1548–1551): 151-163.Google Scholar
- Steen Z, Schwarz M: Nesting and life cycle of the Australian green carpenter bees Xylocopa (Lestis) aeratus Smith and Xylocopa (Lestis) bombylans (Fabricius) (Hymenoptera: Apidae: Xylocopinae). Australian Journal of Entomology. 2000, 39: 291-300. 10.1046/j.1440-6055.2000.00195.x.View ArticleGoogle Scholar
- Cronin A: Social flexibility in a primitively social allodapine bee (Hymenoptera: Apidae): results of a translocation experiment. Oikos. 2001, 94 (2): 337-343. 10.1034/j.1600-0706.2001.940214.x.View ArticleGoogle Scholar
- Sakagami S: Sozialpolymorphismus bei Insekten. Probleme der Kastenbildung im Tierreich, Stuttgart: Wiss. Verlagsgesellschaft 1974 chap. Sozialstruktur und Polymorphismus bei Furchen- und Schmalbienen (Halictinae). 1974, 257-293. chap. Sozialstruktur und Polymorphismus bei Furchen- und Schmalbienen (Halictinae)Google Scholar
- Yanega D: Demography and sociality in halictine bees (Hymenoptera: Halictidae). The evolution of social behavior in insects and arachnids. Edited by: Choe JC, Crespi BJ. 1997, Cambridge: Cambridge University Press, 293-315.View ArticleGoogle Scholar
- Knerer G, Plateaux-Quenu C: Sur la production continue ou periodique de couvain chez les Halictinae (Insectes, Hym.). Comptes Rendus de l'AcadÈmie des Sciences Paris D. 1967, 264: 651-653.Google Scholar
- Knerer G: Biology and social behavior of bees of the genus Halictus Latreille (Hymenoptera; Halictidae). Zoologische Jahrbücher: Abteilung für Systematik, Ökologie und Geographie der Tiere. 1980, 107 (4): 511-536.Google Scholar
- Knerer G: The biology and social behaviour of Evylaeus malachurus (K.) (Hymenoptera; Halictidae) in different climatic regions of Europe. Zoologische Jahrbücher: Abteilung für Systematik, Ökologie und Geographie der Tiere. 1992, 119 (2): 261-290.Google Scholar
- Duchateau M, Velthuis H: Development and reproductive strategies in Bombus terrestris colonies. Behavior. 1988, 107 (3–4): 186-207.View ArticleGoogle Scholar
- Willmer P, Stone G: Behavioral, ecological, and physiological determinants of the activity patterns of bees. Advances in the Study of Behavior. 2004, 34 (34): 347-466.View ArticleGoogle Scholar
- Willmer P: Thermal constraints on activity patterns in nectar feeding insects. Ecological Entomology. 1983, 8 (4): 455-469.View ArticleGoogle Scholar
- Richards M: Social responses to changing envirnments: reproductive behavior and reproductive options in a primitively eusocial sweat bee. PhD thesis. 1994, York University, North York, OntGoogle Scholar
- Strohm E, Linsenmair E: Temperature dependence of provisioning behaviour and investment allocation in the European beewolf Philanthus triangulum F. Ecological Entomology. 1998, 23: 330-339. 10.1046/j.1365-2311.1998.00129.x.View ArticleGoogle Scholar
- Stone G, Gilbert F, Willmer P, Potts S, Semida F, Zalat S: Windows of opportunity and the temporal structuring of foraging activity in a desert solitary bee. Ecological Entomology. 1999, 24 (2): 208-221. 10.1046/j.1365-2311.1999.00181.x.View ArticleGoogle Scholar
- Vogt D: Cold-hardiness in bumblebees. American Zoologist. 2000, 40 (6): 1246-Google Scholar
- Wcislo W, Arneson L, Roesch K, Gonzalez V, Smith A, Fernandez H: The evolution of nocturnal behaviour in sweat bees, Megalopta genalis and M. ecuadoria (Hymenoptera: Halictidae): an escape from competitors and enemies?. Biological Journal of the Linnean Society. 2004, 83 (3): 377-387. 10.1111/j.1095-8312.2004.00399.x.View ArticleGoogle Scholar
- Warrant E, Kelber A, Gislen A, Greiner B, Ribi W: Nocturnal vision and landmark orientation in a tropical halictid bee. Current Biology. 2004, 14 (15): 1309-1318. 10.1016/j.cub.2004.07.057.View ArticlePubMedGoogle Scholar
- Stone GN: Activity patterns of females of the solitary bee Anthophora plumipes in relation to temperature, nectar supplies and body size. Ecological Entomology. 1994, 19 (2): 177-189.View ArticleGoogle Scholar
- Knerer G: Photoperiod as cue for voltinism and caste regulation in halictine bees. Chemistry and Biology of Social Insects. Edited by: Eder J, Rembold H. 1987, Munich: Verlag Johann Peperny, 305-Google Scholar
- Oster G, Wilson E: Caste and ecology in the social insects. 1978, Princeton, N.J.: Princeton University PressGoogle Scholar
- Perleson A, Mirmirani M, Oster G: Optimal Strategies in Immunology. I. B-Cell Differentiation and Proliferation. Journal of Mathematical Biology. 1976, 3: 325-367.View ArticleGoogle Scholar
- Ihaka R, Gentleman R: R: A language for data analysis and graphics. Journal of Computational and Graphical Statistics. 1996, 5 (3): 299-314. 10.2307/1390807.Google Scholar
- Strohm E, Bordon-Hauser A: Advantages and disadvantages of large colony size in a halictid bee: the queen's perspective. Behavioral Ecology and Sociobiology. 2003, 14 (4): 546-553.View ArticleGoogle Scholar
- Richards M: Evidence for geographic variation in colony social organization in an obligately social sweat bee, Lasioglossum malachurum Kirby (Hymenoptera; Halictidae). Canadian Journal of Zoology. 2000, 78 (7): 1259-1266. 10.1139/cjz-78-7-1259.View ArticleGoogle Scholar
- Wyman L, Richards M: Colony social organization of Lasioglossum malachurum Kirby (Hymenoptera, Halictidae) in southern Greece. Insectes Sociaux. 2003, 50 (3): 201-211. 10.1007/s00040-003-0647-7.View ArticleGoogle Scholar
- Cronin A, Hirata M: Social polymorphism in the sweat bee Lasioglossum (Evylaeus) baleicum (Cockerell) (Hymenoptera, Halictidae) in Hokkaido, northern Japan. Insectes Sociaux. 2003, 50: 379-386. 10.1007/s00040-003-0693-1.View ArticleGoogle Scholar
- Gilbert N, Raworth DA: Insects and temperature: A general theory. The Canadian Entomologist. 1996, 128: 1-13.View ArticleGoogle Scholar
- Trudgill D, Honek A, Li D, van Straalen N: Thermal time – concepts and utility. Annals of Applied Biology. 2005, 146: 1-14.View ArticleGoogle Scholar
- Whitfield G, Richards K: Temperature-dependent development and survival of immature stages of the alfalfa leafcutter bee Megachile rotundata Hymenoptera; Megachilidae). Apidologie. 1992, 23: 11-23.View ArticleGoogle Scholar
- Kaitala V, Smith B, Getz W: Nesting strategies of primitively eusocial bees – a model of nest usurpation during the solitary state of the nesting cycle. Journal of Theoretical Biology. 1990, 144 (4): 445-472.View ArticleGoogle Scholar
- Bohart R, Menke A: Sphecid wasps of the world. 1976, University of California PressGoogle Scholar
- Knerer G: Periodicity and strategy of the parasites of a social bee Evylaeus malachurus (Apoidea: Halictidae). Zoologischer Anzeiger Leipzig. 1973, 190 (1–2): 41-63.Google Scholar
- Legewie H: Zum Problem des tierischen Parasitismus, I. Teil: Die Lebensweise der Schmarotzerbiene Sphecodes monilicornis K. (= subquadratus) Hymenoptera: Apoidea). Zeitschrift für Morphologic und Ökologie der Tiere. 1925, 4 (3): 430-464. 10.1007/BF00408464.View ArticleGoogle Scholar
- Sick M, Ayasse M, Tengo J, Engels W, Lubke G, Francke W: Host-parasite relationships in six species of Sphecodes bees and their Halictid hosts – nest intrusion, intranidal behavior, and Dufours Gland volatiles (Hymenoptera, Halictidae). Journal of Insect Behavior. 1994, 7: 101-117. 10.1007/BF01989830.View ArticleGoogle Scholar
- Packer L, Knerer G: Social evolution and its correlates in bees of the subgenus Evylaeus (Hymenoptera, Halictidae). Behavioral Ecology and Sociobiology. 1985, 17 (2): 143-150.Google Scholar
- Coelho B: The biology of the primitively eusocial Augochloropsis iris (Schrottky, 1902)(Hymenoptera, Halictidae). Insectes Sociaux. 2002, 49 (2): 181-190. 10.1007/s00040-002-8299-6. [English 0020–1812]View ArticleGoogle Scholar
- Rosenheim J: An evolutionary argument for egg limitation. Evolution. 1996, 50 (5): 2089-2094. 10.2307/2410766.View ArticleGoogle Scholar
- Heimpel G, Rosenheim J: Egg limitation in parasitoids: A review of the evidence and a case study. Biological Control. 1998, 11 (2): 160-168. 10.1006/bcon.1997.0587.View ArticleGoogle Scholar
- Casas J, Nisbet R, Swarbrick S, Murdoch W: Eggload dynamics and oviposition rate in a wild population of a parasitic wasp. Journal of Animal Ecology. 2000, 69 (2): 185-193. 10.1046/j.1365-2656.2000.00376.x.View ArticleGoogle Scholar
- Richards M, Packer L: Annual variation in survival and reproduction of the primitively eusocial sweat bee Halictus ligatus (Hymenoptera: Halictidae). Canadian Journal of Zoology. 1995, 73 (5): 933-941. [Article- English]View ArticleGoogle Scholar
- Hunt J, Amdam G: Bivoltinism as an antecedent to eusociality in the paper wasp genus Polistes. Science. 2005, 308: 264-267. 10.1126/science.1109724.PubMed CentralView ArticlePubMedGoogle 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.