A scenario for an evolutionary selection of ageing

eLife assessment

Through a theoretical approach, this study makes important contributions to our understanding of the evolutionary causes of the ageing process. Using a simple individual-based model and computational simulations, the authors provide convincing evidence that ageing can be a trait under natural selection, opening the door for further discussion in the context of lifespan extension research.

https://doi.org/10.7554/eLife.92914.3.sa0

Significance of the findings:

Important: Findings that have theoretical or practical implications beyond a single subfield

Landmark
Fundamental
Important
Valuable
Useful

Strength of evidence:

Convincing: Appropriate and validated methodology in line with current state-of-the-art

Exceptional
Compelling
Convincing
Solid
Incomplete
Inadequate

During the peer-review process the editor and reviewers write an eLife Assessment that summarises the significance of the findings reported in the article (on a scale ranging from landmark to useful) and the strength of the evidence (on a scale ranging from exceptional to inadequate). Learn more about eLife Assessments

Abstract
eLife digest
Introduction
Results
Discussion
Materials and methods
Appendix 1
Data availability
References
Article and author information
Metrics

Abstract

Signs of ageing become apparent only late in life, after organismal development is finalized. Ageing, most notably, decreases an individual’s fitness. As such, it is most commonly perceived as a non-adaptive force of evolution and considered a by-product of natural selection. Building upon the evolutionarily conserved age-related Smurf phenotype, we propose a simple mathematical life-history trait model in which an organism is characterized by two core abilities: reproduction and homeostasis. Through the simulation of this model, we observe (1) the convergence of fertility’s end with the onset of senescence, (2) the relative success of ageing populations, as compared to non-ageing populations, and (3) the enhanced evolvability (i.e. the generation of genetic variability) of ageing populations. In addition, we formally demonstrate the mathematical convergence observed in (1). We thus theorize that mechanisms that link the timing of fertility and ageing have been selected and fixed over evolutionary history, which, in turn, explains why ageing populations are more evolvable and therefore more successful. Broadly speaking, our work suggests that ageing is an adaptive force of evolution.

eLife digest

It is a question as old as Darwin’s theory of evolution itself: how is ageing affected by natural selection? The prevailing view is that the process of biological ageing is not adaptive and therefore not directly subject to selection pressures. Take for example a gene causing a fatal disease late after an average individual had reproduced, thus being passed on to the next generation despite its detriment to the individual. This suggests that natural selection acts less strongly on such genes, which can therefore accumulate and cause aging if they do not impact an organism’s reproductive fitness earlier in life.

However, many studies have shown that specific genes control an animal’s lifespan and the onset of ageing through evolutionarily conserved mechanisms. For example, in fruit flies, aging can be categorised into two distinct phases determined by the manifestation of the so-called Smurf phenotype associated with accelerated signs of ageing and an increased risk of death. A pattern where the offspring of older parents live less long than those of younger parents has also been observed across species, also known as the Lansing effect. In this case, ageing can affect the reproductive success of future generations and can therefore be subject to selection pressures.

Roget et al. looked at the trade-offs between an individual’s reproduction and homeostasis using a mathematical model to address whether the distinct phases of aging – as seen in the Smurf phenotype – can appear and be maintained throughout evolution. Using a mathematical model, Roget et al. simulated individuals possessing only one copy of two genes. One controls the duration of reproductive ability, and the other defines the age at which the risk of death becomes non-zero.

This revealed that a simple hypothetical haploid and asexually reproducing system can evolve a life history separated into two phases in the computer simulations. Interestingly, the modelled organisms evolved in a way that the duration of reproduction exceeded the homeostatic maintenance duration. This generated a phase where individuals are capable of reproduction with a high risk of death, similar to the previously described Smurf phase. Roget et al. observed that aging populations showed a lower risk of extinction than non-aging ones, as well as an increased genetic variability of the offspring.

The apparent benefits of ageing in this model imply that ageing can be an adaptive force of evolution and subject to positive selection or, at least less negative selection than expected. This minimal model helps explain trade-offs between reproduction and homeostatic maintenance during evolution. Further work may include parameters such as sexual reproduction and multiple gene copies.

Introduction

Ageing is the umbrella term used to describe the processes that take place when an organism’s capacity to thrive diminishes with time. Patterns of ageing vary greatly given the organism, from negligible senescence to post-reproductive death through progressive age-dependent mortality increase (Jones et al., 2014). While ageing, as an observable process, is evident, the evolutionary role of ageing is unclear and conceptually challenging. An ageing individual is less fit, nevertheless, ageing seems to be broadly present through evolutionary time. Our work aims to explore the question, "is this mere chance– is ageing strictly a by-product of other things under selection– or is it somehow adaptive”? Soon after Charles Darwin published his theory of evolution, August Weismann situated ageing within this framework (Weismann, 1882) by theorizing that, “there exists a specific death-mechanism designed by natural selection to eliminate the old, and therefore worn-out members of a population” (Gavrilov and Gavrilova, 2002). Since then, however, it is mostly accepted that “ageing is not adaptive since it reduces reproductive potential” (Kirkwood and Holliday, 1979) and hence, fitness. Weismann’s own theories eventually evolved to more closely represent this current position.

At present, ageing is typically viewed as a ‘side-effect’, or byproduct, of other processes under selection (Fabian, 2011), which implies that ageing, or the mechanisms that cause ageing, are neither selected nor adaptive— precisely as capacities that would prove advantageous for a given population. This view took precedent starting in the 1950s and it is now assumed that the genetics or molecular processes that drive ageing help to explain how ageing has evolved (Gavrilov and Gavrilova, 2002). Peter Medawar’s theory of mutation accumulation defends that ageing is caused by the progressive accumulation of deleterious mutations with effects that show only late in life (Medawar, 1953). Williams’ antagonistic pleiotropy theory goes further than Medawar’s by presupposing the existence of antagonistic genes and mutations: beneficial at an early age, these genes/mutations prove disadvantageous at a later age (Williams, 1957). Evolutionary conserved genes involved in both the regulation of longevity and organismal growth were discovered in the model organism Caenorhabditis elegans (Kenyon et al., 1993) and later shown to be conserved in flies (Clancy et al., 2001), mice (Blüher et al., 2003), and humans (van Heemst et al., 2005). Thus, genetic modulators for longevity exist and express themselves through evolutionarily conserved physiological mechanisms. With genes involved in the onset of longevity, there is a potential substrate for selective pressure to apply. Regardless, it is generally accepted that ageing is neither a programmed nor beneficial trait for species (Kowald and Kirkwood, 2016).

The Smurf phenotype is a simple age-associated intestinal permeability phenotype that was first observed in Drosophila (Rera et al., 2011). Evolutionarily conserved in nematodes, zebrafish (Dambroise et al., 2016). and mice (Cansell et al., 2023), this phenotype allows for the identification of two distinct subpopulations– non-Smurf individuals and Smurf ones– at any time in a given population. All individuals undergo the transition (from non-Smurf to Smurf) prior to death (Rera et al., 2012; Tricoire and Rera, 2015). In flies, the Smurf phase is characterized by multiple physiological marks of ageing such as the high risk of impending death, loss of energy stores, systemic inflammation, reduced motility (Rera et al., 2012), and reduced fertility (Rera et al., 2018). More generally, the transcriptional hallmarks (Frenk and Houseley, 2018) usually associated with ageing are mostly observed in the latter phase (Zane et al., 2023). To summarize, this phenotype allows for the identification of two successive and necessary phases of life with all the age-related changes occurring in the last. Motivated by these biological observations, we recently assessed (Méléard et al., 2019) the possibility of obtaining, over time, such two phases of life. In order to simplify, we decided to consider the evolution of such a process in a bacteria-like organism, through the design and implementation of an asexual and haploid age-structured population mathematical model. We constrained the evolutive trajectory of ageing (within this model) through the Lansing effect– a transgenerational effect impacting longevity. Smurf individuals carry the propensity to demonstrate this effect. The Lansing effect is a transgenerational phenomenon, first described by Albert Lansing in the late 1940s, whereby the ‘progeny of old parents do not live as long as those of young parents’ (Lansing, 1954; Lansing, 1947). This was first observed in rotifers. More recently, it has been shown that older Drosophila females, and to some extent males tend to produce shorter lived offspring (Priest et al., 2002). Older zebra finch males give birth to offspring with shorter telomere lengths and reduced lifespans (Noguera et al., 2018). In humans, ‘older father’s children have lower evolutionary fitness across four centuries and in four populations’ (Arslan et al., 2017). Despite the absence of consensus regarding any underlying mechanism, the Lansing effect is broadly conserved and therefore relevant (Monaghan et al., 2020). We observed, through this Lansing-positive model, that the ageing phase overlaps with the pre-ageing phase in evolutionary time.

Here, we decided to generalize this model to any system able to reproduce and maintain homeostasis, without the necessary constraint of the Lansing effect, and in hopes of understanding how such a two-phase ageing process might have evolved. We thus show the following:

Through time, the end of the healthy phase and the beginning of the senescent phase converge even in the absence of a transgenerational effect (the Lansing effect).
With an equal Malthusian parameter at t₀, Lansing populations are more successful than non-Lansing populations, suggesting that the individual loss of fitness is compensated at the population level.
Ageing (or senescence-carrying) populations are more evolvable than non-ageing populations. We theorize this is because ageing populations are quicker to explore genotypic space.

This is all to suggest that ageing is, as a function, decreasing both reproductive and homeostatic capabilities of an organism, both an attractor configuration and an adaptive force of evolution, in opposition to what is most commonly assumed.

Results

The model (called bd model) and its population dynamics follows those described in Méléard et al., 2019. Briefly, the model delineates an asexual and haploid population, structured by a life-history trait that is defined by a pair of parameters - genes - (x_b, x_d) where x_b defines the fertility span and x_d, the age at which the mortality risk becomes non-null. Here, we generalized the model to any intensities of birth and death denoted (i_b, i_d) as well as to populations without Lansing effect (Figure 1, see also Appendix 1). The selective pressure is enforced by a logistic competition c mimicking a maximum carrying capacity of the environment, thus no explicit adaptive value is given to any particular trait. Additionally, for each reproduction event, a mutation (h) of probability p can affect both genes x_b and x_d independently, following a Gaussian distribution centered on the parental trait. In Figure 1, the different cases are explored, depending on the respective values of x_b and x_d. Individuals in the Figure 1b–c configuration (for x_b≤x_d) will always give progeny with a genotype (x_b, x_d) ∓ (h_b and/or h_d). The evolutionary outcome of individuals carrying a genotype with x_d <x_b (Figure 1a) is slightly more nuanced and depends on the parental age a and whether the parent carries the possibility for a Lansing effect or not (Figure 1d–f). If a<x_d, or if the parent does not carry a Lansing effect, the genotype of the progeny will be as previously described. But if a>x_d, and if the parent carries the Lansing effect, the progeny then inherits a dramatically reduced x_d (here x_d is set to 0), mimicking a strong Lansing effect.

Figure 1

Download asset Open asset

Three typical configurations of the model with i_b >i_d and their effect on progeny’s genotypes as a function of parental age.

(upper panel) Each haploid individual is defined by a parameter x_b defining its fertility span of intensity *i_b* and a parameter x_d defining the time during which it will maintain itself, with an intensity *i_d*. These parameters can be positive or null. (a) ’Too young to die’: it corresponds to configurations satisfying x_d <x_b. (b) ’Now useless’: it corresponds to configurations where x_b = x_d. (c) ‘Menopause’: it corresponds to configurations where x_d >x_b. (lower panel) Each individual may randomly produce a progeny during its fertility span [0; x_b]. (d) In the case of physiologically young parents (a<x_d), the progeny’s genotype is that of its parent ∓ a Gaussian kernel of mutation centered on the parental gene. In the case of the reproduction event occurring after x_d, for configuration (a) above, two cases are observed, (e) if the organism carries a Lansing effect ability, the x_d of its progeny will be strongly decreased. (f) In the absence of the Lansing effect, the default rule applies.

