From stochastic environments to life histories and back

Environmental stochasticity is known to play an important role in life-history evolution, but most general theory assumes a constant environment. In this paper, we examine life-history evolution in a variable environment, by decomposing average individual fitness (measured by the long-run stochastic growth rate) into contributions from average vital rates and their temporal variation. We examine how generation time, demographic dispersion (measured by the dispersion of reproductive events across the lifespan), demographic resilience (measured by damping time), within-year variances in vital rates, within-year correlations between vital rates and between-year correlations in vital rates combine to determine average individual fitness of stylized life histories. In a fluctuating environment, we show that there is often a range of cohort generation times at which the fitness is at a maximum. Thus, we expect ‘optimal’ phenotypes in fluctuating environments to differ from optimal phenotypes in constant environments. We show that stochastic growth rates are strongly affected by demographic dispersion, even when deterministic growth rates are not, and that demographic dispersion also determines the response of life-history-specific average fitness to within- and between-year correlations. Serial correlations can have a strong effect on fitness, and, depending on the structure of the life history, may act to increase or decrease fitness. The approach we outline takes a useful first step in developing general life-history theory for non-constant environments.     

Accounting for stochasticity in demographic compensation along the elevational range of an alpine plant

Demographic compensation arises when vital rates change in opposite directions across populations, buffering the variation in population growth rates, and is a mechanism often invoked to explain the stability of species geographic ranges. However, studies on demographic compensation have disregarded the effects of temporal variation in vital rates and their temporal correlations, despite theoretical evidence that stochastic dynamics can affect population persistence in temporally varying environments. We carried out a seven‐year‐long demographic study on the perennial plant Arabis alpina (L.) across six populations encompassing most of its elevational range. We discovered demographic compensation in the form of negative correlations between the means of plant vital rates, but also between their temporal coefficients of variation, correlations and elasticities. Even if their contribution to demographic compensation was small, this highlights a previously overlooked, but potentially important, role of stochastic processes in stabilising population dynamics at range margins.     

Stochastic LTRE analysis of the effects of herbivory on the population dynamics of a perennial grassland herb

Herbivores can have strong deleterious effects on vital rates (growth, reproduction, and survival) and thus negatively impact the population dynamics of plant species. In practice, however, these effects might be strongly correlated, for example as a result of tradeoffs between vital rates. To get better insights into the effects of herbivory on the population dynamics of the long‐lived grassland plant Primula veris population projection matrices were constructed from demographic data collected between 1999 and 2008 (nine annual transitions). Data were collected in two large grassland populations, each of which was subjected to two treatments (grazing by cattle versus a mowing treatment), yielding a total of 36 matrices. We applied a lower‐level vital rate life table response experiment (LTRE) using the small noise approximation (SNA) of the stochastic population growth rate to disentangle the contributions of changes in mean vital rates, variability in vital rates, correlations between vital rates and vital rate elasticities to the difference in the stochastic growth rate. Stochastic growth rates (a= log λ_S) were significantly lower in grazed than in mown plots (a= 0.0185 and 0.1019, respectively). SNA LTRE analysis showed that contributions of mean vital rates by far made the largest contribution to the observed difference in a between grazed and control plots. In particular, changes in sexual reproduction rates made the largest contributions to lower the stochastic growth rate in grazed plots: both adult flowering probabilities and flower and seed production were importantly lower in grazed populations, but these negative effects were largely buffered by increased establishment and seedling survival rates. Among the stochastic terms of the SNA decomposition, contributions of covariance and correlations between vital rates had the largest impact, whereas contributions of elasticities were smaller. The strongest correlation driver was the association between adult survival and seedling establishment, suggesting that environmental conditions favouring adult survival also are beneficial for seedling establishment. Overall, our results show that herbivory had a strong negative effect on the long‐term population growth rate of P. veris that was primarily mediated by differences in fecundity (flower and seed production) and germination.     

Evaluating the demographic buffering hypothesis with vital rates estimated for Weddell seals from 30 years of mark-recapture data

