-
Phase transition in random adaptive walks on correlated fitness landscapes
Authors:
Su-Chan Park,
Ivan G. Szendro,
Johannes Neidhart,
Joachim Krug
Abstract:
We study biological evolution on a random fitness landscape where correlations are introduced through a linear fitness gradient of strength $c$. When selection is strong and mutations rare the dynamics is a directed uphill walk that terminates at a local fitness maximum. We analytically calculate the dependence of the walk length on the genome size $L$. When the distribution of the random fitness…
▽ More
We study biological evolution on a random fitness landscape where correlations are introduced through a linear fitness gradient of strength $c$. When selection is strong and mutations rare the dynamics is a directed uphill walk that terminates at a local fitness maximum. We analytically calculate the dependence of the walk length on the genome size $L$. When the distribution of the random fitness component has an exponential tail we find a phase transition of the walk length $D$ between a phase at small $c$ where walks are short $(D \sim \ln L)$ and a phase at large $c$ where walks are long $(D \sim L)$. For all other distributions only a single phase exists for any $c > 0$. The considered process is equivalent to a zero temperature Metropolis dynamics for the random energy model in an external magnetic field, thus also providing insight into the aging dynamics of spin glasses.
△ Less
Submitted 15 April, 2015; v1 submitted 20 August, 2014;
originally announced August 2014.
-
Adaptation in tunably rugged fitness landscapes: The Rough Mount Fuji Model
Authors:
Johannes Neidhart,
Ivan G. Szendro,
Joachim Krug
Abstract:
Much of the current theory of adaptation is based on Gillespie's mutational landscape model (MLM), which assumes that the fitness values of genotypes linked by single mutational steps are independent random variables. On the other hand, a growing body of empirical evidence shows that real fitness landscapes, while possessing a considerable amount of ruggedness, are smoother than predicted by the M…
▽ More
Much of the current theory of adaptation is based on Gillespie's mutational landscape model (MLM), which assumes that the fitness values of genotypes linked by single mutational steps are independent random variables. On the other hand, a growing body of empirical evidence shows that real fitness landscapes, while possessing a considerable amount of ruggedness, are smoother than predicted by the MLM. In the present article we propose and analyse a simple fitness landscape model with tunable ruggedness based on the Rough Mount Fuji (RMF) model originally introduced by Aita et al. [Biopolymers 54:64-79 (2000)] in the context of protein evolution. We provide a comprehensive collection of results pertaining to the topographical structure of RMF landscapes, including explicit formulae for the expected number of local fitness maxima, the location of the global peak, and the fitness correlation function. The statistics of single and multiple adaptive steps on the RMF landscape are explored mainly through simulations, and the results are compared to the known behavior in the MLM model. Finally, we show that the RMF model can explain the large number of second-step mutations observed on a highly-fit first step backgound in a recent evolution experiment with a microvirid bacteriophage [Miller et al., Genetics 187:185-202 (2011)].
△ Less
Submitted 1 August, 2014; v1 submitted 13 February, 2014;
originally announced February 2014.
-
Exact Results for Amplitude Spectra of Fitness Landscapes
Authors:
Johannes Neidhart,
Ivan G. Szendro,
Joachim Krug
Abstract:
Starting from fitness correlation functions, we calculate exact expressions for the amplitude spectra of fitness landscapes as defined by P.F. Stadler [J. Math. Chem. 20, 1 (1996)] for common landscape models, including Kauffman's NK-model, rough Mount Fuji landscapes and general linear superpositions of such landscapes. We further show that correlations decaying exponentially with the Hamming dis…
▽ More
Starting from fitness correlation functions, we calculate exact expressions for the amplitude spectra of fitness landscapes as defined by P.F. Stadler [J. Math. Chem. 20, 1 (1996)] for common landscape models, including Kauffman's NK-model, rough Mount Fuji landscapes and general linear superpositions of such landscapes. We further show that correlations decaying exponentially with the Hamming distance yield exponentially decaying spectra similar to those reported recently for a model of molecular signal transduction. Finally, we compare our results for the model systems to the spectra of various experimentally measured fitness landscapes. We claim that our analytical results should be helpful when trying to interpret empirical data and guide the search for improved fitness landscape models.
△ Less
Submitted 15 July, 2013; v1 submitted 9 January, 2013;
originally announced January 2013.
-
Quantitative analyses of empirical fitness landscapes
Authors:
Ivan G. Szendro,
Martijn F. Schenk,
Jasper Franke,
Joachim Krug,
J. Arjan G. M. de Visser
Abstract:
The concept of a fitness landscape is a powerful metaphor that offers insight into various aspects of evolutionary processes and guidance for the study of evolution. Until recently, empirical evidence on the ruggedness of these landscapes was lacking, but since it became feasible to construct all possible genotypes containing combinations of a limited set of mutations, the number of studies has gr…
▽ More
The concept of a fitness landscape is a powerful metaphor that offers insight into various aspects of evolutionary processes and guidance for the study of evolution. Until recently, empirical evidence on the ruggedness of these landscapes was lacking, but since it became feasible to construct all possible genotypes containing combinations of a limited set of mutations, the number of studies has grown to a point where a classification of landscapes becomes possible. The aim of this review is to identify measures of epistasis that allow a meaningful comparison of fitness landscapes and then apply them to the empirical landscapes to discern factors that affect ruggedness. The various measures of epistasis that have been proposed in the literature appear to be equivalent. Our comparison shows that the ruggedness of the empirical landscape is affected by whether the included mutations are beneficial or deleterious and by whether intra- or intergenic epistasis is involved. Finally, the empirical landscapes are compared to landscapes generated with the Rough Mt.\ Fuji model. Despite the simplicity of this model, it captures the features of the experimental landscapes remarkably well.
△ Less
Submitted 17 October, 2012; v1 submitted 20 February, 2012;
originally announced February 2012.
-
Rare events in population genetics: Stochastic tunneling in a two-locus model with recombination
Authors:
Alexander Altland,
Andrej Fischer,
Joachim Krug,
Ivan G. Szendro
Abstract:
We study the evolution of a population in a two-locus genotype space, in which the negative effects of two single mutations are overcompensated in a high fitness double mutant. We discuss how the interplay of finite population size, $N$, and sexual recombination at rate $r$ affects the escape times $t_\mathrm{esc}$ to the double mutant. For small populations demographic noise generates massive flu…
▽ More
We study the evolution of a population in a two-locus genotype space, in which the negative effects of two single mutations are overcompensated in a high fitness double mutant. We discuss how the interplay of finite population size, $N$, and sexual recombination at rate $r$ affects the escape times $t_\mathrm{esc}$ to the double mutant. For small populations demographic noise generates massive fluctuations in $t_\mathrm{esc}$. The mean escape time varies non-monotonically with $r$, and grows exponentially as $\ln t_{\mathrm{esc}} \sim N(r - r^\ast)^{3/2}$ beyond a critical value $r^\ast$.
△ Less
Submitted 22 February, 2011; v1 submitted 8 December, 2010;
originally announced December 2010.
-
1/f^beta noise in a model for weak ergodicity breaking
Authors:
Markus Niemann,
Ivan G. Szendro,
Holger Kantz
Abstract:
In systems with weak ergodicity breaking, the equivalence of time averages and ensemble averages is known to be broken. We study here the computation of the power spectrum from realizations of a specific process exhibiting 1/f^beta noise, the Rebenshtok-Barkai model. We show that even the binned power spectrum does not converge in the limit of infinite time, but that instead the resulting value is…
▽ More
In systems with weak ergodicity breaking, the equivalence of time averages and ensemble averages is known to be broken. We study here the computation of the power spectrum from realizations of a specific process exhibiting 1/f^beta noise, the Rebenshtok-Barkai model. We show that even the binned power spectrum does not converge in the limit of infinite time, but that instead the resulting value is a random variable stemming from a distribution with finite variance. However, due to the strong correlations in neighboring frequency bins of the spectrum, the exponent beta can be safely estimated by time averages of this type. Analytical calculations are illustrated by numerical simulations.
△ Less
Submitted 15 July, 2010; v1 submitted 4 February, 2010;
originally announced February 2010.
-
Structure of characteristic Lyapunov vectors in spatiotemporal chaos
Authors:
Diego Pazó,
Ivan G. Szendro,
Juan M. López,
Miguel A. Rodríguez
Abstract:
We study Lyapunov vectors (LVs) corresponding to the largest Lyapunov exponents in systems with spatiotemporal chaos. We focus on characteristic LVs and compare the results with backward LVs obtained via successive Gram-Schmidt orthonormalizations. Systems of a very different nature such as coupled-map lattices and the (continuous-time) Lorenz `96 model exhibit the same features in quantitative…
▽ More
We study Lyapunov vectors (LVs) corresponding to the largest Lyapunov exponents in systems with spatiotemporal chaos. We focus on characteristic LVs and compare the results with backward LVs obtained via successive Gram-Schmidt orthonormalizations. Systems of a very different nature such as coupled-map lattices and the (continuous-time) Lorenz `96 model exhibit the same features in quantitative and qualitative terms. Additionally we propose a minimal stochastic model that reproduces the results for chaotic systems. Our work supports the claims about universality of our earlier results [I. G. Szendro et al., Phys. Rev. E 76, 025202(R) (2007)] for a specific coupled-map lattice.
△ Less
Submitted 4 August, 2008;
originally announced August 2008.
-
Dynamics of perturbations in disordered chaotic systems
Authors:
Ivan G. Szendro,
Juan M. Lopez,
Miguel A. Rodriguez
Abstract:
We study the time evolution of perturbations in spatially extended chaotic systems in the presence of quenched disorder. We find that initially random perturbations tend to exponentially localize in space around static pinning centers that are selected by the particular configuration of disorder. The spatial structure of typical perturbations, $δu(x,t)$, is analyzed in terms of the Hopf-Cole tra…
▽ More
We study the time evolution of perturbations in spatially extended chaotic systems in the presence of quenched disorder. We find that initially random perturbations tend to exponentially localize in space around static pinning centers that are selected by the particular configuration of disorder. The spatial structure of typical perturbations, $δu(x,t)$, is analyzed in terms of the Hopf-Cole transform, $h(x,t) \equiv\ln|δu(x,t)|$. Our analysis shows that the associated surface $h(x,t)$ self-organizes into a faceted structure with scale-invariant correlations. Scaling analysis of critical roughening exponents reveals that there are three different universality classes for error propagation in disordered chaotic systems that correspond to different symmetries of the underlying disorder. Our conclusions are based on numerical simulations of disordered lattices of coupled chaotic elements and equations for diffusion in random potentials. We propose a phenomenological stochastic field theory that gives some insights on the path for a generalization of these results for a broad class of disordered extended systems exhibiting space-time chaos.
△ Less
Submitted 3 April, 2008;
originally announced April 2008.
-
Spatiotemporal structure of Lyapunov vectors in chaotic coupled-map lattices
Authors:
Ivan G. Szendro,
Diego Pazó,
Miguel A. Rodríguez,
Juan M. López
Abstract:
The spatiotemporal dynamics of Lyapunov vectors (LVs) in spatially extended chaotic systems is studied by means of coupled-map lattices. We determine intrinsic length scales and spatiotemporal correlations of LVs corresponding to the leading unstable directions by translating the problem to the language of scale-invariant growing surfaces. We find that the so-called 'characteristic' LVs exhibit…
▽ More
The spatiotemporal dynamics of Lyapunov vectors (LVs) in spatially extended chaotic systems is studied by means of coupled-map lattices. We determine intrinsic length scales and spatiotemporal correlations of LVs corresponding to the leading unstable directions by translating the problem to the language of scale-invariant growing surfaces. We find that the so-called 'characteristic' LVs exhibit spatial localization, strong clustering around given spatiotemporal loci, and remarkable dynamic scaling properties of the corresponding surfaces. In contrast, the commonly used backward LVs (obtained through Gram-Schmidt orthogonalization) spread all over the system and do not exhibit dynamic scaling due to artifacts in the dynamical correlations by construction.
△ Less
Submitted 20 September, 2007; v1 submitted 12 June, 2007;
originally announced June 2007.
-
Universal Critical Behavior of the Synchronization Transition in Delayed Chaotic Systems
Authors:
Ivan G. Szendro,
Juan M. López
Abstract:
We numerically investigate the critical behavior of the synchronization transition of two unidirectionally coupled delayed chaotic systems. We map the problem to a spatially extended system to show that the synchronization transition in delayed systems exhibits universal critical properties. We find that the synchronization transition is absorbing and generically belongs to the universality clas…
▽ More
We numerically investigate the critical behavior of the synchronization transition of two unidirectionally coupled delayed chaotic systems. We map the problem to a spatially extended system to show that the synchronization transition in delayed systems exhibits universal critical properties. We find that the synchronization transition is absorbing and generically belongs to the universality class of the bounded Kardar-Parisi-Zhang equation, as occurs in the case of extended systems. We also argue that directed percolation critical behavior may emerge for systems with strong nonlinearities
△ Less
Submitted 30 March, 2005;
originally announced March 2005.
-
Scaling properties of growing noninfinitesimal perturbations in space-time chaos
Authors:
Juan M. Lopez,
Cristina Primo,
Miguel A. Rodriguez,
Ivan G. Szendro
Abstract:
We study the spatiotemporal dynamics of random spatially distributed noninfinitesimal perturbations in one-dimensional chaotic extended systems. We find that an initial perturbation of finite size $ε_0$ grows in time obeying the tangent space dynamic equations (Lyapunov vectors) up to a characteristic time $t_{\times}(ε_0) \sim b - (1/λ_{max}) \ln (ε_0)$, where $λ_{max}$ is the largest Lyapunov…
▽ More
We study the spatiotemporal dynamics of random spatially distributed noninfinitesimal perturbations in one-dimensional chaotic extended systems. We find that an initial perturbation of finite size $ε_0$ grows in time obeying the tangent space dynamic equations (Lyapunov vectors) up to a characteristic time $t_{\times}(ε_0) \sim b - (1/λ_{max}) \ln (ε_0)$, where $λ_{max}$ is the largest Lyapunov exponent and $b$ is a constant. For times $t < t_{\times}$ perturbations exhibit spatial correlations up to a typical distance $ξ\sim t^z$. For times larger than $t_{\times}$ finite perturbations are no longer described by tangent space equations, memory of spatial correlations is progressively destroyed and perturbations become spatiotemporal white noise. We are able to explain these results by mapping the problem to the Kardar-Parisi-Zhang universality class of surface growth.
△ Less
Submitted 2 August, 2004;
originally announced August 2004.