Environmental drivers of varying selective optima in a small passerine: A multivariate, multiepisodic approach

A Long-Lived Alpine Perennial Advances Flowering under Warmer Conditions but Not Enough to Maintain Reproductive Success

AbstractAssessing whether phenological shifts in response to climate change confer a fitness advantage requires investigating the relationships among phenology, fitness, and environmental drivers of selection. Despite widely documented advancements in phenology with warming climate, we lack empirical estimates of how selection on phenology varies in response to continuous climate drivers or how phenological shifts in response to warming conditions affect fitness. We leverage an unusual long-term dataset with repeated, individual measurements of phenology and reproduction in a long-lived alpine plant. We analyze phenotypic plasticity in flowering phenology in relation to two climate drivers, snowmelt timing and growing degree days (GDDs). Plants flower earlier with increased GDDs and earlier snowmelt, and directional selection also favors earlier flowering under these conditions. However, reproduction still declines with warming and early snowmelt, even when flowering is early. Furthermore, the steepness of this reproductive decline increases dramatically with warming conditions, resulting in very little fruit production regardless of flowering time once GDDs exceed approximately 225 degree days or snowmelt occurs before May 15. Even though advancing phenology confers a fitness advantage relative to stasis, these shifts are insufficient to maintain reproduction under warming, highlighting limits to the potential benefits of phenological plasticity under climate change.

The magnitude of selection on growth varies among years and increases under warming conditions in a subarctic seabird

Because of ongoing rapid climate change, many ecosystems are becoming both warmer and more variable, and these changes are likely to alter the magnitude and variability of natural selection acting on wild populations. Critically, changes and fluctuations in selection can impact both population demography and evolutionary change. Therefore, predicting the impacts of climate change depends on understanding the magnitude and variation in selection on traits across different life stages and environments. Long-term experiments in wild settings are a great opportunity to determine the impact of environmental conditions on selection. Here we examined variability in the strength of selection on size traits of nestling black-legged kittiwakes (Rissa tridactyla) in a 25-year study including a food supplementation experiment on Middleton Island in the Gulf of Alaska. Using mixed effect models, we examined the annual variability of stage-specific and resource-specific selection gradients across 25 years. We found that (a) larger and heavier hatchlings were the most likely to survive during early ontogeny, (b) non-food supplemented younger nestlings in a brood experienced the strongest selection, and (c) warmer conditions increased the magnitude of selection on nestling mass and affected non-food supplemented and second-hatched nestlings the most. Our results suggested that variable resource dynamics likely caused some of the changes in selection from year to year and that warming conditions increased the strength of selection on subarctic seabird growth. However, our experimental manipulation revealed that local environmental heterogeneity could buffer the selection expected from broader climatic changes. Consequently, understanding the interactive effects of local conditions and general changes in climate seems likely to improve our ability to predict future selection gradients.

Evolutionary rescue from climate change: male indirect genetic effects on lay-dates and their consequences for population persistence

Changes in avian breeding phenology are among the most apparent responses to climate change in free-ranging populations. A key question is whether populations will be able to keep up with the expected rates of environmental change. There is a large body of research on the mechanisms by which avian lay-dates track temperature change and the consequences of (mal)adaptation on population persistence. Often overlooked is the role of males, which can influence the lay-date of their mate through their effect on the prelaying environment. We explore how social plasticity causing male indirect genetic effects can help or hinder population persistence when female genes underpinning lay-date and male genes influencing female's timing of reproduction both respond to climate-mediated selection. We extend quantitative genetic moving optimum models to predict the consequences of social plasticity on the maximum sustainable rate of temperature change, and evaluate our model using a combination of simulated data and empirical estimates from the literature. Our results suggest that predictions for population persistence may be biased if indirect genetic effects and cross-sex genetic correlations are not considered and that the extent of this bias depends on sex differences in how environmental change affects the optimal timing of reproduction. Our model highlights that more empirical work is needed to understand sex-specific effects of environmental change on phenology and the fitness consequences for population dynamics. While we discuss our results exclusively in the context of avian breeding phenology, the approach we take here can be generalized to many different contexts and types of social interaction.

