A quantitative model of temperature-dependent diapause progression.

Winter diapause in insects is commonly terminated through cold exposure, which, like vernalization in plants, prevents development before spring arrives. Currently, quantitative understanding of the temperature dependence of diapause termination is limited, likely because diapause phenotypes are generally cryptic to human eyes. We introduce a methodology to tackle this challenge. By consecutively moving butterfly pupae of the species Pieris napi from several different cold conditions to 20 °C, we show that diapause termination proceeds as a temperature-dependent rate process, with maximal rates at relatively cold temperatures and low rates at warm and extremely cold temperatures. Further, we show that the resulting thermal reaction norm can predict P. napi diapause termination timing under variable temperatures. Last, we show that once diapause is terminated in P. napi, subsequent development follows a typical thermal performance curve, with a maximal development rate at around 31 °C and a minimum at around 2 °C. The sequence of these thermally distinct processes (diapause termination and postdiapause development) facilitates synchronous spring eclosion in nature; cold microclimates where diapause progresses quickly do not promote fast postdiapause development, allowing individuals in warmer winter microclimates to catch up, and vice versa. The unveiling of diapause termination as one temperature-dependent rate process among others promotes a parsimonious, quantitative, and predictive model, wherein winter diapause functions both as an adaptation against premature development during fall and winter and for synchrony in spring.

Methodological artefacts cause counter-intuitive evolutionary conclusions in a simulation study.

In their simulation study, Garcia-Costoya et al. (2023) conclude that evolutionary constraints might aid populations facing climate change. However, we are concerned that this conclusion is largely a consequence of the simulated temperature variation being too small, and, most importantly, that uneven limitations to standing variation disadvantage unconstrained populations.

Urbanization extends flight phenology and leads to local adaptation of seasonal plasticity in Lepidoptera.

Urbanization is gaining force globally, which challenges biodiversity, and it has recently also emerged as an agent of evolutionary change. Seasonal phenology and life cycle regulation are essential processes that urbanization is likely to alter through both the urban heat island effect (UHI) and artificial light at night (ALAN). However, how UHI and ALAN affect the evolution of seasonal adaptations has received little attention. Here, we test for the urban evolution of seasonal life-history plasticity, specifically changes in the photoperiodic induction of diapause in two lepidopterans, Pieris napi (Pieridae) and Chiasmia clathrata (Geometridae). We used long-term data from standardized monitoring and citizen science observation schemes to compare yearly phenological flight curves in six cities in Finland and Sweden to those of adjacent rural populations. This analysis showed for both species that flight seasons are longer and end later in most cities, suggesting a difference in the timing of diapause induction. Then, we used common garden experiments to test whether the evolution of the photoperiodic reaction norm for diapause could explain these phenological changes for a subset of these cities. These experiments demonstrated a genetic shift for both species in urban areas toward a lower daylength threshold for direct development, consistent with predictions based on the UHI but not ALAN. The correspondence of this genetic change to the results of our larger-scale observational analysis of in situ flight phenology indicates that it may be widespread. These findings suggest that seasonal life cycle regulation evolves in urban ectotherms and may contribute to ecoevolutionary dynamics in cities.

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 to seasonal cues at the fronts of two parallel, climate-induced butterfly range expansions.

Climate change allows species to expand polewards, but non-changing environmental features may limit expansions. Daylength is unaffected by climate and drives life cycle timing in many animals and plants. Because daylength varies over latitudes, poleward-expanding populations must adapt to new daylength conditions. We studied local adaptation to daylength in the butterfly Lasiommata megera, which is expanding northwards along several routes in Europe. Using common garden laboratory experiments with controlled daylengths, we compared diapause induction between populations from the southern-Swedish core range and recently established marginal populations from two independent expansion fronts in Sweden. Caterpillars from the northern populations entered diapause in clearly longer daylengths than those from southern populations, with the exception of caterpillars from one geographically isolated population. The northern populations have repeatedly and rapidly adapted to their local daylengths, indicating that the common use of daylength as seasonal cue need not strongly limit climate-induced insect range expansions.

