Adaptation, plasticity, and extinction in a changing environment: towards a predictive theory — TEI-tagged view

Luis-Miguel ChevinpersName, Russell LandepersName, Georgina M. MacepersName (2010date). “Adaptation, plasticity, and extinction in a changing environment: towards a predictive theory.” PLoS Biologytitle·j.

Object of study — population persistenceterm·dom and extinction under sustained directional environmental change

Research question — Which developmental, genetic, and demographic parameters determine the maximum sustained rate of environmental changeterm an isolated population can survive?

Contribution — a predictive evolutionary-demographic model integrating phenotypic plasticityterm·dom and an explicit environmental variable with a research program for estimating its parameters

Full text — tagged

Adaptation, plasticity, and extinction in a changing environment: towards a predictive theory Abstract The authors analyze developmental, genetic, and demographic mechanisms by which populations tolerate changing environments and discuss empirical methods for determining the critical rate of sustained environmental changeterm·dom that causes population extinctionterm·dom.

Summary Many species are experiencing sustained environmental changeterm·dom mainly due to human activities. The unusual rate and extent of anthropogenic alterations of the environment may exceed the capacity of developmental, genetic, and demographic mechanisms that populations have evolved to deal with environmental change. To begin to understand the limits to population persistenceterm·dom, we present a simple evolutionary model for the critical rate of environmental changeterm·dom beyond which a population must decline and go extinct. We use this model to highlight the major determinants of extinction riskterm·dom in a changing environment, and identify research needs for improved predictions based on projected changes in environmental variables. Two key parameters relating the environment to population biology have not yet received sufficient attention. Phenotypic plasticityterm·dom, the direct influence of environment on the development of individual phenotypes, is increasingly considered an important component of phenotypic change in the wild and should be incorporated in models of population persistenceterm·dom. Environmental sensitivity of selectionterm·dom, the change in the optimum phenotypeterm with the environment, still crucially needs empirical assessment. We use environmental tolerance curvesterm·dom and other examples of ecological and evolutionary responses to climate change to illustrate how these mechanistic approaches can be developed for predictive purposes.

Introduction Global climate change, over-exploitation, and habitat alteration are causing sustained and consistent pressures on wild populations [1] . Because of habitat fragmentation, many species will not be able to track their preferred environment in space, and must therefore adapt in situ to avoid extinction. What determines the maximum rate of environmental changeterm that populations can cope with? Understanding this will inform both models and management plans about where critical thresholds lie or what might affect the potential resilience of particular species or ecological communities.

Two main approaches exist for studying the impact of environmental change on species persistence—niche modellingterm·mth and mechanistic population modellingterm·mth. On the one hand, “climate envelope modelsterm·mth” (or “niche models”) are correlative and focused on the environment. Their conceptual background traces to Hutchinson's multidimensional representation of the niche [2] . They use measurements of environmental variables and records of species presence and absence to infer abiotic correlates of the realized niche of a species. By projecting this niche on a map with environmental data, the spatial locations that satisfy the basic requirements of the species are identified. Niche models have been combined with climate projections to predict range shifts and extinction rates [3] . Aside from the methodological caveats of this approach (changing covariance of environmental variables [4] , population demography not in equilibrium with changing climate [5] , spatial scale of the analysis [6] , [7] , and source-sink dynamics allowing individuals to exist outside their niche [8] ), its major drawback is that it does not account for the biological processes underlying adaptation of a species to its environment. It cannot identify or incorporate the biological differences among species that determine whether or not they can persist in situ. While a new generation of “process-based” (or “mechanistic”) niche models is being developed to overcome these limitations [9] – [11] , they currently rely on simplified demographic processes and generally overlook evolution in response to climate change (but see [12] ). Niche modellingterm·mth is thus currently more suited to understand changes in the distribution of suitable environments for a species at a continental scale than to predict population persistenceterm·dom or guide conservation plans.

On the other hand, mechanistic population modellingterm·mth focuses on the biological processes that underlie population persistenceterm·dom. By combining evolutionary genetics and demography, the conditions that allow a population to maintain a positive intrinsic growth rate in the face of environmental change can be predicted theoretically [13] , [14] and studied empirically [15] . The main focus of mechanistic population studies is not environmental variables, but ecologically important phenotypic traits (morphology, physiology, phenology, or behaviour) that affect population growth through their influence on life history. In contrast to niche modellingterm·mth, this approach allows identification of factors that potentially limit population persistenceterm·dom in a changing environment. However, unless environmental variables are included, as well as phenotypic traits, this approach cannot be used to project the fate of populations.

