VIABILITY SELECTION ON SEASONALLY POLYPHENIC TRAITS: WING MELANIN PATTERN IN WESTERN WHITE BUTTERFLIES

Wing morphological responses to latitude and colonisation in a range expanding butterfly

Populations undergoing rapid climate-driven range expansion experience distinct selection regimes dominated both by increased dispersal at the leading edges and steep environmental gradients. Characterisation of traits associated with such expansions provides insight into the selection pressures and evolutionary constraints that shape demographic and evolutionary responses. Here we investigate patterns in three components of wing morphology (size, shape, colour) often linked to dispersal ability and thermoregulation, along latitudinal gradients of range expansion in the Speckled Wood butterfly (Pararge aegeria) in Britain (two regions of expansion in England and Scotland). We measured 774 males from 54 sites spanning 799 km with a 10-year mean average temperature gradient of 4 °C. A geometric morphometric method was used to investigate variation in size and shape of forewings and hindwings; colour, pattern, and contrast of the wings were examined using a measure of lightness (inverse degree of melanism). Overall, wing size increased with latitude by ∼2% per 100 km, consistent with Bergmann’s rule. Forewings became more rounded and hindwings more elongated with history of colonisation, possibly reflecting selection for increased dispersal ability. Contrary to thermal melanism expectations, wing colour was lighter where larvae developed at cooler temperatures and unrelated to long-term temperature. Changes in wing spot pattern were also detected. High heterogeneity in variance among sites for all of the traits studied may reflect evolutionary time-lags and genetic drift due to colonisation of new habitats. Our study suggests that temperature-sensitive plastic responses for size and colour interact with selection for dispersal traits (wing size and shape). Whilst the plastic and evolutionary responses may in some cases act antagonistically, the rapid expansion of P. aegeria implies an overall reinforcing effect between these two mechanisms.

NATURAL SELECTION FOR ENVIRONMENTALLY INDUCED PHENOTYPES IN TADPOLES

Models suggest that phenotypic plasticity is maintained in situations where the optimal phenotype differs through time or space, so that selection acts in different directions in different environments. Some empirical work supports the general premise of this prediction because phenotypes induced by a particular environment sometimes perform better than other phenotypes when tested in that environment. We have extended these results by estimating the targets of selection in Pseudacris triseriata tadpoles in environments without predators and with larval Anax dragonflies. Tadpoles displayed significant behavioral and morphological plasticity when reared in the presence and absence of nonlethal dragonflies for 32 days in cattle tanks. We measured selection in the absence of free predators by regressing growth and survival in the tanks against activity and several measures of tail and body shape. We measured selection in the presence of predators by exposing groups of 10 tadpoles to Anax in overnight predation trials and regressing the average phenotype of survivors against the number of tadpoles killed. Selection in the two environments acted in opposite directions on both tail and body shape, although the affected fitness components were different. In the presence of Anax, tadpoles with shallow and narrow body, deep tail fin, and wide tail muscle survived best. In the absence of free predators, tadpoles with narrow tail muscle grew significantly faster, and those with shallow tail fin and deep body grew somewhat faster. Activity was unrelated to survival or growth in either environment. Developmental plasticity in tail shape closely paralleled selection, because tail fin depth increased after long-term exposure to Anax and tail muscle width tended to increase. In contrast, there was no plasticity in body shape in spite of strong selection for decreasing body depth. Thus, when confronted with a dragonfly predator, P. triseriata tadpoles adjusted their tail shape (but not body shape) almost exactly in the direction of selection imposed by Anax. These results suggest that phenotypic plasticity in some morphological traits, such as tail depth and tail muscle width, has evolved under intermittent selection by dragonflies. Other traits that undergo selection by dragonflies, such as body morphology, appear developmentally rigid, perhaps because of historically strong opposing selection in nature or other constraints.

MEASURING NATURAL SELECTION ON PHENOTYPIC PLASTICITY

