Life-cycle regulation and life history plasticity in the speckled wood butterfly: are reaction norms predictable?

Evolution of butterfly seasonal plasticity driven by climate change varies across life stages

Photoperiod is a common cue for seasonal plasticity and phenology, but climate change can create cue-environment mismatches for organisms that rely on it. Evolution could potentially correct these mismatches, but phenology often depends on multiple plastic decisions made during different life stages and seasons that may evolve separately. For example, Pararge aegeria (Speckled wood butterfly) has photoperiod-cued seasonal life history plasticity in two different life stages: larval development time and pupal diapause. We tested for climate change-associated evolution of this plasticity by replicating common garden experiments conducted on two Swedish populations 30 years ago. We found evidence for evolutionary change in the contemporary larval reaction norm-although these changes differed between populations-but no evidence for evolution of the pupal reaction norm. This variation in evolution across life stages demonstrates the need to consider how climate change affects the whole life cycle to understand its impacts on phenology.

Local adaptation of life cycles in a butterfly is associated with variation in several circadian clock genes

Many insects exhibit geographical variation in voltinism, the number of generations produced per year. This includes high-latitude species in previously glaciated areas, meaning that divergent selection on life cycle traits has taken place during or shortly after recent colonization. Here, we use a population genomics approach to compare a set of nine Scandinavian populations of the butterfly Pararge aegeria that differ in life cycle traits (diapause thresholds and voltinism) along both north-south and east-west clines. Using a de novo-assembled genome, we reconstruct colonization histories and demographic relationships. Based on the inferred population structure, we then scan the genome for candidate loci showing signs of divergent selection potentially associated with population differences in life cycle traits. The identified candidate genes include a number of components of the insect circadian clock (timeless, timeless2, period, cryptochrome and clockwork orange). Most notably, the gene timeless, which has previously been experimentally linked to life cycle regulation in P. aegeria, is here found to contain a novel 97-amino acid deletion unique to, and fixed in, a single population. These results add to a growing body of research framing circadian gene variation as a potential mechanism for generating local adaptation of life cycles.

The genome sequence of the speckled wood butterfly, Pararge aegeria (Linnaeus, 1758)

We present a genome assembly from an individual female Pararge aegeria (the speckled wood butterfly; Arthropoda; Insecta; Lepidoptera; Nymphalidae). The genome sequence is 517 megabases in span. The majority of the assembly (99.68%) is scaffolded into 29 chromosomal pseudomolecules, with the W and Z sex chromosome assembled. Gene annotation of this assembly on Ensembl has identified 12,288 protein coding genes.

Watching the days go by: Asymmetric regulation of caterpillar development by changes in photoperiod

Many insects possess the plastic ability to either develop directly to adulthood, or enter diapause and postpone reproduction until the next year, depending on environmental cues (primarily photoperiod) that signal the amount of time remaining until the end of the growth season. These two alternative pathways often differ in co-adapted life-history traits, for example, with slower development and larger size in individuals headed for diapause. The developmental timing of these differences may be of adaptive importance: If traits diverge early, the potential for phenotypic differences between the pathways is greater, whereas if traits diverge late, the risk may be lower of expressing a maladaptive phenotype if the selective environment changes during development. Here, we explore the effects of changes in photoperiodic information during life on pupal diapause and associated life-history traits in the butterfly Pararge aegeria. We find that both pupal diapause and larval development rate are asymmetrically regulated: While exposure to long days late in life (regardless of earlier experiences) was sufficient to produce nondiapause development and accelerate larval development accordingly, more prolonged exposure to short days was required to induce diapause and slow down prediapause larval development. While the two developmental pathways diverged early in development, development rates could be partially reversed by altered environmental cues. Meanwhile, pathway differences in body size were more inflexible, despite emerging late in development. These results show how several traits may be shaped by the same environmental cue (photoperiod), but along subtly different ontogenies, into an integrated phenotype.

Local adaptation of photoperiodic plasticity maintains life cycle variation within latitudes in a butterfly

