Mathematical model for Zika virus dynamics with sexual transmission route

Zika virus is a flavivirus transmitted to humans primarily through the bite of infected Aedes mosquitoes. In addition to vector-borne spread, however, the virus can also be transmitted through sexual contact. In this paper, we formulate and analyze a new system of ordinary differential equations which incorporates both vector and sexual transmission routes. Theoretical analysis of this model when there is no disease induced mortality shows that the disease-free equilibrium is locally and globally asymptotically stable whenever the associated reproduction number is less than unity and unstable otherwise. However, when we extend this same model to include Zika induced mortality, which have been documented in Latin America, we find that the model exhibits a backward bifurcation. Specifically, a stable disease-free equilibrium co-exists with a stable endemic equilibrium when the associated reproduction number is less than unity. To further explore model predictions, we use numerical simulations to assess the importance of sexual transmission to disease dynamics. This analysis shows that risky behavior involving multiple sexual partners, particularly among male populations, substantially increases the number of infected individuals in the population, contributing significantly to the disease burden in the community.     

A model of bi-mode transmission dynamics of hepatitis C with optimal control

In this paper, we present a rigorous mathematical analysis of a deterministic model for the transmission dynamics of hepatitis C. The model is suitable for populations where two frequent modes of transmission of hepatitis C virus, namely unsafe blood transfusions and intravenous drug use, are dominant. The susceptible population is divided into two distinct compartments, the intravenous drug users and individuals undergoing unsafe blood transfusions. Individuals belonging to each compartment may develop acute and then possibly chronic infections. Chronically infected individuals may be quarantined. The analysis indicates that the eradication and persistence of the disease is completely determined by the magnitude of basic reproduction number R _c. It is shown that for the basic reproduction number R _c < 1, the disease-free equilibrium is locally and globally asymptotically stable. For R _c > 1, an endemic equilibrium exists and the disease is uniformly persistent. In addition, we present the uncertainty and sensitivity analyses to investigate the influence of different important model parameters on the disease prevalence. When the infected population persists, we have designed a time-dependent optimal quarantine strategy to minimize it. The Pontryagin’s Maximum Principle is used to characterize the optimal control in terms of an optimality system which is solved numerically. Numerical results for the optimal control are compared against the constant controls and their efficiency is discussed.     

The effects of human movement on the persistence of vector-borne diseases

With the recent resurgence of vector-borne diseases due to urbanization and development there is an urgent need to understand the dynamics of vector-borne diseases in rapidly changing urban environments. For example, many empirical studies have produced the disturbing finding that diseases continue to persist in modern city centers with zero or low rates of transmission. We develop spatial models of vector-borne disease dynamics on a network of patches to examine how the movement of humans in heterogeneous environments affects transmission. We show that the movement of humans between patches is sufficient to maintain disease persistence in patches with zero transmission. We construct two classes of models using different approaches: (i) Lagrangian models that mimic human commuting behavior and (ii) Eulerian models that mimic human migration. We determine the basic reproduction number R₀ for both modeling approaches. We show that for both approaches that if the disease-free equilibrium is stable (R₀<1) then it is globally stable and if the disease-free equilibrium is unstable (R₀>1) then there exists a unique positive (endemic) equilibrium that is globally stable among positive solutions. Finally, we prove in general that Lagrangian and Eulerian modeling approaches are not equivalent. The modeling approaches presented provide a framework to explore spatial vector-borne disease dynamics and control in heterogeneous environments. As an example, we consider two patches in which the disease dies out in both patches when there is no movement between them. Numerical simulations demonstrate that the disease becomes endemic in both patches when humans move between the two patches.     

Backward bifurcations in dengue transmission dynamics

