1 Introduction

figure a

Kelly’s work on optimal gambling strategies and the value of side information was arguably the first convincing attempt at applying concepts from information theory for analysis in a different field Kelly (1956). This work was the precursor to growth-optimal portfolio theory which has extended the basic ideas to the realm of capital markets (Cover and Thomas 2006). There has recently been a resurge of interest in employing insights from optimal gambling theory in models of adaptive bet-hedging under fluctuating environments, where close analogies between the economic and biological setting have been convincingly made apparent (Bergstrom 2014; Rivoire and Leibler 2011; Donaldson-Matasci et al. 2010).

Biological bet-hedging was originally proposed to explain the observation of un-germinated seeds of annual plants (Cohen 1966). This strategy involves the variable phenotypic expression of a single genotype, rather than a result of genetic polymorphism, although it is difficult to empirically determine whether observed phenotypic diversity in a population arises from randomization by identical genomes or from an underlying polymorphism (Seger and Brockmann 1987). Indeed, evolutionary biologists have long acknowledged that in a stochastically variable environment, natural selection is likely to favor a gene that randomizes its phenotypic expression (Bergstrom 2014). Recent work has revealed a variety of potential instances of bet-hedging populations: delayed germination in desert winter annual plants that meets postulated criteria of adaptive bet-hedging in a variable environment (Gremer and Venable 2014), bacterial persistence in the presence of antibiotics that appears to constitute an adaptation tuned to the distribution of environmental change (Kussell et al. 2005), flowering times in Lobelia inflata which point to flowering being a conservative bet-hedging strategy (Simons and Johnston 2003), or even bet-hedging as a behavioral phenotype, such as the case of nut hoarding in squirrel populations in anticipation of short or long winters (Bergstrom 2014).

Notwithstanding these empirical findings, identifying actual cases of adaptive bet-hedging in the wild remains elusive. As Seger and Brockmann (1987) have noted more than three decades ago, it is in general difficult to determine whether observed diversity of behavior in a population arises from randomization by genetically identical individuals or from genetic heterogeneity within co-located individuals optimized for different environmental conditions. Moreover, phenotypic heterogeneity can arise within genetically homogenous populations as a form of specialization in a stable environment through stochastic gene expression, positive feedback loops, or asymmetrical cell division, all processes where bet-hedging is not at play (Rubin and Doebeli 2017). These difficulties provide further impetus for constructing better and more elaborate models to test against the data.

Of particular note in classic bet-hedging models is the adoption from economic theory of asymptotic growth rate optimality as the target function for fitness maximization strategies, where growth in wealth is analogous to growth in lineage size. Indeed, since evolution proceeds by shifting gene frequencies over generations, with frequency changes being multiplicative, long-term fitness is commonly measured by geometric mean fitness across generations (Hopper 2018). At the same time, it is also widely acknowledged that long-run growth rate is not a valid measure of fitness under fluctuating environments, such as in the case of bet-hedging populations (Lande 2007).

The resulting intrinsic unpredictability has led some researchers to formulate a probabilistic perspective for natural selection that integrates various effects of uncertainty on natural selection (Yoshimura et al. 2009). The applicability of geometric mean fitness has also come into question under finite-population models, where the probability of fixation provides additional and sometimes more suitable information than the geometric mean fitness (Proulx and Day 2001), and in periodically cycling selection regimes, where evolutionary success depends on the length of the cycle and the strength of selection (Ram et al. 2018). Moreover, both gambling and bet-hedging models targeting optimal growth rate implicitly assume an infinite time horizon in formulating the geometric average, thereby ignoring the finiteness of actual horizons over which both economic and evolutionary processes ultimately act. The problem is further amplified when interim extinction risk is taken into account, especially under finite-population models. Lineage growth trajectories which are highly stochastic are at risk of large “drawdowns,” which may pull the population below some extinction threshold, despite possessing a high asymptotic growth rate. Here we aim to incorporate considerations of finite evolutionary horizons and extinction risk in the search for adaptive optimality in bet-hedging models.

1.1 Background: The Standard Model

Most adaptive bet-hedging models are largely based on the classic horse race gambling model associated with Kelly (1956), where the biological counter-part is a lineage apportioning bets on several possible environments. Assume that k horses run in a race, and let horse \(X_i\) win with probability \(p_i\). If horse \(X_i\) wins, the odds are \(o_i\) for 1. A gambler wishes to apportion his bankroll among the horses \(0<f_i\le 1\), such that \(\sum f_i=1\) and participate in indefinitely repeated races \(n\rightarrow \infty \). How to best apportion the bankroll each time? In this setting, wealth is a discrete-time stochastic process over n periods,

$$\begin{aligned} W_n=\prod _{i=1}^n W_i(X) \end{aligned}$$

