Mathematical Assessment of the Role of Early Latent Infections and Targeted Control Strategies on Syphilis Transmission Dynamics.

A new multi-stage deterministic model for the transmission dynamics of syphilis, which incorporates disease transmission by individuals in the early latent stage of syphilis infection and the reversions of early latent syphilis to the primary and secondary stages, is formulated and rigorously analysed. The model is used to assess the population-level impact of preventive (condom use) and therapeutic measures (treatment using antibiotics) against the spread of the disease in a community. It is shown that the disease-free equilibrium of the model is globally-asymptotically stable whenever the associated control reproduction number (denoted by [Formula: see text]) is less than unity. A special case of the model is shown to have a unique and globally-asymptotically stable endemic equilibrium whenever the associated reproduction number (denoted by [Formula: see text]) exceeds unity. Uncertainty and sensitivity analysis of the model, using parameter values and ranges relevant to syphilis transmission dynamics in Nigeria, show that the top three parameters that drive the syphilis infection (with respect to [Formula: see text]) are the disease transmission rate ([Formula: see text]), compliance in condom use (c) and efficacy of condom ([Formula: see text]). Numerical simulations of the model show that the targeted treatment of secondary syphilis cases is more effective than the targeted treatment of individuals in the primary or early latent stage of syphilis infection.

Comments on "A Mathematical Study to Control Visceral Leishmaniasis: An Application to South Sudan".

Deterministic (ordinary differential equation) models for the transmission dynamics of vector-borne diseases that incorporate disease-induced death in the host(s) population(s) are generally known to exhibit the phenomenon of backward bifurcation (where a stable disease-free equilibrium of the model coexists with a stable endemic equilibrium when the associated reproduction number of the model is less than unity). Further, it is well known that, in these models, the phenomenon of backward bifurcation does not occur when the disease-induced death rate is negligible (e.g., if the disease-induced death rate is set to zero). In a recent paper on the transmission dynamics of visceral leishmaniasis (a disease vectored by sandflies), titled "A Mathematical Study to Control Visceral Leishmaniasis: An Application to South Sudan," published in Bulletin of Mathematical Biology, Vol. 79, Pages 1110-1134, 2017, Ghosh et al. (2017) stated that their deterministic model undergoes a backward bifurcation even when the disease-induced mortality in the host population is set to zero. This result is contrary to the well-established theory on the dynamics of vector-borne diseases. In this short note, we illustrate some of the key errors in the Ghosh et al. (2017) study.

Mathematical analysis of a model for AVL-HIV co-endemicity.

A model for the transmission dynamics of Anthroponotic Visceral Leishmaniasis (AVL) and human immunodeficiency virus (HIV) in a population is developed and used to assess the impact of the spread of each disease on the overall transmission dynamics. As for other vector-borne disease models, the AVL component of the model undergoes backward bifurcation when the associated reproduction number of the AVL-only sub-model (denoted by RL) is less than unity. Uncertainty and sensitivity analyzes of the model, using data relevant to the dynamics of the two diseases in Ethiopia, show that the top three parameters that drive the AVL infection (with respect to the associated response function, RL) are the average number of times a sandfly bites humans per unit time (σV), carrying capacity of vectors (KV) and transmission probability from infected humans to susceptible sandflies (ß2). The distribution of RL is RL∈[0.06,3.94] with a mean of RL=1.08. Furthermore, the top three parameters that affect HIV dynamics (with respect to the response function RH) are the transmission rate of HIV (ßH), HIV-induced death rate (δH), and the modification parameter for the increase in infectiousness of AIDS individuals in comparison to HIV infected without clinical symptoms of AIDS (ωH). The distribution of RH is RH∈[0.88,2.79] with a mean of RH=1.46. The dominant parameters that affect the dynamics of the full VL-HIV model (with respect to the associated reproduction number, RLH, as the response function) are the transmission rate of HIV (ßH), the average number of times a sandfly bites humans per unit time (σV), and HIV-induced death rate (δH) (the distribution of RLH is RLH∈[0.88,3.94] with a mean of RLH=1.64). Numerical simulations of the model show that the two diseases co-exist (with AVL dominating, but not driving HIV to extinction) whenever the reproduction number of each disease exceeds unity. It is shown that AVL can invade a population at HIV-endemic state if a certain threshold quantity, known as invasion reproduction number, exceeds unity.