Here, we propose a model of evolution and population growthterm·mth to address how demographic and evolutionary constraints limit the persistence of a geographically isolated population under sustained environmental changeterm·dom. This model generalizes an earlier one by Lynch and LandepersName [13] by including phenotypic plasticityterm·dom (see Box 1–Glossary) and an explicit environmental variable such as temperature. We use it to highlight the determinants of extinction riskterm·dom and to identify research needs for a more predictive approach to population persistenceterm·dom in a changing environment. Research directed towards key determinants of extinction riskterm·dom that are often overlooked in current empirical work may allow population biologists to take advantage of abundant environmental data in order to begin developing predictions about population persistenceterm·dom under alternative scenarios of climate change [16] , which can help to guide conservation policies and management plans.

Box 1. Glossary critical rate of environmental changeterm·dom: maximum rate of sustained environmental changeterm·dom that allows long-term persistence of a population, denoted as η c .

environmental sensitivity of selectionterm·dom: change in the optimum phenotypeterm with the environment. For a linear relationship, it is measured by the slope Bterm .

generation timeterm·dom: the average age of parents of a cohort of newborn individuals, denoted as T . With discrete non-overlapping generations, T is the mean time between successive reproductive episodes in the population.

genetic varianceterm·dom: the genetic component of phenotypic varianceterm σ 2 , or more precisely, “additive genetic varianceterm·dom” h 2 σ 2 , the statistically additive component of phenotypic varianceterm determining the resemblance between offspring and parents, and the genetic response to selection.

heritabilityterm·dom: the proportion of phenotypic varianceterm in a trait due to additive genetic effects, denoted as h 2 .

intrinsic rate of increaseterm·dom: population growth rate in the absence of intra- or inter-specific competition. For a perfectly adapted population with mean phenotypeterm at the optimum, the intrinsic rate of increaseterm·dom is denoted as r max .

norm of reactionterm·dom: the expected phenotype of a given genotype as a function of the environmentterm·dom.

phenotypic plasticityterm·dom: direct influence of the environment on individual phenotypesterm·dom through developmental mechanisms. For a linear norm of reactionterm·dom plasticity is measured by the slope bterm .

quantitative traitterm·dom: continuously distributed phenotypic character, with phenotypic varianceterm in a population determined by multiple polymorphic genes and environmental effects.

stabilizing selectionterm·dom: natural selection such that individual fitness decreases with increasing phenotypic deviation from an optimum. Its strength is measured by γ .

tolerance curveterm·dom: fitness or performance as a function of the environment.

Extinction Riskterm·dom under Sustained Environmental Changeterm·dom There are three mechanisms by which a population can persist when its local environment changes: dispersalterm·dom to track its preferred environment in space, genetic evolutionterm·dom to the new local conditions, or phenotypic plasticityterm·dom [17] , [18] . Using the model in Box 2 and Text S1 , we give an overview of the determinants of extinction riskterm·dom for an isolated population with a restricted geographic range under sustained environmental changeterm·dom, assuming dispersalterm·dom is not possible. We consider cases where all other possible factors influencing extinction riskterm·dom (demographic and environmental stochasticity, inbreeding depression, genetic drift [19] ) can be ignored relative to the impact of sustained environmental changeterm·dom. The determinants of extinction riskterm·dom then belong to three main classes.

Box 2. A Model of Phenotypic Evolution and Persistence under Sustained Environmental Changeterm·dom We assume that a continuous environmental parameterterm ε (e.g., temperature or precipitation) changes at a constant rate η in time. Adaptation to this changing environment is mediated by a quantitative traitterm·dom z that determines fitness. Selection on multiple correlated traits also can be incorporated [20] , but for simplicity we focus on single trait. Population growth is assumed to be density-independentterm, which can be a good approximation even with density dependence, if the population remains well below the carrying capacity because of maladaptation. In contrast to the original model and its recent extensions [13] , [20] , we include phenotypic plasticityterm·dom in the trait z , such that a given change in the environment directly modifies the phenotype of each individual by a constant amount. The rate of environmental changeterm is expressed per unit time (e.g., °C per year) instead of per generation.