To understand natural selection we need to integrate its measure across environments. We present a method for measuring phenotypic selection that combines the potential for both environmental variation and phenotypic plasticity. The method uses path analysis and a measure of selection that is analogous to selection on breeding values. For individuals growing in alternative environments, paths are created that represent potential changes in the environment. The probabilities for these changes are then multiplied by the path coefficients to calculate selection coefficients. Selection on plasticity is measured as the difference in selection within each environment. We illustrate these methods using data on selection in an experimental population of Arabidopsis thaliana. Individuals from 36 families were grown in one of four environments, a factorial combination of shaded/open and early/late shading. For final height of the inflorescence, there was positive selection in both the open and shaded environments and negative selection on plasticity of height. For bolting time, there was also positive selection in both environments, but no selection on plasticity. We show how to use this information to examine how selection would change with changes in environmental frequencies and their transition probabilities. These methods can be expanded to encompass continuous traits and continuous environments as well as other complexities of natural selection.

Prey from the eyes of predators: Color discriminability of aposematic and mimetic butterflies from an avian visual perspective

Predation exerts strong selection on mimetic butterfly wing color patterns, which also serve other functions such as sexual selection. Therefore, specific selection pressures may affect the sexes and signal components differentially. We tested three predictions about the evolution of mimetic resemblance by comparing wing coloration of aposematic butterflies and their Batesian mimics: (a) females gain greater mimetic advantage than males and therefore are better mimics, (b) due to intersexual genetic correlations, sexually monomorphic mimics are better mimics than female-limited mimics, and (c) mimetic resemblance is better on the dorsal wing surface that is visible to predators in flight. Using a physiological model of avian color vision, we quantified mimetic resemblance from predators' perspective, which showed that female butterflies were better mimics than males. Mimetic resemblance in female-limited mimics was comparable to that in sexually monomorphic mimics, suggesting that intersexual genetic correlations did not constrain adaptive response to selection for female-limited mimicry. Mimetic resemblance on the ventral wing surface was better than that on the dorsal wing surface, implying stronger natural and sexual selection on ventral and dorsal surfaces, respectively. These results suggest that mimetic resemblance in butterfly mimicry rings has evolved under various selective pressures acting in a sex- and wing surface-specific manner.

Predator-induced phenotypic plasticity within- and across-generations: a challenge for theory?

Much work has shown that the environment can induce non-genetic changes in phenotype that span multiple generations. Theory predicts that predictable environmental variation selects for both increased within- and across-generation responses. Yet, to the best of our knowledge, there are no empirical tests of this prediction. We explored the relationship between within- versus across-generation plasticity by evaluating the influence of predator cues on the life-history traits of Daphnia ambigua. We measured the duration of predator-induced transgenerational effects, determined when transgenerational responses are induced, and quantified the cues that activate transgenerational plasticity. We show that predator exposure during embryonic development causes earlier maturation and increased reproductive output. Such effects are detectable two generations removed from predator exposure and are similar in magnitude in response to exposure to cues emitted by injured conspecifics. Moreover, all experimental contexts and traits yielded a negative correlation between within- versus across-generation responses. That is, responses to predator cues within- and across-generations were opposite in sign and magnitude. Although many models address transgenerational plasticity, none of them explain this apparent negative relationship between within- and across-generation plasticities. Our results highlight the need to refine the theory of transgenerational plasticity.

Accelerated evolution of morph-biased genes in pea aphids

Phenotypic plasticity, the production of alternative phenotypes (or morphs) from the same genotype due to environmental factors, results in some genes being expressed in a morph-biased manner. Theoretically, these morph-biased genes experience relaxed selection, the consequence of which is the buildup of slightly deleterious mutations at these genes. Over time, this is expected to result in increased protein divergence at these genes between species and a signature of relaxed purifying selection within species. Here we test these theoretical expectations using morph-biased genes in the pea aphid, a species that produces multiple morphs via polyphenism. We find that morph-biased genes exhibit faster rates of evolution (in terms of dN/dS) relative to unbiased genes and that divergence generally increases with increasing morph bias. Further, genes with expression biased toward rarer morphs (sexual females and males) show faster rates of evolution than genes expressed in the more common morph (asexual females), demonstrating that the amount of time a gene spends being expressed in a morph is associated with its rate of evolution. And finally, we show that genes expressed in the rarer morphs experience decreased purifying selection relative to unbiased genes, suggesting that it is a relaxation of purifying selection that contributes to their faster rates of evolution. Our results provide an important empirical look at the impact of phenotypic plasticity on gene evolution.

