In this paper we describe how to efficiently record the entire genetic history of a population in forwards-time, individual-based population genetics simulations with arbitrary breeding models, population structure and demography. This approach dramatically reduces the computational burden of tracking individual genomes by allowing us to simulate only those loci that may affect reproduction (those having non-neutral variants). The genetic history of the population is recorded as a succinct tree sequence as introduced in the software package msprime, on which neutral mutations can be quickly placed afterwards. Recording the results of each breeding event requires storage that grows linearly with time, but there is a great deal of redundancy in this information. We solve this storage problem by providing an algorithm to quickly ‘simplify’ a tree sequence by removing this irrelevant history for a given set of genomes. By periodically simplifying the history with respect to the extant population, we show that the total storage space required is modest and overall large efficiency gains can be made over classical forward-time simulations. We implement a general-purpose framework for recording and simplifying genealogical data, which can be used to make simulations of any population model more efficient. We modify two popular forwards-time simulation frameworks to use this new approach and observe efficiency gains in large, whole-genome simulations of one to two orders of magnitude. In addition to speed, our method for recording pedigrees has several advantages: (1) All marginal genealogies of the simulated individuals are recorded, rather than just genotypes. (2) A population of N individuals with M polymorphic sites can be stored in O(N log N + M) space, making it feasible to store a simulation’s entire final generation as well as its history. (3) A simulation can easily be initialized with a more efficient coalescent simulation of deep history. The software for recording and processing tree sequences is named tskit.

Species across a wide range of taxa and habitats are shifting phenological events in response to climate change. While advances are common, shifts vary in magnitude and direction within and among species, and the basis for this variation is relatively unknown. We examine previously suggested patterns of variation in phenological shifts in order to understand the cue–response mechanisms that underlie phenological change. Here, we review what is known about the mechanistic basis for nine factors proposed to predict phenological change (latitude, elevation, habitat type, trophic level, migratory strategy, ecological specialization, species’ seasonality, thermoregulatory mode, and generation time). We find that many studies either do not identify a specific underlying mechanism or do not evaluate alternative mechanistic hypotheses, limiting the ability of scientists to predict future responses to global change with accuracy. We present a conceptual framework that emphasizes a critical distinction between environmental (cue‐driven) and organismal (response‐driven) mechanisms causing variation in phenological shifts and discuss how this distinction can reduce confusion in the field and improve predictions of future phenological change.

Phenotypic plasticity and its evolution may help evolutionary rescue in a novel and stressful environment, especially if environmental novelty reveals cryptic genetic variation that enables the evolution of increased plasticity. However, the environmental stochasticity ubiquitous in natural systems may alter these predictions, because high plasticity may amplify phenotype-environment mismatches. Although previous studies have highlighted this potential detrimental effect of plasticity in stochastic environments, they have not investigated how it affects extinction risk in the context of evolutionary rescue and with evolving plasticity. We investigate this question here by integrating stochastic demography with quantitative genetic theory in a model with simultaneous change in the mean and predictability (temporal autocorrelation) of the environment. We develop an approximate prediction of long-term persistence under the new pattern of environmental fluctuations, and compare it with numerical simulations for short-and long-term extinction risk. We find that reduced predictability increases extinction risk and reduces persistence because it increases stochastic load during rescue. This understanding of how stochastic demography, phenotypic plasticity, and evolution interact when evolution acts on cryptic genetic variation revealed in a novel environment can inform expectations for invasions, extinctions, or the emergence of chemical resistance in pests.

Many natural populations exhibit temporal fluctuations in abundance that are consistent with external forcing by a randomly changing environment. As fitness emerges from an interaction between the phenotype and the environment, such demographic fluctuations probably include a substantial contribution from fluctuating phenotypic selection. We study the stochastic population dynamics of a population exposed to random (plus possibly directional) changes in the optimum phenotype for a quantitative trait that evolves in response to this moving optimum.We derive simple analytical predictions for the distribution of log-population size over time, both transiently and at stationarity under Gompertz density regulation. These predictions are well matched by population-and individual-based simulations. The log-population size is approximately reverse gamma distributed, with a negative skew causing an excess of low relative to high population sizes, thus increasing extinction risk relative to a symmetric (e.g. normal) distribution with the same mean and variance. Our analysis reveals how the mean and variance of log population size change with the variance and autocorrelation of deviations of the evolving mean phenotype from the optimum. We apply our results to the analysis of evolutionary rescue in a stochastic environment, and show that random fluctuations in the optimum can substantially increase extinction risk, by both reducing the expected growth rate and increasing the variance of population size by several orders of magnitude.