where \(W(X)=f(X)O(X)\) is the random factor by which the gambler’s wealth is multiplied when horse X wins. More explicitly,

$$\begin{aligned} W_n (f)=\prod _{i=1}^k (f_i o_i )^{H_i}, \quad \text {where, } H\sim \mathrm{Multinomial}(n,k,[p_1,\ldots ,p_k]). \end{aligned}$$

Kelly’s first insight was that choosing to simply maximize expected wealth (for any time horizon n) gives \({{\,\mathrm{\mathrm{arg\,max}}\,}}_f E[W_n (f)]=1\), with the implication that one bets everything on a single horse (the one with the highest \(p_i\)) and a consequent chance of total ruin once that horse loses a race. Therefore, Kelly proposed maximizing the asymptotic growth rate [the rigorous justification provided by Breiman (1961)]. By the law of large numbers random wealth may be expressed as,

$$\begin{aligned} W_n (f) \doteq 2^{n E[\log W(X)]} \end{aligned}$$

where,

$$\begin{aligned} E[\log W(X)]=\sum _{i=1}^k p_i \log f_i o_i \end{aligned}$$

is the asymptotic exponential growth rate. If the gambler stakes his entire wealth each time, i.e., \(\sum f_i=1\), then

$$\begin{aligned} E[\log W(X)]=\sum _{i=1}^k p_i \log o_i -H(p)-D(p||f) \end{aligned}$$

is maximized (convex nonlinear optimization) at “proportional gambling” \(f=p\) where \(D(p\Vert f)\) is minimized, without regard the actual odds provided by the bookie.

Indeed, the notion of proportional gambling, made famous by Kelly’s treatment, has found its way into classic models of diversified bet-hedging. In such models often assumed that “appropriate phenotypes are produced in proportion to the likelihood of each environment” (Hopper 2018) and that consequently “the classical bet-hedging prediction [is] that the optimum probability for employing a strategy is approximately equal to the probability that the strategy will be useful” (King and Masel 2007). Here we follow recent approaches that extend the standard model to non-lethal environments via a full-fitness matrix, such that this notion is no longer directly applicable.

Breiman (1961) was first to show that the Kelly solution is optimal in two convincing ways: [a] that given a Kelly strategy \(\phi ^*\) and any other “essentially different” strategy \(\phi \) (not necessarily a fixed fractional betting strategy),

$$\begin{aligned} \lim _{n\rightarrow \infty } \frac{W_n (\phi ^*)}{W_n (\phi )}= \infty \quad a.s \end{aligned}$$

and [b] that it minimizes the expected time to reach asymptotically large wealth goals. Moreover, this strategy is myopic in the sense that at each iteration of the race one only needs to consider the presently given parameters (Hakansson 1971). However, Kelly strategies may also yield tremendous drawdowns a problem widely recognized in the gambling community, such that optimal Kelly is often viewed as “too risky”; in practice gamblers and investors use “fractional Kelly” which deviates from the optimal solution but reduces the effective variance of the stochastic growth (Fig. 1). In the biological framework, this can lead to abrupt extinction events in finite (especially small) populations with highly stochastic lineage growth trajectories. A further complication is that the underlying probability distributions are merely estimated from past data and model assumptions, leading often to over-betting and increased risk (MacLean et al. 2011).

Fig. 1
figure 1

The asymptotic exponential growth rate and its (finite-horizon) variance for either the two-horse racing model or the classic bet-hedging model (two environments). Note that the strategy (f on the x-axis) that maximizes the growth rate is far from the locus of minimal variance (Color figure online)

In this work, we extend the existing models to incorporate both interim extinction risk and finite evolutionary time horizons within a bet-hedging framework. This requires re-conceptualizing geometric mean fitness for such highly stochastic growth scenarios. We ultimately derive fitness functions that better account for such conditions where the fluctuating environment is strongly coupled to both long and short-term growth and locate optimal stable equilibria.

2 Methods

2.1 The Full-Fitness Matrix Model

We assume environments are i.i.d random events across generations, multinomially distributed (with some results generalized to non-identically distributed environments). Individuals within lineages have a static full-fitness matrix \([O_{ij}]\) in which non-lethal environments have low but generally nonzero fitness (Donaldson-Matasci et al. 2010; Rivoire and Leibler 2011). We adopt a finite-population model where lineages start off with some initial population size \(W_0\), implicitly assumed higher than some bet-hedging evolutionary threshold (King and Masel 2007). Lineages then evolve strategies to randomize individual phenotypes toward maximizing growth across finite horizons in the face of interim extinction threats. More formally, with k environments and phenotypes,

(1)

the general model of lineage growth trajectory across n generations under strategy f is a random process,

