Introduction

Phenotypic variation in life-history traits observed among individuals or populations of the same species can be a consequence of genetic differentiation, or may be due to the effects of different environments on the expression of the phenotype (Schmalhausen, 1949; Endler, 1986). The latter source of variation may be merely a biochemical or physiological interaction of the organism with its environment, or may be adaptive if phenotypic variation increases fitness in the environments encountered (Schmalhausen, 1949; Stearns, 1989). A comprehensive and synthetic understanding of variation requires the examination of both environmental and genetic sources, as well as their interaction (Bernardo, 1993; Sinervo and Svensson, 1998; Niewiarowski, 2001).

A particularly important and widespread mediator of phenotypic variation in ectothermic animals is temperature. The thermal environment experienced by an individual during development has a direct effect on final size, with body size generally increasing as temperature decreases and vice versa (Atkinson, 1994; Partridge and French, 1996; Chown and Gaston, 1999). Moreover, evolutionary responses, with larger individuals in low-temperature environments, have been found in the wild (latitudinal and altitudinal size clines; eg Imasheva et al, 1994; James et al, 1997; Billerbeck et al, 2000) and in the laboratory (Partridge et al, 1994, 1995). Likewise, egg size in ectotherms was commonly found to increase in colder regions and colder seasons (Azevedo et al, 1996; Yampolski and Scheiner, 1996; Blanckenhorn, 2000), and insects lay larger eggs at lower temperatures under laboratory conditions (eg Avelar, 1993; Crill et al, 1996; Sheader, 1996; Ernsting and Isaaks, 1997; Blanckenhorn, 2000). The underlying mechanisms shaping this temperature-size rule as well as its adaptive significance are largely unknown (Azevedo et al, 1996; Crill et al, 1996; Blanckenhorn, 2000; Fox and Czesak, 2000).

For the tropical nymphalid butterfly Bicyclus anynana (Butler, 1879), a previous study showed that the prevailing temperature during oviposition induces a plastic response in egg size, with females kept at a lower oviposition temperature laying larger but fewer eggs than those kept at a higher temperature (Fischer et al, in press). By cross-transferring the experimentally manipulated eggs between temperatures, we showed that the larger eggs produced at a lower temperature had a higher hatching success and yielded larger hatchlings with a slightly higher probability of reaching maturity compared to the smaller eggs produced at a higher temperature. Overall, the lower temperature was slightly more detrimental for smaller than for larger eggs. The slight differences in survival across temperatures together with more substantial ones in fecundity indicate that it could pay off to produce fewer but larger offspring (with increased fitness) at a lower temperature, but more and smaller offspring at a higher temperature where the offspring survival was generally high.

This is in agreement with predictions based on theory (Yampolski and Scheiner, 1996). If this were indeed the case and if, in consequence, the temperature-mediated plasticity in egg size were adaptive (in this species), one might expect that there is genetic variation in the plastic response to temperature for natural selection to act upon (as long as the variation has not been completely eroded; Via and Lande, 1985; Scheiner, 1993).

To address this issue, we use here a full-sib design with sibs being divided between different larval and pupal/egg-laying temperatures, enabling us to draw on several of the questions raised above. First, we will investigate the genetic (family) and environmental (temperature) effects on egg size, as well as potential interactions between the two. The latter would be indicative of the existence of genetic variation in the plastic response of egg size to temperature, and hence that there is potential for short-term evolutionary change. Second, we will explore the effects of larval (including early pupal) rearing temperature and adult (including late pupal) temperature on egg size. Since plastic responses in egg size can be induced by rearing larvae and maintaining the resulting adults at different temperatures, as well as exclusively by different oviposition temperatures, this experiment will reveal information on the timing and hence the sensitive period for the induction of a plastic response. Third, we will estimate the heritability of egg size by means of variance component estimates (using the full-sib families) and parent–offspring regression.

Material and methods

Study organism and experimental population