Variation in shape and consistency of selection between populations of the threatened Hihi (Notiomystis cincta)

The shape and intensity of natural selection can vary between years, potentially resulting in a chronic reduction of fitness as individuals need to track a continually changing optimum of fitness (i.e., a "lag load"). In endangered species, often characterized by small population size, the lack of genetic diversity is expected to limit the response to this constant need to adjust to fluctuating selection, increasing the fitness burden and thus the risk of extinction. Here, we use long-term monitoring data to assess whether the type of selection for a key fitness trait (i.e., lay date) differs between two reintroduced populations of a threatened passerine bird, the hihi (Notiomystis cincta). We apply recent statistical developments to test for the presence or absence of fluctuation in selection in both the Tiritiri Mātangi Island and the Kārori sanctuary populations. Our results support the presence of stabilizing selection in Tiritiri Mātangi with a potential moving optimum for lay date. In Kārori our results favour a regime of directional selection. Although the shape of selection may differ, for both populations an earlier lay date generally increases fitness in both environments. Further, the moving optimum models of lay date on Tiritiri Mātangi, suggesting that selection varies between years, imply a substantial lag load in addition to the fitness burden caused by the population laying too late. Our results highlight the importance of characterizing the form and temporal variation of selection for each population to predict the effects of environmental change and to inform management.

Does mast seeding shape mating time in wild boar? A comparative study

In seasonal environments, the timing of reproduction often matches with the peak of food resources. One well-known effect of global warming is an earlier phenology of resources, leading to a possible mismatch between the timing of reproduction for consumers and food peak. However, global warming may also change the dynamics of food resources, such as the intensity and frequency of pulsed mast seeding. How quantitative changes in mast seeding influence the timing of reproduction of seed consumers remains unexplored. Here, we assess how yearly variation in mast seeding influences mating time in wild boar (Sus scrofa), a widespread seed consumer species. We took advantage of the intensive monitoring of both female reproduction (1636 females) and acorn production over 6 consecutive years across 15 populations of wild boar in the wild. We found that mating time occurs earlier when acorn production increases in most but not all populations. In two out of 15 populations, heavy females mated earlier than light ones. Our findings demonstrate that mast seeding advances the mating time in some populations, which could perhaps impact how boars respond to climate change.

Natural selection fluctuates at an extremely fine spatial scale inside a wild population of snapdragon plants

Spatial variation in natural selection is expected to shape phenotypic variation of wild populations and drive their evolution. Although evidence of phenotypic divergence across populations experiencing different selection regimes is abundant, investigations of intrapopulation variation in selection pressures remain rare. Fine-grained spatial environmental heterogeneity can be expected to influence selective forces within a wild population and thereby alter its fitness function by producing multiple fitness optima at a fine spatial scale. Here, we tested this hypothesis in a wild population of snapdragon plants living on an extremely small island in southern France (about 7500 m² ). We estimated the spline-based fitness function linking individuals' fitness and five morphological traits in interaction with three spatially variable ecological drivers. We found that selection acting on several traits varied both in magnitude and direction in response to environmental variables at the scale of a meter. Our findings illustrate how different phenotypes can be selected at different locations within a population in response to environmental variation. Investigating spatial variation in selection within a population, in association with ecological conditions, represents an opportunity to identify putative ecological drivers of selection in the wild.

Directional selection and the evolution of breeding date in birds, revisited: Hard selection and the evolution of plasticity

