The effect of demographic and environmental variability on disease outbreak for a dengue model with a seasonally varying vector population.

Seasonal changes in temperature, humidity, and rainfall affect vector survival and emergence of mosquitoes and thus impact the dynamics of vector-borne disease outbreaks. Recent studies of deterministic and stochastic epidemic models with periodic environments have shown that the average basic reproduction number is not sufficient to predict an outbreak. We extend these studies to time-nonhomogeneous stochastic dengue models with demographic variability wherein the adult vectors emerge from the larval stage vary periodically. The combined effects of variability and periodicity provide a better understanding of the risk of dengue outbreaks. A multitype branching process approximation of the stochastic dengue model near the disease-free periodic solution is used to calculate the probability of a disease outbreak. The approximation follows from the solution of a system of differential equations derived from the backward Kolmogorov differential equation. This approximation shows that the risk of a disease outbreak is also periodic and depends on the particular time and the number of the initial infected individuals. Numerical examples are explored to demonstrate that the estimates of the probability of an outbreak from that of branching process approximations agree well with that of the continuous-time Markov chain. In addition, we propose a simple stochastic model to account for the effects of environmental variability on the emergence of adult vectors from the larval stage.

Stochastic two-group models with transmission dependent on host infectivity or susceptibility.

Stochastic epidemic models with two groups are formulated and applied to emerging and re-emerging infectious diseases. In recent emerging diseases, disease spread has been attributed to superspreaders, highly infectious individuals that infect a large number of susceptible individuals. In some re-emerging infectious diseases, disease spread is attributed to waning immunity in susceptible hosts. We apply a continuous-time Markov chain (CTMC) model to study disease emergence or re-emergence from different groups, where the transmission rates depend on either the infectious host or the susceptible host. Multitype branching processes approximate the dynamics of the CTMC model near the disease-free equilibrium and are used to estimate the probability of a minor or a major epidemic. It is shown that the probability of a major epidemic is greater if initiated by an individual from the superspreader group or by an individual from the highly susceptible group. The models are applied to Severe Acute Respiratory Syndrome and measles.

The effect of delay in viral production in within-host models during early infection.

Delay in viral production may have a significant impact on the early stages of infection. During the eclipse phase, the time from viral entry until active production of viral particles, no viruses are produced. This delay affects the probability that a viral infection becomes established and timing of the peak viral load. Deterministic and stochastic models are formulated with either multiple latent stages or a fixed delay for the eclipse phase. The deterministic model with multiple latent stages approaches in the limit the model with a fixed delay as the number of stages approaches infinity. The deterministic model framework is used to formulate continuous-time Markov chain and stochastic differential equation models. The probability of a minor infection with rapid viral clearance as opposed to a major full-blown infection with a high viral load is estimated from a branching process approximation of the Markov chain model and the results are confirmed through numerical simulations. In addition, parameter values for influenza A are used to numerically estimate the time to peak viral infection and peak viral load for the deterministic and stochastic models. Although the average length of the eclipse phase is the same in each of the models, as the number of latent stages increases, the numerical results show that the time to viral peak and the peak viral load increase.

Probability of a disease outbreak in stochastic multipatch epidemic models.

Environmental heterogeneity, spatial connectivity, and movement of individuals play important roles in the spread of infectious diseases. To account for environmental differences that impact disease transmission, the spatial region is divided into patches according to risk of infection. A system of ordinary differential equations modeling spatial spread of disease among multiple patches is used to formulate two new stochastic models, a continuous-time Markov chain, and a system of stochastic differential equations. An estimate for the probability of disease extinction is computed by approximating the Markov chain model with a multitype branching process. Numerical examples illustrate some differences between the stochastic models and the deterministic model, important for prevention of disease outbreaks that depend on the location of infectious individuals, the risk of infection, and the movement of individuals.

Extinction thresholds in deterministic and stochastic epidemic models.

The basic reproduction number, ℛ(0), one of the most well-known thresholds in deterministic epidemic theory, predicts a disease outbreak if ℛ(0)>1. In stochastic epidemic theory, there are also thresholds that predict a major outbreak. In the case of a single infectious group, if ℛ(0)>1 and i infectious individuals are introduced into a susceptible population, then the probability of a major outbreak is approximately 1-(1/ℛ(0))( i ). With multiple infectious groups from which the disease could emerge, this result no longer holds. Stochastic thresholds for multiple groups depend on the number of individuals within each group, i ( j ), j=1, …, n, and on the probability of disease extinction for each group, q ( j ). It follows from multitype branching processes that the probability of a major outbreak is approximately [Formula: see text]. In this investigation, we summarize some of the deterministic and stochastic threshold theory, illustrate how to calculate the stochastic thresholds, and derive some new relationships between the deterministic and stochastic thresholds.

Stochastic models for competing species with a shared pathogen.

The presence of a pathogen among multiple competing species has important ecological implications. For example, a pathogen may change the competitive outcome, resulting in replacement of a native species by a non-native species. Alternately, if a pathogen becomes established, there may be a drastic reduction in species numbers. Stochastic variability in the birth, death and pathogen transmission processes plays an important role in determining the success of species or pathogen invasion. We investigate these phenomena while studying the dynamics of deterministic and stochastic models for n competing species with a shared pathogen. The deterministic model is a system of ordinary differential equations for n competing species in which a single shared pathogen is transmitted among the n species. There is no immunity from infection, individuals either die or recover and become immediately susceptible, an SIS disease model. Analytical results about pathogen persistence or extinction are summarized for the deterministic model for two and three species and new results about stability of the infection-free state and invasion by one species of a system of n-1 species are obtained. New stochastic models are derived in the form of continuous-time Markov chains and stochastic differential equations. Branching process theory is applied to the continuous-time Markov chain model to estimate probabilities for pathogen extinction or species invasion. Finally, numerical simulations are conducted to explore the effect of disease on two-species competition, to illustrate some of the analytical results and to highlight some of the differences in the stochastic and deterministic models.