B. anynana is a tropical, fruit-feeding butterfly ranging from southern Africa to Ethiopia (Larsen, 1991). The species exhibits striking phenotypic plasticity (two seasonal morphs), which is thought to function as an adaptation to alternative wet–dry seasonal environments and the associated changes in resting background and predation (Brakefield and Larsen, 1984; Brakefield and Reitsma, 1991; Windig et al, 1994). A laboratory stock population of B. anynana was established at Leiden University in 1988 from over 80 gravid females collected at a single locality in Malawi. Several hundred adults are reared in each generation, maintaining high levels of heterozygosity (Saccheri and Bruford, 1993). Butterflies from the stock population were used for this study.

Experimental design

Butterflies of the parental generation were reared in a climate room at 27°C, high humidity, and a photoperiod of L12:D12 till pupation. Larvae were fed on young maize plants. On day 2 after pupation, female pupae were randomly divided among low (20°C, forming group A) and high (27°C, forming group B) temperatures (with high humidity and L12:D12 throughout). The temperatures chosen are similar to the ones experienced by the butterflies during the dry and wet seasons, respectively, in the field (Brakefield and Mazzotta, 1995; Brakefield and Reitsma, 1991). Following adult eclosion, females were individually marked and mated (at the respective temperature, they were assigned to as pupae) to randomly chosen stock males. In order to account for differences in physiological age, females were given 1 day prior to mating (ie separated from males) at 27°C, and were then put together with males for 2 days for mating, whereas at 20°C females had a premating period of 2 days, followed by 3 days for mating. The use of different premating and mating periods was motivated by the fact that the larval development time is about half as long at 27°C, as compared to 20°C.

Following mating (ie on days 4 and 6 of adult life at 27 and 20°C, respectively), females were placed individually in translucent plastic pots (1 l) containing a fresh cutting of maize to oviposit on. Eggs were collected for 1 (27°C) and 2 (20°C) days, respectively, and subsequently measured (see below). Hence, all eggs included in this analysis were the first ones laid within the females' adult life span, collected within one or two 12 h light periods, thus effectively controlling for any confounding effects of female age (eg Karlsson and Wiklund, 1984; Brakefield et al, 1994; Braby and Jones, 1995). Over the following days, two groups of eggs were collected per female and again divided between both temperatures for rearing (see Figure 1 for illustration of the experimental design). The offspring was reared in families till pupation, with female pupae afterwards being once again randomly divided between 20 and 27°C. Division of pupae took place on day 2 after pupation for individuals being reared at 27°C, and on day 4 for those reared at 20°C. This design results in eight treatment groups with different thermal histories (ie 20–20–20 for pupal/adult temperature of the mother–larval/pupal temperature–pupal/adult temperature of the daughter, 20–20–27, 20–27–20, 20–27–27, all of which belong to group A; and 27–20–20, 27–20–27, 27–27–20, 27–27–27, belonging to group B). F1 female butterflies were marked, mated and individually set up for egg laying in the way outlined above. Throughout the experiment, butterflies had access to moist banana for adult feeding.

Figure 1
figure 1

Schematic overview of experimental protocol. Parental generation females were reared as larvae at 27°C, and on day 2 after pupation randomly divided between 20 and 27°C. Per female butterfly, two groups of eggs were collected at the respective temperatures they were assigned to, and subsequently divided between 20 and 27°C for larval rearing. F1 generation female pupae were once again divided between 20 and 27°C, resulting in eight treatment groups with different thermal histories. The offspring of the parental generation females transferred as pupae to 20 or 27°C form groups A and B, respectively.

Data collection and analysis

For all F1 females, we measured pupal mass (measured 1 day after pupation), larval and pupal development time, and egg size. As the eggs of B. anynana are nearly perfectly spherical, egg size was measured as the cross-sectional area (mm2), using a digital camera (Leica DC200) connected to a binocular microscope. The resulting images were analysed using Scion Image public software (Scion Corporation, 2000). Tight correlations between egg area (applying image analysis) and egg mass as well as hatchling size confirm that this method is highly reliable to measure egg size in B. anynana (Fischer et al, 2002). The mean of about 10 eggs was used to calculate the egg size for individual females. Previous experiments showed that the means did not change substantially above a critical minimum number of 7–8 eggs (data not shown).

Statistics