The seasonal cycle varies geographically and organisms are under selection to express life cycles that optimally exploit their spatiotemporal habitats. In insects, this often means producing an annual number of generations (voltinism) appropriate to the local season length. Variation in voltinism may arise from variation in environmental factors (e.g., temperature or photoperiod) acting on a single reaction norm shared across populations, but it may also result from local adaptation of reaction norms. However, such local adaptation is poorly explored at short geographic distances, especially within latitudes. Using a combination of common-garden rearing and life cycle modeling, we have investigated the causal factors behind voltinism variation in Swedish populations of the butterfly Pararge aegeria, focusing on a set of populations that lie within a single degree of latitude but nonetheless differ in season length and voltinism. Despite considerable differences in ambient temperature between populations, modeling suggested that the key determinant of local voltinism was in fact interpopulation differences in photoperiodic response. These include differences in the induction thresholds for winter diapause, as well as differences in photoperiodic regulation of larval development, a widespread but poorly studied phenomenon. Our results demonstrate previously neglected ways that photoperiodism may mediate insect phenological responses to temperature, and emphasize the importance of local adaptation in shaping phenological patterns in general, as well as for predicting the responses of populations to changes in climate.

Genetic variation underlying local adaptation of diapause induction along a cline in a butterfly

Diapause is a life history strategy allowing individuals to arrest development until favourable conditions return, and it is commonly induced by shortened day length that is latitude specific for local populations. Although understanding the evolutionary dynamics of a threshold trait like diapause induction provides insights into the adaptive process and adaptive potential of populations, the genetic mechanism of variation in photoperiodic induction of diapause is not well understood. Here, we investigate genetic variation underlying latitudinal variation in diapause induction and the selection dynamics acting upon it. Using a genomewide scan for divergent regions between two populations of the butterfly Pararge aegeria that differ strongly in their induction thresholds, we identified and investigated the patterns of variation in those regions. We then tested the association of these regions with diapause induction using between-population crosses, finding significant SNP associations in four genes present in two chromosomal regions, one with the gene period, and the other with the genes kinesin, carnitine O-acetyltransferase and timeless. Patterns of allele frequencies in these two regions in population samples along a latitudinal cline suggest strong selection against heterozygotes at two genes within these loci (period, timeless). Evidence for additional loci modifying the diapause decision was found in patterns of allelic change in relation to induction thresholds over the cline, as well as in backcross analyses. Taken together, population-specific adaptations of diapause induction appear to be due to a combination of alleles of larger and smaller effect size, consistent with an exponential distribution of effect sizes involved in local adaption.

A comparison of the life-history traits between diapause and direct development individuals in the cotton bollworm, Helicoverpa armigera

In order to understand the differences of life-history traits between diapause and direct development individuals in the cotton bollworm, Helicoverpa armigera (Hübner) (Lepidoptera: Noctuidae), the development time, body size, growth rate, and adult longevity were investigated between the two populations, which were induced under 12:12 L:D and 16:8 L:D photoperiods, respectively, at 20, 22, and 25°C. The results indicated that the larval development time, pupal weight, adult weight, and growth rate were significantly different between diapause and direct developing individuals. The diapause developing individuals had a significantly higher pupal and adult weight and a longer larval time compared with direct developing individuals. However, the growth rate in diapause developing individuals was lower than that in the direct developing individuals. Analysis by GLM showed that larval time, pupal and adult weight, and growth rate were significantly influenced by both temperature and developmental pathway. The pupal and adult weights were greater in males than females in both developmental pathways, exhibiting sexual size dimorphism. The dimorphism in adult weight was more pronounced than in pupal weight because female pupae lost more weight at metamorphosis compared to male pupae. Protogyny was observed in both developmental pathways. However, the protogyny phenomenon was more pronounced at lower temperatures in direct developing individuals, whereas it was more pronounced in diapause developing individuals when they experienced higher temperatures in their larval stage and partial pupal period. The adult longevity of diapause developing individuals was significantly longer than that of direct developing individuals. The results reveal that the life-history strategy was different between diapause and direct developing individuals.

The evolution of alternative developmental pathways: footprints of selection on life-history traits in a butterfly