A deterministic model for the transmission dynamics of a strain of dengue disease, which allows transmission by exposed humans and mosquitoes, is developed and rigorously analysed. The model, consisting of seven mutually-exclusive compartments representing the human and vector dynamics, has a locally-asymptotically stable disease-free equilibrium (DFE) whenever a certain epidemiological threshold, known as the basic reproduction number(R(0)) is less than unity. Further, the model exhibits the phenomenon of backward bifurcation, where the stable DFE coexists with a stable endemic equilibrium. The epidemiological consequence of this phenomenon is that the classical epidemiological requirement of making R(0) less than unity is no longer sufficient, although necessary, for effectively controlling the spread of dengue in a community. The model is extended to incorporate an imperfect vaccine against the strain of dengue. Using the theory of centre manifold, the extended model is also shown to undergo backward bifurcation. In both the original and the extended models, it is shown, using Lyapunov function theory and LaSalle Invariance Principle, that the backward bifurcation phenomenon can be removed by substituting the associated standard incidence function with a mass action incidence. In other words, in addition to establishing the presence of backward bifurcation in models of dengue transmission, this study shows that the use of standard incidence in modelling dengue disease causes the backward bifurcation phenomenon of dengue disease.     

Global dynamics of a staged-progression model with amelioration for infectious diseases

We analyze the global dynamics of a mathematical model for infectious diseases that progress through distinct stages within infected hosts with possibility of amelioration. An example of such diseases is HIV/AIDS that progresses through several stages with varying degrees of infectivity; amelioration can result from a host's immune action or more commonly from antiretroviral therapies, such as highly active antiretroviral therapy. For a general n-stage model with constant recruitment and bilinear incidence that incorporates amelioration, we prove that the global dynamics are completely determined by the basic reproduction number R(0). If R(0)≤1, then the disease-free equilibrium P(0) is globally asymptotically stable, and the disease always dies out. If R(0)>1, P(0) is unstable, a unique endemic equilibrium P* is globally asymptotically stable, and the disease persists at the endemic equilibrium. Impacts of amelioration on the basic reproduction number are also investigated.     

On the dynamics of SEIRS epidemic model with transport-related infection

Transportation amongst cities is found as one of the main factors which affect the outbreak of diseases. To understand the effect of transport-related infection on disease spread, an SEIRS (Susceptible, Exposed, Infectious, Recovered) epidemic model for two cities is formulated and analyzed. The epidemiological threshold, known as the basic reproduction number, of the model is derived. If the basic reproduction number is below unity, the disease-free equilibrium is locally asymptotically stable. Thus, the disease can be eradicated from the community. There exists an endemic equilibrium which is locally asymptotically stable if the reproduction number is larger than unity. This means that the disease will persist within the community. The results show that transportation among regions will change the disease dynamics and break infection out even if infectious diseases will go to extinction in each isolated region without transport-related infection. In addition, the result shows that transport-related infection intensifies the disease spread if infectious diseases break out to cause an endemic situation in each region, in the sense of that both the absolute and relative size of patients increase. Further, the formulated model is applied to the real data of SARS outbreak in 2003 to study the transmission of disease during the movement between two regions. The results show that the transport-related infection is effected to the number of infected individuals and the duration of outbreak in such the way that the disease becomes more endemic due to the movement between two cities. This study can be helpful in providing the information to public health authorities and policy maker to reduce spreading disease when its occurs.     

Mathematical Study of a Staged-Progression HIV Model with Imperfect Vaccine

A staged-progression HIV model is formulated and used to investigate the potential impact of an imperfect vaccine. The vaccine is assumed to have several desirable characteristics such as protecting against infection, causing bypass of the primary infection stage, and offering a disease-altering therapeutic effect (so that the vaccine induces reversal from the full blown AIDS stage to the asymptomatic stage). The model, which incorporates HIV transmission by individuals in the AIDS stage, is rigorously analyzed to gain insight into its qualitative features. Using a comparison theorem, the model with mass action incidence is shown to have a globally-asymptotically stable disease-free equilibrium whenever a certain threshold, known as the vaccination reproduction number, is less than unity. Furthermore, the model with mass action incidence has a unique endemic equilibrium whenever this threshold exceeds unity. Using the Li-Muldowney techniques for a reduced version of the mass action model, this endemic equilibrium is shown to be globally-asymptotically stable, under certain parameter restrictions. The epidemiological implications of these results are that an imperfect vaccine can eliminate HIV in a given community if it can reduce the reproduction number to a value less than unity, but the disease will persist otherwise. Furthermore, a future HIV vaccine that induces the bypass of primary infection amongst vaccinated individuals (who become infected) would decrease HIV prevalence, whereas a vaccine with therapeutic effect could have a positive or negative effect at the community level.     

