Limits on the evolutionary rates of biological traits

This paper focuses on the maximum speed at which biological evolution can occur. I derive inequalities that limit the rate of evolutionary processes driven by natural selection, mutations, or genetic drift. These rate limits link the variability in a population to evolutionary rates. In particular, high variances in the fitness of a population and of a quantitative trait allow for fast changes in the trait’s average. In contrast, low variability makes a trait less susceptible to random changes due to genetic drift. The results in this article generalize Fisher’s fundamental theorem of natural selection to dynamics that allow for mutations and genetic drift, via trade-off relations that constrain the evolutionary rates of arbitrary traits. The rate limits can be used to probe questions in various evolutionary biology and ecology settings. They apply, for instance, to trait dynamics within or across species or to the evolution of bacteria strains. They apply to any quantitative trait, e.g., from species’ weights to the lengths of DNA strands.


I. INTRODUCTION
Fisher's theorem of natural selection relates the rate of change in the average fitness of a population with the variability in fitness.It holds for evolutionary processes driven by natural selection [1].Fisher's result suggests that variability serves as a resource by enabling fast evolution.However, the result is of rather limited validity: it does not apply to types that mutate or in the presence of genetic drift [2][3][4].Moreover, Fisher's focus was on the change in fitness.However, in countless instances biologists are interested in other quantitative attributes-or traits-of individuals in a population.Here, I extend Fisher's results by studying (a) the rates of arbitrary biological traits and (b) general evolutionary processes that incorporate mutations and genetic drift.
This paper focuses on the rate d⟨A⟩ dt at which the average ⟨A⟩ of a quantitative trait A changes.A can represent any measurable trait in a population.For example, A could be breeds' weights within a species, the maximum CO2 concentration at which different species can survive, or DNA lengths in bacteria strains.The results apply to a range of settings in evolutionary biology and ecology where one is interested in evolution of traits.
Traits' evolution rates have been studied extensively in quantitative biology.Reference [5] focuses on the biological factors that influence maximum growth rates.References [6,7] and [8,9] study how genetic variance and a population's structure affects trait evolution, respectively.There's also extensive data-based work on traits' evolutionary rates.As one example, Ref. [10] studies the maximum growth rates in mammals.
In the field of applied mathematics, extensions of Fisher's theorem have also been considered.Reference [11], for example, revises Fisher's results by studying the effect of mutations on the change in the average * lpgp@lanl.govfitness of a population (note, though, that it does not focus on other traits).Reference [12] includes an interesting mathematical generalization of Fisher's theorem to arbitrary traits and dynamics beyond natural selection.References [13,14] rely on uncertainty relations from stochastic thermodynamics to study biological processes.
To derive limits on biological evolutionary rates, I will leverage techniques that have proven useful to study the maximum speed of physical processes.The most related results have appeared in Refs.[15,16].Both articles rely on information theory to bound the evolutionary rates of arbitrary traits.The results in Ref. [15] hold for arbitrary processes.However, they do not discriminate how different evolutionary forces affect rates.In contrast, the results in this article isolate the contributions of natural selection, mutations, and genetic drift to evolutionary rates.The results in Ref. [16] separate contributions from natural selection and mutations, but do not account for genetic drift.Moreover, the results in Refs.[15,16] involve information-theoretic quantities (versions of the Fisher information) that can be hard to evaluate unless one possesses enough knowledge about the dynamics of the system.In contrast, the main results in this article depend on averages and standard deviations that are often more accessible from experimental data.
The main outcome of this paper is a set of inequalities that constrain the evolution rate of any quantitative biological trait A (e.g., of a particular phenotype) in terms of simple properties of the system of interest.Specifically, knowledge of expectation values and variances of A and of the fitness f of a population suffices to evaluate the inequalities (see details in Secs.IV and V).In this way, slowly evolving traits can be discriminated from those that can rapidly change without the need to exactly solve the complex dynamics of the system, as pictorially illustrated in Fig. 1.
In Section II, I review general limits on traits' evolutionary rates that hold for any biological system.These results are very general but can sometimes be hard to evaluate.In Sections III, IV, and V, I derive rate limits that hold for systems driven by (a) natural selection, (b) natural selection and mutations, and (c) natural selection, mutations and stochastic forces (genetic drift), respectively.The results in Secs.III-V are tailored to evolutionary processes and can be evaluated with knowledge of averages and standard deviations.I conclude in Sections VI and VII with potential connections to experiments and a discussion.
Non-technical summary.How fast can biological evolution occur in nature?What properties of a population enable drastic changes?To what extent are evolutionary rates affected by different driving forces such as natural selection, mutation, or random processes?I address these questions mathematically by deriving inequalities that limit the speed of evolutionary processes.The inequalities take the form of trade-off relations that relate the rate of change of a biological trait with its variance and with the variance in the fitness of a population.In this way, slowly evolving traits can be discriminated from those that can rapidly change.Evaluating the inequalities can be done without exactly solving the complex dynamics of the system.In broad terms, the results in this paper quantify the extent by which variability in a trait can lead to faster evolution.

