 Research article
 Open Access
 Published:
Effect of stochasticity on coinfection dynamics of respiratory viruses
BMC Bioinformatics volume 20, Article number: 191 (2019)
Abstract
Background
Respiratory viral infections are a leading cause of mortality worldwide. As many as 40% of patients hospitalized with influenzalike illness are reported to be infected with more than one type of virus. However, it is not clear whether these infections are more severe than single viral infections. Mathematical models can be used to help us understand the dynamics of respiratory viral coinfections and their impact on the severity of the illness. Most models of viral infections use ordinary differential equations (ODE) that reproduce the average behavior of the infection, however, they might be inaccurate in predicting certain events because of the stochastic nature of viral replication cycle. Stochastic simulations of single virus infections have shown that there is an extinction probability that depends on the size of the initial viral inoculum and parameters that describe viruscell interactions. Thus the coinfection dynamics predicted by the ODE might be difficult to observe in reality.
Results
In this work, a continuoustime Markov chain (CTMC) model is formulated to investigate probabilistic outcomes of coinfections. This CTMC model is based on our previous coinfection model, expressed in terms of a system of ordinary differential equations. Using the Gillespie method for stochastic simulation, we examine whether stochastic effects early in the infection can alter which virus dominates the infection.
Conclusions
We derive extinction probabilities for each virus individually as well as for the infection as a whole. We find that unlike the prediction of the ODE model, for similar initial growth rates stochasticity allows for a slower growing virus to outcompete a faster growing virus.
Background
With the advent of molecular diagnostic techniques, respiratory tract specimens from patients with influenzalike illness (ILI) are now being recognized as having multiple viruses [1–4]. Around 40% of the hospitalized patients with ILI have coinfections with influenza A virus (IAV), influenza B virus (IBV), respiratory syncytial virus (RSV), human rhinovirus (hRV), adenovirus (AdV), human enterovirus (hEV), human metapneumovirus (hMPV), coronavirus (CoV), parainfluenza virus (PIV), human bocavirus (hBoV) and many others [5–9]. These patients are reported to suffer from heterogeneous disease outcomes such as enhanced [10–12], reduced [13, 14] and unaltered [14–16] severity compared to patients with single virus infections. However, it is not clear how the virusvirus and virushost interactions influence disease severity and lead to these varied outcomes. Two or more virus agents can interact in diverse ways which may arise from the consequences of their inoculation order, inter exposure time, initial inoculums, different combinations of viruses, number of coinfecting viruses and host immune state [17, 18]. Thus, coinfections pose a combinatorial problem which can be challenging to study in a laboratory set up alone.
Coinfection can be better understood using mathematical modeling. While mathematical modeling of single virus infections at the cellular level has proven crucial for finding answers where laboratory experiments are impossible, impractical or expensive [19–23], little has been done in viral coinfection modeling. A few studies [24–26] have used within host models considering interactions of two different strains of the same virus. Among them, Pinilla et al. [24] and Petrie et al. [25] used their models to study competitive mixedinfection experiments of pandemic A/H1N1 influenza with its H275Y mutant strain and Simeonov et al. [26] considered a spatiotemporal model to explain in vitro cellular susceptibility due to the simultaneous presence of RSV A2 and RSV B. Pinky and Dobrovolny [27] proposed a two virus coinfection model to investigate viral interference observed in an experimental study of IAVRSV coinfection (Shinjoh et al. [28]) where they concluded that distinct viruses interact through resource competition. In further investigations [29, 30], they used the model to quantify the impact of resource availability, finding the possibility of chronic single infection if constant cellular regeneration was considered and chronic coinfection if both cellular regeneration and superinfection were considered. However, the majority of the two virus models studied so far have focused on the deterministic approach that reproduces the average behavior of infection kinetics. The exceptions are Dobrovolny et al. [31] and Deecke et al. [32] who investigated two strains of the same virus (wildtype and drug resistant mutant) using a stochastic model to determine mechanisms driving the emergence of drugresistant mutants during the course of a single infection. Since in real life viral infections are stochastic and discrete events, stochastic simulations of infection models will provide further insight into coinfection dynamics.
For example, stochastic simulations of single virus infections have shown that there is an extinction probability that depends on the size of the initial viral inoculum and parameters that describe viruscell interactions [33]. Similarly, experimental studies of viral infections in animals have shown that viruses do not always establish infection in every animal under study [34]. Although the causative phenomenon is still unidentified, there are some possible factors suggested by researchers such as host defense mechanisms, spatial heterogeneity in the target cell population, and the stochastic nature of the virus life cycle [34]. Moreover, evaluation of this quantity can be useful in many situations where the viral dynamics cannot be explained with a simple deterministic model. Numerous stochastic models have been developed to study various aspects of the single viral infection process such as virus release strategies (i.e. budding and bursting) for HIV [33, 34], impact of initial viral dose [35], length of eclipse and infectious phases [33, 34], impact of the immune response [34, 35], and how ongoing proliferation of immune cells acts to decrease the emergence probability of mutated strains [36]. These models have been studied using Monte Carlo simulations of the multitype branching process [37, 38], or by simulating solutions to stochastic differential equations where processes involved in the virus life cycle are diffusion processes (stochasticity is represented by noise terms in the equations) [35, 39, 40].
Of particular interest for stochastic models is the probability of extinction, a feature that ODE models cannot capture. In stochastic models, analytic expression of extinction probability is formulated by keeping track of the number of infected cells [41], the number of virions [42] and both [33] in single virus models, mostly for HIV infection. Yan et al. [34] used a similar method to calculate the extinction probability that includes time dependent immune responses in a single influenza virus model. Stochastic extinction could be a factor in coinfection dynamics since one virus could have a higher extinction probability, even if the two viruses have the same initial viral inoculum or initial growth rate, making it possible for one virus to go extinct while the other viral infection grows. Thus the coinfection outcomes predicted by the ODE model might be difficult to observe in reality.
In this work, we implement a stochastic counterpart of our previously published ODE coinfection model [27], in the form of a continuoustime Markov chain (CTMC) model. Trajectories for the CTMC model are simulated using Gillespie’s tauleap algorithm. In order to investigate how stochastic effects early in the infection impact coinfection, we vary the initial growth rate and compare to predictions from the ODE model. We also derive the extinction coefficient analytically for the model using multitype branching method. While the ODE model found that the virus with a higher growth rate consumes more target cells and produces higher peak viral load compared to the slower growing virus, we find that stochasticity can allow slower growing viruses to consume more target cells and produce more virus than the faster growing virus.
Results
Derivation of extinction coefficient
Stochastic extinction is most relevant during the initial stage of infection. At this stage the number of target cells is small. We can consider that the target cells are constant or equal to the initial number of target cells (T ≈ T_{0}). As a result, the states become decoupled making the stochastic events independent of each other. In addition, each event produces progeny over a lifetime which is also independent of the lifetime of all other events. More details on how to derive a branching process from a CTMC can be found in [43]. Under these conditions, the CTMC model becomes a multitype branching process where the reduced state vectors now represent \(\vec {m}\)=\((n_{E_{1}}\), \(n_{I_{1}}, n_{V_{1}}\), \(n_{E_{2}}, n_{I_{2}}, n_{V_{2}})\), where \(n_{E_{1}}\) and \(n_{E_{2}}\) are the numbers of eclipse cells, \(n_{I_{1}}\) and \(n_{I_{2}}\) are the infected cells, and \(n_{V_{1}}\) and \(n_{V_{2}}\) are the virions of both viruses. Including the assumption of a constant number of target cells, the reduced model is
Thus the continuoustime Markov chain becomes a multitype branching process that describes the dynamics of a population of individuals having birth and death independently according to the specified (in this case exponential) probability mass function. If a timehomogeneous CTMC is a branching process, the only absorbing state is \(\vec {0}\). For this model we defined the absorbing state as \(\vec {0}\) and the probability to reach this state from, say, \(\vec {m}\), is \(\xi (\vec {m})\). This probability is referred to as the extinction probability. Biologically, the extinction probability is defined as the probability that the two types of viruses and all their infected cells are completely eliminated from the host. Once a transition occurs, the current state \(\vec {m}\) is incremented by one of the transition vectors given below.
If the rate of the i^{th} reaction is defined as a_{i} such that a_{1}=β_{1}TV_{1},a_{2}=β_{2}TV_{2}, a_{3}=k_{1}E_{1},a_{4}=k_{2}E_{2}, a_{5}=δ_{1}I_{1},a_{6}=δ_{2}I_{2}, a_{7}=p_{1}I_{1},a_{8}=p_{2}I_{2}, a_{9}=c_{1}V_{1},a_{10}=c_{2}V_{2}, then the probability that the i^{th} reaction is the next reaction is given by
and n_{max} is the number of transitions involved in the model and is equal to 10. The time of the next reaction is a random variable with distribution \(Z(\vec {m}) \exp (Z(\vec {m})t)\) with mean \(\frac {1}{Z(\vec {m})}\) (according to the Gillespie algorithm). The probability that a simultaneous exposure to both viruses eventually evolves to extinction, or reaches the absorbing state, (0,0,0,0,0,0), from state \(\vec {m}\) or the extinction coefficient, \(\xi (\vec {m})\), is
Substituting the expressions for \(P_{i}(\vec {m})\) and \(\xi (\vec {m}+d\vec {m_{i}})\) in Eq. (1), the extinction coefficient becomes:
Although the general solution of this expression is intractable, the CTMC assumption of independent events means that the functional equation of \(\xi (\vec {m})\) can be reduced to an algebraic equation. Thus the extinction probability from a given state is the product of the extinction probabilities from each of the constituents of that state [44], so we can write
where \(\rho ^{n_{V_{1}}}_{V_{1}}\) is the probability that the virus, V_{1}, initiates a process with \(n_{V_{1}}\) number of virus particles that results in extinction. In a similar manner, \(\rho ^{n_{E_{1}}}_{E_{1}}, \rho ^{n_{I_{1}}}_{I_{1}}\) and others are probabilities for eclipse cell, E_{1}, or infected cell, I_{1} and so on. Eq. (3) is recognizable as the fixedpoint equation \(\vec {\varepsilon }= P(\vec {\varepsilon })\), where \(\vec {\varepsilon }= [\varepsilon _{1}, \ldots, \varepsilon _{J}]\) and \(P(\vec {\varepsilon })\) is the probability generating function of the progeny distributions. Now substituting Eq. (3) in Eq. (2), we get
where \(\rho _{V_{i}}\), \(\rho _{I_{i}}\) and \(\rho _{E_{i}}\) are the extinction probabilities when the processes are initiated with a single virus particle or eclipse cell or infectious cell. Solving for each probability, we get \(\rho _{V_{i}}=1\) and \(\rho _{V_{i}}= \frac {c_{i}(p_{i}+\delta _{i})}{p_{i}(c_{i}+\beta _{i}T)}\), \(\rho _{I_{i}}=1\) and \(\rho _{I_{i}}=\frac {\delta _{i}(c_{i}+\beta _{i}T)}{\beta _{i}T(p_{i}+\delta _{i})}\), and \(\rho _{E_{i}}=\rho _{I_{i}}\) where i=1,2. Since probability has to be less than or equal to 1, we can write the solutions of the extinction probabilities as:
Probability of virus extinction Since the extinction of each event is independent, we can write for the probability that both viruses go extinct if the simultaneous infection is initiated with a single virus of each type by the expression \(\rho _{V_{1}}\rho _{V_{2}}\),
Stochastic dynamics of identical viruses
While the probability of virus extinction is an important feature of stochastic models, we are also interested in understanding if stochasticity affects the predicted dynamics of coinfections that survive. Previously in our ODE model [27], we found that the virus with the higher growth rate always outcompetes the slower growing virus. While ODEs can give us the average behaviors of the coinfection process, in real systems the biological processes are stochastic. The randomness associated with births and deaths during the initial infection process may lead to virus extinction even in an exponentially growing virus population [45]. Yan et al. [34] reported that the invasion of viral infection is dependent on the initial viral dose and growth rate of each virus. Here, we are interested in knowing how the coinfection dynamics change with change in growth rates of each virus. First, we will observe the dynamics of coinfection with identical viruses.
Keeping all the initial conditions and transition rates for both viruses equal, we examine the time course of coinfection by plotting number of virus over time. 1000 sample stochastic trajectories of viral load curve for coinfection with identical viruses are shown in Fig. 1. We find that both viruses have peaks above the threshold of detection (100 virions) 88% of the time and 12% of the time one of the viruses experiences extinction. Among 120 (12%) extinctions, virus 1 and virus 2 experience extinction 49 and 65 times out of 1000 simulations, respectively. In another words, there is a 4.9% chance that what begins as a coinfection will result in a single virus infection with virus 2 or 6.5% chance with virus 1.
The ODE model predicts that when all parameters are equal, both viruses will have the same time course, splitting available target cells equally. In the stochastic model, we find that despite having identical growth rates, one virus outcompetes the other virus in particular realizations of the model. Virus 1 has a higher peak viral titer 513 times within 1000 simulations, while virus 2 has the higher peak viral titer 487 times. So while a particular realization of the model will have a clear dominant virus, on average, the viruses are equivalent, in agreement with the ODE model. Additional file 1 includes additional figures examining the distributions when viruses differ. To characterize the viral time course, we calculate peak viral load, time of peak for each virus, as well as duration of coinfection (Fig. 2). The median time of peak for virus 1 is found to be 2.374 ±0.64 days and for virus 2, it is 2.375 ±0.65 days. The median of peak viral load for virus 1 and 2 are (4.0 ±2.6) ×10^{7} and (4.1 ±2.6) ×10^{7}, respectively. From the distributions (Fig. 2), we see that even if the viruses behave differently for a particular model realization, on average, they tend to behave identically. Finally, the distribution of coinfection duration is given in Fig. 2 where the median coinfection duration is found to be 5.730 ±0.059 days. Despite fluctuations in the time course of each virus, the coinfection duration does not vary much.
Stochastic dynamics for different viruses
Since growth rate determines which virus is the stronger competitor [27], we investigate how differences in growth rate between the two viruses changes stochastic infections. Unfortunately, growth rate is not a parameter in the model, so we need to determine which model parameter(s) to change to systematically vary the growth rate. We use the expression for growth rate derived by Smith et al. [46] and determine how growth rate depends on different model parameters (Fig. 3). We find that growth rate varies approximately linearly with virus production rate, p, over a large range of p (p>1), so we will systematically alter p for one virus to alter its growth rate.
For ease of interpretation, we define the relative viral production rate \(r = \frac {p_{1}}{p_{2}}\). We first examine how competition between the viruses changes as the relative growth rates change. Here variation is introduced for virus 1 keeping virus 2 fixed for a range, r=1×10^{−1}×10^{2}. We count the number of times, out of 1000 simulated infections, a particular virus has a higher viral titer peak than the other virus. The results are shown in Fig. 4. When the viruses have identical growth rates, there is a 50% chance that a particular virus will have the higher peak titer, as seen in the previous section. The probability of having a higher peak viral load increases rapidly as the production rate of a virus increases, reaching 90% with a less than 2fold change in viral production. Note that the probability of having the higher peak viral titer never quite reaches 100%, even when there are large differences in growth rate. This indicates that early stochastic events can significantly alter the time course of the infection.
In Fig. 5, we compare coinfection dynamics for the ODE and CTMC models, looking at peak viral load, time of viral peak, and coinfection duration. ODEs predict that if the growth rate of one virus is higher than the other it will always have a higher peak viral load (Fig. 5 (top left)). For the CTMC model, the transition from one virus dominating to the other dominating is not as sharp. Unlike the predictions of ODEs, the CTMC allows for the slower growing virus to dominate infection dynamics. In fact, the median peak viral loads for virus 1 and virus 2 cross closer to a relative viral production rate of 10^{1} rather than 10^{0} as seen in the ODE model. Stochastic variability in the peak viral load (as indicated by the shaded area) for both viruses overlaps for a wide range of relative viral production, indicating that the viruses can have similar peak viral loads.
The time of viral peak also shows some differences between the ODE and CTMC models. For the ODE model, the time of viral peak is similar for both viruses when the relative viral production rate is greater than 10^{0}, although the time of peak decreases as the relative viral production rate increases. This is because the viral production rate of virus 1 is increased over its baseline value causing an earlier time of peak. This drives the earlier time of peak of virus 2, which is the weaker competitor in this case. The decline in time of viral peak is not as sharp in the CTMC model since stochasticity can temper the effect of the increased production rate of virus 2 by allowing virus 1 to still have an opportunity to infect some cells.
Finally, we compare the predicted duration of coinfection variation for ODE and stochastic models (Fig. 5 (bottom row)). Viruses do not coexist for more than about a week in either model. The longest coinfection durations are seen, for both models, when the two viruses have the same growth rates. This is because the faster growing virus outcompetes the slower growing virus leading to short infections for the slower growing virus.
One feature of viral infections that cannot be captured by ODE models is extinction of the infection. Therefore, we simulate the probability of virus extinction, defined as the fraction of times when one virus does not grow above the virus detection limit (detection limit is equal to 100 virus particles), when the coinfection is initiated with a single virus of each type (Fig. 6). Note that this is slightly different than the definition for the probability of extinction calculated in “Derivation of extinction coefficient” section which requires that virus, along with infectious and eclipse cells, all go to zero. The probability of both viruses growing to detectable levels is highest for viruses with similar relative production rates. When relative viral production rates are very different (about 10–100 fold difference), there is a high probability that one virus becomes extinct. When the viruses have very different production rates, the virus with a higher production rate will outcompete the one with a low production rate driving it to extinction. However, since one virus (in this case virus 1) experiences decreased production rate from the base value but initiates infection with the same amount of virus, the probability of extinction reaches close to 100% in a quicker manner for lower relative production rate than that of the higher relative rates.
Discussion
The dynamics of coinfection were previously modeled deterministically in several studies [24, 25, 29]. However, ODE models do not capture the very earliest dynamics of infection where stochastic effects may play an important role. The stochastic model presented here indicates that stochastic effects can dramatically alter the time course of the infection. Our previous ODE coinfection model [27] could not distinguish between two identical/similar viruses, as the predicted time courses are identical. Simulations of the stochastic model, however, indicate that for a particular realization of the model, two identical viruses can have very different time courses, with ∼12% of infections initiated with two viruses resulting in infections with only one detectable virus. When viruses have different growth rates, the ODE model predicts that the virus with the higher growth rate will have a higher peak viral titer. This is not the case for the CTMC where early stochastic effects can allow a slower growing virus to infect more target cells than the faster growing virus, giving the slower virus a competitive advantage that continues over the course of the infection.
The ODE coinfection model resulted in a simple rule for determining which virus would be dominant in a coinfection — the virus with a higher growth rate. Replication of the slowergrowing virus is suppressed due to lack of accessibility to target cells. This simple rule suggests that we can easily determine which viruses will be suppressed in coinfections. For example, application of the ODE model to several respiratory viruses indicated that parainfluenza virus (PIV) replication is substantially reduced during coinfection with other respiratory viruses [27], suggesting that it should be difficult to detect PIV in coinfections. However, PIV is detected in coinfections from 30–80% of the time [15, 47–50]. Some of this unexpectedly high detection rate might be due to stochasticity. PIV detection in coinfection is, however, lower than what is observed for two identical viruses as described in the previous paragraph. The slow growth rate of PIV means that most viruses will outcompete PIV more often than viruses with identical growth rates.
Stochasticity also impacts our ability to use viral interference as a possible mechanism for treating or preventing more serious infections. If we cannot guarantee that a fastgrowing virus will suppress growth of a slowgrowing virus, then this strategy might be risky. For example, some have suggested using defective interfering particles (DIPs) as a possible method for blocking infections [51–55]. DIPs cannot replicate on their own, but have a high growth rate when fullyfunctioning virus is present. Our results indicate that even when there is a large difference in the viral growth rate, there is a nonzero probability that the slowergrowing virus (in this case the fullyfunctional virus) will rise to a higher peak than the faster growing virus, suggesting that use of DIPs for treatment will not be completely effective.
While our extension of the simple coinfection model has provided insight into how stochasticity might affect coinfections, this simple model does not capture all biological processes during the infection. More complex ODE models that include cell regeneration [29] and superinfection [30] have been proposed and reproduce a broader range of behaviors observed during viral coinfections. Stochastic versions of these models can also be developed in the future to examine how stochasticity affects behaviors such as chronic coinfections. Other limitations include the lack of an explicit immune response, which will likely increase the probability of extinction of the coinfection [34], and the inclusion of realistic delays to account for intracellular replication [56]. Despite these shortcomings, this stochastic implementation of a viral coinfection model has shown the extent of variability in the time course of coinfections when stochasticity is introduced.
Conclusions
While ODE models are useful for giving a broad picture of possible dynamical behaviors of infection, in reality each infection is distinct with disease outcome dependent on early stochastic events. This is particular important when considering interactions between viruses during coinfection since stochasticity can drive one or both viruses to extinction before the infection has time to take hold. Our models show that for viral coinfections, this sometimes leads to a less fit virus outcompeting a more fit virus.
Methods
Continuoustime Markov chain model
The previously [27] proposed ODE coinfection model considers the mean concentrations of viruses and cells in a large population. Here, we formulate the probabilistic counterpart of the ODE model, a timehomogeneous CTMC model of two competing viruses with particular account for stochastic effects in the early infection processes. This model considers variability in each viral replication event (for example, infectivity of target cells, transition to eclipse phase, activation of infectious phase and its lifespan, virus production and virus clearance) and takes values on a set of states collectively known as the state space Ω. The states of the full system are defined as \(\vec {m}=(n_{T}, n_{E_{1}}, n_{I_{1}}, n_{V_{1}}, n_{E_{2}}, n_{I_{2}}, n_{V_{2}})\) where the state vectors denote the integer numbers of target cells, eclipse cells, infected cells, virions for virus 1 and 2 respectively. The states are discrete, and the stochastic process is timehomogeneous. The CTMC model that we implement is similar to that of Pearson et al. [33]. Figure 7 illustrates the model diagram. The model is
where T is the number of susceptible target cells, E_{1} and E_{2} are the number of eclipse cells or noninfectious infected cells, I_{1} and I_{2} are the number of active infectious cells and V_{1} and V_{2} are the number of virus particles. Viruses of each type infect the target cells, which are limited, at infection rates, β_{1} and β_{2}. Each type of infected cell transitions into eclipse phases where viruses undertake intracellular processes for progeny virus production. After the time durations of \(\frac {1}{k_{1}}\) and \(\frac {1}{k_{2}}\), eclipse cells become productive infectious cells, I_{1} and I_{2}, which produce viruses at production rates p_{1} and p_{2} over their lifespans of \(\frac {1}{\delta _{1}}\) and \(\frac {1}{\delta _{2}}\) respectively. While viruses are being produced by the infectious cells, some of the cell free viruses are being cleared with clearance rates c_{1} and c_{2}. Thus the number of transitions involved in the model is 10 with respective transition rates (propensities) defined in Table 1. In our model assumption, we have ignored virus particle loss due to absorption into the cells since this amount is negligible compared to free virus particles produced. Finally, no specific immune interactions are considered in this model.
It has been shown that, the stochastic representations of chemical reactions converge to the differential equations as the number of particles goes to infinity when we can assume that the probability of a reaction depends on the density of the reactants [57–59]. We make a similar assumption for the “reactions” involved in viral replication where infection of a cell, for example, depends on the density of both cells and virus. Since biological processes, particularly at the microscopic level, are really a series of chemical reactions, there is an inherent stochasticity to the system that is not simply averaged out because we are not specifically considering the detailed chemical reactions in the model. For instance, the infection of a cell in this model includes binding of the virus to the cell receptor, fusion of the virus with the cell membrane, and opening of the virus membrane to release the contents, among other steps. These are all chemical reactions that can be assumed to occur with probability proportional to density of the reacting chemicals. It seems reasonable then to assume that the overall infection process is also dependent on the density of the larger entities (viruses and cells) that contain these chemicals and we can expect a similar convergence of the Markov chain to the differential equation when there are large numbers of viruses and cells.
Stochastic simulation algorithm
The direct method popularized by Gillespie [60], in general, is used for solving trajectories of timehomogeneous CTMC models. Bartlett [61] first applied this method for epidemic modeling of measles. However, since the computing time of the direct method scales linearly with the initial number of target population [34], the direct method becomes infeasible to simulate viral infection models with realistic number of target cells, i.e. of order 1×10^{8}. Due to the increased simulation efficiency with some accuracy trade off, Gillespie tauleap algorithms are getting more attention. In the tauleap method a small time interval is chosen such that the number of times each transition occurs in this interval is drawn from a Poisson distribution with mean equal to the expected number of transitions during that interval. The time step is fixed for the most basic tauleap method. However, the time step should be small enough that the rate at which transitions occur remains approximately constant during the fixed interval. For numerical implementation of the CTMC model, we use the Gillespie tauleaping method with time step equal to 10^{−3} day. The transitions involved in the stochastic process and their rates are summarized in Table 1. The parameter values for numerical simulation are taken from [19] and are given in Table 2. Viral load is usually measured as concentration with units such as TCID_{50}/mL, PFU/mL or EID_{50}/mL rather than as the total number of virus in the host, while the CTMC model uses discrete values for each of its state variables. In order to convert viral concentration measurements to the number of infectious virus particles, studies have used a conversion factor, α [31, 62, 63]. Although there is no standard value for α, previous estimates suggest that 1 TCID_{50}/mL of nasal wash corresponds to 1×10^{2}−1×10^{5} [62] or 3×10^{4}−3×10^{5}[63] virus particles at the site of infection. So we take α equal to 1×10^{2} to convert the concentrations of virus into numbers of virus particles according to the method mentioned in [31].
Abbreviations
 AdV:

Adenovirus
 CoV:

Coronavirus
 CTMC:

Continuous time Markov chain
 DIP:

Defective interfering particle
 EID:

Egg infectious dose
 hBoV:

Human bocavirus
 hEV:

Human enterovirus
 HIV:

Human immunodeficiency virus
 hMPV:

Human metapneumovirus
 hRV:

Human rhinovirus
 IAV:

Influenza A virus
 IBV:

Influenza B virus
 ILI:

Influenzalike illness
 ODE:

Ordinary differential equation
 PFU:

Plaque forming unit
 PIV:

Parainfluenza virus
 RSV:

Respiratory syncytial virus
 TCID:

Tissue culture infectious dose
References
 1
Goto H, Ihira H, Morishita K, Tsuchiya M, Ohta K, Yumine N, et al.Enhanced growth of influenza A virus by coinfection with human parainfluenza virus type 2. Med Microbiol Immunol. 2015; 205(3):209–18.
 2
Opatowski L, Baguelin M, Eggo RM. Influenza interaction with cocirculating pathogens and its impact on surveillance, pathogenesis, and epidemic profile: A key role for mathematical modelling. PLoS Pathog. 2018; 14(2):e1006770.
 3
Jartti T, SoderlundVenermo M, Hedman K, Ruuskanen O, Makela MJ. New molecular virus detection methods and their clinical value in lower respiratory tract infections in children. Paediatr Resp Rev. 2013; 14:38–45.
 4