Optimal control of a vectored plant disease model for a crop with continuous replanting

ABSTRACT

Vector-transmitted diseases of plants have had devastating effects on agricultural production worldwide, resulting in drastic reductions in yield for crops such as cotton, soybean, tomato, and cassava. Plant-vector-virus models with continuous replanting are investigated in terms of the effects of selection of cuttings, roguing, and insecticide use on disease prevalence in plants. Previous models are extended to include two replanting strategies: frequencyreplanting and abundance-replanting. In frequency-replanting, replanting of infected cuttings depends on the selection frequency parameter ε, whereas in abundance-replanting, replanting depends on plant abundance via a selection rate parameter also denoted as ε. The two models are analysed and new thresholds for disease elimination are defined for each model. Parameter values for cassava, whiteflies, and African cassava mosaic virus serve as a case study. A numerical sensitivity analysis illustrates how the equilibrium densities of healthy and infected plants vary with parameter values. Optimal control theory is used to investigate the effects of roguing and insecticide use with a goal of maximizing the healthy plants that are harvested. Differences in the control strategies in the two models are seen for large values of ε. Also, the combined strategy of roguing and insecticide use performs better than a single control.     

Dynamical analysis of a multi-strain model of HIV in the presence of anti-retroviral drugs

One major drawback associated with the use of anti-retroviral drugs in curtailing HIV spread in a population is the emergence and transmission of HIV strains that are resistant to these drugs. This paper presents a deterministic HIV treatment model, which incorporates a wild (drug sensitive) and a drug-resistant strain, for gaining insights into the dynamical features of the two strains, and determining effective ways to control HIV spread under this situation. Rigorous qualitative analysis of the model reveals that it has a globally asymptotically stable disease-free equilibrium whenever a certain epidemiological threshold (R t 0) is less than unity and that the disease will persist in the population when this threshold exceeds unity. Further, for the case where R t 0 > 1, it is shown that the model can have two co-existing endemic equilibria, and competitive exclusion phenomenon occurs whenever the associated reproduction number of the resistant strain (R t r) is greater than that of the wild strain (R t w). Unlike in the treatment model, it is shown that the model without treatment can have a family of infinitely many endemic equilibria when its associated epidemiological threshold (R(0)) exceeds unity. For the case when [Formula in text], it is shown that the widespread use of treatment against the wild strain can lead to its elimination from the community if the associated reduction in infectiousness of infected individuals (treated for the wild strain) does not exceed a certain threshold value (in this case, the use of treatment is expected to make R t w < R t r.     

Modelling interference between vectors of non-persistently transmitted plant viruses to identify effective control strategies

Aphids are the primary vector of plant viruses. Transient aphids, which probe several plants per day, are considered to be the principal vectors of non-persistently transmitted (NPT) viruses. However, resident aphids, which can complete their life cycle on a single host and are affected by agronomic practices, can transmit NPT viruses as well. Moreover, they can interfere both directly and indirectly with transient aphids, eventually shaping plant disease dynamics. By means of an epidemiological model, originally accounting for ecological principles and agronomic practices, we explore the consequences of fertilization and irrigation, pesticide deployment and roguing of infected plants on the spread of viral diseases in crops. Our results indicate that the spread of NPT viruses can be i) both reduced or increased by fertilization and irrigation, depending on whether the interference is direct or indirect; ii) counter-intuitively increased by pesticide application and iii) reduced by roguing infected plants. We show that a better understanding of vectors’ interactions would enhance our understanding of disease transmission, supporting the development of disease management strategies.     

Impact of chemo-therapy on optimal control of malaria disease with infected immigrants