Developmental pathways may evolve to optimize alternative phenotypes across environments. However, the maintenance of such adaptive plasticity under relaxed selection has received little study. We compare the expression of life-history traits across two developmental pathways in two populations of the butterfly Pararge aegeria where both populations express a diapause pathway but one never expresses direct development in nature. In the population with ongoing selection on both pathways, the difference between pathways in development time and growth rate was larger, whereas the difference in body size was smaller compared with the population experiencing relaxed selection on one pathway. This indicates that relaxed selection on the direct pathway has allowed life-history traits to drift towards values associated with lower fitness when following this pathway. Relaxed selection on direct development was also associated with a higher degree of genetic variation for protandry expressed as within-family sexual dimorphism in growth rate. Genetic correlations for larval growth rate across sexes and pathways were generally positive, with the notable exception of correlation estimates that involved directly developing males of the population that experienced relaxed selection on this pathway. We conclude that relaxed selection on one developmental pathway appears to have partly disrupted the developmental regulation of life-history trait expression. This in turn suggests that ongoing selection may be responsible for maintaining adaptive developmental regulation along alternative developmental pathways in these populations.

Evolution of aquatic insect behaviours across a gradient of disturbance predictability

Natural disturbance regimes--cycles of fire, flood, drought or other events--range from highly predictable (disturbances occur regularly in time or in concert with a proximate cue) to highly unpredictable. While theory predicts how populations should evolve under different degrees of disturbance predictability, there is little empirical evidence of how this occurs in nature. Here, we demonstrate local adaptation in populations of an aquatic insect occupying sites along a natural gradient of disturbance predictability, where predictability was defined as the ability of a proximate cue (rainfall) to signal a disturbance (flash flood). In controlled behavioural experiments, populations from predictable environments responded to rainfall events by quickly exiting the water and moving sufficiently far from the stream to escape flash floods. By contrast, populations from less predictable environments had longer response times and lower response rates, reflecting the uncertainty inherent to these environments. Analysis with signal detection theory showed that for 13 out of 15 populations, observed response times were an optimal compromise between the competing risks of abandoning versus remaining in the stream, mediated by the rainfall-flood correlation of the local environment. Our study provides the first demonstration that populations can evolve in response to differences in disturbance predictability, and provides evidence that populations can adapt to among-stream differences in flow regime.

Phenotypic and genetic variations and correlations in multitrait developmental instability: a multivariate Bayesian model applied to Speckled Wood butterfly (Pararge aegeria) wing measurements

Phenotypic and genetic variation and covariation in developmental instability (DI) have been the subject of many debates. In this paper we develop and apply a statistical model in a Bayesian context to analyse different traits simultaneously in a multivariate model of DI. We apply the model to measurements of yellow spots on the front wing of the Speckled Wood butterfly (Pararge aegeriaL.) in a full-sib breeding experiment. We estimated the posterior distribution of the broad-sense heritability of DI averaged across the five yellow spots, which had a median of 0·19 and a 95% credibility interval ranging between 0·04 and 0·64. Phenotypic and genetic correlations in DI could not be estimated accurately with the present sample size. Yellow spots 4 and 5 appeared to show some degree of developmental integration. The importance of this model and its possible extensions are discussed.

The cost of melanization: butterfly wing coloration under environmental stress

Evolutionary studies typically focus on adaptations to particular environmental conditions, thereby often ignoring the role of possible constraints. Here we focus on the case of variation in dorsal wing melanization in a satyrine butterfly Pararge aegeria. Because melanin is a complex polymer, its synthesis may be constrained if ambient conditions limit the resource budget. This hypothesis was tested by comparing melanization among butterflies that fed as larvae on host grasses experiencing different drought-stress treatments. Treatment differences were validated both at the level of the host plant (nitrogen, carbonate, and water content) and of the butterfly (life-history traits: survival, development time, and size at maturity). Melanization rate was measured as average gray value of the basal dorsal wing area. This area, close to the thorax, is known to be functionally significant for basking in order to thermoregulate. Individuals reared on drought-stressed host plants developed paler wings, and development of darker individuals was slower and less stable as estimated by their level of fluctuating asymmetry. These results provide evidence that melanin is indeed costly to synthesize, and that differences in environmental quality can induce phenotypic variation in wing melanization. Therefore, studies dealing with spatial and/or temporal patterns of variation in wing melanization should not focus on adaptive explanations alone, but rather on a cost-benefit balance under particular sets of environmental conditions.