10.5061/DRYAD.HK64F
Pagliarini, Silvia
University of Verona
Korobeinikov, Andrei
Centre de Recerca Matemàtica
Data from: A mathematical model of marine bacteriophage evolution
Dryad
dataset
2018
virus dynamics
Viral evolution
Darwinian fitness
Phenotype space
Mathematical model
2018-02-07T17:58:13Z
en
https://doi.org/10.1098/rsos.171661
5311 bytes
1
CC0 1.0 Universal (CC0 1.0) Public Domain Dedication
To explore how particularities of a cell-virus system affects viral
evolution, we formulate a mathematical model of marine bacteriophage
evolution. The intrinsic simplicity of real-life phage-bacteria systems
allows to have a reasonably simple model. The model constructed in this
paper is based upon Beretta-Kuang model of bacteria-phage interaction.
Compared to the Beretta-Kuang model, the model assumes the existence of a
multitude of viral variants which correspond to continuously distributed
phenotypes. It is noteworthy that this model does not include any explicit
law or mechanism of evolution; instead it is assumed, in agreement to the
principles of Darwinian evolution, that evolution in this system can occur
as a result of random mutations and natural selection. Simulations with a
leaner fitness landscape (which is chosen for the convenience of
demonstration only) show that a pulse-type traveling wave moving towards
increasing Darwinian fitness appears in the phenotype space. This implies
that the overall fitness of a viral quasispecies steadily increasing in
time. That is, the simulations demonstrate that for an uneven fitness
landscape random mutations combined with a mechanism of natural selection
lead to the Darwinian evolution. It is noteworthy that in this system the
speed of propagation of this wave (and hence the rate of evolution) is not
constant but varies, depending on the current viral fitness and the
abundance of susceptible bacteria.
Model of bacteriophage evolutionImplementation of the model in Matlab. The
model is a 3-equations model, combined with appropriated boundary and
initial conditions (based on biological background and on the
bibliography). Then, the solution is derived using an explicit
method.BPvect.mBPvect_RobinMatlab code. Evolution of marine bacteriophage
modeled by a 3-equation system, solved explicitely applying Robin boundary
conditions for infected cell population and Neumann boundary conditions on
the virus.