The co-existence of multiple oak leaf flushes contributes to the large within-tree variation in chemistry, insect attack and pathogen infection.

Many plant species produce multiple leaf flushes during the growing season, which might have major consequences for within-plant variation in chemistry and species interactions. Yet, we lack a theoretical or empirical framework for how differences among leaf flushes might shape variation in damage by insects and diseases. We assessed the impact of leaf flush identity on leaf chemistry, insect attack and pathogen infection on the pedunculate oak Quercus robur by sampling leaves from each leaf flush in 20 populations across seven European countries during an entire growing season. The first leaf flush had higher levels of primary compounds, and lower levels of secondary compounds, than the second flush, whereas plant chemistry was highly variable in the third flush. Insect attack decreased from the first to the third flush, whereas infection by oak powdery mildew was lowest on leaves from the first flush. The relationship between plant chemistry, insect attack and pathogen infection varied strongly among leaf flushes and seasons. Our findings demonstrate the importance of considering differences among leaf flushes for our understanding of within-tree variation in chemistry, insect attack and disease levels, something particularly relevant given the expected increase in the number of leaf flushes with climate change.

Extensive transcriptomic profiling of pupal diapause in a butterfly reveals a dynamic phenotype.

Diapause is a common adaptation for overwintering in insects that is characterized by arrested development and increased tolerance to stress and cold. While the expression of specific candidate genes during diapause have been investigated, there is no general understanding of the dynamics of the transcriptional landscape as a whole during the extended diapause phenotype. Such a detailed temporal insight is important as diapause is a vital aspect of life cycle timing. Here, we performed a time-course experiment using RNA-Seq on the head and abdomen in the butterfly Pieris napi. In both body parts, comparing diapausing and nondiapausing siblings, differentially expressed genes are detected from the first day of pupal development and onwards, varying dramatically across these formative stages. During diapause there are strong gene expression dynamics present, revealing a preprogrammed transcriptional landscape that is active during the winter. Different biological processes appear to be active in the two body parts. Finally, adults emerging from either the direct or diapause pathways do not show large transcriptomic differences, suggesting the adult phenotype is strongly canalized.

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.

Spring phenology and pathogen infection affect multigenerational plant attackers throughout the growing season.

Climate change has been shown to advance spring phenology, increase the number of insect generations per year (multivoltinism) and increase pathogen infection levels. However, we lack insights into the effects of plant spring phenology and the biotic environment on the preference and performance of multivoltine herbivores and whether such effects extend into the later part of the growing season. To this aim, we used a multifactorial growth chamber experiment to examine the influence of spring phenology on plant pathogen infection, and how the independent and interactive effects of spring phenology and plant pathogen infection affect the preference and performance of multigenerational attackers (the leaf miner Tischeria ekebladella and the aphid Tuberculatus annulatus) on the pedunculate oak in the early, mid and late parts of the plant growing season. Pathogen infection was highest on late phenology plants, irrespective of whether inoculations were conducted in the early, mid or late season. The leaf miner consistently preferred to oviposit on middle and late phenology plants, as well as healthy plants, during all parts of the growing season, whereas we detected an interactive effect between spring phenology and pathogen infection on the performance of the leaf miner. Aphids preferred healthy, late phenology plants during the early season, healthy plants during the mid season, and middle phenology plants during the late season, whereas aphid performance was consistently higher on healthy plants during all parts of the growing season. Our findings highlight that the impact of spring phenology on pathogen infection and the preference and performance of insect herbivores is not restricted to the early season, but that its imprint is still present - and sometimes equally strong - during the peak and end of the growing season. Plant pathogens generally negatively affected herbivore preference and performance, and modulated the effects of spring phenology. We conclude that spring phenology and pathogen infection are two important factors shaping the preference and performance of multigenerational plant attackers, which is particularly relevant given the current advance in spring phenology, pathogen outbreaks and increase in voltinism with climate change.

Thermal performance under constant temperatures can accurately predict insect development times across naturally variable microclimates.