II. RATE LIMITS ON BIOLOGICAL PROCESSES
Consider a set of types denoted by indexes {j} with an evolving number of individuals n j = n j (t) and a total (evolving) population N = j n j (I will also analyze results in terms of frequencies of types in the next section).Let r denote the growth rate of the population, with components r j := ṅj /n j for type j [17,18].Then, the rate of change in the expectation value ⟨A⟩ = j a j n j /N of a quantitative trait A satisfies The covariance cov(A, B) := ⟨AB⟩−⟨A⟩⟨B⟩ characterizes the correlations between two quantities A and B [17].See the Supplementary Material for the proof of Eq. ( 1).Throughout this work, I use da/dt or ȧ interchangeably to denote time derivatives.
In the context of evolutionary biology, Eq. ( 1) is known as the (time-continuous) Price equation [17,[19][20][21].Equation ( 1) is a mathematical identity that holds under very general assumptions (e.g., differentiability) about the underlying dynamics of the system.Similar equations describe the dynamics of classical systems with evolving probability distributions [22] and of open quantum systems [23].Reference [24] discusses the connections between the Price equation in biology and statistical physics, while Ref. [12] contains a detailed analysis of the way that different biological factors contribute to each term in the Price equation.
The first term in the Price equation ( 1), ⟨ Ȧ⟩ = j ȧj n j /N , describes changes in ⟨A⟩ due to explicit time dependence in the values a j of the trait.The second term in the Price equation, j a j d dt (n j /N ) = cov(A, r), corresponds to the change in the average ⟨A⟩ of a trait due to the population changes.Using that the covariance is bounded by the product of standard deviations [17,25], the latter term is constrained by Here, are the variances of the trait A and of the growth rate r, respectively.(The standard deviation σ A is the square root of the variance.) The inequality in Eq. ( 2) sets a fundamental tradeoff between the rate of change of a biological trait A and (i) the variability σ A in the trait and (ii) the variability σ r in the growth rate r with which the populations {n j } change: fast evolution requires a variable trait and variable population growth rate.I equate the variability of a trait A with its standard deviation σ A = ⟨A 2 ⟩ − ⟨A⟩ 2 [25].Note that, from Eq. (1), one concludes that the sign of d⟨A⟩−⟨ Ȧ⟩ depends on whether the trait A is positively or negatively correlated with the growth rate r.
Trade-off relations like Eq. ( 2), typically referred to as speed limits in physics, also constrain the dynamics of quantum [23,26,27] and classical [22] physical systems.In the latter setting, the variance of r is replaced by the Fisher information I F .The Fisher information is a measure of the speed with which a probability distribution evolves.For a time-dependent probability distribution p j , it is given by I F := j p j ṗj pj 2 [28].
Defining p j := n j /N as the frequency of occurrence of type j, I prove in the Supplementary Material that indeed σ r = √ I F .Related limits to biological systems in terms of the Fisher information have been considered in Refs.[15,16].
While the constraint in Eq. ( 2) is extremely general, its practical usefulness may be hindered by the difficulty in relating the variability in the growth rate σ r (or, equivalently, the Fisher information) to the relevant parameters that govern the dynamics of a concrete system.My goal is to derive bounds (mathematical inequalities) on the rate of change of biological quantities -or rate limits for short-tailored to evolutionary processes.The rate limits apply to populations that evolve under mutations, natural selection, and random genetic drift.Mutations take a preponderant role in biology, enabling evolutionary changes that can lead to new types.Such new types can have different fitness than their ancestors.Natural selection is the process by which the frequencies of fitter types tend to increase relative to their peers in a given environment (in the illustration, turtles became fitter to survive in a drastically-changed environment.)A population's evolution can also be affected by chance, known as genetic drift (in the illustration, a random event reduced the population of an ant colony.)The results in this Letter imply that, under natural selection, a trait A with high variability σA can evolve faster than a trait B with small variability σB.At the same time, the trait with higher variability is more susceptible to the effect of random fluctuations in populations due to genetic drift.In this way, the rate limits can be used to discriminate quantifiable traits in terms of their maximum evolutionary rates and their responsiveness to different evolutionary forces.

III. LIMITS TO REPLICATOR EVOLUTIONARY PROCESSES
Under the assumption that mutation rates between types are negligible, the replicator equation, can be used to model population dynamics [17,29,30] (see Refs. [31][32][33][34] for applications of the replicator equation to various other fields).Here, p j := n j /N is the frequency of occurrence of type j, and the fitness f j ≡ f j ({p k }, t) characterizes whether the frequency p j of a type increases or decreases: the populations of types with positive excess fitness, f j ≥ ⟨f ⟩, tend to grow relative to their peers [35].Using Eq. ( 4), it holds that r j = ṗj /p j + Ṅ /N = f j − ⟨f ⟩ + Ṅ /N , which in turn leads to σ r = σ f .Thus, Eq. ( 2) implies that any evolutionary process that can be modeled by the replicator equation is constrained by where I used that the covariance is invariant under the addition of uniform functions [25].The rate of change in natural selection processes is thus limited by the variability of the fitness of the population and the variability of the quantity of interest.
The rate limit in Eq. ( 5) for the replicator equation implies constraints on the dynamics of arbitrary quantitative traits A. If A has no explicit time dependence, i.e., if the a j 's are constant, then Eq. ( 5) becomes a bound on the evolutionary rates d⟨A⟩/dt, discriminating slowly evolving traits from those that can change rapidly.In plain terms, the inequality says that evolution is slow for systems with homogeneous fitness functions, where σ f ≈ 0 (neutral selection regime).In contrast, evolution can be faster on systems with a diverse population such that σ f is large (natural selection regime).This mathematically formalizes and quantifies the common understanding that diversity serves as an evolutionary resource [6,7,36], in this case by enabling fast evolution whenever a trait and fitness have variability across a population.These variabilities can only occur in sufficiently diverse populations.Equation ( 5) involves terms that may be reminiscent of the breeder's equation.In the breeder's equation, ∆⟨A⟩ = Sh 2 , the net change ∆⟨A⟩ in a trait is governed by a measure of heritability (h) and the selection coefficient S. The selection coefficient measures covariance between fitness and a trait, so S = cov(A, f ) [37].Then, Eq. ( 5) relates S in the breeder's equation to the trait and fitness variabilities: |S| ≤ σ A σ f .This may be useful in scenarios where the heritability is known but the selection coefficient is not, or, possibly, to study trait changes beyond the regime of applicability of the breeder's equa-tion [38].(Note, in particular, that the breeder's equation is less general than Price's [37].) It is natural to wonder whether the left-hand and righthand sides of the inequality (5) are similar (i.e., whether the bound is saturated).When this happens, knowledge of the standard deviations σ A and σ r suffices to estimate the evolutionary rate.This happens whenever A has a linear relationship with the growth rate r, i.e., a j ∝ r j + c where c is independent of j, in which case cov(A, r) = σ A σ r [22].This is the case for the fitness function under replicator dynamics, so Eq. ( 1) yields This corollary of the general rate limit (5) provides a simple proof of Fisher's fundamental theorem of natural selection [1][2][3][4].It shows that Fisher's claim is exact for (i) evolutionary processes modeled by a replicator equation with (ii) fitness functions that are independent of time, in which case d⟨f ⟩ dt = σ 2 f .In situations with more general fitness functions f j = f j ({p k }, t), Eq. ( 6) provides a generalized version of Fisher's theorem whereby the velocity with which fitness changes due to changes in population frequencies equals the fitness variance.

IV. LIMITS TO EVOLUTIONARY PROCESSES WITH MUTATIONS
Mutations are a crucial driving force in realistic evolutionary processes [39][40][41].Mutations between types can be described by the replicator-mutator, or quasispecies model: Q kj ≥ 0 is a dimensionless transition matrix that models mutations between types, which satisfies j Q kj = 1 [42][43][44][45].The replicator equation ( 4) is recovered when the mutation matrix is the identity, Q kj = δ kj .
The general inequality (2) holds in this case, too.Note, though, that while for the replicator dynamics the variance in the growth rate equals the variance in fitness, σ 2 r = σ 2 f , this is no longer the case for dynamics with mutations.Relating the variance in the growth rate (or equivalently, the Fisher information) to biologically relevant quantities in concrete settings remains an interesting problem to be explored.
Alternatively, I define the mutation-driven distribution I interpret Π as the frequency with which a given type would hypothetically occur in the future if evolution were only driven by mutations, or, possibly more biologically relevant, in regimes where strong mutation dominate over natural selection processes [39].Note that Π = p in the mutation-less regime.Then, I prove in the Supplementary Material one of the main results of this work: ⟨A⟩ Π := j Π j a j and σ Π A := ⟨A 2 ⟩ Π − ⟨A⟩ 2 Π are the average and the standard deviation evaluated in the mutation-driven distribution Π.I emphasize that the general result in Eq. ( 2) does not imply the rate limit (9).The replicator-mutator evolutionary model, given by Eq. ( 7), was crucial to derive the latter bound.
The rate of change of any trait A is thus constrained by the quantity's standard deviation evaluated in the mutation-driven distribution Π, and the standard deviation of the fitness of the system.As in the mutation-less setting, diversity in the population is seen to give rise to less constrained evolution rates.
One could be puzzled by the appearance of a term in Eq. ( 9) that directly depends on the average fitness ⟨f ⟩ of the population and not just on relative fitness values.After all, shifting the fitness by an additive constant f j → f j + c in the replicator equation (4) does not affect population dynamics.However, this is not the case in the mutator-replicator equation (7), where an additive constant c on the fitness function leads to a change in ṗj of c(Π j − p j ), which depends on how the distribution p and mutation-driven distribution Π differ.The absolute values of the fitness function play a dynamical role in a system with mutations, and this is manifested in the rate limit (9).
To simplify interpretation, let us momentarily consider the case when A does not explicitly depend on time, i.e., the a j 's are constant.Then, Eq. ( 9) and the triangle and reverse triangle inequalities imply upper and lower rate limits, Here, one can identify two distinct sources that contribute to the evolutionary rates of the system.One of the sources, σ Π A σ f , involves the standard deviation of the fitness function, and the standard deviation of the trait of interest evaluated in the mutation-driven distribution.The remaining term depends on the averages of fitness and the quantity of interest.This is somewhat reminiscent of the speed limits for open quantum systems, where two distinct sources to the dynamics of a system lead to additive contributions to the ultimate speed with which a quantity can evolve, which in turn allows to derive lower bounds on speed [23].A different bound on evolutionary rates under natural selection and mutations was derived in Ref. [16] in terms of contributions to the Fisher information (i.e., more in the spirit of Sec.II).In contrast, evaluating Eqs. ( 9) and (10) only requires knowledge of averages and standard deviations.
The results in this article can be used to discriminate traits with fast or slow evolutionary rates and how different evolutionary forces affect them.Given knowledge of expectation values {⟨f ⟩, ⟨A⟩, ⟨A⟩ Π } and standard deviations {σ f , σ A , σ Π A }, the right-hand side of Eqs.(10a) and (10b) constrain the evolutionary rates of a trait A. The rate limits allow mathematically identifying two extremes of evolutionary regimes: (i) The maximum evolution rate of a trait that is robust to mutations, with ⟨A⟩ Π ≈ ⟨A⟩, is bounded by the variabilities in the trait and fitness; (ii) Meanwhile, a population with uniform fitness f j ≈ f k , for which σ f ≈ 0, leads to constrained dynamics where |d⟨A⟩/dt| ≈ ⟨f ⟩ ⟨A⟩ Π − ⟨A⟩ .Natural selection dominates in case (i), often identified as the weak or neutral selection regime [46], while mutations dominate in case (ii).I illustrate this in Fig. 2. (Also see the Supplementary Material for a model illustrating these claims and the inequalities derived in this article.) We can use Eqs.( 9) and (10) to study the evolution rates of quantities often considered to characterize populations.The entropy S := − j p j ln p j , also known as the Shannon-Wiener index in ecology, has been used to measure the diversity in a population [47].While S ≈ 0 if only one type k occurs with p k ≈ 1, one has S = ln N if N types are equally likely to occur.Taking {a j ≡ I j := − ln p j } and using that Ṡ = − j ṗj ln p j from conservation of probability, Eq. ( 9) implies that Here, the relative entropy (sometimes called the Kullback-Leibler divergence) S p∥Π := − j p j ln Πj pj serves as a proxy for the distance between the distribution p and the mutation-driven distribution Π [48][49][50][51].When S p∥Π is small it is hard to distinguish p from Π.Note that σ Π I ≥ 0 but it is unbounded from above.In cases with an homogeneous population (σ f ≈ 0) and finite σ Π I , mutations drives variability since the entropy only evolves due to the mismatch between the two distributions, with a rate dS dt = ⟨f ⟩ S p∥Π .In the mutationless case in which p = Π, Theorem 8 of Ref. [52] implies Then, the maximum entropy rate scales as max Ṡ ≈ ln N σ f /2 for large N .
As a second example, consider the rate of change of the average fitness function for dynamics that incorporate mutations.Equation (9) becomes which coincides with Eq. ( 6) in the mutation-less regime or for neutral mutations that do not influence the fitness landscape, since then ⟨f ⟩ Π = ⟨f ⟩ and σ Π f = σ f .
Equation (12) imposes the most stringent constraints on evolution when the fitness is completely certain when evaluated in the distribution p or in the mutation-driven distribution Π, i.e., when σ Π f = 0 or σ f = 0.This is the neutral selection strong mutation regime (center column of Fig 2 ), when types have comparable fitness and natural selection is not a strong driver of evolution.In these cases, dynamics are due to the difference in fitness between the two distributions, and d⟨f ⟩/dt − ḟ = ⟨f ⟩ ⟨f ⟩ Π − ⟨f ⟩ , that is, mutations dominate.This also illustrates that beneficial and deleterious mutations can be naturally characterized by ⟨f ⟩ Π − ⟨f ⟩ ≥ 0 and ⟨f ⟩ Π − ⟨f ⟩ ≤ 0, respectively, depending on the change in the average fitness of the population.Note that the contribution to the fitness rate depends on the population's average fitness but not on its standard deviation-Fisher's theorem of natural selection does not hold in this regime.
In contrast to the neutral selection regime, faster evolutionary rates are possible for systems with variable fitness landscapes.In this way, rate limits can be used to mathematically formalize evolutionary regimes according to the dynamical forces in action.I will explore this further at the end of the next section.

V. LIMITS TO STOCHASTIC EVOLUTIONARY PROCESSES
Equations ( 9) and ( 12) provide generalizations of Fisher's fundamental theorem of natural selection to replicator-mutator dynamics in terms of universal constraints on evolutionary rates.However, the replicator and replicator-mutator equations are only simplified deterministic models for the dynamics of infinite populations.In realistic scenarios, stochastic forces typically known as genetic drift play a preponderant role in evolutionary processes [53,54].
To account for genetic drift, I consider a stochastic replicator-mutator equation, + p j γ j dW j − l γ l p l dW l , as a model for stochastic evolutionary processes.The first two terms coincide with the replicator-mutator Eq. ( 7) and describe natural selection and mutation dynamics.The last term models genetic drift driven by noise terms dW j .
In the stochastic difference equation ( 13), the population changes by dp j during a time-step dt, over which a stochastic Wiener noise dW j randomly affects the system.The strength of the genetic drift of population j is characterized by the 'volatility parameter' γ j [55], while the term l γ l p l dW l ensures that the frequencies  9) implies that d⟨A⟩/dt − ⟨ Ȧ⟩ ≈ ⟨f ⟩(⟨A⟩Π −⟨A⟩).
---[(c) -Genetic drift regime] The variability of a trait determines its susceptibility to changes due to genetic drift.In the regime ∥γ∥∞σA ≫ σ Π A σ f , ⟨f ⟩|⟨A⟩Π −⟨A⟩| where the stochastic force dominates, the total change in a trait over a time τ is bounded by the integrated variance of the trait and the maximum drift strength.
p j remain normalized.Following the rules of Itô calculus, the zero-mean Wiener noises dW j satisfy dW 2 j = dt and dW j dW k = δ jk dt, where F represents the average of a function F (p) over realizations of the stochastic noise [56].
Dynamics are thus driven by natural selection, mutations, and random genetic drift.Note that the latter dominate dynamics for infinitesimally short times [56].Therefore, in this case, it is more meaningful to focus on integrated changes 1 τ τ 0 d⟨A⟩ − ⟨dA⟩ during a time interval τ rather than on rates of change.The prefactor 1 τ , which makes the quantity a time-average, ensures the same units as in Eqs. ( 5) and (9).
Then, I prove in the Supplementary Material that, where ∥γ∥ ∞ := max j {γ j }.Equation (14) shows that the noise-averaged change in ⟨A⟩ with respect to the change τ 0 jk p k Q kj f k −⟨f ⟩ a j dt due to natural selection and mutations, is bounded by the variability σ A and by the maximum strength ∥γ∥ ∞ of the stochastic forces.(One can apply the same techniques I used to derive Eq. ( 14) to other Wiener-noise dynamical models [53,57,58], see the Supplementary Material for one such example.) In plain terms, Eq. ( 14) provides a criteria to discriminate situations in which biological quantities are resilient against genetic drift.A trait's high variability can make its dynamics more susceptible to changes due to stochastic forces.In contrast, traits with small variability σ A ≈ 0 over a period τ evolve as if natural selection and mutations were the only driving evolutionary forces.
The final result of this paper singles out the contributions of natural selection, mutations, and genetic drift to an inequality on the noise-averaged change of a trait over a time interval τ , There are many ways in which evolutionary regimes can be (often phenomenologically) classified depending on their predominant driving forces on a given model [46,[59][60][61][62]. Equation ( 15) yields a formal way to mathematically demarcate such regimes by comparing the relative strengths between the sources that dominate the maximum evolutionary rate: (a) strong selection regime In words, the characteristic regimes depend on the values of three quantities: the product of the fitness and trait standard deviations (σ Π A σ f ), the mutation-driven change in a trait (⟨f ⟩ ⟨A⟩ Π −⟨A⟩ ), and the genetic drift's intensity times the trait standard deviation (∥γ∥ ∞ σ A ). (a) The strong selection regime occurs when the product of the fitness and trait standard deviations are significantly larger than the other two factors.Under this regime, mutation and genetic drift can be ignored.(b) The strong mutation regime occurs when the mutation-driven change in a trait (weighed by the average fitness) is the larger term.In it, natural selection and genetic drift are negligible.(c) Finally, the genetic drift regime is identified by the intensity of the genetic drift times the trait's standard deviation being the significantly larger term.In it, the trait's evolution is solely due to random events.
In each of the three regimes detailed above, dynamics of a trait is provably due to the corresponding biological drive.Note that these are trait-dependent criteria -as conveyed by Fig. 2, certain traits are more responsive than others to a given evolutionary force.

VI. POTENTIAL IMPLICATIONS TO EXPERIMENTS
In this section, I discuss how the theoretical inequalities derived in Secs.II-V could be contrasted with the evolution of concrete biological systems.I hope to bridge the results in this article, which borrow ideas from physics, to questions of interest in biology.
In any experiment where one has access to the frequencies p j = n j /N with which types occur and to the variability of a trait, evaluating the results in Sec.II is straightforward.In this case, evaluating the left-and right-hand sides of Eq. ( 2) from experimental data would inform whether a trait evolves as fast as possible.A trait that saturates the inequality is exploiting all biological resources to evolve as fast as nature allows.
Evaluating the rate limits on systems driven predominantly by natural selection and mutations (Sec.IV) will involve purposely designed experiments.Some of the quantities in Eq. ( 9) are likely to be directly available from experimental data (e.g., ⟨A⟩ and ⟨f ⟩).However, ⟨A⟩ Π and σ Π a are the trait average and standard deviation evaluated in the mutation-driven frequency distribution Π [Eq.( 8)].Evaluating them would require an experiment where the intensity of mutations can be "cranked up" relative to that of natural selection.Experiments with microbial evolution and growth arena (MEGA) plates, introduced in Ref. [63], could be used for this purpose.Reference [64] discusses another method to influence the fitness landscape and favour certain mutations.
Finally, Sec.V includes the effect of random stochastic influences (genetic drift).Climate change is believed to affect the occurrence and impact of extreme weather events [65,66].Let us consider modeling the effects of extreme events on a population by stochastic forces, as in Sec.V.Then, we can ask: What properties of a biological system make it resilient against the effects of climate change?Relying on the conclusions drawn after Eq. ( 14), one can argue that traits with less variability will remain unaffected by extreme events for longer.This is because the right-hand side of Eq. ( 14) is small for traits with low variability.Traits with high variability can be affected more.Note that this doesn't necessarily mean that the former case yields more resilience against climate change than the latter: it is possible that populations will need to evolve fast to survive and, if so, Eq. ( 14) suggests that variability may help [67,68].

VII. DISCUSSION
Recent works, mostly within the fields of quantum physics [23,27,69,70] and classical statistical mechanics [71][72][73][74][75][76], but also biology [14][15][16], derived trade-off relations that constrain the dynamics of observables.I find it remarkable that constraints that bound speed, while being extremely general, are saturated in certain paradigmatic cases.In stochastic thermodynamics, for example, the rate at which heat is exchanged with a system and the rate at which a system's entropy changes saturates their speed limits for Gibbs states with (arbitrarily) timedependent temperature [22].In an unrelated setting, certain quantum algorithms have been shown to compute as fast as allowed by speed limits on quantum annealing [77].
Here, replicator dynamics-a model often used to describe mutation-less population dynamics in evolutionary biology-have been found to saturate the rate limit for the fitness of a population (the left and right-hand sides of Eq. ( 5) coincide for A = f ).Higher variability in fitness leads to higher evolutionary rates.This last fact was known by Fisher, who connected the rate of change of the average fitness to the variability in the fitness of a population.
More generally, I have shown that connections between the variability in fitness and evolutionary rates hold for dynamics that incorporate mutations and genetic drift.Moreover, these connections can be made for the evolution of arbitrary traits.
Then, the main strengths of our results are that (i) Eq. ( 9) generalizes the constraints on evolutionary rates to arbitrary quantifiable traits and not just the fitness function, (ii) the inequality in Eq. ( 12) concisely shows how fitness variability influences dynamics, in the spirit of Fisher's theorem, and that (iii) Eq. ( 14) accounts for stochastic evolutionary forces.See Refs.[11,15,16] for related work generalizing Fisher's theorem of natural selection, and deriving distinct speed limits on evolutionary processes driven by natural selection and mutations.Speed limits that incorporate stochastic dynamics, however, are mostly unexplored.
There is an extensive literature devoted to understanding the factors that limit biological evolutionary rates [5][6][7][8][9][10].However, the complex nature of this question implies that most of such work is devoted to observations on particular traits within particular types or species, based on experimental observations, and sometimes phenomenological.In contrast, the results in this paper take the role of general mathematical theorems that hold for any biological system whose dynamics can be described by the models considered [Eqs.( 4), (7), or (13)].The techniques introduced here open a path to derive rate limits for modifications of these stochastic dynamical models, too.
It is worth emphasizing that the main results in this work, which constrain rates for dynamics with mutations [Eqs.(9)(10)(11)(12)] and which discriminate quantities that are resilient against stochastic driving forces [Eqs.(14,15)], are not implied by the general speed limits derived in Ref. [22] nor by recent works generalizing the fundamental theorem of natural selection [15,16].Making use of the structure of the differential equations that model biological evolutionary processes was crucial in deriving the new rate limits.In particular, speed limits on stochastic noisy dynamics such as Eq. ( 13) are mostly unexplored.It would be interesting to study applications of the methods used here to other models of evolutionary processes that account for genetic drift and finite size effects [54,[78][79][80], other stochastic drives [81,82], or to quantify the environment's effect on trait plasticity [83].At the same time, while I have argued and illustrated in Fig. 2 how the rate limits are useful in discriminating observable traits in terms of their evolutionary rates, it would be insightful to test these ideas in biological systems from experimental data.
The overarching aim of this work was to mathematically formalize the common knowledge that biological diversity influences evolutionary processes.The results shown here suggest and quantify ways in which variability, as measured by the standard deviations of fitness and of the biological traits of interest, can serve as a resource by allowing for faster evolution.It is tempting to ponder about the implications to related fields [84][85][86].
where γ j characterizes the strength of the stochastic driving forces for population j, and dW j are Wiener noises, which satisfy dW j dW k = δ jk dt and dW 2 j = dt.Let us focus on the noise-averaged change in an expectation value over a time τ relative to the change in ⟨A⟩ due to state changes from natural selection and mutations, as modeled by the replicator-mutator equation.It is given by where I used that δa j := a j − ⟨a⟩, that probability is conserved, j dp j = 0, and Eq.(C1).Here, In order to prevent confusion, I explicitly include the time-dependence of the noise terms.For Wiener processes, it holds that dW t j dW t ′ j = δ jk δ tt ′ , and that the noise terms are independent from all other functions (in the Itô picture), so that p t j dW t j = p t j dW t j = 0.The integrals that lack a line element dt, proportional to a noise term, correspond to Itô integrals [56].
Using the rules of Itô calculus described above and that j δa j p j = 0 gives that Combining Eq. ( C3) and (C5) results in This follows from γ j ≤ ∥γ∥ ∞ := max j {γ j } and p 2 j ≤ p j .This proves Eq. ( 14) in the main text.

Limits under other stochastic dynamics
An alternative stochastic replicator equation of the form was derived in Ref. [57].For simplicity, I assume no mutations in this sub-appendix, but it is easy to generalize the results that follow to account for them.Here, γ j characterizes the strength of the stochastic driving forces for population j, and dW j are Wiener noises, which satisfy dW j dW k = δ jk dt and dW 2 j = dt.Consider the noise-averaged change in an expectation value over a time τ , relative to the change d⟨A⟩ − ⟨dA⟩ replicator = cov(A, f )dt in ⟨A⟩ due state changes from natural selection as modeled by the replicator equation.
Following similar calculations as above, one finds that all cross terms proportional to a single noise term vanish upon averaging, and therefore Then, using the rules of Itô calculus described above and that j δa j p j = 0 gives that the second term is The second line follows from the Cauchy-Schwarz inequality and from γ j ≤ ∥γ∥ ∞ := max j {γ j }, and the third line holds because p 2 j ≤ p j .Note that the first term dominates the upper bound for long times τ ≫ ∥γ∥ 2 ∞ .A looser bound holds by using that ⟨γ 4 ⟩ ≤ ∥γ∥ FIG.1.Evolutionary rate limits.-Inthis work, I derive inequalities that bound the evolutionary rates of quantifiable traits.The rate limits apply to populations that evolve under mutations, natural selection, and random genetic drift.Mutations take a preponderant role in biology, enabling evolutionary changes that can lead to new types.Such new types can have different fitness than their ancestors.Natural selection is the process by which the frequencies of fitter types tend to increase relative to their peers in a given environment (in the illustration, turtles became fitter to survive in a drastically-changed environment.)A population's evolution can also be affected by chance, known as genetic drift (in the illustration, a random event reduced the population of an ant colony.)The results in this Letter imply that, under natural selection, a trait A with high variability σA can evolve faster than a trait B with small variability σB.At the same time, the trait with higher variability is more susceptible to the effect of random fluctuations in populations due to genetic drift.In this way, the rate limits can be used to discriminate quantifiable traits in terms of their maximum evolutionary rates and their responsiveness to different evolutionary forces.

FIG. 2 .
FIG. 2. Discriminating evolutionary regimes.-Therate limits derived in this work constrain the evolution of quantifiable traits of populations driven by natural selection, mutations, or genetic drift.The inequalities depend on (i) the expectation values {⟨f ⟩, ⟨A⟩, ⟨A⟩Π} and (ii) standard deviations {σ f , σA, σ Π A } of the trait of interest and of the fitness profile of the population, and (iii) on the maximum strength ∥γ∥∞ of genetic drift stochastic forces.Here, ⟨A⟩Π and σ Π A are evaluated in the distribution Πj defined in Eq. (8), which characterizes type's frequencies if the system only evolved due to mutations.Equation (15) can be used to mathematically delineate three distinct evolutionary regimes: ---[(a) -Strong selection regime] When σ Π A σ f ≫ ⟨f ⟩|⟨A⟩Π −⟨A⟩|, ∥γ∥∞σA , the contributions of mutations or genetic drift to the changes in a trait are negligible.This can be identified as the regime where natural selection is the dominant evolutionary force.In it, the maximum rate of a trait is constrained by its standard deviation σ Π A evaluated in the mutation-driven distribution Π: smaller trait variability implies smaller maximum rates.---[(b) -Strong mutation regime] Natural selection and genetic drift contributions are negligible when ⟨f ⟩|⟨A⟩Π − ⟨A⟩| ≫ σ Π A σ f , ∥γ∥∞σA .Then, the rate of change of a trait is determined by the change induced by the mutation-driven distribution, and Eq.(9) implies that d⟨A⟩/dt − ⟨ Ȧ⟩ ≈ ⟨f ⟩(⟨A⟩Π −⟨A⟩).---[(c)-Genetic drift regime] The variability of a trait determines its susceptibility to changes due to genetic drift.In the regime ∥γ∥∞σA ≫ σ Π A σ f , ⟨f ⟩|⟨A⟩Π −⟨A⟩| where the stochastic force dominates, the total change in a trait over a time τ is bounded by the integrated variance of the trait and the maximum drift strength.

γ 2 l p 2 l= δf j − γ 2 j p j + l γ 2 l p 2 l
j δa k (dp j − p j δf j dt) (dp k − p k δf k dt ′ ) j δa k (p j X j − p j δf j ) (p k X k − p k δf k )dtdt ′ j δa k p j p k Y j Y k , (C8)by using that dp j = p j X j dt + p j Y j , with the notationX j := f j − ⟨f ⟩ − γ 2 j p j + l (C9) Y j := γ j dW j − l γ l p l dW l .(C10)

δa j δa k p j p k Y j Y k = τ 0 τ 0 jk δa j δa k p j p k γ j γ k dW t j dW t ′ k − 2 τ 0 τ 0 jk δa j δa k p j p k γ j dW t j l p l γ l dW t ′ l + τ 0 τ 0 jk δa j δa k p j p k l p l γ l dW t l m p m γ m dW t ′ m = τ 0 j (δa j ) 2 γ 2 j p 2 j dt − 2 τ 0 jk δa j δa k p j p k γ j p j γ j dt + τ 0 jk δa j δa k p j p k l p 2 l γ 2 δa j δa k p j p k −γ 2 j p j + l γ 2 l p 2 l
δa j δa k (p j X j − p j δf j ) (p k X k − p k δf k )dtdt ′ that j p j δa j =0 in the third line.Combining the previous two equations gives 1 τ τ 0 d⟨A⟩ − ⟨dA⟩ − cov(A, f )dt

4 ∞ , 1 τ τ 0 d⟨A⟩j p j + l γ 2 l p 2 l
− ⟨dA⟩ − cov(A, f )dt 2 equation analogous to Eq. (14) in the main text for the stochastic replicator equation (C7).An analogous constraint holds when mutations are included,dp j = p j k p k f k Q kj − p j ⟨f ⟩ − γ 2 dt + p j γ j dW j − l γ l p l dW l , f k − ⟨f ⟩)Q kj a j dt