Differential characteristics of primary infection and re-infection can cause backward bifurcation in HCV transmission dynamics.

Backward bifurcation, a phenomenon typically characterized by the co-existence of multiple stable equilibria when the associated reproduction number of the model is less than unity, has been observed in numerous disease transmission models. This study establishes, for the first time, the presence of this phenomenon in the transmission dynamics of hepatitis C virus (HCV) within an IDU population. It is shown that the phenomenon does not exist under four scenarios, namely (i) in the absence of re-infection, (ii) in the absence of differential characteristics of HCV infection (with respect to infectivity, progression, treatment and recovery) between re-infected individuals and primary-infected individuals, (iii) when re-infected and treated individuals do not transmit HCV infection and (iv) when the average infectivity-adjusted duration of re-infection is less than that of primary infection. This study identifies, using sensitivity analysis, five parameters of the model that have the most influence on the disease transmission dynamics, namely: effective contact rate, progression rate from acute to chronic infection, recovery rate from acute infection, natural death rate and the relative infectiousness of chronically-infected individuals. Numerical simulations of the model show that the re-infection of recovered individuals has marginal effect on the HCV burden (as measured in terms of the cumulative incidence and the prevalence of the disease) in the IDU community. Furthermore, treatment of infected IDUs, even for small rate (such as 4%), offers significant impact on curtailing HCV spread in the community.

Dynamics of Mycobacterium and bovine tuberculosis in a human-buffalo population.

A new model for the transmission dynamics of Mycobacterium tuberculosis and bovine tuberculosis in a community, consisting of humans and African buffalos, is presented. The buffalo-only component of the model exhibits the phenomenon of backward bifurcation, which arises due to the reinfection of exposed and recovered buffalos, when the associated reproduction number is less than unity. This model has a unique endemic equilibrium, which is globally asymptotically stable for a special case, when the reproduction number exceeds unity. Uncertainty and sensitivity analyses, using data relevant to the dynamics of the two diseases in the Kruger National Park, show that the distribution of the associated reproduction number is less than unity (hence, the diseases would not persist in the community). Crucial parameters that influence the dynamics of the two diseases are also identified. Both the buffalo-only and the buffalo-human model exhibit the same qualitative dynamics with respect to the local and global asymptotic stability of their respective disease-free equilibrium, as well as with respect to the backward bifurcation phenomenon. Numerical simulations of the buffalo-human model show that the cumulative number of Mycobacterium tuberculosis cases in humans (buffalos) decreases with increasing number of bovine tuberculosis infections in humans (buffalo).

Qualitative dynamics of lowly- and highly-pathogenic avian influenza strains.

A new deterministic model for the transmission dynamics of the lowly- and highly-pathogenic avian influenza (LPAI and HPAI) strains is designed and rigorously analyzed. The model exhibits the phenomenon of backward bifurcation, where a stable disease-free equilibrium co-exists with a stable endemic equilibrium whenever the associated reproduction number is less than unity. It is shown that the re-infection of birds infected with the LPAI strain causes the backward bifurcation phenomenon. In the absence of such re-infection, the disease-free equilibrium of the model is globally-asymptotically stable when the associated reproduction number is less than unity. Using non-linear Lyapunov functions of Goh-Volterra type, the LPAI-only and HPAI-only boundary equilibria of the model are shown to be globally-asymptotically stable when they exist. A special case of the model is shown to have a continuum of co-existence equilibria whenever the associated reproduction numbers of the two strains are equal and exceed unity. Furthermore, numerical simulations of the model suggest that co-existence or competitive exclusion of the two strains can occur when the respective reproduction numbers of the two strains exceed unity.

Mathematical Analysis of a Model for Assessing the Impact of Antiretroviral Therapy, Voluntary Testing and Condom Use in Curtailing the Spread of HIV.

This paper presents a deterministic model for evaluating the impact of anti-retroviral drugs (ARVs), voluntary testing (using standard antibody-based and a DNA-based testing methods) and condom use on the transmission dynamics of HIV in a community. Rigorous qualitative analysis of the model show that it has a globally-stable disease-free equilibrium whenever a certain epidemiological threshold, known as the effective reproduction number , is less than unity. The model has an endemic equilibrium whenever . The endemic equilibrium is shown to be locally-asymptotically stable for a special case. Numerical simulations of the model show that the use of the combined testing and treatment strategy is more effective than the use of the standard ELISA testing method with ARV treatment, even for the use of condoms as a singular strategy. Furthermore, the universal strategy (which involves the use of condoms, the two testing methods and ARV treatment) is always more effective than the combined use of the standard ELISA testing method and ARVs.

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.