$$\begin{aligned} W_n=\prod _{i=1}^k \Big (\sum _{j=1}^k f_j o_{ij} \Big )^{H_i} \end{aligned}$$
(2)

where,

$$\begin{aligned} H\sim \mathrm{Multinomial}(n,k,[p_1,\ldots ,p_k ]) \end{aligned}$$

with off-diagonal values reflecting the lower fitness for non-matching environments,

$$\begin{aligned} o_{ii}>o_{ij} \ge 0 \text { and } o_{ii}>1 \end{aligned}$$

and where all individuals in a lineage are bet-hedging,

$$\begin{aligned} \sum _{i=1}^k f_i =1. \end{aligned}$$

And finally, using a straightforward formulation of the growth rate, \(W_n^{1/n}\), a random variable for any finite horizon.

We first derive the asymptotic growth rate optimal “Kelly” solution for this setting (\(f^{\mathrm{Kelly}}\)) with a corresponding bet-hedging region of the environment simplex (“Appendix A”). Relaxing the assumption of i.i.d environments, we derive the static Kelly solution for the case of nonstationary environments—where environments are independent but not identically distributed across generations (“Appendix B”). While under nonstationary environments an optimal growth rate is reached with a dynamic myopic strategy, we focus here on a static strategy since adaptations effectively stabilize across time spans much higher than single generations, such that from evolutionary considerations dynamic strategies are not likely to emerge. Alternative models of fluctuating environments such as Markov chains with underlying switching probabilities (e.g., Li et al. 2017) are not pursued here and left for future work. Finally, we identify a “reference” strategy that admits deterministic growth trajectories, namely the “Dutch book” solution (where the variance of the finite-time growth rate is zero) and characterize the consequent loss of growth incurred by exchanging opportunity for certainty (Appendix C).

2.2 Relative Fitness Payoff Function

We now wish to go beyond the standard approach of targeting the optimization of the asymptotic growth rate as undertaken in the previous section—to incorporate finite evolutionary horizons and extinction risk considerations. For the sake of simplicity, we confine our model here to the case of \(k=2\) environments and phenotypes (so that the two environments occur with probability p and \(1-p\)). To motivate the shift to a finite-horizon framework, we first highlight an important property of our stochastic growth model, known also in portfolio theory (Markowitz 2006). We prove that for any two essentially different strategies, the maximal time \(n_0\) one lineage “dominates” the other is finite for every realization of lineage trajectory pair (“Appendix D”). The exponentially diminishing histogram of last intersection times of given two growth strategies in Fig. 2b (with a single instance of two trajectories for illustration in Fig. 2a) demonstrates this phenomenon.

Fig. 2
figure 2

Two perspectives on growth trajectories \(W_n\) and their last “intersection” point for two lineages with different strategies (a log-optimal Kelly and a suboptimal “fractional Kelly”). a The optimal Kelly strategy eventually departs from a suboptimal strategy without further intersections, “the last intersection.” | b A histogram demonstrating that the last intersection of the trajectories of any two growth strategies occurs at some finite time (here, somewhat above 500 generations) (Color figure online)

The sustained variance and high skewness of the growth rate distribution under any finite horizon necessitates a comparative approach in formulating a fitness payoff function (in fact, the growth rate is asymptotically log-normal as shown in “Appendix E”). Consider a relative fitness measure for two different lineage strategies f and g: The probability that a random trajectory of a lineage with strategy f exceeds the random trajectory of a lineage with strategy g (given time horizon n),

$$\begin{aligned} h(f,g)=P(W_n (f)>W_n (g)) \end{aligned}$$
(3)

with an induced relation defined by,

$$\begin{aligned} W_n (f)\ge W_n (g): P(W_n (f)>W_n (g)) \ge P(W_n (f)<W_n (g)). \end{aligned}$$
(4)

We may interpret this probabilistic relation between two strategies as relative fitness. Note that since realizations of \(W_n (f)\) and \(W_n (g)\) stem from the same underlying stochastic environmental sequence, they will generally be highly correlated (with the corresponding logarithmic growth rates in fact perfectly correlated, as shown in “Appendix F”). Consequently, the probability in Eq. (3) must be derived from their joint distribution rather than simply from marginal distributions. Figure 3 depicts realizations of the log growth rates of \(W_n (f)\) and \(W_n (g)\) as histogram distributions for some choice of strategies f and g, and some finite evolutionary horizon n. Asymptotically with time horizon n, such distributions approach normality with variance going to zero (“Appendix E”).

Fig. 3
figure 3

The distributions of log growth rates for two competing strategies here illustrated as histograms (the x-axis is the log growth rate). In this simple case shown, the strategy generating the higher growth rates right distribution (blue) is evidently of higher fitness than the one generating the left distribution (green). For clarity, only the marginal distributions are shown (the complete picture is in the joint distribution, due to the correlation between any two growth trajectories that share the same fluctuating environment regime) (Color figure online)

