Horizontal gene transfer: numerical comparison between stochastic and deterministic approaches

Horizontal gene Transfer (HT) denotes the transmission of genetic material between two living organisms, while the vertical transmission refers to a DNA transfer from parents to their offspring. Consistent experimental evidence report that this phenomenon plays an essential role in the evolution of certain bacterias. In particular, HT is believed to be the main instrument of developing the antibiotic resistance. In this work, we consider several models which describe this phenomenon: a stochastic jump process (individual-based) and the deterministic nonlinear integrod-ifferential equation obtained as a limit for large populations. We also consider a Hamilton-Jacobi equation, obtained as a limit of the deterministic model under the assumption of small mutations. The goal of this paper is to compare these models with the help of numerical simulations. More specifically, our goal is to understand to which extent the Hamilton-Jacobi model reproduces the qualitative behavior of the stochastic model and the phenomenon of evolutionary rescue in particular.

[1]  Sepideh Mirrahimi,et al.  A class of Hamilton-Jacobi equations with constraint: uniqueness and constructive approach , 2015 .

[2]  G. Barles.  Solutions de viscosité des équations de Hamilton-Jacobi , 1994 .

[3]  V. Tran,et al.  Stochastic dynamics of adaptive trait and neutral marker driven by eco-evolutionary feedbacks , 2013, Journal of mathematical biology.

[4]  Dexi Liu,et al.  Advances in Gene Delivery Systems , 2011, Pharmaceutical Medicine.

[5]  G. Barles,et al.  Dirac concentrations in Lotka-Volterra parabolic PDEs , 2007, 0708.3720.

[6]  P. Collet,et al.  The effect of competition and horizontal trait inheritance on invasion, fixation, and polymorphism. , 2015, Journal of theoretical biology.

[7]  Viet Chi Tran,et al.  Stochastic and deterministic models for age-structured populations with genetically variable traits , 2008, 0809.3767.

[8]  Nicolas Fournier,et al.  A microscopic probabilistic description of a locally regulated population and macroscopic approximations , 2004, math/0503546.

[9]  J. Kingman.  On the genealogy of large populations , 1982, Journal of Applied Probability.

[10]  Axel Klar.  An Asymptotic Preserving Numerical Scheme for Kinetic Equations in the Low Mach Number Limit , 1999 .

[11]  P. Lions.  Generalized Solutions of Hamilton-Jacobi Equations , 1982 .

[12]  R. Ferrière,et al.  From Individual Stochastic Processes to Macroscopic Models in Adaptive Evolution , 2008 .

[13]  D. Posada,et al.  Coalescent Simulation of Intracodon Recombination , 2010, Genetics.

[14]  B. Perthame,et al.  Dirac Mass Dynamics in Multidimensional Nonlocal Parabolic Equations , 2010, 1011.1768.

[15]  A. Klar.  An Asymptotic-Induced Scheme for Nonstationary Transport Equations in the Diffusive Limit , 1998 .

[16]  B. Perthame,et al.  The dynamics of adaptation: an illuminating example and a Hamilton-Jacobi approach. , 2005, Theoretical population biology.

[17]  D. Posada,et al.  Simulation of Genome-Wide Evolution under Heterogeneous Substitution Models and Complex Multispecies Coalescent Histories , 2014, Molecular biology and evolution.

[18]  Shi Jin,et al.  Efficient Asymptotic-Preserving (AP) Schemes For Some Multiscale Kinetic Equations , 1999, SIAM J. Sci. Comput..

[19]  Stochastic dynamics for adaptation and evolution of microorganisms , 2016, 1610.00983.

[20]  J. Hadfield,et al.  Horizontally Transmitted Symbionts and Host Colonization of Ecological Niches , 2013, Current Biology.

[21]  N. Britton,et al.  The Persistence of Parasitic Plasmids , 2007, Genetics.

[22]  P. Lions,et al.  Two approximations of solutions of Hamilton-Jacobi equations , 1984 .