Spalt expression and the development of melanic color patterns in pierid butterflies

Background

Little is currently known about wing pattern development in the butterfly family Pieridae, which consists mostly of black melanized elements on white or yellow/orange backgrounds. A single transcription factor, Spalt (Sal), has been previously associated with the development of some pattern elements in Pieris rapae, but it is unclear to what extent Sal is associated with patterns in other pierid species.

Results

We use immunohistochemistry targeting Sal proteins across several pierids and show that Sal is associated with dense patches of melanization across species but is not associated with vein-melanization or diffuse melanization on the wing. In addition, Sal is expressed along cross-veins and wing compartment midlines that do not develop melanization. Male and female P. rapae spots are sexually dimorphic in size and this dimorphism is also present in the domains of Sal expression. Finally, by disrupting cells positioned in the center of the anterior black spots of P. rapae, before and during the time of Sal expression, spot size was reduced.

Conclusions

Our results suggest, but do not conclusively show, that pierid spots may develop in a manner similar to that of nymphalid eyespots, that is, containing a group of signaling cells at the center of the pattern responsible for the differentiation of the complete spot, and that spots and eyespots share at least one signal-response gene in common, the transcription factor Sal. We propose that focal differentiation and focal signaling mechanisms evolved prior to the split of the nymphalid and pierid lineages.

Phenotypic plasticity and integration in the mangrove rivulus (Kryptolebias marmoratus): a prospectus

The mangrove rivulus (Kryptolebias marmoratus) is a small fish native to mangrove ecosystems in Florida, the Caribbean, Central America, and South America. This species is one of only two self-fertilizing, hermaphroditic vertebrates capable of producing offspring that are genetically identical to both the parent and all siblings. Long bouts of selfing result in individuals with completely homozygous genotypes, effectively allowing for the production of "clones." Rivulus is also extremely sensitive to environmental change, both during development and adulthood. Life-history traits, behavior, physiology, morphology, and even sexual phenotype are shaped to a large extent by the interaction of genes with the environment, and many of these traits appear to co-vary. True reaction norms can be generated for this species in much the same way as has been done for clonally reproducing invertebrates and plants that have contributed immensely to our understanding of the evolution of phenotypic plasticity. That is, rivulus provides the opportunity to place individuals with identical genotypes in many different environments at any point during ontogeny or adulthood. In addition, rivulus populations are characterized by high genotypic diversity, a luxury not afforded by many clonal vertebrates, which allows us to evaluate variation among genotypes in the shape of reaction norms and in patterns of covariance among traits. We provide background information on phenotypic plasticity and phenotypic integration, coupled with a description of characteristics that we feel qualify rivulus as a potentially powerful model in which to study the evolution of reaction norms and covariance among traits.

The redder the better: wing color predicts flight performance in monarch butterflies

The distinctive orange and black wings of monarchs (Danaus plexippus) have long been known to advertise their bitter taste and toxicity to potential predators. Recent work also showed that both the orange and black coloration of this species can vary in response to individual-level and environmental factors. Here we examine the relationship between wing color and flight performance in captive-reared monarchs using a tethered flight mill apparatus to quantify butterfly flight speed, duration and distance. In three different experiments (totaling 121 individuals) we used image analysis to measure body size and four wing traits among newly-emerged butterflies prior to flight trials: wing area, aspect ratio (length/width), melanism, and orange hue. Results showed that monarchs with darker orange (approaching red) wings flew longer distances than those with lighter orange wings in analyses that controlled for sex and other morphometric traits. This finding is consistent with past work showing that among wild monarchs, those sampled during the fall migration are darker in hue (redder) than non-migratory monarchs. Together, these results suggest that pigment deposition onto wing scales during metamorphosis could be linked with traits that influence flight, such as thorax muscle size, energy storage or metabolism. Our results reinforce an association between wing color and flight performance in insects that is suggested by past studies of wing melansim and seasonal polyphenism, and provide an important starting point for work focused on mechanistic links between insect movement and color.