The mismatch between when individuals breed and when we think they should breed has been a long‐standing problem in evolutionary ecology. Price et al. is a classic theory paper in this field and is mainly cited for its most obvious result: if individuals with high nutritional condition breed early, then the advantage of breeding early may be overestimated when information on nutritional condition is absent. Price at al.'s less obvious result is that individuals, on average, are expected to breed later than the optimum. Here, we provide an explanation of their non‐intuitive result in terms of hard selection, and go on to show that neither of their results are expected to hold if the relationship between breeding date and nutrition is allowed to evolve. By introducing the assumption that the advantage of breeding early is greater for individuals in high nutritional condition, we show that their most cited result can be salvaged. However, individuals, on average, are expected to breed earlier than the optimum, not later. More generally, we also show that the hard selection mechanisms that underpin these results have major implications for the evolution of plasticity: when environmental heterogeneity becomes too great, plasticity is selected against, prohibiting the evolution of generalists.

Frequency dependence and the predictability of evolution in a changing environment

Frequency‐dependent (FD) selection, whereby fitness and selection depend on the genetic or phenotypic composition of the population, arises in numerous ecological contexts (competition, mate choice, crypsis, mimicry, etc.) and can strongly impact evolutionary dynamics. In particular, negative frequency‐dependent selection (NFDS) is well known for its ability to potentially maintain stable polymorphisms, but it has also been invoked as a source of persistent, predictable frequency fluctuations. However, the conditions under which such fluctuations persist are not entirely clear. In particular, previous work rarely considered that FD is unlikely to be the sole driver of evolutionary dynamics when it occurs, because most environments are not static but instead change dynamically over time. Here, we investigate how FD interacts with a temporally fluctuating environment to shape the dynamics of population genetic change. We show that a simple metric introduced by Lewontin, the slope of frequency change against frequency near equilibrium, works as a key criterion for distinguishing microevolutionary outcomes, even in a changing environment. When this slope D is between 0 and –2 (consistent with the empirical examples we review), substantial fluctuations would not persist on their own in a large population occupying a constant environment, but they can still be maintained indefinitely as quasi‐cycles fueled by environmental noise or genetic drift. However, such moderate NFDS buffers and temporally shifts evolutionary responses to periodic environments (e.g., seasonality). Stronger FD, with slope D < –2, can produce self‐sustained cycles that may overwhelm responses to a changing environment, or even chaos that fundamentally limits predictability. This diversity of expected outcomes, together with the empirical evidence for both FD and environment‐dependent selection, suggests that the interplay of internal dynamics with external forcing should be investigated more systematically to reach a better understanding and prediction of evolution.

Phenotypic plasticity drives phenological changes in a Mediterranean blue tit population

Earlier phenology induced by climate change, such as the passerines' breeding time, is observed in many natural populations. Understanding the nature of such changes is key to predict the responses of wild populations to climate change. Genetic changes have been rarely investigated for laying date, though it has been shown to be heritable and under directional selection, suggesting that the trait could evolve. In a Corsican blue tit population, the birds' laying date has significantly advanced over 40 years, and we here determine whether this response is of plastic or evolutionary origin, by comparing the predictions of the breeder's and the Robertson-Price (STS) equations, to the observed genetic changes. We compare the results obtained for two fitness proxies (fledgling and recruitment success), using models accounting for their zero inflation. Because the trait appears heritable and under directional selection, the breeder's equation predicts that genetic changes could drive a significant part of the phenological change observed. We, however, found that fitness proxies and laying date are not genetically correlated. The STS, therefore, predicts no evolution of the breeding time, predicting correctly the absence of trend in breeding values. Our results also emphasize that when investigating selection on a plastic trait under fluctuating selection, part of the fitness-trait phenotypic covariance can be due to within individual covariance. In the case of repeated measurements, splitting within and between individual covariance can shift our perspective on the actual intensity of selection over multiple selection episodes, shedding light on the potential for the trait to evolve.

Using ecological context to interpret spatiotemporal variation in natural selection