Under sustained environmental changeterm·dom, assuming constant genetic and phenotypic varianceterm and strength of stabilizing selectionterm, the rate of phenotypic evolution eventually reaches a steady state where the mean phenotypeterm lags a constant distance behind the changing optimum phenotypeterm ( Text S1 ). If this distance is such that the population has a negative growth rate, extinction occurs. The maximum rate of environmental changeterm that allows long-term persistence of the population (or critical rate of environmental changeterm·dom) is The original model [13] revealed how the phenotypic varianceterm σ 2 and heritabilityterm·dom h 2 of the trait, the strength of stabilizing selectionterm γ , and the maximum intrinsic rate of population growthterm r max determined the critical rate of environmental changeterm·dom in the absence of phenotypic plasticityterm·dom.

Three other parameters included here emphasize important factors affecting persistence when the directional environmental change is expressed in environmental units per unit time. First is the generation timeterm·dom T . The remaining two parameters are the environmental sensitivity of selectionterm·dom B , which measures how changes in the environment influence the optimum phenotypeterm, and phenotypic plasticityterm·dom b , which quantifies the direct impact of the environment on development of individual phenotypes. The critical rate of environmental changeterm·dom for long-term persistence increases with decreasing absolute difference between the environmental sensitivity of selectionterm·dom and phenotypic plasticityterm·dom. Although plasticity causes weaker natural selection on the trait and smaller genetic response to selection, this is more than compensated by the plastic phenotypic change that brings the mean phenotypeterm closer to the optimum ( Text S1 ).

Plasticity may also entail a fitness cost [59] . The cost of plasticityterm·dom decreases r max , which can then be formulated as where r 0 is the per-generation growth rate of a population with optimum phenotypeterm and no plasticity, and quantifies the magnitude of the cost for a given plasticity b . By decreasing r max , the cost of plasticityterm·dom decreases the critical rate of environmental changeterm·dom, thus opposing the beneficial effect of plasticity on η c . The overall effect of phenotypic plasticityterm·dom on the critical rate of environmental changeterm·dom is shown in Figure 2 for several values of the cost of plasticityterm·dom. When the cost of plasticityterm·dom is above a threshold ( ), there is an intermediate value of plasticity that maximizes the critical rate of environmental changeterm·dom.

(1) Demographic propertiesterm of the population . These are components of the population growth rate that are not affected by adaptation to the changing environment. They include the generation timeterm·dom, T , and the maximum intrinsic rate of increaseterm·dom r max (per unit time) of a population with optimum mean phenotypeterm. Even if these parameters are not affected by the environmental change, they do constrain the rate of adaptation. Populations with longer generations must evolve faster per generation to adapt to a given rate of environmental changeterm, and populations with a low r max will reach extinction before they can adapt to rapid environmental change [13] – [15] , [20] . The specific influence of each of these variables on extinction riskterm·dom across taxa also depends on their covariation. Broad-scale studies of life-history traits indicate that the per-generation intrinsic rate of increaseterm·dom is of the same order of magnitude across several taxa, whereas r max per unit time is roughly inversely proportional to the generation timeterm·dom [21] – [23] . As a consequence, the critical rate of environmental changeterm·dom should be generally lower for species with longer generation timesterm·dom, which will be at greater extinction riskterm·dom.

(2) Evolutionary potentialterm . This measures how fast genetic evolutionterm·dom occurs for a given deviation of the mean phenotypeterm from the optimum. The strength of stabilizing selectionterm γ measures how the mean fitness decreases with the squared deviation from the optimum and also determines the strength of directional selection for a given deviation from the optimum. Stronger stabilizing selectionterm·dom (larger γ ), although causing increased mortality, allows faster environmental change to be tolerated by causing faster evolutionterm·dom [13] . The genetic varianceterm·dom, the product of phenotypic varianceterm σ 2 and heritabilityterm·dom h 2 , determines how much genetic evolutionterm·dom is caused by a given strength of selection. Higher genetic varianceterm·dom allows persistence under stronger environmental changeterm·dom, because it allows the population to track its phenotypic optimum more closely. Although we focus on a single trait for simplicity, this may actually be a linear combination of measurable traits. In this case, the genetic variancesterm·dom of the original traits and their genetic covariances affect the magnitude of the genetic response to selection in any phenotypic direction [24] , [25] .

