The relationship between eyespot shape and wing shape in the butterfly Bicyclus anynana: A genetic and morphometrical approach

The African butterfly, Bicyclus anynana, normally possesses circular eyespots on its wings. Artificial selection lines, which express ellipsoidal eyespots on the dorsal surface of the forewing, were used to investigate correlated changes in wing shape. Morphometric analysis of linear wing measurements and wing scale counts provided evidence that eyespot shape was correlated with localised shape changes in the corresponding wing-cell, with overall shape changes in the wing, and with the density/arrangement of scales around the eyespot area.     

Reaction norms and the genetic basis of phenotypic plasticity in the wing pattern of the butterfly Bicyclus anynana

The genetic basis of the dry-wet season polyphenism of wing pattern in response to temperature shown by Bicyclus anynana was studied, using a split-family design over four temperatures. Reaction norms crossed, but were only linear in the three highest temperatures, and only when larval development time was used as the environmental axis. Significant full-sib additive variances (V_A) and heritabilities (h²) for plasticity were found using slopes of reaction norms in a bootstrap procedure. Heritabilities were lower in intermediate temperatures, mainly due to differences in the residual variances (V_R). There was no clear trend in V_A across temperatures, contrary to the expectation that V_A would have been depleted by natural selection at the extreme temperatures and not depleted at the intermediate temperatures which occur less frequently in the field. Unpredictability in the onset of the following season at intermediate temperatures might lead to selection for diverse flresponses resulting in relatively high V_Rs. Theoretical models linking reaction norms to genetic parameters in separate environments were difficult to apply in this study, particularly because they are based on the assumption that V_Rs are constant. However the reaction norm approach combined with quantitative genetics provided a valuable insight into the evolution of the observed polyphenism.     

Temperature and humidity acclimation increase desiccation resistance in the butterfly Bicyclus anynana

Desiccation resistance, that is, the ability to reduce water loss, is an ecologically important trait relevant to all terrestrial organisms, which may constrain species distributions. Nevertheless, relatively few studies have investigated plastic capacities in desiccation resistance. We here investigate plastic responses in body mass change, used as a proxy of desiccation resistance, to variation in temperature and relative humidity in the tropical butterfly Bicyclus anynana (Butler) (Lepidoptera: Nymphalidae). Our results indicate that butterflies acclimated to a higher (27 °C) compared with a lower temperature (18 °C) and a lower (50%) compared with a higher (90%) relative humidity displayed a decreased loss of body mass, and therefore likely a loss of body water (27 °C: 11%, 18 °C: 15%; 50% r.h.: 14%, 90% r.h.: 18%). Thus, mass loss was reduced under conditions indicating increased desiccation risk, suggesting adaptive phenotypic plasticity. Effects were most pronounced during the first 24 h after acclimation, indicating quick and transient responses to environmental conditions. As anthropogenic climate change is predicted to increase the magnitude and frequency of heat and drought periods, we argue that more studies on plastic capacities in traits relating to desiccation resistance are needed to better understand species responses.     

Selection for abdominal tergite pigmentation and correlated responses in the trident: a case study in Drosophila melanogaster

In Drosophila melanogaster, abdominal tergite pigmentation and the appearance of a trident‐shaped thoracic pattern exhibit similar biogeographical variation and sensitivity to temperature. These pigmentation traits may be under common selection pressure in natural populations or may be genetically correlated. To investigate the nature of this interaction, replicated populations of D. melanogaster were selected for increased or decreased melanization of the abdominal tergites for 40 generations. Selection for abdominal tergite pigmentation leads to correlated changes in trident formation. Although selection was performed only on female flies, male pigmentation also responded to selection. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 106 , 287–294.     

Pedigree analysis on small laboratory populations of the butterfly Bicyclus anynana: The effects of selection on inbreeding and fitness