External conditions can drive biological rates in ectotherms by directly influencing body temperatures. While estimating the temperature dependence of performance traits such as growth and development rate is feasible under controlled laboratory settings, predictions in nature are difficult. One major challenge lies in translating performance under constant conditions to fluctuating environments. Using the butterfly Pieris napi as model system, we show that development rate, an important fitness trait, can be accurately predicted in the field using models parameterized under constant laboratory temperatures. Additionally, using a factorial design, we show that accurate predictions can be made across microhabitats but critically hinge on adequate consideration of non-linearity in reaction norms, spatial heterogeneity in microclimate and temporal variation in temperature. Our empirical results are also supported by a comparison of published and simulated data. Conclusively, our combined results suggest that, discounting direct effects of temperature, insect development rates are generally unaffected by thermal fluctuations.

Parental effects influence life history traits and covary with an environmental cline in common frog populations.

Across latitudinal clines, the juvenile developmental rates of ectotherms often covary with the length of the growing season, due to life-history trade-offs imposed by the time-constrained environments. However, as the start of the growing season often varies substantially across years, adaptive parental effects on juvenile developmental rates may mediate the costs of a delayed season. By employing a meta-analysis, we tested whether larval developmental rates across a latitudinal cline of the common frog (Rana temporaria) are affected by fluctuating onsets of breeding, across years. We predicted that larval developmental rate will be inversely related to the onset of breeding, and that northern populations will be more prone to shorten their developmental rate in response to late breeding, as the costs of delayed metamorphosis should be highest in areas with a shorter growing season. We found that the larval period of both northern and southern populations responded to parental environmental conditions to a similar degree in absolute terms, but in different directions. In northern populations, a late season start correlated with decreased development time, suggesting that the evolution of parental effects aids population persistence in time-constrained environments. In southern populations, late season start correlated with increased development time, which could potentially be explained as a predator avoidance strategy. Our findings suggest that local ecological variables can induce adaptive parental effects, but responses are complex, and likely trade-off with other ecological factors.

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.

Developmental plasticity in metabolism but not in energy reserve accumulation in a seasonally polyphenic butterfly.

The evolution of seasonal polyphenisms (discrete phenotypes in different annual generations) associated with alternative developmental pathways of diapause (overwintering) and direct development is favoured in temperate insects. Seasonal life history polyphenisms are common and include faster growth and development under direct development than in diapause. However, the physiological underpinnings of this difference remain poorly known despite its significance for understanding the evolution of polyphenisms. We measured respiration and metabolic rates through the penultimate and final larval instars in the butterfly Pieris napi and show that directly developing larvae grew and developed faster and had a higher metabolic rate than larvae entering pupal diapause. The metabolic divergence appeared only in the final instar, that is, after induction of the developmental pathway that takes place in the penultimate instar in P. napi. The accumulation of fat reserves during the final larval instar was similar under diapause and direct development, which was unexpected as diapause is predicted to select for exaggerated reserve accumulation. This suggests that overwinter survival in diapause does not require larger energy reserves than direct development, likely because of metabolic suppression in diapause pupae. The results, nevertheless, demonstrate that physiological changes coincide with the divergence of life histories between the alternative developmental pathways, thus elucidating the proximate basis of seasonal life history polyphenisms.

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.

Elucidating mechanisms for insect body size: partial support for the oxygen-dependent induction of moulting hypothesis.