In our previous work (Méléard et al., 2019), we formally and numerically showed the long-time evolution of the model to converge towards (x_b - x_d)=0 in the case of individuals carrying a Lansing effect. Here, we explore the convergence of (x_b - x_d) without the strong transgenerational effect of ageing. We implemented a new version of the model, devoid of the Lansing effect, and simulated its evolution for a viable - that is allowing the production of at least one progeny - trait (x_b = 1.2, x_d = 1.6). Surprisingly, we still observe a convergence of (x_b - x_d) in finite time. The dynamics of the trait (x_b, x_d) is described by the canonical equation of adaptive dynamics, which depends on the Malthusian parameter and its gradient (Appendix 1). The Malthusian parameter can be interpreted as the age-specific strength of selection (Hamilton, 1966). The speed at which x_b and x_d evolves, decreases with time, just as in the previous form of the model (Méléard et al., 2019), allowing us to recover the well-observed, age-related decrease in the strength of selection (Haldane, 1941; Hamilton, 1966; Medawar, 1953). Simulations of the generalized bd model presented here show that the x_b - x_d distance (the time separating the end of fertility from the increasing risk of death) converges, for any initial trait, towards a positive constant. Thus, the long-term evolution of such a system is a configuration similar to Figure 2a (x_d <x_b). The formal analysis of the generalized bd model confirms that the long-time limit of the traits (x_b - x_d) is the positive constant (defined by the formula in Figure 2b, mathematical analysis presented in Appendix 1), reached after a few dozen simulated generations (Figure 2c). Although we formally demonstrate the long-time limit for any i_b and i_d, all our simulations are run using i_b=i_d=1, in order to limit the number of conditions to assess and report. Surprisingly, the limit value of the trait is not affected by x_b or x_d values - the fertility span and mortality per se - but only by their respective intensities i_b and i_d. These intensities can be interpreted as the instant mortality risk i_d and the probability to give a progeny i_b. Interestingly, the long-time limit values for any i_b and i_d shows a significantly stronger sensitivity to the increasing mortality risk i_d than to reproduction by almost two orders of magnitude (Figure 2d). In addition, for extremely low values of i_b and i_d - that is below 0.01 - the apparent time correlation of the fertility span and mortality onset is almost nonexistent; this is because (x_b - x_d) is large. Biologically, this would appear to an observer as the mortality onset occurring long before the exhaustion of reproductive capacity. Such an organism would be thus characterized as having no significant fertility decrease during the ageing process. On the other end, for individuals showing either a high instant mortality risk or a high probability to give a progeny, the (x_b - x_d) trait is close to 0, meaning that fertility and organismal integrity maintenance are visibly - that is observable by an experimenter - correlated. It is important to note that this mathematical study concerns individuals for which the mean number of descendants per individual is large enough, allowing us to define a viability set of traits (x_b, x_d) (see data availability statement). Because of these mathematical properties, a tradeoff emerges between i_b, i_d, x_b and x_d. Let’s consider an organism - for both the Lansing and non-Lansing cases - with a low reproductive intensity i_b=0.01 and i_d=1. For this organism to propagate, the product i_b * x_b has to be strictly superior to 1, hence here x_b ⩾ 100 (see data availability statement). In this example, the long-time limit of the trait (x_b - x_d) is equal to log(2), thus x_b and x_d are of the same order of magnitude. With the same reasoning, the long-time evolution lower limit of (x_b - x_d), of an organism that is significantly more fertile (with i_b=1, i_d=1), is 1/√3. This model thus allows an elegant explanation for the apparent negative correlation previously described between longevity and fertility without the need of implementing energy trade-offs or relative efficiency of energy allocation between maintenance and reproduction (see data availability statement - examples).

Figure 2

Download asset Open asset

The bd model shows a convergence of x_b - x_d towards a positive value.

Dynamics of the individual-based model shows a convergence of x_b - x_d towards a positive constant value in the absence of the Lansing effect. (a) The generalized b-d model shows a convergence of (x_b - x_d) for any i_b and i_d towards a positive value given by (b) (Annexe 4.3, Figure 2). (c) Simulation of 1000 individuals with initial trait (x_b=1.2, *x_d*=1.6) of intensities *i_b*=*i_d*=1, a competition c=0.0009 and a mutation kernel (P=0.1, σ=0.05) show that the two parameters co-evolvetowards x_b - x_d ≅ 0.55 that is log(3)/2. (d) Landscape of solutions (x_b - x_d) as a function of i_b and i_d (colors separate ranges of 50 units on the z-axis).

In our model, regardless of the initial trait (x_b, x_d) in the viability set, evolution leads to a configuration of the trait such that the risk of mortality starts to increase before the fertility span ends. Similar to biochemical reactions involved in a given pathway that are evolutionarily optimized (e.g. through tunneled reactions and gated electron transfers), we hypothesize here that such a configuration, caused by simple mathematical constraints, creates the conditions for the apparition, selection, and maintenance of a molecular mechanism coupling x_b and x_d. Such a coupling mechanism could thus be the so-called Lansing effect— the only described age-related decline in progeny’s fitness that seems to affect numerous iteroparous species (Lansing, 1947; Monaghan et al., 2020). We assessed the likelihood of survival of an organism carrying such a non-genetic and pro-senescence mechanism when in competition with a population devoid of such a mechanism. To do so, we examined a population divided into two sub-populations: one made of individuals subject to the Lansing effect and the other made up of individuals not subject to the effect. We assume, as before, that each individual is under the same competitive pressure. The two initial sub-populations have the same Darwinian fitness approximated by their Malthusian parameter (see data availability statement, Figure 3—figure supplement 1). Their traits are thus (1.5; 1.3)_Lansing and (1.5; 0.83)_non-Lansing. In order to simplify the analysis, both the birth and death intensities are as follows: i_b=i_d=1 (the model is nevertheless generalized to any (i_b; i_d), see data availability statement). We simulated the evolution of such mixed populations for discrete pairs of mutation rate (p) and competition (c) parameters. Three indexes were calculated for each set of simulation: (Table 1) the ratio of Lansing and non-Lansing populations that collapsed (“-” indicates that all survived), (Table 1) the ratio of total number of progenies produced during the simulation by each population and (Table 1) the relative proportion of the Lansing population at the end of the simulation (Table 1—source data 1). Our 1200 simulations, each with 2.10⁵ birth-death events, summarized in Table 1, show that the Lansing populations survive at least as well as non-Lansing ones (Table 1) especially for a moderate competition parameter (c=9.10^–4) and low (in our simulations) mutation rate (P=0.1). With such conditions, Lansing populations show almost half the risk of disappearance of non-Lansing ones (Table 1), producing nearly three times as many descendants as non-Lansing populations (Table 1), for up to a 20% faster growing population (Table 1). The plots for a single. Thus, although the Lansing effect gives way to a significant proportion of progeny with an extremely low fitness (x_d = 0), pro-ageing populations show a decrease in the risk of collapse. Moreover, we observe a slightly better growth of the population, independent of the magnitude of the Lansing effect (Figure 3—figure supplement 2).

Table 1

Populations with Lansing effect are favorably selected under logistic competition when the mutation rate is non-zero.

p is the mutation rate and c the intensity of the logistic competition. For each couple (p, c), 100 independent simulations were run with 500 individuals per population at t0 of which traits are (1.5; 1.3) Lansing and (1.5; 0.83)non-Lansing so that their respective Malthusian parameters are equal. Each simulation corresponds to 2.105 events of birth or death. Table (a) shows the ratio of Lansing and non-Lansing populations (out of 100 simulations in each case) that did collapse by the end of the simulation. For the lowest competition, none of the populations collapsed within the timeframe of simulations (-). For an intermediate value of competition, approximately less than half of Lansing populations disappear, relative to non-Lansing ones. Table (b) shows the ratio of the number of individuals generated between Lansing and non-Lansing populations. On average, Lansing populations generate approximately twice as many individuals as non-Lansing ones. (c) On average, Lansing populations grow 20% more than the non-Lansing. Values highlighted in green are discussed further below.

		Mutation probability
		0	0.1	0.5	1	0	0.1	0.5	1	0	0.1	0.5	1
		a) Lansing/non-Lansing collapsed population				b) Lansing/non-Lansing number of individuals				c) Lansing/total population size
Competition	9.10^–5	-	-	-	-	1.30	1.38	1.39	1.35	0.57	0.64	0.62	0.59
	9.10^–4	1.02	0.62	0.56	0.66	1.70	2.84	3.58	3.20	0.49	0.62	0.60	0.55
	9.10^–3	1.00	1.05	1.13	1.03	1.05	1.31	1.56	1.84	-	0.43	0.44	0.49

Table 1—source data 1 The magnitude of the Lansing effect does not influence the outcome of evolution. 100 independent simulations were run for each Lansing effect magnitude ranging from 0 (no Lansing effect) to 1 (progeny from parents age € [xd; xb] have xd = 0), starting with 500 Lansing (1.5; 1.3) and 500 non-Lansing (1.5; 0.83) individuals. We plot here the distribution density of xb - xd at the end of the simulation (individuals born in the time interval [990; 1000]), for Lansing populations (blue) and non-Lansing ones (red). Surprisingly, the magnitude of the Lansing effect does not seem to affect the optimal xb - xd solution value.: https://cdn.elifesciences.org/articles/92914/elife-92914-table1-data1-v1.zip
Download elife-92914-table1-data1-v1.zip

In order to understand the evolutionary success of a characteristic that seems to decrease an organism’s fitness, we computed the average Malthusian parameter of each population through time. We had previously identified that this intermediate set of c and p was associated with the highest success rate of Lansing bearing populations and presented the results for this set (highlighted in green, Table 1). First, we observe that, on average, Lansing populations (blue) grow while non-Lansing ones (red) decrease in size (Figure 3a - blue and red curves represent deciles 1, 5, and 9). In the simulations where both populations coexist, the higher fitness of the Lansing population is marginal, with these populations growing 20% more than the non-Lansing population (Figure 3b). This higher success rate seems to be driven by a faster and broader exploration of the Malthusian parameter space in the Lansing population (Figure 3c). This maximization of the Malthusian parameter is not associated with any significant difference of individual lifespan (time of death - time of birth) distributions of either population (Figure 3d). Although subjected to the same competition c, the distribution of the progeny from non-Lansing populations is essentially that of the parental trait in the first 5 generations, while Lansing progenies (not affected by the Lansing effect; we excluded progeny with x_d=0 for the comparison) explore a broader part of the trait space (Figure 3e). Interestingly, low fitness progeny (x_d=0) represents up to 10% of the population for a significant amount of time (Figure 3f). As a consequence, Lansing populations reach the equilibrium trait faster than the non-Lansing ones (Figure 3g). Thus, the relatively higher success rate of Lansing bearing populations seems to be associated with a higher genotypic diversity. This, in theory, leads to a broader range of fitness types. The ‘optimal’ fitness is therefore achieved earlier (or more easily), thus explaining the relative success of the population. This is an example of a population that demonstrates a greater ability to evolve (i.e. the population ‘possesses’ the attribute termed ‘evolvability’).

Figure 3 with 2 supplements see all

Download asset Open asset

The Lansing effect maximizes populational survival by increasing its evolvability.

100 independent simulations were run with a competition intensity of 9.10^–4 and a mutation rate p=0.1 on a mixed population made of 500 non-Lansing individuals and 500 individuals subjected to such effect. At t₀, the population size exceeds the maximum load of the medium thus leading to a population decline at start. At t₀, all individuals are of age 0. Here, we plotted a subset of the 100.10⁶ plus individuals generated during the simulations. Each individual is represented by a segment between its time of birth and its time of death. In each graph, blue and red curves represent deciles 1, 5, and 9 of the distribution at any time for each population type. (a) The higher success rate of Lansing bearing populations does not seem to be associated with a significantly faster population growth but with a lower risk of collapse. (b) For cohabitating populations, the Lansing bearing population (blue) is overgrowing by only 10% the non-Lansing one (red). (c) This higher success rate is associated with a faster and broader exploration of the Malthusian parameter - surrogate for fitness - space in Lansing bearing populations (d) that is not associated with significant changes in the lifespan distribution (e) but a faster increase in genotypic variability within the [0; 10] time interval. (f) This occurs although progeny from physiologically old parents can represent up to 10% of the Lansing bearing population and leads to it reaching the theoretical optimum within the timeframe of simulation (g) with the exception of Lansing progenies. (**e–g**) Horizontal lines represent the theoretical limits for (x_b - x_d) in Lansing (blue) and non-Lansing (red) populations.