1. Life-history theory predicts that those vital rates that make larger contributions to population growth rate ought to be more strongly buffered against environmental variability than are those that are less important. Despite the importance of the theory for predicting demographic responses to changes in the environment, it is not yet known how pervasive demographic buffering is in animal populations because the validity of most existing studies has been called into question because of methodological deficiencies. 2. We tested for demographic buffering in the southern-most breeding mammal population in the world using data collected from 5558 known-age female Weddell seals over 30 years. We first estimated all vital rates simultaneously with mark-recapture analysis and then estimated process variance and covariance in those rates using a hierarchical Bayesian approach. We next calculated the population growth rate's sensitivity to changes in each of the vital rates and tested for evidence of demographic buffering by comparing properly scaled values of sensitivity and process variance in vital rates. 3. We found evidence of positive process covariance between vital rates, which indicates that all vital rates are affected in the same direction by changes in annual environment. Despite the positive correlations, we found strong evidence that demographic buffering occurred through reductions in variation in the vital rates to which population growth rate was most sensitive. Process variation in vital rates was inversely related to sensitivity measures such that variation was greatest in breeding probabilities, intermediate for survival rates of young animals and lowest for survival rates of older animals. 4. Our work contributes to a small but growing set of studies that have used rigorous methods on long-term, detailed data to investigate demographic responses to environmental variation. The information from these studies improves our understanding of life-history evolution in stochastic environments and provides useful information for predicting population responses to future environmental change. Our results for an Antarctic apex predator also provide useful baselines from a marine ecosystem when its top- and middle-trophic levels were not substantially impacted by human activity.     

Does environmental stochasticity matter? Analysis of red deer life-histories on Rum

Summary Most life-history theory assumes that short-term variation in an organism's environment does not affect the survivorships and fecundities of the organisms. This assumption is rarely met. Here we investigate the population and evolutionary biology of red deer,Cervus elephas, to see if relaxation of this assumption is likely to make significant differences to the predicted evolutionary biology of this species. To do this we used 21 years of data from a population of deer on Rum, Western Isles, Scotland. Population growth rates in a stochastic environment were estimated using Tuljapurkar's small noise approximation, confirmed by bootstrap simulation. Numerical differentiation was used to see if the selection pressures (i.e. sensitivities of population growth rate to changes in the vital rates) differ between the stochastic and deterministic cases. The data also allow the costs of reproduction to be estimated. These costs, incorporated as trade-offs into the sensitivity analysis, allow investigation of evolutionary benefits of different life-history tactics. Environmentally induced stochastic variation in the red deer vital rates causes a slight reduction ( 1%) in the predicted population growth rate and has little impact on the estimated selection pressures on the deer's life-history. We thus conclude that, even though density-independent stochastic effects on the population are marked, the deer's fitness is not markedly affected by these and they are adapted to the average conditions they experience. However, the selected life-history is sensitive to the trade-offs between current fecundity, survivorship and future fecundity and it is likely that the environmental variance will affect these trade-offs and, thus, affect the life-history favoured by selection. We also show that the current average life-history is non-optimal and suggest this is a result of selection pressures exerted by culling and predation, now much reduced. As the use of stochastic or deterministic methods provide similar estimates in this case, the use of the latter is justified. Thus,r (the annual per capita rate of population growth) is an appropriate measure of fitness in a population with stochastic numerical fluctuations. In a population of constant size lifetime reproductive success is the obvious measure of fitness to use.     

Selection against demographic stochasticity in age-structured populations

It has been shown that differences in fecundity variance can influence the probability of invasion of a genotype in a population; i.e., a genotype with lower variance in offspring number can be favored in finite populations even if it has a somewhat lower mean fitness than a competitor. In this article, Gillespie's results are extended to population genetic systems with explicit age structure, where the demographic variance (variance in growth rate) calculated in the work of Engen and colleagues is used as a generalization of "variance in offspring number" to predict the interaction between deterministic and random forces driving change in allele frequency. By calculating the variance from the life-history parameters, it is shown that selection against variance in the growth rate will favor a genotypes with lower stochasticity in age-specific survival and fertility rates. A diffusion approximation for selection and drift in a population with two genotypes with different life-history matrices (and therefore different mean growth rates and demographic variances) is derived and shown to be consistent with individual-based simulations. It is also argued that for finite populations, perturbation analyses of both the mean and the variance in growth rate may be necessary to determine the sensitivity of fitness to changes in the life-history parameters.     