We derived and analyzed rigorously a mathematical model that describes the dynamics of malaria infection with the recruitment of infected immigrants, treatment of infectives and spray of insecticides against mosquitoes in the population. Both qualitative and quantitative analysis of the deterministic model are performed with respect to stability of the disease free and endemic equilibria. It is found that in the absence of infected immigrants disease-free equilibrium is achievable and is locally asymptotically stable. Using Pontryagin's Maximum Principle, the optimal strategies for disease control are established. Finally, numerical simulations are performed to illustrate the analytical results.     

To Cut or Not to Cut: A Modeling Approach for Assessing the Role of Male Circumcision in HIV Control

A recent randomized controlled trial shows a significant reduction in women-to-men transmission of HIV due to male circumcision. Such development calls for a rigorous mathematical study to ascertain the full impact of male circumcision in reducing HIV burden, especially in resource-poor nations where access to anti-retroviral drugs is limited. First of all, this paper presents a compartmental model for the transmission dynamics of HIV in a community where male circumcision is practiced. In addition to having a disease-free equilibrium, which is locally-asymptotically stable whenever a certain epidemiological threshold is less than unity, the model exhibits the phenomenon of backward bifurcation, where the disease-free equilibrium coexists with a stable endemic equilibrium when the threshold is less than unity. The implication of this result is that HIV may persist in the population even when the reproduction threshold is less than unity. Using partial data from South Africa, the study shows that male circumcision at 60% efficacy level can prevent up to 220,000 cases and 8,200 deaths in the country within a year. Further, it is shown that male circumcision can significantly reduce, but not eliminate, HIV burden in a community. However, disease elimination is feasible if male circumcision is combined with other interventions such as ARVs and condom use. It is shown that the combined use of male circumcision and ARVs is more effective in reducing disease burden than the combined use of male circumcision and condoms for a moderate condom compliance rate.     

Global dynamics of an epidemiological model with age of infection and disease relapse

In this paper, an epidemiological model with age of infection and disease relapse is investigated. The basic reproduction number for the model is identified, and it is shown to be a sharp threshold to completely determine the global dynamics of the model. By analysing the corresponding characteristic equations, the local stability of a disease-free steady state and an endemic steady state of the model is established. By means of suitable Lyapunov functionals and LaSalle's invariance principle, it is verified that if the basic reproduction number is less than unity, the disease-free steady state is globally asymptotically stable, and hence the disease dies out; if the basic reproduction number is greater than unity, the endemic steady state is globally asymptotically stable and the disease becomes endemic.     

Dynamics of an epidemic model with non-local infections for diseases with latency over a patchy environment

In this paper, with the assumptions that an infectious disease in a population has a fixed latent period and the latent individuals of the population may disperse, we formulate an SIR model with a simple demographic structure for the population living in an n-patch environment (cities, towns, or countries, etc.). The model is given by a system of delay differential equations with a fixed delay accounting for the latency and a non-local term caused by the mobility of the individuals during the latent period. Assuming irreducibility of the travel matrices of the infection related classes, an expression for the basic reproduction number R₀{\mathcal{R}_0} is derived, and it is shown that the disease free equilibrium is globally asymptotically stable if R₀ < 1{\mathcal{R}_0 < 1} , and becomes unstable if ${\mathcal{R}_0 > 1}${\mathcal{R}_0 > 1} . In the latter case, there is at least one endemic equilibrium and the disease will be uniformly persistent. When n = 2, two special cases allowing reducible travel matrices are considered to illustrate joint impact of the disease latency and population mobility on the disease dynamics. In addition to the existence of the disease free equilibrium and interior endemic equilibrium, the existence of a boundary equilibrium and its stability are discussed for these two special cases.     

Role of incidence function in vaccine-induced backward bifurcation in some HIV models