(3) Biological impact of the environmentterm . This links the environment to the biology of individuals in the species. Phenotypic plasticityterm·dom describes the direct impact of the environment on the development of individual phenotypes. It may involve morphological, physiological or behavioural responses, which can occur on different time scales. For continuous environmental variables, plasticity usually is modelled using reaction norms, where the phenotype of a given genotype is plotted as a function of the environment [26] . We focus on linear reaction normsterm for simplicity, although reaction norms can be non-linear [27] , [28] (see below). With linear reaction normsterm, the slope bterm quantifies the degree of plasticity. The environmental sensitivity of selectionterm·dom, B , measures how the optimum phenotypeterm changes with the environment, which for simplicity we also assume is a linear relationship. With no cost of plasticityterm·dom, populations with b closer to B are likely to persist under higher rates of environmental change.

A Research Program This model can be combined with environmental projectionsterm·mth to ask whether future rates of environmental change will threaten the persistence of particular populations or species. Application of the model requires multiple steps, which previously have been undertaken in isolation or in combination, although all have rarely been completed together.

Identify Ecologically Important Traitsterm·dom Affecting Population Growth Ecological investigation in the field is needed to identify traits that potentially determine adaptation to a specific environment. For instance, for ectotherms such as insects and reptiles, thermal adaptation may occur not only through physiological traits governing energy metabolism, but also through behavioural and morphological traits involved in movement between shaded and sunny patches [29] . For many bird species, adaptation to global warming involves adjusting their breeding date so that reproduction coincides with a peak in prey density [30] .

The strength (and direction) of natural selection γ has been estimated in many plant and animal taxa by regressing the fitness of individuals on their phenotypes [31] , [32] . Most of these studies focus on relative fitness, which determines the evolutionary dynamics. However, absolute fitnessterm·dom and r max are needed to understand how phenotypic traits influence population dynamics. Analysis of natural selection has been extended to age- or stage-structured populations, allowing identification of parts of the life cycle where phenotypes most strongly influence demography [33] – [36] .

Estimate Genetic Varianceterm·dom and Plasticity The question of whether genetic variation ( h 2 σ 2 ) along the direction of natural selection limits the rate of evolution has motivated many studies. For instance, in the North American prairie plant Chamaecrista fasciculate , small genetic varianceterm·dom for a linear combination of traits has been predicted to limit the response to selection caused by climate change [25] . Similarly, it was suggested that lack of genetic variation for desiccation and cold resistance restricts the geographic range of several Drosophila species [37] and would limit the response to projected climate change in one of them [38] . Changes in the genetic and phenotypic variation with the environment may also be important. For example, the response to selection on body weight in wild Soay sheep is reduced because of lower heritabilityterm·dom in environments where selection is stronger [39] .

The importance of phenotypic plasticityterm·dom for rapid evolution in response to directional environmental change has only recently been fully appreciated. Many authors have observed adaptive phenotypic evolution without being able to firmly attribute it to genetic change (reviewed in [40] ), but relatively few studies have progressed from considering plasticity as a nuisance parameter to measuring norms of reaction, and their slopes b , in the context of climate change [41] . This was done recently for reproductive timing in mammals [42] and birds [43] , [44] .

Measure the Environmental Sensitivity of Selectionterm·dom Few investigators have measured the environmental sensitivity of phenotypic selection as described by parameter B in Box 2. The most relevant studies have measured selection in a small number of discrete environments. This approach has been applied to birds to describe the influence of changing seasonality on selection on breeding date [45] , and in plants to examine the impact of drought on selection on physiological traits [46] , [47] . It has also been combined with experimental manipulation of the environment, mainly in plants, using either abiotic [48] or biotic factors, such as crowding [49] , [50] or natural enemies [51] , [52] . Other studies have employed ecological criteria to identify changes in the optimum phenotypeterm without actually measuring selection. For instance, the optimum egg-laying date for birds can be deduced from the temporal peak in food availability determined by the population dynamics of insects on which they feed [53] .