Anderson TP, Werno AM, Barratt K, Mahagamasekera P, Murdoch DR, Jennings LC. Comparison of four multiplex PCR assays for the detection of viral pathogens in respiratory specimens. J Virol Meth. 2013; 191(2):118–21.
 5
Kenmoe S, Tchendjou P, Vernet MA, Tetang SM, Mossus T, Ripa MN, et al.Viral etiology of severe acute respiratory infections in hospitalized children in Cameroon, 2011–2013. Influenza and Other Res Vir. 2016; 10(5):386–393.
 6
Scotta M, Chakr V, de Moura A, Becker R, de Souza AP, Jones M, et al.Respiratory viral coinfection and disease severity in children: A systematic review and metaanalysis. J Clin Virol. 2016; 80:45–56.
 7
Asner SA, Rose W, Petrich A, Richardson S, Tran DJ. Is virus coinfection a predictor of severity in children with viral respiratory infections?Clin Microbiol Infect. 2015; 21(3):264.e1–6.
 8
Pretorius MA, Madhi SA, Cohen C, Naidoo D, Groome M, Moyes J, et al.Respiratory Viral Coinfections Identified by a 10Plex RealTime ReverseTranscription Polymerase Chain Reaction Assay in Patients Hospitalized With Severe Acute Respiratory IllnessSouth Africa, 20092010. J Infect Dis. 2012; 206(S1):S159–65.
 9