Body size is a key life history trait, and knowledge of its mechanistic basis is crucial in life history biology. Such knowledge is accumulating for holometabolous insects, whose growth is characterised and body size affected by moulting. According to the oxygen-dependent induction of moulting (ODIM) hypothesis, moult is induced at a critical mass at which oxygen demand of growing tissues overrides the supply from the tracheal respiratory system, which principally grows only at moults. Support for the ODIM hypothesis is controversial, partly because of a lack of proper data to explicitly test the hypothesis. The ODIM hypothesis predicts that the critical mass is positively correlated with oxygen partial pressure (PO2 ) and negatively with temperature. To resolve the controversy that surrounds the ODIM hypothesis, we rigorously test these predictions by exposing penultimate-instar Orthosia gothica (Lepidoptera: Noctuidae) larvae to temperature and moderate PO2  manipulations in a factorial experiment. The relative mass increment in the focal instar increased along with increasing PO2 , as predicted, but there was only weak suggestive evidence of the temperature effect. Probably owing to a high measurement error in the trait, the effect of PO2  on the critical mass was sex specific; high PO2  had a positive effect only in females, whereas low PO2  had a negative effect only in males. Critical mass was independent of temperature. Support for the ODIM hypothesis is partial because of only suggestive evidence of a temperature effect on moulting, but the role of oxygen in moult induction seems unambiguous. The ODIM mechanism thus seems worth considering in body size analyses.

Metabolome dynamics of diapause in the butterfly Pieris napi: distinguishing maintenance, termination and post-diapause phases.

Diapause is a deep resting stage facilitating temporal avoidance of unfavourable environmental conditions, and is used by many insects to adapt their life cycle to seasonal variation. Although considerable work has been invested in trying to understand each of the major diapause stages (induction, maintenance and termination), we know very little about the transitions between stages, especially diapause termination. Understanding diapause termination is crucial for modelling and predicting spring emergence and winter physiology of insects, including many pest insects. In order to gain these insights, we investigated metabolome dynamics across diapause development in pupae of the butterfly Pieris napi, which exhibits adaptive latitudinal variation in the length of endogenous diapause that is uniquely well characterized. By employing a time-series experiment, we show that the whole-body metabolome is highly dynamic throughout diapause and differs between pupae kept at a diapause-terminating (low) temperature and those kept at a diapause-maintaining (high) temperature. We show major physiological transitions through diapause, separate temperature-dependent from temperature-independent processes and identify significant patterns of metabolite accumulation and degradation. Together, the data show that although the general diapause phenotype (suppressed metabolism, increased cold tolerance) is established in a temperature-independent fashion, diapause termination is temperature dependent and requires a cold signal. This revealed several metabolites that are only accumulated under diapause-terminating conditions and degraded in a temperature-unrelated fashion during diapause termination. In conclusion, our findings indicate that some metabolites, in addition to functioning as cryoprotectants, for example, are candidates for having regulatory roles as metabolic clocks or time-keepers during diapause.

Phenological synchrony between a butterfly and its host plants: Experimental test of effects of spring temperature.

Climate-driven changes in the relative phenologies of interacting species may potentially alter the outcome of species interactions. Phenotypic plasticity is expected to be important for short-term response to new climate conditions, and differences between species in plasticity are likely to influence their temporal overlap and interaction patterns. As reaction norms of interacting species may be locally adapted, any such climate-induced change in interaction patterns may vary among localities. However, consequences of spatial variation in plastic responses for species interactions are understudied. We experimentally explored how temperature affected synchrony between spring emergence of a butterfly, Anthocharis cardamines, and onset of flowering of five of its host plant species across a latitudinal gradient. We also studied potential effects on synchrony if climate-driven northward expansions would be faster in the butterflies than in host plants. Lastly, to assess how changes in synchrony influence host use we carried out an experiment to examine the importance of the developmental stage of plant reproductive structures for butterfly oviposition preference. In southern locations, the butterflies were well-synchronized with the majority of their local host plant species across temperatures, suggesting that thermal plasticity in butterfly development matches oviposition to host plant development and that thermal reaction norms of insects and plants result in similar advancement of spring phenology in response to warming. In the most northern region, however, relative phenology between the butterfly and two of its host plant species changed with increased temperature. We also show that the developmental stage of plants was important for egg-laying, and conclude that temperature-induced changes in synchrony in the northernmost region are likely to lead to shifts in host use in A. cardamines if spring temperatures become warmer. Northern expansion of butterfly populations might possibly have a positive effect on keeping up with host plant phenology with more northern host plant populations. Considering that the majority of insect herbivores exploit multiple plant species differing in their phenological response to spring temperatures, temperature-induced changes in synchrony might lead to shifts in host use and changes in species interactions in many temperate communities.