What emerges from this brief outline is that two key areas of study are lacking: the quantification of the environmental sensitivity of selectionterm·dom across continuous environments and the detailed investigation of phenotypic plasticityterm·dom, including its costs and limits under sustained environmental changeterm·dom. Further work, such as that outlined below, would begin to make predictive approaches possible.

Phenotypic Selection along an Environmental Gradient A relationship between an environmental variable and a measure of selection, either observed in the wild or derived from an experiment, has been argued to be the strongest evidence for the cause of natural selection [54] . Quantification of this relationship is required to predict population evolution and persistence in a changing environment. Most of the studies highlighted above compared selective pressures among discrete environmental states, often interpreted as stresses (but see [39] ). Current global change, however, mainly involves continuous environmental variables for which precise time series and projections are available. The environmental sensitivity of selectionterm·dom should therefore be estimated by measuring phenotypic selection along a continuous environmental gradientterm·mth. Apart from measuring a fundamental parameter of the model, this would also allow its basic assumptions to be tested. In particular, our model (like previous ones) assumes that environmental change primarily alters the optimum phenotypeterm with little impact on the width of the fitness function ( Figure 1 ), but this assumption has hardly ever been tested empirically.

Tolerance curvesterm·dom and phenotypic plasticityterm·dom. First row : norms of reaction for three genotypes (colored lines); second row : generalized adaptive landscape depicting fitness as a function of the phenotype and the environment, with the reaction norms projected up onto the fitness surface from the lower plane; third row : environmental tolerance curvesterm·dom representing slices through the adaptive landscape along the lines defined by the reaction norms. The three columns represent alternative scenarios described in the text. In the lower left panel, the dashed green and dashed red tolerance curvesterm·dom include costs of plasticity.

Measurement of selection in several environments may be difficult to carry out in the field because of lack of control over environmental parametersterm. For some species such measurements could be performed in controlled laboratory conditions, provided the relationship between fitness in the laboratory and in the wild can be established. Alternatively, performance can be measured in the laboratory and the relationship between performance and fitness assessed in the wild [55] .

Studies of thermal tolerance curvesterm·dom in ectotherms commonly measure fitness in continuous environments. However, these studies generally relate fitness (or performance) directly to the environment, which makes it difficult to produce evolutionary predictions and to test them experimentally [56] . Modelling tolerance curvesterm·dom with an underlying phenotypic trait subject to natural selection and phenotypic plasticityterm·dom may prove more fruitful, since it should allow predictions about the dynamics of tolerance curveterm·dom evolution (Box 3).

Box 3. Tolerance Curvesterm·dom and Phenotypic Plasticityterm·dom Environmental tolerance curvesterm·dom relating fitness (or performance) to abiotic environmental variables such as temperature [56] or salinity [77] have been studied for a wide variety of organisms, ranging from reptiles and amphibians to insects and plants [56] . They are generally characterized by a few descriptive variables, such as critical environments (where fitness or performance vanishes), optimum environmentterm·dom (where fitness is maximized), and skewness. Previous models investigated the benefits of broad versus narrow tolerance curvesterm·dom under several forms of environmental change [61] , [62] , [78] , but they failed to produce predictions about the dynamics of and constrains on tolerance curveterm·dom evolution. Recently there have been calls for new theoretical developments with clearer empirically testable outcomes [56] .

We propose that incorporating phenotypically plastic traits, which are under stabilizing selectionterm·dom through their influence on lifetime fitness [33] , [34] , [79] , can clarify how tolerance curvesterm·dom emerge and evolve ( Figure 1 ). With linear reaction normsterm and constant strength of stabilizing selectionterm across environments (ridge with constant width in the adaptive landscape), the reaction norm slope determines tolerance breadthterm·dom, while the reaction norm intercept (or elevation) determines the optimum environmentterm·dom (first column in Figure 1 ). With genetic variation in the slope of reaction norms, the tolerance breadthterm·dom may evolve as a consequence of evolving phenotypic plasticityterm·dom in the underlying trait. Including a cost of plasticityterm·dom, where fitness decreases with reaction norm slope regardless of the trait value [80] , [81] , produces a generalist–specialist trade-off between tolerance breadth and maximum fitnessterm·dom, corresponding to the intuitive idea that a “Jack-of-all-trades is a master of none.” This is illustrated in the bottom left panel of Figure 1 —with cost an increasing function of plasticity (dashed lines), the red dashed tolerance curveterm·dom has a lower maximum fitness than the green dashed curve, because the former reaction norm has greater plasticity. This approach can be used to make testable predictions in specific cases based on measures of the variability, inheritance, and plasticity of phenotypic traits.