Our model explains, in mathematical terms, why the mortality onset is evolutionarily linked to reproductive mechanisms (or fertility). Nevertheless, the numerical exploration of our model’s behavior has been limited so far to initial conditions, where the competing populations were of equal Malthusian parameters. The low number of generations involved suggests that the conditions for the development, selection, and maintenance of mechanisms of ageing (Lemoine, 2021) occurrs early on in evolutionary history, in a population of mixed individuals. As such, we decided to test the evolution of the trait (x_b - x_d) in Lansing and non-Lansing bearing individuals of uniformly distributed traits on [–10;+10] (Figure 4 - left panel). We chose to plot one (Figure 4 - central panel) of the hundred simulations we made. This simulation is representative of the general results. Simulations show, in over 110 million individuals, an early counter-selection of extreme trait values, typically (x_b - x_d)>4. Interestingly, the whole space of (x_b - x_d) trait is not explored evenly and the positive part of the trait space represents approximately 2/3 of the individuals (although the branched evolution process led to both the positive (‘Too young to die’ – Figure 1a) and negative (‘Menopause’ – Figure 1c) sides of the trait space). Both the Lansing and non-Lansing bearing populations manage to co-exist until the end of the simulation, each reaching a distribution centered on their respective theoretical solutions (Figure 4 - right panel): 0 for the Lansing (Méléard et al., 2019) and log(3)/2 for the non-Lansing. In this context, where the initial condition does not restrict the competition to individuals of identical Malthusian parameters, the Lansing bearing population is significantly less successful than the non-Lansing one (representing only one third of the final population size). As such, the evolution of a mixed population of individuals with a trait (x_b - x_d) initially uniformly distributed on [–10;+10], with or without a strong inter-generational effect, will lead to a mixed solution of individuals carrying a trait that converges towards the theoretical solution (such as x_d ≲ x_b), thus allowing the maximization of fertility without cluttering the environment with non-fertile individuals. This result is very similar to Weismann’s first intuition (Weismann, 1882). Nevertheless, this interpretation seems somehow finalist (i.e. presumes that the effects necessitate the causes) and does not yet discriminate why the Lansing population is evolutionarily successful in comparison to the non-Lansing population. Thus, we next explore the parameter of evolvability further, which leads us to yet again conceptualize ageing as an adaptive trait.

Figure 4

Download asset Open asset

Mixed populations lead to (x_b - x_d) theoretical limit in a limited time and cohabitation of Lansing and non-Lansing populations.

Starting with a homogenous population of 5000 Lansing bearing and 5000 non-Lansing individuals with traits uniformly distributed from –10 to +10 (left panel), we ran 100 independent simulations on time in [0; 1000]. (center panel) Plotting the trait (x_b - x_d) as a function of time for one simulation shows a rapid elimination of extreme traits and branching evolution. (right panel) The final distribution of traits in each population type is centered on the theoretical convergence limit for each. N_total ≅ 110 millionindividuals, c=9.10^–4, p=0.1.

Populations that consist of Individuals who can transmit ageing ‘information’ to the next generation are relatively more successful, within the framework of our model. Thus, to understand the origin of this pattern, we examined the differential landscape of the Malthusian parameters as a function of the trait (x_b, x_d) for both Lansing and non-Lansing populations. We built this landscape numerically using the Newton method (see data availability statement). First, it is interesting to notice that, from the equations, we have derived the maximum rate of increase for Malthusian parameters, this being 1/i_d with a maximum fitness value capped by i_b (data availability statement). Consistent with our previous characterization of the Trait Substitution Sequence in populations with Lansing effect (Méléard et al., 2019), Lansing individuals have a symmetrical fitness landscape (Figure 5, blue lines) centered on the diagonal x_b = x_d (Figure 5, green diagonal). Along the latter, we can directly observe what is responsible for a ‘selection shadow’. As x_b and x_d increase, a mutation of the same magnitude has smaller and smaller effects on the fitness, thus allowing the accumulation of mutations (Figure 5, blue arrows). The graphical representation of non-Lansing individuals is asymmetric— the rupture of symmetry occurs on the x_b = x_d diagonal. For x_d>x_b (Figure 5, upper diagonal), fitness isoclines of the two types of individuals fully overlap, thus showing an equal response of both Lansing and non-Lansing fitness to mutations. In addition, as expected, the fitness of Lansing individuals is equal to that of non-Lansing ones for a given trait. On the lower part of the graph, corresponding to x_d <x_b, non-Lansing fitness isoclines separate from that of Lansing individuals, making the fitness of non-Lansing individuals higher to that of Lansing ones for a given trait. Nevertheless, the fitness gradient is significantly stronger for Lansing individuals as represented within Figure 5 by the yellow arrow and associated yellow area. For an individual of trait (x_b = 2.45; x_d = 1.05), a mutation making a non-Lansing individual 0.1 in fitness (isocline 0.7 to isocline 0.8) will make a Lansing individual increase its own by 0.42 (isocline 0.1 to above isocline 0.5). With a fourfold difference, the Lansing population produces four times as many individuals as the non-Lansing ones for a given mutation probability. But this reasoning can be extended to any trait (x_b, x_d) with or without Lansing effect. Organisms ageing rapidly - that is with low x_b and x_d - will see their fitness significantly more affected by a given mutation h than individuals with slower ageing affected by the same mutation. As such, because ageing favors the emergence of genetic variants, ageing populations are therefore more evolvable.

Figure 5

Download asset Open asset

The Lansing effect is associated with an increased fitness gradient.

We were able to derive Lansing and non-Lansing Malthusian parameters from the model’s equations (see Annexe 1–2.3 and 1–5) and plot them as a function of the trait (x_b, x_d). The diagonal x_b = x_d is drawn in light green. The corresponding isoclines are overlapping above the diagonal but significantly differ below, with non-Lansing fitness (red lines) being higher than that of Lansing’s (light blue lines). In addition, the distance between two consecutive isoclines is significantly more important in the lower part of the graph for non-Lansing than Lansing bearing populations. As such, a mutation leading a non-Lansing individual’s fitness going from 0.7 to 0.8 (yellow arrow) corresponds to a Lansing individual’s fitness going from 0.1 to 0.52. Finally, Hamilton’s decreasing force of selection with age can be observed along the diagonal with a growing distance between two consecutive fitness isoclines as x_b and x_d continue increasing.

Discussion

Ageing is, despite its phenotypic diversity (Jones et al., 2014), an evolutionarily conserved phenomenon. How ageing evolved, however, is presently debated. Although early theories (Weismann, 1882) conceive ageing as adaptive, ageing is now generally viewed as a side-effect, or byproduct, of diminished selective pressure and therefore not adaptive.

The mathematical model we have presented here allows us to propose an alternative theory: ageing necessarily emerges for any system showing the two minimal properties of life (Trifonov, 2011), namely (a) reproduction with variation (x_b) and (b) organismal maintenance (x_d). We formally show that a haploid and asexual organism with these two properties will rapidly evolve, within a few dozen generations, towards a solution such that (x_b - x_d) is strictly positive, meaning that the risk of mortality starts to increase before the end of the fertility span. Importantly, the time separating both parameters is independent from their absolute values and only depends on the rate of each, respectively i_b for x_b and i_d for x_d. This explains the observed trade-offs (Kirkwood, 2005; Lemaître et al., 2015; Rodrigues and Flatt, 2016) between the fertility of an organism and its lifespan. Thus, our work addresses outstanding questions outlined in the disposable-soma theory (Kirkwood and Holliday, 1979) — why and how a highly fertile organism either dies or ages earlier. Indeed, the lower limit condition for the production of descendants by an individual in our model is x_b * i_b >1. As such, an organism with low fertility (i_b <<1) will obtain a progeny only if fertile longer (x_b >>1). Conversely, a highly fertile organism will evolve towards its minimum viable condition, requiring only a small x_d. The apparent trade-off between fertility and longevity is thus solely a consequence of x_b * i_b >1 and lim_+∞(x_b - x_d)_t. Our model need not implement any constraint on resource allocations or other tradeoffs for this effect to occur.

Because x_b and x_d converge, this favors the onset of a period in which an individual’s fertility drops while its risk of dying becomes non-zero; this is the organism entering the ‘senescence phase’ corresponding to the Smurf phase described in Rera et al., 2012. This necessary convergence of fertility’s end and senescence’s start would thus facilitate the selection of any molecular mechanism that couples the two processes (Echave, 2021). Additionally, and in opposition to what is suggested in Stearns, 1989, we observe that any two genes that are not functionally linked can be co-selected.

While the Lansing effect somewhat decreases the fitness of individuals within a population, the probability of survival of a population is significantly greater in Lansing populations when in competition with a non-Lansing population of equal Malthusian parameter at t₀. We observed, numerically, that this slight increase in survival is mediated by an increase in the genetic variability generated within the population. Thus, we propose that such an active mechanism of ageing can be selected during evolution through its ability to increase an organism’s evolvability. As mentioned above, evolvability is understood as the ‘the capacity to generate heritable selectable phenotypic variation’ (Kirschner and Gerhart, 1998). It is an interesting concept as it allows for a trait that has no direct effect on fitness - even a negative one (Maynard Smith, 1971) - to be under strong selection, given its ability to generate genetic or phenotypic variation. Furthermore, such a two-phase mechanism would be of great advantage in a constantly varying environment. Indeed, when environmental conditions become less permissive, x_d might be affected and individuals would be pushed to enter the [x_d; x_b] space earlier, thus increasing the evolvability of the population. This is what we observe in the laboratory where individuals submitted to harsh conditions will enter the Smurf phase earlier than the control conditions (Rera et al., 2012). Regarding the nature of the transgenerational effect, our model is agnostic and the mere transmission of any negative effect would be sufficient to exert the function.

Because we, without fail, observe the convergence of the end of fertility and the start of senescence, our generalized model - supported by a formal analysis - predicts a high degree of conservation of ageing, specifically as something that can be selected. This gives rise to organisms that lose homeostatic capacities amidst and during the period of fertility. We have identified a mathematical constraint that explains the biphasic pattern of ageing proposed in Tricoire and Rera, 2015, allowing for the positive selection of ageing through evolutionary time. More importantly, the negative impacts of ageing on individuals’ fitness seem to be fully compensated at the population level. Our work, at large, thus demonstrates the following: (1) fertility and senescence always converge if an organism is both fertile and homeostatic, (2) ageing populations are more successful through time, and (3) more evolvable. Therefore, we defend that ageing can, in theory, be re-conceptualized as adaptive.

This two phase model is very simple, yet able to describe all types of ageing observed in the wild, including a rapid post-reproductive onset of mortality, a menopause-like mortality plateau, and what we have identified as a two-phase Smurf-like process. The strong mathematical constraint between x_b, x_d, i_b, and i_d limits the possible configurations. Additionally, our mathematical model of ageing, as a two-phase process (Tricoire and Rera, 2015), shows that the mortality rate of the second phase of life is considerably constant across Drosophila lines of significantly different life expectancies, ranging from 15 to 70 days. In these conditions, if i_d is a constant parameter, can we experimentally affect x_d by acting on i_b and/or x_d? Experimental evolution using only Drosophila progeny conceived later in the life of the parent has shown that the onset of mortality, within these progeny, occurs rather late, sometimes even after the end of the fertility period (Burke et al., 2016; Rose et al., 2002). Although the authors report previous studies of their own with divergent results, other independent experiments have led to results suggesting an increase of x_d following an artificial increase of x_b (Luckinbill and Clare, 1985; Sgrò et al., 2000) as well as the reverse (Stearns et al., 2000).

Without the need to implement resource allocation constraints, pleiotropic antagonistic functions nor late-life accumulation of mutations, our model is able to predict the evolution of ageing while encompassing phenomena that previously led to the two above-mentioned theories (mutation accumulation and antagonistic pleiotropy). More importantly, our model suggests a central role of ageing in evolution, as the mathematical constraint we show is likely to apply to any function affecting fertility and homeostasis. Could this broader application of constraints be responsible for the stereotyped gene expression changes - reminiscent of the so-called hallmarks of ageing - we recently described in Smurfs (Zane et al., 2023)? Although this model helps us to see the conditions under which ageing is an evolutionarily adaptive force, it is still a toy model. The mortality and fertility functions we used are binary and we are now developing more complex versions of the model, notably to assess the interactions existing between i_b, i_d, x_b, and x_d but more importantly to assess their co-evolution with maturation, sex, ploidy, or varying environmental conditions.

Materials and methods

See data availability statement for code, packages and the software used.

Appendix 1

1 The mathematical individual-based bd model

We model an haploid and asexual population of individuals with evolving life-histories by a stochastic individual-based model, similar to the one introduced in Méléard et al., 2019 and a particular case of Ferrière and Tran, 2009. Each individual is characterized by its age and by a life-history trait $x = (x_{b}, x_{d}) \in R_{+}^{2}$ that describes for each individual the age x_b at the end of reproduction and the age x_d when mortality becomes positive. The trait can change through time, by mutations occuring continuously in time.

More precisely, the Markovian dynamics of the population process is defined as follows. The individuals reproduce and die independently. An individual with trait $(x_{b}, x_{d})$ reproduces at rate i_b as long as it is younger than x_b. Further, he cannot die as long as it is younger than x_d and has a natural death rate i_d after age x_d.