The effect of small population size and gene flow on the rate ofinbreeding and loss in fitness in Bicyclus anynana populationswas quantified by means of a pedigree analysis. Laboratorymetapopulations each consisted of four subpopulations with breeding sizeof N = 6 or N = 12 and migration rate of m = 0 or m= 0.33. Pedigrees were established by individually marking about35,000 butterflies. The increase in inbreeding coefficients(F-coefficients) over time was compared to that of simulated populationswith similar N and m. In the seventh generation, the level of inbreedingin larger subpopulations did not deviate significantly from the expectedvalues, but smaller subpopulations were less inbred than expected.Individuals in the small populations still showed considerableinbreeding depression, indicating that only a small proportion of therecessive deleterious alleles had been purged by selection. Two opposingprocesses potentially affected the rate of inbreeding and fitness: (1)Inbreeding depression increased the variance in family size and reducedthe effective population size. This will accelerate the rate ofinbreeding and is expected to selectively purge deleterious recessivealleles. (2) Variance in reproductive success of families was reducedbecause individuals which had a large number of siblings in thepopulation were more likely to mate with a full-sib than individualswith a smaller number of siblings. Subsequent inbreeding depressionreduced the number of viable offspring produced by these full-sibmatings. As a consequence, natural selection purged only some of thedeleterious alleles from the butterfly populations during sevengenerations with inbreeding. These findings emphasise the potentialproblems of using only small numbers of breeding individuals (N10) incaptive populations for conservation purposes.     

Reduced Mating Vigor in Selection Lines of the Butterfly Bicyclus anynana

Despite its acknowledged value for testing evolutionary theory, artificial selection is vulnerable to artifacts that are difficult to control. Here, we show that lines of the butterfly Bicyclus anynana, being artificially selected for differences in egg size, show considerably reduced mating vigor and male mating success than control lines, which is attributed to presumed differences in the levels of inbreeding. Besides selection line, neither forewing length, body dry mass or fat content of males had a detectable impact on male mating performance. Because artificial selection may (e.g. due to inbreeding) yield inadvertent effects even on traits that seem not to be related to the target trait, their results should be interpreted with caution.     

Simultaneous selection on two fitness-related traits in the butterfly Bicyclus anynana

Abstract. Theory about the role of constraints in evolution is abundant, but few empirical data exist to describe the consequences a bias in phenotypic variation has for micro evolution. Responses to natural selection can be severely hampered by a genetic correlation among a suite of traits. Constraints can be studied using antagonistic selection experiments, that is, two-trait selection in opposition to this correlation. The two traits studied here were development time and wing pattern (eyespot size) in the butterfly Bicyclus anynana , both of which have a clear adaptive significance. Rates of response were higher for eyespot size than for development time, but were independent of the concurrent selection (either in the same direction as the correlation or perpendicular to it). Regimes differed in both traits in all directions after 11 generations of selection. The uncoupling lines had higher relative responses than the synergistic lines in development time and equal relative responses in eyespot size. The patterns for eyespot size (reaction norms) were consistent across different rearing temperatures. Differences in lines selected for fast and slow development time were more pronounced at lower temperatures, irrespective of the direction of joint wing pattern selection. Furthermore, correlated responses in pupal weight and growth rate were observed; lines selected for a slower development had higher pupal weights, especially at lower temperatures. The response of the uncoupling lines was not hampered by a lack of selectable genetic variation, and the relative response in the development time was larger than expected based on response in the coupled direction and quantitative genetic predictions. This suggests that the structure of the genetic architecture does not constrain the short-term, independent evolution of both wing pattern and development time.     

Ecdysteroids control eyespot size and wing color pattern in the polyphenic butterfly Bicyclus anynana (Lepidoptera: Satyridae)

The strongly polyphenic African butterfly, Bicyclus anynana, shows conspicuous ventral eyespots and a transverse band in the wet-season form and small eyespots and no band in the dryseason form. These forms are produced when larvae are reared at high and low temperatures, respectively. Truncation selection was applied to a stock population (UNSELECTED-LINE) to produce lines which, at a constant intermediate temperature of 20 °C, always produced the dry season form (LOW-LINE) and the wet-season form (HIGH-LINE) in addition to a line of fast development (FAST-LINE). A relationship between wing pattern and development time was apparent: the FAST-LINE displayed larger eyespots and HIGH-LINE pupae developed faster (mean = 12.5 days) than LOW-LINE pupae (14.1 days). Differences were found among the lines in ecdysteroid titers after pupation. Hemolymph ecdysteroids in HIGH-LINE pupae increased earlier and reached twice the level of those in LOW-LINE pupae during the first 3 days after pupation. FAST-LINE pupae developed faster (11.7 days) than UNSELECTED-LINE pupae (12.8 days) and ecdysteroids in the FAST-LINE increased more quickly and reached higher levels. In the four LINES, ecdysteroid titers in 3 day old pupae were in the order UNSELECT ≈ LOW ⪡ FAST ⪡ HIGH. Thereafter the titers overlapped.An injection of 20-hydroxyecdysone (20E) inhibited pupal development at a dose between 2.5 and 5 μg when it was injected into pupae within 24 h after pupation. At lower doses (0.25–0.5 μg 20E) 22–100% of the pupae in different experimental groups in the LOW- as well as in the HIGH-LINE developed successfully. The pupal stage was significantly shortened, especially in the LOW-LINE. Additionally, 0.25 and 0.5 μg 20E injected into 0–12 h LOW-LINE pupae shifted the wing color pattern towards the wet season form: eyespots increased in size and the transverse wing band appeared in the more conspicuous pattern characteristic of the wet season form. The results demonstrate that ecdysteroids appearing early in the young pupa produce the wet season form of the wings. The same hormonal system mediates both developmental time and wing pattern determination.     

