ORGANISMAL COMPLEXITY AND THE POTENTIAL FOR EVOLUTIONARY DIVERSIFICATION

We present two theoretical approaches to investigate whether organismal complexity, defined as the number of quantitative traits determining fitness, and the potential for adaptive diversification are correlated. The first approach is independent of any specific ecological model and based on curvature properties of the fitness landscape as a function of the dimension of the trait space. This approach indeed suggests a positive correlation between complexity and diversity. An assumption made in this first approach is that the potential for any pair of traits to interact in their effect on fitness is independent of the dimension of the trait space. In the second approach, we circumvent making this assumption by analyzing the evolutionary dynamics in an explicit consumer‐resource model in which the shape of the fitness landscape emerges from the underlying mechanistic ecological model. In this model, consumers are characterized by several quantitative traits and feed on a multidimensional resource distribution. The consumer's feeding efficiency on the resource is determined by the match between consumer phenotype and resource item. This analysis supports a positive correlation between the complexity of the evolving consumer species and its potential to diversify with the additional insight that also increasing resource complexity facilitates diversification.     

The interplay between behavior and morphology in the evolutionary dynamics of resource specialization

We analyze the consequences of diet choice behavior for the evolutionary dynamics of foraging traits by means of a mathematical model. The model is characterized by the following features. Consumers feed on two different substitutable resources that are distributed in a fine-grained manner. On encounter with a resource item, consumers decide whether to attack it so as to maximize their energy intake. Simultaneously, evolutionary change occurs in morphological traits involved in the foraging process. The assumption here is that evolution is constrained by a trade-off in the consumer's ability to forage on the alternative resources. The model predicts that flexible diet choice behavior can guide the direction of evolutionary change and mediate coexistence of different consumer types. Such polymorphisms can evolve from a monomorphic population at evolutionary branching points and also at points where a small genetic change in a trait can provoke a sharp instantaneous and nongenetic change in choice behavior. In the case of weak trade-offs, the evolutionary dynamics of a dimorphic consumer population can lead to alternative evolutionarily stable communities. The robustness of these predictions is checked with individual-based simulations and by relaxing the assumption of optimally foraging consumers.     

An evolutionary quantitative genetics model for phenotypic (co)variances under limited dispersal,with an application to socially synergistic traits

Darwinian evolution consists of the gradual transformation of heritable traits due to natural selection and the input of random variation by mutation. Here, we use a quantitative genetics approach to investigate the coevolution of multiple quantitative traits under selection, mutation, and limited dispersal. We track the dynamics of trait means and of variance–covariances between traits that experience frequency‐dependent selection. Assuming a multivariate‐normal trait distribution, we recover classical dynamics of quantitative genetics, as well as stability and evolutionary branching conditions of invasion analyses, except that due to limited dispersal, selection depends on indirect fitness effects and relatedness. In particular, correlational selection that associates different traits within‐individuals depends on the fitness effects of such associations between‐individuals. We find that these kin selection effects can be as relevant as pleiotropy for the evolution of correlation between traits. We illustrate this with an example of the coevolution of two social traits whose association within‐individuals is costly but synergistically beneficial between‐individuals. As dispersal becomes limited and relatedness increases, associations between‐traits between‐individuals become increasingly targeted by correlational selection. Consequently, the trait distribution goes from being bimodal with a negative correlation under panmixia to unimodal with a positive correlation under limited dispersal.     

Adaptive change in the resource-exploitation traits of a generalist consumer: the evolution and coexistence of generalists and specialists

Mathematical models of consumer-resource systems are used to explore the evolution of traits related to resource acquisition in a generalist consumer species that is capable of exploiting two resources. The analysis focuses on whether evolution of traits determining the capture rates of two resources by a consumer species produce one generalist, two specialists, or all three types, when all types are characterized by a common fitness function. In systems with a stable equilibrium, evolution produces one generalist or two specialists, depending on the second derivative of the trade-off relationship. When there are sustained population fluctuations, the nature of the trade-off between the consumer's capture rates of the two resources still plays a key role in determining the evolutionary outcome. If the trade-off is described by a choice variable between zero and one that is raised to a power n, polymorphic states are possible when n > 1, which implies a positive second derivative of the curve. These states are either dimorphism, with two relatively specialized consumer types, or trimorphism, with a single generalist type and two specialists. Both endogenously driven consumer-resource cycles, and fluctuations driven by an environmental variable affecting resource growth are considered. Trimorphic evolutionary outcomes are relatively common in the case of endogenous cycles. In contrast to a previous study, these trimorphisms can often evolve even when new lineages are constrained to have phenotypes very similar to existing lineages. Exogenous cycles driven by environmental variation in resource growth rates appear to be much less likely to produce a mixture of generalists and specialists than are endogenous consumer-resource cycles.     

