Dynamics of adaptation in an anisotropic phenotype-fitness landscape

F. Hamel (1) , Florian Lavigne (1, 2, 3, 4) , Guillaume Martin (5) , L. Roques (2)

1 I2M - Institut de Mathématiques de Marseille
2 BioSP - Biostatistique et Processus Spatiaux
3 UMR ISEM - Institut des Sciences de l'Evolution de Montpellier
4 LMRS - Laboratoire de Mathématiques Raphaël Salem
5 UMR ISEM - Institut des Sciences de l'Evolution de Montpellier

We consider a diffusive 'replicator-mutator' equation to describe the adaptation of a large asexual population in a n-dimensional phenotypic space, under anisotropic mutation and selection effects. Though this equation has been studied previously when n = 1 or under isotropy assumptions, the n-dimensional anisotropic case remained unexplored. We prove here that the equation admits a unique solution, which is interpreted as the phenotype distribution, and we propose a new and general framework to the study of the quantitative behavior of this solution. When the initial distribution is Gaussian, the equation admits an explicit solution which remains normally distributed at all times, albeit with dynamic mean vector and variance-covariance matrix. In the general case, we derive a degenerate nonlocal parabolic equation satisfied by the distribution of the 'fitness components', and a nonlocal transport equation satisfied by the cumulant generating function of the joint distribution of these components. This last equation can be solved analytically and we then get a general formula for the trajectory of the mean fitness and all higher cumulants of the fitness distribution, over time. Such mean fitness trajectory is the typical outcome of empirical studies of adaptation by experimental evolution, and can thus be compared to empirical data. In sharp contrast with the known results based on isotropic models, our results show that the trajectory of mean fitness may exhibit (n−1) plateaus before it converges. It may thus appear 'non-saturating' for a transient but possibly long time, even though a phenotypic optimum exists. To illustrate the empirical relevance of these results, we show that the anisotropic model leads to a very good fit of Escherichia coli long-term evolution experiment, one of the most famous experimental dataset in experimental evolution. The two 'evolutionary epochs' that have been observed in this experiment have long puzzled the community: we propose that the pattern may simply stem form a climbing hill process, but in an anisotropic fitness landscape.

Mots clés: Adaptation; Mutation; Nonlocal equation; Degenerate PDE; Anisotropy

Domaines: Equations aux dérivées partielles [math.AP]; Génétique des populations [q-bio.PE]

F. Hamel, Florian Lavigne, Guillaume Martin, L. Roques. Dynamics of adaptation in an anisotropic phenotype-fitness landscape. Nonlinear Analysis: Real World Applications, 2020, 54 (103107), pp.103107. ⟨10.1016/j.nonrwa.2020.103107⟩. ⟨hal-02115130⟩