The fitness costs of developmental canalization and plasticity

Organisms are capable of an astonishing repertoire of phenotypic responses to the environment, and these often define important adaptive solutions to heterogeneous and unpredictable conditions. The terms 'phenotypic plasticity' and 'canalization' indicate whether environmental variation has a large or small effect on the phenotype. The evolution of canalization and plasticity is influenced by optimizing selection-targeting traits within environments, but inherent fitness costs of plasticity may also be important. We present a meta-analysis of 27 studies (of 16 species of plant and 7 animals) that have measured selection on the degree of plasticity independent of the characters expressed within environments. Costs of plasticity and canalization were equally frequent and usually mild; large costs were observed only in studies with low sample size. We tested the importance of several covariates, but only the degree of environmental stress was marginally positively related to the cost of plasticity. These findings suggest that costs of plasticity are often weak, and may influence phenotypic evolution only under stressful conditions.

The strength of phenotypic selection in natural populations

How strong is phenotypic selection on quantitative traits in the wild? We reviewed the literature from 1984 through 1997 for studies that estimated the strength of linear and quadratic selection in terms of standardized selection gradients or differentials on natural variation in quantitative traits for field populations. We tabulated 63 published studies of 62 species that reported over 2,500 estimates of linear or quadratic selection. More than 80% of the estimates were for morphological traits; there is very little data for behavioral or physiological traits. Most published selection studies were unreplicated and had sample sizes below 135 individuals, resulting in low statistical power to detect selection of the magnitude typically reported for natural populations. The absolute values of linear selection gradients |beta| were exponentially distributed with an overall median of 0.16, suggesting that strong directional selection was uncommon. The values of |beta| for selection on morphological and on life-history/phenological traits were significantly different: on average, selection on morphology was stronger than selection on phenology/life history. Similarly, the values of |beta| for selection via aspects of survival, fecundity, and mating success were significantly different: on average, selection on mating success was stronger than on survival. Comparisons of estimated linear selection gradients and differentials suggest that indirect components of phenotypic selection were usually modest relative to direct components. The absolute values of quadratic selection gradients |gamma| were exponentially distributed with an overall median of only 0.10, suggesting that quadratic selection is typically quite weak. The distribution of gamma values was symmetric about 0, providing no evidence that stabilizing selection is stronger or more common than disruptive selection in nature.

A simulation study of mutations in the genetic regulatory hierarchy for butterfly eyespot focus determination

The color patterns on the wings of butterflies have been an important model system in evolutionary developmental biology. A recent computational model tested genetic regulatory hierarchies hypothesized to underlie the formation of butterfly eyespot foci [Evans, T.M., Marcus, J.M., 2006. A simulation study of the genetic regulatory hierarchy for butterfly eyespot focus determination. Evol. Dev. 8, 273-283]. The computational model demonstrated that one proposed hierarchy was incapable of reproducing the known patterns of gene expression associated with eyespot focus determination in wild-type butterflies, but that two slightly modified alternative hierarchies were capable of reproducing all of the known gene expressions patterns. Here we extend the computational models previously implemented in Delphi 2.0 to two mutants derived from the squinting bush brown butterfly (Bicyclus anynana). These two mutants, comet and Cyclops, have aberrantly shaped eyespot foci that are produced by different mechanisms. The comet mutation appears to produce a modified interaction between the wing margin and the eyespot focus that results in a series of comet-shaped eyespot foci. The Cyclops mutation causes the failure of wing vein formation between two adjacent wing-cells and the fusion of two adjacent eyespot foci to form a single large elongated focus in their place. The computational approach to modeling pattern formation in these mutants allows us to make predictions about patterns of gene expression, which are largely unstudied in butterfly mutants. It also suggests a critical experiment that will allow us to distinguish between two hypothesized genetic regulatory hierarchies that may underlie all butterfly eyespot foci.

The rules of engagement: how to defend against combinations of predators

