Mathematical Analysis of a Two Strain HIV/AIDS Model with Antiretroviral Treatment

A two strain HIV/AIDS model with treatment which allows AIDS patients with sensitive HIV-strain to undergo amelioration is presented as a system of non-linear ordinary differential equations. The disease-free equilibrium is shown to be globally asymptotically stable when the associated epidemic threshold known as the basic reproduction number for the model is less than unity. The centre manifold theory is used to show that the sensitive HIV-strain only and resistant HIV-strain only endemic equilibria are locally asymptotically stable when the associated reproduction numbers are greater than unity. Qualitative analysis of the model including positivity, boundedness and persistence of solutions are presented. The model is numerically analysed to assess the effects of treatment with amelioration on the dynamics of a two strain HIV/AIDS model. Numerical simulations of the model show that the two strains co-exist whenever the reproduction numbers exceed unity. Further, treatment with amelioration may result in an increase in the total number of infective individuals (asymptomatic) but results in a decrease in the number of AIDS patients. Further, analysis of the reproduction numbers show that antiretroviral resistance increases with increase in antiretroviral use.     

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.     

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.     

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.     

Analysis of Risk-Structured Vaccination Model for the Dynamics of Oncogenic and Warts-Causing HPV Types

A new deterministic model is designed and used to assess the community-wide impact of mass vaccination of new sexually active individuals on the dynamics of the oncogenic and warts-causing HPV types. Rigorous qualitative analyses of the model, which incorporates the two currently available anti-HPV vaccines, reveal that it undergoes competitive exclusion when the reproduction of one HPV risk type (low/high) exceeds unity, while that of the other HPV risk type is less than unity. For the case when the reproduction numbers of the two HPV risk types (low/high) exceed unity, the two risk types co-exist. It is shown that the sub-model with the low-risk HPV types only has at least one endemic equilibrium whenever the associated reproduction threshold exceeds unity. Furthermore, this sub-model undergoes a re-infection-induced backward bifurcation under certain conditions. In the absence of the re-infection of recovered individuals and cancer-induced mortality in males, the associated disease-free equilibrium of the full (risk-structured) model is shown to be globally asymptotically stable whenever the reproduction number of the model is less than unity (that is, the full model does not undergo backward bifurcation under this setting). It is shown, via numerical simulations, that the use of the Gardasil vaccine could lead to the effective control of HPV in the community if the coverage rate is in the range of 73–95 % (84 %). If 70 % of the new sexually active susceptible females are vaccinated with the Gardasil vaccine, additionally vaccinating 34–56 % (45 %) of the new sexually active susceptible males can lead to the effective community-wide control (or elimination) of the HPV types.     

Endemic threshold results in an age-duration-structured population model for HIV infection

In this paper we consider an age-duration-structured population model for HIV infection in a homosexual community. First we investigate the invasion problem to establish the basic reproduction ratio R(0) for the HIV/AIDS epidemic by which we can state the threshold criteria: The disease can invade into the completely susceptible population if R(0)>1, whereas it cannot if R(0)<1. Subsequently, we examine existence and uniqueness of endemic steady states. We will show sufficient conditions for a backward or a forward bifurcation to occur when the basic reproduction ratio crosses unity. That is, in contrast with classical epidemic models, for our HIV model there could exist multiple endemic steady states even if R(0) is less than one. Finally, we show sufficient conditions for the local stability of the endemic steady states.     

Dynamics Analysis of a Multi-strain Cholera Model with an Imperfect Vaccine

A new two-strain model, for assessing the impact of basic control measures, treatment and dose-structured mass vaccination on cholera transmission dynamics in a population, is designed. The model has a globally-asymptotically stable disease-free equilibrium whenever its associated reproduction number is less than unity. The model has a unique, and locally-asymptotically stable, endemic equilibrium when the threshold quantity exceeds unity and another condition holds. Numerical simulations of the model show that, with the expected 50 % minimum efficacy of the first vaccine dose, vaccinating 55 % of the susceptible population with the first vaccine dose will be sufficient to effectively control the spread of cholera in the community. Such effective control can also be achieved if 50 % of the first vaccine dose recipients take the second dose. It is shown that a control strategy that emphasizes the use of antibiotic treatment is more effective than one that emphasizes the use of basic (non-pharmaceutical) anti-cholera control measures only. Numerical simulations show that, while the universal strategy (involving all three control measures) gives the best outcome in minimizing cholera burden in the community, the combined basic anti-cholera control measures and treatment strategy also has very effective community-wide impact.     

Analysis of SIR epidemic models with nonlinear incidence rate and treatment