Achten NB, Wu P, Bont L, Blanken MO, Gebretsadik T, Chappell JD, et al.Interference Between Respiratory Syncytial Virus and Human Rhinovirus Infection in Infancy. J Infect Dis. 2017; 215:1102–6.
 10
Goka EA, Vallely PJ, Mutton KJ, Klapper PE. Single, dual and multiple respiratory virus infections and risk of hospitalization and mortality. Epidemiol Infect. 2015 January; 143(1):37–47.
 11
Semple MG, Cowell A, Dove W, Greensill J, McNamara PS, Halfhide C, et al.Dual infection of infants by human metapneumovirus and human respiratory syncytial virus is strongly associated with severe bronchiolitis. J Infect Dis. 2005; 191(3):382–6.
 12
Richard N, KomurianPradel F, Javouhey E, Perret M, Rajoharison A, Bagnaud A, et al.The impact of dual viral infection in infants admitted to a pediatric intensive care unit associated with severe bronchiolitis. Pediatr Infect Dis J. 2008; 27(3):213–7.
 13
Martin ET, Kuypers J, Wald A, Englund JA. Multiple versus single virus respiratory infections: viral load and clinical disease severity in hospitalized children. Influenza Other Respir Viruses. 2011 May; 6(1):71–77.
 14
Martin ET, Fairchok MP, Stednick ZJ, Kuypers J, Englund JA. Epidemiology of Multiple Respiratory Viruses in Childcare Attendees. J Infect Dis. 2013 January; 207:982–9.
 15