Effective population size,reproductive success and sperm precedence in the butterfly,Bicyclus anynana,in captivity

A pedigree approach is used to estimate the effective population size yn two population cages of the butterfly, Bicyclus anynana. Each cage was founded with 54 individually marked adults of each sex. Matings were recorded over a 3‐day period. Eggs were then collected from each female over a similar period before the numbers of hatching larvae were counted to assess progeny number. The males showed a higher variance in reproductive success than the females. Since about one‐quarter of all females mated more than once, we also examined the pattern of sperm precedence using molecular markers or, in separate crossing experiments, wing pattern mutants. Both instances of complete first and last male sperm precedence, as well as of sperm mixing, were found. In some crosses a ‘leakiness’ was found in which some of the early eggs laid by a female were fertilized by a male partner which was subsequently completely unsuccessful. However, the estimates of effective population size were largely unaffected by the pattern of sperm precedence. Estimates for N_e : N in each cage were close to 0.60. The possibility of obtaining comparable estimates in selected natural populations of butterflies is discussed.     

Differing selection in alternative competitive environments: shade-avoidance responses and germination timing

Under competitive conditions, stem elongation in plants is thought to enhance fitness by increasing light interception. However, the onset of competition should vary with the species of competitor due to interspecific differences in timing of emergence and plant growth form. The fitness benefits of elongation may therefore depend on the timing of this plastic response. Phenotypic selection analyses and path analysis were used to evaluate selection acting on stem elongation at early and late life-history stages and the combination of germination timing and elongation in an annual plant. Velvetleaf (Abutilon theophrasti) were raised in one of three environments experienced by natural populations (cornfields; soybean fields; and disturbed, weedy sites). Due to the rapid growth rate and high density of plants in disturbed areas, selection to increase seedling-stage elongation was expected in weedy sites. Due to the wide spacing of crop plants, competition for light is initially low in cultivated fields, but intensifies as the season progresses. Selection for increased elongation at later nodes was expected in soybean fields because velvetleaf can often overtop soy and thereby increase leaf exposure. In contrast, selection against late elongation was expected in cornfields because velvetleaf are incapable of overtopping corn. Individuals that elongate would experience the carbon cost of allocating to structural tissue, but fail to experience a carbon return through increased light interception. The phenotypic selection analyses were consistent with these predictions and therefore support the role of stem elongation as an adaptation to interspecific competition. Selection also acted on the combination of germination timing and elongation. In the weedy environment, early emergence in conjunction with enhanced stem elongation conveyed the highest fitness. Reduced elongation was favored among individuals that emerged late, potentially because these individuals were unable to overtop neighbors. The results of this study demonstrate that the timing of stem elongation strongly affects competitive success. Environments that differ in the timing of competition for light select for elongation at different life-history stages, and this selection depends on the timing of emergence.     

Direct and indirect responses to selection on pollen size in Brassica rapa L.

Pollen size varies little within angiosperm species, but differs extensively between species, suggesting the action of strong selection. Nevertheless, the potential for genetic responses of pollen size to selection, as determined by additive genetic variance and genetic correlations with other floral traits, has received little attention. To assess this potential, we subjected Brassica rapa to artificial selection for large and small pollen during three generations. This selection caused significant divergence in pollen diameter, with additive genetic effects accounting for over 30% of the observed phenotypic variation in pollen size. Such heritable genetic variation suggests that natural selection could effect evolutionary change in this trait. Selection on pollen size also elicited correlated responses in pollen number (–), flower size (+), style length (+), and ovule number (+), suggesting that pollen size cannot evolve independently. The correlated responses of pollen number, flower size and ovule number probably reflect the genetically determined and physically constrained pattern of resource allocation in B. rapa. In contrast, the positive correlation between pollen size and style length may represent a widespread gametic‐phase disequilibrium in angiosperms that arises from nonrandom fertilization success of large pollen in pistils with long styles.     

