Skip to main content

Main menu

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
    • Front Matter Portal
    • Journal Club
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Home
Home
  • Log in
  • My Cart

Advanced Search

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
    • Front Matter Portal
    • Journal Club
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
Research Article

Asymmetric ecological conditions favor Red-Queen type of continued evolution over stasis

Jan Martin Nordbotten and Nils C. Stenseth
  1. aDepartment of Mathematics, University of Bergen, N-5020 Bergen, Norway;
  2. bDepartment of Civil and Environmental Engineering, Princeton University, Princeton, NJ 08544;
  3. cCentre for Ecological and Evolutionary Synthesis (CEES), Department of Biosciences, University of Oslo, N-0316 Oslo, Norway

See allHide authors and affiliations

PNAS February 16, 2016 113 (7) 1847-1852; first published February 1, 2016; https://doi.org/10.1073/pnas.1525395113
Jan Martin Nordbotten
aDepartment of Mathematics, University of Bergen, N-5020 Bergen, Norway;
bDepartment of Civil and Environmental Engineering, Princeton University, Princeton, NJ 08544;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Nils C. Stenseth
cCentre for Ecological and Evolutionary Synthesis (CEES), Department of Biosciences, University of Oslo, N-0316 Oslo, Norway
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: n.c.stenseth@ibv.uio.no
  1. Contributed by Nils C. Stenseth, December 23, 2015 (sent for review December 11, 2015; reviewed by Eörs Szathmáry and Richard A. Watson)

  • Article
  • Figures & SI
  • Info & Metrics
  • PDF
Loading

Significance

The work presented here demonstrates analytically for the very first time (to our knowledge) that, within a very general theoretical framework, both Red-Queen type of continuous evolution and evolutionary stasis may be the outcomes of ecological interactions within a multispecies ecological community. Whether or not evolution will cease or continue in an abiotically stable environment (i.e., where there are only biotic forces) has been an unsettled problem within evolutionary biology. Our contribution specifies the ecological conditions for which Red-Queen type of continuous evolution and stasis will result. The new and general eco-evolutionary model provides a profoundly new basis for further theoretical and empirical work within the field of coevolution within multispecies ecological systems.

Abstract

Four decades ago, Leigh Van Valen presented the Red Queen’s hypothesis to account for evolution of species within a multispecies ecological community [Van Valen L (1973) Evol Theory 1(1):1–30]. The overall conclusion of Van Valen’s analysis was that evolution would continue even in the absence of abiotic perturbations. Stenseth and Maynard Smith presented in 1984 [Stenseth NC, Maynard Smith J (1984) Evolution 38(4):870–880] a model for the Red Queen’s hypothesis showing that both Red-Queen type of continuous evolution and stasis could result from a model with biotically driven evolution. However, although that contribution demonstrated that both evolutionary outcomes were possible, it did not identify which ecological conditions would lead to each of these evolutionary outcomes. Here, we provide, using a simple, yet general population-biologically founded eco-evolutionary model, such analytically derived conditions: Stasis will predominantly emerge whenever the ecological system contains only symmetric ecological interactions, whereas both Red-Queen and stasis type of evolution may result if the ecological interactions are asymmetrical, and more likely so with increasing degree of asymmetry in the ecological system (i.e., the more trophic interactions, host–pathogen interactions, and the like there are [i.e., +/− type of ecological interactions as well as asymmetric competitive (−/−) and mutualistic (+/+) ecological interactions]). In the special case of no between-generational genetic variance, our results also predict dynamics within these types of purely ecological systems.

  • coevolution
  • evolution within ecological communities
  • ecosystem structure
  • mathematical modeling
  • mathematical analysis

The major part of any species’ environment is represented by the other (interacting) species in the ecosystem; hence, any evolutionary change made by any species within the ecosystem will, in general, be experienced as an environmental change by the coexisting species. This is the basis for the Red Queen’s hypothesis as presented by Van Valen (1)—a proposition that is very similar to an idea suggested several decades earlier by Fisher (1930) (ref. 2, pp. 44–45) as well as Darwin (3). The Red Queen’s hypothesis continues to attract much attention (3⇓⇓⇓⇓⇓⇓–10). However, within a multispecies ecological system it remains unclear whether evolution will cease or continue in the absence of external abiotic perturbations. Although the Darwinian theory of evolution provides a satisfactory explanation of the mechanisms for evolutionary changes and can be tested against short-term events, it does not in itself provide any predictions regarding large-scale and long-term features of the evolutionary dynamics. Essentially, we do not understand to what extent biotic factors (such as between-species interactions) are important determinants of macroevolution, or rather, what is the relative importance of biotic and abiotic processes in macroevolution (cf. ref. 11). To achieve such an understanding, we need to bring ecological and evolutionary dynamics together into a common model framework.

Here, we address how evolution will proceed within a multispecies ecological system in the absence of external abiotic perturbations. This question was asked by John Maynard Smith† in 1980 in a plenary talk on the rate of extinction (and evolution) at the Second International Congress of Systematic and Evolutionary Biology (ICSEB II).‡ To an audience of several hundred evolutionary biologists, he asked: “How many think evolution will cease if all abiotic perturbations cease?”. The audience was split into two fairly equal parts (see ref. 12, a paper which essentially contains his talk). Our aim is to address this unsolved issue by analyzing a general eco-evolutionary model of a multispecies ecosystem containing all types of ecological interactions (i.e., competitive interactions, trophic interactions, etc.), thereby bringing ecological and evolutionary dynamics together in a common model framework. Much previous work has been done on this topic (see, e.g., refs. 13⇓⇓⇓–17; for a recent review, see ref. 8). Our main methodological contribution is to provide a general analytical framework for addressing this basic question that does not rely on numerical simulations. With this new model framework, we provide new results regarding the conditions for evolutionary stasis and Red-Queen dynamics.

The Model

Definition of the Model.

Define the abundance density function of individuals, n(x,t) (abbreviated simply as n when no confusion can arise), taking values at points x in phenotype space (x∈Ω), and in time t for positive times (t∈ℝ+). The phenotype space is in general multidimensional, representing different phenotypical attributes. Thus, if we consider a model with d traits, then Ω⊂ℝd. The domain Ω (bounded or unbounded) corresponds to admissible phenotype states.

We have chosen to develop a purely deterministic model formulation, of which one important reason is to be able to analyze specific model outcomes (rather than their probability density functions). Our choice is to some extent also influenced by our overall question about whether or not evolution will continue should all abiotic perturbations cease (cf. refs. 18 and 19).

The function n denotes the abundance of individuals, which can be characterized by the phenotype vector x at a given time, t. Because the phenotypes are given as a continuous parameter space, we might think of n as a density measured as [individuals]/[phenotype in phenotype space]. When n(x,t)=0, this corresponds to no individuals of phenotype x. Conversely, when n(x,t) has a local and isolated maximum at x, this intuitively corresponds to a monomorphic species with traits x, whereas having several clustered peaks might be seen as a polymorphic species.