Selection and genetic (co)variance in bighorn sheep

Genetic theory predicts that directional selection should deplete additive genetic variance for traits closely related to fitness, and may favor the maintenance of alleles with antagonistically pleiotropic effects on fitness-related traits. Trait heritability is therefore expected to decline with the degree of association with fitness, and some genetic correlations between selected traits are expected to be negative. Here we demonstrate a negative relationship between trait heritability and association with lifetime reproductive success in a wild population of bighorn sheep (Ovis canadensis) at Ram Mountain, Alberta, Canada. Lower heritability for fitness-related traits, however, was not wholly a consequence of declining genetic variance, because those traits showed high levels of residual variance. Genetic correlations estimated between pairs of traits with significant heritability were positive. Principal component analyses suggest that positive relationships between morphometric traits constitute the main axis of genetic variation. Trade-offs in the form of negative genetic or phenotypic correlations among the traits we have measured do not appear to constrain the potential for evolution in this population.     

Differential effects of herbivory and pathogen infestation on plant population dynamics

Plants often interact with antagonists such as herbivores or pathogens. Negative effects on individual plant performance are widely documented, but less is known about whether such effects translate into effects on population viability. In temperate forests, important herbivores include deer. During 2006–2009, I compared vital rates and population growth rates (calculated using integral projection models) between fenced exclosures and grazed control areas, using the perennial herb Phyteuma spicatum as a model species. Deer caused the largest damage to flowering individuals, removing about 24% of all inflorescences and 13% of the above-ground biomass. Only few vital rates seemed to be negatively affected by deer (mainly seed production) and this did not translate into effects on population growth rate. Contrary to expectations, population growth rates tended to be lower in the fenced exclosures in 1 year. This was likely caused by high-pathogen infestation rates, which negatively affected the probability of adult survival and growth. Population growth rate was more sensitive to changes in these vital rates than to changes in seed production. In summary, the results of this demographic study show that grazing effects may be small for long-lived herbs, and that negative effects on vital rates such as seed production may not always translate into effects on population growth rate. The findings also illustrate that other antagonists such as pathogens may be of greater relative importance for differences in population performance than herbivores.     

The contribution of the mitochondrial genome to sex‐specific fitness variance

Maternal inheritance of mitochondrial DNA (mtDNA) facilitates the evolutionary accumulation of mutations with sex‐biased fitness effects. Whereas maternal inheritance closely aligns mtDNA evolution with natural selection in females, it makes it indifferent to evolutionary changes that exclusively benefit males. The constrained response of mtDNA to selection in males can lead to asymmetries in the relative contributions of mitochondrial genes to female versus male fitness variation. Here, we examine the impact of genetic drift and the distribution of fitness effects (DFE) among mutations—including the correlation of mutant fitness effects between the sexes—on mitochondrial genetic variation for fitness. We show how drift, genetic correlations, and skewness of the DFE determine the relative contributions of mitochondrial genes to male versus female fitness variance. When mutant fitness effects are weakly correlated between the sexes, and the effective population size is large, mitochondrial genes should contribute much more to male than to female fitness variance. In contrast, high fitness correlations and small population sizes tend to equalize the contributions of mitochondrial genes to female versus male variance. We discuss implications of these results for the evolution of mitochondrial genome diversity and the genetic architecture of female and male fitness.     

Fitness consequences of maternal and grandmaternal effects

Transgenerational effects are broader than only parental relationships. Despite mounting evidence that multigenerational effects alter phenotypic and life‐history traits, our understanding of how they combine to determine fitness is not well developed because of the added complexity necessary to study them. Here, we derive a quantitative genetic model of adaptation to an extraordinary new environment by an additive genetic component, phenotypic plasticity, maternal and grandmaternal effects. We show how, at equilibrium, negative maternal and negative grandmaternal effects maximize expected population mean fitness. We define negative transgenerational effects as those that have a negative effect on trait expression in the subsequent generation, that is, they slow, or potentially reverse, the expected evolutionary dynamic. When maternal effects are positive, negative grandmaternal effects are preferred. As expected under Mendelian inheritance, the grandmaternal effects have a lower impact on fitness than the maternal effects, but this dual inheritance model predicts a more complex relationship between maternal and grandmaternal effects to constrain phenotypic variance and so maximize expected population mean fitness in the offspring.     