Stochastic models for virus and immune system dynamics.

New stochastic models are developed for the dynamics of a viral infection and an immune response during the early stages of infection. The stochastic models are derived based on the dynamics of deterministic models. The simplest deterministic model is a well-known system of ordinary differential equations which consists of three populations: uninfected cells, actively infected cells, and virus particles. This basic model is extended to include some factors of the immune response related to Human Immunodeficiency Virus-1 (HIV-1) infection. For the deterministic models, the basic reproduction number, R0, is calculated and it is shown that if R0<1, the disease-free equilibrium is locally asymptotically stable and is globally asymptotically stable in some special cases. The new stochastic models are systems of stochastic differential equations (SDEs) and continuous-time Markov chain (CTMC) models that account for the variability in cellular reproduction and death, the infection process, the immune system activation, and viral reproduction. Two viral release strategies are considered: budding and bursting. The CTMC model is used to estimate the probability of virus extinction during the early stages of infection. Numerical simulations are carried out using parameter values applicable to HIV-1 dynamics. The stochastic models provide new insights, distinct from the basic deterministic models. For the case R0>1, the deterministic models predict the viral infection persists in the host. But for the stochastic models, there is a positive probability of viral extinction. It is shown that the probability of a successful invasion depends on the initial viral dose, whether the immune system is activated, and whether the release strategy is bursting or budding.

A habitat-based model for the spread of hantavirus between reservoir and spillover species.

New habitat-based models for spread of hantavirus are developed which account for interspecies interaction. Existing habitat-based models do not consider interspecies pathogen transmission, a primary route for emergence of new infectious diseases and reservoirs in wildlife and man. The modeling of interspecies transmission has the potential to provide more accurate predictions of disease persistence and emergence dynamics. The new models are motivated by our recent work on hantavirus in rodent communities in Paraguay. Our Paraguayan data illustrate the spatial and temporal overlaps among rodent species, one of which is the reservoir species for Jabora virus and others which are spillover species. Disease transmission occurs when their habitats overlap. Two mathematical models, a system of ordinary differential equations (ODE) and a continuous-time Markov chain (CTMC) model, are developed for spread of hantavirus between a reservoir and a spillover species. Analysis of a special case of the ODE model provides an explicit expression for the basic reproduction number, R(0), such that if R(0)<1, then the pathogen does not persist in either population but if R(0)>1, pathogen outbreaks or persistence may occur. Numerical simulations of the CTMC model display sporadic disease incidence, a new behavior of our habitat-based model, not present in other models, but which is a prominent feature of the seroprevalence data from Paraguay. Environmental changes that result in greater habitat overlap result in more encounters among various species that may lead to pathogen outbreaks and pathogen establishment in a new host.

Stochastic model of an influenza epidemic with drug resistance.

A continuous-time Markov chain (CTMC) model is formulated for an influenza epidemic with drug resistance. This stochastic model is based on an influenza epidemic model, expressed in terms of a system of ordinary differential equations (ODE), developed by Stilianakis, N.I., Perelson, A.S., Hayden, F.G., [1998. Emergence of drug resistance during an influenza epidemic: insights from a mathematical model. J. Inf. Dis. 177, 863-873]. Three different treatments-chemoprophylaxis, treatment after exposure but before symptoms, and treatment after symptoms appear, are considered. The basic reproduction number, R(0), is calculated for the deterministic-model under different treatment strategies. It is shown that chemoprophylaxis always reduces the basic reproduction number. In addition, numerical simulations illustrate that the basic reproduction number is generally reduced with realistic treatment rates. Comparisons are made among the different models and the different treatment strategies with respect to the number of infected individuals during an outbreak. The final size distribution is computed for the CTMC model and, in some cases, it is shown to have a bimodal distribution corresponding to two situations: when there is no outbreak and when an outbreak occurs. Given an outbreak occurs, the total number of cases for the CTMC model is in good agreement with the ODE model. The greatest number of drug resistant cases occurs if treatment is delayed or if only symptomatic individuals are treated.

A comparison of persistence-time estimation for discrete and continuous stochastic population models that include demographic and environmental variability.

A discrete-time Markov chain model, a continuous-time Markov chain model, and a stochastic differential equation model are compared for a population experiencing demographic and environmental variability. It is assumed that the environment produces random changes in the per capita birth and death rates, which are independent from the inherent random (demographic) variations in the number of births and deaths for any time interval. An existence and uniqueness result is proved for the stochastic differential equation system. Similarities between the models are demonstrated analytically and computational results are provided to show that estimated persistence times for the three stochastic models are generally in good agreement when the models satisfy certain consistency conditions.

A comparison of three different stochastic population models with regard to persistence time.

Results are summarized from the literature on three commonly used stochastic population models with regard to persistence time. In addition, several new results are introduced to clearly illustrate similarities between the models. Specifically, the relations between the mean persistence time and higher-order moments for discrete-time Markov chain models, continuous-time Markov chain models, and stochastic differential equation models are compared for populations experiencing demographic variability. Similarities between the models are demonstrated analytically, and computational results are provided to show that estimated persistence times for the three stochastic models are generally in good agreement when the models are consistently formulated. As an example, the three stochastic models are applied to a population satisfying logistic growth. Logistic growth is interesting as different birth and death rates can yield the same logistic differential equation. However, the persistence behavior of the population is strongly dependent on the explicit forms for the birth and death rates. Computational results demonstrate how dramatically the mean persistence time can vary for different populations that experience the same logistic growth.