Brand HK, de Groot R, Galama JMD, Brouwer ML, Teuwen K, Hermans PWM, et al.Infection With Multiple Viruses is not Associated With Increased Disease Severity in Children With Bronchiolitis. Pediatr Pulmonol. 2012 April; 47(4):393–400.
 16
RotzenOstlund M, Eriksson M, Lindell AT, Allander T, Wirgart BZ, Grillner L. Children with multiple viral respiratory infections are older than those with single viruses. Acta Paediatr. 2014; 103(1):100–4.
 17
Smith AM. Hostpathogen kinetics during influenza infection and coinfection: insights from predictive modeling. Immunol Rev. 2018; 285(1):97–112.
 18
Kumar N, Sharma S, Barua S, Tripathi BN, Rouse BT. Virological and immunological outcomes of coinfections. Clin Microbiol Rev. 2018; 31(4):e00111–17.
 19
Baccam P, Beauchemin C, Macken CA, Hayden FG, Perelson AS. Kinetics of influenza A virus infection in humans. J Virol. 2006; 80(15):7590–9.
 20
Perelson AS, Neumann AU, Markowitz M, Leonard JM, Ho DD. HIV1 dynamics in vivo: virion clearance rate, infected cell lifespan, and viral generation time. Science. 1996; 271:1582–6.
 21