Here, we include no spatial dimension; however, it is straightforward to consider n as a function also of spatial position y∈ℝ3, such that n=n(x,y,t). Thus, let us consider the following nonlocal and nonlinear partial differential equation describing the population growth rate, governing the dynamics of the abundance function, n=n(x,t):∂n∂t=rn−bn2+n∫Ωα(x,x′)n(x′)dx′+∇⋅(g∇n).[1]This equation is valid for n>0 and is constrained such that ∂n/∂t≥0 when n=0. The parameter functions r, b, and g are all functions of x, as described below. Eq. 1 may be seen either as a generalization of an ecological population growth model to include evolution (through the last term), or conversely, as a reaction-diffusion–type Fisher–Kolmogorov equation with a nonlocal interaction (the second-last term) (20). The following processes are incorporated (see Methods for justification):

  • i) Linear birth and mortality rates are given by r(x)n(x,t). The function r(x) may take positive and negative values (r(x)<0); however, for present phenotypes (and species) it cannot be negative for all x.

  • ii) Quadratic self-interaction corresponding to self-limitation and is given by b(x)n(x,t)2.

  • iii) Bilinear interaction between species are given by the integral interaction operator n(x,t)∫Ωα(x,x′)n(x′,t)dx′. The function α(x,x′) represents the impact of individuals with phenotype x′ on individuals with phenotype x. We choose the sign convention such that competition is characterized by negative values of α and symbiotic interaction by positive values.

  • iv) Generational genetic mutations are included in the model as a diffusion process, given by ∇⋅(g(x)∇n(x,t)); the function g(x) represents the genetically based variability in phenotype between generations, and as usual ∇⋅ is the divergence operator, whereas ∇ is the gradient, both with respect to x. If g=0, the model will correspond to an ecological model as no new genotypes (and resulting phenotypes) will occur. [Notice that g might be seen as a measure of evolvability (cf. refs. 21 and 22).]

We consider the processes included above as archetypical for the kind of interactions occurring in an eco-evolutionary system (23).

The Parameters.

The time-independent parameter functions r, b, α, and g represent the reproductive and competitive traits of individuals possessing the phenotype vector x and are considered intrinsic properties of the multispecies ecological system. This imposed time independence rests on the underlying fundamental assumption that the phenotype space is sufficiently large so that all relevant traits of individuals are included. The arguments of the terms will be suppressed as long as possible to avoid cluttering the equations.

All of our biological conclusions reported below apply for any nonnegative value of b, including b=0. The explicit inclusion of the self-limiting term, b, may seem redundant, as it could be subsumed by the bilinear interaction term by defining a modified interaction term as α¯(x,x′)=α(x,x′)−b(x)δ(x,x′), where δ is the Dirac measure. However, in this case, the modified interaction term α¯ is infinite at x=x′, whereas it is reasonable to assume that the original term α is finite. For simplicity, we therefore keep the two terms separate. Nevertheless, we will allow for b(x)=0 (i.e., no local self-regulation). In this case, mutualistic interaction risks enabling a pair of “runaway” species with exponential growth (Methods, Lemma 2), and the structure of α becomes essential. The presence of a positive self-limitation term, b, allows for simplified criteria on admissible symbiotic interaction without species pairs (or groups) growing exponentially. We further note that the presence of the diffusive term g (i.e., the between-generation genetic variation) ensures that the solution n is continuous; therefore, a purely local self-limiting term b is a reasonable approximation of self-limitation also due to individuals of neighboring traits. It should be noted that if, for instance, α(x,x′) is constant in x′, then the integral ∫Ωα(x,x′)n(x′)dx′ represents the global density; that is, the integral value of α(x,x′), defined below as A(x), defines the limitation due to the global density of species x.

Clustering of Phenotypes into Species.

Our model has no a priori defined species, and the notion of a species can only be considered a posteriori. This implies that we need not consider specific measures by which a new species appear or existing species branch into two separate species—these concepts emerge naturally from Eq. 1. Indeed, due to the structure of Eq. 1 as a nonlocal Fisher–Kolmogorov equation, the solution n(x) will for many parameter functions tend to involve highly nontrivial “clustering” (see ref. 20 for a thorough review). The clustering in the solution corresponds to separate species, similar to the morphological species concept widely used when delimiting species out in nature and in the fossil record. In real ecological systems, species with low abundance may, due to stochasticity, become extinct; this is, however, beyond the scope of this paper.

The Ecological Structure of the Evolutionary Model.

The ecological interactions are defined by the α functions. The ecological interactions can be decomposed into symmetric and asymmetric components. The “symmetric” component of interactions is defined by α(x,x′)+α(x′,x), and measures competitive (−/−) or mutualistic (+/+) interaction. “Asymmetric” interactions are defined by the remaining component α(x,x′)−α(x′,x), and measure the degree of trophic interactions, host–pathogen interactions, and the like (i.e., +/− type). Note that asymmetric competitive (−/−) and mutualistic (+/+) ecological interactions with different reciprocal strengths have both symmetric and asymmetric components. When we refer to a symmetric system, this is one where asymmetric interactions are zero, and thus, α(x,x′)=α(x′,x).

Linking our Approach to Earlier Model Approaches.

Our basic model considering continuous variation in traits thus represents a generalization of a standard multispecies community model (e.g., ref. 24) with the addition of a term, ∇⋅(g∇n), for new genetic variation occurring continuously and randomly (with respect to whatever the local fitness optimum should be) across generations. A diffusion process, similar to the adaptive dynamics approach (see, e.g., ref. 25), is used in the model to represent evolution. With this contribution of ours, we are able to provide analytically derived results, rather than relying on numerical simulation, thus being able to reach more general conclusions.

The main novelty of this new model formulation is to bring together ecological and evolutionary dynamics within a common analytical model framework:

  • i) If there is no between-generation genetic variation (g=0), the model will be a generalization of fairly standard ecological models (e.g., ref. 24) both by making it continuous as well as reducing the emphasis on individual species as we consider the distribution of phenotypes directly; in this case, there will be no evolution.

  • ii) If we allow for between-generation genetic variation (g>0), the ecological model becomes an eco-evolutionary model. When the between-generation variation is significant relative to the intrinsic growth rate of the ecological system, evolution might occur.

Common with many similar earlier models within the field of ecology and evolution (including the adaptive dynamics approach), our model assumes asexual reproduction.

Our Objectives.

We address the following three specific questions:

  • i) Under what conditions are evolution of n bounded, such that the number of individuals n(x) never grows uncontrollably for any trait x?

  • ii) Under what conditions are the stationary solutions of Eq. 1 stable, such that an evolutionarily stable state without any change in the species composition n(x) may arise?

  • iii) Under what conditions, and in what sense, will the species modeled by Eq. 1 continue to evolve (i.e., Red-Queen type of evolution§) and under what conditions will evolution cease (i.e., achieving stasis¶)?

Our results will fully settle question i, but more importantly provide several new results regarding questions ii and iii.

Results