Thermal tolerance curvesterm·dom generally are strongly skewed, with fitness decreasing steeply as temperature approaches the critical thermal maximum [56] . This could be caused by non-linear reaction normsterm·dom for underlying traits, as illustrated by the second column in Figure 1 . Plastic phenotypic responses that saturate (yellow genotype) or increase exponentially (green genotype) in extreme environments could both induce skewed tolerance curvesterm·dom. Alternatively, skewed tolerance curvesterm·dom could be caused by the strength of stabilizing selectionterm (width of the fitness ridge) changing with the environment (third column in Figure 1 ). Finally, the fitness function acting on the trait may itself be skewed in all environments, the height of the fitness ridge may depend on the environment, or the optimum phenotypeterm may change non-linearly with the environment (not shown). These five possible causes of skewness have different implications for the evolution of thermal tolerance in a changing environment, and need to be tested empirically.

Costs and Limits of Plasticity in a Changing Environment Plasticity in Extreme Environments The model in Box 2 assumes a linear reaction normterm, implying that a given amount of environmental change always produces the same plastic phenotypic change. Plasticity is generally studied in the context of environments that vary in space or fluctuate in time with a stationary distribution, but little is known about plastic responses outside the usual range of variation. Extreme environments may disrupt the plastic response, such that the reaction norm may take any shape in environments that were rarely encountered before [57] . This argument is based on the theory of the evolution of reaction norm in heterogeneous environments, assuming little genetic constraint on reaction norm shape. However, reaction norm shapes are likely to be genetically constrained. For many traits, plastic phenotypic responses should reach a physiological limitterm·dom and “saturate” in extreme environments [58] . For instance, body size or metabolic rate obviously cannot increase or decrease indefinitely under sustained environmental changeterm·dom. For other traits, disruption of homeostasis may cause the plastic response to be amplified in extreme environments. Both of these mechanisms could generate skewed tolerance curvesterm·dom like those observed for response to temperature, even when phenotypic selection is symmetric (Box 3).

Costs of Plasticity and the Generalist versus Specialist Trade-Off Phenotypic plasticityterm·dom may entail several types of fitness costs to the organism independent of the expressed phenotype [59] . Here, we shall only distinguish two kinds of costs of plasticity: constitutive and induced. A constitutive costterm·dom reduces individual fitness depending only on the degree of plasticity. This includes the cost of maintaining physiological machinery that allows phenotypic plasticityterm·dom or acquiring information about the environment [59] . The model in Box 2 focuses on such costs. In contrast an induced costterm·dom is a reduction in fitness that depends on the amount of plastic phenotypic change. This kind of cost is physiological and can be understood in terms of constraints on energy allocation. It applies mainly to labile traits that may change repeatedly during the lifetime of an individual, such as breeding date for birds [60] .

In the literature on tolerance curvesterm·dom, the cost of plasticityterm·dom is generally expressed as a trade-off between tolerance breadthterm·dom and fitness in the optimum environmentterm·dom, corresponding to the intuitive idea that the “Jack-of-all-trades is a master of none.” This has been modelled [61] , [62] , but the empirical evidence is still controversial [56] , [63] and based mainly on comparative data among taxa rather than on within-species variation (but see [64] ). We propose that this trade-off may be a consequence of the cost of plasticityterm·dom for an underlying phenotypic trait (Box 3, Figure 1 ).

In the model of Box 2 and Figure 2 , we show how a constitutive costterm·dom of plasticity can limit the critical rate of environmental changeterm·dom for population persistenceterm·dom. An induced costterm·dom would yield a similar result in this model since the environment changes at a constant rate. The cost of plasticityterm·dom may thus limit population persistenceterm·dom in a changing environment, yet little is known about its importance in natural populations. A recent review showed that such costs may be widespread but weak [65] . Further research is needed to measure costs of plasticity for organisms with long generations (such as large mammals and trees) or with life cycles that depend on seasonal timing (e.g., interaction of temperature and photoperiod [66] ), which rely most on plasticity to persist in a rapidly changing environment.