The life-history of an individual with trait $x = (x_{b}, x_{d})$ is described by the couple of birth and death functions $(B_{x}, D_{x}^{c})$ defined on $R_{+}$ by

\forall a \in R_{+}, B_{x} (a) = i_{b} 1_{a \leq x_{b}}, D_{x}^{c} (a) = i_{d} 1_{a > x_{d}} + c N .

Here, the individual age $a$ is the physical age, $N$ the (varying) population size and $c > 0$ the competition pressure exerted by an individual on another one. The death rate will be extended to

D_{x}^{c} (0) = + \infty, for x_{b} < 0 or x_{d} > 0,

meaning that an individual appearing by mutation will be able to survive only if the two components if its trait are non negative.

Note that the date of birth and lifespan of an individual are stochastic and the law of the lifespan on an individual with trait $x$ born at time $τ$ is given by $f_{x} (s) = D_{x}^{c} (s) \exp (- \int_{τ}^{τ + s} D_{x}^{c} (a) d a)$ .

We also take into account genetic mutations which create phenotypic variation, and which added to competition between individuals, will lead to natural selection.

At each reproduction event, a mutation appears instantaneously on each trait x_b and x_d independently with probability $p \in] 0, 1 [$ . If the trait $x_{b}$ mutates (resp. if $x_{d}$ mutates), the trait of the newborn is $x_{b} + h_{b}$ (resp. $x_{d} + h_{d}$ ). The mutation effect $h_{b}$ (resp. $h_{d}$ ) is distributed following a centered Gaussian law with variance $σ^{2}$ . This Gaussian law is denoted by $k (h) d h$ .

Note that a similar model has been defined in Méléard et al., 2019, including a Lansing effect on the reproductive lineage of "old" individuals.

2 The Malthusian parameter

2.1 The demographic parameters

We now introduce the classical demographic parameters for age-structured (without competition) population, where all individuals have the same trait $x \in R_{+}^{2}$ (Charlesworth, 1994). We are looking for a triplet $(λ (x), N_{x}, ϕ_{x})$ where $λ (x) \in R$ is the Malthusian parameter, $N_{x} (a), a \in R_{+}$ the stable age distribution and $ϕ_{x} (a), a \in R_{+}$ the reproductive value. They describe the asymptotic growth of the population dynamics and measure the fitness of life-histories: $λ (x)$ is the growth rate of the population at its demographic equilibrium, $N_{x}$ the age distribution of the population and $ϕ_{x} (a)$ is the probability that an individual with trait $x$ has a newborn after age $a$ . It is known (Charlesworth, 1994), that $(λ (x), N_{x}, ϕ_{x})$ is solution of the direct and dual eigenvalue problems:

{\begin{cases} - \partial_{a} N_{x} (a) - D_{x} (a) N_{x} (a) = λ (x) N_{x} (a) \\ N_{x} (0) = \int_{0}^{+ \infty} B_{x} (α) N_{x} (α) d α, N_{x} (0) = 1, \end{cases}

{\begin{cases} \partial_{a} ϕ_{x} (a) - D_{x} (a) ϕ_{x} (a) + B_{x} (a) ϕ_{x} (0) = λ (x) ϕ_{x} (a) \\ ϕ_{x} (0) = 1 \end{cases}

where $B_{x} (a) = i_{b} 1_{{a > x_{b}}} a n d D_{x} (a) = i_{d} 1_{{a > x_{d}}}$ .

Proposition 2.1

For all $x \in R_{+}^{2}$ , there exists a unique solution $(λ (x), N_{x}, ϕ_{x}) \in R \times L^{1} (R_{+}) \times L^{\infty} (R_{+})$ of (2) and (3). The Malthusian parameter $λ (x)$ is the unique solution of the equation:

i_{b} \int_{0}^{x_{b}} e^{- i_{d} (a - x_{d})_{+} - λ (x) a} d a = 1.

The stable age distribution $N_{x}$ and the reproductive value $ϕ_{x}$ verify

N_{x} (a) = e^{- i_{d} (a - x_{d})_{+} - λ (x) a}, ϕ_{x} (a) = \frac{i_{b} 1_{a \leq x_{b}}}{N_{x} (a)} \int_{a}^{x_{b}} N_{x} (α) d α .

Proof

The proof is straightforward by solving the first equations in (2) and (3), and then by using the equations satisfied by the boundary conditions. □

Remark 2.2

The quantities $λ (x), N_{x}, ϕ_{x}$ are the eigenelements (Proposition 2.1) associated with the linear operator that generates the dynamics $v_{x} (t, a)$ of a non density dependent population with age structure and birth-death rates given by $(B_{x}, D_{x})$ . More precisely, $v_{x} (t, a)$ satisfies the McKendrick Von-Foerster Equation

{\begin{cases} \partial_{t} v_{x} (t, a) + \partial_{a} v_{x} (t, a) = - D_{x} (a) v_{x} (t, a), t \geq 0, a \geq 0 \\ v_{x} (t, 0) = \int_{R_{+}} B_{x} (α) v_{x} (t, α) . \end{cases}

The use of these quantities as an indicator of fitness is justified by the convergence of $e^{- λ (x) t} v_{x} (t, a)$ to $(\int_{R_{+}} v_{x} (0, α) ϕ_{x} (α) d α) N_{x} (a)$ as $t$ tends to infinity (Perthame, 2007 for example).

2.2 Computation and regularity of the Malthusian parameter

The Malthusian parameter $λ (x)$ is defined as the unique real number such that

i_{b} \int_{0}^{x_{b}} e^{- i_{d} (a - x_{d})_{+} - λ (x) a} d a = 1.

Let us introduce

U_{1} = {x \in R_{+}^{2} : x_{b} < x_{d}}, U_{2} = {x \in R_{+}^{2} : x_{d} < x_{b}}, H = {x \in R_{+}^{2} : x_{b} = x_{d}} .

For all $x \in U_{1} \cup H$ , the Malthusian parameter $λ (x)$ satisfies:

i_{b} \int_{0}^{x_{b}} e^{- λ (x) a} d a = 1 = \frac{i_{b}}{λ (x)} (1 - e^{- x_{b} λ (x)}) .

Then $λ (x)$ can be numerically computed by Newton’s method applied to the function $K_{x_{b}} (λ) = \frac{1}{λ} (1 - e^{- x_{b} λ}) - \frac{1}{i_{b}}$ , since $λ (x)$ is solution of $K_{x_{b}} (λ) = 0$ , .

In the case where $x \in U_{2}$ , we have

i_{b} \int_{0}^{x_{b}} e^{- i_{d} (a - x_{d})_{+} - λ (x) a} d a = i_{b} \int_{0}^{x_{d}} e^{- λ (x) a} d a + i_{b} \int_{x_{d}}^{x_{b}} e^{- i_{d} (a - x_{d}) - λ (x) a} d a

= i_{b} {\frac{1}{λ (x)} (1 - e^{- λ (x) x_{d}}) + \frac{e^{i_{d} x_{d}}}{λ (x) + i_{d}} (e^{- (λ (x) + i_{d}) x_{d}} - e^{- (λ (x) + i_{d}) x_{b}})},

which has to be equal to 1. That involves a function

H_{(x_{b}, x_{d})} (λ) = \frac{1}{λ} (1 - e^{- λ_{d}}) + \frac{e^{i_{d} x_{d}}}{λ + i_{d}} (e^{- (λ + i_{d}) x_{d}} - e^{- (λ + i_{d}) x_{b}}) - \frac{1}{i_{b}} .

Newton’s method still allows to resolve numerically the equation and find $λ (x)$ .

Let us now prove some regularity properties of the Malthusian parameter. We show that its gradient is a simple function of the stable age distribution, the reproductive value and the mean generation time $G$ defined for all $x$ by

G (x) = i_{b} \int_{0}^{x_{b}} a N_{x} (a) d a .

Proposition 2.3

The function $x \in (R_{+}^{*})^{2} \mapsto λ (x)$ is of class $C^{1}$ and we have:

\forall x \in (R_{+}^{*})^{2}, \nabla λ (x) = \frac{1}{G (x)} (i_{b} N_{x} (x_{b}), i_{d} N_{x} (x_{d}) ϕ_{x} (x_{d})) .

Note that the derivatives are positive, meaning that $x_{b} \to λ (x)$ and $x_{d} \to λ (x)$ are non decreasing.

Proof

Coming back to the definition of λ and using the the implicit function theorem, we obtain that λ is differentiable and

\begin{array}{ll} x \in U_{1}, \frac{\partial λ (x)}{\partial x_{b}} & = \frac{i_{b} e^{- λ (x) x_{b}}}{G (x)} = \frac{i_{b} N_{x} (x_{b})}{G (x)}; \frac{\partial λ (x)}{\partial x_{d}} = 0 = \frac{N_{x} (x_{d}) ϕ_{x} (x_{d})}{G (x)}; \\ \forall x \in U_{2}, \frac{\partial λ (x)}{\partial x_{b}} & = \frac{i_{b} e^{- i_{d} (x_{b} - x_{d})} e^{- λ (x) x_{b}}}{G (x)} = \frac{i_{b} N_{x} (x_{b})}{G (x)} \\ \frac{\partial λ (x)}{\partial x_{d}} & = \frac{i_{b} i_{d} e^{i_{d} x_{d}} \int_{x_{d}}^{x_{b}} e^{- (i_{d} + λ (x)) a} d a}{G (x)} = \frac{i_{b} i_{d} N_{x} (x_{d}) ϕ_{x} (x_{d})}{G (x)} . \end{array}

We deduce that $λ$ has continuous partial derivatives, which concludes the proof. □

2.3 Viability set

The viability set is the set $V \subset R_{+}^{2}$ of traits $x = (x_{b}, x_{d})$ such that $λ (x) > 0$ . From Equation 4, $λ (x) > 0$ if and only if the mean number $R (x_{b}, x_{d})$ of descendants per individual is larger than one, i.e if and only if we have:

R (x_{b}, x_{d}) := i_{b} \int_{0}^{x_{b}} e^{- i_{d} (a - x_{d})_{+}} d a > 1.

A precise characterization of the set $V$ is given in Lemma 2.4. In Figure 1, we represent the set $V$ for $i_{b} = 1.5$ and $i_{d} = 2$ .

Appendix 1—figure 1

Download asset Open asset

The set $V = {(x_{b}, x_{d}) \in R_{+}^{2}; R (x_{b}, x_{d}) > 1}$ is the convex set delimited by the black curve with equation $R (x_{b}, x_{d}) = 1$ .

Lemma 2.4

We have:

V = {x \in R_{+}^{2}; x_{b} > x_{d} - \log (i_{d} x_{d} + 1 - (i_{b} / i_{d})) if x_{b} > x_{d}; i_{b} x_{b} > 1 if x_{b} \leq x_{d}},

and for all $x \in V$ , $λ (x) \leq i_{b}$ . Moreover, the map $x \in V \mapsto \nabla λ (x)$ is Lipschitz continuous.

Proof

We are looking for which $x = (x_{b}, x_{d}) \in R_{+}^{2}$ , the mean number of descendants $R (x_{b}, x_{d})$ is greater than 1. Recall that $R (x) = i_{b} \int_{0}^{x_{b}} \exp (- i_{d} \int_{0}^{a} 1_{α > x_{d}} d α) d a$ . For $x \in U_{1}$ (defined in Equation 6), we have $R (x) = i_{b} x_{b}$ and $R (x) > 1$ if and only if $i_{b} x_{b} > 1$ . For $x \in U_{2}$ , we have $R (x) = i_{b} x_{d} + (i_{b} / i_{d}) - (i_{b} / i_{d}) e^{- i_{d} (x_{b} - x_{d})}$ and $R (x) > 1$ if and only if $x_{b} > x_{d} - \log (i_{d} x_{d} + 1 - (i_{b} / i_{d}))$ . We conclude for the first assertion arguing that the map $λ \mapsto i_{b} \int_{0}^{x_{b}} \exp (- i_{d} \int_{0}^{a} 1_{α > x_{d}} i_{d} α - λ a) d a$ is decreasing. Let us now show that $λ (x)$ is upper-bounded by i_b. Assume that there exists $x \in V$ such that $λ (x) > i_{b}$ . Then

1 = i_{b} \int_{0}^{x_{b}} e^{- i_{d} (a - x_{d})_{+} - λ (x) a} d a < i_{b} \int_{0}^{x_{b}} e^{- i_{b} a} d a = 1 - e^{- x_{b}},

which is absurd and allows us to conclude. The next claim is shown arguing that the map $x \in V \mapsto \nabla λ (x)$ is differentiable on $U_{1} \cup U_{2}$ and admits bounded partial derivatives. □

Let us develop different examples:

In the case where $i_{b} = 0.01$ and $i_{d} = 1$ , we obtain

R (x_{b}, x_{d}) > 1 ⟺ x_{d} + (1 - (2.01)^{- \frac{1}{1.01}}) > 100,

which gives essentially that x_d has to be greater than 100.

In the case where $i_{b} = i_{d}$ , the formula is simpler. We obtain

R (x_{b}, x_{d}) > 1 ⟺ x_{d} + \frac{1}{i_{b}} (1 - 3^{- \frac{1}{2}}) > \frac{1}{i_{b}} .

We deduce

R (x_{b}, x_{d}) > 1 ⟺ x_{d} > \frac{1}{i_{b} \sqrt{3}} .

If we assume that $i_{b} = i_{d} = 1$ then we obtain that

R (x_{b}, x_{d}) > 1 ⟺ x_{d} > \frac{1}{\sqrt{3}} = 0.577.

Let us finally note that if we assume to be in the limit of the canonical equation and then to be in the case when $x_{b} - x_{d} = \frac{\log 3}{2 i_{b}}$ , we also obtain a characterization of the viability set using x_b:

R (x_{b}, x_{d}) > 1 ⟺ x_{b} > \frac{1}{i_{b}} (\frac{1}{\sqrt{3}} + \frac{\log 3}{2}) .

For $i_{b} = 1$ , that gives $x_{b} > 1.126$

3 Monomorphic equilibrium

Let us come back to the general case with competition, but for a monomorphic population with trait $x$ (and then without mutation). It can be proved (Méléard and Tran, 2009 Proposition 2.4) that for a large population, the stochastic process converges in probability to the solution of the following Gurtin-MacCamy partial differential equation (see Gurtin and Maccamy, 1974).

{\begin{cases} \partial_{t} n_{x} (t, a) + \partial_{a} n_{x} (t, a) = - (D_{x} (a) + c \int_{R_{+}} n_{x} (t, α) d α) n_{x} (t, a) \\ n_{x} (t, 0) = \int_{0}^{+ \infty} B_{x} (α) n_{x} (t, α) d α, (t, a) \in R_{+}^{2} . \end{cases}

This equation describes the density-dependent dynamics of a large population with trait $x$ (without mutation). The trait $x \in R_{+}^{2}$ being given, let us study the positive equilibria of the equation

For $x \in V$ , Equation 9 admits a unique non-trivial solution:

Proposition 3.1

For all $x \in V$ , there exists a unique globally stable equilibrium ${\bar{n}}_{x}$ to Equation 9, that is a solution of

{\begin{cases} - \partial_{a} {\bar{n}}_{x} (a) - (D_{x} (a) + c \int_{R_{+}} {\bar{n}}_{x} (α) d α) {\bar{n}}_{x} (a) = 0 \\ {\bar{n}}_{x} (0) = \int_{0}^{+ \infty} B_{x} (α) {\bar{n}}_{x} (α) d α, \end{cases}

which satisfies $λ (x) = c \int_{R_{+}} {\bar{n}}_{x} (α) d α$ .

Note that

{\bar{n}}_{x} (0) = \frac{λ (x)}{c \int_{0}^{+ \infty} N_{x} (α) d α} .

Proof

The existence part of the proof is trivial from Cansell et al., 2023 and Proposition 2.1 using that $V = {x \in R_{+}^{2} : λ (x) > 0}$ . The long-time behavior of the solutions of (9) is studied in [Webb, 1985, Section 5.4]. □

4 Canonical equation

4.1 Invasion fitness

We now compute the invasion fitness function associated with the individual-based model. We use the definition of invasion fitness given in Méléard and Tran, 2009. The invasion fitness $1 - z (y, x)$ of a mutant with trait $y$ in a resident population with trait $x$ is defined as the survival probability of an age-structured branching process with birth rates $B_{x} (a)$ and death rates $D_{x} (a) + c \int_{R_{+}} {\bar{n}}_{x} (a) d a$ .

Proposition 4.1

Let $y \in (R_{+}^{*})^{2}$ and $x \in V$ , we have

1 - z (y, x) = [\frac{λ (y) - λ (x)}{i_{b}}]_{+} .

Proof

The proof is a direct application of Equation (3.6) in Méléard and Tran, 2009. □

4.2 Trait Substitution sequence and Canonical equation

For this part, we refer principally to Méléard and Tran, 2009 where the Theory of Adaptive Dynamics is rigorously developed for general age-structured populations.

We introduce the canonical equation describing the evolution of the trait $x = (x_{b}, x_{d})$ at a mutation time-scale, under the assumptions of adaptive dynamics (large population, rare and small mutation, invasion and fixation principle, as well known since Metz et al., Dieckman-Law). In Méléard and Tran, 2009, it is shown that this equation can be obtained as a two-step limit from the individual based model. The first step consists in defining the Trait Substitution Process describing the successive advantageous mutant invasions in monomorphic populations at equilibrium. It is obtained as support dynamics of the measure-valued limit of the rescaled population process (at the mutation time-scale), when mutations are rare (but not small). The measure-valued limiting process is rigorously derived from the individual-based model in Dambroise et al., 2016 Section 3. It jumps from a state $δ_{x} (d z) n_{x} (a) d a$ to a state $δ_{y} (d z) - n_{y} (a) d a$ . The trait support process takes values in $V$ and its dynamics is described as follows.

Definition 4.2

The Trait Substitution Sequence is the càdlàg process $(X_{t}, t \geq 0)$ with values in $V$ whose law is characterized by the infinitesimal generator $L$ defined for all bounded and measurable function $φ : V \to R$ by:

L φ (x) = \int_{R^{2}} (φ (x + (h_{1}, h_{2})) - φ (x)) [\frac{λ (x + (h_{1}, h_{2})) - λ (x)}{i_{b}}]_{+} \frac{λ (x)}{c \int_{0}^{+ \infty} N_{x} (a) d a} μ (d h_{1}, d h_{2}),

where $μ (d h) = \frac{δ_{0} (d h_{2}) k (h_{1}) d h_{1} + δ_{0} (d h_{1}) k (h_{2}) d h_{2}}{2}$ and the distribution $k$ has been defined in Section 1.

Note that since by Proposition 2.3, the partial derivatives of $λ$ are positive, then the increment $λ (x + (h_{1}, h_{2})) - λ (x)$ is non negative if and only if h₁ and h₂ are non negative.

The second step consists in assuming that mutation amplitudes are small and of order $ϵ$ , for $ϵ > 0$ . We then define the rescaled process $X^{ϵ}$ by $X^{ϵ} (t) = ϵ X (\frac{t}{ϵ^{2}})$ . introduce the Canonical Equation that describes the limit behaviour of the Trait Substitution Sequence when mutations are small.

Proposition 4.3

Let $T > 0$ . Assume that $X^{ϵ} (0)$ converges to $x^{0} \in V$ in probability. Then the sequence of processes $(X^{ϵ})_{ϵ}$ converges in law in the Skorohod space $D ([0, T], V)$ to the solution $(x (t), t \geq 0)$ of the ordinary differential equation:

\frac{d x}{d t} = \frac{λ (x)}{4 c \int_{0}^{+ \infty} N_{x} (α) d α} \frac{\nabla λ (x) σ^{2}}{i_{b}}, x \in V \subset R_{+}^{2}

Recall that the Malthusian parameter $λ (x)$ is defined in Fabian, 2011, the stable age distribution $N_{x}$ is defined in Frenk and Houseley, 2018 and $σ^{2} (x)$ denotes the variance of the mutation kernel. Recall that (see Proposition 2.3)

\nabla λ (x) = \frac{1}{G (x)} (i_{b} N_{x} (x_{b}), i_{d} N_{x} (x_{d}) ϕ_{x} (x_{d})) .

It describes the strength of selection at ages x_b and x_d. Hence, this canonical equation allows to interpret the age specific strength of selection at ages x_b and x_d as the evolution speed of the traits x_b and x_d respectively, under the assumptions of adaptive dynamics.

Proof

The proof is classical and can be easily adapted from that of [6, Theorem 4.1]. The canonical equation only charges the set $V \subset R_{+}^{2}$ (defined in Section 2.1) and writes as follows:

\frac{d x}{d t} = - \nabla_{y} z (x, x) \frac{{\bar{n}}_{x} (0)}{2} \int_{R_{+}} h^{2} k (h) d h, x \in V .

The set $V$ is the set of traits that admit a positive stable monomorphic equilibrium ${\bar{n}}_{x}$ in a such way that ${\bar{n}}_{x} (0)$ equals the birth rate of a mutant (see Proposition 3.1); $σ^{2}$ is the variance of the mutations and $1 - z (y, x)$ is the invasion fitness. Computing these parameters gives (Equation 11). □

In Figure 2, we present a simulation of a solution of Equation 11. We observe that the traits x_b and x_d increase with time (Figure 2a, b), with decreasing speed tending to zero. The trait $x_{b} (t) - x_{d} (t)$ converges to some positive number (Figure 2c) that we can rigorously compute. That is the aim of the next section.

Appendix 1—figure 2

Download asset Open asset

Simulation of the canonical equation with $x^{0} = (3.5, 1.3)$ and $i_{b} = i_{d} = 1$ .

(a) Dynamics of x_b. (b) Dynamics of x_d. (c) Dynamics of, $x_{b} - x_{d}$ the black curve has equation. $y = \log (3) / 2$ .

4.3 Long-time behaviour of the canonical equation

In this section we study the long-time behaviour of the solutions of the Canonical Equation 11. We prove the following theorem.

Theorem 4.4

Let $x^{0} \in V$ and let $(x (t), t \geq 0)$ be the solution of (11) started at $x^{0} \in V$ . Then we have:

x_{b} (t) - x_{d} (t) \underset{t \to + \infty}{⟶} \frac{\log (1 + \frac{i_{b} + i_{d}}{i_{d}})}{i_{b} + i_{d}} .

We first prove the following lemma. We always denote $U_{1} = {x \in V : x_{b} < x_{d}}$ , $U_{2} = {x \in V : x_{d} < x_{b}}$ and $H = {x \in V : x_{d} = x_{b}}$ .

Lemma 4.5

Let $x^{0} \in V$ and let $(x (t), t \geq 0)$ be the solution of Equation 11 started at $x^{0} \in V$ . Then we have:

(i) There exists $T > 0$ , such that for all $t \geq T$ , $x (t) \in U_{2} \cup H$ .
(ii) There exists $C > 0$ such that for all $t \geq 0$ , $| x_{b} (t) - x_{d} (t) | < C$ ,
(iii) We have $x_{b} (t)$ increases to $+ \infty$ , $x_{d} (t)$ increases to $+ \infty$ and $λ (x (t)) \to i_{b}$ as $t \to + \infty$ .

Proof

For all $x \in V$ , let us define:

v (x) = \frac{λ (x) \int_{0}^{+ \infty} h^{2} k (h) d h}{2 i_{b} c G (x) \int_{0}^{+ \infty} N_{x} (α) d α} .

and we remark that there exist $\underline{v} (x^{0}), \bar{v} (x^{0}) > 0$ such that $\underline{v} (x^{0}) \leq v (x) \leq \bar{v} (x^{0})$ . (i): Let $T := inf {t \geq 0 : x (t) \in U_{2} \cup H} \in [0, + \infty]$ . We first show that $T < + \infty$ . If $x^{0} \in U_{2} \cup H$ , it is obvious. If $x^{0} \in U_{1}$ , assume that $T = + \infty$ . Then for all $t \geq 0$ , $x_{d} (t) = x_{d}^{0}$ . Indeed, as soon as $x_{b} < x_{d}$ , $ϕ (x_{d}) = 0$ and the trait x_d does not move (see Equation 7). We obtain that

\forall t \geq 0, \frac{d (x_{b} (t) - x_{d} (t))}{d t} = \frac{d x_{b} (t)}{d t} \geq v_{1} (x) i_{b} e^{- λ (x) x_{b}} \geq \underline{v} (x^{0}) i_{b} e^{- λ (x) x_{b}} > 0,

that allows us to obtain the contradiction. So we have $T < + \infty$ . We conclude the proof arguing that for all $t \geq 0$ such that $x (t) \in H$ , $d x_{d} (t) / d t = 0$ and $d x_{b} (t) / d t > 0$ . (ii): By (i), we assume without loss of generality that $(x (t), t \geq 0) \subset U_{2} \cup H$ . By Equation 11, we obtain that:

\begin{aligned} \frac{d (x_{b} (t) - x_{d} (t))}{d t} & = i_{b} v (x (t)) (e^{- i_{d} (x_{b} (t) - x_{d} (t)) - λ (x) (t) x_{b} (t)} - i_{d} \int_{x_{d} (t)}^{x_{b} (t)} e^{- i_{d} (a - x_{d} (t)) - λ (x (t)) a} d a) \\ = i_{b} v (x (t)) (e^{- i_{d} (x_{b} (t) - x_{d} (t)) - λ (x) (t) x_{b} (t)} \\ - \frac{i_{d} e^{i_{d} x_{d} (t)}}{i_{d} + λ (x (t))} (e^{- (i_{d} + λ (x (t))) x_{d} (t)} - e^{- (i_{d} + λ (x (t))) x_{b} (t)})) \\ = \frac{i_{b} v (x (t))}{i_{d} + λ (x (t))} ((2 i_{d} + λ (x (t))) e^{- i_{d} (x_{b} (t) - x_{d} (t)) - λ (x (t)) x_{b} (t)} - i_{d} e^{- λ (x (t)) x_{d} (t)}) \\ (7) & = \frac{i_{b} i_{d} v (x (t)) e^{- λ (x (t)) x_{d} (t)}}{i_{d} + λ (x (t))} (\frac{2 i_{d} + λ (x (t))}{i_{d}} e^{- (i_{d} + λ (x (t))) (x_{b} (t) - x_{d} (t))} - 1) . \end{aligned}

By Equation 14 and using the fact that for $x \in V$ , $0 < λ (x) \leq i_{b}$ (Lemma 2.4), we obtain that

\frac{d (x_{b} (t) - x_{d} (t))}{d t} \leq \frac{i_{b} i_{d} v (x (t)) e^{- λ (x (t)) x_{d} (t)}}{i_{d} + λ (x (t))} (\frac{2 i_{d} + i_{b}}{i_{d}} e^{- i_{d} (x_{b} (t) - x_{d} (t))} - 1) .

From the previous inequality, we deduce that on the set

{t \geq 0 : x_{b} (t) - x_{d} (t) > \frac{\log (\frac{2 i_{d} + i_{b}}{i_{d}})}{i_{d}}},

the quantity $x_{b} (t) - x_{d} (t)$ is decreasing, which allows us to conclude.

(iii): As before and by (i), we assume without loss of generality that $(x (t), t \geq 0) \subset U_{2} \cup H$ . Using (ii) and since $λ (x) \leq i_{b}$ (Lemma 2.4), we obtain that

\begin{array}{ll} \frac{d x_{b} (t)}{d t} & = i_{b} v (x (t)) e^{- i_{d} (x_{b} (t) - x_{d} (t))} e^{- λ (x (t)) x_{b} (t)} \\ \geq i_{b} \underline{v} (x^{0}) e^{- C} e^{- i_{b} x_{b} (t)}, \end{array}

that allows to conclude that $x_{b} (t)$ increases to $+ \infty$ and by (ii) we also have a similar behavior for $x_{d} (t)$ . We now prove that $λ (x (t)) \to i_{b}$ as $t \to + \infty$ . Let us recall that for all $t \geq 0$ , $λ (x (t))$ is the unique solution of

i_{b} \int_{0}^{x_{b} (t)} e^{- i_{d} (a - x_{d} (t))_{+} - λ (x (t)) a} d a = 1,

that we rewrite

i_{b} \int_{0}^{x_{d} (t)} e^{- λ (x (t)) a} d a + i_{b} \frac{(e^{- λ (x (t)) x_{d} (t)} - e^{- i_{d} (x_{b} (t) - x_{d} (t)) - λ (x (t)) x_{b} (t)})}{i_{d} + λ (x (t))} = 1.

The map $t \mapsto λ (x (t))$ is clearly increasing (using Equation 7 and the positivity of $x_{b}^{'} (t)$ and $x_{d}^{'} (t)$ ) and bounded by i_b. So there exists $λ^{*} > 0$ such that $λ (x (t)) \to λ^{*}$ . By taking the limit $t \to + \infty$ in Equation 15 and using the previous part of the proof, we deduce that

i_{b} \int_{0}^{+ \infty} e^{- λ^{*} a} d a = 1,

and $λ^{*} = i_{b}$ that concludes the proof. □

We now prove Theorem 4.4.

Proof of Theorem 4.4

By Lemma 4.5 (i), we assume without loss of generality that $(x (t), t \geq 0) \subset U_{2} \cup H$ , i.e that for all $t \geq 0$ , $x_{b} (t) - x_{d} (t) \geq 0$ . We recall that Equality Equation 14 gives:

\frac{d x_{b} (t) - x_{d} (t)}{d t} = \frac{i_{b} i_{d} v (x (t)) e^{- λ (x (t)) x_{d} (t)}}{i_{d} + λ (x (t))} (\frac{2 i_{d} + λ (x (t))}{i_{d}} e^{- (i_{d} + λ (x (t)) (x_{b} (t) - x_{d} (t)))} - 1) .

We define $f, h : R_{+} \to R$ by:

f (t) = \frac{i_{b} i_{d} v (x (t)) e^{- λ (x (t)) x_{d} (t)}}{i_{d} + λ (x (t))}

and

h (t) = \frac{2 i_{d} + λ (x (t))}{i_{d}} e^{- (i_{d} + λ (x (t))) (x_{b} (t) - x_{d} (t))} - \frac{2 i_{d} + i_{b}}{i_{d}} e^{- (i_{d} + i_{b}) (x_{b} (t) - x_{d} (t))} .

Note that $h (t) \to 0$ as $t \to + \infty$ using Lemma 4.5 (ii). Let us also define $u (t) = x_{b} (t) - x_{d} (t)$ . So Equation 16 rewrites

\frac{d u (t)}{d t} = f (t) (\frac{2 i_{d} + i_{b}}{i_{d}} e^{- (i_{d} + i_{b}) u (t)} - 1 + h (t)) .

We deduce that for all $ϵ > 0$ , there exists $t_{0} > 0$ such that for all $t \geq t_{0}$ ,

f (t) (\frac{2 i_{d} + i_{b}}{i_{d}} e^{- 2 (i_{b} + i_{d}) u (t)} - 1 - ϵ) \leq \frac{d u (t)}{d t} \leq f (t) (\frac{2 i_{d} + i_{b}}{i_{d}} e^{- (i_{d} + i_{b}) u (t)} - 1 + ϵ) .

Let us consider the differential equation

\frac{d w (t)}{d t} = f (t) (\frac{2 i_{d} + i_{b}}{i_{d}} e^{- (i_{b} + i_{d}) w (t)} - 1 + ϵ) .

By using the change of variables $s = e^{(i_{b} + i_{d}) w}$ , we solve the previous equation and we find that there exists a constant $C (x^{0})$ such that

w (t) = \frac{1}{i_{b} + i_{d}} \log (\frac{2 i_{d} + i_{b}}{i_{d} (1 - ϵ)} - \frac{C (x^{0})}{1 - ϵ} \exp (- (i_{b} + i_{d}) (1 - ϵ) \int_{0}^{t} f (s) d s)) .

We conclude by proving that the integral above tends to infinity as $t$ tends to infinity. First, the inequality $x_{b} (t) \geq x_{d} (t)$ implies that

f (t) \geq \frac{i_{b} i_{d} v (x (t))}{i_{d} + λ (x (t))} e^{- λ (x (t)) x_{b} (t)} .

Moreover, Equation 11 gives that

v (x (t)) e^{- λ (x (t)) x_{b} (t)} = e^{i_{d} (x_{b} (t) - x_{d} (t))} x_{b}^{'} (t) .

Since $λ (x (t)) \leq i_{b}$ , we obtain that

f (t) \geq \frac{i_{b} i_{d}}{i_{b} + i_{d}} x_{b}^{'} (t)

and that

\int_{0}^{t} f (s) d s \geq \frac{i_{b} i_{d}}{i_{b} + i_{d}} x_{b} (t) \underset{t \to + \infty}{⟶} + \infty .

By Equation 17, we conclude that for all $ϵ > 0$ :

\frac{1}{i_{b} + i_{d}} \log (\frac{2 i_{d} + i_{b}}{i_{d} (1 + ϵ)}) \leq \underset{t \to + \infty}{lim inf} u (t) \leq \underset{t \to + \infty}{lim sup} u (t) \leq \frac{1}{i_{d} + i_{b}} \log (\frac{2 i_{d} + i_{b}}{i_{d} (1 - ϵ)})

that concludes the proof. □

Appendix 1—figure 3

Download asset Open asset

Optimal configurations as time tends to infinity.

5 On the selection of Lansing effect

In this section, we ask the question of the apparition of a pro-senescence and non-genetic mechanism similar to the Lansing effect (Lansing, 1954; Lansing, 1947). We recall that the Lansing effect is the effect through which the progeny of old parents do not live as long as those of young parents.

We will show that the Lansing effect can represent a selective advantage, as an accelerator of the evolution.

5.1 The bd model with Lansing effect

The bd-model with Lansing effect is defined by modifying the bd-model that we introduced in Section 1. It was introduced and studied in details in Méléard et al., 2019 in the case where $i_{b} = i_{d} = 1$ . The authors show that under the assumptions of the adaptive dynamics theory (large population, rare and small mutations), the evolution of the trait $(x_{b}, x_{d})$ is described by the solutions a differential inclusion which reach the diagonal ${(x_{b}, x_{d}) \in R_{+}^{2} : x_{b} = x_{d}}$ and then stay on it. The formula given here are generalized to the case where $i_{b} \neq i_{d}$ .

The model

We assume that an individual which reproduces after age x_d transmits to its descendant a shorter life-expectancy. If an individual with trait $x = (x_{b}, x_{d})$ reproduces at age $a$ , the trait of its descendant is determined by a two-phases mechanism. The first phase is non-genetic and modifies the trait $x$ : if $a < x_{d}$ we define $\tilde{x} = x$ but if $a > x_{d}$ , $\tilde{x} = (x_{b}, 0)$ . The second phase corresponds to genetic mutations which modify the trait $\tilde{x}$ similarly as in Section 1. Hence, on configurations ${(x_{b}, x_{d}) \in R_{+}^{2} : x_{b} < x_{d}}$ , the dynamics is similar as in the model described in Section 1. Let us note that the population is then composed of two subpopulations, a population with traits ${(x_{b}, x_{d}), x_{b} > 0, x_{d} > 0}$ and a population with traits ${(x_{b}, 0), x_{b} > 0}$ .

Demographic parameters

We now give for the model with Lansing effect, the analogous of the demographic parameters introduced in Section 2. We refer to Méléard et al., 2019 for the justification. We denote by $λ^{ℓ} (x)$ the Malthusian parameter describing the asymptotic growth of the population with Lansing effect. It is solution of

i_{b} \int_{0}^{x_{b} \land x_{d}} e^{- λ^{ℓ} (x) a} d a = 1.

Then it can been easily computed by Newton’s method (as seen in Section 2) and the set of viability $V_{ℓ}$ is simple. It is composed of the traits $x = (x_{b}, x_{d})$ such that

x_{b} \land x_{d} > \frac{1}{i_{b}} .

The associated stable age distribution $N_{x}^{ℓ} (a) = (N_{x}^{ℓ, 1} (a), N_{x}^{ℓ, 2} (a))$ satisfies

N_{x}^{ℓ, 1} (a) = i_{b} e^{- i_{d} (a - x_{d})_{+} - λ^{ℓ} (x) a}, N_{x}^{ℓ, 2} (a) = i_{b} F (λ^{ℓ} (x)) e^{- (i_{d} + λ^{ℓ} (x)) a},

where $F$ is some function that we don’t detail here (Méléard et al., 2019, Proposition 3.5). The functions $N_{x}^{ℓ, 1}$ and $N_{x}^{ℓ, 2}$ describe the stable age distributions for populations with traits $(x_{b}, x_{d})$ and $(x_{b}, 0)$ respectively. The generation time $G^{ℓ} (x)$ is given by

G^{ℓ} (x) = \int_{0}^{x_{b} \land x_{d}} a i_{b} e^{- i_{d} (a - x_{d})_{+} - λ^{ℓ} (x) a} d a = i_{b} \int_{0}^{x_{b} \land x_{d}} a e^{- λ^{ℓ} (x) a} d a .

We observe that the Malthusian parameter $λ^{ℓ} (x)$ and the mean generation time $G^{ℓ} (x)$ only take into account the individuals reproducing before age $x_{b} \land x_{d}$ .

Evolution of the trait with Lansing effect

Let us now describe the behaviour of the trait. On the subset ${x_{b} < x_{d}}$ , the Lansing effect doesn’t act. So, the dynamics is similar as the one described in the above sections. The trait dynamics is described by the differential equation

\frac{d x_{b} (t)}{dt} = \frac{\partial λ^{ℓ} (x)}{\partial x_{b}} \frac{λ^{ℓ} (x)}{2 i_{b} c \int_{R_{+}} (N_{x}^{ℓ, 1} + N_{x}^{ℓ, 2}) (a) d a} σ^{2} (x), \frac{d x_{d} (t)}{dt} = 0.

Thus, the trait x_b increases while the trait x_d stays constant. On the subset ${x_{b} > x_{d}}$ , only individuals breeding before the age x_d will have viable offspring. Thus, there is no selective advantage in extending the reproduction phase by increasing x_b, but only in increasing survival by increasing x_d. More precisely, on ${x_{b} < x_{d}}$ , we have:

\frac{d x_{d} (t)}{dt} = \frac{\partial λ^{ℓ} (x)}{\partial x_{d}} \frac{λ^{ℓ} (x)}{2 i_{b} c \int_{R_{+}} (N_{x}^{ℓ, 1} + N_{x}^{ℓ, 2}) (a) d a} σ^{2} (x), \frac{d x_{b} (t)}{dt} = 0.

Indeed, the derivatives of the fitness are given as follows (see Méléard et al., 2019 Proposition 4.1).

\forall x \in U_{1}, \nabla λ^{ℓ} (x) = (\frac{i_{b} e^{- λ^{ℓ} (x) x_{b}}}{G^{ℓ} (x)}, 0); \forall x \in U_{2}, \nabla λ^{ℓ} (x) = (0, \frac{i_{b} e^{- λ^{ℓ} (x) x_{d}}}{G^{ℓ} (x)}) .

We observe that the trait x_d increases while the trait x_b stays constant. Hence, whatever the initial condition, the trait $x$ reaches in finite time the diagonal ${x_{b} = x_{d}}$ and then stays on it. On this diagonal the trait can evolve at different speeds (the dynamics is not unique): the global behavior of the trait is described by a differential inclusion (Méléard et al., 2019, Theorem 4.17).

5.2 Selection for Lansing effect

Let us first note that for Non Lansing and Lansing populations, as observed in the study of adaptive dynamics, the long time strategy leads to traits x_b and x_d going to infinity, with $x_{b} - x_{d} = \log (1 + \frac{i_{b} + i_{d}}{i_{d}}) / (i_{b} + i_{d})$ in the Non Lansing case and $x_{b} = x_{d}$ in the Lansing case (see Méléard et al., 2019, Theorem 4.17 in that case). It is then easy to deduce that in both cases, the Malthusian parameter, which has been proved to be less than i_b, converges to i_b when $t$ tends to infinity. Therefore the evolution will give the same selective advantage to both populations, making possible the cohabitation of the two populations. In addition, we observe that the partial derivatives of the Malthusian parameters with respect to x_b or x_d (in both cases) are positive, meaning that the convergences are increasing. Let us consider a monotype population with trait $(x_{b}, x_{d}) \in U_{2}$ , then by definition, we obtain that

λ^{n ℓ} (x) > λ^{ℓ} (x)

at time 0. Thus, there are periods where the Lansing fitness will increase much more than the Non Lansing one.

In order to assess the relative evolutionary success of Non Lansing/Lansing populations, we consider a population composed of two sub-monomorphic populations with traits respectively $x^{ℓ} = (x_{b}^{ℓ}, x_{d}^{ℓ})$ and $x^{n ℓ} = (x_{b}^{n ℓ}, x_{d}^{n ℓ})$ , the first one subject to the Lansing effect and the second one which is not affected by this senescence effect, both subjected to the same competitive pressure. The traits have been chosen such that the two sub-populations have the same darwinian fitness $λ^{n ℓ} (x^{n ℓ}) = λ^{ℓ} (x^{ℓ})$ . In each sub-population, the dynamics is described either in Section 1 (without Lansing effect) or in Section 5.1 (with Lansing effect). Let us first note that since $λ^{n ℓ} (x^{n ℓ}) = λ^{ℓ} (x^{ℓ})$ and since by definition,

i_{b} \int_{0}^{x_{d}^{ℓ}} e^{- λ^{ℓ} (x^{ℓ}) a} d a = 1 = i_{b} \int_{0}^{x_{d}^{n ℓ}} e^{- λ^{n ℓ} (x^{n ℓ}) a} d a + i_{b} \int_{x_{d}^{n ℓ}}^{x_{b}^{n ℓ}} e^{- i_{d} (a - x_{d}^{n ℓ})} e^{- λ^{n ℓ} (x^{n ℓ}) a} d a,

we deduce immediately that

x_{b}^{n ℓ} > x_{d}^{ℓ} > x_{d}^{n ℓ} .

We observe the isoclines of $λ^{n ℓ}$ and $λ^{ℓ}$ when they have the same values. Although they are very simple (horizontal or vertical lines) in the Lansing case, and in the region U₁ for the non-Lansing case, they have a more complicated form in the region U₂ for the non-Lansing case (Figure 5 of the main paper).

Let us consider the points $x^{n ℓ} \in U_{2}$ such that $λ^{n ℓ} (x^{n ℓ})$ has a fixed constant value. Using the Implicit Function Theorem, we know the existence of a real-valued smooth function $φ_{n ℓ}$ such that for all these points, $x_{d}^{n ℓ} = φ_{n ℓ} (x_{b}^{n ℓ})$ . Further,

φ_{n ℓ}^{'} (x_{b}^{n ℓ}) = - \frac{\frac{\partial λ^{n ℓ}}{\partial x_{b}^{n ℓ}}}{\frac{\partial λ^{n ℓ}}{\partial x_{d}^{n ℓ}}} (x^{n ℓ}) .

The previous computations showed that the partial derivatives of $λ^{n ℓ}$ are positive, and then that $φ_{n ℓ}^{'} (x_{b}^{n ℓ}) < 0$ , yielding the function $φ_{ℓ}$ to be decreasing on U₂. Moreover, the exact computation gives

φ_{ℓ}^{'} (x_{b}^{n ℓ}) = - \frac{i_{b} e^{- i_{d} (x_{b}^{n ℓ} - x_{d}^{n ℓ})} e^{- λ (x^{n ℓ}) x_{b}^{n ℓ}}}{i_{b} i_{d} e^{i_{d} x_{d}^{n ℓ}} \int_{x_{d}^{n ℓ}}^{x_{b}^{n ℓ}} e^{- (i_{d} + λ^{n ℓ} (x^{n ℓ})) a} d a} \geq - \frac{1}{i_{d} (x_{b}^{n ℓ} - x_{d}^{n ℓ})} .

The last inequality explains the almost vertical tangent observed when $x^{n ℓ}$ is close to the diagonal (see Figure 5 of the main paper).

Data availability

All data are available for download at https://github.com/MichaelRera/EvoAgeing/tree/main/article_sims. Package IBMPopSim (R package IBMPopSim v0.3.1): https://cran.r-project.org/web/packages/IBMPopSim/index.html. Github repository for simulation results and code: https://github.com/MichaelRera/EvoAgeing/tree/main/article_sims (copy archived at Rera, 2024). Environment for simulations using IBMPopSim: https://mybinder.org/v2/gh/MichaelRera/EvoAgeing/HEAD. Exploring parameters for Lansing populations: Modele_Lansing_evo.ipynb. Exploring parameters for non-Lansing populations Modele_nonLansing_evo.ipynb. Lansing / non-Lansing competition for equal Malthusian parameters L_nL_compet_eqMalth.ipynb. Lansing / non-Lansing competition (xb-xd) € [-10; 10] L_nL_compet_heteroPop.ipynb.

References

1. Arslan RC
2. Willführ KP
3. Frans EM
4. Verweij KJH
5. Bürkner PC
6. Myrskylä M
7. Voland E
8. Almqvist C
9. Zietsch BP
10. Penke L
(2017) Older fathers’ children have lower evolutionary fitness across four centuries and in four populations
Proceedings. Biological Sciences 284:20171562.

https://doi.org/10.1098/rspb.2017.1562
- PubMed
- Google Scholar
(2003) Extended longevity in mice lacking the insulin receptor in adipose tissue
Science 299:572–574.

https://doi.org/10.1126/science.1078223
- PubMed
- Google Scholar
1. Burke MK
2. Barter TT
3. Cabral LG
4. Kezos JN
5. Phillips MA
6. Rutledge GA
7. Phung KH
8. Chen RH
9. Nguyen HD
10. Mueller LD
11. Rose MR
(2016) Rapid divergence and convergence of life-history in experimentally evolved Drosophila melanogaster
Evolution; International Journal of Organic Evolution 70:2085–2098.

https://doi.org/10.1111/evo.13006
- PubMed
- Google Scholar
Data
1. Cansell C
2. Goepp V
3. Bain F
4. Todd N
5. Douard V
6. Monnoye M
7. Sanchez C
8. Pietrancosta N
9. Rovere C
10. Denis RG
11. Luquet S
12. Rera M
(authors) (2023) Two phases model of ageing in mice: towards a better identification of age-related and late-life metabolic decline [Registered Report Stage 1 Protocol]
figshare.

https://doi.org/10.6084/m9.figshare.23963208.v1
Book
1. Charlesworth B
(1994) Evolution in Age-Structured Populations
Cambridge University Press.

https://doi.org/10.1017/CBO9780511525711
- Google Scholar
1. Clancy DJ
2. Gems D
3. Harshman LG
4. Oldham S
5. Stocker H
6. Hafen E
7. Leevers SJ
8. Partridge L
(2001) Extension of life-span by loss of CHICO, a Drosophila insulin receptor substrate protein
Science 292:104–106.

https://doi.org/10.1126/science.1057991
- PubMed
- Google Scholar
1. Dambroise E
2. Monnier L
3. Ruisheng L
4. Aguilaniu H
5. Joly JS
6. Tricoire H
7. Rera M
(2016) Two phases of aging separated by the Smurf transition as a public path to death
Scientific Reports 6:23523.

https://doi.org/10.1038/srep23523
- PubMed
- Google Scholar
1. Echave J
(2021) Evolutionary coupling range varies widely among enzymes depending on selection pressure
Biophysical Journal 120:4320–4324.

https://doi.org/10.1016/j.bpj.2021.08.042
- PubMed
- Google Scholar
1. Fabian D
(2011)
The evolution of aging

Nature Education Knowledge 3:9.
- Google Scholar
1. Ferrière R
2. Tran VC
(2009) Stochastic and deterministic models for age-structured populations with genetically variable traits
ESAIM: Proc 27:289–310.

https://doi.org/10.1051/proc/2009033
- Google Scholar
1. Frenk S
2. Houseley J
(2018) Gene expression hallmarks of cellular ageing
Biogerontology 19:547–566.

https://doi.org/10.1007/s10522-018-9750-z
- PubMed
- Google Scholar
1. Gavrilov LA
2. Gavrilova NS
(2002) Evolutionary theories of aging and longevity
TheScientificWorldJournal 2:339–356.

https://doi.org/10.1100/tsw.2002.96
- PubMed
- Google Scholar
Preprint
(2023) Efficient Simulation of Individual-Based Population Models: The R Package IBMPopSim
arXiv.

https://doi.org/10.48550/arXiv.2303.06183
- Google Scholar
1. Gurtin ME
2. Maccamy RC
(1974) Non-linear age-dependent population dynamics
Arch Rational Mech Anal 54:281–300.

https://doi.org/10.1007/BF00250793
- Google Scholar
Book
1. Haldane J
(1941)
New Paths In Genetics

London: George Allen and Unwin, Ltd.
- Google Scholar
1. Hamilton WD
(1966) The moulding of senescence by natural selection
Journal of Theoretical Biology 12:12–45.

https://doi.org/10.1016/0022-5193(66)90184-6
- PubMed
- Google Scholar
(2014) Diversity of ageing across the tree of life
Nature 505:169–173.

https://doi.org/10.1038/nature12789
- PubMed
- Google Scholar
1. Kenyon C
2. Chang J
3. Gensch E
4. Rudner A
5. Tabtiang R
(1993) A C. elegans mutant that lives twice as long as wild type
Nature 366:461–464.

https://doi.org/10.1038/366461a0
- PubMed
- Google Scholar
1. Kirkwood TBL
2. Holliday R
(1979) The evolution of ageing and longevity
Proceedings of the Royal Society of London. Series B, Biological Sciences 205:531–546.

https://doi.org/10.1098/rspb.1979.0083
- PubMed
- Google Scholar
1. Kirkwood TBL
(2005) Understanding the odd science of aging
Cell 120:437–447.

https://doi.org/10.1016/j.cell.2005.01.027
- PubMed
- Google Scholar
1. Kirschner M
2. Gerhart J
(1998) Evolvability
PNAS 95:8420–8427.

https://doi.org/10.1073/pnas.95.15.8420
- PubMed
- Google Scholar
1. Kowald A
2. Kirkwood TBL
(2016) Can aging be programmed? A critical literature review
Aging Cell 15:986–998.

https://doi.org/10.1111/acel.12510
- PubMed
- Google Scholar
1. Lansing AI
(1947) A transmissible, cumulative, and reversible factor in aging
Journal of Gerontology 2:228–239.

https://doi.org/10.1093/geronj/2.3.228
- PubMed
- Google Scholar
1. Lansing AI
(1954) A nongenic factor in the longevity of rotifers
Annals of the New York Academy of Sciences 57:455–464.

https://doi.org/10.1111/j.1749-6632.1954.tb36418.x
- PubMed
- Google Scholar
1. Lemaître JF
2. Berger V
3. Bonenfant C
4. Douhard M
5. Gamelon M
6. Plard F
7. Gaillard JM
(2015) Early-late life trade-offs and the evolution of ageing in the wild
Proceedings. Biological Sciences 282:20150209.

https://doi.org/10.1098/rspb.2015.0209
- PubMed
- Google Scholar
1. Lemoine M
(2021) The evolution of the hallmarks of aging
Frontiers in Genetics 12:693071.

https://doi.org/10.3389/fgene.2021.693071
- PubMed
- Google Scholar
1. Luckinbill LS
2. Clare MJ
(1985) Selection for life span in Drosophila melanogaster
Heredity 55 (Pt 1):9–18.

https://doi.org/10.1038/hdy.1985.66
- PubMed
- Google Scholar
1. Maynard Smith J
(1971) What use is sex?
Journal of Theoretical Biology 30:319–335.

https://doi.org/10.1016/0022-5193(71)90058-0
- Google Scholar
1. Medawar PB
(1953) An unsolved problem of biology
Medical Journal of Australia 1:854–855.

https://doi.org/10.5694/j.1326-5377.1953.tb84985.x
- Google Scholar
1. Méléard S
2. Tran VC
(2009) Trait Substitution Sequence process and Canonical Equation for age-structured populations
Journal of Mathematical Biology 58:881–921.

https://doi.org/10.1007/s00285-008-0202-2
- PubMed
- Google Scholar
(2019) A birth–death model of ageing: from individual-based dynamics to evolutive differential inclusions
Journal of Mathematical Biology 79:901–939.

https://doi.org/10.1007/s00285-019-01382-z
- Google Scholar
(2020) Intergenerational transfer of ageing: parental age and offspring lifespan
Trends in Ecology & Evolution 35:927–937.

https://doi.org/10.1016/j.tree.2020.07.005
- PubMed
- Google Scholar
(2018) Experimental demonstration that offspring fathered by old males have shorter telomeres and reduced lifespans
Proceedings. Biological Sciences 285:20180268.

https://doi.org/10.1098/rspb.2018.0268
- PubMed
- Google Scholar
Book
1. Perthame B
(2007) Transport Equations in Biology
New York: Springer.

https://doi.org/10.1007/978-3-7643-7842-4
- Google Scholar
(2002) The role of parental age effects on the evolution of aging
Evolution 56:927–935.

https://doi.org/10.1111/j.0014-3820.2002.tb01405.x
- Google Scholar
1. Rera M
2. Bahadorani S
3. Cho J
4. Koehler CL
5. Ulgherait M
6. Hur JH
7. Ansari WS
8. Lo T Jr
9. Jones DL
10. Walker DW
(2011) Modulation of longevity and tissue homeostasis by the Drosophila PGC-1 homolog
Cell Metabolism 14:623–634.

https://doi.org/10.1016/j.cmet.2011.09.013
- PubMed
- Google Scholar
(2012) Intestinal barrier dysfunction links metabolic and inflammatory markers of aging to death in Drosophila
PNAS 109:21528–21533.

https://doi.org/10.1073/pnas.1215849110
- Google Scholar
(2018) The Smurf transition: new insights on ageing from end-of-life studies in animal models
Current Opinion in Oncology 30:38–44.

https://doi.org/10.1097/CCO.0000000000000419
- PubMed
- Google Scholar
Software
1. Rera M
(2024) EvoAgeing, version swh:1:rev:a9900ab85b63ef21cbf03f612d318f5d7b816399
Software Heritage.

https://archive.softwareheritage.org/swh:1:dir:b2dcbdd92130e6ee1d9285f136f7e646af5e9176;origin=https://github.com/MichaelRera/EvoAgeing;visit=swh:1:snp:85d0365cc037dc68e039d495a19109db25eb7c8c;anchor=swh:1:rev:a9900ab85b63ef21cbf03f612d318f5d7b816399
1. Rodrigues MA
2. Flatt T
(2016) Endocrine uncoupling of the trade-off between reproduction and somatic maintenance in eusocial insects
Current Opinion in Insect Science 16:1–8.

https://doi.org/10.1016/j.cois.2016.04.013
- Google Scholar
1. Rose MR
2. Drapeau MD
3. Yazdi PG
4. Shah KH
5. Moise DB
6. Thakar RR
7. Rauser CL
8. Mueller LD
(2002) Evolution of late-life mortality in Drosophila melanogaster
Evolution; International Journal of Organic Evolution 56:1982–1991.

https://doi.org/10.1111/j.0014-3820.2002.tb00124.x
- PubMed
- Google Scholar
(2000) Selection on age at reproduction in Drosophila melanogaster: female mating frequency as a correlated response
Evolution; International Journal of Organic Evolution 54:2152–2155.

https://doi.org/10.1111/j.0014-3820.2000.tb01257.x
- PubMed
- Google Scholar
1. Stearns SC
(1989) Trade-Offs in Life-History Evolution
Functional Ecology 3:259.

https://doi.org/10.2307/2389364
- Google Scholar
(2000) Experimental evolution of aging, growth, and reproduction in fruitflies
PNAS 97:3309–3313.

https://doi.org/10.1073/pnas.97.7.3309
- PubMed
- Google Scholar
1. Tricoire H
2. Rera M
(2015) A new, discontinuous 2 phases of aging model: lessons from Drosophila melanogaster
PLOS ONE 10:e0141920.

https://doi.org/10.1371/journal.pone.0141920
- PubMed
- Google Scholar
1. Trifonov EN
(2011) Vocabulary of definitions of life suggests a definition
Journal of Biomolecular Structure & Dynamics 29:259–266.

https://doi.org/10.1080/073911011010524992
- PubMed
- Google Scholar
(2005) Reduced insulin/IGF-1 signalling and human longevity
Aging Cell 4:79–85.

https://doi.org/10.1111/j.1474-9728.2005.00148.x
- PubMed
- Google Scholar
Book
1. Webb GF
(1985)
Theory of Nonlinear Age-Dependent Population Dynamics

CRC Press.
- Google Scholar
Book
1. Weismann A
(1882)
Ueber Die Dauer Des Lebens; Ein Vortrag

Jena: G. Fischer.
- Google Scholar
1. Williams GC
(1957) Pleiotropy, natural selection, and the evolution of senescence
Evolution 11:398.

https://doi.org/10.2307/2406060
- Google Scholar
Preprint
1. Zane F
2. Bouzid H
3. Marmol SS
4. Brazane M
5. Besse S
6. Molina JL
7. Cansell C
8. Aprahamian F
9. Durand S
10. Ayache J
11. Antoniewski C
12. Todd N
13. Rera M
(2023) Smurfness-Based Two-Phase Model of Ageing Helps Deconvolve the Ageing Transcriptional Signature
bioRxiv.

https://doi.org/10.1101/2022.11.22.517330
- Google Scholar

Article and author information

Author details

Tristan Roget

Institut Montpelliérain Alexander Grothendieck (IMAG), Université de Montpellier, Montpellier, France

Contribution
Software, Formal analysis, Investigation

Competing interests
No competing interests declared
Claire Macmurray

Université Paris Cité, Paris, France

Contribution
Writing - review and editing

Competing interests
No competing interests declared
Pierre Jolivet

Sorbonne Université, CNRS, Paris, France

Contribution
Software

Competing interests
No competing interests declared
Sylvie Meleard

Institut Universitaire de France et École Polytechnique, CNRS, Institut polytechnique de Paris, Palaiseau, France

Contribution
Conceptualization, Formal analysis, Supervision, Funding acquisition, Investigation, Methodology, Writing - original draft

Competing interests
No competing interests declared
Michael Rera

Université Paris Cité, Institut Jacques Monod, Paris, France

Contribution
Conceptualization, Resources, Data curation, Software, Supervision, Funding acquisition, Investigation, Methodology, Writing - original draft, Project administration, Writing - review and editing

For correspondence
michael.rera@cnrs.fr

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-6574-6511

Funding

Agence Nationale de la Recherche (ANR-20-CE44-0010)

Michael Rera

ATIP/Avenir

Michael Rera

TGCC@CEA (A0090607519)

Pierre Jolivet

Veolia Environnement - Ecole Polytechnique - Museum d'Histoire Naturelle - Fondation X (Chair "Modélisation Mathématique et Biodiversité")

Sylvie Meleard

European Research Council (SINGER 101054787)

Sylvie Meleard

Centre National de la Recherche Scientifique

Michael Rera

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank Dr. Sarah Kaakai for her help in transposing our initial simulation Python codes into the IBMPopSim (Giorgi et al., 2023) framework, Dr. Allon Weiner for a few hours of 'naive' discussion that helped explore the interpretations of this model’s impact on our perception of ageing, Dr. André Klarsfeld for his numerous useful comments on the manuscript. This work was granted access to the GENCI-sponsored HPC resources of TGCC@CEA under allocation A0090607519. This work has been supported by the Chair 'Modélisation Mathématique et Biodiversité' of Veolia Environnement-Ecole Polytechnique-Museum National d’Histoire Naturelle-Fondation X. SM is funded by the European Union (ERC, SINGER, 101054787). Views and opinions expressed are however those of the author(s) only and do not necessarily reflect those of the European Union or the European Research Council. Neither the European Union nor the granting authority can be held responsible for them. MR is funded by the CNRS. This project was partially funded by the ANR ADAGIO (ANR-20-CE44-0010) and the ATIP/Avenir young group leader program for MR.

Version history

Preprint posted: September 19, 2023
Sent for peer review: October 9, 2023
Reviewed Preprint version 1: January 5, 2024
Reviewed Preprint version 2: April 2, 2024
Version of Record published: November 1, 2024

Cite all versions

You can cite all versions using the DOI https://doi.org/10.7554/eLife.92914. This DOI represents all versions, and will always resolve to the latest one.

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

1,568

views
90

downloads
1

citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Mendeley

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Tristan Roget
Claire Macmurray
Pierre Jolivet
Sylvie Meleard
Michael Rera

(2024)

A scenario for an evolutionary selection of ageing

eLife 13:RP92914.

https://doi.org/10.7554/eLife.92914.3

Share this article

Cite this article

Three typical configurations of the model with ib >id and their effect on progeny’s genotypes as a function of parental age.

The bd model shows a convergence of xb - xd towards a positive value.

Populations with Lansing effect are favorably selected under logistic competition when the mutation rate is non-zero.

Table 1—source data 1

The Lansing effect maximizes populational survival by increasing its evolvability.

Mixed populations lead to (xb - xd) theoretical limit in a limited time and cohabitation of Lansing and non-Lansing populations.

The Lansing effect is associated with an increased fitness gradient.

The set V={(xb,xd)∈R+2;R(xb,xd)>1} is the convex set delimited by the black curve with equation R(xb,xd)=1.

Simulation of the canonical equation with x0=(3.5,1.3) and ib=id=1.

Optimal configurations as time tends to infinity.

Author details

Tristan Roget

Contribution

Competing interests

Claire Macmurray

Contribution

Competing interests

Pierre Jolivet

Contribution

Competing interests

Sylvie Meleard

Contribution

Competing interests

Michael Rera

Contribution

For correspondence

Competing interests

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Categories and tags

Further reading

Three typical configurations of the model with i_b >i_d and their effect on progeny’s genotypes as a function of parental age.

The bd model shows a convergence of x_b - x_d towards a positive value.

Mixed populations lead to (x_b - x_d) theoretical limit in a limited time and cohabitation of Lansing and non-Lansing populations.

The set $V = {(x_{b}, x_{d}) \in R_{+}^{2}; R (x_{b}, x_{d}) > 1}$ is the convex set delimited by the black curve with equation $R (x_{b}, x_{d}) = 1$ .

Simulation of the canonical equation with $x^{0} = (3.5, 1.3)$ and $i_{b} = i_{d} = 1$ .