Analysis of interspecific competition in perennial plants using Life Table Response Experiments

The impact of interspecific competition is usually measured by its effect upon plant growth, neglecting impacts upon other stages of the life cycle such as fecundity which have a direct influence upon individual fitness and the asymptotic population growth rate of a population (λ). We used parameterized matrix models for three perennial plant species grown with and without interspecific competition to illustrate how the methodology of Life Table Response Experiments (LTRE) can be used to link any change in population dynamics to changes in any part of the life cycle. Plants were herbaceous grassland species grown for two years in a field experiment at Rothamsted Experimental Station, England. Interspecific competition reduced λ by over 90% in all species. Survival and growth were slightly affected by competition whereas plant fecundity was greatly reduced. Nearly all of the observed difference in λ between the competition treatments was explained by the fecundity terms, and more precisely by a large difference in the number of seeds, and a high sensitivity of λ to the germination rate. Whereas most competition studies focus on the measurement of change in individual fitness, our study illustrates how informative it is to take account not only of the effect of competition upon vital rates but also of how different vital rates affect population growth rate.     

Variability of important vital rates challenges the demographic buffering hypothesis

Selection is assumed to eliminate life-histories showing high variability in vital rates that have the greatest influence on population performance. Therefore, an inverse variability-importance relationship of vital rates is believed to be a universal pattern for diverse life-histories. We tested for such a relationship using multi-year demographic data on a large number of populations of two perennial plant species. Applying different approaches, we first examined the overall variability-importance relationship for the average main vital rates (survival, growth, retrogression, fecundity) per species, and then separately for each population. We found an overall inverse relationship between temporal variation and importance of the average main vital rates for both study species, but these negative species-level correlations were mainly caused by different scales of the examined vital rates. When variability-importance relationships were examined across individual demographic transitions within populations, the abundance of positive and negative correlations depended largely on the method used, and positive correlations were more common after correcting vital rates for sampling variation than when using uncorrected vital rates. Our results cast doubt on the generality of the demographic buffering hypothesis, suggesting that the inverse variability-importance relationship may not be a universal pattern when vital rates are examined for multiple populations of the same plant species.     

Buffering and plasticity in vital rates of oldfield rodents

1. Under the hypothesis of environmental buffering, populations are expected to minimize the variance of the most influential vital rates; however, this may not be a universal principle. Species with a life span <1 year may be less likely to exhibit buffering because of temporal or seasonal variability in vital rate sensitivities. Further, plasticity in vital rates may be adaptive for species in a variable environment with reliable cues. 2. We tested for environmental buffering and plasticity in vital rates using stage-structured matrix models from long-term data sets in four species of grassland rodents. We used periodic matrices to estimate stochastic elasticity for each vital rate and then tested for correlations with a standardized coefficient of variation for each rate. 3. We calculated stochastic elasticities for individual months to test for an association between increased reproduction and the influence of reproduction, relative to survival, on the population growth rate. 4. All species showed some evidence of buffering. The elasticity of vital rates of Peromyscus leucopus (Rafinesque, 1818), Sigmodon hispidus Say & Ord, 1825 and Microtus ochrogaster (Wagner, 1842) was negatively related to vital rate CV. Elasticity and vital rate CV were negatively related in Peromyscus maniculatus (Wagner, 1845), but the relationship was not statistically significant. Peromyscus leucopus and M. ochrogaster showed plasticity in vital rates; reproduction was higher following months where elasticity for reproduction exceeded that of survival. 5. Our results suggest that buffering is common in species with fast life histories; however, some populations that exhibit buffering are capable of responding to short-term variability in environmental conditions through reproductive plasticity.     

Scaling up phenotypic plasticity with hierarchical population models