A few properties of the order induced by this relation are worth highlighting. [a] it is a complete order since any two \(W_n\) are comparable under the relation, [b] it is transitive for any n and consequently a pre-order, and [c] its maximal element is \(W_n^* (f^{\mathrm{Kelly}})\), such that both the order induced by \(E[\log W_n (f)]\) and the order induced by the payoff \(P(W_n (f)>W_n (g))\) form complete preorders and have the same maximal element (“Appendix G”). Despite these beneficial properties, given any “vanilla” strategy g and time horizon n, the strategy that maximizes the payoff function,

$$\begin{aligned} {{\,\mathrm{\mathrm{arg\,max}}\,}}_f P(W_n (f)>W_n (g)) \end{aligned}$$

will vary as a function of g and n (demonstrated by counterexamples), and in particular will not necessarily be \(f^{\mathrm{Kelly}}\). This implies that a wildtype lineage with strategy g different from \(f^{\mathrm{Kelly}}\) will eventually be overtaken by some mutant invasive lineage with a strategy that maximizes this payoff function, a process that may potentially remain in recurrent flux, with invasive lineages replacing a wildtype lineage.

2.3 Competitive Optimality with Risk

To see whether evolutionary stable optima may also emerge, we develop a game-theoretic approach. Players are lineages with particular bet-hedging strategies and random initial population size. Lineages interact by competing over a common niche subject to the same environmental fluctuations. This setup is in some contrast to more standard evolutionary game theory settings, where agents are organisms rather than lineages and where the notion of an iterated strategy is prominent, but maintains the central aspect of interactions formalized in a payoff function (e.g., Stollmeier and Nagler 2018). A lineage survives the competitive encounter by avoiding extinction (defined in what follows) while exceeding its opponent in size over a given time horizon. This outcome is determined by a game-theoretic deterministic payoff function, modified from Eq. (3) to incorporate an extinction threshold and randomized initial lineage size. Ultimately, we are searching for Nash equilibria.

This approach is motivated by the classic work on time-invariant game-theoretic competitive optimality, within the scope of growth-optimal portfolio theory (Bell and Cover 1980, 1988). Bell and Cover consider a competitive setting for a stock portfolio model under any finite number of investment periods and prove that for any relative wealth payoff \(E[\phi (UW_1/VW_2)]\) and portfolio wealth \(W_1\) and \(W_2\), there are conditions on the function \(\phi \) such that the log-optimal Kelly portfolio is a solution to the game, given initial randomizations U and V (independent and of equal expectation). In particular, \(\phi (x)=\chi _{[1,\infty )}(x)\) results in the payoff \(\mathbb {P}(UW_1\le VW_2)\) with the log-optimal portfolio as a game-theoretic solution, given some initial fair randomizations. This additional fair randomization reduces the effect of small differences in end wealth, thus avoiding unwanted cases where the optimal strategy is beat by a small amount most of the time (Cover and Thomas 2006).

2.4 The Payoff Function in a Game-Theoretic Setting

For any time horizon n and extinction threshold d, we define a (deterministic) payoff function: The probability that a random trajectory of a lineage with strategy f exceeds the random trajectory of a lineage with strategy g without first going extinct (given time horizon n),

$$\begin{aligned}&M_n (f,g)=P(u_0 W_n (f)>v_0 W_n (g) | \text { extinction level } d) :=\nonumber \\&\quad P(u_0 W_n (f)>v_0 W_n (g) \wedge W_i,V_i>d, i:1,\ldots ,n)+\nonumber \\&\quad P(u_0 W_i (f)>d, i:1,\ldots ,n \wedge v_0 W_i (g)\le d \text { some } i) \end{aligned}$$
(5)

with initial population size independent randomizations \(u_0\) and \(v_0\), independent and of same mean but possibly of a different distribution class.

This payoff function induces a symmetric discrete-valued non-constant-sum game setting, although it is conceptually “zero-sum” \(M_n(f,g)+M_n(g,f)<1\) (“Appendix H”). Crucially, our payoff matrix is finite since it reflects the finitely many strategies possible in a finite-population model—there can only be N different sized partitions of a population of size N in betting on two environments (under \(k=2\) environments and phenotypes). A low-resolution toy-model instance of the payoff matrix is depicted in Fig. 4.

Fig. 4
figure 4

An example of a (\(2 \times 2\) strategies) evolutionary payoff matrix for a game of two lineages, with primary lineage payoff in bold (Color figure online)

Fig. 5
figure 5