Our results regarding Eq. 1 are collected in Lemma 1, Lemma 2.1, and Lemma 2.2 (Methods), and Lemma 3.1 and Lemma 3.2 (Supporting Information). Together, they characterize the temporal (ecological and evolutionary) dynamics. We summarize the results here, using the notation that the integral of the positive part of the interaction kernel [i.e., the sum of all positive interactions acting on a species with phenotype x, is defined as A2+(x)=∫Ωmax[α(x,x′),0] dx′]. In this summary, the x-dependent parameters are implicitly treated as constants.

  • i) If the level of self-regulation in the system is such that b<A2+, the system is subjected to blowup, which we denote the “unbounded state.” We consider the blowup of the system as unreasonable; this inequality thus imposes a restriction on admissible parameters. From this, we conclude that there must be sufficient self-limiting competition to dominate over synergetic/mutualistic interaction. Sufficient self-limitation can also be achieved by the structure of α (thus allowing for b=0; Lemma 2.1 and Lemma 2.2).

  • ii) If the system is stable as defined in Lemma 1, this implies that a constant (with respect to x) and stable (with respect to time) evolutionary state exists wherein species cannot be identified (i.e., without any peaks in the phenotype space). We denote this as the “constant state.” The structure of α [describing the nature of and strength of the between-species interactions (trophic interaction, competition, and the like)] that allows this to happen would be where there is very little synergetic/mutualistic interaction. We note that because it is the real part of the Fourier transform of α that determines stability, it is only the symmetric interactions that affect stability. Asymmetry in α(x,x′) may thus impact the time evolution, but not the stability of constant state.

  • iii) The case where the system does not blow up (see i above), although unstable in the constant state, is the most biologically interesting. We call this the “species state.” This is the case where the system provides solutions where species may be identified and seen to coexist in a reasonably well-behaved state. For particular parameter regimes, Lemma 3.1 and Lemma 3.2 show that, in the case of only symmetric interactions, stationary solutions corresponding to stasis (defined as a stationary steady state) will arise. It should be noticed that when +/− ecological interactions are absent/missing, the conditions do not correspond to real ecological systems, but might correspond to artificial laboratory systems restricted to only competitive pairs or mutualistic pairs.

  • iv) Our next step is to characterize characterize systems in the species state where the interaction α is nonsymmetric. For homogeneous systems, Lemma 1 shows that, in the linear regime, the solution tends to settle in a continuously evolving state (i.e., Red-Queen type of evolution) depending on the extent and character of asymmetry in the system. It should be noted that asymmetry may exist in mutualistic and competitive ecological interactions as well as the (obvious) trophic and host–parasite type of ecological interactions, because for any two traits x and x′, there may be an asymmetric component of the interaction, even if both α(x,x′) and α(x′,x) are positive (or negative). On this basis, we conclude that asymmetric ecological interactions favor Red-Queen type of evolution.

  • v) Finally, we show, as a corollary to Lemma 2.1 and Lemma 2.2, that for any stable system with a finite trait space Ω, there must exist a state n0(x) for the time evolution of the solution n(x,t), which the time-dependent solution essentially returns to infinitely many times. This state n0 can thus either represent stasis or be interpreted as a limiting factor on the Red-Queen type of evolution.

Altogether, our biological conclusions are as follows:

  • i) An ecological system with only symmetric interactions, such that α(x,x′)=α(x′,x) (each pair of species is affected by the other symmetrically), will, in the absence of external abiotic perturbations, favor stasis. This may be understood intuitively in the sense that, without asymmetric interaction, there is no notion of a preferred direction for continued evolution of the system.

  • ii) An ecological system with asymmetric interactions (such as predator–prey interactions) may lead to both stasis and Red-Queen type of continued evolution: The stronger (or more dominant) the asymmetric ecological interactions (such as +/− interactions and asymmetric −/− and +/+ interactions in the ecological system) are, the more likely Red-Queen type of evolution will be—fully independent of the between-generation genetic variation (quantified by g).#

  • iii) As a corollary, i and ii imply that ecological continued temporal changes (i.e., only temporal variation in the abundances of the interacting phenotypes/species) and evolutionary dynamics (with both temporal variation in the abundance of the interacting phenotypes/species as well as evolution of new phenotypes and extinction of existing phenotypes) are particularly likely if there is asymmetric ecological interactions, including, for instance, asymmetric competitive interaction between species. For evolution to occur, g>0 is required.

Discussion

General.

Our answer to “Will evolution cease if all abiotic perturbations cease?” is that, if all ecological interactions are symmetric (implying no trophic interactions), evolution will cease in an abiotically stable environment. If, however, there are asymmetric interactions, we have shown that the more asymmetry in the ecological system (in the sense of asymmetry with the same frequency in Fourier space as the interspecies distance) the more likely is evolution to continue in a Red-Queen type of fashion, even in the absence of abiotic perturbations, a conclusion profoundly extending the earlier conclusion reached by Stenseth and Maynard Smith (26).

It is worth noting that, in our analysis, the evolutionary pattern (stasis vs. Red Queen) is primarily a property of the ecological system within which a species finds itself [i.e., ∫Ωα(x,x′)n(x)dx′]. This is explicitly the case for the linearized analysis and holds in the general case, both when there is no between-generation variation (i.e., g=0; no evolution), as well as when there is much between-generation variation (i.e., g≫0). This is a theoretical conclusion of profound importance—to our knowledge, no one before has been able to reach a similarly general conclusion within a generalized eco-evolutionary model. Here, it suffices to point out that our analysis links the Red-Queen evolutionary literature profoundly to the food web literature (27, 28). Indeed, this brings ecology and evolution closer together: the evolutionary dynamics being a direct result of the property of the ecological system the species finds itself in, just as with the adaptive dynamics perspective (see, e.g., refs. 13⇓⇓–16 and 25).

We might assume either g=0 (to make it a purely ecological model), or consider g≫0 (to include evolution in the ecological model). It is, however, the combination of both ecological and evolutionary dynamics that applies in nature. Much of the literature on the Red-Queen hypothesis focuses on (+/−) type of ecological interactions, especially trophic and parasite–host interactions (e.g., refs. 9 and 29⇓⇓–32), which is consistent with our results showing that this is indeed the type of asymmetric ecological interactions that would yield Red-Queen type of evolution. The importance of asymmetry has earlier been pointed out by, e.g., refs. 16 and 33.

Interpretation.

The interpretation of the solution n(x,t) in terms of species, lag, and mean and variance in traits (Supporting Information) allows our model framework to be directly linked to the lag–load concept of Maynard Smith (34). Maynard Smith suggested that genetic evolution of a species i will reduce the evolutionary lag Li in proportion to the current lag—that is to say that the following equation holds with a time-independent coefficient βi, a coefficient being equivalent to the coefficient βii in Stenseth and Maynard Smith (26):dLidt=−βiLi(t).[2]Furthermore, Fisher (2) assumed that fitness of a species i is proportional to the genetic variance Σi2. For our model, it is reasonable to assume that, near equilibrium, fitness may be linearly approximated by the evolutionary lag; thus the following equation holds with a time-independent coefficient γi:dLidt=−γiΣi2(t).[3]It is worth noting that the lag–load model and Fisher's fundamental theorem can hold simultaneously, if and only if Σi2(t)=(βi/γi)Li(t). This will only be the case if there is no genetic variance in the optimal state [for which by definition Li(t)=0]. In our model, it is impossible for all individuals of a species to converge to a single trait, because the presence of the g term (i.e., between-generation genetic variation) will ensure that, for every generation, a slight variability in traits remains.

This leads us to the observation that Fisher’s model is applicable at early stages of evolution, which in practice implies the initial response to external environmental forcing. In contrast, the lag–load model is applicable to the later stages of evolution within an environment with fixed external forcing.

Another key concept within the evolutionary literature should be mentioned—the concept of evolutionarily stable strategies (ESS), originally contributed by Maynard Smith and Price (35); see also Maynard Smith (36, 37); for a mathematical definition of this concept for differential-equation population models, see ref. 38. This concept will only apply to the case of stasis or close to stasis—not Red Queen. Furthermore, our results contradict the results given within the context of ESS, where Rosenzweig et al. (39) conclude that Red-Queen–type evolution can only be supported for unbounded traits (see, e.g., refs. 15, 25, and 30 for earlier work on this topic). Their conclusion is obtained within the context of a lag–load formalism, where the species are a priori defined. A complete analysis of the emergence of ESS within our new model remains to be done. We note, however, an important distinction from our framework: in traditional discussions of ESS, one considers a species with an optimized trait, common for all individuals of that species. As discussed after Eq. 3, this is impossible, thus in our model framework an ESS must be considered within the context where the trait variability is the minimum variability permitted by the between-generation genetic variation, evolvability, g.