Neumann AU, Lam NP, Dahari H, Gretch DR, Wiley TE, Layden TJ, et al.Hepatitis C viral dynamics in vivo and the antiviral efficacy of interferonalpha therapy. Science. 1998; 282:103–7.
 22
GonzálezParra G, Dobrovolny HM, Aranda DF, ChenCharpentier B, Rojas RAG. Quantifying rotavirus kinetics in the REH tumor cell line using in vitro data. Virus Res. 2018; 244:53–63.
 23
GonzálezParra G, Dobrovolny HM. Assessing Uncertainty in A2 Respiratory Syncytial Virus Viral Dynamics. Comput Math Meth Med. 2015; 2015:567589.
 24
Pinilla LT, Holder BP, Abed Y, Boivin G, Beauchemin CAA. The H275Y Neuraminidase Mutation of the Pandemic A/H1N1 Influenza Virus Lengthens the Eclipse Phase and Reduces Viral Output of Infected Cells, Potentially Compromising Fitness in Ferrets. J Virol. 2012 October; 86(19):10651–60.
 25
Petrie SM, Butler J, Barr IG, McVernon J, Hurt AC, McCaw JM. Quantifying relative withinhost replication fitness in influenza virus competition experiments. J Theor Biol. 2015; 382(7):256–71.
 26