Individual specialization in a shorebird population with narrow foraging niche

Individual specialization in resource use is a widespread driver for intra-population trait variation, playing a crucial evolutionary role in free-living animals. We investigated the individual foraging specialization of Black-tailed Godwits (Limosa limosa islandica) during the wintering period. Godwits displayed distinct degrees of individual specialization in diet and microhabitat use, indicating the presence of both generalist and specialist birds. Females were overall more specialist than males, primarily consuming polychaetes. Specialist males consumed mainly bivalves, but some individuals also specialized on gastropods or polychaetes. Sexual dimorphism in bill length is probably important in determining the differences in specialization, as longer-billed individuals have access to deep-buried polychaetes inaccessible to most males. Different levels of specialization within the same sex, unrelated to bill length, were also found, suggesting that mechanisms other traits are involved in explaining individual specialization. Godwits specialized on bivalves achieved higher intake rates than non-specialist birds, supporting the idea that individual foraging choices or skills result in different short-term payoffs within the same population. Understanding whether short-term payoffs are good indicators of long-term fitness and how selection operates to favour the prevalence of specialist or generalist godwits is a major future challenge.     

Why climate change will invariably alter selection pressures on phenology

The seasonal timing of lifecycle events is closely linked to individual fitness and hence, maladaptation in phenological traits may impact population dynamics. However, few studies have analysed whether and why climate change will alter selection pressures and hence possibly induce maladaptation in phenology. To fill this gap, we here use a theoretical modelling approach. In our models, the phenologies of consumer and resource are (potentially) environmentally sensitive and depend on two different but correlated environmental variables. Fitness of the consumer depends on the phenological match with the resource. Because we explicitly model the dependence of the phenologies on environmental variables, we can test how differential (heterogeneous) versus equal (homogeneous) rates of change in the environmental variables affect selection on consumer phenology. As expected, under heterogeneous change, phenotypic plasticity is insufficient and thus selection on consumer phenology arises. However, even homogeneous change leads to directional selection on consumer phenology. This is because the consumer reaction norm has historically evolved to be flatter than the resource reaction norm, owing to time lags and imperfect cue reliability. Climate change will therefore lead to increased selection on consumer phenology across a broad range of situations.     

Quantification and decomposition of environment‐selection relationships

In nature, selection varies across time in most environments, but we lack an understanding of how specific ecological changes drive this variation. Ecological factors can alter phenotypic selection coefficients through changes in trait distributions or individual mean fitness, even when the trait‐absolute fitness relationship remains constant. We apply and extend a regression‐based approach in a population of Soay sheep (Ovis aries) and suggest metrics of environment‐selection relationships that can be compared across studies. We then introduce a novel method that constructs an environmentally structured fitness function. This allows calculation of full (as in existing approaches) and partial (acting separately through the absolute fitness function slope, mean fitness, and phenotype distribution) sensitivities of selection to an ecological variable. Both approaches show positive overall effects of density on viability selection of lamb mass. However, the second approach demonstrates that this relationship is largely driven by effects of density on mean fitness, rather than on the trait‐fitness relationship slope. If such mechanisms of environmental dependence of selection are common, this could have important implications regarding the frequency of fluctuating selection, and how previous selection inferences relate to longer term evolutionary dynamics.     

Evolutionary stability of optimal foraging: Partial preferences in the diet and patch models

In this article the patch and diet choice models of the optimal foraging theory are reanalyzed with respect to evolutionary stability of the optimal foraging strategies. In their original setting these fundamental models consider a single consumer only and the resulting fitness functions are both frequency and density independent. Such fitness functions do not allow us to apply the classical game theoretical methods to study an evolutionary stability of optimal foraging strategies for competing animals. In this article frequency and density dependent fitness functions of optimal foraging are derived by separation of time scales in an underlying population dynamical model and corresponding evolutionarily stable strategies are calculated. Contrary to the classical foraging models the results of the present article predict that partial preferences occur in optimal foraging strategies as a consequence of the ecological feedback of consumer preferences on consumer fitness. In the case of the patch occupation model these partial preferences correspond to the ideal free distribution concept while in the case of the diet choice model they correspond to the partial inclusion of the less profitable prey type in predators diet.     

The adaptive dynamics of niche constructing traits in spatially subdivided populations: evolving posthumous extended phenotypes