Life-history data were analyzed separately for the offspring of groups A and B (ie for parental generation females transferred as pupae to either 20 or 27°C). This was done to investigate the family effects, and because the larvae reared at 20°C had to be divided over two climate cells due to space limitations (ie for 20°C larval rearing, the offspring of group A were reared in a different cell than those of group B). Statistical analyses used three-way analyses of variance (ANOVAs) with larval (plus early pupal; hereafter LARVAL) temperature and (late pupal plus) adult (hereafter ADULT) temperature as fixed factors and family as random factor. By removing the nonsignificant interaction terms, the minimum adequate models were constructed. Only those families with at least two offspring in each of the four temperature combination groups were included in the analyses. Heritability of egg size was estimated by means of parent–offspring regression (by regressing the mean egg size of full-sib sisters, weighted by the number of females in each full-sib family, against the mean egg size of the mother) and variance component estimates of the ANOVAs. All statistical tests were performed using JMP version 4.02.

Results

Parental generation

The randomly selected pupae from the stock population, when transferred to either high or low temperature, did not differ in pupal mass (198.8±18.8 mg vs 202.2±14.9 mg; t56=0.76, P=0.45). ADULT temperature had a marginally significant effect on egg size (t64=2.1, P=0.04), with females kept at 20°C laying larger eggs (0.723±0.043 mm2) than those at 27°C (0.703±0.034 mm2).

F1 generation

In both sets of families (ie offspring of groups A and B), genetic (ie families) and environmental (ie LARVAL temperature) sources affected egg size (Tables 1 and 2; Figure 2). Females reared as larvae at 20°C laid significantly larger eggs than their sisters at 27°C (Table 1). Although there seems to be some variation among families in the plastic response to temperature (Figure 2), this was not substantiated statistically (interaction terms between family and larval temperature, explaining 2–3% of the variance, were nonsignificant and hence removed from the models presented in Table 2; in the full model, the P-values were P=0.54, F1,10=0.94 for group A and P=0.21, F1,15=1.53 for group B). Different ADULT temperatures, however, did not result in different-sized eggs (Table 1, 2). Significant interactions between LARVAL and ADULT temperature indicate that the effect of ADULT temperature depends on LARVAL temperature. When reared as larvae at 20°C, egg size tended to increase at the higher rearing temperature, while at 27°C the situation is reversed. Family effects explained 19–23% (h2) of the total variance in egg size, while temperature effects accounted for 10.1 and 24.4% in groups A and B, respectively (Table 3). Heritabilities, as estimated by means of parent–offspring regression, reached somewhat higher values with 0.39 and 0.40 for groups A and B, respectively (Table 4).

Table 1 Egg size, larval and pupal development times and pupal mass for groups of Bicyclus anynana females with different thermal histories
Table 2 Results of three-way analyses of variance (ANOVAs) for the effects of family (random factor), LARVAL, and ADULT temperature (fixed) on female life-history traits in Bicyclus anynana
Figure 2
figure 2

Reaction norms for egg size in relation to LARVAL rearing temperature in the butterfly, Bicyclus anynana. Each line connects the mean egg size of females reared at 20°C with the mean egg size of their sisters reared at 27°C, and thus represents a single split-family. The two graphs refer to parental generation females transferred as pupae to 20°C (a) and to 27°C (b). As offspring, ADULT temperature did not affect egg size, these treatments were pooled.

Table 3 Proportion (in %) of the total variance in egg size explained by different sources for the butterfly Bicyclus anynana
Table 4 Heritability (h2±1 SE) of egg size in the butterfly Bicyclus anynana

As we have found an effect of ADULT temperature on egg size in the parental, but not in the F1 generation (see above), we repeated this part of the experiment by again dividing randomly chosen stock females of a later generation reared at 27°C between 20 and 27°C. Here, subgroups of females were divided as prepupae, or on days 2 or 4 after pupation between both temperatures. Throughout, the groups transferred as prepupae or pupae to 20°C laid larger eggs than the control groups kept at 27°C (prepupae: 0.692 vs 0.671 mm2, n=13 each; pupae day 2: 0.693 vs 0.670 mm2, n=20 each; pupae day 4: 0.700 vs 0.670 mm2, n=16 each). A two-way ANOVA confirmed the temperature effect (F1,92=7.8, P=0.006), whereas neither the timing of division nor its interaction with temperature was significant (time: F2,92=0.1, P=0.91; interaction term: F2,92=0.1, P=0.89).