This paper deals with the nonlinear dynamics of a susceptible-infectious-recovered (SIR) epidemic model with nonlinear incidence rate, vertical transmission, vaccination for the newborns of susceptible and recovered individuals, and the capacity of treatment. It is assumed that the treatment rate is proportional to the number of infectives when it is below the capacity and constant when the number of infectives reaches the capacity. Under some conditions, it is shown that there exists a backward bifurcation from an endemic equilibrium, which implies that the disease-free equilibrium coexists with an endemic equilibrium. In such a case, reducing the basic reproduction number less than unity is not enough to control and eradicate the disease, extra measures are needed to ensure that the solutions approach the disease-free equilibrium. When the basic reproduction number is greater than unity, the model can have multiple endemic equilibria due to the effect of treatment, vaccination and other parameters. The existence and stability of the endemic equilibria of the model are analyzed and sufficient conditions on the existence and stability of a limit cycle are obtained. Numerical simulations are presented to illustrate the analytical results.     

HIV/AIDS dynamics: Impact of economic classes with transmission from poor clinical settings

We formulate and analyze a nonlinear deterministic HIV/AIDS model with two social classes, namely the poor and the rich including transmission from poor clinical settings with a randomly variable population. Four sub-models are derived from the full model, the disease threshold parameters are computed, and it is shown that the disease will die down if these initial threshold parameters are less than unity and will persist if they exceed unity. The impact of economic classes (along with transmission from poor/inadequate clinical settings) on the disease dynamics is assessed, and we observe that even with a single sexual partner, the reproduction number is slightly greater than unity, implying that the additional transmission can only be from clinical settings. Stability (local and global) of both the disease-free and endemic equilibria are then investigated using various techniques of dynamical systems such as Centre Manifold theory and Lyapunov's second method. Analysis on the bifurcation parameter is carried out to assess the impact of related HIV transmission from poor clinical settings. We estimate some of the model parameter values and numerical simulations of the model are represented graphically. Our results show that the prevalence of HIV in rich communities seems to be higher than that in the poor, but the disease develops faster in impoverished individuals.     

Modelling HIV/AIDS in the presence of an HIV testing and screening campaign

Preventing and managing the HIV/AIDS epidemic in South Africa will dominate the next decade and beyond. Reduction of new HIV infections by implementing a comprehensive national HIV prevention programme at a sufficient scale to have real impact remains a priority. In this paper, a deterministic HIV/AIDS model that incorporates condom use, screening through HIV counseling and testing (HCT), regular testing and treatment as control strategies is proposed with the objective of quantifying the effectiveness of HCT in preventing new infections and predicting the long-term dynamics of the epidemic. It is found that a backward bifurcation occurs if the rate of screening is below a certain threshold, suggesting that the classical requirement for the basic reproduction number to be below unity though necessary, is not sufficient for disease control in this case. The global stabilities of the equilibria under certain conditions are determined in terms of the model reproduction number R₀. Numerical simulations are performed and the model is fitted to data on HIV prevalence in South Africa. The effects of changes in some key epidemiological parameters are investigated. Projections are made to predict the long-term dynamics of the disease. The epidemiological implications of such projections on public health planning and management are discussed.     

Malaria Drug Resistance: The Impact of Human Movement and Spatial Heterogeneity

Human habitat connectivity, movement rates, and spatial heterogeneity have tremendous impact on malaria transmission. In this paper, a deterministic system of differential equations for malaria transmission incorporating human movements and the development of drug resistance malaria in an \(n\) patch system is presented. The disease-free equilibrium of the model is globally asymptotically stable when the associated reproduction number is less than unity. For a two patch case, the boundary equilibria (drug sensitive-only and drug resistance-only boundary equilibria) when there is no movement between the patches are shown to be locally asymptotically stable when they exist; the co-existence equilibrium is locally asymptotically stable whenever the reproduction number for the drug sensitive malaria is greater than the reproduction number for the resistance malaria. Furthermore, numerical simulations of the connected two patch model (when there is movement between the patches) suggest that co-existence or competitive exclusion of the two strains can occur when the respective reproduction numbers of the two strains exceed unity. With slow movement (or low migration) between the patches, the drug sensitive strain dominates the drug resistance strain. However, with fast movement (or high migration) between the patches, the drug resistance strain dominates the drug sensitive strain.     

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.     

Modelling the Transmission Dynamics and Control of the Novel 2009 Swine Influenza (H1N1) Pandemic