Niche construction, by which organisms modify the environment in which they live, has been proposed to affect the evolution of many phenotypic traits. But what about the evolution of a niche constructing trait itself, whose expression changes the pattern of natural selection to which the trait is exposed in subsequent generations? This article provides an inclusive fitness analysis of selection on niche constructing phenotypes, which can affect their environment from local to global scales in arbitrarily spatially subdivided populations. The model shows that phenotypic effects of genes extending far beyond the life span of the actor can be affected by natural selection, provided they modify the fitness of those individuals living in the future that are likely to have inherited the niche construction lineage of the actor. Present benefits of behaviors are thus traded off against future indirect costs. The future costs will generally result from a complicated interplay of phenotypic effects, population demography and environmental dynamics. To illustrate these points, I derive the adaptive dynamics of a trait involved in the consumption of an abiotic resource, where resource abundance in future generations feeds back to the evolutionary dynamics of the trait.     

Evolutionary potential of morphological traits across different life‐history stages

Despite accumulating examples of selection acting on heritable traits in the wild, predicted evolutionary responses are often different from observed phenotypic trends. Various explanations have been suggested for these mismatches. These include within‐individual changes across lifespan that can create important variation in genetic architecture of traits and selection acting on them, but also potential problems with the methodological approach used to predict evolutionary responses of traits. Here, we used an 8‐year data set on tree swallow (Tachycineta bicolor) to first assess the effects of differences among three nestling life‐history stages on the genetic (co)variances of two morphological traits (body mass and primary feather length) and the selection acting on them over three generations. We then estimated the evolutionary potential of these traits by predicting their evolutionary responses using the breeder's equation and the secondary theorem of selection approaches. Our results showed variation in strength and direction of selection and slight changes in trait variance across ages. Predicted evolutionary responses differed importantly between both approaches for half of the trait–age combinations we studied, suggesting the presence of environmentally induced correlations between focal traits and fitness possibly biasing breeder's equation predictions. Our results emphasize that predictions of evolutionary potential for morphological traits are likely to be highly variable, both in strength and direction, depending on the life stage and method used, thus mitigating our capacity to predict adaptation and persistence of wild populations.     

Gillespie eco‐evolutionary models (GEMs) reveal the role of heritable trait variation in eco‐evolutionary dynamics

Heritable trait variation is a central and necessary ingredient of evolution. Trait variation also directly affects ecological processes, generating a clear link between evolutionary and ecological dynamics. Despite the changes in variation that occur through selection, drift, mutation, and recombination, current eco‐evolutionary models usually fail to track how variation changes through time. Moreover, eco‐evolutionary models assume fitness functions for each trait and each ecological context, which often do not have empirical validation. We introduce a new type of model, Gillespie eco‐evolutionary models (GEMs), that resolves these concerns by tracking distributions of traits through time as eco‐evolutionary dynamics progress. This is done by allowing change to be driven by the direct fitness consequences of model parameters within the context of the underlying ecological model, without having to assume a particular fitness function. GEMs work by adding a trait distribution component to the standard Gillespie algorithm – an approach that models stochastic systems in nature that are typically approximated through ordinary differential equations. We illustrate GEMs with the Rosenzweig–MacArthur consumer–resource model. We show not only how heritable trait variation fuels trait evolution and influences eco‐evolutionary dynamics, but also how the erosion of variation through time may hinder eco‐evolutionary dynamics in the long run. GEMs can be developed for any parameter in any ordinary differential equation model and, furthermore, can enable modeling of multiple interacting traits at the same time. We expect GEMs will open the door to a new direction in eco‐evolutionary and evolutionary modeling by removing long‐standing modeling barriers, simplifying the link between traits, fitness, and dynamics, and expanding eco‐evolutionary treatment of a greater diversity of ecological interactions. These factors make GEMs much more than a modeling advance, but an important conceptual advance that bridges ecology and evolution through the central concept of heritable trait variation.     

The adaptive dynamics of function-valued traits

This study extends the framework of adaptive dynamics to function-valued traits. Such adaptive traits naturally arise in a great variety of settings: variable or heterogeneous environments, age-structured populations, phenotypic plasticity, patterns of growth and form, resource gradients, and in many other areas of evolutionary ecology. Adaptive dynamics theory allows analysing the long-term evolution of such traits under the density-dependent and frequency-dependent selection pressures resulting from feedback between evolving populations and their ecological environment. Starting from individual-based considerations, we derive equations describing the expected dynamics of a function-valued trait in asexually reproducing populations under mutation-limited evolution, thus generalizing the canonical equation of adaptive dynamics to function-valued traits. We explain in detail how to account for various kinds of evolutionary constraints on the adaptive dynamics of function-valued traits. To illustrate the utility of our approach, we present applications to two specific examples that address, respectively, the evolution of metabolic investment strategies along resource gradients, and the evolution of seasonal flowering schedules in temporally varying environments.     