Studies of inducible defenses have traditionally examined prey responses to one predator at a time. However, prey in nature encounter combinations of predators that should force them to produce phenotypic compromises. We examined how snails (Helisoma trivolvis) alter their phenotype in the presence of three different predator species that were presented alone and in pairwise combinations. When snails were exposed to each predator alone, they formed predator-specific defenses that reflected the differences in each predator's foraging mode. When snails were exposed to pairwise combinations of predators, their phenotype was dependent on their ability to detect each predator, the risk posed by each predator, and the effectiveness of a given defense against each predator. Consequently, responses to combined predators were typically biased towards one of the predators in the pair. This suggests that prey facing combined predators do not form simple intermediate defenses and, as a result, may experience enhanced mortality risk when they encounter natural predator regimes.

Local and global costs of adaptive plasticity to density in Arabidopsis thaliana

Although phenotypic plasticity is demonstrably adaptive in a range of settings, organisms are not perfectly plastic. Costs of plasticity comprise one factor predicted to counter the evolution of this adaptive strategy, yet evidence of costs is rare. Here, we test the fitness effects of plastic life-history and morphological responses to density and costs of this plasticity in recombinant inbred lines of Arabidopsis thaliana. Several costs of plasticity and homeostasis were detected. Of particular relevance, there was a significant cost of plasticity to active stem-elongation responses, an adaptive trait in many species. There was also a cost of plasticity to apical branch production at both high and low density, which resulted from the greater suppression of basal branching in genotypes with plastic apical branch production relative to genotypes with fixed apical branch production. The presence of a cost in multiple environments (i.e., a global cost) is predicted to counter the evolution of plasticity. Experimental segregating progenies such as the one used here are expected to have higher genetic costs of plasticity than arrays of genotypes sampled from natural populations because selection should remove genotypes with costs resulting from linkage disequilibrium or epistasis. The use of experimental progeny arrays therefore increases the ability to evaluate genetic costs.

Frequency and Microenvironmental Pattern of Selection on Plastic Shade‐Avoidance Traits in a Natural Population ofImpatiens capensis

The frequency and predictability of different selective environments are important parameters in models for the evolution of plasticity but have rarely been measured empirically in natural populations. We used an experimental phytometer approach to examine the frequency, predictability, and environmental determinants of heterogeneous selection on phytochrome-mediated shade-avoidance responses in a natural population of the annual plant Impatiens capensis. The strength and direction of selection on shade-avoidance traits varied substantially on a fine spatial scale. The shade-avoidance phenotype had high relative fecundity in some microsites but was disadvantageous in other microsites. Local seedling density proved to be a surprisingly poor predictor of microenvironmental variation in the strength and direction of selection on stem elongation in this study population. At least some of this unpredictability resulted from microenvironmental variation in water availability; the shade-avoidance phenotype was more costly in dry microsites. Thus, environmental heterogeneity in resource availability can affect the relative costs and benefits of expressing shade-avoidance traits independent of local seedling density, the inductive environmental cue. Theory predicts that these conditions may promote local genetic differentiation in reaction norms in structured populations, as observed in I. capensis.

Correlations between scale structure and pigmentation in butterfly wings

We examined the correlation between color and structure of wing scales in the nymphalid butterflies Bicyclus anynana and Heliconius melpomene. All scales in B. anynana are rather similar in comparison to the clear structural differences of differently pigmented scales in H. melpomene. Where scale structural differences in H. melpomene are qualitative, they seem to be quantitative in B. anynana. There is a "gradient" in the density of some structural elements, the cross ribs, in the scales of B. anynana: black, gold, and brown scales show progressively lower cross rib density within an individual. There is, however, high individual variation in the absolute cross rib densities (i.e., scales with a particular color and cross rib density in one individual may have a different color but similar density in another individual). By ectopically inducing color pattern during early pupal development, we examined whether a scale's color and its microstructure could be uncoupled. The effect of these manipulations appears to be different in B. anynana and H. melpomene. In Bicyclus, "black" scales induced by wing damage at an ectopic location normally containing brown scales acquire both an intermediate structure and color between that of brown and normal black scales. In Heliconius, however, intermediate colors or scale structure were never observed, and scales with an altered color (due to damage) always have the same structure as normal scales with that color. The results are discussed on the basis of gene expression patterns, variability in rates of scale development and pigment, and scale sclerotization pathways.

Plasticity to light cues and resources in Arabidopsis thaliana: testing for adaptive value and costs