The paper presents a deterministic compartmental model for the transmission dynamics of swine influenza (H1N1) pandemic in a population in the presence of an imperfect vaccine and use of drug therapy for confirmed cases. Rigorous analysis of the model, which stratifies the infected population in terms of their risk of developing severe illness, reveals that it exhibits a vaccine-induced backward bifurcation when the associated reproduction number is less than unity. The epidemiological consequence of this result is that the effective control of H1N1, when the reproduction number is less than unity, in the population would then be dependent on the initial sizes of the subpopulations of the model. For the case where the vaccine is perfect, it is shown that having the reproduction number less than unity is necessary and sufficient for effective control of H1N1 in the population (in such a case, the associated disease-free equilibrium is globally asymptotically stable). The model has a unique endemic equilibrium when the reproduction number exceeds unity. Numerical simulations of the model, using data relevant to the province of Manitoba, Canada, show that it reasonably mimics the observed H1N1 pandemic data for Manitoba during the first (Spring) wave of the pandemic. Further, it is shown that the timely implementation of a mass vaccination program together with the size of the Manitoban population that have preexisting infection-acquired immunity (from the first wave) are crucial to the magnitude of the expected burden of disease associated with the second wave of the H1N1 pandemic. With an estimated vaccine efficacy of approximately 80%, it is projected that at least 60% of Manitobans need to be vaccinated in order for the effective control or elimination of the H1N1 pandemic in the province to be feasible. Finally, it is shown that the burden of the second wave of H1N1 is expected to be at least three times that of the first wave, and that the second wave would last until the end of January or early February, 2010.     

Mathematical Analysis of the Transmission Dynamics of HIV Syphilis Co-infection in the Presence of Treatment for Syphilis

The re-emergence of syphilis has become a global public health issue, and more persons are getting infected, especially in developing countries. This has also led to an increase in the incidence of human immunodeficiency virus (HIV) infections as some studies have shown in the recent decade. This paper investigates the synergistic interaction between HIV and syphilis using a mathematical model that assesses the impact of syphilis treatment on the dynamics of syphilis and HIV co-infection in a human population where HIV treatment is not readily available or accessible to HIV-infected individuals. In the absence of HIV, the syphilis-only model undergoes the phenomenon of backward bifurcation when the associated reproduction number (\({\mathcal {R}}_{T}\)) is less than unity, due to susceptibility to syphilis reinfection after recovery from a previous infection. The complete syphilis–HIV co-infection model also undergoes the phenomenon of backward bifurcation when the associated effective reproduction number (\({\mathcal {R}}_{C}\)) is less than unity for the same reason as the syphilis-only model. When susceptibility to syphilis reinfection after treatment is insignificant, the disease-free equilibrium of the syphilis-only model is shown to be globally asymptotically stable whenever the associated reproduction number (\({\mathcal {R}}_{T}\)) is less than unity. Sensitivity and uncertainty analysis show that the top three parameters that drive the syphilis infection (with respect to the associated response function, \({\mathcal {R}}_{T}\)) are the contact rate (\(\beta _S\)), modification parameter that accounts for the increased infectiousness of syphilis-infected individuals in the secondary stage of the infection (\(\theta _1\)) and treatment rate for syphilis-only infected individuals in the primary stage of the infection (\(r_1\)). The co-infection model was numerically simulated to investigate the impact of various treatment strategies for primary and secondary syphilis, in both singly and dually infected individuals, on the dynamics of the co-infection of syphilis and HIV. It is observed that if concerted effort is exerted in the treatment of primary and secondary syphilis (in both singly and dually infected individuals), especially with high treatment rates for primary syphilis, this will result in a reduction in the incidence of HIV (and its co-infection with syphilis) in the population.     

Threshold dynamics of a non-autonomous SEIRS model with quarantine and isolation

A model for assessing the effect of periodic fluctuations on the transmission dynamics of a communicable disease, subject to quarantine (of asymptomatic cases) and isolation (of individuals with clinical symptoms of the disease), is considered. The model, which is of a form of a non-autonomous system of non-linear differential equations, is analysed qualitatively and numerically. It is shown that the disease-free solution is globally-asymptotically stable whenever the associated basic reproduction ratio of the model is less than unity, and the disease persists in the population when the reproduction ratio exceeds unity. This study shows that adding periodicity to the autonomous quarantine/isolation model developed in Safi and Gumel (Discret Contin Dyn Syst Ser B 14:209–231, 2010) does not alter the threshold dynamics of the autonomous system with respect to the elimination or persistence of the disease in the population.     

Mathematical Study of the Role of Gametocytes and an Imperfect Vaccine on Malaria Transmission Dynamics