The phenomenon of backward bifurcation in disease models, where a stable endemic equilibrium co-exists with a stable disease-free equilibrium when the associated reproduction number is less than unity, has important implications for disease control. In such a scenario, the classical requirement of the reproduction number being less than unity becomes only a necessary, but not sufficient, condition for disease elimination. This paper addresses the role of the choice of incidence function in a vaccine-induced backward bifurcation in HIV models. Several examples are given where backward bifurcations occur using standard incidence, but not with their equivalents that employ mass action incidence. Furthermore, this result is independent of the type of vaccination program adopted. These results emphasize the need for further work on the incidence functions used in HIV models.     

An SIS patch model with variable transmission coefficients

In this paper, an SIS patch model with non-constant transmission coefficients is formulated to investigate the effect of media coverage and human movement on the spread of infectious diseases among patches. The basic reproduction number R₀ is determined. It is shown that the disease-free equilibrium is globally asymptotically stable if R₀?1, and the disease is uniformly persistent and there exists at least one endemic equilibrium if R₀>1. In particular, when the disease is non-fatal and the travel rates of susceptible and infectious individuals in each patch are the same, the endemic equilibrium is unique and is globally asymptotically stable as R₀>1. Numerical calculations are performed to illustrate some results for the case with two patches.     

Impact of group mixing on disease dynamics

A general mathematical model is proposed to study the impact of group mixing in a heterogeneous host population on the spread of a disease that confers temporary immunity upon recovery. The model contains general distribution functions that account for the probabilities that individuals remain in the recovered class after recovery. For this model, the basic reproduction number R₀ is identified. It is shown that if R₀<1, then the disease dies out in the sense that the disease free equilibrium is globally asymptotically stable; whereas if R₀>1, this equilibrium becomes unstable. In this latter case, depending on the distribution functions and the group mixing strengths, the disease either persists at a constant endemic level or exhibits sustained oscillatory behavior.     

Rational behavioral response and the transmission of STDs

The susceptible-infected (SI) model is extended by allowing for individual optimal choices of self-protective actions against infection, where agents differ with respect to preferences and costs of self-protection. It is shown that a unique endemic equilibrium prevalence exists when the basic reproductive number of a STD is strictly greater than unity, and that the disease-free equilibrium is the unique steady state equilibrium when the basic reproductive number is less than or equal to one. Unlike in models that take individual behavior as given and fixed, the endemic equilibrium prevalence need not vary monotonically with respect to the basic reproductive number. Specifically, with endogenously determined self-protective behavior, a reduction in the basic reproductive number may in fact increase the endemic equilibrium prevalence. The global stability of the endemic steady state is established for the case of a homogeneous population by showing that, for any non-zero initial disease prevalence, there exists an equilibrium path which converges to the endemic steady state.     

Diseases with chronic stage in a population with varying size

An epidemiological model of hepatitis C with a chronic infectious stage and variable population size is introduced. A non-structured baseline ODE model which supports exponential solutions is discussed. The normalized version where the unknown functions are the proportions of the susceptible, infected, and chronic individuals in the total population is analyzed. It is shown that sustained oscillations are not possible and the endemic proportions either approach the disease-free or an endemic equilibrium. The expanded model incorporates the chronic age of the individuals. Partial analysis of this age-structured model is carried out. The global asymptotic stability of the infection-free state is established as well as local asymptotic stability of the endemic non-uniform steady state distribution under some additional conditions. A numerical method for the chronic-age-structured model is introduced. It is shown that this numerical scheme is consistent and convergent of first order. Simulations based on the numerical method suggest that in the structured case the endemic equilibrium may be unstable and sustained oscillations are possible. Closer look at the reproduction number reveals that treatment strategies directed towards speeding up the transition from acute to chronic stage in effect contribute to the eradication of the disease.     

Global stability of an SEIS epidemic model with recruitment and a varying total population size

This paper considers an SEIS epidemic model that incorporates constant recruitment, disease-caused death and disease latency. The incidence term is of the bilinear mass-action form. It is shown that the global dynamics is completely determined by the basic reproduction number R(0). If R(0)1, a unique endemic equilibrium is globally stable in the interior of the feasible region and the disease persists at the endemic equilibrium.