Plants shaded by neighbors or overhead foliage experience both a reduction in the ratio of red to far red light (R:FR), a specific cue perceived by phytochrome, and reduced photosynthetically active radiation (PAR), an essential resource. We tested the adaptive value of plasticity to crowding and to the cue and resource components of foliage shade in the annual plant Arabidopsis thaliana by exposing 36 inbred families from four natural populations to four experimental treatments: (1) high density, full sun; (2) low density, full sun; (3) low density, neutral shade; and (4) low density, low R:FR-simulated foliage shade. Genotypic selection analysis within each treatment revealed strong environmental differences in selection on plastic life-history traits. We used specific contrasts to measure plasticity to density and foliage shade, to partition responses to foliage shade into phytochrome-mediated responses to the R:FR cue and responses to PAR, and to test whether plasticity was adaptive (i.e., in the same direction as selection in each environment). Contrary to expectation, we found no evidence for adaptive plasticity to density. However, we observed both adaptive and maladaptive responses to foliage shade. In general, phytochrome-mediated plasticity to the R:FR cue of foliage shade was adaptive and counteracted maladaptive growth responses to reduced PAR. These results support the prediction that active developmental responses to environmental cues are more likely to be adaptive than are passive resource-mediated responses. Multiple regression analysis detected a few costs of adaptive plasticity and adaptive homeostasis, but such costs were infrequent and their expression depended on the environment. Thus, costs of plasticity may occasionally constrain the evolution of adaptive responses to foliage shade in Arabidopsis, but this constraint may differ among environments and is far from ubiquitous.

Climate extremes: observations, modeling, and impacts

One of the major concerns with a potential change in climate is that an increase in extreme events will occur. Results of observational studies suggest that in many areas that have been analyzed, changes in total precipitation are amplified at the tails, and changes in some temperature extremes have been observed. Model output has been analyzed that shows changes in extreme events for future climates, such as increases in extreme high temperatures, decreases in extreme low temperatures, and increases in intense precipitation events. In addition, the societal infrastructure is becoming more sensitive to weather and climate extremes, which would be exacerbated by climate change. In wild plants and animals, climate-induced extinctions, distributional and phenological changes, and species' range shifts are being documented at an increasing rate. Several apparently gradual biological changes are linked to responses to extreme weather and climate events.

THE FITNESS OF MANIPULATING PHENOTYPES: IMPLICATIONS FOR STUDIES OF FLUCTUATING ASYMMETRY AND MULTIVARIATE SELECTION

Phenotypic manipulation (or phenotypic engineering) that alters trait distributions provides a way to increase the statistical power of detecting relationships between traits and fitness. Manipulations relying on plastic responses, however, assume a specific relationship between the perturbation and the alteration of the traits when multiple traits are involved. We measured several traits, including condition measured as fluctuating asymmetry, in the ladybird beetle Harmonia axyridis under six different diets to examine how altered environments affected multiple traits and their distributions. Although diet affected fluctuating asymmetry, we found no consistent relationship between degree of asymmetry and other phenotypic measures. As expected, individual traits were altered by our treatments. Contrary to expectation, relationships among traits were not constant among diets. Our results suggest that assumptions about the relationship between condition and trait values, especially fluctuating asymmetry, cannot be made. Further, studies that use manipulated phenotypes to statistically determine the form of selection must first demonstrate that the pattern of the phenotypic correlation matrix is not itself altered by the manipulation. If the phenotypic correlation matrix is not constant, then experimental estimates of selection coefficients may not reflect selection that occurs in the wild.

Development, plasticity and evolution of butterfly eyespot patterns

The developmental and genetic bases for the formation, plasticity and diversity of eyespot patterns in butterflies are examined. Eyespot pattern mutants, regulatory gene expression, and transplants of the eyespot developmental organizer demonstrate that eyespot position, number, size and colour are determined progressively in a developmental pathway largely uncoupled from those regulating other wing-pattern elements and body structures. Species comparisons and selection experiments suggest that the evolution of eyespot patterns can occur rapidly through modulation of different stages of this pathway, and requires only single, or very few, changes in regulatory genes.