Inbreeding depression and genetic load in laboratory metapopulations of the butterfly Bicyclus anynana

Abstract.— We investigated the effects of inbreeding on various fitness components and their genetic load in laboratory metapopulations of the butterfly Bicyclus anynana . Six metapopulations each consisted of four subpopulations with breeding population sizes of N = 6 or N = 12 and migration rate of m = 0 or m = 0.33. Metapopulations were maintained for seven generations during which coancestries and pedigrees were established. Individual inbreeding coefficients at the F₇ were calculated and ranged between 0.01 and 0.51. Even though considerable purging had occurred during inbreeding, the genetic load remained higher than that of many outbreeding species: approximately two lethal equivalents were detected for egg sterility, one for zygote survival, one for juvenile survival, and one for longevity. Severe inbreeding depression occurred after seven generations of inbreeding, which jeopardized the metapopulation survival. This finding suggests that the purging of genetic load by intentional inbreeding cannot be recommended for the genetic conservation of species with a high number of lethal.     

Measuring probabilistic reaction norms for age and size at maturation

We present a new probabilistic concept of reaction norms for age and size at maturation that is applicable when observations are carried out at discrete time intervals. This approach can also be used to estimate reaction norms for age and size at metamorphosis or at other ontogenetic transitions. Such estimations are critical for understanding phenotypic plasticity and life-history changes in variable environments, assessing genetic changes in the presence of phenotypic plasticity, and calibrating size- and age-structured population models. We show that previous approaches to this problem, based on regressing size against age at maturation, give results that are systematically biased when compared to the probabilistic reaction norms. The bias can be substantial and is likely to lead to qualitatively incorrect conclusions; it is caused by failing to account for the probabilistic nature of the maturation process. We explain why, instead, robust estimations of maturation reaction norms should be based on logistic regression or on other statistical models that treat the probability of maturing as a dependent variable. We demonstrate the utility of our approach with two examples. First, the analysis of data generated for a known reaction norm highlights some crucial limitations of previous approaches. Second, application to the northeast arctic cod (Gadus morhua) illustrates how our approach can be used to shed new light on existing real-world data.     

Developmental and genetic mechanisms for evolutionary diversification of serial repeats: eyespot size in Bicyclus anynana butterflies

Serially repeated pattern elements on butterfly wings offer the opportunity for integrating genetic, developmental, and functional aspects towards understanding morphological diversification and the evolution of individuality. We use captive populations of Bicyclus anynana butterflies, an emerging model in evolutionary developmental biology, to explore the genetic and developmental basis of compartmentalized changes in eyespot patterns. There is much variation for different aspects of eyespot morphology, and knowledge about the genetic pathways and developmental processes involved in eyespot formation. Also, despite the strong correlations across all eyespots in one butterfly, B. anynana shows great potential for independent changes in the size of individual eyespots. It is, however, unclear to what extent the genetic and developmental processes underlying eyespot formation change in a localized manner to enable such individualization. We use micromanipulations of developing wings to dissect the contribution of different components of eyespot development to quantitative differences in eyespot size on one wing surface. Reciprocal transplants of presumptive eyespot foci between artificial selection lines and controls suggest that while localized antagonistic changes in eyespot size rely mostly on localized changes in focal signal strength, concerted changes depend greatly on epidermal response sensitivities. This potentially reflects differences between the signal-response components of eyespot formation in the degrees of compartmentalization and/or the temporal pattern of selection. We also report on the phenotypic analysis of a number of mutant stocks demonstrating how single alleles can affect different eyespots in concert or independently, and thus contribute to the individualization of serially repeated traits.     

Fitness costs associated with different frequencies and magnitudes of temperature change in the butterfly Bicyclus anynana

Plastic responses to changes in environmental conditions are ubiquitous and typically highly effective, but are predicted to incur costs. We here investigate the effects of different frequencies and magnitudes of temperature change in the tropical butterfly Bicyclus anynana, considering developmental (Experiment 1) and adult stage plasticity (Experiment 2). We predicted negative effects of more frequent temperature changes on development, immune function and/or reproduction. Results from Experiment 1 showed that repeated temperature changes during development, if involving large amplitudes, negatively affect larval time, larval growth rate and pupal mass, while adult traits remained unaffected. However, results from treatment groups with smaller temperature amplitudes yielded no clear patterns. In Experiment 2 prolonged but not repeated exposure to 39 °C increased heat tolerance, potentially reflecting costs of repeatedly activating emergency responses. At the same time fecundity was more strongly reduced in the group with prolonged heat stress, suggesting a trade-off between heat tolerance and reproduction. Clear effects were restricted to conditions involving large temperature amplitudes or high temperatures.     