Character convergence under competition for nutritionally essential resources

Resource competition is thought to drive divergence in resource use traits (character displacement) by generating selection favoring individuals able to use resources unavailable to others. However, this picture assumes nutritionally substitutable resources (e.g., different prey species). When species compete for nutritionally essential resources (e.g., different nutrients), theory predicts that selection drives character convergence. We used models of two species competing for two essential resources to address several issues not considered by existing theory. The models incorporated either slow evolutionary change in resource use traits or fast physiological or behavioral change. We report four major results. First, competition always generates character convergence, but differences in resource requirements prevent competitors from evolving identical resource use traits. Second, character convergence promotes coexistence. Competing species always attain resource use traits that allow coexistence, and adaptive trait change stabilizes the ecological equilibrium. In contrast, adaptation in allopatry never preadapts species to coexist in sympatry. Third, feedbacks between ecological dynamics and trait dynamics lead to surprising dynamical trajectories such as transient divergence in resource use traits followed by subsequent convergence. Fourth, under sufficiently slow trait change, ecological dynamics often drive one of the competitors to near extinction, which would prevent realization of long-term character convergence in practice.     

An Evolutionary Dynamics Model Adapted to Eusocial Insects

This study aims to better understand the evolutionary processes allowing species coexistence in eusocial insect communities. We develop a mathematical model that applies adaptive dynamics theory to the evolutionary dynamics of eusocial insects, focusing on the colony as the unit of selection. The model links long-term evolutionary processes to ecological interactions among colonies and seasonal worker production within the colony. Colony population dynamics is defined by both worker production and colony reproduction. Random mutations occur in strategies, and mutant colonies enter the community. The interactions of colonies at the ecological timescale drive the evolution of strategies at the evolutionary timescale by natural selection. This model is used to study two specific traits in ants: worker body size and the degree of collective foraging. For both traits, trade-offs in competitive ability and other fitness components allows to determine conditions in which selection becomes disruptive. Our results illustrate that asymmetric competition underpins diversity in ant communities.     

The prerequisites for and likelihood of generalist-specialist coexistence

Mathematical models of three-consumer-two-resource systems are used to explore the possibility of coexistence when one consumer is a generalist utilizing both resources, and the other two are specialists utilizing only one. Such coexistence requires strongly saturating functional or numerical responses in at least one consumer and the presence of sustained asynchronous variation in resource abundances. Given these conditions, the effects of three dichotomous factors on the range of parameters allowing coexistence are examined: flexible versus inflexible resource choice by the generalist, endogenous or exogenous cause of resource cycles, and location of the two resources in a single habitat versus two habitats. Coexistence of all three species is found to be possible for all combinations of these factors except for inflexible choice in a two-habitat environment. Generalists experience frequency-dependent fitness because, when they are abundant, they synchronize resource cycles and/or reduce their amplitude. When the generalist can adaptively adjust its relative foraging on the two resources, coexistence conditions are broadened considerably, and coexistence commonly occurs readily with exogenous variation in resource growth and with resources located in distinct habitats. Adaptive behavior increases the generalist's ability to both synchronize and dampen resource cycles.     

On the evolution of specialization with a mechanistic underpinning in structured metapopulations

We analyze the evolution of specialization in resource utilization in a discrete-time metapopulation model using the adaptive dynamics approach. The local dynamics in the metapopulation are based on the Beverton-Holt model with mechanistic underpinnings. The consumer faces a trade-off in the abilities to consume two resources that are spatially heterogeneously distributed to patches that are prone to local catastrophes. We explore the factors favoring the spread of generalist or specialist strategies. Increasing fecundity or decreasing catastrophe probability favors the spread of the generalist strategy and increasing environmental heterogeneity enlarges the parameter domain where the evolutionary branching is possible. When there are no catastrophes, increasing emigration diminishes the parameter domain where the evolutionary branching may occur. Otherwise, the effect of emigration on evolutionary dynamics is non-monotonous: both small and large values of emigration probability favor the spread of the specialist strategies whereas the parameter domain where evolutionary branching may occur is largest when the emigration probability has intermediate values. We compare how different forms of spatial heterogeneity and different models of local growth affect the evolutionary dynamics. We show that even small changes in the resource dynamics may have outstanding evolutionary effects to the consumers.     

Context dependence in complex adaptive landscapes: frequency and trait‐dependent selection surfaces within an adaptive radiation of Caribbean pupfishes