Further Challenges.

The biological challenges emerging from our study are as follows: (i) interpreting previously published models for the Red Queen within our new model framework; (ii) interpreting the literature on food web coupled with evolution (see, e.g., refs. 40 and 41), determining whether we can understand the different modes of evolution as a consequence of the food web stature; (iii) interpreting paleontological records in light of our current results (see, e.g., refs. 42⇓–44); (iv) linking our model structure to the ecometrics literature within the field of evolutionary biology (see, e.g., refs. 45⇓–47; and (v) undertaking experimental studies using microorganisms to test the conclusion of our analysis—and beyond.

The mathematical challenges emerging from our study are as follows: (i) a further analysis of the case of symmetric ecological interaction, to settle the question as to whether symmetric interaction can preclude Red-Queen type of evolutionary dynamics also for moderate g; (ii) a more complete treatment of the long-term dynamics for the nonlinear evolutionary regime, including the categorization of stationary solutions and limit cycles; (iii) a rigorous understanding of the link between the continuous solution and the identification of individual species; (iv) incorporation of spatial heterogeneity and stochastic environmental forcing (see, e.g., ref. 48) in the model framework; (v) incorporation of sexual reproduction and age/size-structured populations in the model framework; and (vi) efficient computational tools allowing for high-dimensional numerical simulation.

Essentially, our paper contributes to the further development of a theory for understanding large-scale features of the evolutionary dynamics (see also refs. 49⇓⇓–52).

Methods

Justification for the Mathematical Model.

Consider the following multispecies ecological model:∂n∂t=fb(n)−fd(n).[4]

The functional fb(n) describes the rate of birth of individuals with trait x, and similarly the functional fd(n) describes the rate of death. We consider death as proportional to the number of individuals n, and furthermore dependent on interaction and competition, justifying the form:fd(n)=(rd+bn−∫Ωα(x,x′)n(x′)dx′)n.[5]

We consider the number of births from trait x as proportional to n(x) but recognize that the actual newborns may have slightly different traits. If we assume that the traits of the next generation are spread according to a multidimensional Gaussian distribution Nd(x,x′) with covariance matrix d, this leads to the following:fb(n)=rb∫ΩNd(x,x′)n(x′)dx.[6]

However, Gaussian spreading is the solution operator of the diffusion equation; thus, for small variability between generations, we can use a Taylor expansion to obtain the following:fb(n)=rb⁡exp(∇⋅(d∇))n≈rb(n+∇⋅(d∇n)).[7]

By letting g=rbd, and r=rb−rd, we obtain Eq. 1. Note that a generalization to nonlocal, “innovative,” genetic mutations can be incorporated by using a non-Gaussian spreading with fatter tails; in this case, the diffusive term would take the form of a fractal derivative operator, e.g., ∇⋅(d∇νn), for some exponent 0<ν≤1.

We may distinguish two types of phenotypes (represented as dimensions in phenotype space): linear and circular. The linear phenotype is the most intuitive and represents any quantity that is strictly ordered (such as size). Linear phenotypes may be either bounded or infinite. The circular phenotype represents quantities that contain no ends, such as relative hues of color in the red–green–blue scheme, stripe orientation on a zebra, plant growth orientation relative to external factors, or phases of growth relative to the seasons (20).

Due to the presence of a second-order differential term, we assign to the boundary of Ω either Dirichlet or Neumann boundary conditions (53). A Dirichlet boundary condition is equivalent to specifying the species density n at the boundary (typically 0), whereas a Neumann boundary condition is equivalent to specifying the normal component of g(n)∇n (and hence the evolutionary drift, also typically 0) across the boundary. Note that, in the dimensions associated with circular phenotypes, the domain is periodic, and thus the domain does not contain a boundary.

We refer to the ecological system given in Eq. 1 with general parameter functions r(x), b(x), g(x), and α(x,x′) as “heterogeneous”; most of our results hold in this setting. However, to facilitate the presentation and to obtain results in Lemma 1, which rely on the Fourier transform, we will sometimes work with the instructive case of constant coefficients, which we refer to as “homogeneous.” This implies the simplifications r(x)=r, b(x)=b, and g(x)=g. Furthermore, in the homogeneous setting, we assume the function describing the ecological interactions to depend only on signed distance in phenotype space [i.e., α(x,x′)=α(x−x′)]. The assumption of homogeneity implies that we study the interspecies interaction, rather than the impact of the competitive advantages of various locations in the phenotype space Ω.

Stability of Constant Stationary States.

Stationary states for Eq. 1 are obtained by solving for ∂n/∂t=0. These stationary states can be divided into three types: zero, constant, and variable. Here, we define the (negative) integral of the kernel as A(x)=−∫Ωα(x,x′)dx′.

Constant steady-state solutions (with respect to x) are available for homogeneous systems. Then the differential term vanishes, and furthermore the integral term simplifies. We are then left with 0=rn−(b+A)n2. This equation has two solutions, the zero solution, which is of no interest, and the nonzero constant solution, which we denote nC=r/(b+A).

Variable Steady-States Solutions.

Variable steady-states solutions are in general defined as any nontrivial solutions to the nonlinear equation obtained from Eq. 1 by setting ∂n/∂t=0. In general, this equation must be solved numerically.

The zero stationary state will always be unstable, because we assume that r>0 for some x. However, the nonzero constant state nC is less obvious. Indeed, the system has some similarities with the Turing equations (54), and the stability depends on the balance between the nonlocal interaction and diffusion.

We proceed considering the homogeneous system, after which by linearizing [Eq. 1] around nC, we obtain the following equation for the deviation m=n−nC. Here, it is assumed that m is infinitesimal and periodic:∂m∂t=−bnCm+nC(α*m)+gΔm.[8]

The linear coefficient evaluates to zero, since r−(b+A)nC=0; furthermore, we denote the convolution integral by an asterisk (*). Considering the case where we take Neumann boundary conditions, a Fourier transform of Eq. 8 now leads to the following:∂m^(k)∂t=f(k)m^(k).[9]

Here, we have denoted the transformed functions by a hat (i.e., m^ is the Fourier transform of m), and the transform variable is denoted by k. The function f is defined as follows:f(k)≡−bnCδ(k)+nCα^(k)−g|2πk|2.[10]

Based on Eqs. 8 and 9, we now deduce the following stability modes:

  • i) By setting k=0, we see that the system is “locally” unstable if nC(b+A)> 0 (i.e., any constant perturbation will grow exponentially).

  • ii) The system is “nonlocally” unstable if there exists a k>0 such that f(k)>0. This implies that patterns of frequency k will be unstable, even though the system may be locally stable according to i).

  • iii) The system is “globally” stable if nC Re{α^(k)}−g|2πk|2<0 for all k>0.

Note that i is just a special case of ii, where k=0. Furthermore, for convenience, we denote the frequency where f(k) attains its maximum as k*, which we refer to as the “most unstable mode.”

We can deduce more from the linearized analysis. Indeed, the linear structure of Eq. 10 allows us to write the time evolution of the perturbation explicitly as follows:m^(k,t)=m^(k,0)ef(k)(11)

In particular, this expression shows that perturbations decay (and grow) in a stationary way (with respect to motion in the Ω domain) if and only if Im{f(k)}=nCIm{α^(k)}=0. In particular, if Im{α^(k)}≠0, the most unstable (and thus dominant) perturbation will be non-stationary. The following lemma summarizes the situation.