Regarding the other life-history traits, a significant family effect was only present in the larval development time for the offspring of group A (Table 2). Either high LARVAL or ADULT temperatures caused predictable decreases in development times (Table 1, 2; a marginally significant effect of ADULT temperature on larval development time is presumably a chance effect of allocation to treatments). However, LARVAL temperature did not consistently influence pupal weight; only the offspring within group B attained higher pupal weights when reared at 27°C. Significant interactions between LARVAL and ADULT temperature for pupal times are caused by the experimental design, as the individuals were divided among temperatures during the pupal stage.

Discussion

Temperature-mediated plasticity in egg size

As already found in other arthropods (eg Azevedo et al, 1996; Crill et al, 1996; Sheader, 1996; Yampolski and Scheiner, 1996), egg size increases at a lower temperature in the butterfly B. anynana. Such temperature-induced plasticity that results in larger offspring at lower temperatures seems to be a near universal rule in ectotherms (Azevedo et al, 1996; Yampolski and Scheiner, 1996; Fox and Czesak, 2000). However, as most previous studies did not distinguish between the effects of larval versus adult temperature, much less is known about the sensitive periods for the induction of the plastic response (Fox and Czesak, 2000). Likewise, the underlying mechanisms remain obscure (Azevedo et al, 1996; Blanckenhorn, 2000; Ernsting and Isaaks, 2000; Fox and Czesak, 2000). From previous experiments that disentangled developmental plasticity from acclimation in the adult stage, we do know that low oviposition temperatures result in larger eggs (eg Avelar, 1993; Huey et al, 1995; Ernsting and Isaaks, 1997; Blanckenhorn, 2000; Fischer et al, in press). However, we could not find a single study investigating exclusively the effects of different temperatures during preadult stages on egg size (by rearing larvae at different temperatures, but maintaining adults in a common environment).

Since egg maturation in butterflies starts in the pupal or early adult phase (eg Boggs, 1997a, 1997b; Ramaswamy et al, 1997), we predicted that lower temperatures experienced from this stage onwards would cause an increase in egg size. This prediction is supported by data from parental generation females and from the additional experiment in which females were transferred to a lower temperature as prepupae, or on days 2 or 4 after pupation. Here, females reared as larvae at 27°C, but reared as pupae (at least partly) at 20°C, laid larger eggs than the control groups kept at 27°C throughout. However, the results based on the families in the main experiment did not provide similarly conclusive support. Interestingly, for the larvae reared at 27°C and divided as pupae among high and low temperatures, there is at least a trend in the predicted direction, whereas this is not the case for the animals reared at 20°C (resulting in a significant interaction between LARVAL and ADULT temperature; Tables 1, 2). Consequently, all the evidence in support of our prediction comes from animals reared as larvae at 27°C, whereas those developing at 20°C did not show a clear plastic response (Table 1). This could be due to a relatively earlier start of egg maturation (and hence determination of egg size) in these slowly developing animals. Note that, as we used the first eggs laid by individual females in our analysis, we cannot detect any responses to temperature once the size of these early eggs has been determined.

Besides these effects of pupal/adult temperatures on egg size, we have found a clear response to the temperature experienced during larval growth and early pupal development (Table 2). Regardless of late pupal and adult temperature, females reared as larvae at 20°C laid larger eggs than their sisters reared at 27°C. Note that the individuals transferred from 20 to 27°C spent ca. 8 days prior to collection of eggs at the higher temperature, and the ones transferred from 27 to 20°C about 14 days at the lower temperature. To our knowledge, this, for the first time, shows that temperature differences experienced that early in development induce a plastic response in egg size. These data suggest that the sensitive period for the determination of the size of first eggs largely occurs during the pupal stage. Here, low temperatures during the first days only, as well as from the last days onwards (at least when reared at a high temperature), will result in increased egg size. Once the butterfly had enclosed from the pupa, however, the size of the first eggs seems to be fixed. This is not in contradiction to a response to prevailing temperature conditions during oviposition, as has been found in a number of species (see above). However, this is a gradual process not affecting the very first eggs laid by a given female (Fischer et al, in press and in prep.). These findings are particularly interesting as B. anynana is one of those species not having any chorionated eggs at eclosion (ie egg maturation starts in the early adult stage), and as it is not able to lay any eggs without additional nutrients from adult income (unpublished observations). Note that the differences in egg size found in our study can neither be explained as a correlated response to changes in body size nor as the result of delayed oviposition at lower temperatures (see also Ernsting and Isaaks, 1997; Fox and Czesak, 2000).