An example of payoff matrix simulations (a and c darker cells represent higher probabilities) and the resulting Nash equilibria in the maximal element matrices (b and d blue cells for column-maximal primary player payoffs, red for row-maximal opponent payoffs and yellow for symmetric Nash equilibria). Only portions of the matrices around the equilibrium are shown. a+b Low-resolution matrices corresponding to a small population with limited strategies | c+d Higher resolution matrices, which correspond to a larger population and subsequently higher number of strategies (multiple maxima in some adjacent cells of panel D is an effect due to a combination of using finite runs in the simulation of the payoff function, such that computed probabilities are rational values, along with high resolution in the range of strategies). The model uses a \(2\times 2\) fitness matrix: \([o_{11}=3.0;o_{12}=0.2;o_{22}=1.8;o_{21}=0.1]\) with \(p=0.594\) and a resulting log-optimal strategy \(f=0.60\)) (Color figure online)

Our goal would be to identify pure strategy Nash equilibria reflecting the evolutionary solutions to competitive bet-hedging. In particular, we would like to explore the conditions under which a bet-hedging setting admits a symmetric equilibrium and whether it is unique. In “Appendix I,” we prove that for an infinite-size payoff matrix (i.e., continuous strategies) the log-optimal strategy is the solution to this game, invariant with the choice of time horizon. Moreover, any finite matrix representing the N strategies possible for a lineage of finite size N necessarily also admits a solution, as illustrated in Fig. 5. This solution is the strategy closest to the log-optimal strategy under the finite resolution framework, such that it converges to it asymptotically with N (“Appendix L”). Finally, under a nonstationary environment model the log-optimal strategy again emerges as the equilibrium static strategy—even given short time horizons (“Appendix M”).

The effect of lineage-size extinction thresholds on actual rates of extinction of random growth trajectories is illustrated in Fig. 6a. As would be expected, higher thresholds of extinction correspond to higher probabilities of extinction, with extinction rates that converge quickly to asymptotic values (“Appendix N”). Numerical simulations indicate that when incorporating low extinction thresholds that result in low extinction rates, the symmetric Nash equilibrium remains stable at the log-optimal strategy. Higher thresholds may result in a number of scenarios: a shift of the symmetric equilibrium away from the log-optimal solution, complete lack of equilibrium solution, or the emergence of multiple symmetric equilibria; in conjunction, multiple pairs of off-diagonal equilibria may appear (see Fig. 6b for one such scenario).

Fig. 6
figure 6

a The accumulative probability of extinction for log-optimal (Kelly) strategies for different lineage-size extinction thresholds. | b An instance of the maximal element matrix of a payoff matrix (a portion around the solutions) resulting from introducing an extinction threshold: two off-diagonal Nash equilibria (black) along with a symmetric Nash equilibrium shifted from the log-optimal strategy (yellow). The simulations use a \(2 \times 2\) fitness matrix: \([o_{11}=2.5;o_{12}=0.2;o_{22}=1.4;o_{21}=0.1]\) with \(p=0.75\), and with \(n=60\) in the simulation of panel B (Color figure online)

2.5 Minimum Time to Reach a Population Threshold Size

To gain further perspective on optimal strategies under highly stochastic growth, we consider evolutionary competition between lineages, where survival is determined by reaching a certain threshold of lineage size in minimal time [e.g., for \(K-\)selected species, see Reznick et al. (2002)]. In effect, the lineage with growth characteristics that minimize the time to reach a certain population size threshold “wins,” a setting with potential relevance in the context of competitively colonizing a limited niche, as in range expansion scenarios [see Villa Martin et al. (2019) for a bet-hedging population expanding into an unoccupied space]. We follow the classic results of Breiman 1961 on the log-optimal portfolio as the optimal strategy minimizing the expected time to reach an asymptotic target wealth, but instead of an infinite target we base the fitness payoff function on finite targets. Initial insight into the effect of strategy choice on the consequent distributions of minimal time (Fig. 7a) is provided by comparing their expectation, where the optimality of Kelly is already apparent (Fig. 7b).

Fig. 7
figure 7

a An example of the marginal distributions of the minimal time to reach a certain population size threshold, for two strategies | b A simulation of \((E[T(g,c)]-E[T(f^*,c)])\) for several target values c, under the range of vanilla strategies g corresponding to positive asymptotic growth rates, where \(f^*=0.8\) (Color figure online)

Instead of considering expectations of (highly correlated) minimal time distributions, we devise a more informative fitness payoff function based on the joint distribution. Crucially, this payoff will naturally be amenable to a game-theoretic approach, in line with the type of analysis in the previous section with payoff \(M_n (f,g)\). As before, we condition the probability on avoiding an extinction threshold. The payoff captures the probability that a trajectory following strategy f reaches threshold c before a trajectory following strategy g, conditioned on avoiding an extinction threshold d. If both trajectories reach c at the same time (since time is in discrete generations), then the one which overshoots with a greater margin above c ‘wins’. Denote by T(fc) the minimal time distribution given strategy f and target lineage size c,