Lemma 1.

The constant stationary state nC is stable if and only if f(k)<0 for all k. Moreover, nonlocally unstable growth is stationary in Ω if and only if Im {α^(k)}=0 for all k where f(k)>0.

The homogenous setting needed for the Fourier transform analysis leads to results where the presence of the evolution term g is strictly stabilizing. However, for heterogeneous parameters, the opposite case may arise, due to the existence of Turing-like instabilities in the system (54) (see Supporting Information for an example).

Boundedness of time evolution:

We now consider the evolution given by equation (1). In particular, we will ascertain under what conditions we can guarantee that a “blow-up” of the solution n cannot occur. This is a prerequisite for the mathematical well-posedness (existence and uniqueness of solutions, etc.) of the model. Indeed, we define the standard Lq norms by the notation ‖n‖q=(∫Ωnq)1/q. Here it is understood that in the present setting, we simplify the notation since it is known that n≥0. Now, by integrating equation (1), we obtainddt‖n‖1≤r+‖n‖1−b−‖n‖22+∫Ωn(y)∫Ωα(x,x′)n(x′)dx′ dy+∫Ω∇⋅(g∇n)dx.(12)Here, we define r+=maxxr(x) and b−=minxb. The third term can be bounded from above by using that n(x)≥0 and the triangle inequality, to obtain∫Ωn(x)∫Ωα(x,x′)n(x′)dx′ dx≤A1+‖n‖22.(13)In order to allow for b=0 we denoteA1+≡maxn(x)>0∫Ωn(x)∫Ωα(x−x′)n(x′)dx′ dx‖n‖22≤maxx∫Ωmax[0,α(x,x′)]dx′≡maxxA2+(x).(14)Furthermore, the differential term can be integrated by parts, after which the boundary terms can be omitted for Neumann boundary conditions. Thus, equation (12) satisfies the inequality:ddt‖n‖1≤r+‖n‖1−b−‖n‖22+A1+‖n‖22.(15)Now, since ‖n‖1≤‖n‖2, it holds that ‖n‖1 is bounded satisfying‖n‖1≤r+b−−A1+≤r+b−−A2+.(16)

Again, we summarize the result in a lemma:

Lemma 2.1.

The time evolution n(x,t) remains bounded in the L1 norm if b−>A1+.

We will, for simplicity of interpretation, typically consider the condition b−>A1+ in Lemma 2.1 (and in similar instances below) in the two specialized cases where either:

  • i) The self-limitation dominates synergetic interaction pointwise, thus, b(x)>A2+(x); or

  • ii) There is no self-limitation, b=0, but interaction has strictly negative eigenvalues, e.g., A1+<0.

Case i) allows us to obtain a stronger result, not available for b=0: It may be of interest to distinguish the total number of individuals, captured by the L1 norm, from the maximum density of individuals with a single trait, which is captured by the maximum norm, denoted L∞. By considering the point x at which n(x,t) attains a maximum (over x), we see that the maximum can remain bounded only if b(x)>A2+(x) for all x. We summarize this in the second part of the lemma.

Lemma 2.2.

The time evolution n(x,t) remains bounded in the L∞ norm if b(x)>A2+(x) for all x.

Lemmas 2.1 and 2.2 have an interesting corollary. Since the solution n(x,t) lives in a space of bounded functions, if furthermore Ω is finite, Eq. 1 is a nonlinear diffusion equation with bounded coefficients, and the space of solutions is compact. This implies that we may choose a sequence of sampling times t={t1,t2,…}, with ti−ti−1>Δt for some time interval Δt, and corresponding solutions n(x,ti). Since n(x,ti) lies in a compact space, any sequence has a sub-sequence which is a Cauchy sequence [see e.g. (55)]. We express this statement in the sense that there exist for any small number ϵ, a sub-sequence of ti of infinitely many sampling times t={tj1,tj2,…} and a fixed point n0(x), such that ‖n(x,tji)−n0(x)‖1≤ϵ for all ji. This shows that even if Red-Queen evolution occurs, the system will return to essentially the same species distribution infinitely many times (resembling convergent evolution).

Evolution Toward Stasis.

In Lemma 1, Lemma 2.1, and Lemma 2.2, we obtain bounds on the parameter values that are required to have biologically plausible evolution. Furthermore, Lemma 1 shows the time evolution in the linearized regime. However, the general question regarding the behavior of solutions in the fully nonlinear regime remains open. Here, we partially answer this question by showing two different cases in which it can be proved that the solution evolves toward stasis. This is achieved by showing that the governing equations have a so-called gradient flow structure when interactions are symmetric and either

  • i) the evolvability is negligible [e.g., g(x)=0; i.e., only ecological dynamics], or

  • ii) the evolvability is sufficiently large [e.g., g(x)≳r(x); i.e., evolutionary dynamics dominating over ecological interaction].

A gradient flow structure is based on the physical analogy where it is possible to define some “potential energy,” and where the system evolves to minimize that energy; gradient flow theory has been previously used to establish stability of biological systems; see, e.g., Hopfield (56) for an early example. If, furthermore, there is dissipation in the system (e.g., friction, in the physical analogy), this implies that the solution will approach a stationary state, i.e., stasis (see, e.g., ref. 57 for an introduction). We provide the details of this approach through Lemma 3.1 and Lemma 3.2 in Supporting Information.

SI: Turing Instabilities

As an example of how evolvability can act in a destabilizing manner, we consider a pathological example, using that the nonlocal integral term can be used to connect separate parts of the domain, which again can be interpreted as a reaction–diffusion equation. To this end, let Ω=[−1,1], and define the species density for the negative and positive parts as n1(x,t)=n(x,t) for x>0, and n2(x,t)=n(x−1,t) for x<0. Define similarly r1, r2, etc., and denote these as constants (e.g., each parameter function takes two values, one for positive and one for negative x). Now if the interaction term has the form α(x,x′)=α1δ(x,x′−1)+α2δ(x,x′+1), then Eq. 1 can be written as a system of local equations:∂ni∂t=rini−bini2+αin1n2+∇⋅(gi∇ni)  for   i=1,2.[S1]To continue the example, we let α1=r1=b1=−α2=1, r2=b2=−1/2. Then the system has a stationary state at n1=n2=1. Proceeding as above, we obtain in place of Eq. 14 the following:∂v^i(k)∂t=∑jfij(k)v^j(k).[S2]The function f is now a matrix, which takes the following form:f(k)≡(−11−11/2)−|2πk|2(g100g2).[S3]This is exactly the type of system Turing studied in the context of chemotaxis. Now, for g1=g2=0, the stationary state is stable, because the eigenvalues of f(k) both have negative real parts. However, when g1=100 and g2=1, the stationary state is unstable, because for k=(4π)−1 there is a positive eigenvalue.

This example shows that it is possible for Eq. 1, depending on the parameter functions, to be stable for g=0, yet unstable for certain g(x)>0. In other words, the presence of evolvability makes an otherwise-stable ecological system unstable.

SI: Evolution Toward Stasis

We refer to the potential energy with the symbol Φ, and the dissipation function (when explicitly needed) as Ψ. Throughout this section, we use the notation n˙=dn/dt. We collect our results below in two lemmas, both of which hold for heterogeneous parameters.

Lemma 3.1.

