Testing for variation in photoperiodic plasticity in a butterfly: Inconsistent effects of circadian genes between geographic scales.

The genetic components of the circadian clock have been implicated as involved in photoperiodic regulation of winter diapause across various insect groups, thereby contributing to adaptation to adverse seasonal conditions. So far, the effects of within-population variation in these genes have not been well explored. Here, we present an experimental test of the effects of within-population variation at two circadian genes, timeless and period, on photoperiodic responses in the butterfly Pararge aegeria. While nonsynonymous candidate SNPs in both of these genes have previously shown to be associated with diapause induction on a between-population level, in the present experiment no such effect was found on a within-population level. In trying to reconcile these results, we examine sequence data, revealing considerable, previously unknown protein-level variation at both timeless and period across Scandinavian populations, including variants unique to the population studied here. Hence, we hypothesize that these variants may counteract the previously observed diapause-averting effect of the candidate SNPs, possibly explaining the difference in results between the experiments. Whatever the cause, these results highlight how the effects of candidate SNPs may sometimes vary across genetic backgrounds, which complicates evolutionary interpretations of geographic patterns of genetic variation.

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.

Parasitism rate differs between herbivore generations in the univoltine, but not bivoltine, range.

With climate change, plant-feeding insects increase their number of annual generations (voltinism). However, to what degree the emergence of a new herbivore generation affects the parasitism rate has not been explored. We performed a field experiment to test whether the parasitism rate differs between the first and the second generations of a specialist leaf miner (Tischeria ekebladella), both in the naturally univoltine and bivoltine parts of the leaf miner's distribution. We found an interactive effect between herbivore generation and geographical range on the parasitism rate. The parasitism rate was higher in the first compared to the second host generation in the part of the range that is naturally univoltine, whereas it did not differ between generations in the bivoltine range. Our experiment highlights that shifts in herbivore voltinism might release top-down control, with potential consequences for natural and applied systems.

Seasonal specialization drives divergent population dynamics in two closely related butterflies.

Seasons impose different selection pressures on organisms through contrasting environmental conditions. How such seasonal evolutionary conflict is resolved in organisms whose lives span across seasons remains underexplored. Through field experiments, laboratory work, and citizen science data analyses, we investigate this question using two closely related butterflies (Pieris rapae and P. napi). Superficially, the two butterflies appear highly ecologically similar. Yet, the citizen science data reveal that their fitness is partitioned differently across seasons. Pieris rapae have higher population growth during the summer season but lower overwintering success than do P. napi. We show that these differences correspond to the physiology and behavior of the butterflies. Pieris rapae outperform P. napi at high temperatures in several growth season traits, reflected in microclimate choice by ovipositing wild females. Instead, P. rapae have higher winter mortality than do P. napi. We conclude that the difference in population dynamics between the two butterflies is driven by seasonal specialization, manifested as strategies that maximize gains during growth seasons and minimize harm during adverse seasons, respectively.

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.

Time- and temperature-dependent dynamics of prothoracicotropic hormone and ecdysone sensitivity co-regulate pupal diapause in the green-veined white butterfly Pieris napi.

Diapause, a general shutdown of developmental pathways, is a vital adaptation allowing insects to adjust their life cycle to adverse environmental conditions such as winter. Diapause in the pupal stage is regulated by the major developmental hormones prothoracicotropic hormone (PTTH) and ecdysone. Termination of pupal diapause in the butterfly Pieris napi depends on low temperatures; therefore, we study the temperature-dependence of PTTH secretion and ecdysone sensitivity dynamics throughout diapause, with a focus on diapause termination. While PTTH is present throughout diapause in the cell bodies of two pairs of neurosecretory cells in the brain, it is absent in the axons, and the PTTH concentration in the haemolymph is significantly lower during diapause than during post diapause development, indicating that the PTTH signaling is reduced during diapause. The sensitivity of pupae to ecdysone injections is dependent on diapause stage. While pupae are sensitive to ecdysone during early diapause initiation, they gradually lose this sensitivity and become insensitive to non-lethal concentrations of ecdysone about 30 days into diapause. At low temperatures, reflecting natural overwintering conditions, diapause termination propensity after ecdysone injection is precocious compared to controls. In stark contrast, at high temperatures reflecting late summer and early autumn conditions, sensitivity to ecdysone does not return. Thus, here we show that PTTH secretion is reduced during diapause, and additionally, that the low ecdysone sensitivity of early diapause maintenance is lost during termination in a temperature dependent manner. The link between ecdysone sensitivity and low-temperature dependence reveals a putative mechanism of how diapause termination operates in insects that is in line with adaptive expectations for diapause.

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.

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.

Changes in the foliar fungal community between oak leaf flushes along a latitudinal gradient in Europe.