A mathematical model is developed to assess the role of gametocytes (the infectious sexual stage of the malaria parasite) in malaria transmission dynamics in a community. The model is rigorously analysed to gain insights into its dynamical features. It is shown that, in the absence of disease-induced mortality, the model has a globally-asymptotically stable disease-free equilibrium whenever a certain epidemiological threshold, known as the basic reproduction number (denoted by ℛ₀), is less than unity. Further, it has a unique endemic equilibrium if ℛ₀>1. The model is extended to incorporate an imperfect vaccine with some assumed therapeutic characteristics. Theoretical analyses of the model with vaccination show that an imperfect malaria vaccine could have negative or positive impact (in reducing disease burden) depending on whether or not a certain threshold (denoted by ∇) is less than unity. Numerical simulations of the vaccination model show that such an imperfect anti-malaria vaccine (with a modest efficacy and coverage rate) can lead to effective disease control if the reproduction threshold (denoted by ℛ_vac) of the disease is reasonably small. On the other hand, the disease cannot be effectively controlled using such a vaccine if ℛ_vac is high. Finally, it is shown that the average number of days spent in the class of infectious individuals with higher level of gametocyte is critically important to the malaria burden in the community.     

Is HIV infection associated with an increased risk for cholera? Insights from a mathematical model

Cholera, a waterborne gastroenteric infection, remains a significant threat to public health in sub-Saharan Africa, the region most heavily affected by HIV. It is biologically plausible that immune suppression caused by HIV infection predisposes to cholera. In this paper, a simple mathematical model is developed and comprehensively analyzed to assess whether HIV infection is associated with an increased risk for cholera or not. Analytical results of the model show that the quantities R(c) and R(h), which represents the reproductive number for cholera and HIV infection, respectively, provide threshold conditions that determine cumulative new single and dual infection cases. These threshold conditions can be used to gain important insights on the epidemiological consequences of HIV and cholera coexistence in the community. Numerical results are provided to support the analytical findings. The findings suggest that in a cholera-endemic area, HIV infection is associated with an increased risk for cholera.     

Linking immunological and epidemiological dynamics of HIV: the case of super-infection

In this paper, a two-strain model that links immunological and epidemiological dynamics across scales is formulated. On the within-host scale, the two strains eliminate each other with the strain with the larger immunological reproduction persisting. However, on the population scale superinfection is possible, with the strain with larger immunological reproduction number super-infecting the strain with the smaller immunological reproduction number. The two models are linked through the age-since-infection structure of the epidemiological variables. In addition, the between-host transmission and the disease-induced death rate depend on the within-host viral load. The immunological reproduction numbers, the epidemiological reproduction numbers and invasion reproduction numbers are computed. Besides the disease-free equilibrium, there are two population-level strain one and strain two isolated equilibria, as well as a population-level coexistence equilibrium when both invasion reproduction numbers are greater than one. The single-strain population-level equilibria are locally asymptotically stable suggesting that in the absence of superinfection oscillations do not occur, a result contrasting previous studies of HIV age-since-infection structured models. Simulations suggest that the epidemiological reproduction number and HIV population prevalence are monotone functions of the within-host parameters with reciprocal trends. In particular, HIV medications that decrease within-host viral load also increase overall population prevalence. The effect of the immunological parameters on the population reproduction number and prevalence is more pronounced when the initial viral load is lower.     

Theoretical Assessment of Public Health Impact of Imperfect Prophylactic HIV-1 Vaccines with Therapeutic Benefits

This paper presents a number of deterministic models for theoretically assessing the potential impact of an imperfect prophylactic HIV-1 vaccine that has five biological modes of action, namely “take,” “degree,” “duration,” “infectiousness,” and “progression,” and can lead to increased risky behavior. The models, which are of the form of systems of nonlinear differential equations, are constructed via a progressive refinement of a basic model to incorporate more realistic features of HIV pathogenesis and epidemiology such as staged progression, differential infectivity, and HIV transmission by AIDS patients. The models are analyzed to gain insights into the qualitative features of the associated equilibria. This allows the determination of important epidemiological thresholds such as the basic reproduction numbers and a measure for vaccine impact or efficacy. The key findings of the study include the following (i) if the vaccinated reproduction number is greater than unity, each of the models considered has a locally unstable disease-free equilibrium and a unique endemic equilibrium; (ii) owing to the vaccine-induced backward bifurcation in these models, the classical epidemiological requirement of vaccinated reproduction number being less than unity does not guarantee disease elimination in these models; (iii) an imperfect vaccine will reduce HIV prevalence and mortality if the reproduction number for a wholly vaccinated population is less than the corresponding reproduction number in the absence of vaccination; (iv) the expressions for the vaccine characteristics of the refined models take the same general structure as those of the basic model.