The adaptive landscape provides the foundational bridge between micro‐ and macroevolution. One well‐known caveat to this perspective is that fitness surfaces depend on ecological context, including competitor frequency, traits measured, and resource abundance. However, this view is based largely on intraspecific studies. It is still unknown how context‐dependence affects the larger features of peaks and valleys on the landscape which ultimately drive speciation and adaptive radiation. Here, I explore this question using one of the most complex fitness landscapes measured in the wild in a sympatric pupfish radiation endemic to San Salvador Island, Bahamas by tracking survival and growth of laboratory‐reared F2 hybrids. I present new analyses of the effects of competitor frequency, dietary isotopes, and trait subsets on this fitness landscape. Contrary to expectations, decreasing competitor frequency increased survival only among very common phenotypes, whereas less common phenotypes rarely survived despite few competitors, suggesting that performance, not competitor frequency, shapes large‐scale features of the fitness landscape. Dietary isotopes were weakly correlated with phenotype and growth, but did not explain additional survival variation. Nonlinear fitness surfaces varied substantially among trait subsets, revealing one‐, two‐, and three‐peak landscapes, demonstrating the complexity of selection in the wild, even among similar functional traits.     

Natural selection, evolvability and bias due to environmental covariance in the field in an annual plant

Estimates of the form and magnitude of natural selection based on phenotypic relationships between traits and fitness measures can be biased when environmental factors influence both relative fitness and phenotypic trait values. I quantified genetic variances and covariances, and estimated linear and quadratic selection coefficients, for seven traits of an annual plant grown in the field. For replicates of 50 paternal half-sib families, coefficients of selection were calculated both for individual phenotypic values of the traits and for half-sib family mean values. The potential for evolutionary response was supported by significant heritability and phenotypic directional selection for several traits but contradicted by the absence of significant genetic variation for fitness estimates and evidence of bias in phenotypic selection coefficients due to environmental covariance for at least two of the traits analysed. Only studies of a much wider range of organisms and traits will reveal the frequency and extent of such bias.     

Evolutionarily unstable fitness maxima and stable fitness minima of continuous traits

Summary We present models of adaptive change in continuous traits for the following situations: (1) adaptation of a single trait within a single population in which the fitness of a given individual depends on the population's mean trait value as well as its own trait value; (2) adaptation of two (or more) traits within a single population; (3) adaptation in two or more interacting species. We analyse a dynamic model of these adaptive scenarios in which the rate of change of the mean trait value is an increasing function of the fitness gradient (i.e. the rate of increase of individual fitness with the individual's trait value). Such models have been employed in evolutionary game theory and are often appropriate both for the evolution of quantitative genetic traits and for the behavioural adjustment of phenotypically plastic traits. The dynamics of the adaptation of several different ecologically important traits can result in characters that minimize individual fitness and can preclude evolution towards characters that maximize individual fitness. We discuss biological circumstances that are likely to produce such adaptive failures for situations involving foraging, predator avoidance, competition and coevolution. The results argue for greater attention to dynamical stability in models of the evolution of continuous traits.     

COMPARING STRENGTHS OF DIRECTIONAL SELECTION: HOW STRONG IS STRONG?

The fundamental equation in evolutionary quantitative genetics, the Lande equation, describes the response to directional selection as a product of the additive genetic variance and the selection gradient of trait value on relative fitness. Comparisons of both genetic variances and selection gradients across traits or populations require standardization, as both are scale dependent. The Lande equation can be standardized in two ways. Standardizing by the variance of the selected trait yields the response in units of standard deviation as the product of the heritability and the variance-standardized selection gradient. This standardization conflates selection and variation because the phenotypic variance is a function of the genetic variance. Alternatively, one can standardize the Lande equation using the trait mean, yielding the proportional response to selection as the product of the squared coefficient of additive genetic variance and the mean-standardized selection gradient. Mean-standardized selection gradients are particularly useful for summarizing the strength of selection because the mean-standardized gradient for fitness itself is one, a convenient benchmark for strong selection. We review published estimates of directional selection in natural populations using mean-standardized selection gradients. Only 38 published studies provided all the necessary information for calculation of mean-standardized gradients. The median absolute value of multivariate mean-standardized gradients shows that selection is on average 54% as strong as selection on fitness. Correcting for the upward bias introduced by taking absolute values lowers the median to 31%, still very strong selection. Such large estimates clearly cannot be representative of selection on all traits. Some possible sources of overestimation of the strength of selection include confounding environmental and genotypic effects on fitness, the use of fitness components as proxies for fitness, and biases in publication or choice of traits to study.