Aim: Leaves support a large diversity of fungi, which are known to cause plant diseases, induce plant defences or influence leaf senescence and decomposition. To advance our understanding of how foliar fungal communities are structured and assembled, we assessed to what extent leaf flush and latitude can explain the within- and among-tree variation in foliar fungal communities. Location: A latitudinal gradient spanning c. 20 degrees in latitude in Europe. Taxa: The foliar fungal community associated with a foundation tree species, the pedunculate oak Quercus robur. Methods: We examined the main and interactive effects of leaf flush and latitude on the foliar fungal community by sampling 20 populations of the pedunculate oak Quercus robur across the tree's range. We used the ITS region as a target for characterization of fungal communities using DNA metabarcoding. Results: Species composition, but not species richness, differed between leaf flushes. Across the latitudinal gradient, species richness was highest in the central part of the oak's distributional range, and foliar fungal community composition shifted along the latitudinal gradient. Among fungal guilds, the relative abundance of plant pathogens and mycoparasites was lower on the first leaf flush, and the relative abundance of plant pathogens and saprotrophs decreased with latitude. Conclusions: Changes in community composition between leaf flushes and along the latitudinal gradient were mostly a result of species turnover. Overall, our findings demonstrate that leaf flush and latitude explain 5%-22% of the small- and large-scale spatial variation in the foliar fungal community on a foundation tree within the temperate region. Using space-for-time substitution, we expect that foliar fungal community structure will change with climate warming, with an increase in the abundance of plant pathogens and mycoparasites at higher latitudes, with major consequences for plant health, species interactions and ecosystem dynamics.

Sensory Organ Investment Varies with Body Size and Sex in the Butterfly Pieris napi.

In solitary insect pollinators such as butterflies, sensory systems must be adapted for multiple tasks, including nectar foraging, mate-finding, and locating host-plants. As a result, the energetic investments between sensory organs can vary at the intraspecific level and even among sexes. To date, little is known about how these investments are distributed between sensory systems and how it varies among individuals of different sex. We performed a comprehensive allometric study on males and females of the butterfly Pieris napi where we measured the sizes and other parameters of sensory traits including eyes, antennae, proboscis, and wings. Our findings show that among all the sensory traits measured, only antenna and wing size have an allometric relationship with body size and that the energetic investment in different sensory systems varies between males and females. Moreover, males had absolutely larger antennae and eyes, indicating that they invest more energy in these organs than females of the same body size. Overall, the findings of this study reveal that the size of sensory traits in P. napi are not necessarily related to body size and raises questions about other factors that drive sensory trait investment in this species and in other insect pollinators in general.

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.

Variation in developmental rates is not linked to environmental unpredictability in annual killifishes.

Comparative evidence suggests that adaptive plasticity may evolve as a response to predictable environmental variation. However, less attention has been placed on unpredictable environmental variation, which is considered to affect evolutionary trajectories by increasing phenotypic variation (or bet hedging). Here, we examine the occurrence of bet hedging in egg developmental rates in seven species of annual killifish that originate from a gradient of variation in precipitation rates, under three treatment incubation temperatures (21, 23, and 25°C). In the wild, these species survive regular and seasonal habitat desiccation, as dormant eggs buried in the soil. At the onset of the rainy season, embryos must be sufficiently developed in order to hatch and complete their life cycle. We found substantial differences among species in both the mean and variation of egg development rates, as well as species-specific plastic responses to incubation temperature. Yet, there was no clear relationship between variation in egg development time and variation in precipitation rate (environmental predictability). The exact cause of these differences therefore remains enigmatic, possibly depending on differences in other natural environmental conditions in addition to precipitation predictability. Hence, if species-specific variances are adaptive, the relationship between development and variation in precipitation is complex and does not diverge in accordance with simple linear relationships.

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.

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.

A region of the sex chromosome associated with population differences in diapause induction contains highly divergent alleles at clock genes.

Developmental plasticity describes the capacity of individuals with the same genotype to induce permanent change in a phenotype depending on a specific external input. One well-studied example of adaptive developmental plasticity is the induction of facultative diapause in insects. Studies investigating the inheritance of diapause induction have suggested diverse genetic origins. However, only few studies have performed genome-wide scans to identify genes affecting the induction decision. Here we compare two populations of the butterfly Pieris napi that differ in the propensity to enter diapause, and despite showing a low genome-wide divergence, we identify a few genomic regions that show high divergence between populations. We then identified a single genomic region associated with diapause induction by genotyping diapausing and directly developing siblings from backcrosses of these populations. This region is located on the Z chromosome and contained three circadian clock genes, cycle, clock, and period. Additionally, period harbored the largest number of SNPs showing complete fixation between populations. We conclude that the heritable basis of between-population variation in the plasticity that determines diapause induction resides on the Z chromosome, with the period gene being the prime candidate for the genetic basis of adaptive plasticity.

Evolutionary impacts of winter climate change on insects.

Overwintering is a serious challenge for insects, and winters are rapidly changing as climate shifts. The capacity for phenotypic plasticity and evolutionary adaptation will determine which species profit or suffer from these changes. Here we discuss current knowledge on the potential and evidence for evolution in winter-relevant traits among insect species and populations. We conclude that the best evidence for evolutionary shifts in response to changing winters remain those related to changes in phenology, but all evidence points to cold hardiness as also having the potential to evolve in response to climate change. Predicting future population sizes and ranges relies on understanding to what extent evolution in winter-related traits is possible, and remains a serious challenge.

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.