Qualitative assessment of the role of public health education program on HIV transmission dynamics.

This paper presents a non-linear deterministic model for assessing the impact of public health education campaign on curtailing the spread of the human immunodeficiency virus (HIV) pandemic in a population. Rigorous qualitative analysis of the model reveals that it exhibits the phenomenon of backward bifurcation (BB), where a stable disease-free equilibrium coexists with a stable endemic equilibrium when a certain threshold quantity, known as the 'effective reproduction number' ('Reff), is less than unity. The epidemiological implication of BB is that a public health education campaign could fail to effectively control HIV even when the classical requirement of having the associated reproduction number less than unity is satisfied. Furthermore, an explicit threshold value is derived above which such an education campaign could lead to detrimental outcome (increase disease burden) and below which it would have positive population-level impact (reduce disease burden in the community). It is shown that the BB phenomenon is caused by imperfect efficacy of the public health education program. The model is used to assess the potential impact of some targeted public health education campaigns using data from numerous countries.

Protecting residential care facilities from pandemic influenza.

It is widely believed that protecting health care facilities against outbreaks of pandemic influenza requires pharmaceutical resources such as antivirals and vaccines. However, early in a pandemic, vaccines will not likely be available and antivirals will probably be of limited supply. The containment of pandemic influenza within acute-care hospitals anywhere is problematic because of open connections with communities. However, other health care institutions, especially those providing care for the disabled, can potentially control community access. We modeled a residential care facility by using a stochastic compartmental model to address the question of whether conditions exist under which nonpharmaceutical interventions (NPIs) alone might prevent the introduction of a pandemic virus. The model projected that with currently recommended staff-visitor interactions and social distancing practices, virus introductions are inevitable in all pandemics, accompanied by rapid internal propagation. The model identified staff reentry as the critical pathway of contagion, and provided estimates of the reduction in risk required to minimize the probability of a virus introduction. By using information on latency for historical and candidate pandemic viruses, we developed NPIs that simulated notions of protective isolation for staff away from the facility that reduced the probability of bringing the pandemic infection back to the facility to levels providing protection over a large range of projected pandemic severities. The proposed form of protective isolation was evaluated for social plausibility by collaborators who operate residential facilities. It appears unavoidable that NPI combinations effective against pandemics more severe than mild imply social disruption that increases with severity.

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.

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.

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.

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.

Assessing the role of basic control measures, antivirals and vaccine in curtailing pandemic influenza: scenarios for the US, UK and the Netherlands.

An increasing number of avian flu cases in humans, arising primarily from direct contact with poultry, in several regions of the world have prompted the urgency to develop pandemic preparedness plans worldwide. Leading recommendations in these plans include basic public health control measures for minimizing transmission in hospitals and communities, the use of antiviral drugs and vaccination. This paper presents a mathematical model for the evaluation of the pandemic flu preparedness plans of the United States (US), the United Kingdom (UK) and the Netherlands. The model is used to assess single and combined interventions. Using data from the US, we show that hospital and community transmission control measures alone can be highly effective in reducing the impact of a potential flu pandemic. We further show that while the use of antivirals alone could lead to very significant reductions in the burden of a pandemic, the combination of transmission control measures, antivirals and vaccine gives the most 'optimal' result. However, implementing such an optimal strategy at the onset of a pandemic may not be realistic. Thus, it is important to consider other plausible alternatives. An optimal preparedness plan is largely dependent on the availability of resources; hence, it is country-specific. We show that countries with limited antiviral stockpiles should emphasize their use therapeutically (rather than prophylactically). However, countries with large antiviral stockpiles can achieve greater reductions in disease burden by implementing them both prophylactically and therapeutically. This study promotes alternative strategies that may be feasible and attainable for the US, UK and the Netherlands. It emphasizes the role of hospital and community transmission control measures in addition to the timely administration of antiviral treatment in reducing the burden of a flu pandemic. The latter is consistent with the preparedness plans of the UK and the Netherlands. Our results indicate that for low efficacy and coverage levels of antivirals and vaccine, the use of a vaccine leads to the greatest reduction in morbidity and mortality compared with the singular use of antivirals. However, as these efficacy and coverage levels are increased, the use of antivirals is more effective.

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.

An sveir model for assessing potential impact of an imperfect anti-sars vaccine.

The control of severe acute respiratory syndrome (SARS), a fatal contagious viral disease that spread to over 32 countries in 2003, was based on quarantine of latently infected individuals and isolation of individuals with clinical symptoms of SARS. Owing to the recent ongoing clinical trials of some candidate anti-SARS vaccines, this study aims to assess, via mathematical modelling, the potential impact of a SARS vaccine, assumed to be imperfect, in curtailing future outbreaks. A relatively simple deterministic model is designed for this purpose. It is shown, using Lyapunov function theory and the theory of compound matrices, that the dynamics of the model are determined by a certain threshold quantity known as the control reproduction number (R(v)). If R(v) =/< 1, the disease will be eliminated from the community; whereas an epidemic occurs if R(v) > 1. This study further shows that an imperfect SARS vaccine with infection-blocking efficacy is always beneficial in reducing disease spread within the community, although its overall impact increases with increasing efficacy and coverage. In particular, it is shown that the fraction of individuals vaccinated at steady-state and vaccine efficacy play equal roles in reducing disease burden, and the vaccine must have efficacy of at least 75% to lead to effective control of SARS (assuming R(0) = 4). Numerical simulations are used to explore the severity of outbreaks when R(v) > 1.

A mathematical model for assessing control strategies against West Nile virus.

Since its incursion into North America in 1999, West Nile virus (WNV) has spread rapidly across the continent resulting in numerous human infections and deaths. Owing to the absence of an effective diagnostic test and therapeutic treatment against WNV, public health officials have focussed on the use of preventive measures in an attempt to halt the spread of WNV in humans. The aim of this paper is to use mathematical modelling and analysis to assess two main anti-WNV preventive strategies, namely: mosquito reduction strategies and personal protection. We propose a single-season ordinary differential equation model for the transmission dynamics of WNV in a mosquito-bird-human community, with birds as reservoir hosts and culicine mosquitoes as vectors. The model exhibits two equilibria; namely the disease-free equilibrium and a unique endemic equilibrium. Stability analysis of the model shows that the disease-free equilibrium is globally asymptotically stable if a certain threshold quantity (R0), which depends solely on parameters associated with the mosquito-bird cycle, is less than unity. The public health implication of this is that WNV can be eradicated from the mosquito-bird cycle (and, consequently, from the human population) if the adopted mosquito reduction strategy (or strategies) can make R0<1. On the other hand, it is shown, using a novel and robust technique that is based on the theory of monotone dynamical systems coupled with a regular perturbation argument and a Liapunov function, that if R0>1, then the unique endemic equilibrium is globally stable for small WNV-induced avian mortality. Thus, in this case, WNV persists in the mosquito-bird population.

Numerical modelling of the perturbation of HIV-1 during combination anti-retroviral therapy.

A competitive, chaos-free, implicit, finite-difference method is developed and used for a novel deterministic model for the perturbation of HIV by combination antiretroviral therapy. The compartmental model monitors the interaction between HIV and CD4(+) T cells, its principal target and site of replication in vivo, in the presence of reverse transcription inhibitors and protease inhibitors. The model exhibits two steady states, an uninfected (trivial) steady state (with no virus present) and an endemically infected state (with virus and infected T cells present). Stability and bifurcation analyses together with numerical simulations of the resulting dynamical system are reported.

A sequential algorithm for the non-linear dual-sorption model of percutaneous drug absorption.

A sequential algorithm is developed for the non-linear dual-sorption model developed by Chandrasekaran et al. [1,2] which monitors pharmacokinetic profiles in percutaneous drug absorption. In the experimental study of percutaneous absorption, it is often observed that the lag-time decreases with the increase in the donor concentration when two or more donor concentrations of the same compound are used. The dual-sorption model has sometimes been employed to explain such experimental results. In this paper, it is shown that another feature observed after vehicle removal may also characterize the dual-sorption model. Soon after vehicle removal, the plots of the drug flux versus time become straight lines on a semilogarithmic scale as in the linear model, but the half-life is prolonged thereafter when the dual-sorption model prevails. The initial half-life after vehicle removal with a low donor concentration is longer than that with a higher donor concentration. These features, if observed in experiments, may be used as evidence to confirm that the dual-sorption model gives an explanation to the non-linear kinetic behaviour of a permeant.