Idiosyncratic development of sensory structures in brains of diapausing butterfly pupae: implications for information processing.

Diapause is an important escape mechanism from seasonal stress in many insects. A certain minimum amount of time in diapause is generally needed in order for it to terminate. The mechanisms of time-keeping in diapause are poorly understood, but it can be hypothesized that a well-developed neural system is required. However, because neural tissue is metabolically costly to maintain, there might exist conflicting selective pressures on overall brain development during diapause, on the one hand to save energy and on the other hand to provide reliable information processing during diapause. We performed the first ever investigation of neural development during diapause and non-diapause (direct) development in pupae of the butterfly Pieris napi from a population whose diapause duration is known. The brain grew in size similarly in pupae of both pathways up to 3 days after pupation, when development in the diapause brain was arrested. While development in the brain of direct pupae continued steadily after this point, no further development occurred during diapause until temperatures increased far after diapause termination. Interestingly, sensory structures related to vision were remarkably well developed in pupae from both pathways, in contrast with neuropils related to olfaction, which only developed in direct pupae. The results suggest that a well-developed visual system might be important for normal diapause development.

Winter chilling speeds spring development of temperate butterflies.

Understanding and predicting phenology has become more important with ongoing climate change and has brought about great research efforts in the recent decades. The majority of studies examining spring phenology of insects have focussed on the effects of spring temperatures alone. Here we use citizen-collected observation data to show that winter cold duration, in addition to spring temperature, can affect the spring emergence of butterflies. Using spatial mixed models, we disentangle the effects of climate variables and reveal impacts of both spring and winter conditions for five butterfly species that overwinter as pupae across the UK, with data from 1976 to 2013 and one butterfly species in Sweden, with data from 2001 to 2013. Warmer springs lead to earlier emergence in all species and milder winters lead to statistically significant delays in three of the five investigated species. We also find that the delaying effect of winter warmth has become more pronounced in the last decade, during which time winter durations have become shorter. For one of the studied species, Anthocharis cardamines (orange tip butterfly), we also make use of parameters determined from previous experiments on pupal development to model the spring phenology. Using daily temperatures in the UK and Sweden, we show that recent variation in spring temperature corresponds to 10-15 day changes in emergence time over UK and Sweden, whereas variation in winter duration corresponds to 20 days variation in the south of the UK versus only 3 days in the south of Sweden. In summary, we show that short winters delay phenology. The effect is most prominent in areas with particularly mild winters, emphasising the importance of winter for the response of ectothermic animals to climate change. With climate change, these effects may become even stronger and apply also at higher latitudes.

Do respiratory limitations affect metabolism of insect larvae before moulting? An empirical test at the individual level.

Recent data suggest that oxygen limitation may induce moulting in larval insects. This oxygen-dependent induction of moulting (ODIM) hypothesis stems from the fact that the tracheal respiratory system of insects grows primarily at moults, whereas tissue mass increases massively between moults. This may result in a mismatch between oxygen supply and demand at the end of each larval instar because oxygen demand of growing tissues exceeds the relatively fixed supply capacity of the respiratory system. The ODIM hypothesis predicts that, within larval instars, respiration and metabolic rates of an individual larva first increase with increasing body mass but eventually level off once the supply capacity of the tracheal system starts to constrain metabolism. Here, we provide the first individual-level test of this key prediction of the ODIM hypothesis. We use a novel methodology where we repeatedly measure respiration and metabolic rates throughout the penultimate- and final-instar larvae in the butterfly Pieris napi In the penultimate instar, respiration and metabolic rates gradually decelerated along with growth, supporting the ODIM hypothesis. However, respiration and metabolic rates increased linearly during growth in the final instar, contradicting the prediction. Moreover, our data suggest considerable variation among individuals in the association between respiration rate and mass in the final instar. Overall, the results provide partial support for the ODIM hypothesis and suggest that oxygen limitation may emerge gradually within a larval instar. The results also suggest that there may be different moult induction mechanisms in larva-to-larva moults compared with the final metamorphic moult.