Spatiotemporal variation in natural selection is expected, but difficult to estimate. Pollinator-mediated selection on floral traits provides a good system for understanding and linking variation in selection to differences in ecological context. We studied pollinator-mediated selection in five populations of Dalechampia scandens (Euphorbiaceae) in Costa Rica and Mexico. Using a nonlinear path-analytical approach, we assessed several functional components of selection, and linked variation in pollinator-mediated selection across time and space to variation in pollinator assemblages. After correcting for estimation error, we detected moderate variation in net selection on two out of four blossom traits. Both the opportunity for selection and the mean strength of selection decreased with increasing reliability of cross-pollination. Selection for pollinator attraction was consistently positive and stronger on advertisement than reward traits. Selection on traits affecting pollen transfer from the pollinator to the stigmas was strong only when cross-pollination was unreliable and there was a mismatch between pollinator and blossom size. These results illustrate how consideration of trait function and ecological context can facilitate both the detection and the causal understanding of spatiotemporal variation in natural selection.

Hydrology influences breeding time in the white-throated dipper

BACKGROUND

Earlier breeding is one of the strongest responses to global change in birds and is a key factor determining reproductive success. In most studies of climate effects, the focus has been on large-scale environmental indices or temperature averaged over large geographical areas, neglecting that animals are affected by the local conditions in their home ranges. In riverine ecosystems, climate change is altering the flow regime, in addition to changes resulting from the increasing demand for renewable and clean hydropower. Together with increasing temperatures, this can lead to shifts in the time window available for successful breeding of birds associated with the riverine habitat. Here, we investigated specifically how the environmental conditions at the territory level influence timing of breeding in a passerine bird with an aquatic lifestyle, the white-throated dipper Cinclus cinclus. We relate daily river discharge and other important hydrological parameters, to a long-term dataset of breeding phenology (1978-2015) in a natural river system.

RESULTS

Dippers bred earlier when winter river discharge and groundwater levels in the weeks prior to breeding were high, and when there was little snow in the catchment area. Breeding was also earlier at lower altitudes, although the effect dramatically declined over the period. This suggests that territories at higher altitudes had more open water in winter later in the study period, which permitted early breeding also here. Unexpectedly, the largest effect inducing earlier breeding time was territory river discharge during the winter months and not immediately prior to breeding. The territory river discharge also increased during the study period.

CONCLUSIONS

The observed earlier breeding can thus be interpreted as a response to climate change. Measuring environmental variation at the scale of the territory thus provides detailed information about the interactions between organisms and the abiotic environment.

Fluctuating optimum and temporally variable selection on breeding date in birds and mammals

Temporal variation in natural selection is predicted to strongly impact the evolution and demography of natural populations, with consequences for the rate of adaptation, evolution of plasticity, and extinction risk. Most of the theory underlying these predictions assumes a moving optimum phenotype, with predictions expressed in terms of the temporal variance and autocorrelation of this optimum. However, empirical studies seldom estimate patterns of fluctuations of an optimum phenotype, precluding further progress in connecting theory with observations. To bridge this gap, we assess the evidence for temporal variation in selection on breeding date by modeling a fitness function with a fluctuating optimum, across 39 populations of 21 wild animals, one of the largest compilations of long-term datasets with individual measurements of trait and fitness components. We find compelling evidence for fluctuations in the fitness function, causing temporal variation in the magnitude, but not the direction of selection. However, fluctuations of the optimum phenotype need not directly translate into variation in selection gradients, because their impact can be buffered by partial tracking of the optimum by the mean phenotype. Analyzing individuals that reproduce in consecutive years, we find that plastic changes track movements of the optimum phenotype across years, especially in bird species, reducing temporal variation in directional selection. This suggests that phenological plasticity has evolved to cope with fluctuations in the optimum, despite their currently modest contribution to variation in selection.

Fluctuations in lifetime selection in an autocorrelated environment