Let the interaction kernel be symmetric, such that α(x,x′)=α(x′,x), let b−>A1+. Then, if the evolvability is negligible, such that g=0, it holds that (i) Eqs. 1 and 2 are a gradient flow, in the sense of the following:1nn˙=−DΦ(n),[S4]with energy function,Φ(n)=∫Ω−rn+b2n2−n2∫Ωα(x,x′)n(x′) dx′dx.[S5]Furthermore, (ii) all flows will be dissipative and lead to a stationary steady state.

Proof:

By differentiating Ω, we obtain that for symmetric α:DΦ=−r+bn−∫Ωα(x,x′)n(x′) dx′.[S6]We recall that the gradient flow structure guarantees a dissipative process becauseddtΦ(y)=〈∂nΦ(n),n˙〉=−∫Ω1nn˙2dy≤0,[S7]where equality holds if and only if the solution is stationary. The process will converge to a stationary point if the energy function Φ has a global minimum, which is assured if b−>A1+ because the solution cannot take negative values. □

Lemma 3.2.

Let the interaction kernel be symmetric, such that α(x,x′)=α(x′,x), let b−>A1+, and let some part of the boundary have Dirichlet boundary conditions. Then if the evolvability is sufficiently large, such that ming(x)>CΩmaxr(x), where CΩ is the Poincare constant for Ω, it holds that (i) Eqs. 1 and 2 are a gradient flow, in the sense of the following:n˙=∂μΨ∗(n,μ)[S8]μ=−DΦ(n)[S9]with energy functionΦ(n)=n22[S10]and dual dissipation functionalΨ*(n,μ)=∫Ω−r2μ2+g2|∇μ|2+n2μ∫Ω[bδ(x,x′)−α(x,x′)]μ(x′) dx′dx.[S11]Furthermore, (ii) the dual functional is monotone and invertible.

Proof:

It follows from substitution that the dual-gradient structure stated in the lemma corresponds to Eq. 1. The task is to ascertain whether a dissipation functional Ψ(n,ξ) exists such thatΨ∗=ℒ{Ψ},[S12]where ℒ is the Legendre transform. This will be the case if Ψ∗ is strictly monotone, which holds because for the last term b−>A1+, whereas for the first two terms, we obtain from the Poincaré inequality that∫Ωg|∇μ|2dx≥CΩ−1maxg(x)∫Ωμ2dx.[S13]Together with the condition ming(x)>CΩmaxr(x), the lemma holds. □

Lemma 3.1 and Lemma 3.2 together show that symmetry is the key mechanism for stasis, both in the presence of high and low evolvability of the system.

SI: Interpretation of Density n(x) in Terms of Species

When the system is close to its steady state, it is possible to relate our model to common notions such as fitness and genetic variance. We therefore introduce the explicit notion of a species in the current framework. Intuitively, a species i will be associated with a local maximum in the density n. Denote such a local maximum as xi. There are two distinct cases: a local maximum can have (essentially) compact support, that is to say that we can identify a region ωi in phenotype space surrounding xi where n<ϵ on the boundary of the region and no other local maxima exists within xi. Here, ϵ represents the threshold when the density n is so rare that we consider it vanishing for practical purposes. This corresponds to a well-defined species. Conversely, several local maxima may be near each other, making the separation into distinct species more arbitrary.

In the case of well-defined species, we can introduce the following notion of lag for a given time t: let the optimal state ni∞(x,t) denote the species density obtained by evolving Eq. 1 within the species region ωi, while keeping n(x,t) constant for all other species (e.g., for all x not in ωi). We can now define the evolutionary lag for species i (30, 31) as the deviation from the optimal state:Li(t)=‖ni(x,t)−ni∞(x,t)‖1‖ni∞(x,t)‖1.[S14]Here, the double vertical brackets may indicate any integral norm, but, as indicated, most relevant is the L1 norm.

Furthermore, we may define the genetic mean and variability for species i as follows:μi(t)=∫ωixn(x,t)dx  and  Σi2(t)=∫ωi(μi−x)2n(x,t)dx.[S15]Note that the position of the most common trait will be xi, whereas the mean trait will be μi: in general, these two are distinct.

Acknowledgments

Comments provided by Juan Bonachela, Erik Hansson, Steinar Engen, Mikael Fortelius, Jarl Giske, Øistein Haugsten Holen, Lee Hsiang Liow, Stig William Omholt, Eörs Szathmary, Kjetil Lysne Voje, Richard A. Watson, and Meike T. Wortel helped us to sharpen our presentation. Lene Martinsen is thanked for proofreading our manuscript as well as for helping to format the text properly. This work was supported by the core funding of Centre for Ecological and Evolutionary Synthesis, two grants from the Research Council of Norway (227860/Evoque, 214285/Red Queen) (to N.C.S.), and a VISTA grant through the Norwegian Academy of Science and Letters (to J.M.N.).

Footnotes

  • ↵1To whom correspondence should be addressed. Email: n.c.stenseth{at}ibv.uio.no.
  • Author contributions: J.M.N. and N.C.S. designed research, performed research, and wrote the paper.

  • Reviewers: E.S., Parmenides Foundation; and R.A.W., University of Southampton.

  • The authors declare no conflict of interest.

  • ↵†It is worth noticing that Maynard Smith (signed as “A Correspondent”) highlighted the importance of the Van Valen paper in Nature (58) soon after its publication.

  • ↵‡Maynard Smith J, Second International Congress of Systematic and Evolutionary Biology (ICSEB II), July 17–24, 1980, University of British Columbia, Vancouver, BC, Canada.

  • ↵§Red-Queen type of evolution is defined as continuous change in the coexisting species due to selection pressures caused by changes in their biotic environment.

  • ↵¶Stasis is defined as no evolution within any of the coexisting species due to their interactions with their biotic or abiotic environment, but with occasional minimal evolution due to genetic drift.

  • ↵#In the linearized analysis, the asymmetric interactions must have a structure that has frequencies that overlap with the most unstable frequencies. In other words, for the most unstable frequency k*, the imaginary part of the Fourier transform of α cannot be zero if Red Queen is to arise.

  • This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1525395113/-/DCSupplemental.

Freely available online through the PNAS open access option.