Critical rate of environmental changeterm·dom with costly phenotypic plasticityterm·dom. The maximum rate of environmental changeterm allowing long-term persistence of a population, η c , is plotted against plasticity b , for several values of the cost of plasticityterm·dom. For a given plasticity b , the cost increases with decreasing ω b . For each ω b , rates of environmental change higher than the corresponding line cause population extinctionterm·dom. Parameters: r max = 0.140, T = 1, γ = 1/51, B = 2, σ 2 = 1, and h 2 = 0.5.

Limitations of the Model Our model includes a number of simplifying assumptions. We assume a constant shape of the fitness function across environments, allowing only the optimum phenotypeterm to change with the environment. However, the strength of stabilizing selectionterm around the optimum phenotypeterm γ may depend on the environment ( Figure 1 , third column, second row). Furthermore the maximum fitness r max may generally be lower in more extreme environments because of physiological constraints, such that the ridge in Figure 1 would have an intermediate maximum. This was shown to occur in a recent comparative study of thermal tolerance curvesterm·dom among several insect taxa, where the maximum fitness was lower for species that had their optimum at a lower temperature [67] . Other demographic parameters including the generation timeterm·dom T may change with the environment. Genotype × environment interaction (or genetic variance in plasticityterm·dom) can cause the genetic varianceterm·dom ( h 2 σ 2 ) to depend on the environment [68] – [71] . Although these mechanisms are not included in our simple model, they all can be analyzed theoretically and studied empirically, and should be included where they are suspected to be important.

Plasticity itself may evolve if it varies genetically, and new environments can cause directional selection on plasticity [72] , [73] . However, more information is needed about plasticity and its inheritance in extreme environments before the evolution of plasticity can be included in the analysis of persistence under sustained environmental changeterm·dom.

Finally, we considered a population that cannot disperse nor receive migrants from other populations experiencing different environments. We focused on this situation because it is one that is most commonly overlooked in niche modellingterm·mth studies. Our model thus applies best to species with habitats restricted by dispersalterm·dom barriers or species which disperse slowly relative to the rate of environmental changeterm [74] . At the other extreme, a population may be able to persist by following its preferred environment in space without any evolution, as assumed in niche modellingterm·mth, but this is very unlikely for species occupying fragmented habitatsterm·dom. For populations that can both disperse and evolve in response to a changing environment, gene flowterm·dom among populations with different environments may limit local adaptation and restrict the geographic range of the population or even cause extinction [74] , [75] .

Conclusion Our aim was to describe an approach based on evolutionary and demographic mechanisms that can be used to make predictions on population persistenceterm·dom in a changing environment and to highlight the most important variables to measure. While this approach is obviously more costly and time-consuming than niche modellingterm·mth, its results are also likely to be more useful for specific purposes because it explicitly incorporates the factors that limit population response to environmental change.

The feasibility of such a mechanistic approach has been demonstrated by a few recent studies. Deutsch et al. [76] used thermal tolerance curvesterm·dom to predict the fitness consequence of climate change for many species of terrestrial insects across latitudes, but without explicitly considering phenotypic plasticityterm·dom or genetic evolutionterm·dom. Kearney et al. [12] combined biophysical models of energy transfers with measures of heritabilityterm·dom of egg desiccation to predict how climate change would affect the distribution of the mosquito Aedes aegiptii in Australia. Egg desiccation was treated as a threshold trait, but the possibility of phenotypic plasticityterm·dom or evolution of the threshold was not considered. These encouraging efforts call for more empirical studies where genetic evolutionterm·dom and phenotypic plasticityterm·dom are combined with demography to make predictions about population persistenceterm·dom in a changing environment. The simple approach we have outlined is a necessary step towards a more specific and comprehensive understanding of the influence of environmental change on population extinctionterm·dom.

Supporting Information A model of plasticity, evolution, and extinction in a changing environment.

(0.48 MB PDF) Click here for additional data file.

Acknowledgments We thank A.P. Hendry and M.A. Olalla-Tarraga for discussions.

Footnotes References Associated Data Supplementary Materials A model of plasticity, evolution, and extinction in a changing environment.

(0.48 MB PDF) Click here for additional data file.