Most natural environments vary stochastically and are temporally autocorrelated. Previous theory investigating the effects of environmental autocorrelation on evolution mostly assumed that total fitness resulted from a single selection episode. Yet organisms are likely to experience selection repeatedly along their life, in response to possibly different environmental states. We model the evolution of a quantitative trait in organisms with non-overlapping generations undergoing several episodes of selection in a randomly fluctuating and autocorrelated environment. We show that the evolutionary dynamics depends not directly on fluctuations of the environment, but instead on those of an effective phenotypic optimum that integrates the effects of all selection episodes within each generation. The variance and autocorrelation of the integrated optimum shape the variance and predictability of selection, with substantial qualitative and quantitative deviations from previous predictions considering a single selection episode per generation. We also investigate the consequence of multiple selection episodes per generation on population load. In particular, we identify a new load resulting from within-generation fluctuating selection, generating the death of individuals without significance for the evolutionary dynamics. Our study emphasizes how taking into account fluctuating selection within lifetime unravels new properties of evolutionary dynamics, with crucial implications notably with respect to responses to global changes.

Where is the optimum? Predicting the variation of selection along climatic gradients and the adaptive value of plasticity. A case study on tree phenology

Many theoretical models predict when genetic evolution and phenotypic plasticity allow adaptation to changing environmental conditions. These models generally assume stabilizing selection around some optimal phenotype. We however often ignore how optimal phenotypes change with the environment, which limit our understanding of the adaptive value of phenotypic plasticity. Here, we propose an approach based on our knowledge of the causal relationships between climate, adaptive traits, and fitness to further these questions. This approach relies on a sensitivity analysis of the process‐based model phenofit, which mathematically formalizes these causal relationships, to predict fitness landscapes and optimal budburst dates along elevation gradients in three major European tree species. Variation in the overall shape of the fitness landscape and resulting directional selection gradients were found to be mainly driven by temperature variation. The optimal budburst date was delayed with elevation, while the range of dates allowing high fitness narrowed and the maximal fitness at the optimum decreased. We also found that the plasticity of the budburst date should allow tracking the spatial variation in the optimal date, but with variable mismatch depending on the species, ranging from negligible mismatch in fir, moderate in beech, to large in oak. Phenotypic plasticity would therefore be more adaptive in fir and beech than in oak. In all species, we predicted stronger directional selection for earlier budburst date at higher elevation. The weak selection on budburst date in fir should result in the evolution of negligible genetic divergence, while beech and oak would evolve counter‐gradient variation, where genetic and environmental effects are in opposite directions. Our study suggests that theoretical models should consider how whole fitness landscapes change with the environment. The approach introduced here has the potential to be developed for other traits and species to explore how populations will adapt to climate change.

Stabilizing selection and adaptive evolution in a combination of two traits in an arctic ungulate

Stabilizing selection is thought to be common in wild populations and act as one of the main evolutionary mechanisms, which constrain phenotypic variation. When multiple traits interact to create a combined phenotype, correlational selection may be an important process driving adaptive evolution. Here, we report on phenotypic selection and evolutionary changes in two natal traits in a semidomestic population of reindeer (Rangifer tarandus) in northern Finland. The population has been closely monitored since 1969, and detailed data have been collected on individuals since they were born. Over the length of the study period (1969-2015), we found directional and stabilizing selection toward a combination of earlier birth date and heavier birth mass with an intermediate optimum along the major axis of the selection surface. In addition, we demonstrate significant changes in mean traits toward earlier birth date and heavier birth mass, with corresponding genetic changes in breeding values during the study period. Our results demonstrate evolutionary changes in a combination of two traits, which agree closely with estimated patterns of phenotypic selection. Knowledge of the selective surface for combinations of genetically correlated traits are vital to predict how population mean phenotypes and fitness are affected when environments change.

Selective Sweep at a QTL in a Randomly Fluctuating Environment