References

  1. ↵
    1. Van Valen L
    (1973) A new evolutionary law. Evol Theory 1:1–30
    .
    OpenUrl
  2. ↵
    1. Fisher RA
    (1930) The Genetical Theory of Natural Selection (Clarendon, Oxford)
    .
  3. ↵
    1. Darwin C
    (1859) On the Origin of Species (John Murray, London)
    .
  4. ↵
    1. Maynard Smith J
    (1976) A comment on the Red Queen. Am Nat 110(973):325–330
    .
    OpenUrlCrossRef
  5. ↵
    1. Stenseth NC
    (1979) Where have all the species gone? On the nature of extinction and the Red Queen hypothesis. Oikos 33(2):196–227
    .
    OpenUrlCrossRef
  6. ↵
    1. PJ Greenwood,
    2. Harvey PH,
    3. Slatkin M
    1. Stenseth NC
    (1985) Darwinian evolution in ecosystems: the Red Queen viw. Evolution. Essays in Honour of John Maynard Smith, eds PJ Greenwood, Harvey PH, Slatkin M (Cambridge Univ Press, Cambridge, UK), pp 55–72
    .
  7. ↵
    1. Salthe SN
    (2008) Natural selection in relation to complexity. Artif Life 14(3):363–374
    .
    OpenUrlCrossRefPubMed
  8. ↵
    1. Liow LH,
    2. Van Valen L,
    3. Stenseth NC
    (2011) Red Queen: From populations to taxa and communities. Trends Ecol Evol 26(7):349–358
    .
    OpenUrlCrossRefPubMed
  9. ↵
    1. Brockhurst MA, et al.
    (2014) Running with the Red Queen: The role of biotic conflicts in evolution. Proc Roy Soc Biol Sci 281(1797):20141382
    .
    OpenUrlAbstract/FREE Full Text
  10. ↵
    1. Voje KL,
    2. Holen ØH,
    3. Liow LH,
    4. Stenseth NC
    (2015) The role of biotic forces in driving macroevolution: Beyond the Red Queen. Proc Roy Soc Biol Sci 282(1808):20150186
    .
    OpenUrlAbstract/FREE Full Text
  11. ↵
    1. Vrba ES,
    2. Denton GH,
    3. Partridge TC,
    4. Burckle LH
    1. Vrba ES
    (1995) On the connections between paleoclimate and evolution. Paleoclimate and Evolution with Emphasis on Human Origins, eds Vrba ES, Denton GH, Partridge TC, Burckle LH (Yale Univ Press, New Haven, CT), pp 24–45
    .
  12. ↵
    1. Maynard Smith J
    (1989) The causes of extinction. Philos Trans R Soc Lond B Biol Sci 325(1228):241–252
    .
    OpenUrlAbstract/FREE Full Text
  13. ↵
    1. Kisdi E
    (1999) Evolutionary branching under asymmetric competition. J Theor Biol 197(2):149–162
    .
    OpenUrlCrossRefPubMed
  14. ↵
    1. Kisdi E,
    2. Geritz SA
    (2001) Evolutionary disarmament in interspecific competition. Proc Roy Soc Biol Sci 268(1485):2589–2594
    .
    OpenUrlAbstract/FREE Full Text
  15. ↵
    1. Kisdi É,
    2. Jacobs FJA,
    3. Geritz SAH
    (2002) Red Queen evolution by cycles of evolutionary branching and extinction. Selection 2(1-2):161–176
    .
    OpenUrlCrossRef
  16. ↵
    1. Law R,
    2. Marrow P,
    3. Dieckmann U
    (1997) On evolution under asymmetric competition. Evol Ecol 11(4):485–501
    .
    OpenUrlCrossRef
  17. ↵
    1. Khibnik AI,
    2. Kondrashov AS
    (1979) Three mechanisms of Red Queen dynamics. Proc Roy Soc Biol Sci 264:1049–1056
    .
    OpenUrl
  18. ↵
    1. Ferrière R,
    2. Dieckmann U,
    3. Couvet D
    1. Dieckmann U,
    2. Ferrière R
    (2004) Adaptive dynamics and evolving biodiversity. Evolutionary Conservation Biology, eds Ferrière R, Dieckmann U, Couvet D (Cambridge Univ Press, Cambridge, UK), pp 188–224
    .
  19. ↵
    1. Parvinen K
    (2005) Evolutionary suicide. Acta Biotheor 53(3):241–264
    .
    OpenUrlCrossRefPubMed
  20. ↵
    1. Murray JD
    (2003) Mathematical Biology II: Spatial Models and Biomedical Applications (Springer, New York), 3rd Ed
    .
  21. ↵
    1. Houle D
    (1992) Comparing evolvability and variability of quantitative traits. Genetics 130(1):195–204
    .
    OpenUrlAbstract/FREE Full Text
  22. ↵
    1. Hansen TF,
    2. Houle D
    (2008) Measuring and comparing evolvability and constraint in multivariate characters. J Evol Biol 21(5):1201–1219
    .
    OpenUrlCrossRefPubMed
  23. ↵
    1. Roughgarden J
    (1979) Theory of Population Genetics and Evolutionary Ecology: An Introduction (Macmillan, New York)
    .
  24. ↵
    1. May RM
    (1973) Stability and complexity in model ecosystems. Monogr Popul Biol 6:1–235
    .
    OpenUrlPubMed
  25. ↵
    1. Dieckmann U,
    2. Brännström Å,
    3. HilleRisLambers R,
    4. Ito H
    (2007) The adaptive dynamics of community structure. Mathematics for Ecology and Environmental Sciences, Biological and Medical Physics, Biomedical Engineering, eds Takeuchi Y, Iwasa Y, Sato K (Springer, Berlin), pp 145–177
    .
  26. ↵
    1. Stenseth NC,
    2. Maynard Smith J
    (1984) Coevolution in ecosystems: Red Queen evolution or stasis? Evolution 38(4):870–880
    .
    OpenUrlCrossRef
  27. ↵
    1. Pimm SL,
    2. Pimm JW
    (1982) Resource use, competition, and resource availability in Hawaiian honeycreepers. Ecology 63(5):1468–1480
    .
    OpenUrlCrossRef
  28. ↵
    1. Pimm SL
    (1984) The complexity and stability of ecosystems. Nature 307(5949):321–326
    .
    OpenUrlCrossRef
  29. ↵
    1. Dawkins R,
    2. Krebs JR
    (1979) Arms races between and within species. Proc R Soc Lond B Biol Sci 205(1161):489–511
    .
    OpenUrlAbstract/FREE Full Text
  30. ↵
    1. Marrow P,
    2. Law R,
    3. Cannings C
    (1992) The coevolution of predator-prey interactions: ESSs and Red Queen dynamics. Proc Roy Soc Biol Sci 250(1328):133–141
    .
    OpenUrlAbstract/FREE Full Text
  31. ↵
    1. Nuismer SL,
    2. Ridenhour BJ,
    3. Oswald BP
    (2007) Antagonistic coevolution mediated by phenotypic differences between quantitative traits. Evolution 61(8):1823–1834
    .
    OpenUrlCrossRefPubMed
  32. ↵
    1. Rabajante JF, et al.
    (2015) Red Queen dynamics in multi-host and multi-parasite interaction system. Sci Rep 5:10004
    .
    OpenUrlCrossRefPubMed
  33. ↵
    1. Heard MJ,
    2. Sax DF
    (2013) Coexistence between native and exotic species is facilitated by asymmetries in competitive ability and susceptibility to herbivores. Ecol Lett 16(2):206–213
    .
    OpenUrlCrossRefPubMed
  34. ↵
    1. Maynard Smith J
    (1976) What determines the rate of evolution? Am Nat 110(973):331–338
    .
    OpenUrlCrossRef
  35. ↵
    1. Maynard Smith J,
    2. Price GR
    (1973) The logic of animal conflict. Nature 246:15–18
    .
    OpenUrlCrossRef
  36. ↵
    1. Maynard Smith J
    (1982) Evolution and the Theory of Games (Cambridge Univ Press, Cambridge, UK)
    .
  37. ↵
    1. Maynard Smith J
    (1976) Evolution and the theory of games. Am Sci 64(1):41–45
    .
    OpenUrlPubMed
  38. ↵
    1. Reed J,
    2. Stenseth NC
    (1984) Evolution of cannibalism in an age-structured population. Bull Math Biol 46(3):371–377
    .
    OpenUrlCrossRefPubMed
  39. ↵
    1. Rosenzweig ML,
    2. Brown JS,
    3. Vincent TL
    (1987) Red Queens and ESS: The coevolution of evolutionary rates. Evol Ecol 1(1):59–94
    .
    OpenUrlCrossRef
  40. ↵
    1. Ito H,
    2. Shimada M,
    3. Ikegami T
    (2009) Coevolutionary dynamics of adaptive radiation for food-web development. Popul Ecol 51(1):65–81
    .
    OpenUrlCrossRef
  41. ↵
    1. Takahashi D,
    2. Brännström Å,
    3. Mazzucco R,
    4. Yamauchi A,
    5. Dieckmann U
    (2013) Abrupt community transitions and cyclic evolutionary dynamics in complex food webs. J Theor Biol 337:181–189
    .
    OpenUrlCrossRefPubMed
  42. ↵
    1. Barnosky AD
    (2001) Distinguishing the effects of the Red Queen and Court Jester on Miocene mammal evolution in the northern Rocky Mountains. J Vertebr Paleontol 21(1):172–185
    .
    OpenUrlCrossRef
  43. ↵
    1. Benton MJ
    (2009) The Red Queen and the Court Jester: Species diversity and the role of biotic and abiotic factors through time. Science 323(5915):728–732
    .
    OpenUrlAbstract/FREE Full Text
  44. ↵
    1. Benton MJ
    (2015) Exploring macroevolution using modern and fossil data. Proc Roy Soc Biol Sci 282(1810):20150569
    .
    OpenUrlAbstract/FREE Full Text
  45. ↵
    1. Eronen JT, et al.
    (2010) Ecometrics: The traits that bind the past and present together. Integr Zool 5(2):88–101
    .
    OpenUrlCrossRefPubMed
  46. ↵
    1. Polly PD, et al.
    (2011) History matters: Ecometrics and integrative climate change biology. Proc Roy Soc Biol Sci 278(1709):1131–1140
    .
    OpenUrlAbstract/FREE Full Text
  47. ↵
    1. Fritz SA, et al.
    (2013) Diversity in time and space: Wanted dead and alive. Trends Ecol Evol 28(9):509–516
    .
    OpenUrlCrossRefPubMed
  48. ↵
    1. Sæther BE,
    2. Engen S
    (2015) The concept of fitness in fluctuating environments. Trends Ecol Evol 30(5):273–281
    .
    OpenUrlCrossRefPubMed
  49. ↵
    1. Pelletier F,
    2. Garant D,
    3. Hendry AP
    (2009) Eco-evolutionary dynamics. Philos Trans R Soc Lond B Biol Sci 364(1523):1483–1489
    .
    OpenUrlAbstract/FREE Full Text
  50. ↵
    1. Post DM,
    2. Palkovacs EP
    (2009) Eco-evolutionary feedbacks in community and ecosystem ecology: Interactions between the ecological theatre and the evolutionary play. Philos Trans R Soc Lond B Biol Sci 364(1523):1629–1640
    .
    OpenUrlAbstract/FREE Full Text
  51. ↵
    1. Kendall BE
    (2015) Some directions in ecological theory. Ecology 96(12):3117–3125
    .
    OpenUrlCrossRef
  52. ↵
    1. Power DA, et al.
    (2015) What can ecosystems learn? Expanding evolutionary ecology with learning theory. Biol Direct 10(1):69
    .
    OpenUrlCrossRefPubMed
  53. ↵
    1. Evans LC
    (2010) Partial Differential Equations (American Mathematical Society, Providence, RI), 2nd Ed
    .
  54. ↵
    1. Turing AM
    (1952) The chemical basis of morphogenesis. Phil Trans R Soc B 237(641):37–72
    .
    OpenUrlAbstract/FREE Full Text
  55. ↵
    1. Rudin W
    (1987) Real and Complex Analysis (McGraw-Hill, New York), 3rd Ed
    .
  56. ↵
    1. Hopfield JJ
    (1984) Neurons with graded response have collective computational properties like those of two-state neurons. Proc Natl Acad Sci USA 81(10):3088–3092
    .
    OpenUrlAbstract/FREE Full Text
  57. ↵
    1. Mielke A
    (2011) A gradient structure for reaction–diffusion systems and for energy-drift-diffusion systems. Nonlinearity 24(4):1329–1346
    .
    OpenUrlCrossRef
  58. ↵
    A Correspondent (1973) The Red Queen. Nature 245:353–354
    .
