Gulbudak2019.2 - Heterogeneous viral strategies promote coexistence in virus-microbe systems (Chronic)

  public model
Model Identifier
Short description
This is a mathematical model describing describing the population dynamics of microbes infected by chronically infecting viruses.
Related Publication
  • Heterogeneous viral strategies promote coexistence in virus-microbe systems.
  • Gulbudak H, Weitz JS
  • Journal of theoretical biology , 2/ 2019 , Volume 462 , pages: 65-84 , PubMed ID: 30389532
  • Department of Mathematics, University of Louisiana at Lafayette, Lafayette, LA, USA. Electronic address:
  • Viral infections of microbial cells often culminate in lysis and the release of new virus particles. However, viruses of microbes can also initiate chronic infections, in which new viruses particles are released via budding and without cell lysis. In chronic infections, viral genomes may also be passed on from mother to daughter cells during division. The consequences of chronic infections for the population dynamics of viruses and microbes remains under-explored. In this paper we present a model of chronic infections as well as a model of interactions between lytic and chronic viruses competing for the same microbial population. In the chronic only model, we identify conditions underlying complex bifurcations such as saddle-node, backward and Hopfbifurcations, leading to parameter regions with multiple attractors and/or oscillatory behavior. We then utilize invasion analysis to examine the coupled nonlinear system of microbes, lytic viruses, and chronic viruses. In so doing we find unexpected results, including a regime in which the chronic virus requires the lytic virus for survival, invasion, and persistence. In this regime, lytic viruses decrease total cell densities, so that a subpopulation of chronically infected cells experience decreased niche competition. As such, even when chronically infected cells have a growth disadvantage, lytic viruses can, paradoxically, enable the proliferation of both chronically infected cells and chronic viruses. We discuss the implications of our results for understanding the ecology and long-term evolution of heterogeneous viral strategies.
Submitter of the first revision: Johannes Meyer
Submitter of this revision: Johannes Meyer
Modellers: Johannes Meyer

Metadata information

hasProperty (2 statements)
Mathematical Modelling Ontology Ordinary differential equation model
Gene Ontology transmission of virus

Curation status


Connected external resources

SBGN view in Newt Editor

Name Description Size Actions

Model files

Gulbudak2019.2.Chronic.xml SBML L2V4 Representation of Gulbudak2019.2 - Heterogeneous viral strategies promote coexistence in virus-microbe systems (Chronic) 33.43 KB Preview | Download

Additional files

Gulbudak2019.2.Chronic.cps COPASI file of Gulbudak2019.2 - Heterogeneous viral strategies promote coexistence in virus-microbe systems (Chronic) 62.46 KB Preview | Download
Gulbudak2019.2.Chronic.sedml SED-ML file of Gulbudak2019.2 - Heterogeneous viral strategies promote coexistence in virus-microbe systems (Chronic) 2.69 KB Preview | Download

  • Model originally submitted by : Johannes Meyer
  • Submitted: Nov 10, 2019 8:35:33 PM
  • Last Modified: Nov 10, 2019 8:35:33 PM
  • Version: 2 public model Download this version
    • Submitted on: Nov 10, 2019 8:35:33 PM
    • Submitted by: Johannes Meyer
    • With comment: Automatically added model identifier BIOMD0000000846
: Variable used inside SBML models

Species Initial Concentration/Amount

3.32E7 item

C14283 ; C14141
0.0 item

8.3E8 item
Reactions Rate Parameters
V_C => compartment*mu*V_C mu = 0.0866
C => compartment*d_tilde*C d_tilde = 0.05
=> S compartment*r*S*(1-N/K) r = 0.339; N = 8.3E8; K = 8.947E8
S + V_C => C compartment*phi_tilde*S*V_C phi_tilde = 5.0E-12
=> C compartment*r_tilde*C*(1-N/K) r_tilde = 0.2; N = 8.3E8; K = 8.947E8
=> V_C; C compartment*alpha*C alpha = 0.1
S => compartment*d*S d = 0.0416666666666667
Curator's comment:
(added: 10 Nov 2019, 20:35:25, updated: 10 Nov 2019, 20:35:25)
Reproduced plot of Figure E.1 (a) in the original publication. Model simulated using COPASI 4.24 (Build 197), plot produced using Wolfram Mathematica 11.3.