Matching field and laboratory environments: effects of neglecting daily temperature variation on insect reaction norms

The tropical butterfly, Bicyclus anynana, exhibits seasonal polyphenism. The wet season form has large eyespots and a pale band while these characters are much less conspicuous or absent in the dry season form. This plasticity is induced in the laboratory by use of a standard series of constant temperatures in the larval stage yielding a continuous norm of reaction. Butterflies in this study were reared from hatchling larvae in seven regimes which differed with respect to thermoperiod or photoperiod. The effect of rearing treatment on the phenotypic plasticity of the adult wing pattern, on life history traits and on larval feeding rhythms was investigated. Photoperiod had little effect except that constant light produced a higher mortality and tended to produce a longer development time. Thermoperiod had a major effect on the life history traits in comparison to a constant temperature regime with the same daily mean: development time was shorter with higher growth rates. The faster development was associated with a substantial shift in the wing pattern towards the wet season form. Larvae feed mostly at night both under constant and thermoperiod (cool nights) conditions. The results are discussed with respect to the necessity of matching field and laboratory environments in studies of norms of reaction or of life history traits where the adaptive significance of the variation is important. Fluctuating conditions in nature, especially with respect to thermoperiod, must be taken into account.     

Evolutionary significance of genotypic variation in developmental reaction norms for a perennial grass under competitive stress

The developmental reaction norm (DRN) represents the set of ontogenetic trajectories that can be produced by a genotype exposed to different environmental conditions. Genetic variation in the DRN for growth traits and in the patterns of biomass allocation is critical to phenotypic evolution in heterogeneous environments. The DRN and patterns of biomass allocation were investigated in 11 clones of the caespitose, corm-forming, perennial grass Phleum pratense in relation to competitive stress imparted by Lolium perenne in a 16 week glasshouse experiment. A separate experiment assessed the ability of basal buds flanking a corm to sprout and the relationship of corm mass to sprout mass for the same clones. Corm fresh mass varied among clones and was significantly correlated with the dry mass of the tillers that sprouted from basal buds. In the competition experiment, clones in competitive environments varied significantly from those in non-competetive environments in terms of their DRNs for number of tillers and shoot dry mass. Thus, selection of DRNs would favour different genotypes in the two environments and at different times. Significant negative genetic correlations were detected for tiller number and mean tiller mass in the noncompetitive, but not the competitive, environment. Biomass allocation to stem bases was significantly greater for clones under competitive stress. Allocation to storage tissues such as corms may be adaptive if it enhances persistence in the competitive field environments typically occupied by caespitose grasses. Root and shoot allocation showed a significant clone by competition interaction. For P. pratense, genotypic variation in growth trajectories plays an important role in determining variation in individual performance, a condition necessary for the continued evolution of the DRN.     

The evolution of life histories in spatially heterogeneous environments: Optimal reaction norms revisited

Summary Natural populations live in heterogeneous environments, where habitat variation drives the evolution of phenotypic plasticity. The key feature of population structure addressed in this paper is the net flow of individuals from source (good) to sink (poor) habitats. These movements make it necessary to calculate fitness across the full range of habitats encountered by the population, rather than independently for each habitat. As a consequence, the optimal phenotype in a given habitat not only depends on conditions there but is linked to the performance of individuals in other habitats. We generalize the Euler-Lotka equation to define fitness in a spatially heterogeneous environment in which individuals disperse among habitats as newborn and then stay in a given habitat for life. In this case, maximizing fitness (the rate of increase over all habitats) is equivalent to maximizing the reproductive value of newborn in each habitat but not to maximizing the rate of increase that would result if individuals in each habitat were an isolated population. The new equation can be used to find optimal reaction norms for life history traits, and examples are calculated for age at maturity and clutch size. In contrast to previous results, the optimal reaction norm differs from the line connecting local adaptations of isolated populations each living in only one habitat. Selection pressure is higher in good and frequent habitats than in poor and rare ones. A formula for the relative importance of these two factors allows predictions of the habitat in which the genetic variance about the optimal reaction norm should be smallest.