PreviousNext
Back to top
Article Alerts
Email Article

Thank you for your interest in spreading the word on PNAS.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Asymmetric ecological conditions favor Red-Queen type of continued evolution over stasis
(Your Name) has sent you a message from PNAS
(Your Name) thought you would like to see the PNAS web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Red Queen’s hypothesis
Jan Martin Nordbotten, Nils C. Stenseth
Proceedings of the National Academy of Sciences Feb 2016, 113 (7) 1847-1852; DOI: 10.1073/pnas.1525395113

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
Red Queen’s hypothesis
Jan Martin Nordbotten, Nils C. Stenseth
Proceedings of the National Academy of Sciences Feb 2016, 113 (7) 1847-1852; DOI: 10.1073/pnas.1525395113
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Mendeley logo Mendeley

Article Classifications

  • Biological Sciences
  • Evolution
  • Physical Sciences
  • Applied Mathematics
Proceedings of the National Academy of Sciences: 113 (7)
Table of Contents

Submit

Sign up for Article Alerts

Jump to section

  • Article
    • Abstract
    • The Model
    • Results
    • Discussion
    • Methods
    • SI: Turing Instabilities
    • SI: Evolution Toward Stasis
    • SI: Interpretation of Density n(x) in Terms of Species
    • Acknowledgments
    • Footnotes
    • References
  • Figures & SI
  • Info & Metrics
  • PDF

You May Also be Interested in

Setting sun over a sun-baked dirt landscape
Core Concept: Popular integrated assessment climate policy models have key caveats
Better explicating the strengths and shortcomings of these models will help refine projections and improve transparency in the years ahead.
Image credit: Witsawat.S.
Model of the Amazon forest
News Feature: A sea in the Amazon
Did the Caribbean sweep into the western Amazon millions of years ago, shaping the region’s rich biodiversity?
Image credit: Tacio Cordeiro Bicudo (University of São Paulo, São Paulo, Brazil), Victor Sacek (University of São Paulo, São Paulo, Brazil), and Lucy Reading-Ikkanda (artist).
Syrian archaeological site
Journal Club: In Mesopotamia, early cities may have faltered before climate-driven collapse
Settlements 4,200 years ago may have suffered from overpopulation before drought and lower temperatures ultimately made them unsustainable.
Image credit: Andrea Ricci.
Steamboat Geyser eruption.
Eruption of Steamboat Geyser
Mara Reed and Michael Manga explore why Yellowstone's Steamboat Geyser resumed erupting in 2018.
Listen
Past PodcastsSubscribe
Birds nestling on tree branches
Parent–offspring conflict in songbird fledging
Some songbird parents might improve their own fitness by manipulating their offspring into leaving the nest early, at the cost of fledgling survival, a study finds.
Image credit: Gil Eckrich (photographer).

Similar Articles

Site Logo
Powered by HighWire
  • Submit Manuscript
  • Twitter
  • Facebook
  • RSS Feeds
  • Email Alerts

Articles

  • Current Issue
  • Special Feature Articles – Most Recent
  • List of Issues

PNAS Portals

  • Anthropology
  • Chemistry
  • Classics
  • Front Matter
  • Physics
  • Sustainability Science
  • Teaching Resources

Information

  • Authors
  • Editorial Board
  • Reviewers
  • Subscribers
  • Librarians
  • Press
  • Site Map
  • PNAS Updates
  • FAQs
  • Accessibility Statement
  • Rights & Permissions
  • About
  • Contact

Feedback    Privacy/Legal

Copyright © 2021 National Academy of Sciences. Online ISSN 1091-6490