Simeonov I, Gong X, Kim O, Poss M, Chiaromonte F, Fricks J. Exploratory Spatial Analysis of in vitro Respiratory Syncytial Virus Coinfections. Virus. 2010; 2:2782–802.
 27
Pinky L, Dobrovolny HM. Coinfections of the Respiratory Tract: Viral Competition for Resources. PLoS One. 2016; 11(5):e0155589.
 28
Shinjoh M, Omoe K, Saito N, Matsuo N, Nerome K. In vitro growth profiles of respiratory syncytial virus in the presence of influenza virus. Acta Virol. 2000; 44(2):91–7.
 29
Pinky L, Dobrovolny HM. The impact of cell regeneration on the dynamics of viral coinfection. Chaos. 2017; 27(6):063109.
 30
Pinky L, GonzálezParra G, Dobrovolny HM. Superinfection and cell regeneration can lead to chronic viral coinfections. J Theor Biol. 2019 January; 466:24–38.
 31
Dobrovolny HM, Beauchemin CA. Modelling the emergence of influenza drug resistance: The roles of surface proteins, the immune response and antiviral mechanisms. PloS One. 2017; 12(7):e0180582.
 32
Deecke LAC, Dobrovolny HM. Intermittent treatment of severe influenza. J Theor Biol. 2018; 442:129–38.
 33
Pearson JE, Krapivsky P, Perelson AS. Stochastic Theory of Early Viral Infection: Continuous versus Burst Production of Virions. PLoS Comp Biol. 2011; 7(2):e1001058.
 34
Yan AWC, Cao P, McCaw JM. On the extinction probability in models of withinhost infection: the role of latency and immunity. J Math Biol. 2016; 73:787–813.
 35
Yuan Y, Allen LJ. Stochastic models for virus and immune system dynamics. Math Biosci. 2011; 234(2):84–94.
 36
Hartfield M, Alizon S. Withinhost stochastic emergence dynamics of immuneescape mutants. PLoS Comput Biol. 2015; 11(3):e1004149.
 37
Heffernan JM, Wahl LM. Monte Carlo estimates of natural variation in HIV infection. J Theo Bio. 2005; 236(2):137–53.
 38
Allen LJ, Lahodny GE. Extinction thresholds in deterministic and stochastic epidemic models. J Bio Dyna. 2012; 6(2):590–611.
 39
Kirupaharan N, Allen LJ. Coexistence of multiple pathogen strains in stochastic epidemic models with densitydependent mortality. Bull Math Bio. 2004; 66(4):841–64.
 40
Allen LJ. An introduction to stochastic processes with applications to biology. In: CRC Press. Boca Raton: Chapman and Hall/CRC: 2010.
 41
Merrill SJ. The stochastic dance of early HIV infection. J Com App Math. 2005; 184(1):242–57.
 42