Adaptation is mediated by phenotypic traits that are often near continuous, and undergo selective pressures that may change with the environment. The dynamics of allelic frequencies at underlying quantitative trait loci (QTL) depend on their own phenotypic effects, but also possibly on other polymorphic loci affecting the same trait, and on environmental change driving phenotypic selection. Most environments include a substantial component of random noise, characterized both by its magnitude and its temporal autocorrelation, which sets the timescale of environmental predictability. I investigate the dynamics of a mutation affecting a quantitative trait in an autocorrelated stochastic environment that causes random fluctuations of an optimum phenotype. The trait under selection may also exhibit background polygenic variance caused by many polymorphic loci of small effects elsewhere in the genome. In addition, the mutation at the QTL may affect phenotypic plasticity, the phenotypic response of given genotype to its environment of development or expression. Stochastic environmental fluctuations increase the variance of the evolutionary process, with consequences for the probability of a complete sweep at the QTL. Background polygenic variation critically alters this process, by setting an upper limit to stochastic variance of population genetics at the QTL. For a plasticity QTL, stochastic fluctuations also influences the expected selection coefficient, and alleles with the same expected trajectory can have very different stochastic variances. Finally, a mutation may be favored through its effect on plasticity despite causing a systematic mismatch with optimum, which is compensated by evolution of the mean background phenotype.

The evolution of phenotypic plasticity when environments fluctuate in time and space

Most theoretical studies have explored the evolution of plasticity when the environment, and therefore the optimal trait value, varies in time or space. When the environment varies in time and space, we show that genetic adaptation to Markovian temporal fluctuations depends on the between-generation autocorrelation in the environment in exactly the same way that genetic adaptation to spatial fluctuations depends on the probability of philopatry. This is because both measure the correlation in parent-offspring environments and therefore the effectiveness of a genetic response to selection. If the capacity to genetically respond to selection is stronger in one dimension (e.g., space), then plasticity mainly evolves in response to fluctuations in the other dimension (e.g., time). If the relationships between the environments of development and selection are the same in time and space, the evolved plastic response to temporal fluctuations is useful in a spatial context and genetic differentiation in space is reduced. However, if the relationships between the environments of development and selection are different, the optimal level of plasticity is different in the two dimensions. In this case, the plastic response that evolves to cope with temporal fluctuations may actually be maladaptive in space, resulting in the evolution of hyperplasticity or negative plasticity. These effects can be mitigated by spatial genetic differentiation that acts in opposition to plasticity resulting in counter-gradient variation. These results highlight the difficulty of making space-for-time substitutions in empirical work but identify the key parameters that need to be measured in order to test whether space-for-time substitutions are likely to be valid.

Short-term insurance versus long-term bet-hedging strategies as adaptations to variable environments

Understanding how organisms adapt to environmental variation is a key challenge of biology. Central to this are bet‐hedging strategies that maximize geometric mean fitness across generations, either by being conservative or diversifying phenotypes. Theoretical models have identified environmental variation across generations with multiplicative fitness effects as driving the evolution of bet‐hedging. However, behavioral ecology has revealed adaptive responses to additive fitness effects of environmental variation within lifetimes, either through insurance or risk‐sensitive strategies. Here, we explore whether the effects of adaptive insurance interact with the evolution of bet‐hedging by varying the position and skew of both arithmetic and geometric mean fitness functions. We find that insurance causes the optimal phenotype to shift from the peak to down the less steeply decreasing side of the fitness function, and that conservative bet‐hedging produces an additional shift on top of this, which decreases as adaptive phenotypic variation from diversifying bet‐hedging increases. When diversifying bet‐hedging is not an option, environmental canalization to reduce phenotypic variation is almost always favored, except where the tails of the fitness function are steeply convex and produce a novel risk‐sensitive increase in phenotypic variance akin to diversifying bet‐hedging. Importantly, using skewed fitness functions, we provide the first model that explicitly addresses how conservative and diversifying bet‐hedging strategies might coexist.