Individuals respond to different environments by developing different phenotypes, which is generally seen as a mechanism through which individuals can buffer adverse environmental conditions and increase their fitness. To understand the consequences of phenotypic plasticity it is necessary to study how changing a particular trait of an individual affects either its survival, growth, reproduction or a combination of these demographic vital rates (i.e. fitness components). Integrating vital rate changes due to phenotypic plasticity into models of population dynamics allows detailed study of how phenotypic changes scale up to higher levels of integration and forms an excellent tool to distinguish those plastic trait changes that really matter at the population level. A modeling approach also facilitates studying systems that are even more complex: traits and vital rates often co-vary or trade-off with other traits that may show plastic responses over environmental gradients. Here we review recent developments in the literature on population models that attempt to include phenotypic plasticity with a range of evolutionary assumptions and modeling techniques. We present in detail a model framework in which environmental impacts on population dynamics can be followed analytically through direct and indirect pathways that importantly incorporate phenotypic plasticity, trait-trait and trait-vital rate relationships. We illustrate this framework with two case studies: the population-level consequences of phenotypic responses to nutrient enrichment of plant species occurring in nutrient-poor habitats and of responses to changes in flooding regimes due to climate change. We conclude with exciting prospects for further development of this framework: selection analyses, modeling advances and the inclusion of spatial dynamics by considering dispersal traits as well.     

Buffering of life histories against environmental stochasticity: accounting for a spurious correlation between the variabilities of vital rates and their contributions to fitness

Life-history theory predicts vital rates that on average make large contributions to the annual multiplication rate of a lineage should be highly buffered against environmental variability. This prediction has been tested by looking for a negative correlation between the sensitivities (or elasticities) of the elements in a projection matrix and their variances (or coefficients of variation). Here, we show by constructing random matrices that a spurious negative correlation exists between the sensitivities and variances, and between the elasticities and coefficients of variation, of matrix elements. This spurious correlation arises in part because size transition probabilities, which are bounded by 0 and 1, have a limit to their variability that often does not apply to matrix elements representing reproduction. We advocate an alternative analysis based on the underlying vital rates (not the matrix elements) that accounts for the inherent limit to the variability of zero-to-one vital rates, corrects for sampling variation, and tests for a declining upper limit to variability as a vital rate's fitness contribution increases. Applying this analysis to demographic data from five populations of the alpine cushion plant Silene acaulis, we provide evidence of stronger buffering in the vital rates that most influence fitness.     

Heterozygosity-fitness correlations in the gilthead sea bream Sparus aurata using microsatellite loci from unknown and gene-rich genomic locations

Heterozygosity-fitness correlations (HFC) were assessed for a sample of a gilthead sea bream Sparus aurata population. Two hundred and seventy-one fish were genotyped at 22 known and novel microsatellite loci, from which correlations between the multilocus heterozygosity index (I(MLH) ) and various fitness traits (fork length, mass and specific growth rates) were calculated. Significant global HFCs were found in this sample (0·02 ≤r(2) ≤ 0·08). In addition, all the significant correlations found in this work were negative, indicating that heterozygotes had lower fitness than their homozygote counterparts. Marker location could not explain the observed HFCs. Evidence of inbreeding, outbreeding or population and family structuring was not found in this work. The presence of undetected general effects that may lead to the appearance of HFCs, however, cannot be ruled out. These results seem to be best explained by the occurrence of local effects (due to linkage) or even by possible direct locus advantages.     

Extinction in relation to demographic and environmental stochasticity in age-structured models

The demographic variance of an age-structured population is defined. This parameter is further split into components generated by demographic stochasticity in each vital rate. The applicability of these parameters are investigated by checking how an age-structured population process can be approximated by a diffusion with only three parameters. These are the deterministic growth rate computed from the expected projection matrix and the environmental and demographic variances. We also consider age-structured populations where the fecundity at any stage is either zero or one, and there is neither environmental stochasticity nor dependence between individual fecundity and survival. In this case the demographic variance is uniquely determined by the vital rates defining the projection matrix. The demographic variance for a long-lived bird species, the wandering albatross in the southwestern part of the Indian Ocean, is estimated. We also compute estimates of the age-specific contributions to the total demographic variance from survival, fecundity and the covariance between survival and fecundity.     