$$\begin{aligned} M_c (f,g)=P(T(f,c)<T(g,c) | \text { extinction level } d). \end{aligned}$$

More precisely, we denote new trajectories \(\{W^E_k\}_{k=1}^n\) by

$$\begin{aligned} W_0^E \in U[a,b], \end{aligned}$$

and for all \(k=0,\ldots , n-1\)

$$\begin{aligned} \quad W^E_{k+1} = {\left\{ \begin{array}{ll}W^E_{k} \overline{o}_1(f)^{x_{k+1}} \overline{o}_2(f)^{1-x_{k+1}}, \text { if } W^E_{k} \ge E\\ 0, \text { if } W^E_{k} < E. \end{array}\right. } \end{aligned}$$

We denote also by

$$\begin{aligned} T(f,c) := \min \{n: W^E_n(f) \ge c\} \end{aligned}$$

the first time when the trajectory \(\{W^E_k\}_{k=0}^n\) cut the threshold c. \(T(f,C) = \infty \) if and only if this trajectory does not cut the threshold.

Then the payoff matrix \(M_c(f,g)\) is defined by

$$\begin{aligned} M_c(f,g) :=&\mathbb {P}(T(f,C) < T(g,C) | \text {given that at least one of them is finite})\nonumber \\&+ \mathbb {P}(T(f,C) = T(g,C), W^E_{T(f,C)} > V^E_{T(g,C)} | \nonumber \\&\text {given that at least one of them is finite}). \end{aligned}$$
(6)

We then identify pure strategy Nash equilibria reflecting the evolutionary solutions with the new relative payoff \(M_c (f,g)\). In “Appendix J,” we prove that again Kelly is the solution to the game, invariant to the evolutionary “choice” of target population size c and that under a nonstationary environment regime Kelly emerges as the static equilibrium strategy. Finally, we highlight a deep mathematical link of this probabilistic perspective for minimal time optimality to the competitive optimality setting with payoff \(M_n (f,g)\). Formally, \(M_c(f,g)\) can be rewritten as a convex linear combination of \(M_n(f,g)\): \(M_c(f,g) = \sum _{n=0}^{\infty } P( W_0 W_n(f) > V_0 W_n(g), T(f,c) = n)\) (see “Appendix J” for more details).

3 Discussion

In this work, we provide further support for the robustness of the expected log criterion as an optimality solution for biological bet-hedging. We develop a game-theoretic framework inherently invariant to the span of evolutionary horizons while incorporating considerations of interim extinction risk and use multiple optimality criteria to strengthen our results. This approach goes beyond standard models of bet-hedging, which focus on indefinite “long-term” growth rates and that ignore accounting for interim risk. Previous work generally upholds that “phenotypes with the greatest long-term average growth rate will dominate the entire population” as “the basic principle” used in optimization (Yoshimura and Jansen 1996), or that a proxy for the likely outcome of evolution is “to think of organisms as maximizing the long-term growth rate of their lineage” (Donaldson-Matasci et al. 2010).

Nevertheless, some authors have recently acknowledged the importance of accounting for finite time horizons. For instance, Rivoire and Leibler (2011) note in passing that in their model “the growth rate emerges as a unique measure of fitness when considering the long-term limit \(T\rightarrow \infty \), but, if considering a finite ‘horizon’ there may be a different strategy that outperforms [it].” Indeed, as some evolutionists have argued, short-term fitness measures are also needed to achieve a full understanding of how evolution works in variable environments, as geometric mean fitness concerns the long-run evolutionary outcome (Okasha 2018). Moreover, long-term fitness metrics are typically formulated without regard to transient short-term population dynamics, in which lineages might come close to extinction. Under more inclusive models with extinction, selection in a fluctuating environment can also favor bet-hedging strategies that ultimately increase the risk of extinction (Libby and Ratcliff 2019). Given such considerations, the benefit of explicitly incorporating extinction considerations in stochastic growth models is clearly evident.

We have opted to focus on symmetric Nash equilibria rather than evolutionary stable strategies (ESS), which are strategies that cannot be beaten if the fraction of the rival invading mutants in the population is sufficiently small and are generally invoked in settings with iterative match-ups between individuals rather than lineages (Smith and Price 1973). Since the payoff in our game-theoretic setting pits one lineage against another (two different strategies), there is no explicit sense of invading mutants [but see Olofsson et al. (2009) for an ESS approach to bet-hedging]. Moreover, some of the classic aspects of Nash’s theorem do not directly apply within our setting. The theorem states that for every two-person zero-sum game with finitely many strategies, there exists a mixed strategy that solves the game (Nash 1951). While our framework is indeed “two-person”, it is not zero-sum and has finitely many strategies. Crucially, since an implicit goal of theoretical work such as ours may be toward predicting which strategies are likely to evolve, we focus on pure strategies rather than mixed ones, where the uniqueness of the equilibrium solution emerges as especially beneficial (echoing the classic approach of growth rate log-optimality where there is always a unique solution due to convexity).

We are not the first to attempt to model the expected minimal time to reach a finite asymptotic target, an extension of the seminal result of Breiman (1961) on properties of the log-optimal portfolio. Aucamp (1977) derived the first such analysis, given some basic assumptions that concern reaching a wealth target exactly vs. “overshooting” it. More recently, Kardaras and Platen (2010) find that in a continuous time or asset price model where a finite target can be exactly reached with no overshooting, the Kelly solution is still optimal; in a discrete time model Kelly is only approximately optimal, but if “time rebates” are introduced (to compensate overshooting the goal in the last investment period) it becomes exactly optimal. While these results on the expectation of the time distribution are in line with our analysis of stochastic lineage growth optimality, we obtain an even stronger result: Given finite population size targets, the log-optimal strategy emerges as a Nash equilibrium under a payoff function based on the joint distribution of minimal time trajectories.

Interestingly, Kelly (1956) has anticipated the application of his ideas in biological bet-hedging, writing “Although the model adopted here is drawn from the real-life situation of gambling, it is possible that it could apply to certain other economic situations...the essential requirements for the validity of the theory are the possibility of reinvestment of profits and the ability to control or vary the amount of money invested or bet in different categories.” It does not require a leap of the imagination to notice analogies of “economic situations” to evolutionary strategies, of “reinvestment of profits” to biological reproduction and growth, and of the “control” of invested money to evolved adaptive optimality. Of course, it is best appreciated with Shannon’s famous “bandwagon” warning in mind, cautioning over hasty attempts to apply insights from information theory to other fields (Shannon 1956).

3.1 Other Approaches to Optimization Under Finite Horizon and Risk

A seemingly straightforward way of introducing finite (albeit still arbitrary) horizons into optimization settings is by considering the expectation of a finite-horizon growth rate. This is the approach adopted in some recent stock portfolio models for finite horizons (Vince and Zhu 2013; Morgan 2015). Within our formalism from Eq. (2), this amounts to finding,

$$\begin{aligned} {{\,\mathrm{\mathrm{arg\,max}}\,}}_f E\Big [W_n(f)^{\frac{1}{n}}\Big ] = {{\,\mathrm{\mathrm{arg\,max}}\,}}_f \Bigg (\sum _{i=1}^k \Big (\sum _{j=1}^k f_j o_{ij} \Big )^{\frac{1}{n}} p_i \Bigg )^n. \end{aligned}$$

However, this implicitly assumes some arbitrary utility function, in this case the n-th root, the maximization of which requiring some justification. In contrast, Kelly’s focus on \({{\,\mathrm{\mathrm{arg\,max}}\,}}_f E[\log W_n]\) while implicitly assumes logarithmic utility, is equivalent the limit of the above expression, and leads to desired optimality properties as famously laid out by Breiman (1961).

A more convincing approach to maximizing wealth with risk management over finite horizons was proposed in Rujeerapaiboon et al. (2015) for portfolio construction. The authors consider the optimization of a minimum bound for finite-horizon growth,

$$\begin{aligned} {{\,\mathrm{\mathrm{arg\,max}}\,}}_f \Big \{{{\,\mathrm{\mathrm{arg\,max}}\,}}_c \mathbb {P}\Big (\frac{1}{n} \log W_n\ge c\Big )\ge 1-\varepsilon \Big \} \end{aligned}$$

with a degree of freedom corresponding roughly to a risk aversion or a choice of certainty parameter.

The expression above allows deriving the portfolio giving the highest minimum bound for wealth for any level of certainty \(\varepsilon \). While choosing a particular horizon n and a risk aversion parameter is perfectly sensible in an investment setting, the translation to the biological framework is problematic: what would be evolution’s risk aversion in this setting? Or the appropriate time horizon for optimization? Any choice of these two parameters would inescapably be arbitrary in nature. In an alternative approach Rujeerapaiboon et al. (2018) reformulate the Kelly gambling setting in terms of the Conservative Expected Value (CEV), a risk-averse expectation for highly skewed distributions. This amounts essentially to devising a systematic way of constructing fractional Kelly strategies such that it is strongly coupled with the infimum of the finite-horizon growth rate. Here again, there is an implicit arbitrariness in the choice of horizon length if applied in the context of an evolutionary framework, which we seek to avoid.

Other authors have focused on incorporating risk to the standard Kelly gambling setting with an infinite time horizon. For instance, Busseti et al. (2016) develop a systematic way to trade-off growth rate and drawdown risk by formulating a risk-constrained Kelly gambling problem within the standard setting of growth rate maximization under asymptotic horizons. The additional risk constraint limits the probability of a drawdown to a specified level. Nevertheless, for our purposes, percentage drawdown is arguably not a natural metric for representing lineage extinction risks, as compared with explicit extinction thresholds, especially in scenarios of competing finite-size populations (Ashby et al. 2017). Still other approaches may seek to target risk minimization as a primary criterion. In an approach akin to our Dutch book analysis, Wolf et al. (2005) minimize the growth rate variance and consequently the probability of extinction due to “unlucky” environmental trajectories. However, this is at the inevitable expense of achieving high stochastic growth rates, a vital aspect of evolutionary fitness.

3.2 Game-Theoretic Competitive Optimality of Bell and Cover

The results presented here can also be seen as both a special case and an extension of the classic results of Bell and Cover (1980, 1988). There are several important distinctions: [a] their setting is formulated for continuous random variables whereas our environments are discrete events, [b] their payoff implies a zero-sum game whereas our game is nonzero-sum (more accurately, non-constant-sum) due to the effect of extinctions, and [c] their payoff function is a straightforward probability while our payoff is effectively a conditional probability (includes considerations of extinction risk). Moreover, implicit in Bell and Cover’s setting is an infinitely sized payoff matrix, whereas our payoff matrix is finite since it reflects a finite number of strategies possible in a finite population. These distinctions have enabled us to show that, at least given the particular payoff function and discrete framework, the emerging symmetric Nash equilibrium is in fact a strict and unique one.

Some authors have generalized or utilized other aspects of the classic competitive optimality results. Most recently, Garivaltis (2019) has shown that discrete-time results of Bell and Cover (1988) hold equally well for continuous-time rebalanced portfolios in a competitive setting between two investors, each aiming to maximize the expected ratio of one’s own wealth to the other. In an original use of evolutionary ideas in finance, Lo et al. (2017) and Orr (2017) consider a payoff function capturing relative wealth of two competing investors each with some set initial wealth, focusing on finite-period analysis. They analyze optimal strategies of a primary player against a given “vanilla” strategy, a framework consistent with our initial relative payoff non-game-theoretic setting. They find that the particular vanilla strategy chosen plays an important role in the optimal allocation, in conjunction with initial wealth of both players.

Finally, our game-theoretic analysis may hint at a solution to a “coincidence” pointed out in Bell and Cover (1980). They were left perplexed as to why competitive optimality for a finite horizon turned out, by “coincidence,” to have the same solution (namely, Kelly) as in the growth-optimal portfolio: “Finally, it is tantalizing that \(b^*\) arises as the solution to such dissimilar problems [...] The underlying for this coincidence will be investigated.” Their follow-up 1988 paper suggests a “possible reason for the robustness of log-optimal portfolios” or why “log-optimal portfolios behave well in the competitive investment game”: namely that the wealth generated from any portfolio is always within “fair reach” of the wealth from the log-optimal portfolio. Indeed, the Kuhn–Tucker conditions and the consequent bound on the wealth ratio (Cover and Thomas 2006, Theorem 16.2.2) already imply that game-theoretic optimality is the driving force behind the asymptotic dominance. Fair randomization of initial wealth then leads to the game-theoretic solution for any increasing function of the wealth ratio. Our investigation of the payoff matrix suggests another perspective to this “coincidence.” Asymptotically with horizon n, the payoff matrix becomes maximally “contrasted,” with off-diagonal cells converging to probabilities of 0 or 1 (except those on “fault lines”), such that the Nash equilibrium emerges naturally. In effect, the “saddle-point” equilibrium, which has been established as invariant with n, asymptotically attains maximum curvature (“Appendix K”).

4 Conclusion

In this work, we have argued that under fluctuating environments and trait randomization geometric mean fitness should also encompass considerations of stochastic growth and extinction risk under finite evolutionary horizons. We show that for both the relative maximal growth payoff and the relative minimal time payoff there is a unique pure strategy symmetric equilibrium, which is invariant with evolutionary time horizon and robust to low extinction risk. Coinciding with the classic bet-hedging modeling approach, this is the Kelly log-optimal strategy. With higher thresholds of extinction, the equilibrium may shift away from Kelly and possibly branch out to multiple equilibria. Future work will be required to generalize the model to competitive optimality payoffs beyond pairwise lineages, Markovian environmental sequential transitions, random fitness matrices, and to more precisely capture the effect of high extinction thresholds on the optimal evolutionary solutions.