Heritability of egg size and genotype–environment interactions

The parent–offspring regressions (h2 ca. 0.4) as well as the variance component estimates (h2 ca. 0.2) show that heritability of egg size is low to moderate, as could be expected for a life-history trait (Roff, 1997). These values are also in broad agreement with selection experiments using B. anynana, revealing that about 40% of the variation in egg size is due to additive genetic variance (Fischer et al, unpublished data). Generally, there are surprisingly few data available on genetic variation in egg size within insect populations (Fox and Czesak, 2000). In two seed beetles, egg size is highly heritable and the estimates range between 0.22 and 0.91, varying between populations and host plants (Fox et al, 1999). Also laboratory selection experiments have demonstrated that egg size is heritable in insects (Harvey, 1983; Schwarzkopf et al, 1999).

Environmental and genetic sources tended to affect egg size to a similar extent, emphasising the importance of considering both types of variation (eg Sinervo and Svensson, 1998; Niewiarowski, 2001). However, we were not able to demonstrate the existence of genotype–environment interactions, although there seems to be some variation in the slopes of reaction norms (Figure 2). As heritabilities of plasticity are usually lower than trait heritabilities (Scheiner, 1993; Scheiner and Yampolski, 1998), they are relatively difficult to detect and we cannot rule out that the absence of genetic variation in the plastic response is due to insufficient statistical power (ie small sample sizes) in our analysis. Based on our present data, power analyses indicated that the chance of detecting a significant interaction was only about 6%. The lowest number needed to find a significant interaction at the 5% level is 2000–3000 individuals. Apart from insufficient statistical power, our negative results may indicate that temperature-mediated egg size plasticity is caused by physiological constraints, or that there is a strong directional selection favoring a direct response to temperature, thus eroding variation in the plastic response. Both interpretations are in line with the fact that the production of bigger eggs at lower temperatures is widespread in insects (see above).

Although temperature-mediated plasticity in egg and body size is a widespread phenomenon (see above), an adjustment of egg size to the prevailing temperature conditions could be of particular importance for tropical B. anynana in its seasonal environment (Fischer et al, in press). Here, reproduction is largely confined to the wet season when oviposition plants are abundantly available, and where 2–3 generations occur. Towards the end of the wet season, however, there is a marked decrease in temperature, starting several weeks before larval food plants dry out completely and reproduction ceases (Windig et al, 1994; Brakefield, 1997). Females ovipositing during this phase may well benefit from laying larger eggs with enhanced fitness.

Pupal weight and development times

Regarding pupal weight and development time, there was no detectable family effect, with only one exception (larval development time for the offspring of group A). This suggests low heritabilities for these important life-history traits (cf. Roff, 1997). The realized heritabilities in B. anynana are indeed low for development time (0.06, Zijlstra et al, unpublished data); however, they are moderate for pupal weight (ca. 0.25, Zwaan et al, unpublished data). The effects of temperature on development time follow the well-known pattern of a faster development at higher temperatures (Honek and Kocourek, 1990). Pupal weight, however, did not increase at lower temperature, comprising an exception from the above-mentioned temperature-size rule (Atkinson, 1994; Partridge and French, 1996). This may indicate that a constant temperature of 20°C is already suboptimal for growth and development for this tropical butterfly (Brakefield, 1997).

In summary, this study shows that developmental temperature can induce a plastic response in egg size even at a surprisingly early developmental stage (ie larval plus early pupal), long before first eggs actually mature. Further, our results stress the importance of considering both environmental and genetic sources of variation for this type of trait. However, further research is required to demonstrate whether genetic variation occurs for the plastic response of egg size to temperature.