The devil in the details of life-history evolution: Instability and reversal of genetic correlations during selection on<Emphasis Type="Italic">Drosophila</Emphasis> development

The evolutionary relationships between three major components of Darwinian fitness, development rate, growth rate and preadult survival, were estimated using a comparison of 55 distinct populations ofDrosophila melanogaster variously selected for age-specific fertility, environmental-stress tolerance and accelerated development. Development rate displayed a strong net negative evolutionary correlation with weight at eclosion across all selection treatments, consistent with the existence of a size-versus-time tradeoff between these characters. However, within the data set, the magnitude of the evolutionary correlation depended upon the particular selection treatments contrasted. A previously proposed tradeoff between preadult viability and growth rate was apparent only under weak selection for juvenile fitness components. Direct selection for rapid development led to sharp reductions in both growth rates and viability. These data add to the mounting results from experimental evolution that illustrate the sensitivity of evolutionary correlations to (i) genotype-by-environment (G X E) interaction, (ii) complex functional-trait interactions, and (iii) character definition. Instability, disappearance and reversal of patterns of genetic covariation often occur over short evolutionary time frames and as the direct product of selection, rather than some stochastic process. We suggest that the functional architecture of fitness is a rapidly evolving matrix with reticulate properties, a matrix that we understand only poorly.     

Population Dynamics of Breeding Mallards in the Great Lakes States

Abstract: Mallard (Anas platyrhynchos) populations in the United States portion of the Great Lakes region increased through the 1990s but have since declined. To promote sustainable growth of this population, managers need to understand how perturbation of vital rates will affect annual population growth rate (Λ). We developed a stage-based model representing the female mallard population in the Great Lakes using vital rates generated from a landscape-level study documenting reproductive parameters from 2001 to 2003. We conducted perturbation analyses (i.e., sensitivity analyses) to identify vital rates that most influence Λ and variance decomposition analyses to determine the proportion of variation in Λ explained by variation in each vital rate. Perturbation analyses indicated that Λ was most sensitive to changes in nonbreeding survival, duckling survival, and nest success. Therefore, changes in these vital rates would be expected to result in the greatest ΔΛ. Process variation in breeding season parameters accounted for 63% of variation in Λ. Breeding season parameters explaining the most variation were duckling survival (32%) and nest success (16%). Survival of adult females outside the breeding season accounted for 36% of variation in Λ. Harvest derivation, high harvest, and high sensitivity of Λ to nonbreeding survival for Great Lakes female mallards suggests there is a strong potential for managing the Great Lakes mallard population via harvest management. Because Λ was highly sensitive to changes in duckling survival, we suggest programs that emphasize wetland protection, enhancement, and restoration as a management strategy to improve population growth for breeding mallards.     

Environment dependence of mutational parameters for viability in Drosophila melanogaster

The genomic rate of mildly deleterious mutations (U) figures prominently in much evolutionary and ecological theory. In Drosophila melanogaster, estimates of U have varied widely, from <0.1 to nearly 1 per zygote. The source of this variation is unknown, but could include differences in the conditions used for assaying fitness traits. We examined how assay conditions affect estimates of the rates and effects of viability-depressing mutations in two sets of lines with accumulated spontaneous mutations on the second chromosome. In each set, the among-line variance in egg-to-adult viability was significantly greater when viability was assayed using a high parental density than when it was assayed using a low density. In contrast, the proportional decline in viability due to new mutations did not differ between densities. Two other manipulations, lowering the temperature and adding ethanol to the medium, had no significant effects on either the mean decline or among-line variance. Cross-environment genetic correlations in viability were generally close to one, implying that most mutations reduced viability in all environments. Using data from the low-density, lower-bound estimates of U approached the classic, high values of Mukai and Ohnishi; at the high density, U estimates were similar to recently reported low values. The difference in estimated mutation rates, taken at face value, would imply that many mutations affected fitness at low density but not at high density, but this is shown to be incompatible with the observed high cross-environment correlations. Possible reasons for this discrepancy are discussed. Regardless of the interpretation, the results show that assay conditions can have a large effect on estimates of mutational parameters for fitness traits.