Tuckwell HC, Shipman PD, Perelson AS. The probability of HIV infection in a new host and its reduction with microbicides. Math Biosci. 2008; 214(12):81–6.
 43
Pinky L. Viral coinfections of the respiratory tract. UMI Thesis. Fort Worth: TCU Press; 2018.
 44
Allen LJ, Jang SR, Roeger LI. Predicting population extinction or disease outbreaks with stochastic models. Lett Biomath. 2017; 4(1):1–22.
 45
Allen LJ, van den Driessche P. Relations between deterministic and stochastic thresholds for disease extinction in continuousand discretetime infectious disease models. Math Biosci. 2013; 243(1):99–108.
 46
Smith AM, Adler FR, Perelson AS. An accurate twophase approximate solution to an acute viral infection model. J Math Biol. 2010 May; 60(5):711–26.
 47
Cai XY, Wang Q, Lin GY, Cai ZW, Lin CX, Chen PZ, et al.Respiratory Virus Infections Among Children in South China. J Med Virol. 2014; 86:1249–55.
 48
Dierig A, Heron LG, Lambert SB, Yin JK, Leask J, Chow MYK, et al.Epidemiology of respiratory viral infections in children enrolled in a study of influenza vaccine effectiveness. Influenza Other Resp Viruses. 2014; 8(3):293–301.
 49
Martin ET, Kuypers J, Wald A, Englund JA. Multiple versus single virus respiratory infections: viral load and clinical disease severity in hospitalized children. Influenza Other Resp Viruses. 2011; 6(1):71–7.
 50
Zhang G, Hu Y, Wang H, Zhang L, Bao Y, Zhou X. High Incidence of Multiple Viral Infections Identified in Upper Respiratory Tract Infected Children under Three Years of Age in Shanghai, China. PLoS one. 2012; 7(9):e44568.
 51
Dimmock NJ, Easton AJ. Cloned Defective Interfering Influenza RNA and a Possible PanSpecific Treatment of Respiratory Virus Diseases. Viruses. 2015; 2015(7):3768–88.
 52
Zhao H, To KK, Chu H, Ding Q, Zhao X, Li C, et al.Dualfunctional peptide with defective interfering genes effectively protects mice against avian and seasonal influenza. Nature Communications. 2018; 15(9):2358.
 53
Wasik MA, Eichwald L, Genzel Y, Reichl U. Cell culturebased production of defective interfering particles for influenza antiviral therapy. Appl Microbiol Biotech. 2018; 102(3):1167–77.
 54
Baltes A, Akpinar F, Inankur B, Yin J. Inhibition of infection spread by cotransmitted defective interfering particles. Plos One. 2017; 12(9):e0184029.
 55
Smith CM, Scott PD, O’Callaghan C, Easton AJ, Dimmock NJA, Defective Interfering Influenza RNA. Inhibits Infectious Influenza Virus Replication in Human Respiratory Tract Cells: A Potential New Human Antiviral. Viruses. 2016; 8(8):237.
 56
Holder BP, Beauchemin CAA. Exploring the effect of biological delays in kinetic models of influenza within a host or cell culture. BMC Public Health. 2011; 11(S1):S10.
 57
Kurtz T. Limit theorems for sequences of jump Markov processes approximating ordinary differential processes. J Appl Prob. 1971; 8(2):344–56.
 58
Kurtz T. Reationship between stochastic and deterministic models for chemical reactions. J Chem Phys. 1972; 57(7):2976–8.
 59
Kurtz T. Solutions of ordinary differential equations as limits of pure jump Markov processes. J Appl Prob. 1970; 7(1):49–58.
 60
Gillespie DT. Exact stochastic simulation of coupled chemical reactions. J Phys Chem. 1977; 81(25):2340–61.
 61
Bartlett MS. Stochastic processes or the statistics of change. J R Stat Soc: Ser C: Appl Stat. 1953; 2(1):44–64.
 62
Handel A, Longini IM, Antia R. Neuraminidase inhibitor resistance in influenza: assessing the danger of its generation and spread. PLoS Comput Biol. 2007; 3(12):e240.
 63
Perelson AS, Rong L, Hayden FG. Combination antiviral therapy for influenza: predictions from modeling of human infections. Journal Infect Dis. 2012; 205(11):1642–5.
Acknowledgements
Not applicable.
Funding
Not applicable.
Availability of data and materials
All data is contained within the manuscript. All code is available at https://doi.org/http://github.com/hdobrovo/stochastic_coinfection.
Author information
Affiliations
Contributions
LP, GGP and HMD conceived the experiment(s), LP conducted the experiment(s), LP analyzed the results, LP and HMD wrote the manuscript. All authors reviewed the manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file
Additional file1
Pdf file containing additional figures. (PDF 26,531 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Pinky, L., GonzalezParra, G. & Dobrovolny, H.M. Effect of stochasticity on coinfection dynamics of respiratory viruses. BMC Bioinformatics 20, 191 (2019). https://doi.org/10.1186/s1285901927936
Received:
Accepted:
Published:
Keywords
 Viral coinfection
 Respiratory virus
 Withinhost model
 Continuoustime Markov chain
 Multitype branching process
 Extinction probability